diff --git a/Mathlib.lean b/Mathlib.lean
index 3a6ceab1f6a98..e19af77f39161 100644
--- a/Mathlib.lean
+++ b/Mathlib.lean
@@ -3893,6 +3893,7 @@ import Mathlib.NumberTheory.EulerProduct.ExpLog
 import Mathlib.NumberTheory.FLT.Basic
 import Mathlib.NumberTheory.FLT.Four
 import Mathlib.NumberTheory.FLT.MasonStothers
+import Mathlib.NumberTheory.FLT.Polynomial
 import Mathlib.NumberTheory.FLT.Three
 import Mathlib.NumberTheory.FactorisationProperties
 import Mathlib.NumberTheory.Fermat
diff --git a/Mathlib/NumberTheory/FLT/MasonStothers.lean b/Mathlib/NumberTheory/FLT/MasonStothers.lean
index eee4ff81ac138..bd5d8a2c3a537 100644
--- a/Mathlib/NumberTheory/FLT/MasonStothers.lean
+++ b/Mathlib/NumberTheory/FLT/MasonStothers.lean
@@ -16,9 +16,6 @@ Proof is based on this online note by Franz Lemmermeyer http://www.fen.bilkent.e
 which is essentially based on Noah Snyder's paper "An Alternative Proof of Mason's Theorem",
 but slightly different.
 
-## TODO
-
-Prove polynomial FLT using Mason-Stothers theorem.
 -/
 
 open Polynomial UniqueFactorizationMonoid UniqueFactorizationDomain EuclideanDomain
@@ -51,13 +48,25 @@ private theorem abc_subcall {a b c w : k[X]} {hw : w ≠ 0} (wab : w = wronskian
       exact Nat.add_lt_add_right abc_dr_ndeg_lt _
 
 /-- **Polynomial ABC theorem.** -/
-theorem Polynomial.abc {a b c : k[X]} (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0) (hab : IsCoprime a b)
-    (hbc : IsCoprime b c) (hca : IsCoprime c a) (hsum : a + b + c = 0) :
+protected theorem Polynomial.abc
+    {a b c : k[X]} (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0)
+    (hab : IsCoprime a b) (hsum : a + b + c = 0) :
     ( natDegree a + 1 ≤ (radical (a * b * c)).natDegree ∧
       natDegree b + 1 ≤ (radical (a * b * c)).natDegree ∧
       natDegree c + 1 ≤ (radical (a * b * c)).natDegree ) ∨
       derivative a = 0 ∧ derivative b = 0 ∧ derivative c = 0 := by
   set w := wronskian a b with wab
+  have hbc : IsCoprime b c := by
+    rw [add_eq_zero_iff_neg_eq] at hsum
+    rw [← hsum, IsCoprime.neg_right_iff]
+    convert IsCoprime.add_mul_left_right hab.symm 1
+    rw [mul_one]
+  have hsum' : b + c + a = 0 := by rwa [add_rotate] at hsum
+  have hca : IsCoprime c a := by
+    rw [add_eq_zero_iff_neg_eq] at hsum'
+    rw [← hsum', IsCoprime.neg_right_iff]
+    convert IsCoprime.add_mul_left_right hbc.symm 1
+    rw [mul_one]
   have wbc : w = wronskian b c := wronskian_eq_of_sum_zero hsum
   have wca : w = wronskian c a := by
     rw [add_rotate] at hsum
diff --git a/Mathlib/NumberTheory/FLT/Polynomial.lean b/Mathlib/NumberTheory/FLT/Polynomial.lean
new file mode 100644
index 0000000000000..e6beb32165d5a
--- /dev/null
+++ b/Mathlib/NumberTheory/FLT/Polynomial.lean
@@ -0,0 +1,275 @@
+/-
+Copyright (c) 2024 Jineon Baek and Seewoo Lee. All rights reserved.
+Released under Apache 2.0 license as described in the file LICENSE.
+Authors: Jineon Baek, Seewoo Lee
+-/
+import Mathlib.Algebra.Polynomial.Expand
+import Mathlib.Algebra.GroupWithZero.Defs
+import Mathlib.NumberTheory.FLT.Basic
+import Mathlib.NumberTheory.FLT.MasonStothers
+import Mathlib.RingTheory.Polynomial.Content
+
+/-!
+# Fermat's Last Theorem for polynomials over a field
+
+This file states and proves the Fermat's Last Theorem for polynomials over a field.
+For `n ≥ 3` not divisible by the characteristic of the coefficient field `k` and (pairwise) nonzero
+coprime polynomials `a, b, c` (over a field) with `a ^ n + b ^ n = c ^ n`,
+all polynomials must be constants.
+
+More generally, we can prove non-solvability of the Fermat-Catalan equation: there are no
+non-constant polynomial solutions to the equation `u * a ^ p + v * b ^ q + w * c ^ r = 0`, where
+`p, q, r ≥ 3` with `p * q + q * r + r * p ≤ p * q * r` , `p, q, r` not divisible by `char k`,
+and `u, v, w` are nonzero elements in `k`.
+FLT is the special case where `p = q = r = n`, `u = v = 1`, and `w = -1`.
+
+The proof uses the Mason-Stothers theorem (Polynomial ABC theorem) and infinite descent
+(in the characteristic p case).
+-/
+
+open Polynomial UniqueFactorizationMonoid UniqueFactorizationDomain
+
+variable {k R : Type*} [Field k] [CommRing R] [IsDomain R] [NormalizationMonoid R]
+  [UniqueFactorizationMonoid R]
+
+private lemma Ne.isUnit_C {u : k} (hu : u ≠ 0) : IsUnit (C u) :=
+  Polynomial.isUnit_C.mpr hu.isUnit
+
+-- auxiliary lemma that 'rotates' coprimality
+private lemma rot_coprime
+    {p q r : ℕ} {a b c : k[X]} {u v w : k}
+    {hp : 0 < p} {hq : 0 < q} {hr : 0 < r}
+    {hu : u ≠ 0} {hv : v ≠ 0} {hw : w ≠ 0}
+    (heq : C u * a ^ p + C v * b ^ q + C w * c ^ r = 0) (hab : IsCoprime a b) : IsCoprime b c := by
+  have hCu : IsUnit (C u) := Ne.isUnit_C hu
+  have hCv : IsUnit (C v) := Ne.isUnit_C hv
+  have hCw : IsUnit (C w) := Ne.isUnit_C hw
+  rw [← IsCoprime.pow_iff hp hq, ← isCoprime_mul_units_left hCu hCv] at hab
+  rw [add_eq_zero_iff_neg_eq] at heq
+  rw [← IsCoprime.pow_iff hq hr, ← isCoprime_mul_units_left hCv hCw,
+    ← heq, IsCoprime.neg_right_iff]
+  convert IsCoprime.add_mul_left_right hab.symm 1 using 2
+  rw [mul_one]
+
+private lemma ineq_pqr_contradiction {p q r a b c : ℕ}
+    (hp : 0 < p) (hq : 0 < q) (hr : 0 < r)
+    (hineq : q * r + r * p + p * q ≤ p * q * r)
+    (hpa : p * a < a + b + c)
+    (hqb : q * b < a + b + c)
+    (hrc : r * c < a + b + c) : False := by
+  suffices h : p * q * r * (a + b + c) + 1 ≤ p * q * r * (a + b + c) by simp at h
+  calc
+    _ = (p * a) * (q * r) + (q * b) * (r * p) + (r * c) * (p * q) + 1 := by ring
+    _ ≤ (a + b + c) * (q * r) + (a + b + c) * (r * p) + (a + b + c) * (p * q) := by
+      rw [Nat.succ_le]
+      gcongr
+    _ = (q * r + r * p + p * q) * (a + b + c) := by ring
+    _ ≤ _ := by gcongr
+
+private theorem Polynomial.flt_catalan_deriv [DecidableEq k]
+    {p q r : ℕ} (hp : 0 < p) (hq : 0 < q) (hr : 0 < r)
+    (hineq : q * r + r * p + p * q ≤ p * q * r)
+    (chp : (p : k) ≠ 0) (chq : (q : k) ≠ 0) (chr : (r : k) ≠ 0)
+    {a b c : k[X]} (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0)
+    (hab : IsCoprime a b) {u v w : k} (hu : u ≠ 0) (hv : v ≠ 0) (hw : w ≠ 0)
+    (heq : C u * a ^ p + C v * b ^ q + C w * c ^ r = 0) :
+    derivative a = 0 ∧ derivative b = 0 ∧ derivative c = 0 := by
+  have hbc : IsCoprime b c := by apply rot_coprime heq <;> assumption
+  have hca : IsCoprime c a := by
+    rw [add_rotate] at heq; apply rot_coprime heq <;> assumption
+  have hCu := C_ne_zero.mpr hu
+  have hCv := C_ne_zero.mpr hv
+  have hCw := C_ne_zero.mpr hw
+  have hap := pow_ne_zero p ha
+  have hbq := pow_ne_zero q hb
+  have hcr := pow_ne_zero r hc
+  have habp : IsCoprime (C u * a ^ p) (C v * b ^ q) := by
+    rw [isCoprime_mul_units_left hu.isUnit_C hv.isUnit_C]; exact hab.pow
+  have hbcp : IsCoprime (C v * b ^ q) (C w * c ^ r) := by
+    rw [isCoprime_mul_units_left hv.isUnit_C hw.isUnit_C]; exact hbc.pow
+  have hcap : IsCoprime (C w * c ^ r) (C u * a ^ p) := by
+    rw [isCoprime_mul_units_left hw.isUnit_C hu.isUnit_C]; exact hca.pow
+  have habcp := hcap.symm.mul_left hbcp
+
+  -- Use Mason-Stothers theorem
+  rcases Polynomial.abc
+      (mul_ne_zero hCu hap) (mul_ne_zero hCv hbq) (mul_ne_zero hCw hcr)
+      habp heq with nd_lt | dr0
+  · simp_rw [radical_mul habcp, radical_mul habp,
+        radical_mul_of_isUnit_left hu.isUnit_C,
+        radical_mul_of_isUnit_left hv.isUnit_C,
+        radical_mul_of_isUnit_left hw.isUnit_C,
+        radical_pow a hp, radical_pow b hq, radical_pow c hr,
+        natDegree_mul hCu hap,
+        natDegree_mul hCv hbq,
+        natDegree_mul hCw hcr,
+        natDegree_C, natDegree_pow, zero_add,
+        ← radical_mul hab,
+        ← radical_mul (hca.symm.mul_left hbc)] at nd_lt
+
+    obtain ⟨hpa', hqb', hrc'⟩ := nd_lt
+    have hpa := hpa'.trans natDegree_radical_le
+    have hqb := hqb'.trans natDegree_radical_le
+    have hrc := hrc'.trans natDegree_radical_le
+    rw [natDegree_mul (mul_ne_zero ha hb) hc,
+        natDegree_mul ha hb, Nat.add_one_le_iff] at hpa hqb hrc
+    exfalso
+    exact (ineq_pqr_contradiction hp hq hr hineq hpa hqb hrc)
+  · rw [derivative_C_mul, derivative_C_mul, derivative_C_mul,
+        mul_eq_zero_iff_left (C_ne_zero.mpr hu),
+        mul_eq_zero_iff_left (C_ne_zero.mpr hv),
+        mul_eq_zero_iff_left (C_ne_zero.mpr hw),
+        derivative_pow_eq_zero chp,
+        derivative_pow_eq_zero chq,
+        derivative_pow_eq_zero chr] at dr0
+    exact dr0
+
+-- helper lemma that gives a baggage of small facts on `contract (ringChar k) a`
+private lemma find_contract {a : k[X]}
+    (ha : a ≠ 0) (hda : derivative a = 0) (chn0 : ringChar k ≠ 0) :
+    ∃ ca, ca ≠ 0 ∧
+      a = expand k (ringChar k) ca ∧
+      a.natDegree = ca.natDegree * ringChar k := by
+  have heq := (expand_contract (ringChar k) hda chn0).symm
+  refine ⟨contract (ringChar k) a, ?_, heq, ?_⟩
+  · intro h
+    rw [h, map_zero] at heq
+    exact ha heq
+  · rw [← natDegree_expand, ← heq]
+
+variable [DecidableEq k]
+
+private theorem Polynomial.flt_catalan_aux
+    {p q r : ℕ} {a b c : k[X]} {u v w : k}
+    (heq : C u * a ^ p + C v * b ^ q + C w * c ^ r = 0)
+    (hp : 0 < p) (hq : 0 < q) (hr : 0 < r)
+    (hineq : q * r + r * p + p * q ≤ p * q * r)
+    (chp : (p : k) ≠ 0) (chq : (q : k) ≠ 0) (chr : (r : k) ≠ 0)
+    (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0) (hab : IsCoprime a b)
+    (hu : u ≠ 0) (hv : v ≠ 0) (hw : w ≠ 0) :
+    a.natDegree = 0 := by
+  cases' eq_or_ne (ringChar k) 0 with ch0 chn0
+  -- characteristic zero
+  · obtain ⟨da, -, -⟩ := flt_catalan_deriv
+      hp hq hr hineq chp chq chr ha hb hc hab hu hv hw heq
+    have czk : CharZero k := by
+      apply charZero_of_inj_zero
+      intro n
+      rw [ringChar.spec, ch0]
+      exact zero_dvd_iff.mp
+    rw [eq_C_of_derivative_eq_zero da]
+    exact natDegree_C _
+  -- characteristic p > 0
+  · generalize eq_d : a.natDegree = d
+    -- set up infinite descent
+    -- strong induct on `d := a.natDegree`
+    induction d
+      using Nat.case_strong_induction_on
+      generalizing a b c ha hb hc hab heq with
+    | hz => rfl
+    | hi d ih_d => -- have derivatives of `a, b, c` zero
+      obtain ⟨ad, bd, cd⟩ := flt_catalan_deriv
+        hp hq hr hineq chp chq chr ha hb hc hab hu hv hw heq
+      -- find contracts `ca, cb, cc` so that `a(k) = ca(k^ch)`
+      obtain ⟨ca, ca_nz, eq_a, eq_deg_a⟩ := find_contract ha ad chn0
+      obtain ⟨cb, cb_nz, eq_b, eq_deg_b⟩ := find_contract hb bd chn0
+      obtain ⟨cc, cc_nz, eq_c, eq_deg_c⟩ := find_contract hc cd chn0
+      set ch := ringChar k
+      suffices hca : ca.natDegree = 0 by
+        rw [← eq_d, eq_deg_a, hca, zero_mul]
+      by_contra hnca; apply hnca
+      apply ih_d _ _ _ ca_nz cb_nz cc_nz <;> clear ih_d <;> try rfl
+      · apply (isCoprime_expand chn0).mp
+        rwa [← eq_a, ← eq_b]
+      · have _ : ch ≠ 1 := CharP.ringChar_ne_one
+        have hch2 : 2 ≤ ch := by omega
+        rw [← add_le_add_iff_right 1, ← eq_d, eq_deg_a]
+        refine le_trans ?_ (Nat.mul_le_mul_left _ hch2)
+        omega
+      · rw [eq_a, eq_b, eq_c, ← expand_C ch u, ← expand_C ch v, ← expand_C ch w] at heq
+        simp_rw [← map_pow, ← map_mul, ← map_add] at heq
+        rwa [Polynomial.expand_eq_zero (zero_lt_iff.mpr chn0)] at heq
+
+/-- Nonsolvability of the Fermat-Catalan equation. -/
+theorem Polynomial.flt_catalan
+    {p q r : ℕ} (hp : 0 < p) (hq : 0 < q) (hr : 0 < r)
+    (hineq : q * r + r * p + p * q ≤ p * q * r)
+    (chp : (p : k) ≠ 0) (chq : (q : k) ≠ 0) (chr : (r : k) ≠ 0)
+    {a b c : k[X]} (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0) (hab : IsCoprime a b)
+    {u v w : k} (hu : u ≠ 0) (hv : v ≠ 0) (hw : w ≠ 0)
+    (heq : C u * a ^ p + C v * b ^ q + C w * c ^ r = 0) :
+    a.natDegree = 0 ∧ b.natDegree = 0 ∧ c.natDegree = 0 := by
+  -- WLOG argument: essentially three times flt_catalan_aux
+  have hbc : IsCoprime b c := by
+    apply rot_coprime heq hab <;> assumption
+  have heq' : C v * b ^ q + C w * c ^ r + C u * a ^ p = 0 := by
+    rwa [add_rotate] at heq
+  have hca : IsCoprime c a := by
+    apply rot_coprime heq' hbc <;> assumption
+  refine ⟨?_, ?_, ?_⟩
+  · apply flt_catalan_aux heq <;> assumption
+  · rw [add_rotate] at heq hineq
+    rw [mul_rotate] at hineq
+    apply flt_catalan_aux heq <;> assumption
+  · rw [← add_rotate] at heq hineq
+    rw [← mul_rotate] at hineq
+    apply flt_catalan_aux heq <;> assumption
+
+theorem Polynomial.flt
+    {n : ℕ} (hn : 3 ≤ n) (chn : (n : k) ≠ 0)
+    {a b c : k[X]} (ha : a ≠ 0) (hb : b ≠ 0) (hc : c ≠ 0)
+    (hab : IsCoprime a b) (heq : a ^ n + b ^ n = c ^ n) :
+    a.natDegree = 0 ∧ b.natDegree = 0 ∧ c.natDegree = 0 := by
+  have hn' : 0 < n := by linarith
+  rw [← sub_eq_zero, ← one_mul (a ^ n), ← one_mul (b ^ n), ← one_mul (c ^ n), sub_eq_add_neg,
+    ← neg_mul] at heq
+  have hone : (1 : k[X]) = C 1 := by rfl
+  have hneg_one : (-1 : k[X]) = C (-1) := by simp only [map_neg, map_one]
+  simp_rw [hneg_one, hone] at heq
+  apply flt_catalan hn' hn' hn' _
+    chn chn chn ha hb hc hab one_ne_zero one_ne_zero (neg_ne_zero.mpr one_ne_zero) heq
+  have eq_lhs : n * n + n * n + n * n = 3 * n * n := by ring
+  rw [eq_lhs, mul_assoc, mul_assoc]
+  exact Nat.mul_le_mul_right (n * n) hn
+
+theorem fermatLastTheoremWith'_polynomial {n : ℕ} (hn : 3 ≤ n) (chn : (n : k) ≠ 0) :
+    FermatLastTheoremWith' k[X] n := by
+  rw [FermatLastTheoremWith']
+  intros a b c ha hb hc heq
+  obtain ⟨a', eq_a⟩ := gcd_dvd_left a b
+  obtain ⟨b', eq_b⟩ := gcd_dvd_right a b
+  set d := gcd a b
+  have hd : d ≠ 0 := gcd_ne_zero_of_left ha
+  rw [eq_a, eq_b, mul_pow, mul_pow, ← mul_add] at heq
+  have hdc : d ∣ c := by
+    have hn : 0 < n := by omega
+    have hdncn : d ^ n ∣ c ^ n := ⟨_, heq.symm⟩
+
+    rw [dvd_iff_normalizedFactors_le_normalizedFactors hd hc]
+    rw [dvd_iff_normalizedFactors_le_normalizedFactors
+          (pow_ne_zero n hd) (pow_ne_zero n hc),
+        normalizedFactors_pow, normalizedFactors_pow] at hdncn
+    simp_rw [Multiset.le_iff_count, Multiset.count_nsmul,
+      mul_le_mul_left hn] at hdncn ⊢
+    exact hdncn
+  obtain ⟨c', eq_c⟩ := hdc
+  rw [eq_a, mul_ne_zero_iff] at ha
+  rw [eq_b, mul_ne_zero_iff] at hb
+  rw [eq_c, mul_ne_zero_iff] at hc
+  rw [mul_comm] at eq_a eq_b eq_c
+  refine ⟨d, a', b', c', ⟨eq_a, eq_b, eq_c⟩, ?_⟩
+  rw [eq_c, mul_pow, mul_comm, mul_left_inj' (pow_ne_zero n hd)] at heq
+  suffices goal : a'.natDegree = 0 ∧ b'.natDegree = 0 ∧ c'.natDegree = 0 by
+    simp [natDegree_eq_zero] at goal
+    obtain ⟨⟨ca', ha'⟩, ⟨cb', hb'⟩, ⟨cc', hc'⟩⟩ := goal
+    rw [← ha', ← hb', ← hc']
+    rw [← ha', C_ne_zero] at ha
+    rw [← hb', C_ne_zero] at hb
+    rw [← hc', C_ne_zero] at hc
+    exact ⟨ha.right.isUnit_C, hb.right.isUnit_C, hc.right.isUnit_C⟩
+  apply flt hn chn ha.right hb.right hc.right _ heq
+  convert isCoprime_div_gcd_div_gcd _
+  · exact EuclideanDomain.eq_div_of_mul_eq_left ha.left eq_a.symm
+  · exact EuclideanDomain.eq_div_of_mul_eq_left ha.left eq_b.symm
+  · rw [eq_b]
+    exact mul_ne_zero hb.right hb.left