We study the radial extremal function $\varphi(|{\,\cdot\,}|)$ arising in the problem of finding the sharp Nikolskii constant $\mathcal C_d$ in $\mathit{PW}_1^1(\mathbb R^d)$ for arbitrary dimension $d\ge1$. We prove a factorization $\varphi=\Phi_1\Phi_2$, where $\Phi_1$ and $\Phi_2$ are entire functions of exponential type $1/2$ satisfying a functional equation and second-order differential equations with polynomial coefficients. As a result, the original extremal problem is reduced to a one-dimensional spectral problem depending on at most $d+1$ parameters. We also obtain a zeta interpretation of the coefficients of the polynomial appearing in the functional equation and a multiplicative equilibrium condition for the zeros of the extremal function. These results can be used to construct several algorithms for computing $\mathcal C_d$.
Let $\varphi_{d}(|{\,\cdot\,}|)$ be the radial extremal function in the problem for the sharp Nikolskii constant $\mathcal C_d^{-1}=\inf \|f\|_{1}$ over functions $f\in\mathit{PW}\,_{1}^{1}(\mathbb R^{d})$, $f(0)=1$. For every odd dimension, we construct an entire function $\Phi$ of exponential type $1/2$ such that $\varphi_d(z)=\Phi(z)\Phi(-z)$, and $\Phi$ satisfies a quadratic functional equation and a second-order linear differential equation with polynomial coefficients. Thus, the problem of finding the extremal function is reduced to a spectral problem with finitely many parameters. This result extends a recent one-dimensional result, but uses a different method. For example, in dimension $d=3$ it leads to a seven-diagonal spectral scheme that allows us to compute $\mathcal C_3$ to high accuracy. Even dimensions remain open within this approach.