Solution (source code)

= Solution

For a unit <vector> $v$, a <Householder reflection> is $H=I-2vv^T$. Clearly $H^T=H$, and $H^TH=(I-2vv^T)^2=I$.