Solution (source code)

= Solution

Put $Z=[X,W]$. Its ridgeless minimum-norm fit is $Z^T(ZZ^T)^+Y$, so its first $p$ coordinates are
$$
\widetilde\beta_{1:p}=X^T(XX^T+WW^T)^+Y.
$$
The <strong law of large numbers> applied entrywise gives $WW^T\to\lambda I_n$ almost surely as $d\to\infty$. Continuity of inversion then yields
$$
\widetilde\beta_{1:p}\longrightarrow
X^T(XX^T+\lambda I_n)^{-1}Y
=(X^TX+\lambda I_p)^{-1}X^TY=\widehat\beta_\lambda
$$
almost surely, where the middle equality is the <push-through identity>.