diff --git a/Mathlib.lean b/Mathlib.lean
index 3192355b35fd8c..583df52fc94b5c 100644
--- a/Mathlib.lean
+++ b/Mathlib.lean
@@ -3219,6 +3219,7 @@ public import Mathlib.Combinatorics.Enumerative.Schroder
 public import Mathlib.Combinatorics.Enumerative.Stirling
 public import Mathlib.Combinatorics.Extremal.RuzsaSzemeredi
 public import Mathlib.Combinatorics.Graph.Basic
+public import Mathlib.Combinatorics.Graph.GrahamPollak
 public import Mathlib.Combinatorics.HalesJewett
 public import Mathlib.Combinatorics.Hall.Basic
 public import Mathlib.Combinatorics.Hall.Finite
diff --git a/Mathlib/Combinatorics/SimpleGraph/CompleteMultipartite.lean b/Mathlib/Combinatorics/SimpleGraph/CompleteMultipartite.lean
index 0c142464863b3f..f95b185cd4345f 100644
--- a/Mathlib/Combinatorics/SimpleGraph/CompleteMultipartite.lean
+++ b/Mathlib/Combinatorics/SimpleGraph/CompleteMultipartite.lean
@@ -57,7 +57,34 @@ open Finset Fintype
 
 universe u
 namespace SimpleGraph
-variable {α : Type u}
+variable {α : Type u} {ι : Type*}
+
+/-- `G` is `IsCompleteMultipartiteWith f` iff adjacency in `G`
+    is determined by inequality of `f`-values. -/
+def IsCompleteMultipartiteWith (G : SimpleGraph α) (f : α → ι) : Prop :=
+  ∀ u v, G.Adj u v ↔ f u ≠ f v
+
+namespace IsCompleteMultipartiteWith
+
+variable {G : SimpleGraph α} {f : α → ι} {C : G.IsCompleteMultipartiteWith f}
+include C
+
+/-- For a complete multipartite graph labelled via `f`, `u` and `v` are adjacent iff `f u ≠ f v`. -/
+def adj_iff_ne ⦃u v : α⦄ : G.Adj u v ↔ f u ≠ f v := C u v
+
+/-- The neighbors of `v` are vertices which do not match on `f`. -/
+lemma neighborSet_eq (v : α) : G.neighborSet v = {u | f v ≠ f u} := Set.ext <| C v
+
+section finite
+
+variable (v : α) [Fintype (G.neighborSet v)] [Fintype {u | f v ≠ f u}]
+
+lemma neighborFinset_eq : G.neighborFinset v = {u | f v ≠ f u}.toFinset :=
+  Set.toFinset_inj.mpr <| C.neighborSet_eq v
+
+end finite
+
+end IsCompleteMultipartiteWith
 
 /-- `G` is `IsCompleteMultipartite` iff non-adjacency is transitive -/
 def IsCompleteMultipartite (G : SimpleGraph α) : Prop := Transitive (¬ G.Adj · ·)
@@ -66,6 +93,7 @@ theorem bot_isCompleteMultipartite : (⊥ : SimpleGraph α).IsCompleteMultiparti
   simp [IsCompleteMultipartite, Transitive]
 
 variable {G : SimpleGraph α}
+
 /-- The setoid given by non-adjacency -/
 def IsCompleteMultipartite.setoid (h : G.IsCompleteMultipartite) : Setoid α :=
     ⟨(¬ G.Adj · ·), ⟨G.loopless, fun h' ↦ by rwa [adj_comm] at h', fun h1 h2 ↦ h h1 h2⟩⟩
diff --git a/Mathlib/Combinatorics/SimpleGraph/GrahamPollak.lean b/Mathlib/Combinatorics/SimpleGraph/GrahamPollak.lean
new file mode 100644
index 00000000000000..c5a217aff9e8de
--- /dev/null
+++ b/Mathlib/Combinatorics/SimpleGraph/GrahamPollak.lean
@@ -0,0 +1,149 @@
+/-
+Copyright (c) 2026 Julian Berman. All rights reserved.
+Released under Apache 2.0 license as described in the file LICENSE.
+Authors: Julian Berman, Aaron Hill
+-/
+module
+
+public import Mathlib.Analysis.RCLike.Lemmas
+public import Mathlib.Combinatorics.SimpleGraph.Bipartite
+public import Mathlib.Combinatorics.SimpleGraph.CompleteMultipartite
+public import Mathlib.Combinatorics.SimpleGraph.EdgeLabeling
+public import Mathlib.Data.Matrix.ColumnRowPartitioned
+
+/-!
+
+# The Graham-Pollak theorem
+
+-/
+
+@[expose] public section
+
+namespace SimpleGraph
+
+variable {V : Type*} {G : SimpleGraph V} {u v : V} {α ι : Type*} {a : α} {f : V → ι}
+
+def IsCompleteBipartiteWith (left : Set V) : Prop := G.IsCompleteMultipartiteWith (· ∈ left)
+
+namespace IsCompleteBipartiteWith
+
+variable {left : Set V} (C : G.IsCompleteBipartiteWith left)
+include C
+
+lemma isCompleteMultipartiteWith : G.IsCompleteMultipartiteWith (· ∈ left) := C
+
+lemma adj_iff_not_mem (hv : v ∈ left) : G.Adj v u ↔ u ∉ left := by
+  simp [C.isCompleteMultipartiteWith.adj_iff_ne, hv]
+
+lemma adj_iff_mem (hv : v ∉ left) : G.Adj v u ↔ u ∈ left := by
+  simp [C.isCompleteMultipartiteWith.adj_iff_ne, hv]
+
+lemma neighborSet_eq_of_mem_left (hv : v ∈ left) : G.neighborSet v = leftᶜ := by
+  grind [C.isCompleteMultipartiteWith.neighborSet_eq]
+
+lemma neighborSet_eq_of_not_mem_left (hv : v ∉ left) : G.neighborSet v = left := by
+  grind [C.isCompleteMultipartiteWith.neighborSet_eq]
+
+lemma bipartite : G.IsBipartiteWith left leftᶜ := by
+  refine ⟨disjoint_compl_right, fun v u hadj ↦ ?_⟩
+  grind [C.isCompleteMultipartiteWith.adj_iff_ne.mp hadj]
+
+section finite
+
+variable [Fintype ↑left] [Fintype ↑(G.neighborSet v)]
+
+lemma neighborFinset_eq_of_mem_left [Fintype ↑leftᶜ] (hv : v ∈ left.toFinset) :
+    G.neighborFinset v = leftᶜ.toFinset := by
+  grind [neighborFinset_def, neighborSet_eq_of_mem_left]
+
+lemma neighborFinset_eq_of_not_mem_left (hv : v ∉ left.toFinset) :
+    G.neighborFinset v = left.toFinset := by
+  grind [neighborFinset_def, neighborSet_eq_of_not_mem_left]
+
+end finite
+
+end IsCompleteBipartiteWith
+
+variable (G) in
+def IsCompleteBipartite :=
+  ∃ left, G.IsCompleteBipartiteWith left
+
+section finite
+
+variable [Fintype V] [Fintype α] {𝓁 : TopEdgeLabeling V α}
+
+open Finset Fintype LinearMap in
+open scoped Matrix in
+/--
+The Graham-Pollak theorem:
+
+In a complete graph on `|V|` vertices, any edge labeling into complete bipartite subgraphs uses
+at least `|V| - 1` distinct labels.
+-/
+theorem card_le_card_add_one_of_forall_IsCompleteBipartite
+  (completeBipartiteOf : ∀ a, IsCompleteBipartite <| 𝓁.labelGraph a) :
+    card V ≤ card α + 1 := by
+  classical
+  rcases subsingleton_or_nontrivial V
+  · grind [card_le_one_iff_subsingleton]
+  by_contra! h
+  let M : Matrix (Fin 1 ⊕ α) V ℝ := Matrix.fromRows
+    (Matrix.replicateCol V ![1])
+    (Matrix.of fun m n ↦ (completeBipartiteOf m).choose.indicator 1 n)
+  obtain ⟨c, hc, hc_nezero⟩ : ∃ x ∈ ker M.toLin', x ≠ 0 := (ker _).ne_bot_iff.mp <| by
+    apply ker_ne_bot_of_finrank_lt
+    simp only [Module.finrank_fintype_fun_eq_card, card_sum, card_unique]
+    grind
+  have h_sum : ∑ v, c v = 0 := by
+    have : (M *ᵥ c) (.inl 0) = 0 := by simp_all
+    simp only [M, Matrix.mulVec, dotProduct] at this
+    aesop
+  · have h_disjoint (u : V) :
+        ((univ : Finset α) : Set α).PairwiseDisjoint (𝓁.labelGraph · |>.neighborFinset u) := by
+      intro
+      grind [Finset.disjoint_left, mem_neighborFinset, 𝓁.labelGraph_adj]
+    have h_partition (u : V) :
+        univ.erase u = (univ : Finset α).biUnion (𝓁.labelGraph · |>.neighborFinset u) := by
+      ext v
+      simp only [mem_erase, Finset.mem_univ, Finset.mem_biUnion, mem_neighborFinset]
+      exact ⟨fun _ ↦ ⟨𝓁 ⟨s(u, v), by tauto⟩, by tauto⟩,
+             by grind only [TopEdgeLabeling.labelGraph_adj]⟩
+    have h_sq : ∑ v, c v ^ 2 = -∑ u, ∑ v ∈ univ.erase u, c u * c v := by simp [← mul_sum, sq, h_sum]
+    have h_zero : ∑ v, c v ^ 2 = 0 := by
+      rw [h_sq,
+          Finset.sum_congr rfl fun u _ ↦ by rw [h_partition, sum_biUnion <| h_disjoint _],
+          neg_eq_zero,
+          sum_comm]
+      refine sum_eq_zero fun v _ ↦ ?_
+      let cbp := completeBipartiteOf v
+      have h_left : ∑ u ∈ (completeBipartiteOf v).choose, c u = 0 := by
+        suffices ∑ x, ((completeBipartiteOf v).choose.toFinset : Set _).indicator c x = 0 by
+          rwa [sum_indicator_subset _ (by simp)] at this
+        have : (M *ᵥ c) (.inr v) = 0 := by simp_all
+        simp only [M, Matrix.mulVec, dotProduct, Set.indicator_apply] at this
+        aesop
+      let sum_eq (S : Finset V) := ∑ x ∈ S, ∑ i ∈ 𝓁.labelGraph v |>.neighborFinset x, c x * c i = 0
+      have h_L_sum : sum_eq cbp.choose.toFinset := by
+        dsimp [sum_eq]
+        rw [Finset.sum_congr rfl fun _ hv ↦ by
+          rw [cbp.choose_spec.neighborFinset_eq_of_mem_left hv, ← mul_sum],
+            ← sum_mul, h_left, zero_mul]
+      have h_R_sum : sum_eq cbp.chooseᶜ.toFinset := by
+        dsimp [sum_eq]
+        rw [Finset.sum_congr rfl fun _ hx ↦ by
+          rw [cbp.choose_spec.neighborFinset_eq_of_not_mem_left (by grind), ← mul_sum],
+            h_left]
+        simp only [mul_zero, sum_const_zero]
+      rw [show univ = cbp.choose.toFinset ∪ cbp.chooseᶜ.toFinset by simp,
+          sum_union <| Finset.disjoint_left.mpr fun v hvL hvR ↦
+            Set.disjoint_left.mp cbp.choose_spec.bipartite.disjoint
+              (Set.mem_toFinset.mp hvL) (Set.mem_toFinset.mp hvR),
+          h_L_sum,
+          h_R_sum,
+          add_zero]
+    simp only [Finset.sum_eq_zero_iff_of_nonneg (fun _ _ ↦ sq_nonneg _), sq_eq_zero_iff] at h_zero
+    bound
+
+end finite
+
+end SimpleGraph