In my introductory graduate analysis class, we learned the Hardy-Littlewood Sobolev inequality (in dimension $d$): if $p, q > 1$, $0 < \lambda < d$ are such that $$ \frac{1}{p} + \frac{1}{q} + \frac{\lambda}{d} = 2, $$ then for all compactly supported smooth functions on $\mathbb{R}^d$, $$ \int_{\mathbb{R}^d}\int_{\mathbb{R}^d} f(x) |x - y|^{-\lambda}g(y)\, dx\, dy \leq C\|{f}\|_p \|{g}\|_q $$ for some constant $C$. We were then told that Sobolev's inequality, $$ \|u\|_{L^{p^*}} \leq C\|Du\|_{L^p} $$ for some (different than above, perhaps) constant $C$, was a consequence of Hardy-Littlewood-Sobolev. However, I had no idea how to show this. Firstly, the $L^{p^*}$ norm of a function $u$ (here, $p^{*}$ is the Sobolev conjugate of $p$) is a single integral, and I don't understand how we can recover this from the double integral in H-L-S.
My only idea was to somehow relate the $|x - y|$ to the difference quotient of $u$, but I don't see how to do this (or if this is even right). Even with this, I have no idea how to pick the right $f$ and $g$.
Any help in the right direction would be appreciated.