ConceptComplete

Cofactor Expansion

Cofactor expansion (also called Laplace expansion) is the standard recursive method for computing determinants. It expands the determinant along any row or column, reducing an n×nn \times n determinant to a sum of (n1)×(n1)(n-1) \times (n-1) determinants.


Definitions

Definition4.3Minor and cofactor

Let AMn×n(F)A \in M_{n \times n}(F).

The (i,j)(i, j)-minor MijM_{ij} is the determinant of the (n1)×(n1)(n-1) \times (n-1) submatrix obtained by deleting row ii and column jj from AA.

The (i,j)(i, j)-cofactor is

Cij=(1)i+jMij.C_{ij} = (-1)^{i+j} M_{ij}.

The sign pattern (1)i+j(-1)^{i+j} follows a checkerboard:

(+++++)\begin{pmatrix} + & - & + & \cdots \\ - & + & - & \cdots \\ + & - & + & \cdots \\ \vdots & \vdots & \vdots & \ddots \end{pmatrix}

Definition4.4Cofactor expansion

Expansion along row ii:

det(A)=j=1naijCij=j=1n(1)i+jaijMij.\det(A) = \sum_{j=1}^n a_{ij} C_{ij} = \sum_{j=1}^n (-1)^{i+j} a_{ij} M_{ij}.

Expansion along column jj:

det(A)=i=1naijCij=i=1n(1)i+jaijMij.\det(A) = \sum_{i=1}^n a_{ij} C_{ij} = \sum_{i=1}^n (-1)^{i+j} a_{ij} M_{ij}.

Both give the same result regardless of which row or column is chosen.


Examples

Example3 x 3 expansion along row 1

A=(213014521).A = \begin{pmatrix} 2 & 1 & 3 \\ 0 & -1 & 4 \\ 5 & 2 & 1 \end{pmatrix}.

Expanding along row 1:

det(A)=2det(1421)1det(0451)+3det(0152)\det(A) = 2 \cdot \det\begin{pmatrix} -1 & 4 \\ 2 & 1 \end{pmatrix} - 1 \cdot \det\begin{pmatrix} 0 & 4 \\ 5 & 1 \end{pmatrix} + 3 \cdot \det\begin{pmatrix} 0 & -1 \\ 5 & 2 \end{pmatrix}

=2(18)1(020)+3(0+5)=2(9)1(20)+3(5)=18+20+15=17.= 2(-1 - 8) - 1(0 - 20) + 3(0 + 5) = 2(-9) - 1(-20) + 3(5) = -18 + 20 + 15 = 17.

ExampleSame matrix, expansion along column 1

Column 1 has a zero entry at position (2,1)(2, 1), so one term vanishes:

det(A)=2C11+0C21+5C31\det(A) = 2 \cdot C_{11} + 0 \cdot C_{21} + 5 \cdot C_{31}

=2det(1421)+0+5(1)3+1det(1314)= 2 \cdot \det\begin{pmatrix} -1 & 4 \\ 2 & 1 \end{pmatrix} + 0 + 5 \cdot (-1)^{3+1}\det\begin{pmatrix} 1 & 3 \\ -1 & 4 \end{pmatrix}

=2(9)+5(4+3)=18+35=17.= 2(-9) + 5(4 + 3) = -18 + 35 = 17.

Same answer. Choosing a row/column with zeros saves work.

Example4 x 4 with strategic expansion

A=(1003200101200034).A = \begin{pmatrix} 1 & 0 & 0 & 3 \\ 2 & 0 & 0 & 1 \\ 0 & 1 & 2 & 0 \\ 0 & 0 & 3 & 4 \end{pmatrix}.

Expand along column 2 (two zeros): only rows 3 contributes.

det(A)=(1)3+21det(103201034)=det(103201034).\det(A) = (-1)^{3+2} \cdot 1 \cdot \det\begin{pmatrix} 1 & 0 & 3 \\ 2 & 0 & 1 \\ 0 & 3 & 4 \end{pmatrix} = -\det\begin{pmatrix} 1 & 0 & 3 \\ 2 & 0 & 1 \\ 0 & 3 & 4 \end{pmatrix}.

Expand this 3×33 \times 3 along column 2: only row 3 contributes.

=(1)3+23det(1321)=(3)(16)=(3)(5)=15.= -(-1)^{3+2} \cdot 3 \cdot \det\begin{pmatrix} 1 & 3 \\ 2 & 1 \end{pmatrix} = -(- 3)(1 - 6) = -(-3)(-5) = -15.

ExampleTriangular matrix via cofactor

For an upper triangular 3×33 \times 3 matrix, expanding along column 1 repeatedly:

det(a0b00c)=adet(b0c)=abc=abc.\det\begin{pmatrix} a & * & * \\ 0 & b & * \\ 0 & 0 & c \end{pmatrix} = a \cdot \det\begin{pmatrix} b & * \\ 0 & c \end{pmatrix} = a \cdot bc = abc.

Generalizing: the determinant of a triangular matrix is the product of diagonal entries.

ExampleBlock triangular determinant

For a block upper triangular matrix (AB0D)\begin{pmatrix} A & B \\ 0 & D \end{pmatrix} where AA is k×kk \times k and DD is ×\ell \times \ell:

det(AB0D)=det(A)det(D).\det\begin{pmatrix} A & B \\ 0 & D \end{pmatrix} = \det(A) \cdot \det(D).

This generalizes the triangular matrix result.

ExampleVandermonde determinant

The Vandermonde matrix with nodes x1,,xnx_1, \ldots, x_n:

V=(1x1x12x1n11x2x22x2n11xnxn2xnn1)V = \begin{pmatrix} 1 & x_1 & x_1^2 & \cdots & x_1^{n-1} \\ 1 & x_2 & x_2^2 & \cdots & x_2^{n-1} \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & x_n & x_n^2 & \cdots & x_n^{n-1} \end{pmatrix}

has det(V)=1i<jn(xjxi)\det(V) = \prod_{1 \leq i < j \leq n} (x_j - x_i). So VV is invertible iff all xix_i are distinct.

For n=3n = 3: det(V)=(x2x1)(x3x1)(x3x2)\det(V) = (x_2 - x_1)(x_3 - x_1)(x_3 - x_2).


The adjugate matrix

Definition4.5Adjugate (classical adjoint)

The adjugate (or classical adjoint) of AA is the transpose of the cofactor matrix:

adj(A)=(Cij)T,i.e., (adj(A))ij=Cji.\text{adj}(A) = (C_{ij})^T, \quad \text{i.e., } (\text{adj}(A))_{ij} = C_{ji}.

Theorem4.2Adjugate formula for inverse

For any n×nn \times n matrix AA:

Aadj(A)=adj(A)A=det(A)In.A \cdot \text{adj}(A) = \text{adj}(A) \cdot A = \det(A) \cdot I_n.

If AA is invertible:

A1=1det(A)adj(A).A^{-1} = \frac{1}{\det(A)} \text{adj}(A).

ExampleAdjugate of a 2 x 2 matrix

A=(abcd)A = \begin{pmatrix} a & b \\ c & d \end{pmatrix}, adj(A)=(dbca)\text{adj}(A) = \begin{pmatrix} d & -b \\ -c & a \end{pmatrix}, so A1=1adbc(dbca)A^{-1} = \frac{1}{ad - bc}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix}.

ExampleAdjugate of a 3 x 3 matrix

For A=(123014560)A = \begin{pmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 5 & 6 & 0 \end{pmatrix}:

det(A)=1(024)2(020)+3(05)=24+4015=1\det(A) = 1(0-24) - 2(0-20) + 3(0-5) = -24 + 40 - 15 = 1.

Computing all nine cofactors and transposing gives adj(A)=A1\text{adj}(A) = A^{-1} (since detA=1\det A = 1).

ExampleExpansion along 'alien' row

What happens if we expand using entries of row ii with cofactors of row kk (kik \neq i)?

j=1naijCkj=0(ik).\sum_{j=1}^n a_{ij} C_{kj} = 0 \quad (i \neq k).

This is because the sum equals the determinant of a matrix with two identical rows (ii and kk), which is zero. This fact is used to prove Aadj(A)=det(A)IA \cdot \text{adj}(A) = \det(A) \cdot I.

ExampleAdjugate of a singular matrix

If rank(A)=n1\text{rank}(A) = n - 1, then adj(A)\text{adj}(A) has rank 1. If rank(A)n2\text{rank}(A) \leq n - 2, then adj(A)=0\text{adj}(A) = 0. The adjugate captures the "next level" of singularity.

ExampleDeterminant from 2 x 2 minors

For a 3×33 \times 3 matrix, each cofactor is a 2×22 \times 2 determinant. Computing all nine cofactors involves nine 2×22 \times 2 determinants. This is more work than row reduction for large matrices, but the adjugate formula is theoretically important.

ExampleCofactor expansion for skew-symmetric matrices

For a 3×33 \times 3 skew-symmetric matrix A=(0aba0cbc0)A = \begin{pmatrix} 0 & a & b \\ -a & 0 & c \\ -b & -c & 0 \end{pmatrix}, expanding along row 1:

det(A)=0M11adet(acb0)+bdet(a0bc)=a(bc)+b(ac)=0.\det(A) = 0 \cdot M_{11} - a \cdot \det\begin{pmatrix} -a & c \\ -b & 0 \end{pmatrix} + b \cdot \det\begin{pmatrix} -a & 0 \\ -b & -c \end{pmatrix} = -a(bc) + b(ac) = 0.

Every odd-dimensional skew-symmetric matrix has determinant 00.


RemarkComputational complexity

Cofactor expansion has O(n!)O(n!) complexity -- impractical for n>10n > 10. For numerical computation, Gaussian elimination (O(n3)O(n^3)) is vastly superior. Cofactor expansion remains valuable for theoretical proofs, symbolic determinants, and small matrices.

RemarkLooking ahead

The determinant connects to invertibility in the Determinant and Invertibility Theorem, and the adjugate formula leads to Cramer's Rule for explicitly solving linear systems.