Let , , and so on be as described in this previous exercise. Compute the relations matrix of with respect to the generating set for . Conclude that Theorem 21 and the algorithm described on page 480-481 in D&F for determining the invariant factor decomposition of an -module follow from the algorithm described in the series of exercises beginning here.

Recall that , where is the Kronecker delta. So , and thus by definition the relations matrix of is .

By the algorithm we discussed here, using a sequence of elementary row and column operations we can transform this relations matrix into a diagonal matrix whose diagonal entries are ordered by divisibility and are precisely the invariant factors of as an module. This diagonal matrix is the Smith Normal Form of . Recall that each elementary row or column operation can be achieved by left (row) or right (column) multiplying by an appropriate elementary matrix (as described here) and that these elementary matrices are invertible. So we have products and of elementary matrices such that . Now and are also products of elementary matrices, and we have . In particular, and have the same invariant factors, and so is also the Smith Normal Form of . This proves Theorem 21 in D&F.

Now we will verify the correctness of the algorithm given on page 480-481 of D&F. Fix the standard basis on an -dimensional vector space , let be a linear transformation on , and let be the matrix of with respect to the standard basis. By Theorem 21, there are products and of elementary matrices such that is in Smith Normal Form. ( is the row operations used in reverse order, and is the column operations used in order.) With respect to the generators , has relations matrix , and is in Smith Normal Form. Now recall (from here) that elementary column operations on a relations matrix correspond to ‘elementary perturbations’ of the generating set of . In particular, letting be the standard basis of as a matrix, performing these elementary perturbations in order corresponds precisely to the column operations described in part (2) of the algorithm on page 480. The resulting generating set consists of elements of -linear combinations of the , which we convert to -linear combinations of the via (that is, by ‘evaluating at ‘). The nonzero elements of the resulting generating set, , then generate (homomorphically) so that , and moreover the annihilator of is the corresponding invariant factor .

In particular, the elements with are -linearly independent (since is the monic polynomial of least degree which annihilates ), and so by dimension considerations are a basis for the -vector space . Let be the matrix whose columns are these basis vectors in order; note that is invertible by this previous exercise. Now if is the matrix for with respect to the standard bases, then is the matrix for with respect to the basis forming the columns of . (This is the standard ‘change of basis’ procedure.)We can see then (by evaluating on each ) that is precisely the direct sum of the companion matrices of the invariant factors . That is, is in rational canonical form.

[This is kind of rough. I would probably rearrange the material a bit, compared to the book’s ordering of results. No time to think about that now though.]