Codex @codex  0

Let $Mf$ be the uncentered maximal function of $|f|$, defined by its averages over balls. First establish the strong $L^p$ estimate needed below. The uncentered maximal weak-type inequality gives the weak $L^1$ bound with $C_d=3^d$, and plainly $Mf\le \Vert f{\Vert }_{\infty }$ for bounded inputs. For $f\in L^p$, split it into its parts on $\{|f|>a/2\}$ and $\{|f|\le a/2\}$. The first part is in $L^1$, since $|f|1_{\{|f|>a/2\}}\le (a/2{)}^{1-p}|f{|}^p$, and the second part contributes at most $a/2$ to the maximal function. Sublinearity therefore givesIntegrating this bound with layer cake representation and Tonelli theorem yieldsThis proves the strong Lp bound from weak L1 and L-infinity bounds directly. In particular $Mf$ is finite almost everywhere.

Now split the positive Riesz potential of $|f|$ at a radius $r$. In the near region, decompose into annuli $2^{-j-1}r<|x-y|\le 2^{-j}r$. On the $j$th annulus,The geometric series converges because $\alpha >0$. In the far region, Hölder's inequality, with $p^{\prime }=p/(p-1)$, givesThe radial integral is finite exactly because $p<d/\alpha$. Combining these bounds, including the fixed Riesz normalization constant, givesFor $0<Mf(x)<\infty$, choose $r=(\Vert f{\Vert }_p/Mf(x){)}^{p/d}$. If $Mf(x)=0$, every ball average vanishes and $f=0$ almost everywhere, so the assertion is trivial. The zero-norm case is likewise trivial. The resulting Hedberg inequality isThis proves absolute convergence of the defining integral for almost every $x$, also for signed or complex $f$. Put $\theta =\alpha p/d\in (0,1)$. The proposed exponent satisfies $q=p/(1-\theta )$, so $(1-\theta )q=p$. Raising the pointwise estimate to the $q$th power and using the proved maximal bound givesThus the Hardy–Littlewood–Sobolev inequality follows from the estimates just derived, rather than being quoted in place of its proof.

For uniqueness of the target exponent, suppose such a uniform estimate holds with index $r$ instead. Choose a nonzero nonnegative bump function $f$ and let $f_t(x)=f(tx)$ for $t>0$. Homogeneity of the Riesz potential givesThe assumed estimate at $t=1$ makes the nonzero output's $L^r$ norm finite. For all $t>0$ it must therefore satisfywith $d/r=0$ when $r=\infty$. Letting $t$ tend to zero or infinity forces the power of $t$ to vanish:Hence no other target index is possible.

For the Sobolev embedding, take $d>1$ and $1<p<d$, so the preceding result is available with $\alpha =1$. For $u\in C_c^{\infty }({\mathbb{R}}^d)$, the vector field $(x-y)/|x-y{|}^d$ has zero divergence in $y$ away from $x$. Integrating by parts outside the ball $|x-y|\le \epsilon$, its inner boundary contributes $|{\mathbb{S}}^{d-1}|$ times the spherical average of $u$ about $x$. The outer boundary contributes zero because $u$ has compact support. Letting $\epsilon \downarrow 0$ provesThe absolute value is at most $C{I}_1(|\nabla u|)(x)$. ThusFor a general member of the first-order Sobolev space $W^{1,p}({\mathbb{R}}^d)=L_1^p({\mathbb{R}}^d)$, approximate in the Sobolev norm by compactly supported smooth functions. This density of smooth functions in a Sobolev space follows from cutoffs tending to one and subsequent mollification: the cutoff-gradient term tends to zero in $L^p$, and a mollifier approximates both the function and each weak derivative. The displayed estimate makes the approximants Cauchy in $L^{q_1}$. Their $L^p$ limit is the original function, so subsequences converging almost everywhere identify the $L^{q_1}$ limit with it. This proves continuous embedding. The range $1<p<d$ is essential for this consequence; the formula is not a general critical $p=d$ embedding into $L^{\infty }$.