Solution (source code)

= Solution

Since $M^T=I-\mathbf a\mathbf b^T=I-N$, nilpotence already gives $(I+N)(I-N)=I$. Explicitly, in <index notation>,
$$
(LM^T)_{ij}=\sum_k(\delta_{ik}+a_ib_k)(\delta_{jk}-b_ja_k)=\delta_{ij}-b_ja_i+a_ib_j-a_ib_j\sum_ka_kb_k=\delta_{ij}.
$$
Thus $\boxed{LM^T=I}$, or equivalently $M^T=L^{-1}$. Here $\delta_{ij}$ is the <Kronecker delta>.