Solution (source code)

= Solution

The differentiable loss has gradient $X^T(X\beta-Y)/n$. The <subdifferential sum rule> and part (b) give
$$
\partial Q(\beta)=\frac1nX^T(X\beta-Y)+\lambda\sqrt m\,u,
$$
where blockwise
$$
u^{(k)}=\frac{\beta^{(k)}}{\lVert\beta^{(k)}\rVert_2}
\quad\text{if }\beta^{(k)}\ne0,
\qquad
\lVert u^{(k)}\rVert_2\leq1
\quad\text{if }\beta^{(k)}=0.
$$