Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

Graph Theory

322 missions · 117 completed

Missions

Open205Completed117All322
Combinatorics·Captain: Gabewhigham

Conway's 99-graph problemOpen Problem

Motivation

A strongly regular graph with parameters (n,k,λ,μ)(n,k,\lambda,\mu)(n,k,λ,μ) is a finite simple graph on nnn vertices in which every vertex has exactly kkk neighbours, every pair of adjacent vertices has exactly λ\lambdaλ common neighbours, and every pair of non-adjacent vertices has exactly μ\muμ common neighbours. For most parameter tuples the elementary counting and integrality conditions already decide existence; the interesting cases are those that survive every known feasibility test and still resist construction. The tuple (99,14,1,2)(99,14,1,2)(99,14,1,2) is the smallest such case in the family λ=1\lambda = 1λ=1, μ=2\mu = 2μ=2, and its existence has been open for more than fifty years. John Horton Conway offered $1000 for a resolution, as one of five problems posed at the 2014 DIMACS conference on Challenges of Identifying Integer Sequences (Conway, Five $1,000 Problems (Update 2017)).

Timeline of the problem and of what is known about it:

  • 1969/1971 — the parameter set is raised by Norman Biggs in his Southampton lectures (Finite Groups of Automorphisms, LMS Lecture Note Series 6, p. 111).
  • 1973 — Berlekamp, van Lint and Seidel construct a strongly regular graph with parameters (243,22,1,2)(243,22,1,2)(243,22,1,2) as the coset graph of the perfect ternary Golay code, settling one of the five feasible parameter tuples in this family.
  • 1975 — the existence question appears as Problem 7 (attributed to J. J. Seidel) in R. K. Guy's problem list, The Geometry of Metric and Linear Spaces, Springer LNM 490, pp. 237–238; Conway had worked on it by then.
  • 1984 — H. A. Wilbrink, On the (99,14,1,2)(99,14,1,2)(99,14,1,2) strongly regular graph, shows that such a graph cannot be vertex-transitive: no group of automorphisms can act transitively on its 99 vertices.
  • 1988 — Brouwer and Neumaier, A remark on partial linear spaces of girth 5 with an application to strongly regular graphs, Combinatorica 8, 57–61.
  • 2004 — Makhnev and Minakova, On automorphisms of strongly regular graphs with λ=1\lambda=1λ=1, μ=2\mu=2μ=2, Discrete Math. Appl. 14(2), and 2011 — Behbahani and Lam, Strongly regular graphs with non-trivial automorphisms, Discrete Math. 311, 132–144: further restrictions on the possible automorphism groups.
  • 2014/2017 — Conway's prize offer publicises the problem.

No graph with these parameters has been found, and no non-existence proof is known.

Setting

Fix a finite vertex set VVV and a simple graph ggg on VVV (irreflexive, symmetric adjacency Adj\mathrm{Adj}Adj). For vertices v,wv,wv,w write N(v)={u:Adj(v,u)}N(v) = \{u : \mathrm{Adj}(v,u)\}N(v)={u:Adj(v,u)} for the neighbourhood of vvv and N(v)∩N(w)N(v)\cap N(w)N(v)∩N(w) for the set of common neighbours. The graph ggg is strongly regular with parameters (n,k,λ,μ)(n,k,\lambda,\mu)(n,k,λ,μ), written IsSRGWith g n k λ μ\mathrm{IsSRGWith}\ g\ n\ k\ \lambda\ \muIsSRGWith g n k λ μ, when

  • ∣V∣=n|V| = n∣V∣=n;
  • ∣N(v)∣=k|N(v)| = k∣N(v)∣=k for every vertex vvv;
  • ∣N(v)∩N(w)∣=λ|N(v)\cap N(w)| = \lambda∣N(v)∩N(w)∣=λ whenever vvv and www are adjacent;
  • ∣N(v)∩N(w)∣=μ|N(v)\cap N(w)| = \mu∣N(v)∩N(w)∣=μ whenever v≠wv \neq wv=w are non-adjacent.

The case λ=1\lambda = 1λ=1 says that every edge lies in exactly one triangle — equivalently, the neighbourhood of each vertex induces a perfect matching, so such graphs are locally linear. The case μ=2\mu = 2μ=2 says that every non-adjacent pair is the pair of opposite corners of exactly one 444-cycle. Conway's problem asks for (n,k)=(99,14)(n,k) = (99,14)(n,k)=(99,14) with these two local conditions.

Counting paths of length two from a fixed vertex gives k(k−λ−1)=(n−k−1)μk(k-\lambda-1) = (n-k-1)\muk(k−λ−1)=(n−k−1)μ, which for λ=1\lambda=1λ=1, μ=2\mu=2μ=2 reduces to 2n=k2+22n = k^2 + 22n=k2+2; with k=14k = 14k=14 this yields n=99n = 99n=99. Writing AAA for the adjacency matrix, III for the identity and JJJ for the all-ones matrix, strong regularity is equivalent to the matrix identity A2=kI+λA+μ(J−I−A)A^2 = kI + \lambda A + \mu(J - I - A)A2=kI+λA+μ(J−I−A), which for (99,14,1,2)(99,14,1,2)(99,14,1,2) reads A2+A=12I+2JA^2 + A = 12I + 2JA2+A=12I+2J; the eigenvalues of AAA other than k=14k=14k=14 are then 333 and −4-4−4, and integrality of their multiplicities (545454 and 444444) is one of the feasibility conditions that (99,14,1,2)(99,14,1,2)(99,14,1,2) passes.

Formalization targets

Goal

∃ α, ∃ g a simple graph on α,IsSRGWith g 99 14 1 2.\exists\ \alpha,\ \exists\ g \text{ a simple graph on } \alpha,\quad \mathrm{IsSRGWith}\ g\ 99\ 14\ 1\ 2 .∃ α, ∃ g a simple graph on α,IsSRGWith g 99 14 1 2.

The goal is Mathlib's own proof_wanted conway_99 in Mathlib/Combinatorics/SimpleGraph/StronglyRegular.lean, stated verbatim: existence of a finite type carrying a strongly regular graph with parameters (99,14,1,2)(99,14,1,2)(99,14,1,2). A resolution in either direction is welcome — a proof settles the existence half, and a proof of the negation settles the non-existence half; the platform records the two as proof and disproof of the same statement.

Supporting targets

2n=k2+2,k even,k∈{2,4,14,22,112,994}2n = k^2 + 2, \qquad k \text{ even}, \qquad k \in \{2,4,14,22,112,994\}2n=k2+2,k even,k∈{2,4,14,22,112,994}

for every strongly regular graph with λ=1\lambda = 1λ=1, μ=2\mu = 2μ=2: the counting identity, local linearity, and the integrality restriction that cuts the family down to five non-degenerate parameter tuples.

∃ g, IsSRGWith g 9 4 1 2,∃ g, IsSRGWith g 243 22 1 2\exists\, g,\ \mathrm{IsSRGWith}\ g\ 9\ 4\ 1\ 2, \qquad \exists\, g,\ \mathrm{IsSRGWith}\ g\ 243\ 22\ 1\ 2∃g, IsSRGWith g 9 4 1 2,∃g, IsSRGWith g 243 22 1 2

the two members of the family that are known to exist: the 3×33\times 33×3 rook's graph (the Paley graph on 999 vertices) and the Berlekamp–van Lint–Seidel graph.

∣E(g)∣=693,∣{triangles of g}∣=231,A2+A=12I+2J,g not vertex-transitive|E(g)| = 693, \qquad |\{\text{triangles of } g\}| = 231, \qquad A^2 + A = 12I + 2J, \qquad g \text{ not vertex-transitive}∣E(g)∣=693,∣{triangles of g}∣=231,A2+A=12I+2J,g not vertex-transitive

structural consequences for a hypothetical 999999-graph, the last one being Wilbrink's theorem.

Significance

A (99,14,1,2)(99,14,1,2)(99,14,1,2) graph, if it exists, is a locally linear graph of maximal density in its parameter range and a partial linear space of girth 555 with 999999 points and 231231231 lines of size 333; its existence would also produce new association schemes and new examples for the general classification of strongly regular graphs. A non-existence proof would be the first case in this family ruled out by anything other than the classical feasibility conditions, and would say something new about how far local conditions (λ=1\lambda=1λ=1, μ=2\mu=2μ=2) constrain global structure.

Nothing in this mission is presently formalized. Mathlib defines SimpleGraph.IsSRGWith, proves the counting identity IsSRGWith.param_eq, the complement rule IsSRGWith.compl, and the matrix identity IsSRGWith.matrix_eq, and records the 999999-graph problem as a proof_wanted. The supporting targets are of three kinds: results that are proved in the literature and only need formalizing (existence at (9,4,1,2)(9,4,1,2)(9,4,1,2) and (243,22,1,2)(243,22,1,2)(243,22,1,2); Wilbrink's non-vertex-transitivity; the integrality restriction on kkk); routine consequences that supply reusable infrastructure (edge and triangle counts, the spectral identity, evenness of kkk); and the goal itself, which is open mathematics.

Difficulty

The obvious approaches fail for concrete reasons. Exhaustive search is out of range: the graph has 693693693 edges among (992)=4851\binom{99}{2} = 4851(299​)=4851 pairs, and no isomorph-free generation of locally linear graphs on 999999 vertices is feasible. Algebraic constructions are blocked by Wilbrink's theorem — the graph cannot be vertex-transitive, so it is not a Cayley graph and cannot be produced by the group-theoretic constructions that yield most known strongly regular graphs, including the two that work at (9,4,1,2)(9,4,1,2)(9,4,1,2) and (243,22,1,2)(243,22,1,2)(243,22,1,2). On the non-existence side, every classical feasibility test (the counting identity, integrality of the eigenvalue multiplicities, the Krein conditions, the absolute bound) is passed by (99,14,1,2)(99,14,1,2)(99,14,1,2), so a proof of non-existence needs an argument that does not factor through the parameters alone.

Formalization scope

All statements are phrased with Mathlib's SimpleGraph.IsSRGWith on a Fintype vertex type with DecidableRel adjacency, and use Fintype.card, SimpleGraph.edgeFinset, SimpleGraph.cliqueFinset 3 (triangles as 333-cliques), SimpleGraph.adjMatrix over Z\mathbb{Z}Z, and graph isomorphisms g ≃g g for automorphisms. The goal quantifies over α : Type together with a Fintype α instance, so the vertex set is finite by construction and the empty type does not satisfy the cardinality clause; the statement is therefore not vacuously satisfiable. Note that Mathlib's definition constrains λ\lambdaλ only through pairs that are actually adjacent and μ\muμ only through pairs that are actually distinct and non-adjacent, so degenerate small graphs (the one-vertex graph, K3K_3K3​) do satisfy IsSRGWith with λ=1\lambda = 1λ=1, μ=2\mu = 2μ=2; the supporting statements carry the cardinality hypotheses (0<n0 < n0<n, 1<n1 < n1<n) that exclude them where needed, and the degenerate degree k=2k = 2k=2 is listed explicitly in the classification of feasible degrees.

Infrastructure a complete development needs, and which is reusable beyond this mission: interface lemmas for counting common neighbours in a strongly regular graph; the spectral theory of the adjacency matrix (multiplicities of the two non-principal eigenvalues, and their integrality), which is the missing ingredient for the classification of feasible degrees; a Lean construction of the perfect ternary Golay code and its coset graph, for the (243,22,1,2)(243,22,1,2)(243,22,1,2) case; and decision procedures for strong regularity of an explicitly given small graph, for the (9,4,1,2)(9,4,1,2)(9,4,1,2) case. Contributions to any of these are welcome, as are partial non-existence results (for instance, restrictions on automorphisms of prime order) submitted as separate statements.

Selected references

  • N. Biggs, Finite Groups of Automorphisms: Course Given at the University of Southampton, October–December 1969, London Mathematical Society Lecture Note Series 6, Cambridge University Press, 1971, p. 111.
  • E. R. Berlekamp, J. H. van Lint, J. J. Seidel, A strongly regular graph derived from the perfect ternary Golay code, in: A Survey of Combinatorial Theory, North-Holland, 1973, pp. 25–30.
  • R. K. Guy, Problems, in: The Geometry of Metric and Linear Spaces, Springer Lecture Notes in Mathematics 490, 1975, pp. 233–244 (Problem 7, J. J. Seidel, pp. 237–238). doi:10.1007/BFb0081147
  • H. A. Wilbrink, On the (99,14,1,2)(99,14,1,2)(99,14,1,2) strongly regular graph, in: Papers dedicated to J. J. Seidel, EUT Report 84-WSK-03, Eindhoven University of Technology, 1984, pp. 342–355. PDF
  • A. E. Brouwer, A. Neumaier, A remark on partial linear spaces of girth 5 with an application to strongly regular graphs, Combinatorica 8 (1988), 57–61. doi:10.1007/BF02122552
  • A. A. Makhnev, I. M. Minakova, On automorphisms of strongly regular graphs with λ=1\lambda=1λ=1, μ=2\mu=2μ=2, Discrete Mathematics and Applications 14 (2004), no. 2. doi:10.1515/156939204872374
  • M. Behbahani, C. Lam, Strongly regular graphs with non-trivial automorphisms, Discrete Mathematics 311 (2011), 132–144. doi:10.1016/j.disc.2010.10.005
  • J. H. Conway, Five $1,000 Problems (Update 2017), OEIS. PDF
24 thms7 active usersReviewed
Combinatorics·Captain: hao jia

P3-Partitions of Cubic 3-Connected Graphs (OPG-46613)Open Problem

Motivation

A P3P_3P3​-packing in a graph is a collection of pairwise vertex-disjoint paths on three vertices. Determining the largest such packing is NP-hard even in restricted graph classes, so structural hypotheses that force an optimal packing are of independent interest in graph factor theory. The present question asks whether 3-vertex-connectivity and cubicity force the strongest possible packing whenever the vertex count permits a perfect partition.

A. Kelmans attributes the broader packing problem to 1984. In Problem 1.10 of Packing 3-vertex Paths in Cubic 3-connected Graphs, the question is whether every cubic 3-connected graph GGG satisfies λ(G)=⌊∣V(G)∣/3⌋\lambda(G)=\lfloor |V(G)|/3\rfloorλ(G)=⌊∣V(G)∣/3⌋. Theorem 3.1 of that paper proves that the divisible-order factor statement is equivalent to several apparently stronger deletion and prescribed-edge statements; it does not prove the open claim itself. OPG-46613 records the divisible-order form targeted here.

A 2026 candidate analysis in the Vibe Mathing problem repository investigated a tempting sufficient route: find a perfect matching whose complementary 2-factor has every cycle length divisible by three. Candidate C01 explains why that condition would yield a P3P_3P3​-factor. Candidate C02 gives an explicit proposed family HqH_qHq​ of order 18+12q18+12q18+12q that has P3P_3P3​-factors but is claimed not to satisfy the stronger matching condition. These candidate claims have computational and partial Lean checks, but no complete Lean kernel proof; they are milestones here, not declarations that the original problem or the candidate family has already been formally established.

Setting

All graphs are finite and simple. A graph is cubic when every vertex has exactly three neighbors. It is 3-vertex-connected here when it has at least four vertices and deleting any set of at most two vertices leaves a connected induced graph.

A P3P_3P3​-factor is represented by a natural number bbb, together with a bijection

Fin⁡(b)×Fin⁡(3)≃V(G),\operatorname{Fin}(b)\times\operatorname{Fin}(3)\simeq V(G),Fin(b)×Fin(3)≃V(G),

such that, in every block, positions 000 and 111 are adjacent and positions 111 and 222 are adjacent. The path is not required to be induced: an ambient edge between positions 000 and 222 is allowed because the two selected path edges still form a copy of P3P_3P3​.

A 2-factor is a spanning 2-regular subgraph. It is called divisible when every one of its connected components has order divisible by three. A divisible matching complement is a perfect matching MMM such that the relative complement G∖MG\setminus MG∖M is a divisible 2-factor.

The explicit graph HqH_qHq​ is defined on Fin⁡(18+12q)\operatorname{Fin}(18+12q)Fin(18+12q). Its first nine vertices form the fixed Petersen-minus-one-vertex brick from C02; the remaining vertices form the stated cycle-and-opposite-chord brick with three joining edges. The full adjacency relation is part of the Lean definition rather than an external data file.

Formalization targets

Main goal

For every finite simple graph GGG,

(G cubic)∧(G 3-vertex-connected)∧3∣∣V(G)∣⟹G has a P3-factor.\bigl(G\text{ cubic}\bigr)\land \bigl(G\text{ 3-vertex-connected}\bigr)\land 3\mid |V(G)| \quad\Longrightarrow\quad G\text{ has a }P_3\text{-factor}.(G cubic)∧(G 3-vertex-connected)∧3∣∣V(G)∣⟹G has a P3​-factor.

This is the OPG-46613 target. Cubicity forces the order to be even, so within this domain divisibility by three is equivalent to divisibility by six.

Literature and route milestones

The mission also formalizes the (z1)⇔(z8)(z1)\Leftrightarrow(z8)(z1)⇔(z8) part of Kelmans's Theorem 3.1: the divisible-order factor claim is equivalent to the assertion that deleting any specified 3-vertex path leaves a P3P_3P3​-factor. Two route lemmas state that divisible 2-factors split into P3P_3P3​-factors and that, in cubic graphs, divisible 2-factors are equivalent to divisible perfect-matching complements.

Candidate boundary milestones

The C02 milestones ask first for the complete 18-vertex statement and then for the full family:

∀q∈N,Hq is cubic and 3-vertex-connected, has a P3-factor, and has no divisible matching complement.\forall q\in\mathbb N,\quad H_q\text{ is cubic and 3-vertex-connected, has a }P_3\text{-factor, and has no divisible matching complement}.∀q∈N,Hq​ is cubic and 3-vertex-connected, has a P3​-factor, and has no divisible matching complement.

This separates a sufficient method from the root conclusion. It is not a counterexample to OPG-46613 because every HqH_qHq​ in the proposed family explicitly satisfies the desired P3P_3P3​ conclusion.

Significance

A proof of the main goal would settle the divisible-order form of a long-standing path-packing problem. Through Kelmans's equivalences it would also control several deletion and prescribed-edge variants for cubic 3-connected graphs. A disproof would require a graph satisfying all domain hypotheses but lacking a P3P_3P3​-factor; the C02 family does not claim this.

Formalizing the candidate boundary is useful even before the root is resolved. It turns a route exclusion into a checkable theorem and prevents a search campaign from silently assuming that every relevant graph possesses a divisible complementary 2-factor. The definitions of noninduced P3P_3P3​-factors, vertex connectivity by deletion, perfect matchings, 2-factors, and component-order divisibility are intended to be reusable in later graph-factor work.

Difficulty

The perfect-matching route is attractive because the complement of a perfect matching in a cubic graph is 2-regular. The obstruction is that its cycles need not have lengths divisible by three. The C02 candidate family is designed to expose exactly that gap: a persistent 5-cycle is claimed to occur in every complementary 2-factor even though an unrelated P3P_3P3​-factor exists. Consequently, proving the main theorem cannot simply assume that a favorable perfect matching always exists.

The formal difficulty is also semantic. Connectivity must mean vertex connectivity, the complement must be relative to GGG on the same vertex set, component sizes must refer to the 2-factor rather than the ambient graph, and P3P_3P3​ must remain noninduced. Weakening any of these points can create a materially different or vacuous theorem.

Formalization scope

The development targets Lean 4.33.1 and Mathlib revision 0df444a360eaa60ab8c11dca51a86af692955474. Graphs use SimpleGraph on finite vertex types. Degree is the cardinality of the actual neighbor subtype. Three-vertex-connectivity explicitly quantifies over all finite deletion sets of cardinality at most two and includes a four-vertex order guard.

The main theorem is universe-polymorphic and does not hard-code a finite graph enumeration. The HqH_qHq​ family includes q=0q=0q=0. The factor structure uses a bijection, so disjointness and coverage cannot be discharged by duplicate or omitted vertices. Ambient chords do not invalidate a block, while both required consecutive adjacencies must be genuine graph edges. The candidate family statements remain open theorem goals ending in sorry; the shared definition module itself is sorry-free.

Welcome contributions include proofs of the model lemmas, the finite H0H_0H0​ statement, the general C02 family, Kelmans's equivalence, or decompositions of the root theorem into faithful reusable lemmas. Numerical enumeration alone is supporting evidence and should not be presented as a kernel proof.

Selected references

  • A. Kelmans, Packing 3-vertex Paths In Cubic 3-connected Graphs, arXiv:0910.2766v2, 2011, Problem 1.10 (p. 3) and Theorem 3.1 (pp. 7–8). https://arxiv.org/abs/0910.2766v2
  • UnsolvedMath, OPG-46613: P3-partitions of cubic 3-connected graphs. https://www.unsolvedmath.com/problems/OPG-46613
  • Vibe Mathing, C01: divisible-cycle implication and a 30-vertex obstruction, fixed repository revision 14b8dc64ac2d89c98cf3a2bbb2fcba76ced0df6a. https://github.com/vibemathing/problem-opg-46613-cubic-p3-partition/blob/14b8dc64ac2d89c98cf3a2bbb2fcba76ced0df6a/research/artifacts/candidates/opg46613-c01/proof.md
  • Vibe Mathing, C02: an 18-vertex obstruction and an infinite family with P3-factors, fixed repository revision 14b8dc64ac2d89c98cf3a2bbb2fcba76ced0df6a. https://github.com/vibemathing/problem-opg-46613-cubic-p3-partition/blob/14b8dc64ac2d89c98cf3a2bbb2fcba76ced0df6a/research/artifacts/candidates/opg46613-c02/proof.md
16 thms5 active usersReviewed
Computational GeometryDiscrete Geometry·Captain: hao jia

Uniform Obstacle Bounds for Planar Graphs (OPG-37357)Open Problem

Motivation

An obstacle representation turns a graph into a visibility system: vertices are points in the plane, and nonedges are blocked by polygonal obstacles. The obstacle number asks for the minimum number of obstacles needed. OPG-37357 records two different questions for planar graphs. The first asks whether one obstacle can ever be insufficient. The second asks whether some universal constant bounds the ordinary obstacle number of every planar graph.

The status of the two parts is different. Berman, Chappell, Faudree, Gimbel, Hartman, and Williams proved in 2017 that explicit planar graphs, including the icosahedron and their graphs X4X_4X4​ and X6X_6X6​, have ordinary obstacle number two. Thus the first question has a published positive answer. The universal-constant question remains the research target here. A separate invariant called planar or plane obstacle number requires a crossing-free visibility drawing; results for that invariant must not be substituted for the ordinary obstacle number used by this mission.

Setting

A finite simple graph GGG has a kkk-obstacle drawing when its vertices are placed injectively as points in R2\mathbb R^2R2 and there are kkk pairwise disjoint closed connected polygonal obstacles such that

uv∈E(G)⟺[p(u),p(v)] meets no obstacle.uv\in E(G) \quad\Longleftrightarrow\quad [p(u),p(v)]\text{ meets no obstacle}.uv∈E(G)⟺[p(u),p(v)] meets no obstacle.

Graph vertices lie outside every obstacle. The ordinary obstacle number obs⁡(G)\operatorname{obs}(G)obs(G) is the least such kkk. The drawing itself may contain crossings between visible graph edges; planarity is a property of the abstract input graph, not an extra constraint on the obstacle drawing.

The Lean model represents a polygonal obstacle as a connected finite union of closed filled triangles. This gives a compact polygonal region with exact real-coordinate segment incidence. Straight-line planarity of the abstract graph is represented separately.

Formalization targets

The two-part OPG record

The source records both

∃ finite planar G, obs⁡(G)>1\exists\text{ finite planar }G,\ \operatorname{obs}(G)>1∃ finite planar G, obs(G)>1

and

∃k∈N ∀ finite planar H, obs⁡(H)≤k.\exists k\in\mathbb N\ \forall\text{ finite planar }H, \ \operatorname{obs}(H)\le k.∃k∈N ∀ finite planar H, obs(H)≤k.

The first assertion is known in the literature and appears as a published-result milestone. The second is open and is therefore the mission's main theorem. Together they preserve the two-part source without presenting the whole record as unresolved.

Published first part

A milestone formalizes the stronger published statement

∃ finite planar G,obs⁡(G)≤2andobs⁡(G)≰1.\exists\text{ finite planar }G, \qquad \operatorname{obs}(G)\le2 \quad\text{and}\quad \operatorname{obs}(G)\not\le1.∃ finite planar G,obs(G)≤2andobs(G)≤1.

This captures ordinary obstacle number exactly two without hard-coding one graph before its adjacency data and lower-bound certificate are formalized.

Universal bound

The open milestone asks for a single natural number kkk, chosen before the graph, that works for every finite planar graph. The number of obstacle corners is not bounded by this theorem; only the number of connected polygonal obstacles is.

Significance

The published first part establishes that planarity alone does not force a one-obstacle representation. The second part asks whether planar graphs nevertheless have uniformly bounded visibility complexity. A positive answer would produce a common finite obstacle budget independent of graph order; a negative answer would require a family of planar graphs with unbounded ordinary obstacle number.

Formalization is especially useful because several nearby notions differ by one word but have different known bounds: ordinary versus plane obstacle number, arbitrary polygonal versus convex obstacles, and fixed-placement versus freely chosen drawings. The mission's definitions make those choices explicit and provide reusable segment-obstacle semantics for later geometric graph formalizations.

Difficulty

A finite combinatorial graph does not come with a canonical visibility drawing. Even when one starts with an arbitrary connected blocking set, replacing it by one bounded simple polygon requires compactness, component, incidence, and polygonal-neighborhood arguments. Conversely, lower bounds must quantify over every possible placement and obstacle, not merely refute a selected coordinate drawing.

Counting results for unrestricted graphs do not automatically preserve planarity. Bounds for planar obstacle number impose a crossing-free drawing and therefore answer a different question. The known two-obstacle examples close only the existential first part and give no universal kkk.

Formalization scope

All graph vertex types are finite. Obstacles are closed connected polygonal regions represented by finite triangle unions; they are pairwise disjoint and avoid graph vertices. Visibility uses the full closed segment, so tangency or boundary contact blocks a nonedge. The planarity witness is independent of the obstacle drawing. Empty and one-vertex graphs remain in the universal quantifier and should be handled without division or nonemptiness assumptions.

The repository's fixed-placement polygonization argument and finite arrangement code are candidate_only. They may motivate supporting lemmas, but they neither prove the unrestricted obstacle-drawing completeness theorem nor settle the universal bound. Contributions are welcome on exact geometry primitives, the published two-obstacle construction and lower bound, conversions between connected blockers and polygonal obstacles, and the universal root. A proof for the plane invariant, convex invariant, one fixed drawing, or a finite order cutoff must be labeled at that narrower scope.

Selected references

  • L. W. Berman, G. G. Chappell, J. R. Faudree, J. Gimbel, C. Hartman, and G. I. Williams, Graphs with Obstacle Number Greater than One, JGAA 21(6), 2017. https://doi.org/10.7155/jgaa.00452
  • J. Gimbel, P. Ossona de Mendez, and P. Valtr, Obstacle Numbers of Planar Graphs, Graph Drawing 2017. https://arxiv.org/abs/1706.06992
  • M. Balko, S. Chaplick, R. Ganian, S. Gupta, M. Hoffmann, P. Valtr, and A. Wolff, Bounding and Computing Obstacle Numbers of Graphs, SIAM Journal on Discrete Mathematics 38(2), 2024. https://arxiv.org/abs/2206.15414
  • Open Problem Garden / UnsolvedMath, OPG-37357. https://www.unsolvedmath.com/problems/OPG-37357
6 thms4 active usersReviewed
CombinatoricsOptimization·Captain: hao jia

Clique Partitions of Chordal Graphs (Erdos Problem 81)Open Problem

Motivation

An edge partition into cliques compresses the adjacency structure of a graph into complete pieces without allowing any edge to be counted twice. Erdős Problem 81 asks for the asymptotically sharp upper bound on the number of pieces needed when the graph is chordal. Chordal graphs have strong elimination structure, but that structure does not make the partition parameter additive under arbitrary edge deletion, and obtaining a linear error term remains substantially stronger than identifying the leading quadratic coefficient.

Erdős, Ordman, and Zalcstein studied clique partitions of chordal graphs in 1993. Their examples already exhibit the n2/6n^2/6n2/6 scale, while their general upper estimate had a larger quadratic coefficient. Later dense-packing results of Haxell–Rödl and Yuster compare fractional and integer triangle packings with an o(n2)o(n^2)o(n2) gap. The project candidate combines that interface with chordal elimination arguments to formulate a uniform n2/6+o(n2)n^2/6+o(n^2)n2/6+o(n2) milestone. It does not supply the O(n)O(n)O(n) remainder asked for by the root.

Setting

A finite simple graph is chordal when it has no induced cycle of length greater than three. The Lean definition uses the equivalent perfect-elimination form: vertices admit an injective ranking such that the later neighbors of every vertex form a clique.

An edge partition into cliques is a finite family P\mathcal PP of complete vertex sets such that every edge of GGG belongs to exactly one member of P\mathcal PP. Members may share vertices but may not share edges. Write cp⁡(G)\operatorname{cp}(G)cp(G) for the minimum possible number of pieces.

The asymptotic notation

n26+O(n)\frac{n^2}{6}+O(n)6n2​+O(n)

means that there are constants C>0C>0C>0 and n0≥1n_0\ge1n0​≥1, chosen independently of GGG and nnn, such that every chordal nnn-vertex graph with n≥n0n\ge n_0n≥n0​ has a clique partition with at most n2/6+Cnn^2/6+Cnn2/6+Cn pieces.

Formalization targets

Erdős Problem 81

The root theorem is

∃C>0 ∃n0≥1 ∀n≥n0 ∀G chordal on n vertices,cp⁡(G)≤n26+Cn.\exists C>0\ \exists n_0\ge1\ \forall n\ge n_0\ \forall G\text{ chordal on }n\text{ vertices}, \qquad \operatorname{cp}(G)\le \frac{n^2}{6}+Cn.∃C>0 ∃n0​≥1 ∀n≥n0​ ∀G chordal on n vertices,cp(G)≤6n2​+Cn.

The quantifier order is essential: CCC and n0n_0n0​ are universal and cannot depend on the graph.

Leading-coefficient milestone

The supporting target records the weaker uniform statement

∀ε>0 ∃n0 ∀n≥n0 ∀G chordal on n vertices,cp⁡(G)≤(16+ε)n2.\forall\varepsilon>0\ \exists n_0\ \forall n\ge n_0\ \forall G\text{ chordal on }n\text{ vertices}, \qquad \operatorname{cp}(G)\le \left(\frac16+\varepsilon\right)n^2.∀ε>0 ∃n0​ ∀n≥n0​ ∀G chordal on n vertices,cp(G)≤(61​+ε)n2.

This is the precise n2/6+o(n2)n^2/6+o(n^2)n2/6+o(n2) form. It is not equivalent to the root: choosing ε=1/n\varepsilon=1/nε=1/n is invalid because the cutoff may depend on the fixed value of ε\varepsilonε.

Significance

The root would determine the clique-partition extremum for chordal graphs up to a linear remainder, matching the scale of the complete-split examples that motivate the coefficient 1/61/61/6. It would refine a leading-order asymptotic theorem into a uniform estimate strong enough to distinguish second-order behavior.

Formalization creates a clean interface among perfect elimination orderings, exact edge partitions, fractional edge-and-triangle decompositions, and integer triangle packings. It also forces the proof to distinguish a partition from a cover and original graph order from the order of any auxiliary hypergraph. These definitions can support other decomposition problems on chordal and split graphs.

Difficulty

Perfect elimination does not by itself give the sharp partition count. Greedily taking maximal cliques may overlap in edges or accumulate too many singleton pieces. Similarly, a fractional edge-and-triangle partition can achieve the right leading coefficient while integer rounding loses o(n2)o(n^2)o(n2) pieces; the root requires that loss to be only O(n)O(n)O(n).

The dense-packing theorem has quantifiers of the form “for every fixed ε>0\varepsilon>0ε>0 there exists N(ε)N(\varepsilon)N(ε).” It therefore yields a uniform subquadratic error but no linear error. Any proof of the root must add a chordal-specific rounding or extremal reduction rather than treating the general packing theorem as if its ε\varepsilonε could vary with nnn.

Formalization scope

Graphs are finite and simple. Chordality is encoded by existence of a perfect-elimination ranking, including disconnected and edgeless graphs. A clique piece is a finite vertex set that spans a complete subgraph. Exactness means every actual edge occurs in exactly one piece; no nonedge can occur inside a piece. Bounds are compared in R\mathbb RR so the displayed asymptotic expressions retain their conventional form, while the number of parts remains a natural number.

The candidate derivation of the leading coefficient imports finite linear-programming duality and the Haxell–Rödl/Yuster fixed-triangle packing approximation. It is candidate_only, not an admitted result or kernel proof. Contributions may formalize the perfect-elimination lemmas, the fractional compression, the uniform packing interface, complete-split lower examples, or the root linear rounding theorem. A result for edge-and-triangle pieces only, a fractional partition, or one fixed order must not be presented as the unrestricted integer clique-partition theorem.

Selected references

  • P. Erdős, E. T. Ordman, and Y. Zalcstein, Clique Partitions of Chordal Graphs, Combinatorics, Probability and Computing 2(4), 1993. https://doi.org/10.1017/S0963548300000808
  • P. E. Haxell and V. Rödl, Integer and Fractional Packings in Dense Graphs, Combinatorica 21, 2001. https://doi.org/10.1007/s004930170003
  • R. Yuster, Integer and fractional packing of families of graphs, 2003. https://arxiv.org/abs/math/0305350
  • Erdős Problems, Problem 81. https://www.erdosproblems.com/81
16 thms4 active usersReviewed
Information TheoryQuantum Error Correction·Captain: mikedeng1

Quantum Expander Codes: On a Biregular Left-Right Expanding Graph with δA, δB < 1/6, the Greedy Generator Decoder Corrects Every Quantum Error of Weight Below min(γA·nA, γB·nB)/(3(1+ΔB))Research Paper

Motivation

A quantum LDPC code is a CSS code whose two check matrices HXH_XHX​ and HZH_ZHZ​ have rows and columns of bounded weight, so that every check measures only a few qubits. Such codes are the main candidates for fault-tolerant quantum memories, because their checks are local, but a code is only useful together with a decoder: an efficient procedure that recovers the error, up to stabilizers, from its syndrome. For classical LDPC codes built from expander graphs, Sipser and Spielman showed in 1996 that a simple bit-flipping algorithm corrects a constant fraction of the minimum distance in linear time (Sipser–Spielman 1996). In the quantum setting the obvious analogue fails, and before this work no family of quantum LDPC codes was known to admit an efficient decoder correcting a number of adversarial errors proportional to its minimum distance.

Timeline.

  • 1996: Sipser and Spielman, expander codes and their linear-time bit-flipping decoder.
  • 2009: Tillich and Zémor introduce the hypergraph product, giving quantum LDPC codes of positive rate and minimum distance of order n\sqrt nn​ (arXiv:0903.0566).
  • 2015: Leverrier, Tillich and Zémor take the hypergraph product of an expanding biregular graph with itself and prove that a greedy generator-flipping decoder corrects every error of weight up to a constant fraction of the minimum distance (arXiv:1504.00822, FOCS 2015). This mission formalizes that result.

Setting

Let G=(A∪B,E)G = (A \cup B, \mathcal E)G=(A∪B,E) be a finite bipartite graph, nA=∣A∣n_A = |A|nA​=∣A∣, nB=∣B∣n_B = |B|nB​=∣B∣, with nB≤nAn_B \le n_AnB​≤nA​, and write a∼ba \sim ba∼b for adjacency. GGG is (ΔA,ΔB)(\Delta_A, \Delta_B)(ΔA​,ΔB​)-biregular if every vertex of AAA has ΔA\Delta_AΔA​ neighbours and every vertex of BBB has ΔB\Delta_BΔB​. For S⊆AS \subseteq AS⊆A let Γ(S)\Gamma(S)Γ(S) be its neighbourhood. GGG is (γA,δA)(\gamma_A, \delta_A)(γA​,δA​)-left-expanding if ∣Γ(S)∣≥(1−δA)ΔA∣S∣|\Gamma(S)| \ge (1 - \delta_A)\Delta_A |S|∣Γ(S)∣≥(1−δA​)ΔA​∣S∣ for every S⊆AS \subseteq AS⊆A with ∣S∣≤γAnA|S| \le \gamma_A n_A∣S∣≤γA​nA​; right expansion is the same on the BBB side, and left-right expanding means both. Γu(S)\Gamma_u(S)Γu​(S) denotes the unique neighbours of SSS, adjacent to exactly one vertex of SSS.

Let HHH be the nB×nAn_B \times n_AnB​×nA​ binary biadjacency matrix of GGG. The quantum code QG\mathcal Q_GQG​ has n=nA2+nB2n = n_A^2 + n_B^2n=nA2​+nB2​ qubits, indexed by A2∪B2A^2 \cup B^2A2∪B2 (written αa\alpha aαa and bβb\betabβ), and check matrices

HX=(1nA⊗H, HT⊗1nB),HZ=(H⊗1nA, 1nB⊗HT),H_X = (\mathbb 1_{n_A} \otimes H,\ H^T \otimes \mathbb 1_{n_B}), \qquad H_Z = (H \otimes \mathbb 1_{n_A},\ \mathbb 1_{n_B} \otimes H^T),HX​=(1nA​​⊗H, HT⊗1nB​​),HZ​=(H⊗1nA​​, 1nB​​⊗HT),

with rows indexed by A×BA \times BA×B and B×AB \times AB×A. The row bababa of HZH_ZHZ​ has support the Z-generator gba={αa:α∼b}∪{bβ:β∼a}g_{ba} = \{\alpha a : \alpha \sim b\} \cup \{b\beta : \beta \sim a\}gba​={αa:α∼b}∪{bβ:β∼a}. An XXX error is e∈F2A2∪B2e \in \mathbb F_2^{A^2 \cup B^2}e∈F2A2∪B2​, its syndrome is σX(e)=HXeT\sigma_X(e) = H_X e^TσX​(e)=HX​eT, and CZ⊥\mathcal C_Z^\perpCZ⊥​ is the space of sums of Z-generators. The weight w(eX,eZ)w(e_X, e_Z)w(eX​,eZ​) of a quantum error counts the qubits where eXe_XeX​ or eZe_ZeZ​ is nonzero.

The decoder. On a syndrome sss, one iteration picks a vector e1e_1e1​ whose support lies in some generator gbag_{ba}gba​, such that s1=s+σX(e1)s_1 = s + \sigma_X(e_1)s1​=s+σX​(e1​) has ∣s1∣<∣s∣|s_1| < |s|∣s1​∣<∣s∣ and the ratio (∣s∣−∣s1∣)/∣e1∣(|s| - |s_1|)/|e_1|(∣s∣−∣s1​∣)/∣e1​∣ is maximal among all such vectors. It iterates on the current syndrome, fails if it gets stuck at a nonzero syndrome, and otherwise outputs the sum eOe_OeO​ of the flipped vectors. The ZZZ decoder is the same with HXH_XHX​ and HZH_ZHZ​ exchanged.

Formalization targets

Goal: Theorem 2

Assume GGG is (ΔA,ΔB)(\Delta_A, \Delta_B)(ΔA​,ΔB​)-biregular and (γA,δA,γB,δB)(\gamma_A, \delta_A, \gamma_B, \delta_B)(γA​,δA​,γB​,δB​)-left-right expanding with

δA<16,δB<16,\delta_A < \tfrac16, \qquad \delta_B < \tfrac16,δA​<61​,δB​<61​,

and let (eX,eZ)(e_X, e_Z)(eX​,eZ​) be a quantum error with

w(eX,eZ)<w0=13(1+ΔB)min⁡(γAnA,γBnB).w(e_X, e_Z) < w_0 = \frac{1}{3(1 + \Delta_B)} \min(\gamma_A n_A, \gamma_B n_B).w(eX​,eZ​)<w0​=3(1+ΔB​)1​min(γA​nA​,γB​nB​).

Then on every run, the XXX decoder on input σX(eX)\sigma_X(e_X)σX​(eX​) never outputs a failure and its output satisfies eO+eX∈CZ⊥e_O + e_X \in \mathcal C_Z^\perpeO​+eX​∈CZ⊥​. The ZZZ decoder on σZ(eZ)\sigma_Z(e_Z)σZ​(eZ​) behaves symmetrically, with eO+eZ∈CX⊥e_O + e_Z \in \mathcal C_X^\perpeO​+eZ​∈CX⊥​.

Milestones

In the paper's order of proof:

  • Lemma 3, ∣Γu(SA)∣≥(1−2δA)ΔA∣SA∣|\Gamma_u(S_A)| \ge (1 - 2\delta_A)\Delta_A |S_A|∣Γu​(SA​)∣≥(1−2δA​)ΔA​∣SA​∣;
  • Proposition 4, distance at least γAnA\gamma_A n_AγA​nA​ for the classical code;
  • the distance half of eq. (3), dQ≥min⁡(d,dT)d_Q \ge \min(d, d^T)dQ​≥min(d,dT), already posed on the platform;
  • Corollary 5, dQ≥min⁡(γAnA,γBnB)d_Q \ge \min(\gamma_A n_A, \gamma_B n_B)dQ​≥min(γA​nA​,γB​nB​);
  • the column-weight bound ∣σX(e)∣≤ΔB∣e∣|\sigma_X(e)| \le \Delta_B |e|∣σX​(e)∣≤ΔB​∣e∣;
  • Lemma 7, existence of a critical generator;
  • the displays (19), (22), (23), (24) of Appendix B;
  • Lemma 8, a step of ratio at least 1/31/31/3 inside a critical generator;
  • Lemma 10, the XXX half of the goal.

Significance

The result. Theorem 2 gives the first family of quantum LDPC codes with an efficient decoder that corrects a number of adversarial errors proportional to the minimum distance. Here that number is of order n\sqrt nn​. The same analysis yields a robustness property: low-weight errors have syndromes of proportional weight (Corollary 9 of the paper). Later work on quantum expander codes builds on this analysis, for example on their decoding under random noise (arXiv:1711.08351) and on fault tolerance with constant overhead (arXiv:1808.03821).

Formalizing it. The result is proved on paper. As far as is known, it has no machine-checked proof. The mission produces a formal statement of the decoder and its correctness, and it also produces reusable pieces:

  • unique-neighbour expansion;
  • the distance of expander codes;
  • the distance of hypergraph product codes, tied to an open item of the Tillich–Zémor mission;
  • a generic formalization of greedy syndrome-weight decoding.

Difficulty

The obvious quantum analogue of Sipser–Spielman flips single qubits that lower the syndrome weight. It fails because of degeneracy: an error and its sum with a generator have the same syndrome. A small error can therefore look locally like a large one, and no single qubit need be worth flipping. The paper's decoder instead flips arbitrary patterns inside one generator. The difficulty is to show that some pattern always lowers the syndrome weight by a constant fraction of its size. The proof needs a combinatorial structure, the critical generator, whose existence uses unique-neighbour expansion on both sides of the graph, and a case analysis on its six parts. A second difficulty is that the decoder minimizes the syndrome weight, not the error weight, so the intermediate error can grow. Bounding that growth is what limits the guarantee to w0w_0w0​ rather than the distance.

Formalization scope

  • The graph is a relation adj : A → B → Prop on finite types, and HHH = bipMatrix adj : Matrix B A (ZMod 2).
  • The published TillichZemor.HGP.Setting is reused. The paper's HXH_XHX​ is TillichZemor.HGP.hz H H and its HZH_ZHZ​ is hx H H; the X/Z names are swapped relative to the 2009 paper, and this is checked against eq. (5) in a sanity file.
  • Qubits are (A × A) ⊕ (B × B), weights are hammingNorm, and CZ⊥\mathcal C_Z^\perpCZ⊥​ is CodingTheory.dualCode of ker H_Z.
  • Expansion constants are real and cardinalities are compared in R\mathbb RR. Distance statements are per-vector, so no ∞\infty∞-valued minimum distance is cast to a real number.
  • The decoder is defined once (Candidate, Step, Reach, DecodesCorrectly), generic in the syndrome and generator matrices. Its ratio is computed in Q\mathbb QQ.
  • Reach contains every state reachable under any tie-breaking. Correctness requires, at every reachable state, a further step while the syndrome is nonzero, and eO+ee_O + eeO​+e in the row space of the generator matrix once it is zero. Termination is automatic, because each step lowers a natural number.

Additions and omissions.

  • The hypothesis ΔA≥1\Delta_A \ge 1ΔA​≥1 is added throughout. With ΔA=0\Delta_A = 0ΔA​=0 the expansion conditions are vacuous, H=0H = 0H=0, and every error is invisible.
  • In Lemma 8 the vector e1e_1e1​ is required to be nonzero. With e1=0e_1 = 0e1​=0 the printed inequality reads 0≥00 \ge 00≥0.
  • Lemma 8 is posed in its first assertion only. The critical generator is taken with respect to a minimum-weight representative of the coset, as in Appendix B.
  • The clauses "runs in time linear in the code length" and "letting ΔA\Delta_AΔA​, ΔB\Delta_BΔB​ be fixed and allowing nA,nBn_A, n_BnA​,nB​ to grow" are not formalized. The goal holds for every finite graph meeting the hypotheses.

Trivializing formalizations ruled out. The goal names the paper's own decoder, not "there exists a decoder", which a brute-force coset decoder would satisfy. It quantifies over every run, not a single favourable one. It asserts both never-failing and correctness, not partial correctness alone. It uses HXH_XHX​ = hz H H, not the swapped matrix.

Welcome contributions. Proofs of the counting milestones (Lemma 3, Proposition 4, the column-weight bound, (19), (22)–(24)) are independent of one another. The distance bound of Tillich–Zémor, open on the platform, is the deepest dependency. Unique-neighbour expansion and the hypergraph-product distance bound are useful beyond this mission.

Selected references

  • A. Leverrier, J.-P. Tillich, G. Zémor, Quantum expander codes, FOCS 2015; arXiv:1504.00822v1 (2015). https://arxiv.org/abs/1504.00822
  • J.-P. Tillich, G. Zémor, Quantum LDPC codes with positive rate and minimum distance proportional to the square root of the blocklength, IEEE Trans. Inf. Theory 60(2), 2014; arXiv:0903.0566. https://arxiv.org/abs/0903.0566
  • M. Sipser, D. A. Spielman, Expander codes, IEEE Trans. Inf. Theory 42(6), 1996. https://doi.org/10.1109/18.556667
  • O. Fawzi, A. Grospellier, A. Leverrier, Efficient decoding of random errors for quantum expander codes, STOC 2018; arXiv:1711.08351. https://arxiv.org/abs/1711.08351
  • O. Fawzi, A. Grospellier, A. Leverrier, Constant overhead quantum fault-tolerance with quantum expander codes, FOCS 2018; arXiv:1808.03821. https://arxiv.org/abs/1808.03821
16 thms3 active usersReviewed
Combinatorics·Captain: mikedeng1

A Solution to the 1-2-3 Conjecture: Every Finite Graph Without a K2 Component Has an Edge-Weighting from {1,2,3} Whose Weighted Degrees Properly Color the VerticesResearch Paper

Motivation

Many graph-colouring questions ask whether a proper colouring can be produced indirectly, from local data, instead of being assigned vertex by vertex. In the 1-2-3 Conjecture, the colours are produced by weighting the edges of a graph and colouring each vertex with the sum of the weights at it. Karoński, Łuczak and Thomason asked in 2004 whether weights from {1,2,3}\{1,2,3\}{1,2,3} always suffice (J. Combin. Theory Ser. B 91 (2004)). The question started as a local version of the irregularity strength problem, where all vertices, not only adjacent ones, must receive different sums. It became one of the most studied problems on graph labellings. Surveys: Seamone (arXiv:1211.5122) and Grytczuk (DOI 10.1016/j.ejc.2020.103213).

Timeline.

  • 2004: Karoński, Łuczak, Thomason state the conjecture and prove it for 3-colourable graphs.
  • 2007: Addario-Berry, Dalal, McDiarmid, Reed and Thomason give the first finite bound, k=30k=30k=30; later k=16k=16k=16 (Addario-Berry et al.) and k=13k=13k=13 (Wang and Yu).
  • 2010: Kalkowski, Karoński and Pfender reach k=5k=5k=5 (J. Combin. Theory Ser. B 100 (2010)).
  • Before 2023: Przybyło proves the conjecture for ddd-regular graphs with d≥108d\ge 10^8d≥108 and for graphs with minimum degree δ(G)≥Clog⁡Δ(G)\delta(G)\ge C\log\Delta(G)δ(G)≥ClogΔ(G); Zhong proves it for graphs with minimum degree at least 0.99985∣V∣0.99985|V|0.99985∣V∣ (references [20], [21], [31] of the paper).
  • 2023: Keusch reaches k=4k=4k=4 with a maximum-cut and flow argument (Combinatorica 43 (2023), 651–658; reference [18] of the paper).
  • 2024: Keusch proves the conjecture, k=3k=3k=3 (J. Combin. Theory Ser. B (2024), arXiv:2303.02611).

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite simple graph and write N(v)N(v)N(v) for the neighbourhood of a vertex vvv. A kkk-edge-weighting is a function ω:E→{1,…,k}\omega:E\to\{1,\dots,k\}ω:E→{1,…,k}. The weighted degree of vvv is

sω(v)=∑w∈N(v)ω({v,w}).s_\omega(v)=\sum_{w\in N(v)}\omega(\{v,w\}).sω​(v)=w∈N(v)∑​ω({v,w}).

Two adjacent vertices v,wv,wv,w have a coloring conflict if sω(v)=sω(w)s_\omega(v)=s_\omega(w)sω​(v)=sω​(w). The weighting is vertex-coloring if no edge is a coloring conflict, that is, if sωs_\omegasω​ is a proper vertex colouring of GGG.

If GGG has a connected component consisting of a single edge {v,w}\{v,w\}{v,w} (a component isomorphic to K2K_2K2​), then sω(v)=sω(w)=ω({v,w})s_\omega(v)=s_\omega(w)=\omega(\{v,w\})sω​(v)=sω​(w)=ω({v,w}) for every ω\omegaω, so no weighting is vertex-coloring. This is the only obstruction.

The proof uses a few further objects:

  • for W⊆VW\subseteq VW⊆V, deg⁡W(v)=∣N(v)∩W∣\deg_W(v)=|N(v)\cap W|degW​(v)=∣N(v)∩W∣;
  • for disjoint S,T⊆VS,T\subseteq VS,T⊆V, the bipartite subgraph G(S,T)G(S,T)G(S,T) with vertex set S∪TS\cup TS∪T and the edges of GGG between SSS and TTT;
  • a good R-B-partition (Definition 6): a partition V=R∪BV=R\cup BV=R∪B with RRR independent, G(R,B)G(R,B)G(R,B) connected and ∣B∣|B|∣B∣ even;
  • a maximum cut C=(S,T)C=(S,T)C=(S,T) and an auxiliary unit-capacity flow network GC,F,σG_{C,F,\sigma}GC,F,σ​ built from it (Lemma 4).

The Lean names are wdeg G ω v for sω(v)s_\omega(v)sω​(v), IsWeighting G k ω, IsVertexColoring G ω, NoK2Component G and IsGoodPartition G U R B.

Formalization targets

Goal: Theorem 1 (p. 2)

If no connected component of GGG is isomorphic to K2K_2K2​, then there is ω:E→{1,2,3}\omega:E\to\{1,2,3\}ω:E→{1,2,3} with

∑u∈N(v)ω({u,v})≠∑u∈N(w)ω({u,w})for every edge {v,w}∈E.\sum_{u\in N(v)}\omega(\{u,v\})\ne\sum_{u\in N(w)}\omega(\{u,w\})\qquad\text{for every edge }\{v,w\}\in E.u∈N(v)∑​ω({u,v})=u∈N(w)∑​ω({u,w})for every edge {v,w}∈E.

The connected case (p. 14)

Every finite connected graph with ∣V∣≥3|V|\ge 3∣V∣≥3 has a vertex-coloring {1,2,3}\{1,2,3\}{1,2,3}-edge-weighting. This is the statement the case analysis of the proof establishes. Theorem 1 follows by treating each component separately.

Milestones along the proof

In the order the proof uses them:

  • Lemma 2: an independent set RRR with G(R,V∖R)G(R,V\setminus R)G(R,V∖R) connected, alternating along a shortest path.
  • Lemma 4: a flow of value ∣F∣|F|∣F∣ in GC,F,σG_{C,F,\sigma}GC,F,σ​.
  • Lemma 3: a weighting and odd designated colours fff, proper on GGG, with ∣sω+h−f∣=1|s_\omega+h-f|=1∣sω​+h−f∣=1.
  • Lemma 5: a bipartite weighting with prescribed parities on RRR and exact weighted degrees on BBB, plus a second version with extra control along a fixed path.
  • Lemmas 7, 8 and 9: the three basic situations in which a good R-B-partition of GGG, or of GGG minus one vertex, yields a vertex-coloring weighting.
  • Lemma 10: every connected graph of minimum degree at least 222 has an edge whose removal, together with its two ends, leaves a connected graph.

Significance

The result. Theorem 1 settles a conjecture that was open for nineteen years. The constant is optimal: the triangle K3K_3K3​ has no vertex-coloring {1,2}\{1,2\}{1,2}-weighting (two of its three edges get the same weight, and then two vertices tie), and deciding whether a vertex-coloring {1,2}\{1,2\}{1,2}-weighting exists is NP-complete (Dudek and Wajc, 2011). Theorem 1 also closes the sequence of upper bounds 30,16,13,6,5,430,16,13,6,5,430,16,13,6,5,4 on the number of weights needed. Variants remain open: total weightings with vertex weights {1,2}\{1,2\}{1,2} and edge weights {1,2}\{1,2\}{1,2} (the 1-2 Conjecture), and list versions.

Formalizing it. The theorem is proved, in a 16-page paper. As far as is known it has no machine-checked proof. A formal proof would verify the case analysis at the end of the paper (pp. 11–14, about fifteen sub-cases with three auxiliary graphs), which is where errors in such arguments usually hide. The flow lemma is quoted from the 2023 paper and not reproved here. Formalizing it means proving an integral max-flow statement for a specific unit-capacity network, and the result can be reused for other flow-based weighting arguments.

Difficulty

The obvious approach is to fix the weights greedily, edge by edge or vertex by vertex. This fails because changing one edge weight changes two weighted degrees at once, so a repair at one vertex can create a conflict at a neighbour. With only three weights there is no room to absorb this. Every earlier bound above 333 used the extra weight as slack.

The difficulty in the paper is parity. The construction needs the set BBB of blue vertices to have even cardinality, and ∣B∣|B|∣B∣ cannot always be made even. The proof therefore removes one or two vertices, applies one of three base lemmas, and re-inserts the vertices. That re-insertion may need one more round of weight changes along a path. Each sub-case has to keep track of parities, of exact weighted degrees on BBB, and of the connectivity of G(R,B)G(R,B)G(R,B) simultaneously.

Formalization scope

Representation.

  • Graphs are SimpleGraph V on a finite vertex type V with decidable adjacency. The paper's "graph" is a finite simple graph: weighted degrees are finite sums, and the proof uses a maximum cut and vertex orderings.
  • An edge-weighting is a function ω : Sym2 V → ℕ on unordered pairs. IsWeighting G k ω requires ω(e)∈{1,…,k}\omega(e)\in\{1,\dots,k\}ω(e)∈{1,…,k} for every edge eee of GGG; values on non-edges are unconstrained and never used.
  • The weighted degree sums over G.neighborFinset v. Since {u,v}={v,u}\{u,v\}=\{v,u\}{u,v}={v,u} in Sym2, the single sum equals both sides of Theorem 1's displayed inequality.
  • "No component isomorphic to K2K_2K2​" uses Mathlib's connected components and graph isomorphism. A sanity file proves it equivalent to "no edge both of whose ends have degree 1".
  • G(R,B)G(R,B)G(R,B) is Mathlib's G.between R B, induced on R∪BR\cup BR∪B wherever that vertex set is not all of VVV.

Conventions that differ from the page. There are two, both disclosed in the item statements:

  1. The vertex type is finite.
  2. The flow of Lemma 4 is integral. The page says "an sss-ttt-flow"; the paper uses integrality immediately afterwards, and the two versions are equivalent by the integrality theorem.

Lemma 5's condition α(v)∈{2deg⁡(v)−1,2deg⁡(v),2deg⁡(v)+1}\alpha(v)\in\{2\deg(v)-1,2\deg(v),2\deg(v)+1\}α(v)∈{2deg(v)−1,2deg(v),2deg(v)+1} is stated in Z\mathbb ZZ, and Lemma 3's ∣sω+h−f∣=1|s_\omega+h-f|=1∣sω​+h−f∣=1 is computed in Z\mathbb ZZ.

Ruled out. None of the following formalizations is faithful:

  • weights allowed to be 000, which deletes edges;
  • weights constrained off the edges;
  • the K2K_2K2​ hypothesis strengthened to minimum degree at least 222, or to connectivity;
  • "adjacent vertices" replaced by "all pairs of vertices", which is the irregularity-strength problem.

Infrastructure. A complete development needs:

  • spanning trees and parity arguments on trees (Lemma 5);
  • depth-first search trees (Lemma 10);
  • maximum cuts and an integral max-flow or augmenting-path argument for GC,F,σG_{C,F,\sigma}GC,F,σ​ (Lemma 4), together with path decompositions of unit flows (the proof of Lemma 3).

The flow and cut definitions are stated for arbitrary finite graphs and can be reused. Contributions are welcome at every level: proofs of individual lemmas, general flow-decomposition results, and the final case analysis.

Selected references

  • R. Keusch, A Solution to the 1-2-3 Conjecture, J. Combin. Theory Ser. B, 2024. arXiv:2303.02611v4, DOI 10.1016/j.jctb.2024.01.002
  • R. Keusch, Vertex-coloring graphs with 4-edge-weightings, Combinatorica 43 (2023), 651–658. Cited as [18] in arXiv:2303.02611v4.
  • M. Karoński, T. Łuczak, A. Thomason, Edge weights and vertex colours, J. Combin. Theory Ser. B 91 (2004), 151–157. DOI 10.1016/j.jctb.2003.12.001
  • M. Kalkowski, M. Karoński, F. Pfender, Vertex-coloring edge-weightings: towards the 1-2-3-conjecture, J. Combin. Theory Ser. B 100 (2010). DOI 10.1016/j.jctb.2009.06.002
  • A. Dudek, D. Wajc, On the complexity of vertex-coloring edge-weightings, Discrete Math. Theor. Comput. Sci. 13 (2011), 45–50. Cited as [11] in arXiv:2303.02611v4.
  • B. Seamone, The 1-2-3 Conjecture and related problems: a survey, 2012. arXiv:1211.5122
  • J. Grytczuk, From the 1-2-3 Conjecture to the Riemann Hypothesis, European J. Combin. 91 (2021), 103213. DOI 10.1016/j.ejc.2020.103213
14 thms3 active usersReviewed
Combinatorics·Captain: mikedeng1

Generalized Turán Problems for Even Cycles 6: The Maximum Number of (2l+1)-Cycles in a Graph with No Cycle of Length 3, …, 2l or 2k+1 Is Ω(n^{1+1/(2k+2)}) and O(n²)Research Paper

Motivation

Generalized Turán problems ask for the maximum number of copies of a fixed graph HHH in an nnn-vertex graph that avoids a family F\mathcal FF of forbidden subgraphs. The ordinary Turán number counts edges (H=K2H = K_2H=K2​); Alon and Shikhelman (arXiv:1506.03397) started the systematic study of the general function ex(n,H,F)\mathrm{ex}(n, H, \mathcal F)ex(n,H,F), building on earlier results of Erdős on cliques, Bollobás and Győri on pentagons in triangle-free graphs, and Győri and Li on triangles in graphs without a long odd or even cycle. Counting cycles in graphs with forbidden cycles is closely tied to extremal problems for Berge hypergraphs, and the constructions used are sparse hypergraphs of large girth.

Gerbner, Győri, Methuku and Vizer, Generalized Turán problems for even cycles (arXiv:1712.07079, J. Combin. Theory Ser. B, 2020), determine or bound ex(n,Cℓ,F)\mathrm{ex}(n, C_{\ell}, \mathcal F)ex(n,Cℓ​,F) for many cycle families. This mission covers one of their results in §6.2: the number of odd cycles C2l+1C_{2l+1}C2l+1​ in a graph in which every shorter cycle is forbidden and, in addition, one longer odd cycle C2k+1C_{2k+1}C2k+1​ is forbidden. Its companion (Theorem 17, a separate mission of this series) forbids a longer even cycle C2kC_{2k}C2k​ instead, and there the paper obtains a sub-quadratic upper bound; with an odd extra cycle only a quadratic upper bound is proved.

Setting

All graphs are finite and simple. For an integer a≥3a \ge 3a≥3, CaC_aCa​ is the cycle of length aaa. For a set AAA of integers a≥3a \ge 3a≥3, write CA={Ca:a∈A}\mathcal C_A = \{C_a : a \in A\}CA​={Ca​:a∈A}; a graph is CA\mathcal C_ACA​-free if it contains no CaC_aCa​ with a∈Aa \in Aa∈A as a (not necessarily induced) subgraph. The family C2l\mathcal C_{2l}C2l​ consists of all cycles C3,C4,…,C2lC_3, C_4, \dots, C_{2l}C3​,C4​,…,C2l​.

For graphs HHH and GGG, N(H,G)\mathcal N(H, G)N(H,G) is the number of copies of HHH in GGG: subgraphs of GGG isomorphic to HHH, counted without labelling. The generalized Turán number is

ex(n,H,F)=max⁡{N(H,G):G has n vertices and contains no member of F}.\mathrm{ex}(n, H, \mathcal F) = \max\{\mathcal N(H, G) : G \text{ has } n \text{ vertices and contains no member of } \mathcal F\}.ex(n,H,F)=max{N(H,G):G has n vertices and contains no member of F}.

For a vertex vvv of GGG, Ni(v)N_i(v)Ni​(v) is the set of vertices at distance exactly iii from vvv.

A hypergraph is a finite set of hyperedges (vertex sets); it is rrr-uniform if each hyperedge has rrr vertices. A Berge cycle of length j≥2j \ge 2j≥2 is a sequence v1,h1,v2,h2,…,vj,hj,v1v_1, h_1, v_2, h_2, \dots, v_j, h_j, v_1v1​,h1​,v2​,h2​,…,vj​,hj​,v1​ of distinct vertices and distinct hyperedges with vi,vi+1∈hiv_i, v_{i+1} \in h_ivi​,vi+1​∈hi​ and vj,v1∈hjv_j, v_1 \in h_jvj​,v1​∈hj​. The girth of a hypergraph is the length of its shortest Berge cycle; "girth at least sss" means there is no Berge cycle of length 2,…,s−12, \dots, s-12,…,s−1.

Formalization targets

Goal: Theorem 18

For all integers k>l≥2k > l \ge 2k>l≥2, as n→∞n \to \inftyn→∞,

Ω(n1+12k+2)=ex(n,C2l+1,C2l∪{C2k+1})=O(n2).\Omega\big(n^{1+\frac{1}{2k+2}}\big) = \mathrm{ex}\big(n, C_{2l+1}, \mathcal C_{2l} \cup \{C_{2k+1}\}\big) = O(n^2).Ω(n1+2k+21​)=ex(n,C2l+1​,C2l​∪{C2k+1​})=O(n2).

Both the lower and the upper bound are part of the goal; the implied constants may depend on kkk and lll only.

Milestones

  1. Theorem 29 (Nešetřil–Rödl). For r≥2r \ge 2r≥2, s≥3s \ge 3s≥3 and all large nnn there is an rrr-uniform hypergraph on nnn vertices of girth at least sss with at least n1+1/sn^{1+1/s}n1+1/s hyperedges.
  2. The construction. Replacing each hyperedge of a (2l+1)(2l+1)(2l+1)-uniform hypergraph of girth at least 2k+22k+22k+2 by a copy of C2l+1C_{2l+1}C2l+1​ yields a graph with no cycle of length at most 2k+12k+12k+1 other than 2l+12l+12l+1, and with at least as many copies of C2l+1C_{2l+1}C2l+1​ as hyperedges.
  3. Cycles through a vertex. In a C2l∪{C2k+1}\mathcal C_{2l} \cup \{C_{2k+1}\}C2l​∪{C2k+1​}-free graph there is no edge inside Ni(v)N_i(v)Ni​(v) for i<li < li<l, every C2l+1C_{2l+1}C2l+1​ through vvv uses exactly one edge inside Nl(v)N_l(v)Nl​(v), and so the number of C2l+1C_{2l+1}C2l+1​'s through vvv is at most the number of edges inside Nl(v)N_l(v)Nl​(v).
  4. Claim 19. In a C2l∪{C2k+1}\mathcal C_{2l} \cup \{C_{2k+1}\}C2l​∪{C2k+1​}-free graph the number of edges inside Nl(v)N_l(v)Nl​(v) is at most C(k,l) ∣Nl(v)∣≤C(k,l) nC(k,l)\,|N_l(v)| \le C(k,l)\, nC(k,l)∣Nl​(v)∣≤C(k,l)n.

The supporting published theorem erdos_gallai_path_conjecture (the Erdős–Gallai path theorem: a graph on nnn vertices with no path of m+1m+1m+1 edges has at most mn/2mn/2mn/2 edges) is included as a reference item; the paper's proof of Claim 19 invokes it.

Significance

Theorem 18 shows that forbidding all cycles up to length 2l2l2l together with one odd cycle C2k+1C_{2k+1}C2k+1​ still permits polynomially more than nnn copies of C2l+1C_{2l+1}C2l+1​, at least n1+1/(2k+2)n^{1 + 1/(2k+2)}n1+1/(2k+2), and at most O(n2)O(n^2)O(n2). Compared with Theorem 17, where an even cycle C2kC_{2k}C2k​ is forbidden and the count is O(n1+l/(l+1))O(n^{1+l/(l+1)})O(n1+l/(l+1)), it isolates the role of the parity of the extra forbidden cycle. The authors conjecture that the true order is sub-quadratic (Conjecture 34, §8 of the paper), so the gap between the two exponents in Theorem 18 is a stated open question; a formal development of the known bounds is the base on which any improvement would be checked.

The result is proved in the paper; to the best of current knowledge it has no machine-checked proof. Formalizing it produces: a reusable encoding of generalized Turán numbers for cycle families on top of Mathlib's copyCount; Berge cycles and the girth of uniform hypergraphs; distance layers Ni(v)N_i(v)Ni​(v); and a formal statement of the Nešetřil–Rödl existence theorem, which is used as a black box by many extremal constructions.

Difficulty

The upper bound requires careful bookkeeping of distance layers around a vertex in a graph without short cycles. The obvious estimate, bounding the copies of C2l+1C_{2l+1}C2l+1​ through vvv by the number of pairs of vertices of Nl(v)N_l(v)Nl​(v), gives only ∣Nl(v)∣2|N_l(v)|^2∣Nl​(v)∣2 per vertex and hence nothing better than the trivial bound overall; the content of Claim 19 is that forbidding the single cycle C2k+1C_{2k+1}C2k+1​ makes the edge count inside Nl(v)N_l(v)Nl​(v) linear, and the parity of 2k+12k+12k+1 is used essentially.

The lower bound rests on Theorem 29, a cited result that is not proved in the paper, and for which Mathlib has no supporting theory of uniform hypergraphs or Berge cycles. For the construction itself, the paper only asserts that the cycle-replacement graph has no short cycle other than the inserted C2l+1C_{2l+1}C2l+1​'s; a cycle of the graph may wander through many hyperedges, and relating it to the girth of the hypergraph is the step a formal proof has to supply.

Formalization scope

  • Vertices of an nnn-vertex graph are Fin n; graphs are SimpleGraph (Fin n). Cycles are Mathlib's SimpleGraph.cycleGraph, and N(H,G)\mathcal N(H, G)N(H,G) is Mathlib's unlabelled copyCount.
  • C2l∪{C2k+1}\mathcal C_{2l} \cup \{C_{2k+1}\}C2l​∪{C2k+1​} is encoded as the set of lengths Set.Icc 3 (2 * l) ∪ {2 * k + 1}: every length 3,…,2l3, \dots, 2l3,…,2l, not only 2l2l2l.
  • ex(n,H,CA)\mathrm{ex}(n, H, \mathcal C_A)ex(n,H,CA​) is the supremum of a nonempty finite set of natural numbers, hence a true maximum (checked in a sanity file).
  • Ni(v)N_i(v)Ni​(v) uses the extended distance edist, so vertices of other components lie in no layer.
  • Hypergraphs are Finset (Finset V); "girth at least sss" is a predicate (no Berge cycle of length 2,…,s−12, \dots, s-12,…,s−1), not an infimum, so a hypergraph without Berge cycles has every girth.
  • Ω\OmegaΩ and OOO are Mathlib's Asymptotics.IsBigO along atTop, with real exponents; both directions of the goal are stated.
  • Claim 19's O(⋅)O(\cdot)O(⋅) is a single constant chosen before nnn, the graph and the vertex. The milestone on cycles through vvv uses the page's C2l∪{C2k+1}\mathcal C_{2l} \cup \{C_{2k+1}\}C2l​∪{C2k+1​}-freeness assumption.

A formalization that states only the upper bound, uses labelled copies, lets a constant depend on the graph, omits the short cycles C3,…,C2l−1C_3, \dots, C_{2l-1}C3​,…,C2l−1​ from the forbidden family, or defines ex\mathrm{ex}ex by a supremum that can take a junk value, does not prove this theorem.

Contributions welcome: proofs of the milestones, a proof of Theorem 29, and general lemmas on breadth-first layers in graphs of large girth, which are reusable by the companion missions of this series.

Selected references

  • D. Gerbner, E. Győri, A. Methuku, M. Vizer, Generalized Turán problems for even cycles, J. Combin. Theory Ser. B (2020); arXiv:1712.07079v3. https://arxiv.org/abs/1712.07079, https://doi.org/10.1016/j.jctb.2020.05.005
  • N. Alon, C. Shikhelman, Many T copies in H-free graphs, J. Combin. Theory Ser. B 121 (2016). https://arxiv.org/abs/1506.03397
  • J. Nešetřil, V. Rödl, On a probabilistic graph-theoretical method, Proc. Amer. Math. Soc. 72(2), 417–421 (1978). https://www.ams.org/journals/proc/1978-072-02/
  • P. Erdős, T. Gallai, On maximal paths and circuits of graphs, Acta Math. Acad. Sci. Hungar. 10 (1959). https://doi.org/10.1007/BF02024498
7 thms3 active usersReviewed
CombinatoricsOperations Research·Captain: mikedeng1

An Introduction to Timetabling 1: Every Class–Teacher Requirement Matrix Has a Timetable over Any p Days with Every Daily Load and Pair Count Between ⌊r/p⌋ and ⌈r/p⌉Research Paper

Motivation

School and university timetabling is one of the oldest applications of combinatorial optimization. In his invited review An introduction to timetabling (European Journal of Operational Research, 1985), D. de Werra opens with the class–teacher model: classes meet teachers for a prescribed number of one-period lectures, and the lectures must be placed into periods or days so that no class and no teacher is overloaded. Every richer timetabling model of the review (preassignments, unavailabilities, course scheduling) is built on top of this one, and the review keeps returning to its graph-theoretic reading as an edge-colouring problem in a bipartite multigraph.

The model comes in three levels of demand. The daily problem asks for a clash-free timetable. The weekly problem asks for an assignment of lectures to days within daily load limits. The balanced weekly problem asks, in addition, that the load of every class and every teacher be spread as evenly as possible over the week, and that the lectures of each class–teacher pair be spread evenly too. The last one is the practical requirement that a teacher should not have all four lectures with one class on Monday.

Timeline (as cited in the review).

  • 1916. D. König proves that a bipartite multigraph of maximum degree Δ\DeltaΔ can be edge-coloured with Δ\DeltaΔ colours. The review states it as Proposition 2.1, citing Berge's Graphes [1].
  • 1975. D. de Werra, "A few remarks on chromatic scheduling" (reference [37] of the review), gives the balanced version: for every ppp the edges of a bipartite multigraph can be ppp-coloured so that every node and every bundle of parallel edges is balanced. This is Proposition 2.3.
  • 1978. D. de Werra, "Some comments on a note about timetabling" (INFOR, reference [38]), gives the weekly version with daily load limits, Proposition 2.2.
  • 1985. The review collects the three results in matrix form as problems CT1, CT2 and CT3 and sketches a network-flow construction for CT3, citing Krarup [19].

Setting

There are mmm classes c1,…,cmc_1,\dots,c_mc1​,…,cm​ and nnn teachers t1,…,tnt_1,\dots,t_nt1​,…,tn​. The requirement matrix R=(rij)R=(r_{ij})R=(rij​) is an m×nm\times nm×n matrix of nonnegative integers: rijr_{ij}rij​ is the number of lectures class cic_ici​ must have with teacher tjt_jtj​. Write ri⋅=∑jrijr_{i\cdot}=\sum_j r_{ij}ri⋅​=∑j​rij​ for the total load of class cic_ici​ and r⋅j=∑irijr_{\cdot j}=\sum_i r_{ij}r⋅j​=∑i​rij​ for that of teacher tjt_jtj​. For nonnegative integers sss and p≥1p\ge1p≥1, ⌊s/p⌋\lfloor s/p\rfloor⌊s/p⌋ and ⌈s/p⌉\lceil s/p\rceil⌈s/p⌉ are the floor and ceiling of the rational s/ps/ps/p.

A schedule over ppp periods or days is an array x=(xijk)x=(x_{ijk})x=(xijk​) of nonnegative integers (i≤mi\le mi≤m, j≤nj\le nj≤n, k≤pk\le pk≤p) that places every lecture exactly once:

∑k=1pxijk=rijfor all i,j.(1)\sum_{k=1}^p x_{ijk}=r_{ij}\quad\text{for all } i,j. \tag{1}k=1∑p​xijk​=rij​for all i,j.(1)
  • CT1 asks in addition that xijk∈{0,1}x_{ijk}\in\{0,1\}xijk​∈{0,1} and that ∑jxijk≤1\sum_j x_{ijk}\le1∑j​xijk​≤1, ∑ixijk≤1\sum_i x_{ijk}\le1∑i​xijk​≤1: no class and no teacher has two lectures in one period.
  • CT2, given positive integers aia_iai​, bjb_jbj​, asks that ∑jxijk≤ai\sum_j x_{ijk}\le a_i∑j​xijk​≤ai​ and ∑ixijk≤bj\sum_i x_{ijk}\le b_j∑i​xijk​≤bj​ on every day kkk.
  • CT3 asks that on every day kkk
⌊ri⋅p⌋≤∑jxijk≤⌈ri⋅p⌉,⌊r⋅jp⌋≤∑ixijk≤⌈r⋅jp⌉,⌊rijp⌋≤xijk≤⌈rijp⌉.(8–10)\Big\lfloor \tfrac{r_{i\cdot}}{p}\Big\rfloor\le\sum_j x_{ijk}\le\Big\lceil \tfrac{r_{i\cdot}}{p}\Big\rceil,\quad \Big\lfloor \tfrac{r_{\cdot j}}{p}\Big\rfloor\le\sum_i x_{ijk}\le\Big\lceil \tfrac{r_{\cdot j}}{p}\Big\rceil,\quad \Big\lfloor \tfrac{r_{ij}}{p}\Big\rfloor\le x_{ijk}\le\Big\lceil \tfrac{r_{ij}}{p}\Big\rceil. \tag{8–10}⌊pri⋅​​⌋≤j∑​xijk​≤⌈pri⋅​​⌉,⌊pr⋅j​​⌋≤i∑​xijk​≤⌈pr⋅j​​⌉,⌊prij​​⌋≤xijk​≤⌈prij​​⌉.(8–10)

In the Lean development these are IsCT1 R p x, IsCT2 R p a b x and IsCT3 R p x, with the floor and ceiling written lo s p and hi s p and the totals rowSum R i, colSum R j.

Formalization targets

Goal: Proposition 2.3

∀ R∈Nm×n, ∀ p≥1:∃ x,  x solves CT3 for R over p days.\forall\,R\in\mathbb N^{m\times n},\ \forall\,p\ge1:\quad \exists\,x,\ \ x \text{ solves CT3 for } R \text{ over } p \text{ days}.∀R∈Nm×n, ∀p≥1:∃x,  x solves CT3 for R over p days.

The statement fixes no constant: it asserts that perfectly balanced weekly timetables always exist.

Milestones

  1. One balanced day (flow sketch, p. 153): for p≥1p\ge1p≥1 there is a matrix yyy satisfying the bounds (8), (9), (10) at once.
  2. The remaining days stay balanced (same sketch): if p≥2p\ge2p≥2 and ⌊s/p⌋≤y≤⌈s/p⌉\lfloor s/p\rfloor\le y\le\lceil s/p\rceil⌊s/p⌋≤y≤⌈s/p⌉, then ⌊s/p⌋≤⌊(s−y)/(p−1)⌋\lfloor s/p\rfloor\le\lfloor (s-y)/(p-1)\rfloor⌊s/p⌋≤⌊(s−y)/(p−1)⌋ and ⌈(s−y)/(p−1)⌉≤⌈s/p⌉\lceil (s-y)/(p-1)\rceil\le\lceil s/p\rceil⌈(s−y)/(p−1)⌉≤⌈s/p⌉.
  3. Proposition 2.1 (König): CT1 is solvable iff r⋅j≤pr_{\cdot j}\le pr⋅j​≤p and ri⋅≤pr_{i\cdot}\le pri⋅​≤p for all i,ji,ji,j.
  4. Proposition 2.2: CT2 is solvable iff r⋅j≤p bjr_{\cdot j}\le p\,b_jr⋅j​≤pbj​ and ri⋅≤p air_{i\cdot}\le p\,a_iri⋅​≤pai​.
  5. Minimum number of days: CT2 is solvable over ppp days iff p≥max⁡(max⁡j⌈r⋅j/bj⌉,max⁡i⌈ri⋅/ai⌉)p\ge\max\big(\max_j\lceil r_{\cdot j}/b_j\rceil,\max_i\lceil r_{i\cdot}/a_i\rceil\big)p≥max(maxj​⌈r⋅j​/bj​⌉,maxi​⌈ri⋅​/ai​⌉).
  6. Figure 1: the printed timetable for R=(124302)R=\begin{pmatrix}1&2&4\\3&0&2\end{pmatrix}R=(13​20​42​), p=3p=3p=3, solves CT3.

Significance

The result. Proposition 2.3 says that the balancing requirement costs nothing: whatever the requirement matrix, the week can be organised so that each class and each teacher has an essentially constant daily load and each class–teacher pair is spread evenly. It contains König's edge-colouring theorem (take ppp at least the maximum degree, so every daily load is at most one) and the weekly Proposition 2.2 as special cases. Combined with Proposition 2.1 applied to each day, it reduces weekly timetabling to independent daily problems, which is the decomposition the review recommends ("one may start by solving the weekly problem … and then solve separately the resulting daily problems", p. 154). In graph language it is the existence of equitable edge colourings of bipartite multigraphs, a standard tool in scheduling and in the theory of edge colourings.

Formalizing it. All results of this mission are classical and proved in the literature; none of them is formalized on the platform. Mathlib contains Hall's marriage theorem but not König's edge-colouring theorem for bipartite multigraphs, and not balanced or equitable colourings. A complete development gives a machine-checked König edge-colouring theorem in matrix form, the matrix-rounding step with simultaneous row, column and entry bounds, and the balanced-colouring theorem itself.

Difficulty

Each family of constraints alone is easy: spreading a single number rijr_{ij}rij​ evenly over ppp days, or a single row total, is a matter of division with remainder. The difficulty is that one array must satisfy (1), (8), (9) and (10) simultaneously, so that the per-entry roundings must add up to correct roundings of every row total and every column total on every day. Rounding each entry independently breaks the row and column constraints, and greedy day-by-day choices can leave a remainder that no longer fits the bounds. The crux is integrality: the rational average rij/pr_{ij}/prij​/p satisfies every bound on every day, but the problem asks for integers, and an integral one-day slice must round all entries, all row totals and all column totals consistently.

Formalization scope

Classes, teachers and days are Fin m, Fin n, Fin p; RRR and xxx are natural-number valued, which is the paper's "xijk≥0x_{ijk}\ge0xijk​≥0 integer". Constraint (4) of CT1 is written xijk≤1x_{ijk}\le1xijk​≤1. Floors and ceilings are Nat.floor and Nat.ceil of the rational quotient. In (8) and (9) they are applied to the row and column totals, never summed entrywise. The minimum number of days uses Finset.sup, whose empty value is 000, and is stated as an equivalence ("solvable over ppp days iff p≥pmin⁡p\ge p_{\min}p≥pmin​") rather than as an infimum.

Disclosed choices:

  • The goal and the one-day milestone assume p≥1p\ge1p≥1; the paper's "for any ppp" means any number of days, and for p=0p=0p=0 Lean's division returns 000.
  • The remaining-days milestone assumes p≥2p\ge2p≥2 so that p−1p-1p−1 days remain.
  • Propositions 2.1, 2.2 and the minimum-days statement hold for every ppp, including p=0p=0p=0, and assume nothing more than the page: ai,bja_i,b_jai​,bj​ positive.
  • The edge-colouring and network-flow readings of the paper are stated in prose only; the formal statements are in matrix form, which is how the propositions are stated.

The formalization does not admit a trivial reading: CT3 requires all four constraint families at once, (1) is an equality, and the sanity checks show both that the Figure 1 data satisfy CT3 and that CT1 is not satisfiable for R=(2)R=(2)R=(2) and p=1p=1p=1.

Useful infrastructure, reusable beyond this mission: integral flows with lower bounds, or total unimodularity of bipartite incidence matrices; König's edge-colouring theorem for bipartite multigraphs; floor and ceiling arithmetic for natural-number quotients. Proofs of any milestone, alternative proofs of the goal, and a sorry-free proof of the Figure 1 instance are all welcome.

Selected references

  • D. de Werra, An introduction to timetabling, European Journal of Operational Research 19 (1985) 151–162. https://doi.org/10.1016/0377-2217(85)90167-5
  • D. König, Über Graphen und ihre Anwendung auf Determinantentheorie und Mengenlehre, Mathematische Annalen 77 (1916) 453–465. https://doi.org/10.1007/BF01456961
  • C. Berge, Graphes, Gauthier-Villars, Paris, 1983 (reference [1] of the review).
  • D. de Werra, A few remarks on chromatic scheduling, in B. Roy (ed.), Combinatorial Programming: Methods and Applications, Reidel, Dordrecht, 1975, 337–342.
  • D. de Werra, Some comments on a note about timetabling, INFOR 16 (1978) 90–92.
  • J. Krarup, Chromatic optimisation: Limitations, objectives, uses, references, European Journal of Operational Research 11 (1982) 1–19.
8 thms3 active usersReviewed
Dynamic ProgrammingOperations ResearchOptimization·Captain: mikedeng1

The Complexity of Markov Decision Processes III: The Optimal Discounted Cost of a Deterministic Process Is the Least Discounted Cost of a SigmaResearch Paper

Motivation

A Markov decision process with finitely many states and decisions, a discount factor β∈(0,1)\beta\in(0,1)β∈(0,1) and the objective of minimizing expected discounted cost is the basic model of sequential decision making in operations research and reinforcement learning. It is solvable in polynomial time by linear programming. Papadimitriou and Tsitsiklis (Math. Oper. Res. 12(3), 1987) asked the finer question of whether it can be solved fast in parallel, and showed a sharp split: the general (stochastic) problem is P-complete (their Theorem 1), so it is unlikely to admit a polylogarithmic parallel algorithm, while the deterministic special case, where every transition probability is 0 or 1, is in NC (their Theorem 4).

The deterministic case is a shortest-path problem with a twist. Every policy traces a walk in a graph, and the cost of the walk counts its ttt-th arc with weight βt\beta^tβt. The paper's argument identifies the optimal discounted cost with the cost of a finite combinatorial object, a sigma, and computes it from discounted min-plus matrix products. This mission formalizes that identity.

Setting

A deterministic stationary Markov decision process consists of a finite set SSS of nnn states, an initial state s0∈Ss_0\in Ss0​∈S, and for each state sss a finite nonempty set DsD_sDs​ of decisions. Decision i∈Dsi\in D_si∈Ds​ costs c(s,i)∈Rc(s,i)\in\mathbb Rc(s,i)∈R and leads with certainty to the state next(s,i)\mathrm{next}(s,i)next(s,i). Equivalently, it is a directed graph on SSS with an arc s→next(s,i)s\to\mathrm{next}(s,i)s→next(s,i) of weight c(s,i)c(s,i)c(s,i) for each decision; parallel arcs and loops are allowed.

A policy δ\deltaδ chooses a decision δ(s,t)∈Ds\delta(s,t)\in D_sδ(s,t)∈Ds​ for every state sss and time t=0,1,2,…t=0,1,2,\dotst=0,1,2,…; it is stationary if δ(s,t)\delta(s,t)δ(s,t) does not depend on ttt. From s0s_0s0​ it produces the states st+1=next(st,δ(st,t))s_{t+1}=\mathrm{next}(s_t,\delta(s_t,t))st+1​=next(st​,δ(st​,t)) and the discounted cost

Jβ(δ)=∑t=0∞βtc(st,δ(st,t)).J_\beta(\delta)=\sum_{t=0}^{\infty}\beta^t c\bigl(s_t,\delta(s_t,t)\bigr).Jβ​(δ)=t=0∑∞​βtc(st​,δ(st​,t)).

The optimal discounted cost is Jβ∗(s0)=inf⁡δJβ(δ)J^*_\beta(s_0)=\inf_\delta J_\beta(\delta)Jβ∗​(s0​)=infδ​Jβ​(δ), the infimum over all policies δ(s,t)\delta(s,t)δ(s,t).

A sigma is a walk w0=s0,w1,…,wNw_0=s_0,w_1,\dots,w_Nw0​=s0​,w1​,…,wN​ along decisions et∈Dwte_t\in D_{w_t}et​∈Dwt​​ whose nodes w0,…,wN−1w_0,\dots,w_{N-1}w0​,…,wN−1​ are distinct and whose last node repeats an earlier one, wN=wjw_N=w_jwN​=wj​ with j<Nj<Nj<N: the walk from s0s_0s0​ up to the first repetition of a node. It is a prefix of jjj arcs followed by a cycle of l=N−j≥1l=N-j\ge1l=N−j≥1 arcs. Its discounted cost cβ(P)c_\beta(P)cβ​(P) is the discounted cost of the infinite walk that follows it and repeats its cycle forever.

In the min-plus semiring (addition is min⁡\minmin, multiplication is +++, zero is +∞+\infty+∞), let AuvA_{uv}Auv​ be the least cost of a decision from uuu to vvv (+∞+\infty+∞ if none), let rArArA multiply every finite entry by rrr, and let

B0=I,Bk=A⊗βA⊗⋯⊗βk−1A,B_0=I,\qquad B_k=A\otimes\beta A\otimes\cdots\otimes\beta^{k-1}A ,B0​=I,Bk​=A⊗βA⊗⋯⊗βk−1A,

with III the min-plus identity. The entry (Bk)uv(B_k)_{uv}(Bk​)uv​ is the least discounted cost ∑t<kβtct\sum_{t<k}\beta^t c_t∑t<k​βtct​ of a kkk-arc walk from uuu to vvv.

Formalization targets

Goal

For 0<β<10<\beta<10<β<1,

Jβ∗(s0)=min⁡P sigmacβ(P)=min⁡{(Bk)s0u+βk(Bl)uu1−βl  :  u∈S, 0≤k≤n, 1≤l≤n, both entries finite},J^*_\beta(s_0)=\min_{P\text{ sigma}}c_\beta(P)=\min\Bigl\{(B_k)_{s_0u}+\frac{\beta^{k}(B_l)_{uu}}{1-\beta^{l}}\;:\;u\in S,\ 0\le k\le n,\ 1\le l\le n,\ \text{both entries finite}\Bigr\},Jβ∗​(s0​)=P sigmamin​cβ​(P)=min{(Bk​)s0​u​+1−βlβk(Bl​)uu​​:u∈S, 0≤k≤n, 1≤l≤n, both entries finite},

with both minima attained. This is the claim on pp. 446–447 on which Theorem 4 rests. The algebraic expression is the corrected form of the paper's (see Formalization scope).

Milestones

  1. The closed form of the cost of a sigma (p. 446): cβ(P)=∑t<jβtct+βj1−βl∑t<lβtcj+tc_\beta(P)=\sum_{t<j}\beta^tc_t+\frac{\beta^j}{1-\beta^l}\sum_{t<l}\beta^tc_{j+t}cβ​(P)=∑t<j​βtct​+1−βlβj​∑t<l​βtcj+t​.
  2. A stationary optimal policy exists (§2, p. 444; cited on p. 446). The general finite stochastic statement is on the platform as BertsekasDP.discounted_main_theorem (Proved), and the deterministic case is stated in this mission's model.
  3. The path of a stationary policy, cut at its first repeated node, is a sigma with the same discounted cost, and every sigma arises this way (p. 446).
  4. The products BjB_jBj​ are the cheapest discounted jjj-arc walks (p. 447).

Significance

The identity reduces an optimization over infinite objects, the time-dependent policies, to a minimum over finitely many closed-form expressions in matrix entries. Each BkB_kBk​ is a product of kkk min-plus matrices, and matrix products parallelize well; that is the source of the NC algorithm. The same reduction, with the discounting removed, underlies the minimum-mean-cycle characterization of the average-cost problem (mission II of this series), and the sigma structure reappears in the stationary finite-horizon analysis (mission IV).

On the formalization side, the result has been proved in the paper only as a short paragraph, which cites the existence of stationary optimal policies and leaves the algebra to the reader. Its printed formulas contain index slips (below). A machine-checked version fixes the exact constants. To our knowledge none of these statements has a formal proof. The existence of stationary optimal policies for finite discounted problems is formalized on the platform in Bertsekas's general stochastic setting.

Difficulty

The obvious argument says: an optimal policy may be taken stationary, a stationary deterministic policy produces an eventually periodic path, and an eventually periodic path is a sigma. The first step is the substantial one: the optimum is an infimum over all time-dependent policies δ(s,t)\delta(s,t)δ(s,t), and a non-stationary policy can follow any walk at all, so this step is a statement about the whole dynamic program, not a combinatorial fact about a single walk. The second part of the goal adds a matching argument in the other direction: the expression ranges over arbitrary kkk-arc walks and closed lll-walks, which need not be simple, so one must check that none of them beats the best sigma, and that the best sigma has k≤nk\le nk≤n and l≤nl\le nl≤n. Finally, the discount factors must be tracked exactly; the printed formulas are off by one power of β\betaβ in two places.

Formalization scope

The Lean development lives in the namespace MDPComplexity.Discounted, with one definition file and five theorems. Conventions committed to:

  • The model is stationary and deterministic. A DetMDP S has a finite state type, a nonempty finite decision type at each state (the paper says "a finite set DsD_sDs​"; nonemptiness is implicit since a policy must choose a decision), a successor map and real costs.
  • Policies are maps δ(s,t)\delta(s,t)δ(s,t), time is N\mathbb NN from 000, and the optimum optDisc is the infimum over all such policies. Defining it over stationary policies or over sigmas would make the goal true by definition. That shortcut is ruled out, and the stationarity step is part of the content.
  • Costs are tsums, which are junk (000) for non-summable series; every statement that sums costs assumes 0<β<10<\beta<10<β<1, under which the series converge absolutely.
  • A sigma is encoded as a lasso, wN=wjw_N=w_jwN​=wj​ with j<Nj<Nj<N. This includes j=0j=0j=0, a sigma whose cycle passes through s0s_0s0​, which the printed form (u0,…,uk,v1,…,vl,v1)(u_0,\dots,u_k,v_1,\dots,v_l,v_1)(u0​,…,uk​,v1​,…,vl​,v1​) omits but the definition "a path from u0u_0u0​ until the first repetition of a node" includes. The paper's kkk is j−1j-1j−1.
  • Corrected slips: (i) in the sigma display the cycle sum carries βj−1\beta^{j-1}βj−1, not βj\beta^jβj, for its jjj-th arc, and the indices are made consistent; (ii) the algorithm's expression is λ+βkμ/(1−βl)\lambda+\beta^k\mu/(1-\beta^l)λ+βkμ/(1−βl) with l≥1l\ge1l≥1, not λ+βk+1μ/(1+βl)\lambda+\beta^{k+1}\mu/(1+\beta^l)λ+βk+1μ/(1+βl) (one state with a loop of cost 111 has optimum 1/(1−β)1/(1-\beta)1/(1−β), which the printed form misses); (iii) B0B_0B0​, used but not defined, is the min-plus identity.
  • Matrix entries are in Tropical (WithTop ℝ), so +∞+\infty+∞ is ⊤ and the product is Mathlib's matrix product over the tropical semiring.

Not formalized: membership in NC, the processor count n4n^4n4, the log⁡2n+3log⁡n\log^2 n+3\log nlog2n+3logn parallel time, and the P-completeness of the general problem. Mathlib has no parallel complexity classes. Only the identity that makes the algorithm correct is a target.

Contributions welcome: proofs of any milestone. A bridge from DetMDP to BertsekasSSPModel (taking decisions as controls and 0/1 transition rows) would let milestone 2 be derived from the platform's proved Proposition 7.3.1. A direct proof of the deterministic case is also welcome. The lasso and walk lemmas are reusable for the other missions of the series.

Selected references

  • C. H. Papadimitriou and J. N. Tsitsiklis, The Complexity of Markov Decision Processes, Mathematics of Operations Research 12(3), 441–450, 1987. https://doi.org/10.1287/moor.12.3.441
  • D. P. Bertsekas, Dynamic Programming and Optimal Control, Vol. I, 3rd ed., Athena Scientific, 2005, §7.3 (discounted problems, Proposition 7.3.1).
  • M. L. Puterman, Markov Decision Processes: Discrete Stochastic Dynamic Programming, Wiley, 1994, Ch. 6 (existence of stationary optimal policies for discounted models). https://doi.org/10.1002/9780470316887
7 thms3 active usersReviewed
Complexity TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

Complexity of Machine Scheduling Problems 6: DIRECTED HAMILTON PATH Reduces to No-Wait Flow Shop Makespan and Total Completion TimeResearch Paper

Motivation

In a no-wait flow shop every job passes through the machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​ in the same order, and once it has started it may never wait between two machines. The constraint comes from processes in which the material changes state if it is left standing, such as hot metal rolling, chemical and food processing, and some pharmaceutical lines; see the survey of Hall and Sriskandarajah (doi:10.1287/opre.44.3.510). The question of how hard it is to schedule such a shop well is the question this mission formalizes.

Brucker, Lenstra and Rinnooy Kan settled it for the case in which the number of machines is part of the input. In their report Complexity of Machine Scheduling Problems (Mathematisch Centrum, 1975; journal version in Annals of Discrete Mathematics 1, 1977, doi:10.1016/S0167-5060(08)70743-X), Theorem 5 reduces DIRECTED HAMILTON PATH to the no-wait flow shop. Both minimizing the makespan Cmax⁡C_{\max}Cmax​ and minimizing the total completion time ∑Cj\sum C_j∑Cj​ are thereby NP-complete.

Timeline.

  • 1964: Gilmore and Gomory solve the two-machine no-wait flow shop with makespan in polynomial time (doi:10.1287/opre.12.5.655).
  • 1972: Wismer (doi:10.1287/opre.20.3.689) and Reddi and Ramamoorthy (doi:10.1057/jors.1972.52) show that no-wait makespan minimization is a travelling-salesman problem with arc weights computed from the processing times.
  • 1972: Karp proves DIRECTED HAMILTON CIRCUIT NP-complete (doi:10.1007/978-1-4684-2001-2_9).
  • 1975: Brucker, Lenstra and Rinnooy Kan reduce DIRECTED HAMILTON CIRCUIT to DIRECTED HAMILTON PATH (their Theorem 2(d)), and DIRECTED HAMILTON PATH to n∣m∣F,no wait∣Cmax⁡n|m|F,\textit{no wait}|C_{\max}n∣m∣F,no wait∣Cmax​ and to n∣m∣F,no wait,wj=1∣∑wjCjn|m|F,\textit{no wait},w_j=1|\sum w_jC_jn∣m∣F,no wait,wj​=1∣∑wj​Cj​ (Theorem 5).
  • 1984: Röck shows that the problem stays NP-hard for three machines (doi:10.1145/62.65).

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines. Job JℓJ_\ellJℓ​ needs processing time pℓi∈Np_{\ell i}\in\mathbb Npℓi​∈N on machine MiM_iMi​. Write

qℓi=∑r=1ipℓr,qℓ0=0,q_{\ell i}=\sum_{r=1}^{i}p_{\ell r},\qquad q_{\ell 0}=0,qℓi​=r=1∑i​pℓr​,qℓ0​=0,

for the time job JℓJ_\ellJℓ​ spends on its first iii machines. Because a job never waits, a schedule is determined by the start times Bℓ∈NB_\ell\in\mathbb NBℓ​∈N. The operation of JℓJ_\ellJℓ​ on MiM_iMi​ occupies [Bℓ+qℓ,i−1, Bℓ+qℓi)[B_\ell+q_{\ell,i-1},\,B_\ell+q_{\ell i})[Bℓ​+qℓ,i−1​,Bℓ​+qℓi​), and job JℓJ_\ellJℓ​ completes at Cℓ=Bℓ+qℓmC_\ell=B_\ell+q_{\ell m}Cℓ​=Bℓ​+qℓm​. A schedule is feasible when no two distinct jobs occupy the same machine at the same time. The delay

cjk=max⁡1≤i≤m{qji−qk,i−1}c_{jk}=\max_{1\le i\le m}\{q_{ji}-q_{k,i-1}\}cjk​=1≤i≤mmax​{qji​−qk,i−1​}

is the least gap Bk−BjB_k-B_jBk​−Bj​ that lets JkJ_kJk​ follow JjJ_jJj​ on every machine.

A directed graph G=(V,A)G=(V,A)G=(V,A) on V={0,…,n−1}V=\{0,\dots,n-1\}V={0,…,n−1} has a Hamilton path if its vertices can be ordered σ(0),…,σ(n−1)\sigma(0),\dots,\sigma(n-1)σ(0),…,σ(n−1) with every (σ(i),σ(i+1))∈A(\sigma(i),\sigma(i+1))\in A(σ(i),σ(i+1))∈A.

From GGG the paper builds an instance with nnn jobs and m=n(n−1)+2m=n(n-1)+2m=n(n−1)+2 machines. Each ordered pair (j,k)(j,k)(j,k) of distinct jobs is assigned a middle machine ι(j,k)∈{2,…,m−1}\iota(j,k)\in\{2,\dots,m-1\}ι(j,k)∈{2,…,m−1}, and the partial sums qℓiq_{\ell i}qℓi​ are perturbed from iμi\muiμ by ±λ\pm\lambda±λ or ±(λ+1)\pm(\lambda+1)±(λ+1) on the machines ι(ℓ,⋅)\iota(\ell,\cdot)ι(ℓ,⋅) and just before the machines ι(⋅,ℓ)\iota(\cdot,\ell)ι(⋅,ℓ), depending on whether the pair is an arc. The parameters satisfy λ≥1\lambda\ge1λ≥1 and μ≥2λ+3\mu\ge2\lambda+3μ≥2λ+3.

Formalization targets

Goal (Theorem 5)

DIRECTED HAMILTON PATH ∝ n∣m∣F,no wait∣Cmax⁡andDIRECTED HAMILTON PATH ∝ n∣m∣F,no wait,wj=1∣∑wjCj,\text{DIRECTED HAMILTON PATH}\ \propto\ n|m|F,\textit{no wait}|C_{\max}\quad\text{and}\quad \text{DIRECTED HAMILTON PATH}\ \propto\ n|m|F,\textit{no wait},w_j=1|\textstyle\sum w_jC_j,DIRECTED HAMILTON PATH ∝ n∣m∣F,no wait∣Cmax​andDIRECTED HAMILTON PATH ∝ n∣m∣F,no wait,wj​=1∣∑wj​Cj​,

where ∝\propto∝ is polynomial-time many-one reducibility between the binary-coded recognition languages.

Milestones, in the order the proof uses them

  1. Eq. (9): cjkc_{jk}cjk​ is the least gap between BjB_jBj​ and BkB_kBk​ for which JkJ_kJk​ follows JjJ_jJj​ on every machine.
  2. The travelling-salesman reformulation: if all processing times are positive, a feasible schedule with Cmax⁡≤yC_{\max}\le yCmax​≤y (resp. ∑Cℓ≤y\sum C_\ell\le y∑Cℓ​≤y) exists iff some job order has path length (resp. summed completion times along the path) at most yyy.
  3. An admissible ordering ι\iotaι exists for every n≠2n\ne2n=2.
  4. All processing times of the construction are at least 111.
  5. The delays of the construction: cjk=μ+2λc_{jk}=\mu+2\lambdacjk​=μ+2λ if (j,k)∈A(j,k)\in A(j,k)∈A, and μ+2λ+2\mu+2\lambda+2μ+2λ+2 otherwise.
  6. Theorem 5(a), the equivalence: GGG has a Hamilton path   ⟺  \iff⟺ Cmax⁡≤(n−1)(μ+2λ)+mμC_{\max}\le(n-1)(\mu+2\lambda)+m\muCmax​≤(n−1)(μ+2λ)+mμ.
  7. Theorem 5(b), the equivalence: GGG has a Hamilton path   ⟺  \iff⟺ ∑jCj≤12n(n−1)(μ+2λ)+nmμ\sum_jC_j\le\tfrac12n(n-1)(\mu+2\lambda)+nm\mu∑j​Cj​≤21​n(n−1)(μ+2λ)+nmμ.
  8. Theorem 2(d), off the goal's path: G′G'G′ has a Hamilton circuit iff the graph obtained by splitting a vertex v′v'v′ into v′v'v′ and a new sink v′′v''v′′ has a Hamilton path.

Significance

The result. Theorem 5, together with the NP-completeness of DIRECTED HAMILTON PATH, places both no-wait criteria among the NP-complete problems once mmm is part of the input. The Gilmore–Gomory algorithm for two machines therefore cannot be extended to arbitrary mmm unless P = NP, and heuristics and exact exponential methods for no-wait shops are justified. The reduction also gives a structural fact of independent use: every directed graph can be realized, up to two arc-weight values, as the delay matrix of a no-wait flow shop.

Formalizing it. The result has been proved for fifty years; to our knowledge no machine-checked version exists. This mission produces a checked no-wait flow shop model, a checked travelling-salesman reformulation of it, the correctness of the paper's construction, and a polynomial-time reduction in an explicit Turing-machine model. It also records a gap in the printed proof: the ordering ι\iotaι the paper calls easy to construct does not exist for n=2n=2n=2.

Difficulty

Three steps are not routine.

  1. The passage from schedules to job orders assumes that the order in which jobs start is the order in which they visit every machine. With zero processing times this fails, so the reformulation needs the positivity milestone.
  2. Computing the delays requires a case analysis over all machines and all pairs of perturbations. The property of ι\iotaι is exactly what keeps the "+++" and "−-−" cases from colliding, and the constants λ≥1\lambda\ge1λ≥1, μ≥2λ+3\mu\ge2\lambda+3μ≥2λ+3 are tight enough that the comparison must be done carefully.
  3. The goal asks for an actual Turing machine with a polynomial step bound. It has to build ι\iotaι for every n≠2n\ne2n=2 and handle n≤2n\le2n≤2 separately. The equivalences alone do not give this.

Formalization scope

Jobs and machines are 000-based (Fin n, Fin m); cum p ℓ i is the paper's qℓiq_{\ell i}qℓi​, and the machine indices ι(j,k)\iota(j,k)ι(j,k) are the paper's 111-based indices. Start times are natural numbers, and release dates are 000. Section 3 computes times from processing orders on integer data, and every criterion is regular, so real start times would not change the yes-instances. A zero-length operation occupies the empty interval. Cmax⁡≤yC_{\max}\le yCmax​≤y is stated as Cℓ≤yC_\ell\le yCℓ​≤y for every job. The delays, the partial sums of the construction and its processing times are computed in Z\mathbb ZZ. The instance uses their conversion to N\mathbb NN, which is exact by milestone 4. The threshold of Theorem 5(b) is compared in Q\mathbb QQ, as printed.

The paper's loose phrases are made explicit as follows:

  • "scheduled directly after" means "precedes on every machine";
  • "equivalent to solving the TRAVELLING SALESMAN problem" means the two threshold equivalences of milestone 2, with the ∑Cj\sum C_j∑Cj​ version as the reading of "constructed as in (a)";
  • "such an ordering can easily be constructed" is stated for n≠2n\ne2n=2, the only case in which it is true.

Directed graphs are Boolean adjacency matrices. A graph with 000 or 111 vertices has a Hamilton path, and a one-vertex graph has a Hamilton circuit iff it has a loop.

The goal is CookPvsNP.PolyReducible from the published definition CookPvsNP_defs, with polynomial time measured on Cook's one-tape Turing machines. Instances are coded with the alphabet BSym and binary code encNats of the published ProjSchedTW.Complexity.Encoding (Neumann, Schwindt and Zimmermann). A graph is coded as nnn followed by its adjacency matrix; a flow-shop instance as nnn, mmm, the matrix (pℓi)(p_{\ell i})(pℓi​) and the threshold yyy.

Stating only the equivalences 6–7 and calling the result "reducible" would drop polynomiality; the goal therefore asserts PolyReducible. The target languages contain exactly the no-wait flow-shop instances, with no waiting allowed, so the reduction cannot land in a looser problem. The parameters λ,μ\lambda,\muλ,μ always carry λ≥1\lambda\ge1λ≥1, μ≥2λ+3\mu\ge2\lambda+3μ≥2λ+3.

Contributions are welcome at every level:

  • the general travelling-salesman reformulation, which is reusable for any no-wait flow-shop result;
  • the combinatorial existence of ι\iotaι;
  • the arithmetic of the construction;
  • Turing-machine infrastructure for computing arithmetic list transformations in polynomial time, which every reduction in this series needs.

Selected references

  • P. Brucker, J. K. Lenstra, A. H. G. Rinnooy Kan, Complexity of Machine Scheduling Problems, Mathematisch Centrum Report BW 43/75, Amsterdam, 1975; Annals of Discrete Mathematics 1 (1977) 343–362. doi:10.1016/S0167-5060(08)70743-X
  • R. M. Karp, Reducibility among combinatorial problems, in Complexity of Computer Computations, Plenum, 1972, 85–103. doi:10.1007/978-1-4684-2001-2_9
  • P. C. Gilmore, R. E. Gomory, Sequencing a one state-variable machine: a solvable case of the traveling salesman problem, Operations Research 12 (1964) 655–679. doi:10.1287/opre.12.5.655
  • D. A. Wismer, Solution of the flowshop-scheduling problem with no intermediate queues, Operations Research 20 (1972) 689–697. doi:10.1287/opre.20.3.689
  • S. S. Reddi, C. V. Ramamoorthy, On the flow-shop sequencing problem with no wait in process, Operational Research Quarterly 23 (1972) 323–331. doi:10.1057/jors.1972.52
  • H. Röck, The three-machine no-wait flow shop is NP-complete, Journal of the ACM 31 (1984) 336–345. doi:10.1145/62.65
  • N. G. Hall, C. Sriskandarajah, A survey of machine scheduling problems with blocking and no-wait in process, Operations Research 44 (1996) 510–525. doi:10.1287/opre.44.3.510
  • S. A. Cook, The complexity of theorem-proving procedures, STOC 1971, 151–158. doi:10.1145/800157.805047
14 thms3 active usersReviewed
Dynamic ProgrammingOperations ResearchOptimization·Captain: mikedeng1

Send-and-Split Method for Minimum-Concave-Cost Network Flows II: Subproblem Minimum Costs Are the Greatest Solution of the Send-and-Split Equations, Unique if Circulations Cost Positively (Theorem 2)Research Paper

Motivation

Many network-design and production-planning problems ask for a cheapest way to route a commodity from supply points to demand points when the cost of an arc grows less than proportionally with the amount sent: setup charges, economies of scale, and fixed-plus-linear costs all give concave arc costs. Minimizing a concave function over the flow polytope is NP-hard in general, and the classical linear-cost machinery (potentials, negative-cycle tests) no longer applies. Erickson, Monma and Veinott (Math. Oper. Res. 12 (1987)) gave the send-and-split method, a dynamic program over subsets of the demand nodes that solves the uncapacitated problem in time polynomial in the number of nodes and arcs and exponential only in the number of demand nodes. It contains as special cases the Dreyfus–Wagner recursion for Steiner trees in graphs (Networks 1 (1971)), and the paper applies it to concave-cost production, inventory and network-design problems (§6).

Timeline. Zangwill (1968) solved minimum-concave-cost flows on special networks by dynamic programming. Dreyfus and Wagner (1971) gave the subset recursion for Steiner trees with positive arc lengths. Erickson, Monma and Veinott (1987) extended the subset recursion to arbitrary uncapacitated networks with additive concave arc costs and arbitrary real demands, proved that the minimum costs of the subproblems form the greatest solution of the resulting equations (Theorem 2), and showed it is the only solution when every simple circulation has positive cost.

Setting

A network consists of a directed graph G=(N,A)G = (N, A)G=(N,A) with nodes N={0,…,n−1}N = \{0, \dots, n-1\}N={0,…,n−1} and arcs AAA, ordered pairs of distinct nodes, together with a real demand vector r=(ri)r = (r_i)r=(ri​). A preflow is a nonnegative matrix x=(xij)x = (x_{ij})x=(xij​) supported on AAA; a flow for rrr is a preflow with inflow minus outflow equal to rir_iri​ at every node. Each arc has a cost function cijc_{ij}cij​, concave on [0,∞)[0, \infty)[0,∞) with cij(0)=0c_{ij}(0) = 0cij​(0)=0, and the cost of a preflow is c(x)=∑(i,j)∈Acij(xij)c(x) = \sum_{(i,j)\in A} c_{ij}(x_{ij})c(x)=∑(i,j)∈A​cij​(xij​). A minimum-cost flow is a flow whose cost is at most that of every flow. A simple circulation is a circulation equal to some θ>0\theta > 0θ>0 on the arcs of a simple directed circuit and 000 elsewhere.

The demand nodes are D={i:ri≠0}D = \{ i : r_i \neq 0 \}D={i:ri​=0}. For ∅⊂I⊆D\emptyset \subset I \subseteq D∅⊂I⊆D let rI=∑j∈Irjr_I = \sum_{j \in I} r_jrI​=∑j∈I​rj​. The subproblem i→Ii \to Ii→I keeps the demands of the nodes in III, sets all others to zero, and subtracts rIr_IrI​ at node iii; CiIC_{iI}CiI​ is its minimum cost, +∞+\infty+∞ when it has no flow. Let AI=AA_I = AAI​=A if rI>0r_I > 0rI​>0 and the reversed arcs if rI<0r_I < 0rI​<0, with sending cost cij(rI)c_{ij}(r_I)cij​(rI​), read as cji(−rI)c_{ji}(-r_I)cji​(−rI​) when rI<0r_I < 0rI​<0. The send-and-split equations for an array C′C'C′ are

CiI′=Cj,I∖{j}′ (j∈I)if rI=0,(1)C'_{iI} = C'_{j, I\setminus\{j\}} \ (j \in I) \quad \text{if } r_I = 0, \qquad (1)CiI′​=Cj,I∖{j}′​ (j∈I)if rI​=0,(1) CiI′=min⁡(i,j)∈AI[cij(rI)+CjI′]∧BiI′if rI≠0,(2)C'_{iI} = \min_{(i,j)\in A_I} \bigl[ c_{ij}(r_I) + C'_{jI} \bigr] \wedge B'_{iI} \quad \text{if } r_I \neq 0, \qquad (2)CiI′​=(i,j)∈AI​min​[cij​(rI​)+CjI′​]∧BiI′​if rI​=0,(2)

with BiI′=min⁡∅⊂J⊂I[CiJ′+Ci,I∖J′]B'_{iI} = \min_{\emptyset\subset J\subset I} [C'_{iJ} + C'_{i,I\setminus J}]BiI′​=min∅⊂J⊂I​[CiJ′​+Ci,I∖J′​] for ∣I∣>1|I| > 1∣I∣>1, and BiI′=0B'_{iI} = 0BiI′​=0 if I={i}I = \{i\}I={i}, +∞+\infty+∞ otherwise, for ∣I∣=1|I| = 1∣I∣=1 (3).

Formalization targets

Goal: Theorem 2 (Dynamic-Programming Equations)

If there is a minimum-cost flow for rrr, then

C solves (1)–(3),CjI′≤CjI  for every +∞-or-real solution C′,C \text{ solves (1)–(3)}, \qquad C'_{jI} \le C_{jI} \ \text{ for every } +\infty\text{-or-real solution } C',C solves (1)–(3),CjI′​≤CjI​  for every +∞-or-real solution C′,

for all j∈Nj \in Nj∈N and ∅⊂I⊆D\emptyset \subset I \subseteq D∅⊂I⊆D. If moreover every simple circulation has positive cost, then every +∞+\infty+∞-or-real solution equals CCC.

Milestones

In attack order: the subproblem alternative (a minimum-cost flow or no flow, p. 640); equation (1); inequalities (4), (5), (6) (p. 641); CiI=BiIC_{iI} = B_{iI}CiI​=BiI​ for i∈Ii \in Ii∈I (p. 641); CCC satisfies (1) and (2) (p. 642); and two generic statements about Bellman's equations (2)′ for minimum-cost chains in a graph with nonnegative, respectively positive, simple circuits: the minimum chain costs form the greatest solution and are nondecreasing in the arc costs, and with positive circuits the solution is unique (p. 642).

Significance

Theorem 2 is the correctness theorem of the send-and-split method: it says that solving equations (1)–(3) by increasing ∣I∣|I|∣I∣, with a shortest-chain computation for each III, yields the true subproblem minimum costs, in particular the optimum CiDC_{iD}CiD​ of the original problem. Its uniqueness clause identifies when the equations alone determine the answer, and the counterexample on p. 641 (zero costs on a strongly connected graph) shows that the positivity hypothesis cannot simply be dropped. The same structure underlies the Steiner-tree recursion.

The theorem is proved in the paper; to our knowledge none of it is formalized. This mission supplies Lean definitions of uncapacitated concave-cost flows with extended-real minimum costs and the send-and-split equations, and poses the correctness theorem as a proof goal. The two Bellman milestones are general facts about minimum-cost chains with +∞+\infty+∞-or-real values and are reusable well beyond this paper. Related items on the platform are credited, not reused: the Dreyfus–Wagner Steiner recursion DreyfusWagner.Steiner.steinerLength_recurrence and DreyfusWagner.Steiner.optimal_decomposition (the special case of undirected positive lengths), and BellmanRouting.PolicySpace.routing_equation_unique (uniqueness of the routing equation on a complete graph with positive real times).

Difficulty

The equations (2) mix two minima: sending the whole demand rIr_IrI​ across one arc, and splitting III at the current node. Showing that CCC satisfies them requires that some optimal flow of every subproblem has a tree-like structure, which in turn rests on the existence theory for concave-cost flows (an optimum is attained at an extreme flow whose support is a forest). Showing that CCC is the greatest solution requires an induction on ∣I∣|I|∣I∣ in which, for fixed III, equation (2) is a shortest-chain system whose arc costs to an auxiliary node depend on the solution at smaller sets; comparison then needs monotonicity of minimum chain costs in the arc costs, including arcs whose cost drops from +∞+\infty+∞. The naive attempt to prove uniqueness by iterating (2) fails without the positivity hypothesis, because zero-cost circuits admit spurious solutions.

Formalization scope

Nodes are Fin n; arcs are a Finset (Fin n × Fin n) with no loops; arc costs are functions R→R\mathbb{R} \to \mathbb{R}R→R assumed concave on Set.Ici 0 with cij(0)=0c_{ij}(0) = 0cij​(0)=0 at arcs (the paper's standing assumption, stated as a hypothesis). Minimum costs, sending costs, splitting terms and solutions take values in EReal; the empty infimum is +∞+\infty+∞, as on the page.

Explicit readings of the paper's phrases:

  • "there is a minimum-cost flow" means a flow for rrr whose cost is at most that of every flow for rrr;
  • "+∞+\infty+∞ or real-valued" means ≠−∞\neq -\infty=−∞, required on every ∅⊂I⊆D\emptyset \subset I \subseteq D∅⊂I⊆D;
  • "greatest solution" is two claims: CCC is a solution, and every solution is ≤C\le C≤C entrywise;
  • "each simple circulation has positive cost" quantifies over every simple circuit and every θ>0\theta > 0θ>0;
  • "CiIC_{iI}CiI​ is finite" is rendered as CiI=c(x)C_{iI} = c(x)CiI​=c(x) for a minimum-cost flow xxx of the subproblem;
  • "minimum-cost chain" is the infimum of walk costs, attained or +∞+\infty+∞; arc costs of +∞+\infty+∞ in (2)′ are missing arcs.

The definition of CiIC_{iI}CiI​ is the infimum over flows, never the equations themselves; a formalization defining CCC by (1)–(3), or omitting the ≠−∞\neq -\infty=−∞ condition on solutions, or dropping the clause that CCC is itself a solution, would make the theorem trivial or false and is ruled out. Arrays are compared only on admissible sets. The graph structure of (2)′ is the platform definition BertsekasSPGraph. Running-time bounds, the planar refinements of §4 (Theorem 3) and the capacitated reduction of §5 are out of scope. Proofs of any milestone, and a formalization of the paper's Theorem 1 that the proofs rely on, are welcome.

Selected references

  • R. E. Erickson, C. L. Monma, A. F. Veinott, Jr., Send-and-Split Method for Minimum-Concave-Cost Network Flows, Mathematics of Operations Research 12(4), 1987, 634–664. https://doi.org/10.1287/moor.12.4.634
  • S. E. Dreyfus, R. A. Wagner, The Steiner Problem in Graphs, Networks 1(3), 1971, 195–207. https://doi.org/10.1002/net.3230010302
  • R. Bellman, On a Routing Problem, Quarterly of Applied Mathematics 16(1), 1958, 87–90. https://doi.org/10.1090/qam/102435
  • W. M. Hirsch, A. J. Hoffman, Extreme Varieties, Concave Functions, and the Fixed Charge Problem, Communications on Pure and Applied Mathematics 14(3), 1961, 355–369. https://doi.org/10.1002/cpa.3160140312
  • W. I. Zangwill, Minimum Concave Cost Flows in Certain Networks, Management Science 14(7), 1968, 429–450. https://doi.org/10.1287/mnsc.14.7.429
15 thms3 active usersReviewed
CombinatoricsOperations Research·Captain: mikedeng1

The Strong Perfect Graph Theorem I: A Graph Is Perfect If and Only If It Is BergeResearch Paper

Motivation

A perfect graph is one whose coloring problem has a particularly sharp answer on every induced subgraph: the fewest colors needed is exactly the size of its largest clique. This makes a local obstruction, a clique, certify the optimum number of colors throughout the graph. Claude Berge proposed in 1961 that perfection could be recognized by the absence of two kinds of induced odd cycles, one in the graph and one in its complement. The equivalence became known as the strong perfect graph conjecture. Chudnovsky, Robertson, Seymour, and Thomas proved it in their 2006 paper, which also proves a structural decomposition of the graphs under study. The paper connects this question to graph coloring, Shannon capacity, and linear and integer programming. Chudnovsky et al., pp. 51–54

The earlier complement theorem was proved by Lovász in 1972 and appears as Theorem 1.1 in the paper. The strong conjecture remained unresolved for roughly four decades; the authors' Theorem 1.2 settles it. Their proof places a second result, Theorem 1.3, beside the equivalence: a graph with no forbidden odd hole or antihole must either belong to a basic class or admit one of several specified decompositions. The graph classes and decompositions are therefore part of the statement of the route to the main result, not merely vocabulary for a proof. Chudnovsky et al., pp. 52–56

Setting

All graphs here are finite and simple. The complement G‾\overline GG has the same vertices as GGG, and two distinct vertices are adjacent in G‾\overline GG exactly when they are not adjacent in GGG. For a vertex set XXX, the notation G∣XG|XG∣X means the induced subgraph on XXX. A clique is a set of pairwise adjacent vertices. Its largest possible size in a graph HHH is ω(H)\omega(H)ω(H), and χ(H)\chi(H)χ(H) is the minimum number of colors in a proper vertex coloring of HHH.

A hole is an induced cycle of length at least four. An antihole of GGG is a hole in G‾\overline GG. A graph is Berge if every hole and antihole has even length. Thus a perfect graph requires χ(G∣X)=ω(G∣X)\chi(G|X)=\omega(G|X)χ(G∣X)=ω(G∣X) for every X⊆V(G)X\subseteq V(G)X⊆V(G), while a Berge graph satisfies a restriction on induced cycles in both GGG and G‾\overline GG. “Induced” matters: a cycle with a chord is not a hole. Chudnovsky et al., pp. 51–52

For the structural milestones, a basic graph is a bipartite graph, the complement of one, a line graph of a bipartite graph, the complement of such a line graph, or a double split graph. The latter consists of paired vertices ai,bia_i,b_iai​,bi​ and cj,djc_j,d_jcj​,dj​ with the within-pair and between-pair adjacencies specified on pp. 52–53. A proper 2-join partitions the vertices into two sides with two prescribed complete cross-edge blocks, connected-component conditions on both sides, and a special odd-path condition. A proper homogeneous pair is a pair of vertex sets whose outside vertices split into four nonempty adjacency classes. A balanced skew partition has one side disconnected and the other disconnected in the complement, together with parity restrictions on induced paths and antipaths. Chudnovsky et al., pp. 52–54

Formalization targets

Theorem 1.2: perfection and the Berge property

The goal is the exact equivalence for every finite simple graph:

G is perfect⟺G is Berge.G\text{ is perfect}\quad\Longleftrightarrow\quad G\text{ is Berge}.G is perfect⟺G is Berge.

No order bound, chosen graph class, or decomposition hypothesis is attached to the goal. Chudnovsky et al., p. 52, 1.2

Structural and reduction milestones

Theorem 1.1 says that GGG perfect implies G‾\overline GG perfect. Theorem 1.5 says that a minimum imperfect graph, a Berge nonperfect graph with the smallest vertex count among all such graphs, cannot admit a balanced skew partition. Theorem 13.5 says that a recalcitrant graph—a Berge graph with the listed line-graph, double-split, 2-join, homogeneous-pair, and balanced-skew outcomes absent—has GGG or G‾\overline GG bipartite. Theorem 1.3 states the decomposition conclusion:

G Berge⟹G basic ∨ G or G‾ has a proper 2-join ∨ G has a proper homogeneous pair ∨ G has a balanced skew partition.G\text{ Berge}\Longrightarrow G\text{ basic}\ \lor\ G\text{ or }\overline G\text{ has a proper 2-join}\ \lor\ G\text{ has a proper homogeneous pair}\ \lor\ G\text{ has a balanced skew partition}.G Berge⟹G basic ∨ G or G has a proper 2-join ∨ G has a proper homogeneous pair ∨ G has a balanced skew partition.

The milestone order records the two reduction results, the later structural capstone, and the decomposition statement it yields. The paper's other section results that establish 13.5 are posed in the remaining missions of this series. Chudnovsky et al., pp. 52, 54–55, 154

Significance

Theorem 1.2 gives a forbidden-induced-subgraph characterization of perfect graphs. Its cycle condition is intrinsic to the graph and its complement; its coloring condition quantifies over every induced subgraph. Together with Theorem 1.3, it ties a numerical property of colorings to explicit graph structures and separations. The complement theorem and the exclusion of decompositions for a minimum imperfect graph explain why the structural alternatives have the strength needed for the equivalence. Chudnovsky et al., pp. 52–56

The mathematical theorem was proved in the cited paper. This mission poses its statements in Lean and seeks machine-checked proofs; the draft theorem declarations are open targets. The definition layer is useful beyond this mission: induced holes, Berge graphs, perfection, balanced skew partitions, and the decomposition predicates can support the later missions without changing what each source statement means. No machine-checked proof of these draft targets is claimed here.

Difficulty

The forward implication can be tested on induced odd cycles, but that observation does not settle the converse. Excluding odd holes and antiholes does not give an immediate coloring of an arbitrary induced subgraph. The paper instead establishes a detailed account of what a Berge graph can look like when it is not in a basic class. The delicate point in turning this account into Theorem 1.2 is that each decomposition outcome must be incompatible with a minimum imperfect graph. Ordinary skew partitions are too broad for that role; the balanced parity conditions are part of the statement. The structural conclusion 13.5 collects restrictions established across many later sections, so formalizing its prerequisites is a substantial graph-theoretic task. Chudnovsky et al., pp. 52–56, 154

Formalization scope

Lean uses SimpleGraph V with finite vertices, decidable vertex equality, and the Mathlib complement, induced subgraph, chromatic number, clique number, bipartiteness, and line graph. A hole is a list in cyclic order whose adjacency relation agrees exactly with the cycle edges; a path is likewise listed in one orientation with exactly its consecutive edges. Antiholes and antipaths use the complement graph. The empty vertex set is connected, matching p. 54. The double split partition is encoded by an equivalence from the four indexed parts to the whole vertex type, making disjointness and coverage explicit. A minimum imperfect graph is globally minimal by vertex count, among all finite Berge graphs.

No hypothesis beyond the paper's finite, simple graph convention is added to the goal or numbered milestones. In particular, “Berge” includes holes in both GGG and G‾\overline GG; “perfect” ranges over every induced subgraph; and the complement occurs only in those decomposition outcomes where the paper places it. A non-induced cycle or a missing complement condition would make the formal target different. The definitions and Lean proofs of the graph classes, their boundary cases, and the structural milestones are welcome contributions. Later missions pose the paper's intervening numbered results rather than duplicating them here.

Selected references

  • Maria Chudnovsky, Neil Robertson, Paul Seymour, and Robin Thomas, The strong perfect graph theorem, Annals of Mathematics 164 (2006), 51–229. DOI 10.4007/annals.2006.164.51.
16 thms3 active usersReviewed
CombinatoricsLinear algebra·Captain: mikedeng1

On a Conjecture of Spectral Extremal Problems: If the Extremal Graphs for F Are Turán Graphs Plus a Fixed Number of Edges, Every F-Free Graph of Maximum Spectral Radius Is ExtremalResearch Paper

Motivation

Extremal graph theory asks how many edges a graph on nnn vertices can have without containing a fixed graph FFF. The answer, the Turán number ex(n,F)\mathrm{ex}(n,F)ex(n,F), and the graphs attaining it, the set Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) of extremal graphs, are known precisely only for special FFF; the Erdős–Stone–Simonovits theorem gives ex(n,F)=(1−1χ(F)−1+o(1))n22\mathrm{ex}(n,F) = (1 - \frac{1}{\chi(F)-1} + o(1))\frac{n^2}{2}ex(n,F)=(1−χ(F)−11​+o(1))2n2​, where χ(F)\chi(F)χ(F) is the chromatic number.

Spectral extremal graph theory asks the same question with the number of edges replaced by the spectral radius λ(G)\lambda(G)λ(G), the largest eigenvalue of the adjacency matrix. Since λ(G)≥2e(G)/n\lambda(G) \ge 2e(G)/nλ(G)≥2e(G)/n, a spectral bound implies an edge bound, and spectral extremal results are usually stronger than their edge versions. Nikiforov (Linear Algebra Appl. 427, 2007) showed that the Turán graph Tn,rT_{n,r}Tn,r​ maximises λ\lambdaλ among Kr+1K_{r+1}Kr+1​-free graphs, so for F=Kr+1F = K_{r+1}F=Kr+1​ the spectral and the edge extremal graphs coincide. Whether this happens for other FFF is the subject of the paper.

Timeline:

  • 1941 Turán: Tn,rT_{n,r}Tn,r​ is the unique extremal graph for Kr+1K_{r+1}Kr+1​.
  • 2003 Chen, Gould, Pfender, Wei (J. Combin. Theory Ser. B 89): ex(n,Fk,r+1)=e(Tn,r)+O(1)\mathrm{ex}(n, F_{k,r+1}) = e(T_{n,r}) + O(1)ex(n,Fk,r+1​)=e(Tn,r​)+O(1) for kkk copies of Kr+1K_{r+1}Kr+1​ sharing one vertex.
  • 2007 Nikiforov (Linear Algebra Appl. 427): spectral Turán theorem for Kr+1K_{r+1}Kr+1​. 2009 Nikiforov (J. Graph Theory 62): spectral stability for large forbidden subgraphs, the source of Lemma 2.5.
  • 2020 Cioabă, Feng, Tait, Zhang (Electron. J. Combin. 27): the spectral extremal graph for the friendship graph FkF_kFk​ lies in Ex(n,Fk)\mathrm{Ex}(n, F_k)Ex(n,Fk​).
  • 2022 Cioabă, Desai, Tait (European J. Combin. 99) conjecture: if the graphs in Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) are Turán graphs plus O(1)O(1)O(1) edges, then the spectral extremal graphs lie in Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) for large nnn. Before the general proof it was known for Kr+1K_{r+1}Kr+1​, the friendship graphs FkF_kFk​, the graphs Hs,kH_{s,k}Hs,k​ (Li, Peng) and the intersecting cliques Fk,rF_{k,r}Fk,r​ (Desai, Kang, Li, Ni, Tait, Wang, arXiv:2108.03587).
  • 2022/2023 Wang, Kang, Xue (arXiv:2203.10831; J. Combin. Theory Ser. B, 2023): the conjecture holds in general. This mission formalizes their Theorem 1.2.

Setting

All graphs are finite and simple. For a graph GGG on nnn vertices, A(G)A(G)A(G) is its 0/10/10/1 adjacency matrix and λ(G)\lambda(G)λ(G) is the largest eigenvalue of A(G)A(G)A(G). A graph GGG is FFF-free if no subgraph of GGG is isomorphic to FFF. The Turán number ex(n,F)\mathrm{ex}(n,F)ex(n,F) is the maximum number of edges e(G)e(G)e(G) over FFF-free graphs GGG on nnn vertices, and Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) is the set of FFF-free nnn-vertex graphs with ex(n,F)\mathrm{ex}(n,F)ex(n,F) edges. The Turán graph Tn,rT_{n,r}Tn,r​ is the complete rrr-partite graph on nnn vertices with parts of sizes ⌊n/r⌋\lfloor n/r\rfloor⌊n/r⌋ or ⌈n/r⌉\lceil n/r\rceil⌈n/r⌉.

The hypothesis on FFF is: for fixed integers r≥2r \ge 2r≥2 and a≥0a \ge 0a≥0 and all large nnn, ex(n,F)=e(Tn,r)+a\mathrm{ex}(n,F) = e(T_{n,r}) + aex(n,F)=e(Tn,r​)+a and every graph in Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) contains a spanning copy of Tn,rT_{n,r}Tn,r​, i.e. is Tn,rT_{n,r}Tn,r​ plus aaa edges. The paper states "adding O(1)O(1)O(1) edges" and fixes the constant at the start of Section 3 (p. 4): "We may assume that the graphs in Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) are obtained from Tn,rT_{n,r}Tn,r​ by adding aaa edges." The mission follows that reading. Examples: Kr+1K_{r+1}Kr+1​ with a=0a = 0a=0, the friendship graphs, and the intersecting cliques Fk,r+1F_{k,r+1}Fk,r+1​.

A graph GGG is spectral extremal for FFF if it is FFF-free and λ(G)≥λ(G′)\lambda(G) \ge \lambda(G')λ(G)≥λ(G′) for every FFF-free G′G'G′ on the same nnn vertices.

Formalization targets

Goal: Theorem 1.2

For r≥2r \ge 2r≥2, a≥0a \ge 0a≥0 and FFF satisfying the hypothesis, there is NNN such that for all n≥Nn \ge Nn≥N every spectral extremal graph GGG for FFF on nnn vertices satisfies

e(G)=ex(n,F),i.e.G∈Ex(n,F).e(G) = \mathrm{ex}(n,F), \qquad\text{i.e.}\qquad G \in \mathrm{Ex}(n,F).e(G)=ex(n,F),i.e.G∈Ex(n,F).

The statement contains no numerical constant, and NNN depends only on FFF, rrr, aaa.

Milestones

The milestones follow the proof, in order: strict monotonicity of λ\lambdaλ under proper subgraphs of a connected graph (Lemma 2.3); connectivity of GGG (Lemma 3.1); the bound λ(G)≥(1−1r)n−r4n+2an\lambda(G) \ge (1-\frac1r)n - \frac{r}{4n} + \frac{2a}{n}λ(G)≥(1−r1​)n−4nr​+n2a​ (Lemma 3.2); spectral stability for χ(F)=r+1\chi(F) = r+1χ(F)=r+1 (Corollary 2.6); for every maximum rrr-cut V1∪⋯∪VrV_1 \cup \dots \cup V_rV1​∪⋯∪Vr​, ∑ie(Vi)≤εn2\sum_i e(V_i) \le \varepsilon n^2∑i​e(Vi​)≤εn2 and ∣Vi∣=(1r±3ε)n|V_i| = (\frac1r \pm 3\sqrt\varepsilon)n∣Vi​∣=(r1​±3ε​)n (Lemma 3.3); a counting inequality for intersections (Lemma 2.8); e(G[Vi])≤ae(G[V_i]) \le ae(G[Vi​])≤a and minimum degree above (1−1r−3rε1/3)n(1 - \frac1r - 3r\varepsilon^{1/3})n(1−r1​−3rε1/3)n (Lemma 3.6); at most 2a2a2a vertices of each part have a neighbour in that part, and all others see every other part completely (Lemma 3.7); Perron entries xu≥1−20a2r2/nx_u \ge 1 - 20a^2r^2/nxu​≥1−20a2r2/n for a≥1a \ge 1a≥1 (Lemma 3.8); e(Gin)−e(Gout)≤ae(G_{in}) - e(G_{out}) \le ae(Gin​)−e(Gout​)≤a (Lemma 3.9); balancing two parts of a complete multipartite graph increases λ\lambdaλ (Lemma 2.7); and the maximum partition is balanced, ∣ni−nj∣≤1|n_i - n_j| \le 1∣ni​−nj​∣≤1 (Lemma 3.10).

Significance

The theorem settles the Cioabă–Desai–Tait conjecture: for every FFF whose extremal graphs are Turán graphs plus a bounded number of edges, the spectral extremal problem reduces to the edge extremal problem for large nnn. This recovers the earlier cases (friendship graphs, the graphs Hs,kH_{s,k}Hs,k​, intersecting cliques) at once, and it turns any future determination of Ex(n,F)\mathrm{Ex}(n,F)Ex(n,F) of this type into a spectral result with no further work.

The result is proved on paper; it has no machine-checked proof. Mathlib has Turán's theorem (extremalNumber_top, uniqueness of turanGraph) and the definition of extremalNumber, but no spectral extremal graph theory: no Perron–Frobenius theorem for graphs, no spectral Turán theorem, no stability theorem. The formalization would supply these, and the milestone statements are independently reusable: strict monotonicity of λ\lambdaλ (Lemma 2.3), the spectral comparison of complete multipartite graphs (Lemma 2.7) and spectral stability (Corollary 2.6) are standard tools of the area.

Difficulty

The natural first idea, comparing GGG with an extremal graph HHH by Rayleigh quotients, fails at the start: λ(G)≥λ(H)\lambda(G) \ge \lambda(H)λ(G)≥λ(H) gives only e(G)≥e(Tn,r)−o(n2)e(G) \ge e(T_{n,r}) - o(n^2)e(G)≥e(Tn,r​)−o(n2), far from ex(n,F)\mathrm{ex}(n,F)ex(n,F). Closing an additive gap of O(1)O(1)O(1) edges requires control of the Perron vector to within O(1/n)O(1/n)O(1/n) at every vertex and exact control of the part sizes. The second is the delicate step: an imbalance of one vertex between two parts costs Θ(1/n)\Theta(1/n)Θ(1/n) in λ\lambdaλ, while the aaa extra edges contribute only O(1/n2)O(1/n^2)O(1/n2) beyond the Turán graph, so the two effects must be compared at different scales. Further, the proof needs the deep spectral stability theorem of Nikiforov (Lemma 2.5), whose own proof is long.

Formalization scope

Graphs on nnn vertices are SimpleGraph (Fin n), matching Mathlib's extremalNumber n F; FFF is a graph on any finite type. λ(G)\lambda(G)λ(G) is the largest eigenvalue of G.adjMatrix ℝ (index 000 of Mathlib's decreasingly sorted eigenvalues₀), not an absolute value and not a norm. "FFF-free" is F.Free G (no copy of FFF, not necessarily induced). "Sufficiently large nnn" is ∃ N, ∀ n ≥ N with NNN chosen after F,r,aF, r, aF,r,a and before GGG; "sufficiently small ε\varepsilonε" is ∃ ε₀ > 0, ∀ ε ∈ (0, ε₀). Partitions V1∪⋯∪VrV_1 \cup \dots \cup V_rV1​∪⋯∪Vr​ are labellings Fin n → Fin r whose parts may be empty; the Section 3 lemmas hold for every partition maximising the number of crossing edges. Lemma 3.8 carries the extra hypothesis a≥1a \ge 1a≥1, because as printed it is false for a=0a = 0a=0 (for F=Kr+1F = K_{r+1}F=Kr+1​ and r∤nr \nmid nr∤n the Perron vector of Tn,rT_{n,r}Tn,r​ has entries below 111); the goal does not assume it.

Trivializing formalizations are ruled out: λ\lambdaλ is not defined from the edge count (which would make spectral and edge extremality the same), the hypothesis on FFF does not contain the conclusion and is satisfiable (a sorry-free check for F=Kr+1F = K_{r+1}F=Kr+1​, a=0a = 0a=0 was compiled), the conclusion is e(G)=ex(n,F)e(G) = \mathrm{ex}(n,F)e(G)=ex(n,F) and not a weaker bound, and the threshold is not chosen after GGG.

A complete development needs the Perron–Frobenius theorem for irreducible nonnegative symmetric matrices, the Rayleigh quotient characterisation of λ\lambdaλ, the spectrum of complete multipartite graphs, Nikiforov's spectral stability lemma, and max-cut partition arguments. All of these are reusable beyond this mission, and proofs of any milestone, of the cited Lemmas 2.1, 2.2 and 2.5, or of Nikiforov's spectral Turán theorem are welcome.

Selected references

  • J. Wang, L. Kang, Y. Xue, On a conjecture of spectral extremal problems, J. Combin. Theory Ser. B, 2023; arXiv:2203.10831v1 (2022). https://arxiv.org/abs/2203.10831
  • S. Cioabă, D. N. Desai, M. Tait, The spectral radius of graphs with no odd wheels, European J. Combin. 99 (2022) 103420.
  • S. Cioabă, L. H. Feng, M. Tait, X. D. Zhang, The maximum spectral radius of graphs without friendship subgraphs, Electron. J. Combin. 27(4) (2020) P4.22.
  • G. Chen, R. J. Gould, F. Pfender, B. Wei, Extremal graphs for intersecting cliques, J. Combin. Theory Ser. B 89 (2003) 159–171.
  • V. Nikiforov, Bounds on graph eigenvalues II, Linear Algebra Appl. 427 (2007) 183–189.
  • V. Nikiforov, Stability for large forbidden subgraphs, J. Graph Theory 62(4) (2009) 362–368.
  • D. N. Desai, L. Kang, Y. Li, Z. Ni, M. Tait, J. Wang, Spectral extremal graphs for intersecting cliques, arXiv:2108.03587v2 (2021). https://arxiv.org/abs/2108.03587
18 thms3 active usersReviewed
CombinatoricsLinear OptimizationOperations Research·Captain: mikedeng1

The Matroids with the Max-Flow Min-Cut Property: Binary Mengerian Clutters and the Q6 MinorResearch Paper

Motivation

Several classical theorems of combinatorial optimization say that a family of sets arising from a graph packs: the maximum number of pairwise disjoint members equals the minimum size of a set meeting every member. König's theorem on bipartite graphs, Menger's theorem, the max-flow min-cut theorem of Ford and Fulkerson, Edmonds' branching theorem and the Lucchesi–Younger theorem all have this form (Seymour 1977, (1.1)–(1.5)). In the capacitated version (weights on elements, integral flows) the max-flow min-cut theorem says more: the packing property survives every deletion and replication of elements. Clutters with this stronger property are called Mengerian. For each 000–111 matrix they are exactly the systems whose covering linear program and its dual have integral optima for every integral weight vector, which is why the notion matters to integer programming and polyhedral combinatorics.

Seymour's paper answers the question for the class of binary clutters, the clutters coming from binary matroids, which includes path collections, cut collections and odd-circuit collections of graphs. Earlier, Gallai's theorem implied that ports of regular matroids are Mengerian (Seymour 1977, p. 200); combined with Tutte's excluded-minor characterization of regular matroids, this showed that binary clutters without Q6Q_6Q6​ or b(Q6)b(Q_6)b(Q6​) minors are Mengerian. Seymour shows that the second excluded minor is unnecessary, so a single small clutter is the only obstruction.

Setting

All sets are finite. A clutter L\mathbf LL is a finite collection of finite sets, no member of which is contained in another; ∅\emptyset∅ and {∅}\{\emptyset\}{∅} are the two trivial clutters. Its ground set is E(L)=⋃A∈LAE(\mathbf L)=\bigcup_{A\in\mathbf L}AE(L)=⋃A∈L​A. The blocker b(L)b(\mathbf L)b(L) is the collection of minimal subsets of E(L)E(\mathbf L)E(L) that meet every member of L\mathbf LL, and τ(L)\tau(\mathbf L)τ(L) is the minimum cardinality of a member of b(L)b(\mathbf L)b(L).

L\mathbf LL is Mengerian if L={∅}\mathbf L=\{\emptyset\}L={∅}, or if for every weight map w:E(L)→Z+w:E(\mathbf L)\to\mathbb Z^+w:E(L)→Z+ there is an integral packing q:L→Z+q:\mathbf L\to\mathbb Z^+q:L→Z+ with ∑A∋xq(A)≤w(x)\sum_{A\ni x}q(A)\le w(x)∑A∋x​q(A)≤w(x) for each x∈E(L)x\in E(\mathbf L)x∈E(L) and

∑A∈Lq(A)=min⁡B∈b(L)∑x∈Bw(x).\sum_{A\in\mathbf L}q(A)=\min_{B\in b(\mathbf L)}\sum_{x\in B}w(x).A∈L∑​q(A)=B∈b(L)min​x∈B∑​w(x).

For a set ZZZ, the deletion is L∖Z={A∈L:A∩Z=∅}\mathbf L\setminus Z=\{A\in\mathbf L:A\cap Z=\emptyset\}L∖Z={A∈L:A∩Z=∅} and the contraction L/Z\mathbf L/ZL/Z is the collection of minimal members of {A−Z:A∈L}\{A-Z:A\in\mathbf L\}{A−Z:A∈L} (minimal, not minimal nonempty). A minor of L\mathbf LL is any clutter obtained by a finite sequence of deletions and contractions.

A clutter is binary if ∣A∩B∣|A\cap B|∣A∩B∣ is odd for all A∈LA\in\mathbf LA∈L and B∈b(L)B\in b(\mathbf L)B∈b(L); this is condition (3.2)(ii) of the paper, which is equivalent to being a port of a binary matroid. Finally

Q6={{1,3,5},{1,4,6},{2,3,6},{2,4,5}},Q_6=\{\{1,3,5\},\{1,4,6\},\{2,3,6\},\{2,4,5\}\},Q6​={{1,3,5},{1,4,6},{2,3,6},{2,4,5}},

the triangles of K4K_4K4​ with its edges labelled 1,…,61,\dots,61,…,6.

For the structure theory, a circuit of a binary clutter is a minimal nonempty C⊆E(L)C\subseteq E(\mathbf L)C⊆E(L) with ∣C∩B∣|C\cap B|∣C∩B∣ even for every B∈b(L)B\in b(\mathbf L)B∈b(L); xxx and yyy are parallel when {x,y}\{x,y\}{x,y} is a circuit, and the point ⟨x⟩\langle x\rangle⟨x⟩ is the parallel class of xxx. With mb(L)={B∈b(L):∣B∣=τ(L)}mb(\mathbf L)=\{B\in b(\mathbf L):|B|=\tau(\mathbf L)\}mb(L)={B∈b(L):∣B∣=τ(L)}, L\mathbf LL is critical if E(mb(L))=E(L)E(mb(\mathbf L))=E(\mathbf L)E(mb(L))=E(L). In a critical binary clutter, x→yx\to yx→y means that every member of mb(L)mb(\mathbf L)mb(L) containing xxx contains yyy while y∉⟨x⟩y\notin\langle x\rangley∈/⟨x⟩, and yyy is initial if no xxx has x→yx\to yx→y. MBC abbreviates "Mengerian binary clutter".

Formalization targets

Goal: Seymour's theorem (p. 209)

For every binary clutter L\mathbf LL,

L is Mengerian  ⟺  L has no minor isomorphic to Q6.\mathbf L\ \text{is Mengerian}\iff \mathbf L\ \text{has no minor isomorphic to } Q_6 .L is Mengerian⟺L has no minor isomorphic to Q6​.

Milestones

In the order the proof uses them:

  • (2.3) Every minor of a Mengerian clutter is Mengerian.
  • Section 1, p. 193. Q6Q_6Q6​ is not Mengerian. With (2.3) this is the "only if" direction.
  • (3.6)(i) Circuits of a binary clutter have at least two elements.
  • (3.6)(iii) If Z⊆E(L)Z\subseteq E(\mathbf L)Z⊆E(L) meets every member of b(L)b(\mathbf L)b(L) evenly, then ZZZ is a disjoint union of circuits. If it meets every member oddly, then ZZZ is a disjoint union of circuits and one member of L\mathbf LL.
  • (4.3) In a critical MBC, x→yx\to yx→y implies y↛xy\not\to xy→x.
  • (4.4) In a critical MBC, x→yx\to yx→y gives a circuit C∋x,yC\ni x,yC∋x,y with ∣C∣≥3|C|\ge3∣C∣≥3, z→yz\to yz→y for z∈C−{y}z\in C-\{y\}z∈C−{y}, and ∣B−(C−{y})∣≥τ(L)−1|B-(C-\{y\})|\ge\tau(\mathbf L)-1∣B−(C−{y})∣≥τ(L)−1 for B∈b(L)B\in b(\mathbf L)B∈b(L).
  • (4.5) In a critical MBC, a non-initial xxx lies on a circuit CCC with ∣C∣≥3|C|\ge3∣C∣≥3 whose other elements are initial and point to xxx, and ∣B∩(C−{x})∣≤1|B\cap(C-\{x\})|\le1∣B∩(C−{x})∣≤1 for B∈mb(L)B\in mb(\mathbf L)B∈mb(L).
  • (4.6) A nontrivial critical MBC has a member consisting of initial elements.
  • (5.1) A binary clutter with six elements x1,y1,x2,y2,x3,y3x_1,y_1,x_2,y_2,x_3,y_3x1​,y1​,x2​,y2​,x3​,y3​ whose only circuits are the three sets {xi,yi,xj,yj}\{x_i,y_i,x_j,y_j\}{xi​,yi​,xj​,yj​}, together with a member AAA that meets each pair {xi,yi}\{x_i,y_i\}{xi​,yi​} once and satisfies a minimality condition, has a Q6Q_6Q6​ minor.

Significance

The theorem is an excluded-minor characterization of the max-flow min-cut property. For binary clutters it decides exactly when the covering system Mx≥1Mx\ge1Mx≥1, x≥0x\ge0x≥0 has integral optimal primal and dual solutions for every integral cost vector, and it identifies Q6Q_6Q6​ as the single obstruction. Its matroid form (the Corollary, p. 220) states that for a matroid MMM the port Ω(M)\Omega(M)Ω(M) is Mengerian for every element Ω\OmegaΩ if and only if MMM is binary and has no F7∗F_7^*F7∗​ minor. Consequences discussed in the paper include the two-commodity setting of (3.5): the clutter of minimal edge sets joining sss to s′s's′ or ttt to t′t't′ is Mengerian exactly when the graph does not reduce to the configuration of its Figure 2. The theorem is also a basis for later work on ideal and Mengerian clutters, such as Cornuéjols' book Combinatorial Optimization: Packing and Covering (SIAM, 2001).

The result has been proved since 1977. To our knowledge no machine-checked proof exists. Mathlib at the pinned revision has matroids but no clutters, blockers, clutter minors, or matroids representable over GF(2). This mission builds that layer. The minor-closedness of the Mengerian property (2.3), the parity decomposition (3.6)(iii) and the structure theory of critical Mengerian binary clutters (4.3)–(4.6) are results in their own right and are useful beyond the main theorem.

Difficulty

The "only if" direction is short: minors of Mengerian clutters are Mengerian, and Q6Q_6Q6​ fails with unit weights. The "if" direction is, in the author's words, "very much harder". A natural first idea is to show directly, by LP duality, that the covering polyhedron of a Q6Q_6Q6​-free binary clutter is integral. This does not work: integrality of the polyhedron is the weak max-flow min-cut property, and Q6Q_6Q6​ itself has that property while not being Mengerian, so no argument that sees only fractional optima can separate the two cases. The paper's proof works with a minimal counterexample and derives the Q6Q_6Q6​ minor from the structure of critical Mengerian binary clutters in Section 4; its intermediate claims (5.2)–(5.39) hold only for that minimal counterexample, which is why they are not milestones here.

Formalization scope

Elements form a type α with decidable equality. A clutter is L : Finset (Finset α) with the clutter axiom as a hypothesis, E(L)E(\mathbf L)E(L) is the union of members, and deletion and contraction take an arbitrary finite set ZZZ. Weights www and packings qqq are N\mathbb NN-valued. The minimum in the Mengerian condition is expressed as "some B∈b(L)B\in b(\mathbf L)B∈b(L) of least weight has weight equal to the packing value", never as an infimum. {∅}\{\emptyset\}{∅} is Mengerian by the paper's convention, and τ({∅})\tau(\{\emptyset\})τ({∅}), which the paper leaves undefined, has the junk value 000 in Lean; every item reading τ\tauτ excludes {∅}\{\emptyset\}{∅} or is vacuous there. "Minor" is the reflexive–transitive closure of single deletions and contractions. "Has a Q6Q_6Q6​ minor" means that some minor equals the image of Q6Q_6Q6​ (on Fin 6, with the paper's labels shifted down by one) under an injective relabelling Fin 6 ↪ α. Binary clutters are defined by (3.2)(ii); the paper defines them as ports of binary matroids and quotes (3.2) [15, 28] for the equivalence, and Mathlib has no GF(2)-representable matroids at this revision. Circuits are defined intrinsically, which makes (3.6)(ii) hold by definition.

Four readings would change the theorem and are ruled out: real-valued packings qqq (the weak max-flow min-cut property, which Q6Q_6Q6​ has, so the goal would be false), a non-minimal blocker or one not restricted to E(L)E(\mathbf L)E(L), dropping the {∅}\{\emptyset\}{∅} exception, and reading "Q6Q_6Q6​ minor" as literal equality instead of isomorphism.

A complete development needs the blocker calculus ((2.1), (2.2), cited from [28] with proofs omitted), the parity theory of binary clutters, and the replication operation Lw\mathbf L_wLw​. The clutter layer (blocker, minors, Mengerian, binary, circuits) is reusable for later work on ideal clutters, Lehman's theorem and the Corollary's matroid form. Proofs of any milestone, of the helper facts b(b(L))=Lb(b(\mathbf L))=\mathbf Lb(b(L))=L, (2.1) and (2.2), and of the equivalences in (3.2) are welcome.

Selected references

  • P. D. Seymour, The Matroids with the Max-Flow Min-Cut Property, J. Combin. Theory Ser. B 23 (1977) 189–222. https://doi.org/10.1016/0095-8956(77)90031-4
  • J. Edmonds and D. R. Fulkerson, Bottleneck extrema, J. Combin. Theory 8 (1970) 299–306. https://doi.org/10.1016/S0021-9800(70)80083-7
  • L. R. Ford and D. R. Fulkerson, Maximal flow through a network, Canad. J. Math. 8 (1956) 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • G. Cornuéjols, Combinatorial Optimization: Packing and Covering, CBMS-NSF Regional Conf. Ser. in Appl. Math. 74, SIAM, 2001. https://doi.org/10.1137/1.9780898717105
30 thms3 active usersReviewed
Quantum Error CorrectionQuantum Information·Captain: mikedeng1

Quantum Tanner Codes 1: On a Ramanujan Left-Right Cayley Complex with Robust Dual Tensor Codes, the Quantum Tanner Code Has Rate at Least (1−2ρ)² and Distance at Least δn/(4Δ^{3/2+ε})Research Paper

Why quantum Tanner codes matter

An error-correcting code protects information by adding redundancy. A quantum code must also account for two kinds of errors, commonly called XXX and ZZZ errors. Quantum low-density parity-check codes use checks that each touch only a bounded number of qubits, while each qubit participates in only a bounded number of checks. Such locality matters when the number of qubits grows. For years it was difficult to make both the number of encoded qubits and the minimum distance proportional to the block length while retaining bounded-weight checks. Leverrier and Zémor's quantum Tanner construction gives a route through local codes placed on an expanding square complex; its main distance theorem is the target of this mission Leverrier and Zémor, 2022.

The construction combines two earlier ingredients. Left-right Cayley complexes provide the squares and expansion properties used by the code Dinur, Evra, Livne, Lubotzky, and Mozes, 2021. Robustness of a dual tensor code, including after puncturing coordinates, controls local words Panteleev and Kalachev, 2021. The paper states the local reduction property as Theorem 1 and derives the rate and distance bounds as Theorem 17. This mission formalizes that implication under the paper's stated component-code and expansion conditions; the separate random-code result that supplies robust examples belongs to the second mission of this series.

The square complex and its codes

Let GGG be a finite group and let A,B⊆GA,B\subseteq GA,B⊆G be inverse-closed sets, each of size Δ\DeltaΔ. The total no-conjugacy condition says ag≠gbag\ne gbag=gb for every a∈Aa\in Aa∈A, b∈Bb\in Bb∈B, and g∈Gg\in Gg∈G. Vertices are two copies of the group, V=G×{0,1}V=G\times\{0,1\}V=G×{0,1}. A square associated with v=(g,i)v=(g,i)v=(g,i) and (a,b)∈A×B(a,b)\in A\times B(a,b)∈A×B has vertices vvv, (ag,1−i)(ag,1-i)(ag,1−i), (gb,1−i)(gb,1-i)(gb,1−i), and (agb,i)(agb,i)(agb,i). The set of squares is QQQ, and a binary word on QQQ has length n=∣Q∣n=|Q|n=∣Q∣. At a vertex vvv, the incident squares are indexed by A×BA\times BA×B; this indexing makes the local view xvx_vxv​ of a word x∈F2Qx\in\mathbb F_2^Qx∈F2Q​ an A×BA\times BA×B matrix Leverrier and Zémor, §§3.2 and 4.1.

A binary linear code CAC_ACA​ has coordinates indexed by AAA, and CBC_BCB​ has coordinates indexed by BBB. Their tensor code CA⊗CBC_A\otimes C_BCA​⊗CB​ consists of matrices whose columns lie in CAC_ACA​ and whose rows lie in CBC_BCB​. The dual tensor code CA⊗F2B+F2A⊗CBC_A\otimes\mathbb F_2^B+\mathbb F_2^A\otimes C_BCA​⊗F2B​+F2A​⊗CB​ is the sum of the column and row codes. In Definition 5, it is www-robust when every member of weight at most www is supported on a union of few rows and columns, with the bounds determined by the two component minimum distances. Definition 8 requires this property after deleting up to ppp coordinates from each side. These are conditions on all relevant codewords and all permitted puncturings, including the unpunctured pair Leverrier and Zémor, pp. 15–16.

The two Tanner codes are C0=T(G0□,(CA⊗CB)⊥)\mathcal C_0=T(\mathcal G^\square_0,(C_A\otimes C_B)^\perp)C0​=T(G0□​,(CA​⊗CB​)⊥) and C1=T(G1□,(CA⊥⊗CB⊥)⊥)\mathcal C_1=T(\mathcal G^\square_1,(C_A^\perp\otimes C_B^\perp)^\perp)C1​=T(G1□​,(CA⊥​⊗CB⊥​)⊥). Here Gi□\mathcal G^\square_iGi□​ is the square-adjacency graph on the vertices in copy iii. Together they define the quantum CSS code Q=(C0,C1)\mathcal Q=(\mathcal C_0,\mathcal C_1)Q=(C0​,C1​). Its dimension counts encoded qubits. Its distance is the smaller of the minimum weights in C0∖C1⊥\mathcal C_0\setminus\mathcal C_1^\perpC0​∖C1⊥​ and C1∖C0⊥\mathcal C_1\setminus\mathcal C_0^\perpC1​∖C0⊥​ Leverrier and Zémor, pp. 14, 18.

Formalization targets

The goal is Theorem 17. Fix 0<ρ<1/20<\rho<1/20<ρ<1/2, 0<ε<1/20<\varepsilon<1/20<ε<1/2, 1/2+ε<γ<11/2+\varepsilon<\gamma<11/2+ε<γ<1, and δ>0\delta>0δ>0. For sufficiently large Δ\DeltaΔ, take a left-right Cayley complex satisfying total no-conjugacy and the paper's Ramanujan conditions. Require 0<dim⁡CA≤ρΔ0<\dim C_A\le\rho\Delta0<dimCA​≤ρΔ, dim⁡CB=Δ−dim⁡CA\dim C_B=\Delta-\dim C_AdimCB​=Δ−dimCA​, minimum distance at least δΔ\delta\DeltaδΔ for CA,CBC_A,C_BCA​,CB​ and both duals, and Δ3/2−ε/2\Delta^{3/2-\varepsilon/2}Δ3/2−ε/2 robustness with Δγ\Delta^\gammaΔγ resistance to puncturing for both relevant dual tensor codes. The target is

k≥(1−2ρ)2n,d≥δn4Δ3/2+ε.k\ge (1-2\rho)^2 n,\qquad d\ge \frac{\delta n}{4\Delta^{3/2+\varepsilon}}.k≥(1−2ρ)2n,d≥4Δ3/2+εδn​.

The milestones mirror the paper's route to this result. Lemma 4 bounds the nontrivial eigenvalues of the union and square graphs. Propositions 6 and 7 turn dual-tensor robustness into a decomposition and a 3/23/23/2 closeness bound. Theorem 1 supplies a weight-reducing local codeword for every sufficiently light nonzero Tanner-code word. Display (5) supplies the dimension lower bound. Theorem 17 combines the two distance directions with that rate estimate Leverrier and Zémor, pp. 12, 15–16, 18–23.

What the result supplies

The lower bound on kkk ensures a positive fraction of qubits can carry logical information when ρ<1/2\rho<1/2ρ<1/2. The lower bound on ddd rules out small undetectable logical errors. For a fixed allowable degree Δ\DeltaΔ, both bounds scale linearly with nnn. The paper also shows that random component codes can satisfy the robustness assumptions; that probabilistic result is a separate formalization target Leverrier and Zémor, Theorems 9 and 17.

A checked development would make explicit how the square indexing, duals, punctured codes, and CSS quotient fit together. The source proves the mathematical result; the statements here remain open Lean theorems. Reusable contributions include the tensor-code interface, spectral bounds for the graph construction, and lemmas connecting local views to the global Tanner codes. A proof of the existence of infinite Ramanujan complex families is outside this mission's stated target.

Where the proof is difficult

A light global codeword need not have a local view that is itself close to a tensor-code word. Some square-graph vertices have unusually large local weight, and puncturing is needed to isolate coordinates shared with those exceptional vertices. The paper's small-set expansion estimates and robustness conditions control those interactions. Simply averaging local weights does not yield the weight-reducing local codeword asserted by Theorem 1 Leverrier and Zémor, §4.2.

Formalization scope

Words are functions into ZMod 2, codes are submodules, and a matrix is a function on A×BA\times BA×B with AAA indexing rows. The square coordinates are four-element vertex sets. The construction uses V=G×Fin⁡(2)V=G\times\operatorname{Fin}(2)V=G×Fin(2) and counts nnn as the number of squares rather than building n=Δ2∣G∣/2n=\Delta^2|G|/2n=Δ2∣G∣/2 into a definition. The minimum distance of the zero code is infinity, and robustness uses inequalities involving that extended value. Powers of Δ\DeltaΔ and the parameters w,pw,pw,p are real; the paper uses integer w,pw,pw,p in Definition 8, while the theorem's power thresholds need real comparisons. The formal robustness predicate includes Definition 5's range 0≤w≤Δ20\le w\le\Delta^20≤w≤Δ2; the theorem thresholds satisfy it for sufficiently large Δ\DeltaΔ. These conventions are recorded with each item.

The goal quantifies over every finite complex satisfying inverse closure, total no-conjugacy, and Ramanujan conditions, with a uniform threshold Δ0\Delta_0Δ0​ chosen before the group and codes. This pins the complexes in the cited infinite family without formalizing their existence. In particular, the goal contains no auxiliary set or local reduction assumption that could make the distance bound immediate. The published CodingTheory and TillichZemor.HGP.Setting definitions supply dual codes and CSS dimension; the paper-specific setting supplies the square complex, local views, Tanner codes, and robustness. Contributions that keep these conventions and close any milestone theorem are welcome.

Selected references

  • A. Leverrier and G. Zémor, Quantum Tanner codes, arXiv preprint, 2022. arXiv:2202.13641v3.
  • I. Dinur, S. Evra, R. Livne, A. Lubotzky, and S. Mozes, Locally Testable Codes with constant rate, distance, and locality, arXiv preprint, 2021. arXiv:2111.04808.
  • P. Panteleev and G. Kalachev, Asymptotically good quantum and locally testable classical LDPC codes, arXiv preprint, 2021. arXiv:2111.03654.
9 thms2 active usersReviewed
OptimizationQuantum Information·Captain: mikedeng1

A Quantum Approximate Optimization Algorithm II: On v-Regular Graphs the Variance of the Cut Value in the QAOA State Is at Most 2((v−1)^(2p+2)−1)/(v−2) Times the Number of EdgesResearch Paper

Motivation

The Quantum Approximate Optimization Algorithm (QAOA) of Farhi, Goldstone and Gutmann (arXiv:1411.4028, 2014) is a family of shallow quantum circuits for combinatorial optimization, indexed by a depth ppp and 2p2p2p real angles. It has become one of the standard candidate algorithms for near-term quantum hardware, and its MaxCut instance is the usual benchmark on which it is analysed.

The algorithm is run by preparing a quantum state, measuring it in the computational basis to obtain a bit string zzz, and evaluating the classical objective C(z)C(z)C(z). Its analysis in the paper bounds the mean FpF_pFp​ of C(z)C(z)C(z). A mean alone does not say what a single run returns. Section III of the paper supplies the missing piece for MaxCut on bounded-degree graphs: the variance of C(z)C(z)C(z) grows at most linearly in the number of edges, so the measured cut values are concentrated around FpF_pFp​ and a modest number of repetitions reproduces the mean. This mission formalizes that variance bound, equation (30) of the paper, together with the locality and counting facts it rests on.

Setting

Let GGG be a simple graph on vertices {0,…,n−1}\{0,\dots,n-1\}{0,…,n−1} with edge set {⟨jk⟩}\{\langle jk\rangle\}{⟨jk⟩} of size mmm, each edge listed once. Vertices are qubits; a computational basis state is a bit string z∈{0,1}nz\in\{0,1\}^nz∈{0,1}n.

  • The clause operator of an edge is C⟨jk⟩=12(−σjzσkz+1)C_{\langle jk\rangle}=\tfrac12(-\sigma^z_j\sigma^z_k+1)C⟨jk⟩​=21​(−σjz​σkz​+1). It is diagonal: C⟨jk⟩(z)=1C_{\langle jk\rangle}(z)=1C⟨jk⟩​(z)=1 if zj≠zkz_j\ne z_kzj​=zk​ and 000 otherwise.
  • The MaxCut objective is C=∑⟨jk⟩C⟨jk⟩C=\sum_{\langle jk\rangle}C_{\langle jk\rangle}C=∑⟨jk⟩​C⟨jk⟩​; C(z)C(z)C(z) is the number of edges cut by zzz.
  • The mixer is B=∑jσjxB=\sum_j\sigma^x_jB=∑j​σjx​, and U(C,γ)=e−iγCU(C,\gamma)=e^{-i\gamma C}U(C,γ)=e−iγC, U(B,β)=e−iβBU(B,\beta)=e^{-i\beta B}U(B,β)=e−iβB.
  • The initial state is ∣s⟩=2−n/2∑z∣z⟩|s\rangle=2^{-n/2}\sum_z|z\rangle∣s⟩=2−n/2∑z​∣z⟩.
  • For p≥1p\ge1p≥1 and angles γ=(γ1,…,γp)\boldsymbol\gamma=(\gamma_1,\dots,\gamma_p)γ=(γ1​,…,γp​), β=(β1,…,βp)\boldsymbol\beta=(\beta_1,\dots,\beta_p)β=(β1​,…,βp​), the QAOA state is
∣γ,β⟩=U(B,βp) U(C,γp)⋯U(B,β1) U(C,γ1) ∣s⟩.|\boldsymbol\gamma,\boldsymbol\beta\rangle=U(B,\beta_p)\,U(C,\gamma_p)\cdots U(B,\beta_1)\,U(C,\gamma_1)\,|s\rangle .∣γ,β⟩=U(B,βp​)U(C,γp​)⋯U(B,β1​)U(C,γ1​)∣s⟩.
  • Measuring it gives zzz with probability P(z)=∣⟨z∣γ,β⟩∣2P(z)=|\langle z|\boldsymbol\gamma,\boldsymbol\beta\rangle|^2P(z)=∣⟨z∣γ,β⟩∣2. The mean is Fp(γ,β)=⟨γ,β∣C∣γ,β⟩=∑zP(z)C(z)F_p(\boldsymbol\gamma,\boldsymbol\beta)=\langle\boldsymbol\gamma,\boldsymbol\beta|C|\boldsymbol\gamma,\boldsymbol\beta\rangle=\sum_zP(z)C(z)Fp​(γ,β)=⟨γ,β∣C∣γ,β⟩=∑z​P(z)C(z) and the variance is
⟨γ,β∣C2∣γ,β⟩−⟨γ,β∣C∣γ,β⟩2.\langle\boldsymbol\gamma,\boldsymbol\beta|C^2|\boldsymbol\gamma,\boldsymbol\beta\rangle-\langle\boldsymbol\gamma,\boldsymbol\beta|C|\boldsymbol\gamma,\boldsymbol\beta\rangle^2 .⟨γ,β∣C2∣γ,β⟩−⟨γ,β∣C∣γ,β⟩2.
  • The light cone (the paper's subgraph g(j,k)g(j,k)g(j,k)) of an edge is the set of qubits at graph distance at most ppp from jjj or from kkk.

In Lean these are edgeCut, cutValue, qaoaState, expectCut and edgeTerm in the shared FarhiQAOA.Cubic setting, and edgePairs, cutVariance, edgePairTerm, edgeCov and lightCone in FarhiQAOA.Concentration, all built on the published QAlgorithms definitions of qubit registers, measurement probabilities, the uniform superposition, the problem Hamiltonian and the QAOA unitary.

Formalization targets

Goal: the variance bound (30)

For a vvv-regular graph with v≥3v\ge3v≥3, every depth p≥1p\ge1p≥1 and all real angles,

⟨γ,β∣C2∣γ,β⟩−⟨γ,β∣C∣γ,β⟩2≤2[(v−1)2p+2−1(v−1)−1]⋅m.\langle\boldsymbol\gamma,\boldsymbol\beta|C^2|\boldsymbol\gamma,\boldsymbol\beta\rangle-\langle\boldsymbol\gamma,\boldsymbol\beta|C|\boldsymbol\gamma,\boldsymbol\beta\rangle^2\le2\left[\frac{(v-1)^{2p+2}-1}{(v-1)-1}\right]\cdot m .⟨γ,β∣C2∣γ,β⟩−⟨γ,β∣C∣γ,β⟩2≤2[(v−1)−1(v−1)2p+2−1​]⋅m.

The constant is the paper's, as printed.

Milestones

  1. (27)–(28). The variance is the double sum, over ordered pairs of edges, of the covariances ⟨C⟨jk⟩C⟨j′k′⟩⟩−⟨C⟨jk⟩⟩⟨C⟨j′k′⟩⟩\langle C_{\langle jk\rangle}C_{\langle j'k'\rangle}\rangle-\langle C_{\langle jk\rangle}\rangle\langle C_{\langle j'k'\rangle}\rangle⟨C⟨jk⟩​C⟨j′k′⟩​⟩−⟨C⟨jk⟩​⟩⟨C⟨j′k′⟩​⟩.
  2. Disjoint light cones. If g(j,k)g(j,k)g(j,k) and g(j′,k′)g(j',k')g(j′,k′) share no qubit, the covariance of the two edges is 000.
  3. (26). For maximum degree v≥3v\ge3v≥3, ∣g(j,k)∣≤2[((v−1)p+1−1)/((v−1)−1)]|g(j,k)|\le2\big[((v-1)^{p+1}-1)/((v-1)-1)\big]∣g(j,k)∣≤2[((v−1)p+1−1)/((v−1)−1)].
  4. (29). For maximum degree v≥3v\ge3v≥3, at most 2[((v−1)2p+2−1)/((v−1)−1)]2\big[((v-1)^{2p+2}-1)/((v-1)-1)\big]2[((v−1)2p+2−1)/((v−1)−1)] edges ⟨j′k′⟩\langle j'k'\rangle⟨j′k′⟩ have g(j′,k′)g(j',k')g(j′,k′) meeting g(j,k)g(j,k)g(j,k).
  5. Summand bound. Each covariance is at most 111 in absolute value.

Companions, printed in the paper for v=2v=2v=2: ∣g(j,k)∣≤2p+2|g(j,k)|\le2p+2∣g(j,k)∣≤2p+2 and at most 4p+44p+44p+4 meeting edges.

Significance

The result. For fixed vvv and ppp the bound makes the standard deviation of the cut value O(m)O(\sqrt m)O(m​), while the mean is of order mmm. The algorithm's output is therefore concentrated: a single measurement already returns a cut close to FpF_pFp​ in relative terms, and the paper uses this to argue that the sample mean of polynomially many runs estimates FpF_pFp​ within 111. The same locality structure, that each edge term depends only on a bounded neighbourhood, is what makes FpF_pFp​ classically computable for fixed ppp and underlies later analyses of QAOA on bounded-degree and random regular graphs.

Formalizing it. The bound is proved in the paper in a few lines; none of it is machine-checked. A formal proof has to make the paper's operator locality argument precise (the conjugated clause operator is supported on the light cone, and the initial state is a product state), connect it to the computational-basis probabilities, and carry out the graph counting with explicit constants. The resulting locality lemma for QAOA circuits on graphs is reusable for any statement about QAOA at fixed depth.

Difficulty

The algebra of (28) and the bound by 111 are routine. The substance is in two places.

The vanishing of a covariance requires showing that U†C⟨jk⟩UU^\dagger C_{\langle jk\rangle}UU†C⟨jk⟩​U, for the full depth-ppp QAOA unitary UUU, acts only on the qubits of g(j,k)g(j,k)g(j,k), and that for two operators acting on disjoint sets of qubits the expectation in the product state ∣s⟩|s\rangle∣s⟩ factorizes. Both are statements about operators on (C2)⊗n(\mathbb C^2)^{\otimes n}(C2)⊗n, whereas the objects in the statement are probabilities of computational basis states; the formal proof must pass between the two. Matrix exponentials of commuting sums have to be split into commuting factors, which in Lean means working with NormedSpace.exp on matrices indexed by bit strings.

The counting (26) and (29) is a ball-growth bound in a graph of bounded degree. The paper's count, via (26) with ppp replaced by 2p+12p+12p+1, bounds vertices, while (29) counts edges; the passage between the two, and the extended graph distance on disconnected graphs, must be made exact with the printed constant.

Formalization scope

  • Qubits are indexed by Fin n, bit strings are Fin n → Bool with true the bit 111, and states live in EuclideanSpace ℂ (Fin n → Bool).
  • All expectations are sums ∑zP(z) f(z)\sum_zP(z)\,f(z)∑z​P(z)f(z) over the computational basis. This equals the operator expectation in the paper because every operator involved is diagonal.
  • Each edge is counted once (pairs j<kj<kj<k), in mmm, in CCC, in the double sum of (28) and in the count (29).
  • Graph distance is the extended distance SimpleGraph.edist (+∞+\infty+∞ between components).
  • The goal is stated for vvv-regular graphs with v≥3v\ge3v≥3, the domain of the printed formula. The milestones (26) and (29) assume only maximum degree ≤v\le v≤v, as printed there.
  • The angles range over all of Rp\mathbb R^pRp. The paper restricts them to [0,2π][0,2\pi][0,2π] and [0,π][0,\pi][0,π] by periodicity; the bound does not depend on that restriction.
  • The QAOA state is the state itself, not an arbitrary vector: normalisation of the measurement probabilities follows from unitarity and is part of what must be proved.

A variance defined with two different states, or a probability vector that is not normalised, would make the goal vacuous or false; the definitions use the single state ∣γ,β⟩|\boldsymbol\gamma,\boldsymbol\beta\rangle∣γ,β⟩ throughout.

Out of scope: the O(m)O(\sqrt m)O(m​) standard-deviation reading and the sampling claim after (30), which carry no explicit constant; the measurement-and-repetition procedure, running time, circuit depth and the classical search for the angles.

Welcome contributions: a general light-cone lemma for conjugation by products of exponentials of local Hamiltonians on qubits; factorization of product-state expectations over disjoint qubit sets; ball-size bounds in graphs of bounded degree.

Selected references

  • E. Farhi, J. Goldstone, S. Gutmann, A Quantum Approximate Optimization Algorithm, arXiv preprint, 2014. https://arxiv.org/abs/1411.4028 (v1; Section III, pp. 6–7)
  • E. Farhi, J. Goldstone, S. Gutmann, M. Sipser, Quantum Computation by Adiabatic Evolution, arXiv preprint, 2000. https://arxiv.org/abs/quant-ph/0001106
14 thms2 active usersReviewed
OptimizationQuantum Information·Captain: mikedeng1

A Quantum Approximate Optimization Algorithm I: At p = 1 the QAOA Expected Cut Is at Least 0.6924 Times the Maximum Cut on Every Connected 3-Regular GraphResearch Paper

Motivation

The quantum approximate optimization algorithm (QAOA) of Farhi, Goldstone and Gutmann (arXiv:1411.4028, 2014) is a family of shallow quantum circuits, indexed by an integer depth p≥1p \ge 1p≥1, that prepare states whose computational-basis measurements are biased toward good solutions of a combinatorial optimization problem. It became the standard candidate for optimization on near-term quantum hardware, and its performance on MaxCut is the benchmark against which later analyses and experiments are measured.

The original paper proves one concrete performance guarantee: at depth p=1p = 1p=1, on 3-regular graphs, the algorithm achieves an approximation ratio of at least 0.69240.69240.6924. Short timeline:

  • 2014, Farhi–Goldstone–Gutmann: the p=1p = 1p=1 bound 0.69240.69240.6924 on 3-regular graphs, by reducing the expected cut to three subgraph functions and a numerical minimization (arXiv:1411.4028).
  • 2018, Wang, Hadfield, Jiang and Rieffel: closed-form expressions for the p=1p = 1p=1 MaxCut expectation on arbitrary graphs (arXiv:1706.02998, Phys. Rev. A 97, 022304).
  • Classically, Goemans and Williamson's semidefinite rounding achieves 0.8780.8780.878 on all graphs (J. ACM 42, 1995), so the p=1p = 1p=1 guarantee is, as the paper itself notes, below the best classical algorithms.

Setting

A graph GGG has vertex set {0,…,n−1}\{0, \dots, n-1\}{0,…,n−1}, and 3-regular means every vertex has exactly three neighbours. Each vertex carries a qubit; a computational basis state is a bit string z=z0⋯zn−1z = z_0 \cdots z_{n-1}z=z0​⋯zn−1​. For an edge ⟨jk⟩\langle jk\rangle⟨jk⟩ the clause C⟨jk⟩(z)C_{\langle jk\rangle}(z)C⟨jk⟩​(z) is 111 if zj≠zkz_j \neq z_kzj​=zk​ and 000 otherwise, and the cut objective is

C(z)=∑⟨jk⟩C⟨jk⟩(z),C(z) = \sum_{\langle jk\rangle} C_{\langle jk\rangle}(z),C(z)=⟨jk⟩∑​C⟨jk⟩​(z),

the number of edges cut by the bipartition zzz defines. As an operator, CCC is diagonal in the computational basis.

Let B=∑jσjxB = \sum_j \sigma^x_jB=∑j​σjx​ and let ∣s⟩=2−n/2∑z∣z⟩|s\rangle = 2^{-n/2}\sum_z |z\rangle∣s⟩=2−n/2∑z​∣z⟩. For angles γ,β∈R\gamma, \beta \in \mathbb Rγ,β∈R, the depth-one QAOA state and its expected cut are

∣γ,β⟩=e−iβB e−iγC ∣s⟩,F1(γ,β)=⟨γ,β∣ C ∣γ,β⟩=∑z∣⟨z∣γ,β⟩∣2 C(z).|\gamma,\beta\rangle = e^{-i\beta B}\, e^{-i\gamma C}\, |s\rangle, \qquad F_1(\gamma,\beta) = \langle\gamma,\beta|\,C\,|\gamma,\beta\rangle = \sum_z |\langle z|\gamma,\beta\rangle|^2\, C(z).∣γ,β⟩=e−iβBe−iγC∣s⟩,F1​(γ,β)=⟨γ,β∣C∣γ,β⟩=z∑​∣⟨z∣γ,β⟩∣2C(z).

The angles may be restricted to γ∈[0,2π]\gamma \in [0, 2\pi]γ∈[0,2π], β∈[0,π]\beta \in [0, \pi]β∈[0,π] by periodicity.

Two kinds of small subgraphs govern the analysis. A crossed square is a set of four vertices spanning exactly five edges (a 4-cycle with one diagonal). An isolated triangle is a triangle whose three outgoing edges end on three distinct vertices. SSS and TTT denote their numbers in GGG. Finally, for an edge ⟨jk⟩\langle jk\rangle⟨jk⟩ the subgraph type is the shape of the edge together with its adjacent edges; in a 3-regular graph there are three types g4,g5,g6g_4, g_5, g_6g4​,g5​,g6​, according as jjj and kkk have two, one or no common neighbours. Each type ggg defines a function fg(γ,β)f_g(\gamma, \beta)fg​(γ,β), the expectation of the clause of the central edge in the depth-one QAOA state of the subgraph ggg alone.

Formalization targets

Goal: the depth-one ratio on cubic graphs

For every connected 3-regular graph GGG there are angles γ∈[0,2π]\gamma \in [0, 2\pi]γ∈[0,2π], β∈[0,π]\beta \in [0, \pi]β∈[0,π] with

F1(γ,β)  ≥  0.6924⋅C(z)for every bit string z,F_1(\gamma, \beta) \;\ge\; 0.6924 \cdot C(z) \qquad \text{for every bit string } z,F1​(γ,β)≥0.6924⋅C(z)for every bit string z,

that is, M1=max⁡γ,βF1≥0.6924⋅max⁡zC(z)M_1 = \max_{\gamma,\beta} F_1 \ge 0.6924 \cdot \max_z C(z)M1​=maxγ,β​F1​≥0.6924⋅maxz​C(z). The constant is exactly the paper's 0.69240.69240.6924, and K4K_4K4​ is included.

Milestones

The paper's argument passes through the following claims, each a milestone:

  1. At p=1p = 1p=1 each edge term of a 3-regular graph equals fg4f_{g_4}fg4​​, fg5f_{g_5}fg5​​ or fg6f_{g_6}fg6​​ according to its subgraph type (p. 5).
  2. In a connected 3-regular graph other than K4K_4K4​, there are SSS edges of type g4g_4g4​ and 4S+3T4S + 3T4S+3T of type g5g_5g5​ (pp. 8–9).
  3. Equation (32): F1=Sfg4+(4S+3T)fg5+(3n/2−5S−3T)fg6F_1 = S f_{g_4} + (4S + 3T) f_{g_5} + (3n/2 - 5S - 3T) f_{g_6}F1​=Sfg4​​+(4S+3T)fg5​​+(3n/2−5S−3T)fg6​​ (p. 9).
  4. 3T+4S≤n3T + 4S \le n3T+4S≤n (p. 8).
  5. Every cut has at most 3n/2−S−T3n/2 - S - T3n/2−S−T edges (p. 9).
  6. For s,t≥0s, t \ge 0s,t≥0 with 4s+3t≤14s + 3t \le 14s+3t≤1, the scaled ratio (35) M1(1,s,t)/(3/2−s−t)M_1(1, s, t)/(3/2 - s - t)M1​(1,s,t)/(3/2−s−t) is at least 0.69240.69240.6924 (p. 9).
  7. On K4K_4K4​ the depth-one ratio exceeds 0.69240.69240.6924 (p. 8).

Significance

The result is the first rigorous worst-case guarantee for the QAOA, and the method behind it (the expected objective is a sum of local terms, each depending only on a bounded neighbourhood of an edge) is the template for every later analysis of fixed-depth QAOA on bounded-degree graphs. It shows that a depth-one circuit, of depth independent of nnn, achieves a constant approximation ratio uniformly over an infinite class of graphs.

The paper's proof is complete apart from one numerical step: the minimization of (35) is reported as a numerical evaluation. No machine-checked version of the result exists. A formal proof requires an exact evaluation of the three subgraph functions at suitable angles, an exact combinatorial count of edge types, and a certified inequality in place of the numerical minimization. The margin is small: the exact minimum is 12+133≈0.692450\tfrac12 + \tfrac{1}{3\sqrt 3} \approx 0.69245021​+33​1​≈0.692450, about 5⋅10−55 \cdot 10^{-5}5⋅10−5 above the paper's constant.

Difficulty

The light-cone reduction (milestone 1) is the conceptual step: one must show that conjugating a two-qubit clause by e−iγCe^{-i\gamma C}e−iγC and e−iβBe^{-i\beta B}e−iβB leaves an operator that involves only the qubits adjacent to the edge, and that the remaining qubits contribute a factor of one. In matrix form on 2n2^n2n-dimensional space this is a statement about tensor-product structure that the definitions do not make visible. The combinatorial count (milestones 2–4) needs a careful case analysis of how triangles sit in a cubic graph, including the exceptional graph K4K_4K4​, where the count fails. The final inequality (milestone 6) has a margin of order 10−410^{-4}10−4, so crude bounds on the trigonometric functions fgf_{g}fg​ do not suffice.

Formalization scope

  • Representation. Bit strings are Qubits n = Fin n → Bool from the published QAlgorithms layer, with true the bit 111. States are vectors in EuclideanSpace ℂ (Qubits n); ∣s⟩|s\rangle∣s⟩, CCC as a diagonal operator, BBB and the QAOA unitary are the published definitions uniformState, problemHamiltonian, qaoaMixer, qaoaUnitary. The latter takes mixer angles before phase angles.
  • Conventions. Graphs are SimpleGraph (Fin n) with decidable adjacency; "3-regular" is IsRegularOfDegree 3. The cut counts each edge once. F1F_1F1​ is defined as ∑zPr⁡[z] C(z)\sum_z \Pr[z]\,C(z)∑z​Pr[z]C(z), which equals ⟨γ,β∣C∣γ,β⟩\langle\gamma,\beta|C|\gamma,\beta\rangle⟨γ,β∣C∣γ,β⟩ because CCC is diagonal. The constant 0.69240.69240.6924 is the exact rational 1731/25001731/25001731/2500.
  • Added hypotheses. Milestones 2 and 3 assume GGG connected and n≠4n \ne 4n=4, which excludes K4K_4K4​, the case the paper sets aside on p. 8. No other hypothesis is added.
  • Out of scope. The goal concerns the expected cut at good angles, which is what the paper's analysis bounds. Measuring the state, repeating, and the claim that a cut of that size is then found with high probability are not formalized; nor are running time, circuit depth, or the search for optimal angles.
  • Ruled-out trivializations. A cut objective counting each edge twice, or an unnormalized state, would leave the ratio statement formally intact but about the wrong objects; the setting counts each edge once and uses the normalized uniform superposition. The goal quantifies over all bit strings, so it cannot be satisfied by bounding a particular cut.
  • Contributions welcome. Reusable lemmas on the light-cone locality of conjugated Pauli operators, closed forms of the depth-one MaxCut expectation, and counting lemmas for triangles in cubic graphs all serve beyond this mission.

Selected references

  • E. Farhi, J. Goldstone, S. Gutmann, A Quantum Approximate Optimization Algorithm, arXiv preprint, 2014. https://arxiv.org/abs/1411.4028
  • Z. Wang, S. Hadfield, Z. Jiang, E. G. Rieffel, Quantum approximate optimization algorithm for MaxCut: A fermionic view, Phys. Rev. A 97, 022304, 2018. https://arxiv.org/abs/1706.02998
  • M. X. Goemans, D. P. Williamson, Improved approximation algorithms for maximum cut and satisfiability problems using semidefinite programming, J. ACM 42(6), 1995. https://doi.org/10.1145/227683.227684
14 thms2 active usersReviewed
Control TheoryDynamical SystemsLinear algebra·Captain: mikedeng1

Synchronization via Pinning Control on General Complex Networks II: Componentwise LMIs on the Strongly Connected Components Globally Synchronize a Pinned Network with a Directed Spanning TreeResearch Paper

Motivation

Networks of identical dynamical systems coupled along a directed graph model power grids, sensor and robot swarms, neural populations and coupled oscillators. A basic control question is how to steer every node of such a network onto a prescribed reference trajectory by acting directly on only a few nodes. This is pinning control: feedback is injected at a small set of pinned vertices, and the coupling has to carry the reference signal to the rest of the network.

For undirected or strongly connected networks, criteria for pinning synchronization go back to Wang and Chen (2002) and Chen, Liu and Lu (2007), and the first mission of this series treats the strongly connected case of Yu, Chen, Lü and Kurths (2013). Many networks of interest are not strongly connected: information flows down a hierarchy of groups, and the minimal requirement for the reference to reach every node is that the network, with the reference added as a virtual leader, contains a directed spanning tree rooted at the leader. Section 4 of that paper gives a criterion for this general case that splits into one condition per strongly connected component. This mission formalizes that criterion, Theorem 4.5, with the lemmas its proof uses and three corollaries.

Setting

There are NNN vertices with states xi(t)∈Rnx_i(t)\in\mathbb R^nxi​(t)∈Rn. The coupling configuration matrix G∈RN×NG\in\mathbb R^{N\times N}G∈RN×N has Gij≥0G_{ij}\ge0Gij​≥0 for i≠ji\ne ji=j, where Gij>0G_{ij}>0Gij​>0 records a connection from vertex jjj to vertex iii, and zero row sums, Gii=−∑j≠iGijG_{ii}=-\sum_{j\ne i}G_{ij}Gii​=−∑j=i​Gij​. Given a function f:Rn×R+→Rnf:\mathbb R^n\times\mathbb R_+\to\mathbb R^nf:Rn×R+​→Rn, a coupling strength ccc, an inner coupling matrix Γ∈Rn×n\Gamma\in\mathbb R^{n\times n}Γ∈Rn×n and pinning gains di≥0d_i\ge0di​≥0 (vertex iii is pinned when di>0d_i>0di​>0), the controlled network is

x˙i=f(xi,t)+c∑j=1NGijΓxj−c di Γ(xi−s),i=1,…,N,\dot x_i=f(x_i,t)+c\sum_{j=1}^N G_{ij}\Gamma x_j-c\,d_i\,\Gamma(x_i-s),\qquad i=1,\dots,N,x˙i​=f(xi​,t)+cj=1∑N​Gij​Γxj​−cdi​Γ(xi​−s),i=1,…,N,

where the reference sss solves s˙=f(s,t)\dot s=f(s,t)s˙=f(s,t). The network is globally synchronized when xi(t)−s(t)→0x_i(t)-s(t)\to0xi​(t)−s(t)→0 for every vertex and every pair of solutions.

The nonlinearity is controlled by Assumption 1: there is a matrix KKK with (x−y)T(f(x,t)−f(y,t))≤(x−y)TKΓ(x−y)(x-y)^{\mathsf T}(f(x,t)-f(y,t))\le(x-y)^{\mathsf T}K\Gamma(x-y)(x−y)T(f(x,t)−f(y,t))≤(x−y)TKΓ(x−y) for all x,yx,yx,y and t≥0t\ge0t≥0. Throughout §4 the paper also assumes that KKK and Γ\GammaΓ commute, that Γ\GammaΓ is positive definite, and puts θ=λmax⁡((K+KT)/2)≥0\theta=\lambda_{\max}((K+K^{\mathsf T})/2)\ge0θ=λmax​((K+KT)/2)≥0.

The vertices split into strongly connected components V1,…,VMV_1,\dots,V_MV1​,…,VM​. In the Frobenius normal form the components are ordered so that no edge goes from a later component to an earlier one, which makes GGG block lower triangular. For each component jjj the mission uses three matrices:

  • G~j\widetilde G_jGj​, the diagonal block of G−diag(d)G-\mathrm{diag}(d)G−diag(d) on VjV_jVj​;
  • AjA_jAj​, the within-component part of GGG with zero row sums;
  • Π~j=diag(Πj)\widetilde\Pi_j=\mathrm{diag}(\Pi_j)Πj​=diag(Πj​), where Πj\Pi_jΠj​ is a positive vector with ΠjTAj=0\Pi_j^{\mathsf T}A_j=0ΠjT​Aj​=0.

In Lean, component membership is a map blk : Fin N → Fin M, and Fin M carries the Frobenius order.

Formalization targets

Goal: Theorem 4.5

If, for every component jjj,

2θΠ~j+c(Π~jG~j+G~jTΠ~j)<0,(4.8)2\theta\widetilde\Pi_j+c\bigl(\widetilde\Pi_j\widetilde G_j+\widetilde G_j^{\mathsf T}\widetilde\Pi_j\bigr)<0, \tag{4.8}2θΠj​+c(Πj​Gj​+GjT​Πj​)<0,(4.8)

then the controlled network with a directed spanning tree is globally synchronized.

Milestones

  1. The θ bound (p. 1402). Commutation and Γ>0\Gamma>0Γ>0 turn Assumption 1 into (x−y)T(f(x,t)−f(y,t))≤θ(x−y)TΓ(x−y)(x-y)^{\mathsf T}(f(x,t)-f(y,t))\le\theta(x-y)^{\mathsf T}\Gamma(x-y)(x−y)T(f(x,t)−f(y,t))≤θ(x−y)TΓ(x−y).
  2. Lemma 2.13. The strict Schur complement test for a symmetric block matrix.
  3. Lemma 4.4. If Π~jG~j+G~jTΠ~j<0\widetilde\Pi_j\widetilde G_j+\widetilde G_j^{\mathsf T}\widetilde\Pi_j<0Πj​Gj​+GjT​Πj​<0 on every diagonal block of a block lower-triangular G‾\overline GG, then positive block weights Δj\Delta_jΔj​ exist with ΔΠ~G‾+G‾TΠ~Δ<0\Delta\widetilde\Pi\overline G+\overline G^{\mathsf T}\widetilde\Pi\Delta<0ΔΠG+GTΠΔ<0.
  4. (4.10). Along solutions, V=∑jΔje~jTΠ~je~jV=\sum_j\Delta_j\widetilde e_j^{\mathsf T}\widetilde\Pi_j\widetilde e_jV=∑j​Δj​ejT​Πj​ej​ satisfies V˙≤eT{[ΔΠ~(cG‾+θIN)+(cG‾+θIN)TΠ~Δ]⊗Γ}e\dot V\le e^{\mathsf T}\{[\Delta\widetilde\Pi(c\overline G+\theta I_N)+(c\overline G+\theta I_N)^{\mathsf T}\widetilde\Pi\Delta]\otimes\Gamma\}eV˙≤eT{[ΔΠ(cG+θIN​)+(cG+θIN​)TΠΔ]⊗Γ}e.

Companions

  • Corollary 4.7: the same conclusion under θΠ~j+c(A‾j−Π~jCj−Π~jD‾j)<0\theta\widetilde\Pi_j+c(\overline A_j-\widetilde\Pi_jC_j-\widetilde\Pi_j\overline D_j)<0θΠj​+c(Aj​−Πj​Cj​−Πj​Dj​)<0, which separates the outer in-degrees CjC_jCj​ from the pinning gains D‾j\overline D_jDj​.
  • Corollary 4.8 (i): an explicit coupling-strength threshold under Assumption 2.
  • Corollary 4.9: partial synchronization of the first components.

Significance

Theorem 4.5 replaces a single Nn×NnNn\times NnNn×Nn matrix inequality for the whole network with MMM inequalities, one of the size of each component. In practice the test runs component by component down the Frobenius order: decide which components need pinning, choose gains for each one separately, and get synchronization of the whole network. Corollary 4.7 makes this explicit in terms of the outer in-degrees: a component with no incoming edges from earlier components has to be pinned (Remark 4), while a component that receives enough weight from upstream needs no control. These conditions underlie the pinning scheme algorithm of §5 of the paper.

On the formal side, as far as a search of the platform and of Mathlib shows, none of these results is machine-checked. Mathlib has the positive semidefinite Schur complement and irreducible nonnegative matrices, but not the strict Schur complement, block-diagonal Lyapunov scaling, or Lyapunov arguments for networked ODEs with Kronecker-structured couplings. The mission produces, beyond the paper's theorem:

  • a reusable strict Schur complement lemma;
  • a block-triangular diagonal-stability lemma;
  • a template for Lyapunov estimates on coupled networks in Lean.

Difficulty

The obvious argument is to take one Lyapunov function weighted by a positive left null vector of the whole coupling matrix, as in the strongly connected case. Here it fails: the coupling matrix of a network that is not strongly connected has no positive left null vector, and condition (4.8) is local to each component. The central step is the passage from blockwise to global negative definiteness (Lemma 4.4). Condition (4.8) says nothing about the off-diagonal blocks that couple a component to earlier ones, and a block lower-triangular matrix whose symmetrized diagonal blocks are negative definite need not have a negative definite symmetric part; the component weights Δj\Delta_jΔj​ have to be found so that those couplings are dominated, which is a quantitative statement about Schur complements. Turning the resulting matrix inequality into convergence of the trajectories then needs a differential inequality on [0,∞)[0,\infty)[0,∞), a Kronecker-product definiteness argument with Γ\GammaΓ, and the commutation of KKK and Γ\GammaΓ to reach the θ bound.

Formalization scope

States are Fin n → ℝ and matrices are indexed by Fin N, Fin n, or the subtype of one component's vertices. Solutions are functions on ℝ that satisfy their ODE on [0,∞)[0,\infty)[0,∞) through HasDerivWithinAt … (Set.Ici 0). "For any initial conditions" means for every pair of solutions (s,x)(s,x)(s,x); existence of solutions is not claimed. Convergence is componentwise, which agrees with every norm on Rn\mathbb R^nRn. "Positive definite" is Mathlib's Matrix.PosDef, which includes symmetry, and X<0X<0X<0 is (-X).PosDef. λmax⁡\lambda_{\max}λmax​ is the least Rayleigh bound, IsLamMax, never a supremum of a possibly empty set.

The vertices are not permuted into Frobenius form. The hypotheses say instead that:

  • every component is nonempty;
  • every component is strongly connected inside itself;
  • every edge between components points forward in the order;
  • every vertex is reachable from a pinned vertex.

The weights Πj\Pi_jΠj​ are tied to the paper's object by ΠjTAj=0\Pi_j^{\mathsf T}A_j=0ΠjT​Aj​=0. The standing assumptions of p. 1402 (Assumption 1, KΓ=ΓKK\Gamma=\Gamma KKΓ=ΓK, θ=λmax⁡((K+KT)/2)≥0\theta=\lambda_{\max}((K+K^{\mathsf T})/2)\ge0θ=λmax​((K+KT)/2)≥0, Γ>0\Gamma>0Γ>0) are hypotheses of every goal-level statement. The goal does not assume the weights Δj\Delta_jΔj​, the bound (4.10), or definiteness of a whole-network matrix; it assumes only (4.8) on each component, so it cannot be trivialized by assuming its own conclusion. For n=0n=0n=0 the λmax⁡\lambda_{\max}λmax​ hypothesis cannot be met, so that degenerate dimension is excluded.

The pages contain three misprints. (4.10) has "−c-c−c" in its third line and "++++++" in its last; the last line, with a single +++, is formalized. Corollary 4.9 ranges "i=1,…,Ni=1,\dots,Ni=1,…,N" where the component range is meant.

A complete development needs a Gronwall-type decay lemma on [0,∞)[0,\infty)[0,∞), the strict Schur complement, definiteness of Kronecker products, and the square root of a positive definite matrix commuting with matrices that commute with it. Each of these is reusable well beyond this mission, and contributions to any of them are welcome, as are proofs of the milestones in any order.

Selected references

  • W. Yu, G. Chen, J. Lü, J. Kurths, Synchronization via pinning control on general complex networks, SIAM J. Control Optim. 51(2) (2013), 1395–1416. https://doi.org/10.1137/100781699
  • X. F. Wang, G. Chen, Pinning control of scale-free dynamical networks, Physica A 310 (2002), 521–531. https://doi.org/10.1016/S0378-4371(02)00772-0
  • T. Chen, X. Liu, W. Lu, Pinning complex networks by a single controller, IEEE Trans. Circuits Syst. I 54(6) (2007), 1317–1326. https://doi.org/10.1109/TCSI.2007.895383
  • S. Boyd, L. El Ghaoui, E. Feron, V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, SIAM, 1994. https://doi.org/10.1137/1.9781611970777
  • R. A. Brualdi, H. J. Ryser, Combinatorial Matrix Theory, Cambridge University Press, 1991. https://doi.org/10.1017/CBO9781107325708
7 thms2 active usersReviewed
CombinatoricsProbability·Captain: mikedeng1

Algorithms for #BIS-Hard Problems on Expander Graphs 4: On Almost Every Δ-Regular Bipartite Graph with Δ ≥ Cq² log² q, Proper q-Colorings Reduce to Kotecký–Preiss Polymer ModelsResearch Paper

Motivation

Counting the proper qqq-colorings of a graph, and sampling one uniformly at random, are basic problems of approximate counting with a long history in statistical physics, where proper colorings are the zero-temperature antiferromagnetic Potts model. On bipartite graphs these problems are tied to #BIS, the problem of approximately counting independent sets in a bipartite graph, whose complexity is unresolved. Galanis, Štefankovič, Vigoda and Yang showed that approximately counting proper qqq-colorings is #BIS-hard on bipartite graphs of maximum degree Δ\DeltaΔ whenever q≤Δ/(2log⁡Δ)q \le \Delta/(2\log\Delta)q≤Δ/(2logΔ) (SIAM J. Comput. 45(6), 2016). A positive result in that regime therefore has to use structure beyond bipartiteness and bounded degree.

Jenssen, Keevash and Perkins (SIAM J. Comput. 49(4) 2020) use expansion. Their Theorem 5 says that for Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q there is an FPTAS and a polynomial-time sampler for proper qqq-colorings on almost every Δ\DeltaΔ-regular bipartite graph. This answered Conjecture 1 of the authors' SODA 2019 extended abstract. Liao, Lin, Lu and Mao (APPROX/RANDOM 2019) independently proved a similar result under a slightly stronger condition on Δ\DeltaΔ. The organizing idea, that most colorings of a strongly expanding bipartite graph follow one of a few patterns, comes from the work of Peled and Spinka on rigidity of colorings of Zd\mathbb Z^dZd (arXiv:1808.03597).

This mission is the fourth of a series formalizing the four main results of that paper; it covers Theorem 5.

Setting

Fix integers q≥3q \ge 3q≥3, Δ≥3\Delta \ge 3Δ≥3 and m≥1m \ge 1m≥1. The vertex set is V=O⊔EV = \mathcal O \sqcup \mathcal EV=O⊔E, two sides of mmm vertices each (so n=2mn = 2mn=2m). Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) is the finite set of labelled simple graphs on VVV whose edges all join O\mathcal OO to E\mathcal EE and in which every vertex has degree Δ\DeltaΔ. A property holds for almost every Δ\DeltaΔ-regular bipartite graph if the fraction of graphs in Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) having it tends to 111 as m→∞m \to \inftym→∞.

For S⊆VS \subseteq VS⊆V the vertex boundary ∂S\partial S∂S is the set of vertices outside SSS with a neighbour in SSS, and S+=S∪∂SS^+ = S \cup \partial SS+=S∪∂S. The graph is a bipartite (σ,ρ)(\sigma,\rho)(σ,ρ)-expander if ∣∂S∣≥ρ∣S∣|\partial S| \ge \rho|S|∣∂S∣≥ρ∣S∣ for every SSS inside one side with ∣S∣≤σm|S| \le \sigma m∣S∣≤σm.

A proper qqq-coloring is a map f:V→[q]f : V \to [q]f:V→[q] with f(u)≠f(v)f(u) \ne f(v)f(u)=f(v) on every edge; Xq(G)\mathcal X_q(G)Xq​(G) is their set and ZG(q)=∣Xq(G)∣Z_G(q) = |\mathcal X_q(G)|ZG​(q)=∣Xq​(G)∣. A pattern is a pair (A,B)(A,B)(A,B) of disjoint sets with A∪B=[q]A \cup B = [q]A∪B=[q]; there are 2q2^q2q of them. A coloring agrees with (A,B)(A,B)(A,B) at vvv if v∈Ov \in \mathcal Ov∈O and f(v)∈Af(v) \in Af(v)∈A, or v∈Ev \in \mathcal Ev∈E and f(v)∈Bf(v) \in Bf(v)∈B. For S⊆VS \subseteq VS⊆V, χA,B(S)\chi_{A,B}(S)χA,B​(S) is the set of proper colorings that disagree with (A,B)(A,B)(A,B) at every vertex of SSS and agree everywhere else.

G3G^3G3 joins two vertices at GGG-distance at most 333. A set is little if ∣S∣≤4qlog⁡ΔΔm|S| \le 4q\frac{\log\Delta}{\Delta}m∣S∣≤4qΔlogΔ​m. A polymer is a little set that is connected in G3G^3G3; C(G)\mathcal C(G)C(G) is the set of polymers, and two polymers are compatible if their GGG-distance exceeds 333. For a pattern (A,B)(A,B)(A,B), the weight and the polymer partition function are

wA,B(γ)=∣χA,B(γ)∣∣A∣m∣B∣m,ΞA,B(G)=∑Γ∏γ∈ΓwA,B(γ),w_{A,B}(\gamma) = \frac{|\chi_{A,B}(\gamma)|}{|A|^m|B|^m}, \qquad \Xi_{A,B}(G) = \sum_{\Gamma} \prod_{\gamma\in\Gamma} w_{A,B}(\gamma),wA,B​(γ)=∣A∣m∣B∣m∣χA,B​(γ)∣​,ΞA,B​(G)=Γ∑​γ∈Γ∏​wA,B​(γ),

the sum over all families Γ⊆C(G)\Gamma \subseteq \mathcal C(G)Γ⊆C(G) of pairwise compatible polymers, the empty family included. Finally Z^\hat ZZ^ is an ε\varepsilonε-relative approximation to ZZZ if e−εZ^≤Z≤eεZ^e^{-\varepsilon}\hat Z \le Z \le e^{\varepsilon}\hat Ze−εZ^≤Z≤eεZ^.

Formalization targets

Goal: the deterministic core of Theorem 5

There is an absolute constant C>0C > 0C>0 such that for all q≥3q \ge 3q≥3 and Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q, almost every Δ\DeltaΔ-regular bipartite graph GGG satisfies

e−εZ~G(q)≤ZG(q)≤eεZ~G(q),Z~G(q)=∑(A,B)∣A∣m∣B∣m ΞA,B(G),ε=e−m/(4q),e^{-\varepsilon}\tilde Z_G(q) \le Z_G(q) \le e^{\varepsilon}\tilde Z_G(q), \qquad \tilde Z_G(q) = \sum_{(A,B)} |A|^m|B|^m\,\Xi_{A,B}(G), \quad \varepsilon = e^{-m/(4q)},e−εZ~G​(q)≤ZG​(q)≤eεZ~G​(q),Z~G​(q)=(A,B)∑​∣A∣m∣B∣mΞA,B​(G),ε=e−m/(4q),

and, for every pattern (A,B)(A,B)(A,B) and every polymer γ∈C(G)\gamma \in \mathcal C(G)γ∈C(G), the Kotecký–Preiss condition

∑γ′∈C(G): d(γ′,γ)≤3wA,B(γ′) e∣γ′∣+g(γ′)≤∣γ∣,g(γ)=Δ10q2log⁡Δ∣γ∣.\sum_{\gamma' \in \mathcal C(G):\ d(\gamma',\gamma)\le 3} w_{A,B}(\gamma')\,e^{|\gamma'| + g(\gamma')} \le |\gamma|, \qquad g(\gamma) = \frac{\Delta}{10q^2\log\Delta}|\gamma|.γ′∈C(G): d(γ′,γ)≤3∑​wA,B​(γ′)e∣γ′∣+g(γ′)≤∣γ∣,g(γ)=10q2logΔΔ​∣γ∣.

The constant CCC is left existential, as in the paper.

Milestones

In the order of the proof:

  1. Lemma 23, the expansion of almost every regular bipartite graph;
  2. Lemma 27, the bound on ∣χA,B(S)∣|\chi_{A,B}(S)|∣χA,B​(S)∣;
  3. Lemma 28, the covering of Xq(G)\mathcal X_q(G)Xq​(G) by χA,B(S)\chi_{A,B}(S)χA,B​(S) with SSS little;
  4. the product formula after (13);
  5. Claim 30, the size bound for sparse sets;
  6. the binomial estimate, an existing platform theorem;
  7. Lemma 29, the reduction;
  8. the weight bound of §5.1;
  9. the per-vertex estimate (15).

Significance

The goal is the whole mathematical content of Theorem 5. The paper's algorithm (its Theorem 8, adapted from Helmuth, Perkins and Regts, Algorithmic Pirogov–Sinai theory) takes exactly these two facts as input. Given the Kotecký–Preiss condition, the cluster expansion of each log⁡ΞA,B(G)\log\Xi_{A,B}(G)logΞA,B​(G) converges and truncating it gives the estimator; the reduction transfers the estimate to ZG(q)Z_G(q)ZG​(q). With the goal proved, the remaining distance to the FPTAS is algorithmic bookkeeping, not combinatorics.

The result has been proved since 2019 and is not formalized anywhere. A formal proof would add a machine-checked polymer-model reduction for a spin system with exponentially many ground states, a reusable model of uniformly random regular bipartite graphs, and a checked instance of the "expansion implies few ground states" argument used throughout this literature.

Difficulty

The obvious first idea is a one-ground-state expansion, as for the hard-core model, where configurations are close to "all even" or "all odd". For colorings there are 2q2^q2q patterns. A coloring can lie close to several of them, so the reduction must show that the sets χA,B(S)\chi_{A,B}(S)χA,B​(S) with SSS little overlap only negligibly. That needs an inclusion–exclusion step whose error is exponentially small in m/qm/qm/q. A second difficulty is the size of polymers: little sets are qqq times larger than the sets of the hard-core analysis. Claim 30 is what keeps the union of compatible polymers from becoming a large set that the expansion no longer controls. Finally, the constants must be uniform in qqq. The single condition Δ≥Cq2log⁡2q\Delta \ge Cq^2\log^2 qΔ≥Cq2log2q has to cover the expansion lemma, the reduction and the Kotecký–Preiss estimate at once.

Formalization scope

The FPTAS and the sampling algorithm are not formalized. Mathlib has no machine model in which "polynomial time" could be stated, and the running times carry unquantified constants. The goal is the mathematics the algorithm rests on: the reduction and the Kotecký–Preiss condition, with the paper's constants e−m/(4q)e^{-m/(4q)}e−m/(4q), 4qlog⁡Δ/Δ4q\log\Delta/\Delta4qlogΔ/Δ and Δ/(10q2log⁡Δ)\Delta/(10q^2\log\Delta)Δ/(10q2logΔ). Neither half may be dropped. Replacing e−m/(4q)e^{-m/(4q)}e−m/(4q) by an existential error, or summing the Kotecký–Preiss condition over compatible rather than incompatible polymers, would make the goal trivial or wrong.

Conventions committed to in Lean:

  • Vertices are Fin m ⊕ Fin m, and Gbip(2m,Δ)\mathcal G^{\mathrm{bip}}(2m,\Delta)Gbip(2m,Δ) is the set of labelled Δ\DeltaΔ-regular graphs with edges only between the two fixed sides. The paper does not fix a bipartition; this fixed-sides model is a disclosed choice, matching §5's ∣O∣=∣E∣=m|\mathcal O| = |\mathcal E| = m∣O∣=∣E∣=m.
  • Distances are Mathlib's extended distance edist, so vertices in different components are at distance ∞\infty∞ and never incompatible.
  • A pattern is encoded by A⊆[q]A \subseteq [q]A⊆[q] with B=[q]∖AB = [q]\setminus AB=[q]∖A. Patterns with an empty side are kept: their weights are 000 under Lean's x/0=0x/0 = 0x/0=0, so they contribute 000 to Z~G(q)\tilde Z_G(q)Z~G​(q), as on the page for large mmm.
  • log is the natural logarithm; the binomial estimate is used in nats.

Added hypotheses:

  • A,B≠∅A, B \ne \emptysetA,B=∅ in Lemma 27 and in the product formula, which divide by ∣A∣|A|∣A∣ and ∣B∣|B|∣B∣;
  • Δ≥3\Delta \ge 3Δ≥3 in Lemma 28 and Claim 30, the paper's standing convention (§1.3);
  • m≥1m \ge 1m≥1 in Lemma 29, which fails on the empty vertex set.

The goal needs none of these, since its limit is in mmm and its constant absorbs Δ\DeltaΔ.

A complete development needs a counting theory of random regular bipartite graphs (for Lemma 23, through Bassalygo's expansion theorem), counts of connected sets in bounded-degree graphs, and inclusion–exclusion over patterns. The random-graph model and the polymer-model definitions are reusable across this series and beyond. Proofs of any milestone, and alternative routes to Lemma 23, are welcome.

Selected references

  • M. Jenssen, P. Keevash and W. Perkins, Algorithms for #BIS-hard problems on expander graphs, SIAM J. Comput. 49(4) (2020) 681–710. https://doi.org/10.1137/19M1286669
  • A. Galanis, D. Štefankovič, E. Vigoda and L. Yang, Ferromagnetic Potts model: refined #BIS-hardness and related results, SIAM J. Comput. 45(6) (2016) 2004–2065.
  • R. Peled and Y. Spinka, Rigidity of proper colorings of Zd\mathbb Z^dZd, arXiv:1808.03597 (2018). https://arxiv.org/abs/1808.03597
  • C. Liao, J. Lin, P. Lu and Z. Mao, Counting independent sets and colorings on random regular bipartite graphs, APPROX/RANDOM 2019, Schloss Dagstuhl–Leibniz-Zentrum für Informatik (2019).
  • T. Helmuth, W. Perkins and G. Regts, Algorithmic Pirogov–Sinai theory, Probab. Theory Related Fields (2020); extended abstract at STOC 2019.
  • R. Kotecký and D. Preiss, Cluster expansion for abstract polymer models, Comm. Math. Phys. 103(3) (1986) 491–498.
  • L. A. Bassalygo, Asymptotically optimal switching circuits, Problems Inform. Transmission 17(3) (1981) 206–211.
13 thms2 active usersReviewed
Mathematical PhysicsProbability·Captain: mikedeng1

Algorithms for #BIS-Hard Problems on Expander Graphs 3: On Almost Every Δ-Regular Bipartite Graph, the Hard-Core Model at λ > 50 log²Δ/Δ Reduces to Two Kotecký–Preiss Polymer ModelsResearch Paper

Motivation

The hard-core model assigns weight to independent vertex sets of a graph. Its partition function counts those sets when the activity is one and weights larger sets more heavily when the activity increases. Counting independent sets on bipartite graphs is a standard benchmark for approximate counting: the general problem is difficult, while random regular graphs have additional expansion that can make approximation possible. Jenssen, Keevash and Perkins prove an approximation and sampling theorem for almost every sufficiently high degree regular bipartite graph when the activity is above an explicit threshold Jenssen–Keevash–Perkins, Theorem 2.

The question is specific to a random graph model. It does not assert the same conclusion for every regular bipartite graph. The authors use a criterion attributed to Bassalygo to show that almost every graph in the model has the expansion needed by their argument Jenssen–Keevash–Perkins, Theorem 22 and Lemma 23. Liao, Lin, Lu and Mao independently obtained related counting results under somewhat stronger parameter conditions Liao–Lin–Lu–Mao, 2019.

Setting

Let GGG be a simple graph with two labelled sides O\mathcal OO and E\mathcal EE, each containing mmm vertices, and let every vertex have degree Δ\DeltaΔ. An independent set III contains no adjacent pair. For activity λ>0\lambda>0λ>0, the hard-core partition function is

ZG(λ)=∑I independentλ∣I∣.Z_G(\lambda)=\sum_{I\text{ independent}}\lambda^{|I|}.ZG​(λ)=I independent∑​λ∣I∣.

For a vertex set SSS, its vertex boundary ∂S\partial S∂S consists of vertices outside SSS with a neighbor in SSS. The graph is a bipartite (σ,ρ)(\sigma,\rho)(σ,ρ)-expander when every subset of either side with size at most σm\sigma mσm has boundary of size at least ρ\rhoρ times its size. The expansion parameters used here are σ=4log⁡Δ/Δ\sigma=4\log\Delta/\Deltaσ=4logΔ/Δ and ρ=Δ/(4log⁡Δ)−1/2\rho=\Delta/(4\log\Delta)-1/2ρ=Δ/(4logΔ)−1/2. A set is tiny when its size is at most σm\sigma mσm Jenssen–Keevash–Perkins, Definition 21 and §4.4.

The graph G2G^2G2 joins distinct vertices whose distance in GGG is at most two. On either side, a polymer is a nonempty, tiny, G2G^2G2-connected vertex set. Two polymers on the same side are compatible if every pair of their vertices is more than distance two apart in GGG. A polymer γ\gammaγ has weight wγ=λ∣γ∣/(1+λ)∣∂γ∣w_\gamma=\lambda^{|\gamma|}/(1+\lambda)^{|\partial\gamma|}wγ​=λ∣γ∣/(1+λ)∣∂γ∣. Summing products of weights over mutually compatible polymer families gives the partition functions ΞE(G)\Xi^{\mathcal E}(G)ΞE(G) and ΞO(G)\Xi^{\mathcal O}(G)ΞO(G) Jenssen–Keevash–Perkins, §§4.1, 4.4.

“Almost every” means that, for fixed Δ\DeltaΔ, the fraction of labelled Δ\DeltaΔ-regular bipartite graphs satisfying a property tends to one as mmm tends to infinity. The random graph is chosen uniformly from graphs on the two fixed sides Jenssen–Keevash–Perkins, p. 3.

Formalization targets

The first target is the authors' expansion premise: for all sufficiently large Δ\DeltaΔ, almost every graph is a bipartite (4log⁡Δ/Δ,Δ/(4log⁡Δ)−1/2)(4\log\Delta/\Delta,\Delta/(4\log\Delta)-1/2)(4logΔ/Δ,Δ/(4logΔ)−1/2)-expander. The main target combines two deterministic claims on almost every such graph. For every λ>50log⁡2Δ/Δ\lambda>50\log^2\Delta/\Deltaλ>50log2Δ/Δ, write n=2mn=2mn=2m and

Z^G(λ)=(1+λ)m(ΞE(G)+ΞO(G)),ε0=(1+λ)−n/4.\widehat Z_G(\lambda)=(1+\lambda)^m\bigl(\Xi^{\mathcal E}(G)+\Xi^{\mathcal O}(G)\bigr),\qquad \varepsilon_0=(1+\lambda)^{-n/4}.ZG​(λ)=(1+λ)m(ΞE(G)+ΞO(G)),ε0​=(1+λ)−n/4.

First, Z^G\widehat Z_GZG​ is an ε0\varepsilon_0ε0​-relative approximation in the paper's exponential convention:

e−ε0Z^G(λ)≤ZG(λ)≤eε0Z^G(λ).e^{-\varepsilon_0}\widehat Z_G(\lambda)\le Z_G(\lambda)\le e^{\varepsilon_0}\widehat Z_G(\lambda).e−ε0​ZG​(λ)≤ZG​(λ)≤eε0​ZG​(λ).

Second, each of the two polymer models obeys the Kotecký–Preiss condition with decay g(γ)=∣γ∣Δlog⁡(1+λ)/(10log⁡Δ)g(\gamma)=|\gamma|\Delta\log(1+\lambda)/(10\log\Delta)g(γ)=∣γ∣Δlog(1+λ)/(10logΔ):

∑γ′ incompatible with γwγ′e∣γ′∣+g(γ′)≤∣γ∣.\sum_{\gamma'\text{ incompatible with }\gamma}w_{\gamma'}e^{|\gamma'|+g(\gamma')}\le |\gamma|.γ′ incompatible with γ∑​wγ′​e∣γ′∣+g(γ′)≤∣γ∣.

The same graph event must support every activity above the threshold; the activity quantifier is inside the limiting probability. The intermediate targets are Bassalygo's entropy criterion, the high probability expansion lemma, the tiny-side lemma for independent sets, the approximation lemma at its weaker activity threshold, a published binomial-entropy estimate, and the per-vertex polymer sum bound Jenssen–Keevash–Perkins, Theorem 22, Lemmas 23–25 and §4.5.

Significance

The reduction expresses a weighted independent-set count through two finite polymer partition functions whose weights decay under the same expansion condition. The KP bound supplies the convergence condition used by the paper's general polymer approximation theorem. Together these are the mathematical inputs to its algorithmic Theorem 2 Jenssen–Keevash–Perkins, Theorems 2, 8.

A formal proof would make the expansion assumptions, distance convention, graph probability limit, and error normalization explicit. The target is a formalization of a published result, not an open conjecture. The proposal contains statements and definitions; their proofs remain to be supplied. The binomial estimate is already machine checked on Prove2Me and is included by reference.

Difficulty

An independent set can occupy vertices on both sides of a bipartite graph. Merely taking one side as the dominant phase double counts sets that are sparse on both sides. The tiny-side lemma and the polymer decomposition identify and control that overlap. The random graph step also requires an expansion estimate strong enough to apply to every small subset of either side, rather than only to a typical subset Jenssen–Keevash–Perkins, §§4.4–4.5.

Formalization scope

Lean uses finite simple graphs on Fin(m)⊕Fin(m)\mathrm{Fin}(m)\oplus\mathrm{Fin}(m)Fin(m)⊕Fin(m), with the left and right summands serving as the odd and even sides. This fixes a labelled bipartition for the uniform random graph model. Distances use extended natural graph distance, so vertices in different components have infinite distance. Polymer families include the empty family, and graph connectivity excludes the empty polymer. Sizes are natural numbers; thresholds involving logarithms are real inequalities. Binary entropy in the cited expansion theorem is measured in bits, while Mathlib's binary entropy is converted from natural units. The approximation uses exponential relative error exactly as in Definition 10 of the paper.

The FPTAS, sampling algorithm, and running-time assertion are outside this mission because Mathlib has no matching machine and cost model. The goal formalizes both the partition-function reduction and its KP condition, the mathematical inputs consumed by the paper's algorithm theorem. It quantifies over all activities above the threshold within the almost-every event, preventing a separate favorable graph set for each activity. The fixed-sides model and the explicit domain of Bassalygo's entropy formula are disclosed representation choices. The tiny-side milestone (Lemma 24) carries a graph-size threshold depending on the degree: as printed, for every graph size, it fails on small dense graphs (the complete bipartite graph minus disjoint complete blocks), and the main result uses it only as the graph size grows. Reusable contributions include finite graph expansion lemmas, polymer partition-function identities, and bounds on connected vertex sets.

Selected references

  • M. Jenssen, P. Keevash and W. Perkins, Algorithms for #BIS-hard problems on expander graphs, SIAM Journal on Computing 49(4), 2020, 681–710. DOI.
  • L. Bassalygo, Asymptotically optimal switching circuits, Problems of Information Transmission 17(3), 1981, 206–211. Publisher record.
  • C. Liao, J. Lin, P. Lu and Z. Mao, Counting independent sets and colorings on random regular bipartite graphs, APPROX/RANDOM, 2019. DOI.
8 thms2 active usersReviewed
Combinatorics·Captain: mikedeng1

Cycle Lengths in Expanding Graphs 2: Every α-Expander on n Vertices Has Ω(α³/log(1/α))·n Distinct Cycle LengthsResearch Paper

Motivation

An expander graph has many edges leaving each small vertex set. This condition is useful when a graph must remain connected after local losses, and it also constrains the cycles that can occur. Knowing that one long cycle exists leaves open how many different lengths its cycles can have. For example, a bipartite graph has no odd cycle, so expansion alone does not make every length available. The question here is quantitative: how large must the set of available lengths be in terms of the graph's order and expansion parameter? Friedman and Krivelevich answered it with a lower bound whose dependence on the expansion parameter is polynomial up to a logarithm. Their Theorem 2 is the target of this mission. Friedman and Krivelevich, Cycle lengths in expanding graphs, pp. 2–3

The same paper separates two questions. Theorem 1 says that cycle lengths are well spread over a range, but its resulting lower bound on their number has a much smaller dependence on the expansion parameter. Theorem 2 obtains a stronger count without claiming that those lengths form an interval. This distinction matters: counting many available lengths and prescribing an arbitrary length are different demands on a graph. The paper's third theorem studies the latter demand under a different graph condition. Friedman and Krivelevich, pp. 2–4

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite simple graph with n=∣V∣n=|V|n=∣V∣ vertices. A simple graph has no loops or parallel edges. For a vertex set U⊆VU\subseteq VU⊆V, its external neighborhood is

NG(U)={v∈V∖U:v is adjacent to some vertex in U}.N_G(U)=\{v\in V\setminus U:\text{$v$ is adjacent to some vertex in $U$}\}.NG​(U)={v∈V∖U:v is adjacent to some vertex in U}.

The graph is an α\alphaα-expander if α>0\alpha>0α>0 and, for every U⊆VU\subseteq VU⊆V with ∣U∣≤⌈n/2⌉|U|\le\lceil n/2\rceil∣U∣≤⌈n/2⌉,

∣NG(U)∣≥α∣U∣.|N_G(U)|\ge\alpha|U|.∣NG​(U)∣≥α∣U∣.

The ceiling is part of the paper's formal Definition on p. 2. It matters at odd nnn, where the largest tested set has (n+1)/2(n+1)/2(n+1)/2 vertices. A related (k,α)(k,\alpha)(k,α)-expander uses the same inequality for every UUU with ∣U∣≤k|U|\le k∣U∣≤k; here kkk may be real and the size of UUU is still an integer. This broader notion appears in the supporting long-path lemma. Friedman and Krivelevich, pp. 1–2, 4

A cycle is a closed path with no repeated vertex apart from its starting and ending vertex. Its length is its number of edges, hence at least three. Write

L(G)={ℓ∈N:G contains a cycle of length ℓ}.L(G)=\{\ell\in\mathbb N:\text{$G$ contains a cycle of length $\ell$}\}.L(G)={ℓ∈N:G contains a cycle of length ℓ}.

The quantity ∣L(G)∣|L(G)|∣L(G)∣ counts lengths, rather than the number of cycles. Two cycles of the same length contribute just one element. Since a simple cycle visits distinct vertices, L(G)L(G)L(G) is finite and its cardinality is an ordinary natural number. Friedman and Krivelevich, p. 2

Formalization targets

The goal is the polynomial lower bound from Theorem 2. There is one absolute constant c>0c>0c>0 such that, for each 0<α≤10<\alpha\le10<α≤1 and every sufficiently large nnn, every α\alphaα-expander GGG on nnn vertices satisfies

∣L(G)∣≥c α3log⁡2(2/α) n.|L(G)|\ge c\,\frac{\alpha^3}{\log_2(2/\alpha)}\,n.∣L(G)∣≥clog2​(2/α)α3​n.

The paper writes Ω(α3/log⁡(1/α))n\Omega(\alpha^3/\log(1/\alpha))nΩ(α3/log(1/α))n. The denominator above records the endpoint convention used in this mission; it agrees up to an absolute factor for α≤1/2\alpha\le1/2α≤1/2. The constant ccc is uniform across all α\alphaα, graph orders, and graphs. The threshold for nnn may depend on α\alphaα. Friedman and Krivelevich, p. 3, Theorem 2

Four intermediate targets reproduce claims used by the paper. Lemma 2.5 gives a path of a prescribed length from a prescribed start vertex when every kkk-vertex set has enough external neighbors and the start vertex's component is large enough. In the proof of Theorem 2, deleting ⌊αn/4⌋\lfloor\alpha n/4\rfloor⌊αn/4⌋ vertices leaves a path of length ⌈αn/4⌉−1\lceil\alpha n/4\rceil-1⌈αn/4⌉−1. The chosen integer

k(α)=2⌈log⁡2(3/α)log⁡2(1+α/2)⌉+1k(\alpha)=2\left\lceil\frac{\log_2(3/\alpha)}{\log_2(1+\alpha/2)}\right\rceil+1k(α)=2⌈log2​(1+α/2)log2​(3/α)​⌉+1

has order at most log⁡2(2/α)/α\log_2(2/\alpha)/\alphalog2​(2/α)/α with an absolute constant. Finally, the paper's set X0X_0X0​ of path vertices connected to the deleted set by short paths satisfies ∣X0∣≥α2n/12|X_0|\ge\alpha^2n/12∣X0​∣≥α2n/12 once the rounding-sensitive small cases are excluded. These are the milestone statements, in the order in which they enter the target. Friedman and Krivelevich, pp. 6, 13–14

Significance

The theorem guarantees linearly many distinct cycle lengths for every fixed positive expansion parameter, and gives an explicit polynomial-scale dependence on that parameter. A long cycle alone proves only that one length occurs. Even a statement placing cycles near many target lengths can lose substantial strength in its dependence on α\alphaα. The count in Theorem 2 supplies a separate guarantee for graph families where the expansion parameter varies. It does not imply that every length in an interval occurs. Friedman and Krivelevich, p. 3

The result is proved in the cited paper. The formalization target is a machine-checked proof of that known statement and its selected intermediate claims. The reusable output includes an exact interface for external vertex expansion, simple cycle lengths, and the prescribed-start path lemma. The local Lean statements currently compile as open targets; they do not supply machine-checked proofs of the claims.

Difficulty

Expansion gives large neighborhoods and forces connectedness, but these facts do not by themselves distinguish cycle lengths. A tempting first argument finds a long cycle and then treats its subpaths as shorter cycles; the subpaths do not close into cycles without additional edges. A second temptation is to count cycles rather than lengths: a graph may have many cycles while all of them share only a few lengths. The central challenge is to obtain a family of cycles whose lengths can be certified as different while keeping the graph's expansion parameter visible in the count. The proof also has integer rounding at the scale αn\alpha nαn, where claims written for large graphs can fail on small ones. Friedman and Krivelevich, pp. 13–14

Formalization scope

Graphs are SimpleGraph (Fin n) for the goal and proof-specific milestones; the general path lemma allows any finite vertex type. Vertex sets are Set V, with Set.ncard for finite cardinality. An α\alphaα-expander tests every set up to ⌈n/2⌉=(n+1)/2\lceil n/2\rceil=(n+1)/2⌈n/2⌉=(n+1)/2 in natural-number arithmetic. Cycles use SimpleGraph.Walk.IsCycle, and paths use Walk.IsPath; their length is the number of edges. This rules out a two-edge closed walk being counted as a cycle. The formal development also uses the explicit k(α)k(\alpha)k(α) and the paper's X0X_0X0​, including its condition on the internal vertices of a short connecting path.

All logarithms are binary. The shifted log⁡2(2/α)\log_2(2/\alpha)log2​(2/α) replaces log⁡2(1/α)\log_2(1/\alpha)log2​(1/α) only inside asymptotic bounds, because the latter vanishes at α=1\alpha=1α=1. The constant in the lower bound precedes α\alphaα in the quantifier order. The goal includes n≥n0(α)n\ge n_0(\alpha)n≥n0​(α), a necessary disclosed addition: the two-vertex complete graph is an α\alphaα-expander for α≤1\alpha\le1α≤1 but has no cycle. The X0X_0X0​ milestone adds αn≥12\alpha n\ge12αn≥12 to make the paper's suppressed rounding valid; otherwise the initial vertex set can be empty. The deletion milestone states the path claim for any deleted set of the indicated size, since the displayed estimate uses only that size. A formal proof needs finite graph paths, component sizes, vertex neighborhoods, and cardinality estimates. Contributions to these general lemmas are useful beyond this mission. No claim is discharged by letting an asymptotic constant depend on α\alphaα, by counting closed walks, or by exploiting division by zero.

Selected references

  • L. Friedman and M. Krivelevich, Cycle lengths in expanding graphs, arXiv:1912.11011v2, 2020; published in Combinatorica 41 (2021). arXiv; DOI.
7 thms2 active usersReviewed
CombinatoricsDiscrete Geometry·Captain: mikedeng1

Twin-width I: Tractable FO Model Checking 2: Every Subgraph of a Unit d-Dimensional Ball Graph with Clique Number k Has Twin-width at Most (3⌈√d⌉)^d·kResearch Paper

Motivation

Twin-width is a graph parameter introduced by Bonnet, Kim, Thomassé and Watrigant in Twin-width I: Tractable FO Model Checking (J. ACM 69(1), Article 3, 2021; arXiv:2004.14789). Its main theorem says that first-order model checking is fixed-parameter tractable on every class of graphs of bounded twin-width, provided a witnessing contraction sequence is given. This turns every class shown to have bounded twin-width into a class on which a whole logic becomes algorithmically tractable, so the question "which natural classes have bounded twin-width?" has direct algorithmic content.

Section 4 of the paper answers that question for several classes. This mission covers the geometric part of Section 4: ddd-dimensional grids, grids with diagonals, and unit ddd-dimensional ball graphs with bounded clique number. Geometric intersection graphs are a standard source of hard instances in parameterized algorithms. Width parameters designed for dense graphs (rank-width, clique-width, boolean-width) are unbounded already on the planar n×nn\times nn×n grid, whereas twin-width stays bounded on grids of every fixed dimension and transfers from there to ball graphs.

Setting

All graphs are finite and simple. Two vertex sets X,YX,YX,Y of a graph GGG are homogeneous if every pair (x,y)∈X×Y(x,y)\in X\times Y(x,y)∈X×Y is an edge or no pair is. For a partition P\mathcal PP of V(G)V(G)V(G), the red degree of a part XXX is the number of other parts not homogeneous to XXX; P\mathcal PP is a ddd-partition if every red degree is at most ddd. The graph GGG has twin-width at most ddd, tww⁡(G)≤d\operatorname{tww}(G)\le dtww(G)≤d, if there is a sequence P0,…,PN\mathcal P_0,\dots,\mathcal P_NP0​,…,PN​ of ddd-partitions of V(G)V(G)V(G) that starts with the partition into singletons, ends with at most one part, and in which each partition arises from the previous one by merging two parts. This is the paper's partition form of its definition by contractions of trigraphs (graphs with black and red edges, §3).

For the all-red trigraph Hr=(V,∅,E(H))H^r=(V,\emptyset,E(H))Hr=(V,∅,E(H)), in which every edge of HHH is red, the paper's contraction rule never produces a black edge, and two parts are red-adjacent exactly when some edge of HHH joins them. The red twin-width tww⁡(Hr)≤d\operatorname{tww}(H^r)\le dtww(Hr)≤d is defined as above with this adjacency.

For d,n≥1d,n\ge 1d,n≥1, the ddd-dimensional nnn-grid PndP^d_nPnd​ has vertex set [n]d[n]^d[n]d, with x∼yx\sim yx∼y iff ∑i∣xi−yi∣=1\sum_{i}|x_i-y_i|=1∑i​∣xi​−yi​∣=1. The grid with diagonals Kn,dK_{n,d}Kn,d​ has the same vertices, with distinct x∼yx\sim yx∼y iff max⁡i∣xi−yi∣≤1\max_i|x_i-y_i|\le 1maxi​∣xi​−yi​∣≤1, and Kn,drK^r_{n,d}Kn,dr​ is its all-red trigraph; RndR^d_nRnd​ is the all-red trigraph of PndP^d_nPnd​.

Given centres c(v)∈Rdc(v)\in\mathbb R^dc(v)∈Rd for vvv in a finite set VVV, the unit ball graph GGG joins distinct u,vu,vu,v iff the closed unit balls around c(u)c(u)c(u) and c(v)c(v)c(v) meet, that is ∥c(u)−c(v)∥2≤2\|c(u)-c(v)\|_2\le 2∥c(u)−c(v)∥2​≤2. Its clique number is the largest size of a set of pairwise adjacent vertices. A subgraph of GGG is obtained by deleting vertices and edges.

Formalization targets

Goal: Theorem 4.5 (p. 3:15)

For every d,kd,kd,k, every subgraph HHH of a unit ddd-dimensional ball graph GGG with clique number at most kkk satisfies

tww⁡(H)≤(3⌈d ⌉)d k.\operatorname{tww}(H)\le (3\lceil\sqrt d\,\rceil)^d\,k .tww(H)≤(3⌈d​⌉)dk.

Milestones, in the order the paper uses them

  1. Red paths have twin-width at most 222 (§3, p. 3:12).
  2. Induced subgraphs do not increase twin-width, for graphs and for all-red trigraphs (§4.1, p. 3:12).
  3. tww⁡(Rnd)≤3d\operatorname{tww}(R^d_n)\le 3dtww(Rnd​)≤3d for positive d,nd,nd,n (the statement proved inside Theorem 4.3, p. 3:15).
  4. Theorem 4.3: tww⁡(Pnd)≤3d\operatorname{tww}(P^d_n)\le 3dtww(Pnd​)≤3d for positive d,nd,nd,n.
  5. Every subgraph of PndP^d_nPnd​ has twin-width at most 3d3d3d (p. 3:15).
  6. Lemma 4.4: every subgraph of Kn,drK^r_{n,d}Kn,dr​ has twin-width at most 2(3d−1)2(3^d-1)2(3d−1).
  7. In a unit ball graph without a clique of size k+1k+1k+1, each half-open cell of side 2/d2/\sqrt d2/d​ contains at most kkk centres (proof of Theorem 4.5, p. 3:16).

A companion item, Theorem 4.1 (p. 3:13, graph case), states that adding one vertex joined to an arbitrary set at most changes twin-width from ttt to 2(t+1)2(t+1)2(t+1).

Significance

The goal theorem places unit ball graphs of bounded clique number, in every fixed dimension, among the classes of bounded twin-width. Combined with the paper's main theorem, and with the paper's remark that a contraction sequence can be computed in polynomial time from a geometric representation, it gives fixed-parameter tractable first-order model checking on these graphs. The clique-number hypothesis is needed: unit disk graphs without a clique bound have unbounded twin-width (p. 3:16), as shown in Twin-width II. Theorem 4.3 and Lemma 4.4 are of independent use as twin-width bounds for grids, and the subgraph statements show that grids are an exception to the general fact that bounded twin-width is not preserved by taking (non-induced) subgraphs (p. 3:13).

All results here are proved in the paper. To our knowledge none of them, and no notion of twin-width with red degree, is machine-checked in Lean or Mathlib. The work is to formalize the paper's arguments: the parallel contraction of grid layers, the supercell contraction, and the geometric counting.

Difficulty

Twin-width is not monotone under subgraphs, so the obvious route, "a unit ball graph of bounded clique number is a subgraph of a grid-like graph, hence has bounded twin-width", fails as stated. The paper repairs it by working with all-red trigraphs throughout: for an all-red trigraph, removing edges can only remove red adjacencies, and this is why Theorem 4.3 is proved in its red form and Lemma 4.4 is stated for Kn,drK^r_{n,d}Kn,dr​. A second difficulty is bookkeeping: the bound 3d3d3d comes from contracting the nnn layers of PndP^d_nPnd​ "in parallel" and tracking the red degree of every vertex at each step, and Lemma 4.4 is asserted "by the arguments of Theorem 4.3" without a separate proof, so its argument has to be reconstructed.

Formalization scope

Lean represents twin-width by the predicates TwinWidthLE G d and RedTwinWidthLE H d on finite vertex types: an explicit sequence of Finpartitions of the vertex set from ⊥ (singletons) to at most one part, consecutive partitions differing by one merge of two distinct parts. Twin-width is never a number, so no infimum over an empty set enters. [n]d[n]^d[n]d is Fin d → Fin n; centres live in EuclideanSpace ℝ (Fin d) with the Euclidean distance; balls are closed of radius 111, so centres at distance exactly 222 are adjacent. ⌈d ⌉\lceil\sqrt d\,\rceil⌈d​⌉ is the natural-number ceiling of the real square root. "Clique number kkk" is encoded as CliqueFree (k+1) (clique number at most kkk), which is the same theorem because the bound is monotone in kkk; subgraphs are a finite vertex set SSS together with a graph H≤G[S]H\le G[S]H≤G[S] on it.

The hypotheses 1≤d1\le d1≤d, 1≤n1\le n1≤n appear in Theorem 4.3, its red form, and the grid-subgraph statement, as on the page ("positive integers"). No other hypothesis is added; Lemma 4.4, the cell bound and the goal hold as stated also for d=0d=0d=0 or n=0n=0n=0. The red-degree count uses "some edge joins the two parts", which is derived from the paper's contraction rule for trigraphs without black edges and is used only for all-red trigraphs; the twin-width of graphs counts non-homogeneous parts. This rules out the trivializing readings: a twin-width defined as an infimum that defaults to 000, a sequence allowed to jump to one part, homogeneity counted as adjacency for graphs, or a bound depending on the graph. The polynomial-time sentence of Theorem 4.5 is not formalized.

Reusable infrastructure: the partition-form twin-width predicates (shared with the other missions of this series), grids and grids with diagonals in coordinates, and unit ball graphs. Contributions are welcome on the general lemmas (induced subgraphs, red-to-black transfer, small explicit sequences) as well as on the grid theorems themselves.

Selected references

  • É. Bonnet, E. J. Kim, S. Thomassé, R. Watrigant, Twin-width I: Tractable FO Model Checking, Journal of the ACM 69(1), Article 3, 2021. https://doi.org/10.1145/3486655 (preprint: https://arxiv.org/abs/2004.14789)
  • É. Bonnet, C. Geniet, E. J. Kim, S. Thomassé, R. Watrigant, Twin-width II: small classes, Combinatorial Theory 2(2), 2022. https://arxiv.org/abs/2006.09877
10 thms2 active usersReviewed
Combinatorics·Captain: mikedeng1

Twin-width I: Tractable FO Model Checking 1: Every Graph with Boolean-width k Has Twin-width at Most 2^(k+1) − 1Research Paper

Motivation

Twin-width was introduced by Bonnet, Kim, Thomassé and Watrigant (J. ACM 69(1), Article 3, 2021) as a graph parameter under which first-order model checking is fixed-parameter tractable when a contraction sequence is given. Its interest lies in the classes it covers at once: planar and, more generally, KtK_tKt​-minor-free graphs, bounded-degree grids, unit ball graphs of bounded clique number, and the classes of bounded rank-width or clique-width. These are the classical dense-graph width measures, for which FO (indeed MSO1_11​) model checking was already known to be tractable through Courcelle–Makowsky–Rotics.

This mission formalizes the link to the dense width measures, stated in the paper through boolean-width (Bui-Xuan, Telle and Vatshelle, Theor. Comput. Sci. 412(39), 2011). It is known that boolw⁡(G)≤modw⁡(G)≤cw⁡(G)≤2rw⁡(G)+1−1\operatorname{boolw}(G)\le\operatorname{modw}(G)\le\operatorname{cw}(G)\le 2^{\operatorname{rw}(G)+1}-1boolw(G)≤modw(G)≤cw(G)≤2rw(G)+1−1 (Vatshelle's PhD thesis, Chapter 4), so a bound on twin-width in terms of boolean-width shows that bounded twin-width captures bounded rank-width, clique-width, module-width and boolean-width simultaneously.

Setting

Let GGG be a finite simple graph on a vertex set VVV.

Two vertex sets X,Y⊆VX,Y\subseteq VX,Y⊆V are homogeneous if either all pairs (x,y)∈X×Y(x,y)\in X\times Y(x,y)∈X×Y are edges or none is. For a partition P\mathcal PP of VVV, a part XXX has red degree equal to the number of other parts Y≠XY\ne XY=X that are not homogeneous to XXX; P\mathcal PP is a ddd-partition if every part has red degree at most ddd. The graph GGG has twin-width at most ddd, written tww⁡(G)≤d\operatorname{tww}(G)\le dtww(G)≤d, if there are partitions P0,…,PN\mathcal P_0,\dots,\mathcal P_NP0​,…,PN​ of VVV with P0\mathcal P_0P0​ the partition into singletons, PN\mathcal P_NPN​ having at most one part, each Pi+1\mathcal P_{i+1}Pi+1​ obtained from Pi\mathcal P_iPi​ by merging two parts, and every Pi\mathcal P_iPi​ a ddd-partition. This is the paper's partition form (p. 3:32) of its definition through trigraphs and red edges (pp. 3:11–3:12).

A decomposition tree of GGG is a rooted binary tree TTT, every internal node having exactly two children, whose leaves are in one-to-one correspondence with VVV. Each edge eee of TTT lies above exactly one rooted subtree and induces a cut Pe=(Ae,Be)P_e=(A_e,B_e)Pe​=(Ae​,Be​) of VVV: AeA_eAe​ is the set of leaf labels of that subtree and Be=V∖AeB_e=V\setminus A_eBe​=V∖Ae​. For disjoint A,BA,BA,B write

νG(A,B)=∣{N(S)∩B:S⊆A}∣,\nu_G(A,B)=\bigl|\{N(S)\cap B: S\subseteq A\}\bigr|,νG​(A,B)=​{N(S)∩B:S⊆A}​,

the number of different neighbourhoods in BBB of subsets of AAA, the empty subset included. The boolean-width of the cut is log⁡2νG(A,B)\log_2\nu_G(A,B)log2​νG​(A,B), the boolean-width of TTT is its maximum over the edges of TTT, and boolw⁡(G)\operatorname{boolw}(G)boolw(G) is the minimum over all decomposition trees. In Lean, nbhdCount G A B is νG(A,B)\nu_G(A,B)νG​(A,B), and BoolWidthLE G N states that some decomposition tree has νG(Ae,Be)≤N\nu_G(A_e,B_e)\le NνG​(Ae​,Be​)≤N on every edge, so boolw⁡(G)≤k\operatorname{boolw}(G)\le kboolw(G)≤k is BoolWidthLE G (2^k).

Formalization targets

Goal: Theorem 4.2 (p. 3:14)

Every graph with boolean-width kkk has twin-width at most 2k+1−12^{k+1}-12k+1−1. With N=2kN=2^{k}N=2k, the integer neighbourhood count,

BoolWidthLE⁡(G,N) ⟹ tww⁡(G)≤2N−1.\operatorname{BoolWidthLE}(G,N)\ \Longrightarrow\ \operatorname{tww}(G)\le 2N-1 .BoolWidthLE(G,N) ⟹ tww(G)≤2N−1.

Milestones (proof of Theorem 4.2, p. 3:14)

  1. A graph with at most NNN vertices has tww⁡(G)≤N\operatorname{tww}(G)\le Ntww(G)≤N.
  2. A rooted binary tree with more than 2N2N2N leaves (N≥1N\ge1N≥1) has a rooted subtree below some edge with between N+1N+1N+1 and 2N2N2N leaves.
  3. If νG(A,B)≤N\nu_G(A,B)\le NνG​(A,B)≤N and ∣A∣≥N+1|A|\ge N+1∣A∣≥N+1, two distinct vertices of AAA have the same neighbourhood in BBB.
  4. If (A,B)(A,B)(A,B) partitions VVV, ∣A∣≤2N|A|\le 2N∣A∣≤2N, and u≠v∈Au\ne v\in Au=v∈A have the same neighbourhood in BBB, then contracting uuu and vvv yields a (2N−1)(2N-1)(2N−1)-partition.

Significance

The theorem places the classes of bounded boolean-width, and therefore of bounded clique-width and rank-width, inside the classes of bounded twin-width, with a single-exponential bound. Combined with the paper's model-checking algorithm, it recovers tractability of FO model checking on these classes once a contraction sequence is computed, and it situates twin-width as a common generalization of the dense width measures and of sparse classes such as KtK_tKt​-minor-free graphs.

No machine-checked version of twin-width or boolean-width is known to exist; neither notion is in Mathlib. The mission produces Lean definitions of twin-width in partition form, of decomposition trees and of the neighbourhood count behind boolean-width, together with the four combinatorial steps the paper's proof states explicitly. The definitions are shared, in identical form, with the other missions of this series on the same paper.

Difficulty

The first three milestones are elementary: a counting bound on contraction sequences, a walk down a binary tree, and a pigeonhole argument. The fourth is a direct count of non-homogeneous parts.

The difficulty is the iteration. After the first contraction the leaves of the updated tree are contracted vertices, so "the same neighbourhood in BeB_eBe​" must be read as homogeneity to all of BeB_eBe​ for a set of vertices, and the boolean-width bound of the cuts must be re-established for sets rather than single vertices. The paper asserts the invariant that the contractions never create a red component of more than 2k+12^{k+1}2k+1 vertices without proving it. A complete proof of the goal therefore has to supply an invariant, carried through all ∣V∣−1|V|-1∣V∣−1 contractions, that keeps every part's red degree below 2N−12N-12N−1. The naive repetition of the first step does not by itself control red edges between previously contracted parts.

Formalization scope

Twin-width is the predicate TwinWidthLE G d, never a numerical infimum, so a junk value of an empty or unbounded infimum cannot make the goal trivial. A contraction sequence must merge exactly two parts at each step and start from singletons; it cannot jump to one part. Homogeneity counts all pairs between two parts, in both of its alternatives. Boolean-width is bounded on every edge of the decomposition tree, not only at the root, and the count runs over neighbourhoods of all subsets of AeA_eAe​, not only of single vertices. The bound 2N−12N-12N−1 depends only on NNN, not on the graph.

Boolean-width is kept as the integer count N=2boolw⁡(G)N=2^{\operatorname{boolw}(G)}N=2boolw(G); since boolw⁡(G)=k\operatorname{boolw}(G)=kboolw(G)=k means that the least admissible count is 2k2^k2k, the Lean goal for this NNN is exactly the printed theorem. Subtractions 2N−12N-12N−1 are in N\mathbb NN and exact where used. Graphs live on a Fintype vertex type with decidable equality and decidable adjacency. A graph on the empty vertex set has no decomposition tree; the goal is silent about it, and it has twin-width 000. No hypothesis is added to any statement relative to the page except N≥1N\ge 1N≥1 in milestone 2, which the page has implicitly (2k≥12^k\ge12k≥1) and without which the conclusion is impossible.

The definitions needed are the setting module TwinWidthI.BoolWidth.Setting: merge steps, homogeneity, ddd-partitions, twin-width, decomposition trees (DTree), the count νG\nu_GνG​ and BoolWidthLE. Proofs of the milestones, a proof of the goal, and lemmas towards the iteration invariant are welcome.

Selected references

  • É. Bonnet, E. J. Kim, S. Thomassé and R. Watrigant, Twin-width I: Tractable FO Model Checking, J. ACM 69(1), Article 3, 2021. https://doi.org/10.1145/3486655
  • B.-M. Bui-Xuan, J. A. Telle and M. Vatshelle, Boolean-width of graphs, Theoretical Computer Science 412(39), 5187–5204, 2011. https://doi.org/10.1016/j.tcs.2011.05.022
  • M. Vatshelle, New Width Parameters of Graphs, PhD thesis, University of Bergen, 2012 (reference [39] of Bonnet et al.; source of the chain boolw ≤ modw ≤ cw ≤ 2^{rw+1} − 1).
6 thms2 active usersReviewed
Combinatorics·Captain: mikedeng1

Steinberg's Conjecture is false: There Exists a Planar Graph with No Cycles of Length Four or Five That Is Not 3-ColorableResearch Paper

Motivation

Deciding whether a planar graph is 3-colorable is NP-complete, so a large part of the theory of planar graph coloring looks for simple local conditions that guarantee 3-colorability. Grötzsch's theorem (1959) is the model result: every triangle-free planar graph is 3-colorable. Forbidding triangles is not the only possibility, and in 1976 Richard Steinberg asked whether forbidding cycles of two other lengths suffices.

Steinberg's Conjecture (1976). Every planar graph with no cycles of length four or five is 3-colorable.

The conjecture drove four decades of work on planar coloring. It was surveyed by Borodin as one of the central problems of the area (Borodin 2013, §7), and it was one of the six graph theory problems ranked with the highest importance in the Open Problem Garden.

Timeline:

  • 1976. Steinberg poses the conjecture (Steinberg 1993 records its history).
  • 1991. Erdős suggests the relaxation: find the smallest kkk such that every planar graph with no cycles of lengths 4,…,k4, \dots, k4,…,k is 3-colorable. Abbott and Zhou show k≤11k \le 11k≤11 (Ars Combin. 32, 1991).
  • 1995. Sanders and Zhao improve this to k≤9k \le 9k≤9 (Graphs Combin. 11, 1995).
  • 2005. Borodin, Glebov, Raspaud and Salavatipour prove k≤7k \le 7k≤7: planar graphs without cycles of lengths 4 to 7 are 3-colorable (JCTB 93).
  • 2007. Voigt shows that the list-coloring analogue fails: there is a planar graph without 4- and 5-cycles that is not 3-choosable (Discrete Math. 307, 2007).
  • 2016. Cohen-Addad, Hebdige, Král', Li and Salgado construct a counterexample, so the conjecture is false (arXiv:1604.05108; J. Combin. Theory Ser. B, 2017).

Setting

All graphs are finite and simple. A cycle of length kkk in a graph HHH is a sequence of k≥3k \ge 3k≥3 pairwise distinct vertices v0,…,vk−1v_0, \dots, v_{k-1}v0​,…,vk−1​ with vivi+1v_i v_{i+1}vi​vi+1​ an edge for each iii (indices mod kkk). HHH has no cycles of length four or five if no cycle of HHH has length 444 or 555; triangles and cycles of length at least six are allowed. A 3-coloring of HHH assigns one of three colors to each vertex so that adjacent vertices get different colors, and HHH is 3-colorable if it has one. The distance dH(u,v)d_H(u,v)dH​(u,v) is the number of edges of a shortest uuu–vvv path. HHH is planar if it can be drawn in the plane without crossing edges.

The paper builds its counterexample in three steps, from three concrete graphs given as drawings.

  1. G1G_1G1​ (Figure 1) has 15 vertices a,…,oa, \dots, oa,…,o and 23 edges. Its contact vertices are aaa, bbb, ccc.
  2. G2G_2G2​ (Figure 2) consists of three copies of G1G_1G1​ glued along contact vertices around a central triangle defdefdef. It has 42 vertices and 72 edges, and its contact vertices are aaa, bbb, ccc.
  3. GGG (Figure 3) consists of four copies of G2G_2G2​ glued at a central vertex aaa, together with the extra vertices bbb, eee, e′e'e′ and twelve extra edges. It has 166 vertices and 300 edges.

In the Lean development these are the definitions G1, G2 and G. The hypothesis "no cycles of length four or five" is the definition NoFourFiveCycle.

Formalization targets

Goal: Theorem 3 (p. 4)

∃ H finite, planar:H has no cycle of length 4 or 5  and  H is not 3-colorable.\exists\, H \text{ finite, planar}:\quad H \text{ has no cycle of length } 4 \text{ or } 5 \ \text{ and } \ H \text{ is not 3-colorable}.∃H finite, planar:H has no cycle of length 4 or 5  and  H is not 3-colorable.

This is the paper's main theorem and exactly the negation of Steinberg's Conjecture. It leaves the witness free, so any counterexample proves it.

Milestones

  1. Lemma 1 (p. 2): G1G_1G1​ has no 4- or 5-cycle; no 3-coloring gives a,b,ca, b, ca,b,c the same color; and d(a,b)=d(a,c)=3d(a,b) = d(a,c) = 3d(a,b)=d(a,c)=3, d(b,c)=4d(b,c) = 4d(b,c)=4.
  2. Lemma 2 (p. 3): G2G_2G2​ has no 4- or 5-cycle; no 3-coloring gives a,b,ca, b, ca,b,c the same color; and the contacts are pairwise at distance 444.
  3. GGG has no cycle of length four or five (proof of Theorem 3, p. 5).
  4. GGG has no 3-coloring (proof of Theorem 3, p. 5).
  5. GGG is planar (§2, p. 2, and Figure 3).
  6. G1G_1G1​ and G2G_2G2​ are planar (§2, p. 2, and Figures 1–2).

Milestones 3–5 together give the goal, with GGG as the witness.

Significance

The result. Theorem 3 settles a 40-year-old problem negatively. It shows that forbidding 4- and 5-cycles is not a local condition that forces 3-colorability, which redirects the Erdős relaxation: the threshold kkk for "no cycles of lengths 4,…,k4, \dots, k4,…,k" lies in {6,7}\{6, 7\}{6,7}, since k≤7k \le 7k≤7 is known and k=5k = 5k=5 is now excluded. The paper notes that the same graph also disproves the Strong Bordeaux Conjecture and the Novosibirsk 3-Color Conjecture.

Formalizing it. The theorem is proved, and the proof is short. As far as is known it has no machine-checked version. Every claim in the proof is a finite statement about explicit graphs, so the work here is to verify the three gadgets and the gluing arguments in Lean, including planarity of an explicit 166-vertex graph, which is asserted in the paper only through its figures. Planarity, cycle structure and coloring of a concrete graph of this size are a useful benchmark for computational verification in Lean.

Difficulty

The obvious approach is brute force: GGG is a fixed finite graph, so in principle each claim can be decided. In practice decide is not workable at this size. Excluding 4- and 5-cycles means quantifying over walks of a 166-vertex graph. Non-3-colorability is a search over 31663^{166}3166 assignments unless it is reduced, as the paper does, through the gadget lemmas. The planarity witness needs explicit coordinates for 166 points and a check of all (3002)\binom{300}{2}(2300​) pairs of segments, in exact arithmetic. The gluing arguments ("a cycle through two copies has length at least six") are where the paper's distance claims are used. Turning those claims into statements about Mathlib walks on unions of embedded copies is the conceptual work.

Formalization scope

Lean representation and conventions:

  • Graphs are Mathlib SimpleGraphs. G1 lives on Fin 15, with a,…,o↦0,…,14a, \dots, o \mapsto 0, \dots, 14a,…,o↦0,…,14. G2 lives on Fin 42, with a,b,c,d,e,f↦0,…,5a, b, c, d, e, f \mapsto 0, \dots, 5a,b,c,d,e,f↦0,…,5, and G lives on Fin 166, with a,b,c,c′,e,e′,d,f,d′,f′↦0,…,9a, b, c, c', e, e', d, f, d', f' \mapsto 0, \dots, 9a,b,c,c′,e,e′,d,f,d′,f′↦0,…,9.
  • Copies. G2 is literally (⨆ k, G1.map (copy2 k)) ⊔ the triangle defdefdef, with explicit embeddings copy2 k : Fin 15 ↪ Fin 42. G is (⨆ k, G2.map (copyG k)) ⊔ the twelve extra edges, with copyG k : Fin 42 ↪ Fin 166. The contact tables are in each definition's statement. A local check confirms the edge counts 23, 72 and 300, and that G agrees edge for edge with an independent transcription of Figure 3.
  • Cycles are Mathlib walks p : H.Walk v v with p.IsCycle, and their length is p.length.
  • 3-colorings are H.Coloring (Fin 3), and 3-colorability is H.Colorable 3. "No 3-coloring gives a,b,ca, b, ca,b,c one color" is ∀ C, ¬ (C a = C b ∧ C b = C c). Colorings in which two of the three share a color are still allowed, as in the paper.
  • Distance is SimpleGraph.dist. Its junk value 000 for unreachable pairs is never hit, because every asserted distance is positive.
  • Planarity is the published straight-line predicate OPG37357.IsPlanar: an injective placement of the vertices in R2\mathbb R^2R2 in which no vertex lies inside an edge segment and edges with four distinct endpoints have disjoint segments. For finite simple graphs it is equivalent to planarity by Fáry's theorem, and it implies topological planarity.
  • The goal's witness is a graph on Fin n. This loses no generality for finite graphs, and the empty graph cannot serve as a trivial witness because it is 3-colorable.

The following trivializing formalizations are ruled out. Closed walks, instead of cycles, would exclude every graph with an edge. Girth is not used: GGG contains triangles. Lemmas 1 and 2 do not assert that aaa and bbb get different colors, which would be false for G1G_1G1​. An infinite or unspecified vertex type is not allowed for the witness.

Welcome contributions: decision procedures or certificates for short cycles and colorings of explicit finite graphs; a straight-line embedding of GGG with rational coordinates; lemmas on cycles and distances in graphs glued from embedded copies, which are reusable for other gadget constructions; and proofs of Lemmas 1–2 by computation or by the paper's arguments.

Selected references

  • V. Cohen-Addad, M. Hebdige, D. Král', Z. Li, E. Salgado, Steinberg's Conjecture is false, J. Combin. Theory Ser. B 122 (2017), 452–456; preprint arXiv:1604.05108v2. https://arxiv.org/abs/1604.05108v2
  • O. V. Borodin, Colorings of plane graphs: A survey, Discrete Math. 313 (2013), 517–539. https://doi.org/10.1016/j.disc.2012.11.011
  • O. V. Borodin, A. N. Glebov, A. Raspaud, M. R. Salavatipour, Planar graphs without cycles of length from 4 to 7 are 3-colorable, J. Combin. Theory Ser. B 93 (2005), 303–331. https://doi.org/10.1016/j.jctb.2004.11.001
  • D. P. Sanders, Y. Zhao, A note on the three color problem, Graphs Combin. 11 (1995), 91–94. (Reference [11] of arXiv:1604.05108v2.)
  • R. Steinberg, The state of the three color problem, in: Quo Vadis, Graph Theory?, Ann. Discrete Math. 55 (1993), 211–248. (Reference [12] of arXiv:1604.05108v2.)
  • M. Voigt, A non-3-choosable planar graph without cycles of length 4 and 5, Discrete Math. 307 (2007), 1013–1015. (Reference [13] of arXiv:1604.05108v2.)
13 thms2 active usersReviewed
Complexity TheoryLinear algebraTheoretical Computer Science·Captain: mikedeng1

Most Tensor Problems Are NP-Hard 2: Graph 3-Colorability Reduces to Tensor Bilinear Feasibility over the Reals and over the Complex NumbersResearch Paper

Motivation

A matrix A∈Cn×nA \in \mathbb C^{n\times n}A∈Cn×n is singular exactly when the pair of linear systems x⊤A=0\mathbf x^\top A = 0x⊤A=0, Ay=0A\mathbf y = 0Ay=0 has a solution with both x\mathbf xx and y\mathbf yy nonzero, and this is decided in polynomial time by Gaussian elimination. The natural analogue for a 3-tensor replaces the two linear systems by three families of bilinear equations, one for each way of contracting the tensor against two vectors. For tensors whose dimensions satisfy the Gelfand–Kapranov–Zelevinsky condition, solvability of that system is equivalent to the vanishing of the hyperdeterminant, and for any tensor it is equivalent to 000 being a singular value in the sense of Lim (2005). The same system appears in multiview geometry (Cohen and Tomasi 1997).

C. J. Hillar and L.-H. Lim, Most Tensor Problems Are NP-Hard (J. ACM 60(6), 2013; arXiv:0911.1393), show that, unlike the matrix case, this feasibility question is NP-hard over R\mathbb RR and over C\mathbb CC (their Theorem 3.7). The proof is a reduction from graph 3-colorability through the polynomial color encoding of a graph, an idea the authors trace to Bayer's thesis (1982) and to Lovász (1994) and De Loera (1995). This mission formalizes that reduction and its consequence for tensor singular values (Theorem 6.2).

Setting

A 3-tensor over a field FFF is an array A=[ ⁣[aijk] ⁣]∈Fl×m×n\mathcal A = [\![a_{ijk}]\!] \in F^{l\times m\times n}A=[[aijk​]]∈Fl×m×n. Its slices are the matrices Ai(j,k)=Bj(i,k)=Ck(i,j)=aijkA_i(j,k) = B_j(i,k) = C_k(i,j) = a_{ijk}Ai​(j,k)=Bj​(i,k)=Ck​(i,j)=aijk​, all read off the same array.

Tensor bilinear feasibility (Problem 3.1) asks, for a rational tensor A∈Ql×m×n\mathcal A \in \mathbb Q^{l\times m\times n}A∈Ql×m×n and F=RF = \mathbb RF=R or C\mathbb CC, whether the system (9)

v⊤Aiw=0 (i≤l),u⊤Bjw=0 (j≤m),u⊤Ckv=0 (k≤n)\mathbf v^\top A_i\mathbf w = 0\ (i \le l),\qquad \mathbf u^\top B_j \mathbf w = 0\ (j \le m),\qquad \mathbf u^\top C_k \mathbf v = 0\ (k \le n)v⊤Ai​w=0 (i≤l),u⊤Bj​w=0 (j≤m),u⊤Ck​v=0 (k≤n)

has a solution u∈Fl\mathbf u \in F^lu∈Fl, v∈Fm\mathbf v\in F^mv∈Fm, w∈Fn\mathbf w \in F^nw∈Fn with each of the three vectors nonzero. In Lean this is TBF A.

The color encoding CGC_GCG​ of a graph G=(V,E)G = (V,E)G=(V,E) on vvv vertices (Definition 2.9) is the set of 4v4v4v quadratics in x1,…,xv,y1,…,yv,zx_1,\dots,x_v,y_1,\dots,y_v,zx1​,…,xv​,y1​,…,yv​,z:

xiyi−z2,yiz−xi2,xiz−yi2,∑j:{i,j}∈E(xi2+xixj+xj2),i=1,…,v.x_iy_i - z^2,\quad y_iz - x_i^2,\quad x_iz - y_i^2,\quad \sum_{j:\{i,j\}\in E}(x_i^2 + x_ix_j + x_j^2),\qquad i = 1,\dots,v .xi​yi​−z2,yi​z−xi2​,xi​z−yi2​,j:{i,j}∈E∑​(xi2​+xi​xj​+xj2​),i=1,…,v.

The tensor AG∈Zl×m×n\mathcal A_G \in \mathbb Z^{l\times m\times n}AG​∈Zl×m×n of the proof of Theorem 3.7 has m=n=2v+1m = n = 2v+1m=n=2v+1 and l=v(2v+5)l = v(2v+5)l=v(2v+5) slices: v(2v+1)v(2v+1)v(2v+1) slices whose bilinear forms are the 2×22\times 22×2 minors of the matrix [v w][\mathbf v\ \mathbf w][v w], then 3v+v3v + v3v+v slices whose bilinear forms restrict, on w=v\mathbf w = \mathbf vw=v, to the polynomials of CGC_GCG​ (Lean: AG G).

The tensor B\mathcal BB of Lemma 3.6 doubles a real tensor A\mathcal AA into R2l×2m×2n\mathbb R^{2l\times 2m\times 2n}R2l×2m×2n through the slices Bi=[Ai00−Ai]B_i = \left[\begin{smallmatrix} A_i & 0 \\ 0 & -A_i\end{smallmatrix}\right]Bi​=[Ai​0​0−Ai​​], Bl+i=[0AiAi0]B_{l+i} = \left[\begin{smallmatrix} 0 & A_i \\ A_i & 0\end{smallmatrix}\right]Bl+i​=[0Ai​​Ai​0​] (Lean: doubling A).

An ℓ2\ell^2ℓ2-singular value σ\sigmaσ of A\mathcal AA (Definition 6.1) admits nonzero u,v,w\mathbf u,\mathbf v,\mathbf wu,v,w with ∑j,kaijkvjwk=σui\sum_{j,k}a_{ijk}v_jw_k = \sigma u_i∑j,k​aijk​vj​wk​=σui​ and the two analogous equations; an ℓ3\ell^3ℓ3-singular value has σui2\sigma u_i^2σui2​, σvj2\sigma v_j^2σvj2​, σwk2\sigma w_k^2σwk2​ on the right. At σ=0\sigma = 0σ=0 both reduce to (9).

Formalization targets

Goal: Theorem 3.7

3-COLORABILITY≤pTBFRand3-COLORABILITY≤pTBFC,\textsf{3-COLORABILITY} \le_p \textsf{TBF}_{\mathbb R} \quad\text{and}\quad \textsf{3-COLORABILITY} \le_p \textsf{TBF}_{\mathbb C},3-COLORABILITY≤p​TBFR​and3-COLORABILITY≤p​TBFC​,

polynomial-time many-one reductions between languages of binary codes. NP-hardness of tensor bilinear feasibility then follows from the NP-completeness of 3-colorability (Karp 1972), which is cited and not part of the statement.

Milestones, in the order the proof uses them

  1. Lemma 2.10. CGC_GCG​ has a nonzero complex zero iff GGG is 3-colorable.
  2. The minors (proof of Theorem 3.7). For nonzero v,w∈CN\mathbf v, \mathbf w \in \mathbb C^Nv,w∈CN, all 2×22\times22×2 minors of [v w][\mathbf v\ \mathbf w][v w] vanish iff v=cw\mathbf v = c\mathbf wv=cw for some c∈Cc \in \mathbb Cc∈C.
  3. Correctness over C\mathbb CC (proof of Theorem 3.7). For v≥1v \ge 1v≥1, TBFC(AG)  ⟺  G\mathrm{TBF}_{\mathbb C}(\mathcal A_G) \iff GTBFC​(AG​)⟺G is 3-colorable.
  4. Lemma 3.6. TBFR(B)  ⟺  TBFC(A)\mathrm{TBF}_{\mathbb R}(\mathcal B) \iff \mathrm{TBF}_{\mathbb C}(\mathcal A)TBFR​(B)⟺TBFC​(A) for real A\mathcal AA.
  5. Theorem 6.2 (further result). 3-colorability reduces to deciding whether 000 is an ℓ2\ell^2ℓ2-singular value, and whether it is an ℓ3\ell^3ℓ3-singular value, over R\mathbb RR and over C\mathbb CC.

Significance

Theorem 3.7 is the hub of the paper's bilinear results: Section 6 derives from it the NP-hardness of deciding whether 000 is a tensor singular value (Theorem 6.2), and the triple bilinear problem (Theorem 3.8) is handled by the same techniques. It separates tensors from matrices in a precise sense: the matrix analogue of (9) is decided by a determinant, while the tensor version encodes a classical NP-complete problem. It does not settle the hardness of the hyperdeterminant, since the dimensions l=v(2v+5)l = v(2v+5)l=v(2v+5), m=n=2v+1m = n = 2v+1m=n=2v+1 of AG\mathcal A_GAG​ violate the GKZ condition; the paper leaves that as a conjecture.

The result is proved in the paper; none of it is formalized on Prove2Me or, to our knowledge, elsewhere. A complete development produces a machine-checked polynomial-time reduction in Cook's Turing-machine model between two concrete algebraic problems, the color-encoding lemma for 3-colorability (reusable by any reduction that goes through polynomial systems, including Theorems 1.3 and 2.6 of the same paper), and the realification of bilinear systems.

Difficulty

The algebra of the reduction is short; two places need care. First, the bilinear system does not see CGC_GCG​ directly: (9) involves two independent vectors v,w\mathbf v, \mathbf wv,w, and the reduction recovers the quadratic system only because the minors force v=cw\mathbf v = c\mathbf wv=cw, which requires both vectors to be nonzero. Second, the converse direction must produce a nonzero u\mathbf uu satisfying the second and third families, a homogeneous linear system of 4v+24v+24v+2 equations in l=v(2v+5)l = v(2v+5)l=v(2v+5) unknowns; this works only when l>4v+2l > 4v+2l>4v+2, which fails for the empty graph, so the reduction has to treat v=0v = 0v=0 separately.

The goal also asks for polynomial-time computability in a Turing-machine model: decoding a graph code, rejecting malformed words, and writing the binary code of AG\mathcal A_GAG​ (or of its doubling over R\mathbb RR). This bookkeeping, not the algebra, is the bulk of a complete proof.

Formalization scope

  • Tensors are functions I → J → K → F over finite index types, 0-based: the paper's aijka_{ijk}aijk​ is A (i-1) (j-1) (k-1). Problem 3.1 requires u\mathbf uu, v\mathbf vv, w\mathbf ww each nonzero, not merely (u,v,w)≠0(\mathbf u, \mathbf v, \mathbf w) \neq 0(u,v,w)=0, and the three equation families come from one tensor.
  • The doubled tensor of Lemma 3.6 is indexed by I ⊕ I, J ⊕ J, K ⊕ K, with Sum.inr i standing for the paper's index l+il+il+i. Lemma 3.6 is stated for this explicit B\mathcal BB; a bare existential would be witnessed by any tensor with the right truth value. Only the equivalence of existence is stated, not the paper's "one-to-one" correspondence (which pairs u′=(p,q)\mathbf u' = (\mathbf p,\mathbf q)u′=(p,q) with the conjugate-type vector p−iq\mathbf p - i\mathbf qp−iq).
  • The paper does not print the matrices of the last 4v4v4v slices of AG\mathcal A_GAG​. The mission fixes the bilinearizations xiy^i−zz^x_i\hat y_i - z\hat zxi​y^​i​−zz^, yiz^−xix^iy_i\hat z - x_i\hat x_iyi​z^−xi​x^i​, xiz^−yiy^ix_i\hat z - y_i\hat y_ixi​z^−yi​y^​i​ and ∑j∼i(xix^i+xix^j+xjx^j)\sum_{j\sim i}(x_i\hat x_i + x_i\hat x_j + x_j\hat x_j)∑j∼i​(xi​x^i​+xi​x^j​+xj​x^j​); each restricts to the printed polynomial on w=v\mathbf w = \mathbf vw=v.
  • "Nontrivial zero" in the minors claim is read as both v\mathbf vv and w\mathbf ww nonzero; the correctness milestone over C\mathbb CC assumes v≥1v \geq 1v≥1, since for v=0v = 0v=0 the printed tensor has no slices and its system (9) is infeasible.
  • "Polynomially reducible" is formalized as Cook–Karp many-one reducibility (CookPvsNP.PolyReducible), which the paper's proof constructs and which implies the paper's oracle notion of §1.4. "NP-hard" in Theorem 6.2 is rendered the same way, and "(ℓ2\ell^2ℓ2 or ℓ3\ell^3ℓ3)" as the claim for each notion. Rational inputs are coded as numerator and denominator in binary; tensors as l,m,nl, m, nl,m,n and then their entries in lexicographic order.
  • The singular-value equations (20)–(21) are taken literally over C\mathbb CC, without conjugates and without normalization; only σ=0\sigma = 0σ=0 is used, where they are scale-free.
  • NP-completeness of 3-colorability is cited, not formalized; no statement mentions P, NP-membership, or hyperdeterminants.
  • A trivializing formalization is ruled out: the goal is the polynomial-time reduction, not the bare equivalence of yes-instances, and the reduction must handle every input word.

The mission reuses the published definitions CookPvsNP_defs (Turing machines, polynomial-time computability, many-one reducibility) and ProjSchedTW.Complexity.Encoding (the binary alphabet and number codes), and Mathlib's SimpleGraph.Colorable. Contributions welcome: the algebraic milestones (Lemma 2.10, the minors, Lemma 3.6) are independent and self-contained; a library of polynomial-time computable list and arithmetic operations on BSym codes would serve both this goal and Theorem 6.2.

Selected references

  • C. J. Hillar, L.-H. Lim, Most Tensor Problems Are NP-Hard, J. ACM 60(6), Art. 45, 2013. https://doi.org/10.1145/2512329 (arXiv:0911.1393v5, https://arxiv.org/abs/0911.1393)
  • R. M. Karp, Reducibility among combinatorial problems, in Complexity of Computer Computations, Plenum, 1972. https://doi.org/10.1007/978-1-4684-2001-2_9
  • L.-H. Lim, Singular values and eigenvalues of tensors: a variational approach, Proc. IEEE CAMSAP, 2005. https://arxiv.org/abs/math/0607648
  • I. M. Gelfand, M. M. Kapranov, A. V. Zelevinsky, Discriminants, Resultants, and Multidimensional Determinants, Birkhäuser, 1994. https://doi.org/10.1007/978-0-8176-4771-1
  • S. A. Cook, The P versus NP problem, Clay Mathematics Institute problem description, 2000. https://www.claymath.org/wp-content/uploads/2022/06/pvsnp.pdf
14 thms2 active usersReviewed
PreviousPage 1 of 9Next

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me