# Prove that a square matrix with integer entries can be decomposed into a product of square matrices with integer matrixes

Let $A$ be an $n ×n$ matrix with integer entries and $b_1,...,b_k$ are integers satisfying $det A = b_1•••b_k$. Show that there exist $n × n$ matrices $B_1,..., B_k$ with integer entries such that $A = B_1•••B_k$ and $det B_i = b_i$ for all $i = 1,...,k.$

It suffices to prove the case for $k=2$ and the generalized result follows from induction. I am not sure how to proceed here. We can surely decompose $A$ into a product of elementary matrices $A = C_1•••C_j$, but how can we ensure if they contain a consecutive product of $C_1,...,C_l$'s that satisfy $det C_1•••C_l = b_1$ and $det C_{l+1}•••C_j = b_2$? Also, there's no guarantee that those will be matrices with integer entries...

• This is maybe overkill, but are you familiar with Smith normal form (en.wikipedia.org/wiki/Smith_normal_form)? Sep 26, 2017 at 23:00
• no... and I am mot sure if I can easily comprehend the wiki article... Sep 26, 2017 at 23:05
• Take $B_1=A$ and $B_2=\dots = B_k=I_n$. Oct 1, 2017 at 22:13
• that doesn't work, as then $det B_i = 1$ for all $i >1$. What if $det A = 3*7*11*47$? Oct 1, 2017 at 22:18
• Then $\det B_1=3\times 7\times 11\times 47\in\mathbb{N}$ and for $i>1$, $\det B_i=1$. Maybe you want the $b_i$ to be a prime? Oct 2, 2017 at 3:27

$\def\diag{\operatorname{diag}}$ The suggestion of Qiaochu Yuan is just what you need. You can do this using Smith normal form, and you are not likely to be able to do this without using Smith normal form. I assume you are only interested in the case where the determinant is non zero. In that case, the theorem regarding Smith normal form is that there exist integer matrices $U$ and $V$ with determinant equal to $\pm 1$, and a diagonal integer matrix $D = \operatorname{diag}(d_1, d_2, \dots, d_n)$ such that $A = U D V$. (Moreover, the diagonal entries satisfy $d_j$ divides $d_{j+1}$ for all $j$.) Thus $\det(A) = d_1 d_2 \cdots d_n$. Now you can get your result by factoring the diagonal matrix $D$. You can factor $D$ into (mutually commuting) diagonal matrices, each with only one diagonal entry different from 1, and with that special entry equal to some prime factor of $\det(A)$. Then you can regroup these factors to obtain the desired matrix factors with determinants equal to your chosen factors $b_j$. For example, with $D = \diag(1, 4, 12)$, we have $\det(D) = 48$ and $D$ factors as $$D = \diag(1,2,1)\diag(1,2,1) \diag(1,1,2) \diag(1,1,2) \diag(1,1,3).$$ If your desired factors are $b_1 = 8$ and $b_2 = 6$, then you can re-group the factors above as $$D = \diag(1, 4, 2) \diag(1, 1, 6).$$ So finally, you would take $B_1 = U \diag(1,4,2)$ and $B_2 = \diag(1,1,6)V$. A small amount of additional fiddling might be necessary to take care of unwanted signs $\pm 1$, ask me if you are concerned about that. If you need to implement this method computationally, you can do that, because there is an algorithm for computing Smith normal form, i.e. for computing $U, V, D$.