In this question, we just consider square matrices.
The cofactor matrix $\mathrm{C}(\mathrm{A}) = (c_{ij})$ of a $n$-by-$n$ matrix $\mathrm{A} = (a_{ij})$ is a $n$-by-$n$ matrix ($n > 0$) with
$$c_{ij} = (-1)^{i+j} \cdot \det{(\mathrm{A}_{ij})} \; \forall i,j,$$
where $\mathrm{A}_{ij}$ is the matrix remaining after removing the $i$-th row and the $j$-th column from $\mathrm{A}$.
Problem: Let $\mathrm{A},\mathrm{B}$ be $n$-by-$n$ matrices. Prove that $\mathrm{C}(\mathrm{AB}) = \mathrm{C}(\mathrm{A})\cdot\mathrm{C}(\mathrm{B})$.
I start learning Linear Algebra recently, and this property is one of the problems giving me hard time proving. After giving up, I begin searching for a solution but not succeeding (my effort might not be enough). Therefore, I hope I can find a solution or a link to a solution here, and it's better that this solution uses no linear space and mapping knowledge. For further information about cofactor matrices and its relatives, you can read the Wikipedia article linked.
Thanks in advance!
Let $C(A)$ be matrix with entries $\alpha_{ij}$, and likewise $B(A)=(\beta_{ij})$. Then the entries $\gamma_{ij}$ of $C(A)C(B)$ can be written as $$ \gamma_{ij} = \sum_{k=1}^n \alpha_{ik} \beta_{kj} = \sum_{k=1}^n(-1)^{i+j+2k} \det(A_{ik}) \det(B_{kj}) = (-1)^{i+j}\sum_{k=1}^n\det(A_{ik}) \det(B_{kj}). $$ One can check that $$ (AB)_{ij} = A_{i\cdot}B_{\cdot j}, $$ where $A_{i\cdot}$ and $B_{\cdot j}$ are submatrices of $A$ and $B$ without the $i$th row and $j$th columns respectively, $A_{i\cdot}\in\mathbb R^{n-1,n}$, $B_{\cdot j}\in\mathbb R^{n,n-1}$ (or any other field instead of $\mathbb R$).
By the Cauchy-Binet formula, $$ \det(A_{i\cdot}B_{\cdot j})= \sum_{k=1}^n\det(A_{ik}) \det(B_{kj}). $$ This proves $$ \gamma_{ij} = (-1)^{i+j}\sum_{k=1}^n\det(A_{ik}) \det(B_{kj}) =(-1)^{i+j}\det((AB)_{ij} ), $$ hence the entries of $C(AB)$ and $C(A)C(B)$ conincide.