Solution (source code)

= Solution

The reverse-regression estimate is $\hat b=\sum X_iY_i/\sum Y_i^2$. Therefore
$$
\hat b\hat\beta=\frac{(\sum X_iY_i)^2}{(\sum X_i^2)(\sum Y_i^2)}\le1
$$
by Cauchy–Schwarz. Equality means the two data <vectors> are proportional, which is exactly when both fitted residual sums, and hence both variance MLEs, vanish.