Huang's sensitivity theorem

A formalization of Hao Huang's sensitivity theorem: in the hypercube of dimension n ≥ 1, if one colors more than half the vertices then at least one vertex has at least √n colored neighbors.

A fun summer collaboration by Reid Barton, Johan Commelin, Jesse Han, Chris Hughes, Robert Y. Lewis, and Patrick Massot, based on Don Knuth's account of the story (https://www.cs.stanford.edu/~knuth/papers/huang.pdf), using the Lean theorem prover (https://leanprover.github.io/), by Leonardo de Moura at Microsoft Research, and his collaborators (https://leanprover.github.io/people/), and using Lean's user maintained mathematics library (https://github.com/leanprover-community/mathlib).

The project was developed at https://github.com/leanprover-community/lean-sensitivity and is now archived at https://github.com/leanprover-community/mathlib/blob/master/archive/sensitivity.lean

The next two lines assert we do not want to give a constructive proof, but rather use classical logic.

We also want to use the notation ∑ for sums.

The hypercube

Notations:

• ℕ denotes natural numbers (including zero).
• Fin n = {0, ⋯ , n - 1}.
• Bool = {true, false}.

def Sensitivity.Q (n : ℕ) : Type

The hypercube in dimension n.

Equations

• Sensitivity.Q n = (Fin n → Bool)

instance Sensitivity.instInhabitedQ (n : ℕ) : Inhabited (Sensitivity.Q n)

Equations

• Sensitivity.instInhabitedQ n = inferInstanceAs (Inhabited (Fin n → Bool))

instance Sensitivity.instFintypeQ (n : ℕ) : Fintype (Sensitivity.Q n)

Equations

• Sensitivity.instFintypeQ n = inferInstanceAs (Fintype (Fin n → Bool))

def Sensitivity.π {n : ℕ} : Sensitivity.Q n.succ → Sensitivity.Q n

The projection from Q n.succ to Q n forgetting the first value (ie. the image of zero).

Equations

• Sensitivity.π p = p ∘ Fin.succ

theorem Sensitivity.Q.ext {n : ℕ} {f : Sensitivity.Q n} {g : Sensitivity.Q n} (h : ∀ (x : Fin n), f x = g x) : f = g

n will always denote a natural number.

instance Sensitivity.Q.instUniqueOfNatNat : Unique (Sensitivity.Q 0)

Q 0 has a unique element.

Equations

• Sensitivity.Q.instUniqueOfNatNat = { default := fun (x : Fin 0) => true, uniq := Sensitivity.Q.instUniqueOfNatNat.proof_1 }

theorem Sensitivity.Q.card (n : ℕ) : Fintype.card (Sensitivity.Q n) = 2 ^ n

Q n has 2^n elements.

Until the end of this namespace, n will be an implicit argument (still a natural number).

theorem Sensitivity.Q.succ_n_eq {n : ℕ} (p : Sensitivity.Q n.succ) (q : Sensitivity.Q n.succ) : p = q ↔ p 0 = q 0 ∧ Sensitivity.π p = Sensitivity.π q

def Sensitivity.Q.adjacent {n : ℕ} (p : Sensitivity.Q n) : Set (Sensitivity.Q n)

The adjacency relation defining the graph structure on Q n: p.adjacent q if there is an edge from p to q in Q n.

Equations

• p.adjacent = {q : Sensitivity.Q n | ∃! i : Fin n, p i ≠ q i}

theorem Sensitivity.Q.not_adjacent_zero (p : Sensitivity.Q 0) (q : Sensitivity.Q 0) : q ∉ p.adjacent

In Q 0, no two vertices are adjacent.

theorem Sensitivity.Q.adj_iff_proj_eq {n : ℕ} {p : Sensitivity.Q n.succ} {q : Sensitivity.Q n.succ} (h₀ : p 0 ≠ q 0) : q ∈ p.adjacent ↔ Sensitivity.π p = Sensitivity.π q

If p and q in Q n.succ have different values at zero then they are adjacent iff their projections to Q n are equal.

theorem Sensitivity.Q.adj_iff_proj_adj {n : ℕ} {p : Sensitivity.Q n.succ} {q : Sensitivity.Q n.succ} (h₀ : p 0 = q 0) : q ∈ p.adjacent ↔ Sensitivity.π q ∈ (Sensitivity.π p).adjacent

If p and q in Q n.succ have the same value at zero then they are adjacent iff their projections to Q n are adjacent.

theorem Sensitivity.Q.adjacent.symm {n : ℕ} {p : Sensitivity.Q n} {q : Sensitivity.Q n} : q ∈ p.adjacent ↔ p ∈ q.adjacent

The vector space

def Sensitivity.V : ℕ → Type

The free vector space on vertices of a hypercube, defined inductively.

Equations

• Sensitivity.V 0 = ℝ
• Sensitivity.V n.succ = (Sensitivity.V n × Sensitivity.V n)

theorem Sensitivity.V.ext {n : ℕ} {p : Sensitivity.V n.succ} {q : Sensitivity.V n.succ} (h₁ : p.1 = q.1) (h₂ : p.2 = q.2) : p = q

V n is a real vector space whose equality relation is computable.

instance Sensitivity.V.instDecidableEq (n : ℕ) : DecidableEq (Sensitivity.V n)

Equations

• Sensitivity.V.instDecidableEq n = Nat.recAux (id inferInstance) (fun (n : ℕ) (a : DecidableEq (Sensitivity.V n)) => id inferInstance) n

instance Sensitivity.V.instAddCommGroup (n : ℕ) : AddCommGroup (Sensitivity.V n)

Equations

• Sensitivity.V.instAddCommGroup n = Nat.recAux (id inferInstance) (fun (n : ℕ) (a : AddCommGroup (Sensitivity.V n)) => id inferInstance) n

instance Sensitivity.V.instModuleReal (n : ℕ) : Module ℝ (Sensitivity.V n)

Equations

• Sensitivity.V.instModuleReal n = Nat.recAux (id inferInstance) (fun (n : ℕ) (a : Module ℝ (Sensitivity.V n)) => id inferInstance) n

noncomputable def Sensitivity.e {n : ℕ} : Sensitivity.Q n → Sensitivity.V n

The basis of V indexed by the hypercube, defined inductively.

Equations

• Sensitivity.e = fun (x : Sensitivity.Q 0) => 1
• Sensitivity.e = fun (x : Sensitivity.Q n.succ) => bif x 0 then (Sensitivity.e (Sensitivity.π x), 0) else (0, Sensitivity.e (Sensitivity.π x))

@[simp]

theorem Sensitivity.e_zero_apply (x : Sensitivity.Q 0) : Sensitivity.e x = 1

noncomputable def Sensitivity.ε {n : ℕ} : Sensitivity.Q n → Sensitivity.V n →ₗ[ℝ] ℝ

The dual basis to e, defined inductively.

Equations

• One or more equations did not get rendered due to their size.
• Sensitivity.ε x_2 = LinearMap.id

theorem Sensitivity.duality {n : ℕ} (p : Sensitivity.Q n) (q : Sensitivity.Q n) : (Sensitivity.ε p) (Sensitivity.e q) = if p = q then 1 else 0

theorem Sensitivity.epsilon_total {n : ℕ} {v : Sensitivity.V n} (h : ∀ (p : Sensitivity.Q n), (Sensitivity.ε p) v = 0) : v = 0

Any vector in V n annihilated by all ε p's is zero.

theorem Sensitivity.dualBases_e_ε (n : ℕ) : Module.DualBases Sensitivity.e Sensitivity.ε

e and ε are dual families of vectors. It implies that e is indeed a basis and ε computes coefficients of decompositions of vectors on that basis.

We will now derive the dimension of V, first as a cardinal in dim_V and, since this cardinal is finite, as a natural number in finrank_V

theorem Sensitivity.dim_V {n : ℕ} : Module.rank ℝ (Sensitivity.V n) = 2 ^ n

instance Sensitivity.instFiniteDimensionalRealV {n : ℕ} : FiniteDimensional ℝ (Sensitivity.V n)

Equations

• ⋯ = ⋯

theorem Sensitivity.finrank_V {n : ℕ} : FiniteDimensional.finrank ℝ (Sensitivity.V n) = 2 ^ n

The linear map

noncomputable def Sensitivity.f (n : ℕ) : Sensitivity.V n →ₗ[ℝ] Sensitivity.V n

The linear operator $f_n$ corresponding to Huang's matrix $A_n$, defined inductively as a ℝ-linear map from V n to V n.

Equations

• Sensitivity.f 0 = 0
• Sensitivity.f n.succ = ((Sensitivity.f n).coprod LinearMap.id).prod (LinearMap.id.coprod (-Sensitivity.f n))

The preceding definition uses linear map constructions to automatically get that f is linear, but its values are somewhat buried as a side-effect. The next two lemmas unbury them.

@[simp]

theorem Sensitivity.f_zero : Sensitivity.f 0 = 0

theorem Sensitivity.f_succ_apply {n : ℕ} (v : Sensitivity.V n.succ) : (Sensitivity.f n.succ) v = ((Sensitivity.f n) v.1 + v.2, v.1 - (Sensitivity.f n) v.2)

In the next statement, the explicit conversion (n : ℝ) of n to a real number is necessary since otherwise n • v refers to the multiplication defined using only the addition of V.

theorem Sensitivity.f_squared {n : ℕ} (v : Sensitivity.V n) : (Sensitivity.f n) ((Sensitivity.f n) v) = ↑n • v

We now compute the matrix of f in the e basis (p is the line index, q the column index).

theorem Sensitivity.f_matrix {n : ℕ} (p : Sensitivity.Q n) (q : Sensitivity.Q n) : |(Sensitivity.ε q) ((Sensitivity.f n) (Sensitivity.e p))| = if p ∈ q.adjacent then 1 else 0

noncomputable def Sensitivity.g (m : ℕ) : Sensitivity.V m →ₗ[ℝ] Sensitivity.V m.succ

The linear operator $g_m$ corresponding to Knuth's matrix $B_m$.

Equations

• Sensitivity.g m = (Sensitivity.f m + (↑m + 1).sqrt • LinearMap.id).prod LinearMap.id

In the following lemmas, m will denote a natural number.

Again we unpack what are the values of g.

theorem Sensitivity.g_apply {m : ℕ} (v : Sensitivity.V m) : (Sensitivity.g m) v = ((Sensitivity.f m) v + (↑m + 1).sqrt • v, v)

theorem Sensitivity.g_injective {m : ℕ} : Function.Injective ⇑(Sensitivity.g m)

theorem Sensitivity.f_image_g {m : ℕ} (w : Sensitivity.V m.succ) (hv : ∃ (v : Sensitivity.V m), (Sensitivity.g m) v = w) : (Sensitivity.f m.succ) w = (↑m + 1).sqrt • w

The main proof

In this section, in order to enforce that n is positive, we write it as m + 1 for some natural number m.

dim X will denote the dimension of a subspace X as a cardinal.

fdim X will denote the (finite) dimension of a subspace X as a natural number.

Span S will denote the ℝ-subspace spanned by S.

Card X will denote the cardinal of a subset of a finite type, as a natural number.

In the following, ⊓ and ⊔ will denote intersection and sums of ℝ-subspaces, equipped with their subspace structures. The notations come from the general theory of lattices, with inf and sup (also known as meet and join).

theorem Sensitivity.exists_eigenvalue {m : ℕ} (H : Set (Sensitivity.Q m.succ)) (hH : H.toFinset.card ≥ 2 ^ m + 1) : ∃ y ∈ Submodule.span ℝ (Sensitivity.e '' H) ⊓ LinearMap.range (Sensitivity.g m), y ≠ 0

If a subset H of Q (m+1) has cardinal at least 2^m + 1 then the subspace of V (m+1) spanned by the corresponding basis vectors non-trivially intersects the range of g m.

theorem Sensitivity.huang_degree_theorem {m : ℕ} (H : Set (Sensitivity.Q m.succ)) (hH : H.toFinset.card ≥ 2 ^ m + 1) : ∃ q ∈ H, (↑m + 1).sqrt ≤ ↑(H ∩ q.adjacent).toFinset.card

Huang sensitivity theorem also known as the Huang degree theorem