Prove that for any $n, k\geq2$ there exist nonsingular nondiagonal matrices $A_1, A_2, \dots, A_k \in M_n(\mathbb{R})$ such that $$A_1^{-1} + A_2^{-1} + \dots + A_k^{-1} = \left( A_1 + A_2 + \dots + A_k \right)^{-1}$$
For $n=2$, if we have $A^{-1} + B^{-1} = (A+B)^{-1}$, if I am not mistaken, we can prove that $$\det(A)=\det(B)=\det(A+B)$$ if $A,B\in M_n(\mathbb{R})$, so it would be natural to consider that $\det(A_1)=\det(A_2)=\dots=\det(A_k)$, but from here I don't have any idea what should I do. What should we do?
If $k$ is odd, then it is easy: take $A_{2m} = -A_{2m+1}$ for every $m\ge 1$ so that the whole expression reduces to $A_1^{-1} = A_1^{-1}$ that is true for every $A_1$ invertible.
If $k$ is even, then using the same trick as before, what's left is $$ A_1^{-1} + A_2^{-1} = (A_1+A_2)^{-1} $$ So if you solve it for $k=2$, you solve it for every $k$.
Let's prove that $n=3$ and $k=2$ produces an absurd. It turns out that $$ 2I + A_1A_2^{-1} + A_2A_1^{-1} = I $$ and if $X = A_1A_2^{-1}$ then $X^{-1} = A_2A_1^{-1}$ and $$ X^2 + X + I = 0. $$ This means that the minimal polynomial of $X$ is $x^2+x+1$, so the characteristic polynomial must be $(x^2+x+1)^m$ for some $m$, but then $n=2m=3$, that is impossible.