Cayley-Hamilton theorem

At a Glance

Why This Chapter Matters

Cayley-Hamilton delivers 10–20 marks across three different question styles: compute a high matrix power, find the inverse via the characteristic equation, or evaluate a matrix polynomial by polynomial division. All three reduce to the same workflow — compute the characteristic polynomial, apply Cayley-Hamilton to reduce powers, then evaluate numerically. The 3×33\times 3 matrix with repeated eigenvalue (appearing in both 2019 and 2023) adds one step: the double-root derivative condition. A student who has practised this template can solve any variant in under 15 minutes.

Minimum Theory

Cayley-Hamilton theorem. Every n×nn\times n matrix AA satisfies its own characteristic equation: if p(λ)=det⁡(λI−A)p(\lambda)=\det(\lambda I-A) then p(A)=Op(A)=O. This means every matrix power AkA^k (for k≥nk\ge n) can be expressed as a polynomial in AA of degree ≤n−1\le n-1.

Power reduction via Lagrange interpolation. To find ANA^N for large NN, write λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) where deg⁡r<n\deg r<n. Since p(A)=Op(A)=O, AN=r(A)A^N=r(A). Determine the nn coefficients of rr by substituting each eigenvalue (with multiplicity): a simple root λ0\lambda_0 gives λ0N=r(λ0)\lambda_0^N=r(\lambda_0); a double root λ1\lambda_1 gives both λ1N=r(λ1)\lambda_1^N=r(\lambda_1) and Nλ1N−1=r′(λ1)N\lambda_1^{N-1}=r'(\lambda_1) (differentiate both sides of λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) and use p(λ1)=p′(λ1)=0p(\lambda_1)=p'(\lambda_1)=0).

Inverse via Cayley-Hamilton. For 2×22\times 2: from A2−(tr A)A+(det⁡A)I=OA^2-(\text{tr}\,A)A+(\det A)I=O, rearrange to isolate II on one side and factor out AA: A−1=((tr A)I−A)/det⁡AA^{-1}=\bigl((\text{tr}\,A)I-A\bigr)/\det A.

Pipeline: from matrix A to A^n via the Cayley-Hamilton power-reduction technique

Question Archetypes

ArchetypeYou are seeing this when…
cayley-hamilton-power”Use the Cayley-Hamilton theorem to find ANA^N” for large NN
cayley-hamilton-application”Verify CH; find A−1A^{-1}; evaluate a matrix polynomial p(A)p(A)“

cayley-hamilton-power (2 question(s); 2019, 2023)

Recognition Cues

Solution Template

  1. State the theorem. Every n×nn\times n matrix satisfies its characteristic equation p(A)=Op(A)=O.
  2. Compute p(λ)=det⁡(λI−A)p(\lambda)=\det(\lambda I-A). Expand and factorise.
  3. Write λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) with rr of degree ≤2\le 2 (i.e.\ r=aλ2+bλ+cr=a\lambda^2+b\lambda+c). So AN=aA2+bA+cIA^N=aA^2+bA+cI.
  4. Find a,b,ca,b,c from eigenvalue conditions. Simple root: substitute; double root: substitute and differentiate.
  5. Compute A2A^2 entry-by-entry. Assemble AN=aA2+bA+cIA^N=aA^2+bA+cI.

Worked Example

2019 Paper 1, 2019-P1-Q4a (15 marks)

State the Cayley-Hamilton theorem. Use it to find A100A^{100} for

A=(100101010).A=\begin{pmatrix}1&0&0\\1&0&1\\0&1&0\end{pmatrix}.

Step 1 — State theorem. Every square matrix satisfies its own characteristic equation: p(A)=Op(A)=O.

Step 2 — Characteristic polynomial.

p(λ)=det⁡(λ−100−1λ−10−1λ)=(λ−1)(λ2−1)=(λ−1)2(λ+1).p(\lambda)=\det\begin{pmatrix}\lambda-1&0&0\\-1&\lambda&-1\\0&-1&\lambda\end{pmatrix}=(\lambda-1)(\lambda^2-1)=(\lambda-1)^2(\lambda+1).

So p(λ)=λ3−λ2−λ+1p(\lambda)=\lambda^3-\lambda^2-\lambda+1 and by CH, A3=A2+A−IA^3=A^2+A-I.

Step 3 — Write remainder. λ100=q(λ)(λ−1)2(λ+1)+(aλ2+bλ+c)\lambda^{100}=q(\lambda)(\lambda-1)^2(\lambda+1)+(a\lambda^2+b\lambda+c), so A100=aA2+bA+cIA^{100}=aA^2+bA+cI.

Step 4 — Eigenvalue conditions.

λ=1\lambda=1 (value): 1=a+b+c.(1)1=a+b+c. \qquad(1)

λ=1\lambda=1 (derivative, double root): 100=2a+b.(2)100=2a+b. \qquad(2)

λ=−1\lambda=-1 (simple root): (−1)100=1=a−b+c.(3)(-1)^{100}=1=a-b+c. \qquad(3)

From (1)−(3)(1)-(3): 2b=02b=0, so b=0b=0. From (2)(2): a=50a=50. From (1)(1): c=1−50=−49c=1-50=-49.

Step 5 — Compute A2A^2 and assemble.

A2=(100110101).A^2=\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}.

A100=50A2−49I=50(100110101)−49(100010001)=(10050105001).A^{100}=50A^2-49I=50\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}-49\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}=\begin{pmatrix}1&0&0\\50&1&0\\50&0&1\end{pmatrix}.

  A100=(10050105001).  \boxed{\;A^{100}=\begin{pmatrix}1&0&0\\50&1&0\\50&0&1\end{pmatrix}.\;}


2023 Paper 1, 2023-P1-Q3a (20 marks)

(i) Verify CH for A=(100101010)A=\begin{pmatrix}1&0&0\\1&0&1\\0&1&0\end{pmatrix}. (ii) Show An=An−2+A2−IA^n=A^{n-2}+A^2-I for n≥3n\ge 3 and find A40A^{40}.

Part (i) — Verify CH.

Characteristic polynomial: p(λ)=(λ−1)2(λ+1)=λ3−λ2−λ+1p(\lambda)=(\lambda-1)^2(\lambda+1)=\lambda^3-\lambda^2-\lambda+1.

Compute A2A^2 and A3A^3 entry-by-entry:

A2=(100110101),A3=(100201110).A^2=\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix},\qquad A^3=\begin{pmatrix}1&0&0\\2&0&1\\1&1&0\end{pmatrix}.

A3−A2−A+I=(1−1−1+1002−1−1+00−1−0+11−0−1+01−1−0+01−0−1+00−1−0+1)=(000000000).  ✓A^3-A^2-A+I=\begin{pmatrix}1-1-1+1&0&0\\2-1-1+0&0-1-0+1&1-0-1+0\\1-1-0+0&1-0-1+0&0-1-0+1\end{pmatrix}=\begin{pmatrix}0&0&0\\0&0&0\\0&0&0\end{pmatrix}.\;\checkmark

Part (ii) — Recurrence and A40A^{40}.

From CH: A3=A2+A−IA^3=A^2+A-I, i.e.\ A3−A=A2−IA^3-A=A^2-I, i.e.\ A3−A1=A2−IA^3-A^1=A^2-I.

Claim: An−An−2=A2−IA^n-A^{n-2}=A^2-I for all n≥3n\ge 3.

Proof by induction. Base n=3n=3: A3−A=A2−IA^3-A=A^2-I follows directly from CH. Inductive step: assume Ak−Ak−2=A2−IA^k-A^{k-2}=A^2-I. Multiply both sides on the left by AA: Ak+1−Ak−1=A3−A=A2−IA^{k+1}-A^{k-1}=A^3-A=A^2-I (using the base case again). This proves Ak+1−Ak−1=A2−IA^{k+1}-A^{k-1}=A^2-I.

Telescope to A40A^{40}: sum the identity An−An−2=A2−IA^n-A^{n-2}=A^2-I for n=4,6,8,…,40n=4,6,8,\ldots,40 (19 terms, all even steps):

A40−A2=19(A2−I)⇒A40=20A2−19I.A^{40}-A^2=19(A^2-I)\qquad\Rightarrow\qquad A^{40}=20A^2-19I.

A40=20(100110101)−19I=(10020102001).A^{40}=20\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}-19I=\begin{pmatrix}1&0&0\\20&1&0\\20&0&1\end{pmatrix}.

  A40=(10020102001).  \boxed{\;A^{40}=\begin{pmatrix}1&0&0\\20&1&0\\20&0&1\end{pmatrix}.\;}

Common Traps


cayley-hamilton-application (1 question(s); 2014)

Recognition Cues

Solution Template

  1. Characteristic polynomial. p(λ)=λ2−(tr A)λ+det⁡Ap(\lambda)=\lambda^2-(\text{tr}\,A)\lambda+\det A.
  2. Verify CH. Compute A2A^2 directly and check A2−(tr A)A+(det⁡A)I=OA^2-(\text{tr}\,A)A+(\det A)I=O.
  3. Inverse. Rearrange CH: A(A−(tr A)I)=−(det⁡A)IA(A-(\text{tr}\,A)I)=-(\det A)I, so A−1=((tr A)I−A)/det⁡AA^{-1}=((\text{tr}\,A)I-A)/\det A.
  4. Higher powers. Use A2=(tr A)A−(det⁡A)IA^2=(\text{tr}\,A)A-(\det A)I as the recurrence; compute A3,A4,A5,…A^3,A^4,A^5,\ldots iteratively as αkA+βkI\alpha_k A+\beta_k I.
  5. Sum the polynomial. Collect all AA-coefficients and all II-coefficients. Write the result as a single matrix.

Worked Example

2014 Paper 1, 2014-P1-Q2b-ii (10 marks)

Verify CH for A=(1423)A=\begin{pmatrix}1&4\\2&3\end{pmatrix}; find A−1A^{-1}; evaluate A5−4A4−7A3+11A2−A−10IA^5-4A^4-7A^3+11A^2-A-10I.

Step 1 — Characteristic polynomial. tr A=4\text{tr}\,A=4, det⁡A=3−8=−5\det A=3-8=-5. So p(λ)=λ2−4λ−5p(\lambda)=\lambda^2-4\lambda-5.

Step 2 — Verify CH. A2=(916817)A^2=\begin{pmatrix}9&16\\8&17\end{pmatrix}. Check: A2−4A−5I=(9−4−516−16−08−8−017−12−5)=OA^2-4A-5I=\begin{pmatrix}9-4-5&16-16-0\\8-8-0&17-12-5\end{pmatrix}=O. ✓\checkmark

Step 3 — Inverse. From A2−4A−5I=OA^2-4A-5I=O: A(A−4I)=5IA(A-4I)=5I, so

A−1=A−4I5=15(−342−1).A^{-1}=\frac{A-4I}{5}=\frac{1}{5}\begin{pmatrix}-3&4\\2&-1\end{pmatrix}.

Step 4 — Reduce higher powers. Recurrence: A2=4A+5IA^2=4A+5I.

A3=4A2+5A=4(4A+5I)+5A=21A+20I.A^3=4A^2+5A=4(4A+5I)+5A=21A+20I.

A4=A⋅A3=21A2+20A=21(4A+5I)+20A=104A+105I.A^4=A\cdot A^3=21A^2+20A=21(4A+5I)+20A=104A+105I.

A5=A⋅A4=104A2+105A=104(4A+5I)+105A=521A+520I.A^5=A\cdot A^4=104A^2+105A=104(4A+5I)+105A=521A+520I.

Step 5 — Evaluate polynomial. Collect:

TermAA-coeffII-coeff
A5A^5521521520520
−4A4-4A^4−416-416−420-420
−7A3-7A^3−147-147−140-140
11A211A^244445555
−A-A−1-100
−10I-10I00−10-10
Sum1155

Result: A+5I=(6428)A+5I=\begin{pmatrix}6&4\\2&8\end{pmatrix}.

  A5−4A4−7A3+11A2−A−10I=(6428).  \boxed{\;A^5-4A^4-7A^3+11A^2-A-10I=\begin{pmatrix}6&4\\2&8\end{pmatrix}.\;}

(Alternative: polynomial division of λ5−4λ4−7λ3+11λ2−λ−10\lambda^5-4\lambda^4-7\lambda^3+11\lambda^2-\lambda-10 by λ2−4λ−5\lambda^2-4\lambda-5 gives remainder λ+5\lambda+5, so p(A)=A+5Ip(A)=A+5I.)

Common Traps


Marks-Aware Writing

10-mark question (2014): Verification (Steps 1–2), inverse (Step 3), and polynomial evaluation (Steps 4–5) each carry roughly 3+3+43+3+4 marks. An answer that verifies CH but miscomputes the inverse loses 3 marks; missing the polynomial evaluation entirely loses 4.

15-mark question (2019): Stating CH (2 marks), characteristic polynomial (3 marks), setting up and solving the a,b,ca,b,c system (5 marks including the double-root derivative condition), computing A2A^2 and the final matrix (5 marks).

20-mark question (2023): Verification (10 marks) and recurrence + A40A^{40} (10 marks). For verification: A2A^2 and A3A^3 computed correctly (5 marks), substitution and confirmation of zero (5 marks). For part (ii): proving the recurrence by induction (5 marks), telescoping correctly to get A40=20A2−19IA^{40}=20A^2-19I (3 marks), final matrix (2 marks).

Practice Set

(No additional practice items beyond the three worked examples.)

Ready to drill what you just read?

Daily Practice turns these patterns into one adaptive set a day — practised daily until they're automatic, free for everyone.

See Daily Practice →

This chapter is part of the Maths Coverage Map — 14 years, mapped. Get the take-away PDF free.