Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Graph Theory

79 missions · 47 completed

Missions

Open32Completed47All79
🏆Completed
CombinatoricsNumber Theory·Captain: xiangyazi24

Proofs from THE BOOKTextbook

Proofs from THE BOOK: verified results and open formalization tasks

This Textbook project develops a reusable Lean library around Martin Aigner and Günter M. Ziegler's Proofs from THE BOOK. It combines results imported from the existing proof_in_the_book repository with precise contribution targets from the sixth edition (2018). The aim is to preserve mathematical meaning, reuse existing proofs, and make the remaining work accessible to other contributors.

What is already verified

The original import contains 156 distinct platform-accepted results. Euclid, the original main theorem, is retained as a completed milestone when the project goal moves to the sixth-edition extension. Every one of the repository's 40 chapter topics has accepted results. Each result certifies its actual Lean statement, including its hypotheses; this does not certify every argument or every theorem in a chapter. Some proofs reuse Mathlib, while others were developed in the repository. Their source and proof notes retain that distinction.

The imported source snapshot is 873d52e0c88cd351f594221e70c3c5b3559777a9. Imported results use Lean 4.30.0 and Mathlib c5ea00351c28e24afc9f0f84379aa41082b1188f. Immutable public source links are used only where the linked source matches the verified artifact. Compatibility changes, unsuccessful attempts, and verification evidence are retained in the integration project.

The live goal is the explicit conjunction of the 21 linked sixth-edition extension targets. Its reduction connects these targets to the goal, so proving the remaining children advances the project. This goal is deliberately narrower than “every theorem and every proof in the book”; the unlinked topology tasks below are additional formalization work.

Sixth-edition contribution targets

New milestones explicitly marked 6th ed. cover Chapters 7 (spectral theorem and determinants), 15 (round circles and links), 35 (finite Kakeya), 37 (permanents and entropy), and 45 (probabilistic counting). They include the precise definitions and boundary conditions needed to state the results. Compiled Open targets are requests for proofs, not proved results. The spectral theorem has a direct Mathlib proof; community results are reused under their actual statements and with attribution.

Two Chapter 15 tasks intentionally remain unlinked mathematical milestones: the full non-equivalence assertion for the depicted Borromean, Tait, and trivial links, and the Fox-coloring invariance bridge for equivalent link diagrams. These invite formalization of the diagrams and the topology bridge as well as proof. The separate modular Fox calculations do not by themselves establish ambient non-equivalence.

The crossing-lemma target is the universal good-drawing form: actual injective edge arcs and exact finite intersection records appear in its interface. It does not assume the desired crossing bound. The Ramsey target preserves the real exponent for odd k. The related public Erdős–Ramsey result with a rounded exponent is identified as a supporting result, not as proof of that full target.

Chapter numbering and statement scope

Older milestones use the repository's chapter labels. Repository Chapters 1–21 match the bundled fourth edition; Chapter 22 inserts Van der Waerden's permanent theorem, and Chapters 23–40 correspond to fourth-edition Chapters 22–39. The sixth edition has 45 chapters, so these organizational labels are not sixth-edition chapter numbers. New milestones give sixth-edition numbers and printed source pages explicitly.

Some existing formalizations preserve narrower statements or additional premises. Examples include repository Chapter 13's dihedral-angle conclusion, Chapter 28's Dilworth lower-bound result, and the geometric premises in Chapter 36. Read the actual linked theorem and its description before reusing it. A chapter title or the former Euclid main theorem is not a completion certificate for the collection.

How to contribute

Choose an Open linked theorem and inspect its definitions, exact binders, Mathlib revision, and prior attempts. Reuse a compatible existing result when it proves that statement; preserve the original contributor's attribution. Submit a matching proof for verification. For an unlinked milestone, first formalize and review the source statement and its definitions. These are known textbook results awaiting formalization or proof in this project, rather than claims of new unresolved mathematics.

Source repository: https://github.com/xiangyazi24/proof_in_the_book

Book: Aigner and Ziegler, Proofs from THE BOOK, Sixth Edition (2018), https://doi.org/10.1007/978-3-662-57265-8

209 thms11 active users
🏆Completed
CombinatoricsComplexity TheoryOperations Research+1·Captain: mikedeng1

Scheduling Subject to Resource Constraints: Classification and Complexity II: Unit-Time Jobs on Two Uniform Machines with Unit Resources Are Strongly NP-hardResearch Paper

Motivation

Machine scheduling asks how to assign jobs to machines over time. In many applications a job also needs additional scarce resources while it runs: a tool, a skilled operator, a memory bank, a channel. Adding such resources can turn a problem with a polynomial algorithm into an NP-hard one. Błażewicz, Lenstra and Rinnooy Kan (DAM 1983) extended the standard three-field classification α ∣ β ∣ γ\alpha\,|\,\beta\,|\,\gammaα∣β∣γ of scheduling problems (Graham, Lawler, Lenstra and Rinnooy Kan 1979) with a resource field resλσρres\lambda\sigma\rhoresλσρ. They then drew the complete borderline between easy and hard problems for unit-time jobs, identical or uniform machines and the makespan criterion. Their Fig. 2 marks each problem type as polynomially solvable or NP-hard.

This mission formalizes the two hardness results of that classification that come from graph partition problems (Theorems 2 and 3, p. 15). Two identical machines are easy under any resource constraints (Theorem 1, after Garey and Johnson 1975). Theorems 2 and 3 show that a third identical machine, or two machines of different speeds, already makes the problem strongly NP-hard, once the number of unit resources is part of the input.

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​. Each machine processes at most one job at a time, and each job runs on one machine without interruption. Machine MiM_iMi​ has a speed qi>0q_i>0qi​>0, and every job has unit execution requirement pj=1p_j=1pj​=1, so it takes time 1/qi1/q_i1/qi​ on MiM_iMi​. Identical machines (PPP) are the case qi=1q_i=1qi​=1; uniform machines (QQQ) allow arbitrary speeds.

There are lll resources R1,…,RlR_1,\dots,R_lR1​,…,Rl​. Resource RhR_hRh​ has a positive integer size shs_hsh​, the amount available at any time. Job JjJ_jJj​ needs a nonnegative integer amount rhjr_{hj}rhj​ of RhR_hRh​ throughout its execution. A schedule assigns each job a machine μ(j)\mu(j)μ(j) and a start time Sj≥0S_j\ge 0Sj​≥0. Its completion time is Cj=Sj+1/qμ(j)C_j=S_j+1/q_{\mu(j)}Cj​=Sj​+1/qμ(j)​, and it is being executed at every time ttt with Sj≤t<CjS_j\le t<C_jSj​≤t<Cj​. A schedule is feasible when:

  • jobs on the same machine do not overlap in time;
  • at every time ttt, the set StS_tSt​ of jobs being executed satisfies
∑j∈Strhj≤sh(h=1,…,l).\sum_{j\in S_t} r_{hj}\le s_h\qquad(h=1,\dots,l).j∈St​∑​rhj​≤sh​(h=1,…,l).

The makespan is Cmax⁡=max⁡jCjC_{\max}=\max_j C_jCmax​=maxj​Cj​.

The resource type res⋅11res{\cdot}11res⋅11 means three things: the number lll of resources is part of the input, every size is sh=1s_h=1sh​=1, and every requirement satisfies rhj≤1r_{hj}\le1rhj​≤1. A unit resource is therefore a conflict: two jobs that both need it can never run at the same time. The problems here have no precedence constraints. Pm ∣ res⋅11, pj=1 ∣ Cmax⁡Pm\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Pm∣res⋅11,pj​=1∣Cmax​ and Qm ∣ res⋅11, pj=1 ∣ Cmax⁡Qm\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Qm∣res⋅11,pj​=1∣Cmax​ ask for a feasible schedule of minimum makespan. Their decision versions ask, for a threshold yyy, whether a feasible schedule with Cmax⁡≤yC_{\max}\le yCmax​≤y exists.

The source problems are two graph problems on a graph G=(V,E)G=(V,E)G=(V,E) with ∣V∣=3t|V|=3t∣V∣=3t:

  • PARTITION INTO TRIANGLES: can VVV be partitioned into ttt triples of pairwise adjacent vertices?
  • PARTITION INTO PATHS OF LENGTH 2: can VVV be partitioned into ttt triples, each with at most one nonadjacent pair, that is, each spanning a path of length 2?

Both are NP-complete (Garey and Johnson 1979, problems GT11 and GT13).

The construction of p. 15 introduces one job per vertex and one unit resource R{j,k}R_{\{j,k\}}R{j,k}​ per nonadjacent pair {j,k}\{j,k\}{j,k}, required by JjJ_jJj​ and JkJ_kJk​ only.

Formalization targets

Goal: Theorem 3

Q2 ∣ res⋅11, pj=1 ∣ Cmax⁡ is NP-hard in the strong sense.Q2\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}\ \text{is NP-hard in the strong sense.}Q2∣res⋅11,pj​=1∣Cmax​ is NP-hard in the strong sense.

Formally: if the language of PARTITION INTO PATHS OF LENGTH 2 is NP-hard, then the language of unary codes of yes-instances of the decision version of Q2 ∣ res⋅11, pj=1 ∣ Cmax⁡Q2\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Q2∣res⋅11,pj​=1∣Cmax​ is NP-hard. The two speeds are arbitrary positive integers.

Milestones

  1. The construction's key property (p. 15). In the constructed instance, two distinct jobs can be executed simultaneously if and only if their vertices are adjacent.
  2. The triangle equivalence (proof of Theorem 2). GGG has a partition into triangles if and only if the constructed instance on three identical machines has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t.
  3. Theorem 2. P3 ∣ res⋅11, pj=1 ∣ Cmax⁡P3\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}P3∣res⋅11,pj​=1∣Cmax​ is NP-hard in the strong sense, given the NP-hardness of PARTITION INTO TRIANGLES.
  4. The paths equivalence (proof of Theorem 3). GGG has a partition into paths of length 2 if and only if the constructed instance on two uniform machines with speeds q1=2q_1=2q1​=2, q2=1q_2=1q2​=1 has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t.

Significance

The results. Theorems 2 and 3 are two of the minimal NP-hard problems in the paper's classification. Together with Theorem 1 they place the borderline exactly: with unit resources whose number is part of the input, two identical machines are polynomial, while three identical machines, or two machines of different speeds, are strongly NP-hard. Strong NP-hardness rules out pseudo-polynomial algorithms unless P = NP, and it carries over to every more general resource type and machine environment in Fig. 1 and Fig. 2. Section 4.1 of the paper also derives hardness for ∑Cj\sum C_j∑Cj​ and Lmax⁡L_{\max}Lmax​ from these instances.

Formalizing it. The results are classical and proved on paper, but the paper's proofs are one sentence each ("Clearly", "It is easily seen"). No machine-checked proof exists, and the platform has no model of resource-constrained scheduling with real-valued time. This mission produces that model. It also produces a precise statement of strong NP-hardness on top of Cook's Turing-machine definitions, and the first formal NP-hardness reductions from graph partition problems to scheduling.

Difficulty

The scheduling half of each equivalence depends on the real-time model. On two uniform machines of speeds 2 and 1, jobs take time 12\tfrac1221​ and 111, so jobs on the fast machine start at half-integers or anywhere else. The resource constraint must hold at every real time, not at a finite set of checkpoints. An argument that treats time as integer slots applies to the triangle case but does not transfer to the paths case.

The complexity half needs polynomial-time computability of the construction on Cook's one-tape Turing machines, on encoded strings that include malformed inputs. It also needs closure of polynomial-time reductions under composition, which the imported complexity layer states but does not prove.

Formalization scope

  • Time and schedules. Start times are nonnegative reals, execution intervals are half-open [Sj,Cj)[S_j,C_j)[Sj​,Cj​), and the resource constraints are imposed at every real time. Schedules are nonpreemptive.
  • Indices. Jobs, machines and resources are 0-based (Fin n, Fin m, Fin l), so q1,q2q_1,q_2q1​,q2​ are q 0, q 1.
  • Decision versions. "NP-hard" refers to the decision version with a threshold yyy. Thresholds are natural numbers and the Q2Q2Q2 speeds are positive integers. This restricted problem is a subproblem of the one with rational data, so its hardness is the stronger statement.
  • Encodings and strong NP-hardness. Instances are strings over a two-letter alphabet with every number in unary. Graphs are ttt in unary followed by the 3t×3t3t\times 3t3t×3t adjacency matrix, so ∣V∣=3t|V|=3t∣V∣=3t is part of the instance. Languages contain only codes of yes-instances. Strong NP-hardness is NP-hardness of the unary code language. With unary numbers, Max(I)≤Length(I)\mathrm{Max}(I)\le\mathrm{Length}(I)Max(I)≤Length(I), so this is equivalent to Garey and Johnson's definition. The complexity layer is the published module CookPvsNP_defs.
  • Cited hypothesis. Each hardness theorem takes as its only hypothesis the NP-hardness of its source problem, which the paper cites from Garey and Johnson rather than proves. The hypothesis is a true statement about a nonempty, non-universal language. The statements are not weakened to a reduction between languages, and they assume nothing about P versus NP.
  • Source problems. The paper's phrase "three vertices, at most two of which are nonadjacent" is read as "at most one nonadjacent pair", which is Garey and Johnson's GT13. Reading it as "at most two nonadjacent pairs" would admit triples with a single edge and change the problem. PARTITION INTO PATHS OF LENGTH 2 reuses the published definition CubicP3Partition.P3Factor, a spanning non-induced P3P_3P3​-factor.
  • Construction. Resources are indexed by the nonadjacent pairs j<kj<kj<k in lexicographic order, one per unordered pair and none for a pair {j,j}\{j,j\}{j,j}. A diagonal resource would make every job infeasible.
  • Not trivial. A model that checks resources only at integer times, or only at start times, would make the paths equivalence false. A hypothesis on the target problem would make the goal circular. The definitions rule out both.

Welcome contributions: proofs of the two equivalences, polynomial-time computability of the construction on Cook's machines, and a general composition lemma for polynomial-time reductions. The composition lemma is reusable for every hardness mission built on CookPvsNP_defs.

Selected references

  • J. Błażewicz, J. K. Lenstra, A. H. G. Rinnooy Kan, Scheduling subject to resource constraints: classification and complexity, Discrete Applied Mathematics 5 (1983) 11–24. https://doi.org/10.1016/0166-218X(83)90012-4
  • M. R. Garey, D. S. Johnson, Complexity results for multiprocessor scheduling under resource constraints, SIAM Journal on Computing 4 (1975) 397–411. https://doi.org/10.1137/0204035
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, San Francisco, 1979, ISBN 0-7167-1045-5.
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Annals of Discrete Mathematics 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • S. Cook, The P versus NP problem, Clay Mathematics Institute Millennium Problems. https://www.claymath.org/wp-content/uploads/2022/06/pvsnp.pdf
41 thms7 active usersReviewed
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
🏆Completed
CombinatoricsTheoretical Computer Science·Captain: hao jia

Immune High-Girth Bipartite Graphs (Feghali-Lucke-Paulusma-Ries 2025)Research Paper

Motivation

A matching cut is a vertex bipartition whose crossing edges form a matching. The property was introduced under the name decomposability and has links to graph algorithms, stable cutsets in line graphs, and several graph-labeling problems. An Open Problem Garden question asked whether sufficiently large girth forces a matching cut once average degree is bounded.

Feghali, Lucke, Paulusma, and Ries answered that question negatively. Their conference paper appeared at ISAAC 2023, and the version of record was published in Algorithmica in 2025. The paper proves NP-completeness for bipartite graphs of arbitrarily prescribed girth and bounded maximum degree. A central input, Lemma 5, is a stronger structural existence statement: for every girth threshold there is an immune 141414-regular bipartite graph of at least that girth, and it has a perfect matching.

This is therefore a ResearchPaper mission, not a new open-problem mission. Its goal is to formalize the published theorem and its graph-theoretic consequence. Repository candidate constructions and finite arithmetic audits remain separate and are not credited as solving the problem.

Setting

For a finite simple graph GGG and a vertex set A⊆V(G)A\subseteq V(G)A⊆V(G), the associated cut consists of all edges with one endpoint in AAA and one in V(G)∖AV(G)\setminus AV(G)∖A. The cut is nontrivial when both shores are nonempty. It is a matching cut when each vertex is incident with at most one crossing edge. A graph is called immune in the cited paper when it has no matching cut.

The girth is the length of a shortest simple cycle; forests have infinite girth. A graph is 141414-regular when every vertex has exactly fourteen neighbors. Bipartiteness is witnessed by a partition into two independent sides. A perfect matching pairs every vertex with one adjacent partner.

The original OPG wording has a literal one-vertex boundary ambiguity: with nonempty shores required, K1K_1K1​ has no matching cut, average degree zero, and infinite girth. The research-paper target avoids that vacuity by constructing connected graphs with at least two vertices, exact degree fourteen, and arbitrarily large finite girth.

Formalization targets

Lemma 5 — immune high-girth graphs

The main theorem follows the paper's structural lemma:

∀g≥3 ∃G,G is finite, connected, bipartite, and 14-regular,\forall g\ge3\ \exists G, \quad G\text{ is finite, connected, bipartite, and $14$-regular}, ∀g≥3 ∃G,G is finite, connected, bipartite, and 14-regular, girth⁡(G)≥g,G has no matching cut,G has a perfect matching. \operatorname{girth}(G)\ge g, \qquad G\text{ has no matching cut}, \qquad G\text{ has a perfect matching}.girth(G)≥g,G has no matching cut,G has a perfect matching.

The graph may depend on ggg. The existence quantifier does not request an efficient algorithm or a numerical order bound.

Negative OPG consequence

A supporting theorem removes the perfect-matching and bipartite fields and records the direct substantive counterexample family:

∀g≥3 ∃G,d‾(G)=14<15,girth⁡(G)≥g,G has no matching cut.\forall g\ge3\ \exists G, \qquad \overline d(G)=14<15, \quad \operatorname{girth}(G)\ge g, \quad G\text{ has no matching cut}.∀g≥3 ∃G,d(G)=14<15,girth(G)≥g,G has no matching cut.

Thus choosing d=15d=15d=15 refutes the intended universal assertion that some girth threshold works for every graph of average degree below ddd.

Significance

The theorem shows that large girth and bounded degree do not force matching cuts. The examples are highly nontrivial: they are connected, regular, bipartite, and can have arbitrarily large girth. This separates local tree-like structure from the global expansion that prevents a matching cut.

Within the paper, the immune graphs serve as gadgets for hardness reductions. The journal theorem states that, for every g≥3g\ge3g≥3, Matching Cut is NP-complete even for bipartite graphs of girth at least ggg and maximum degree at most 606060. Formalizing Lemma 5 supplies the graph-theoretic core needed to reconstruct that result without forcing this mission to formalize an entire complexity-theory reduction in its first stage.

Difficulty

Large girth alone makes bounded neighborhoods look like trees, and trees have many matching cuts. Immunity must therefore come from global expansion rather than short local cycles. The paper obtains the required family from Lubotzky–Phillips–Sarnak Ramanujan graphs and combines spectral and isoperimetric bounds to show that every nontrivial cut has too many crossing incidences to be a matching.

A formal proof must bridge several exact interfaces: existence of suitable primes, the finite Cayley-graph construction, bipartiteness and regularity, the girth lower bound, the spectral-to-isoperimetric inequality, and Hall's theorem for the perfect matching. None of these can be replaced by a finite sample or an asymptotic slogan.

Formalization scope

Graphs are finite and simple. Connectedness is nonempty mutual graph reachability. A simple cycle is a cyclic list of at least three distinct vertices; girth at least ggg means every such cycle has length at least ggg, so forests satisfy every threshold. A matching cut requires both shores nonempty and is encoded by the condition that every vertex has at most one crossing neighbor. A perfect matching is represented by an adjacent involution.

The main theorem explicitly requires at least two vertices, although exact 14-regularity already forces nontrivial order; the redundant bound documents exclusion of the K1K_1K1​ ambiguity. The mission does not claim that the frozen OPG contract was well-posed at order one. It formalizes the paper's substantive counterexample family and the consequence for the intended question.

Candidate files in the associated repository explore alternative bounded-degree constructions and integer counts. They are candidate_only and are not proof dependencies. Contributions should follow the published Lemma 5 and its cited inputs, or provide a separately sourced proof of the same declaration. The later maximum-degree-60 NP-completeness theorem is welcome as a future extension after the finite complexity framework is fixed.

Selected references

  • C. Feghali, F. Lucke, D. Paulusma, and B. Ries, Matching Cuts in Graphs of High Girth and H-Free Graphs, Algorithmica 87 (2025), 1199–1221. https://doi.org/10.1007/s00453-025-01318-8
  • C. Feghali, F. Lucke, D. Paulusma, and B. Ries, ISAAC 2023 version. https://drops.dagstuhl.de/entities/document/10.4230/LIPIcs.ISAAC.2023.31
  • A. Lubotzky, R. Phillips, and P. Sarnak, Ramanujan graphs, Combinatorica 8 (1988), 261–277. https://doi.org/10.1007/BF02126799
  • Open Problem Garden, Matching cut and girth. https://www.openproblemgarden.org/op/matching_cut_and_girth
20 thms6 active usersReviewed
🏆Completed
Combinatorics·Captain: hao jia

Formalizing an 8-Vertex Candidate Counterexample to the Geodesic-Cycle Assignment Problem (OPG-500)Open Problem

[VM-STATUS-20260908-R05-PROVED]

Status update (2026-09-08): The root theorem OPG500Counterexample.eight_vertex_counterexample is now Proved by an accepted Prove2Me submission. All six milestones are proved and the root has zero open leaves. The accepted proof has also been independently rebuilt with Lean 4.33.1 / Mathlib 0df444a360eaa60ab8c11dca51a86af692955474. The historical text below describes the mission as it stood before formal closure.


Motivation and historical context

Peripheral cycles occupy a distinguished place in structural graph theory. A cycle is peripheral when it is induced and does not separate the graph after its vertices are removed. Tutte proved in 1963 that the peripheral cycles of a finite 3-connected graph generate its binary cycle space. This theorem links a local, visibly embedded kind of cycle to the global algebraic structure of all cycles.

Weighted geodesic cycles provide a different generating family. Georgakopoulos and Sprüssel proved in 2009 that, for every finite graph with positive edge lengths, every cycle is a binary sum of weighted geodesic cycles whose lengths do not exceed the length of the original cycle. In the same paper they posed Problem 3: can the edges of every finite 3-connected graph be assigned positive lengths so that every weighted geodesic cycle is peripheral? A positive answer would recover Tutte's generation theorem through metric structure.

The present target tests the opposite possibility on one explicitly specified graph with eight vertices. A candidate argument and finite certificates are available in the frozen OPG-500 research repository, but those artifacts are explicitly marked candidate_only: they are neither a published counterexample nor a machine-checked resolution. The purpose of the formal target is to determine whether the proposed universal obstruction survives complete definition, proof, and statement-faithfulness checks.

Setting

Let GGG be a finite simple graph. A positive edge-length assignment is a function

ℓ:E(G)⟶R\ell:E(G)\longrightarrow \mathbb Rℓ:E(G)⟶R

such that ℓ(e)>0\ell(e)>0ℓ(e)>0 for every edge eee. The length of a finite path or cycle is the sum of the lengths of its edges.

A simple cycle CCC is ℓ\ellℓ-geodesic when, for every pair of vertices x,yx,yx,y on CCC, at least one of the two xxx–yyy arcs of CCC has length equal to the shortest-path distance between xxx and yyy in GGG. Equivalently, there is no xxx–yyy path in GGG whose length is strictly smaller than both xxx–yyy arcs of CCC. The definition concerns vertices of the cycle and permits ties between shortest paths.

A simple cycle is peripheral when it is induced and deleting all of its vertices leaves a connected graph or the empty graph. This is vertex deletion, not edge deletion.

Fix the graph HHH on vertices 0,1,…,70,1,\ldots,70,1,…,7. The vertices 0,1,2,30,1,2,30,1,2,3 induce K4K_4K4​. For each i∈{0,1,2,3}i\in\{0,1,2,3\}i∈{0,1,2,3}, set yi=7−iy_i=7-iyi​=7−i and join yiy_iyi​ to exactly the three core vertices other than iii. The four vertices yiy_iyi​ are pairwise nonadjacent. Thus the frozen edge set is

{01,02,03,04,05,06,12,13,14,15,17,23,24,26,27,35,36,37}.\{01,02,03,04,05,06,12,13,14,15,17,23,24,26,27,35,36,37\}.{01,02,03,04,05,06,12,13,14,15,17,23,24,26,27,35,36,37}.

The labels and edge set are part of the statement and are not interchangeable with earlier candidate labelings without an explicit isomorphism.

Formalization targets

Main target: the universal eight-vertex obstruction

Formalize the following statement for the fixed graph HHH:

H is 3-connectedand∀ℓ:E(H)→R>0,  ∃C,  C is an ℓ-geodesic simple cycle of H and is not peripheral.H\text{ is 3-connected}\quad\text{and}\quad \forall\ell:E(H)\to\mathbb R_{>0},\; \exists C,\; C\text{ is an $\ell$-geodesic simple cycle of $H$ and is not peripheral}.H is 3-connectedand∀ℓ:E(H)→R>0​,∃C,C is an ℓ-geodesic simple cycle of H and is not peripheral.

The existential cycle may depend on ℓ\ellℓ. The universal quantifier includes all strictly positive real assignments, including assignments with tied shortest paths. This is the stable target; finite samples and rational specializations are subordinate checks rather than replacements for it.

Supporting targets

The development should also formalize the finite weighted geodesic-cycle generation theorem of Georgakopoulos and Sprüssel, the exact 3-connectivity and peripheral-cycle classification of HHH, the required shortest-path and tight-subgraph statements, the finite cycle-space rank statements, and the finite minimum/descent principle used to select a cycle outside a closed binary span. These targets should remain separate declarations so that their assumptions and reuse boundaries are visible.

Significance

A proof of the main target would give a negative answer to the finite problem by exhibiting a 3-connected graph for which no positive edge weighting can make all geodesic cycles peripheral. It would not contradict Tutte's theorem: peripheral cycles may still generate the cycle space even though they cannot be made to contain every geodesic cycle for any weighting. The distinction between these two generation mechanisms is part of the mathematical content.

A formal development would add more than a checked final sentence. It would provide reusable definitions for positively weighted finite graphs and vertex-geodesic cycles, a precise treatment of the two arcs between cycle vertices, explicit deletion semantics for peripheral cycles, and finite cycle-space infrastructure. It would also separate purely finite graph facts from statements quantified over arbitrary real weights. The candidate repository currently supplies finite enumeration and abstract Lean fragments, but no existing artifact checks this full dependency chain.

Until the complete main theorem is verified, the eight-vertex graph remains a candidate obstruction and the original problem remains unresolved by this development.

Difficulty

The central difficulty is the universal quantification over real edge lengths. Testing many integer or rational vectors cannot cover it. Shortest paths need not be unique, so an argument that silently perturbs the weights or assumes unique geodesics can change which cycles are geodesic. Every strict and weak inequality must therefore agree with the source definition, including tie cases.

The graph is small but the semantic boundary is not. A formal cycle representation must expose the two cycle arcs for every vertex pair without admitting malformed or repeated-vertex objects. The peripheral predicate must combine inducedness with connectivity after vertex deletion and must classify all cycles, not only a selected family of triangles. Finally, finite cycle-space computations and rank inequalities must be connected to actual paths and weighted geodesicity; a propositional or enumerative certificate alone does not establish that bridge.

Formalization scope

The Lean development will use Fin 8 for the vertices of HHH and a SimpleGraph representation for adjacency. Weights will be functions on the edge subtype, so values on nonedges cannot affect the theorem. All weights are real and strictly positive. Paths and cycles are finite and simple; arbitrary walks do not count as target witnesses. Geodesicity is vertex-based and includes tied shortest paths. Peripheral cycles use inducedness and vertex deletion, with a connected-or-empty remainder.

The main theorem must retain the quantifier order “for every weighting, there exists a cycle.” It may not be weakened to rational weights, finitely many tested assignments, nonnegative weights, one selected weighting, edge-geodesicity, or the assertion that only four named core triangles fail to be peripheral. Definitions must be sorry-free, and nontrivial mathematical claims must be theorem declarations with separately checked proofs.

Reusable contributions include finite weighted-path length, shortest-path attainment in finite positive graphs, the equivalence of the two geodesic formulations, cycle-arc APIs, vertex-deletion connectivity, binary edge-vector encodings, and finite descent outside a closed span. Graph-specific finite certificates are welcome only when their checker is represented in Lean or their conclusions are otherwise proved in the kernel.

Selected references

  • A. Georgakopoulos and P. Sprüssel, Geodetic topological cycles in locally finite graphs, Electronic Journal of Combinatorics 16 (2009), R144. Section 3.1, Theorem 3.1; Section 5, Problem 3. https://arxiv.org/abs/0911.3999v1
  • Open Problem Garden, Geodesic cycles and Tutte's Theorem, problem statement and vertex-based definition. https://www.openproblemgarden.org/op/geodesic_cycles_and_tuttes_theorem
  • W. T. Tutte, How to draw a graph, Proceedings of the London Mathematical Society 13 (1963), 743–768. Cited as reference [18] by Georgakopoulos and Sprüssel for peripheral-cycle generation.
  • Vibe Mathing, frozen OPG-500 candidate repository at commit a41fe59b4535851ea55f6e868e938b9aaf81e924. https://github.com/vibemathing/problem-opg-500-geodesic-cycles/tree/a41fe59b4535851ea55f6e868e938b9aaf81e924
9 thms6 active usersReviewed
🏆Completed
Combinatorics·Captain: mikedeng1

Applied Combinatorics I: Graph Theory and Dirac's Hamiltonicity TheoremTextbook

Motivation

Graphs are the most basic combinatorial model of pairwise relations: road networks, frequency interference between radio stations, schedules, circuit layouts. Chapter 5 of Keller and Trotter's Applied Combinatorics (appliedcombinatorics.org, CC BY-SA 4.0) introduces the vocabulary of graph theory and proves its first structural theorems: when a graph can be traversed edge by edge (Euler, 1736), when it can be toured vertex by vertex (Dirac, 1952), when it can be colored with two colors, and how far the chromatic number can drift from the size of the largest clique.

Hamiltonicity is the vertex analogue of Euler's edge-traversal problem and behaves very differently: Euler's problem has a simple parity characterization, while deciding whether a graph has a hamiltonian cycle is NP-complete (Karp, 1972). Sufficient conditions are therefore the main tool, and the minimum-degree condition of Dirac (1952) is the first and most cited of them; Ore's condition (1960) and the Bondy–Chvátal closure (1976) refine it.

Setting

A graph G=(V,E)G = (V, E)G=(V,E) consists of a finite vertex set VVV and a set EEE of 2-element subsets of VVV, the edges; xy∈Exy \in Exy∈E means xxx and yyy are adjacent. The degree deg⁡G(v)\deg_G(v)degG​(v) is the number of neighbours of vvv. The mission uses Mathlib's SimpleGraph V with [Fintype V]; "a graph on nnn vertices" means Fintype.card V = n.

  • A cycle is a sequence (x1,…,xn)(x_1, \dots, x_n)(x1​,…,xn​) of n≥3n \ge 3n≥3 distinct vertices with xixi+1∈Ex_i x_{i+1} \in Exi​xi+1​∈E for i<ni < ni<n and x1xn∈Ex_1 x_n \in Ex1​xn​∈E; its length is nnn. A graph is acyclic if it has no cycle, and a tree if it is connected and acyclic. A leaf of a tree is a vertex of degree 111.
  • A hamiltonian cycle is a sequence (x1,…,xn)(x_1, \dots, x_n)(x1​,…,xn​) in which every vertex appears exactly once, xixi+1∈Ex_i x_{i+1} \in Exi​xi+1​∈E for i<ni < ni<n, and x1xn∈Ex_1 x_n \in Ex1​xn​∈E. A graph is hamiltonian if it has one (AppliedComb.Graphs.IsHamiltonian).
  • An eulerian circuit is a sequence (x0,…,xt)(x_0, \dots, x_t)(x0​,…,xt​), repetition allowed, with x0=xtx_0 = x_tx0​=xt​, consecutive entries adjacent, and every edge equal to xixi+1x_i x_{i+1}xi​xi+1​ for exactly one i<ti < ti<t. A graph without isolated vertices is eulerian if it has one (IsEulerian).
  • A proper coloring assigns colors to vertices so that adjacent vertices differ; the chromatic number χ(G)\chi(G)χ(G) is the least number of colors in a proper coloring, and the clique number ω(G)\omega(G)ω(G) is the largest size of a set of pairwise adjacent vertices.
  • An interval graph is the intersection graph of closed intervals [av,bv]⊂R[a_v, b_v] \subset \mathbb R[av​,bv​]⊂R, v∈Vv \in Vv∈V: distinct u,vu, vu,v are adjacent iff their intervals meet (IsIntervalGraph).

Formalization targets

Goal: Dirac's theorem (Theorem 5.18)

If ∣V∣=n≥1 and deg⁡G(v)≥⌈n2⌉ for all v∈V, then G is hamiltonian.\text{If } |V| = n \ge 1 \text{ and } \deg_G(v) \ge \left\lceil \tfrac n2 \right\rceil \text{ for all } v \in V, \text{ then } G \text{ is hamiltonian.}If ∣V∣=n≥1 and degG​(v)≥⌈2n​⌉ for all v∈V, then G is hamiltonian.

Milestones

In the book's order:

  1. Proposition 5.11. A tree on n≥2n \ge 2n≥2 vertices has at least two leaves.
  2. Theorem 5.13 (Euler). A graph without isolated vertices is eulerian if and only if it is connected and every degree is even.
  3. Theorem 5.21. χ(G)≤2\chi(G) \le 2χ(G)≤2 if and only if GGG contains no odd cycle.
  4. Proposition 5.24 (generalized pigeonhole). If f:X→Yf : X \to Yf:X→Y and ∣X∣≥(m−1)∣Y∣+1|X| \ge (m-1)|Y| + 1∣X∣≥(m−1)∣Y∣+1, some fibre of fff contains mmm distinct elements.
  5. Proposition 5.25. For every t≥3t \ge 3t≥3 there is a finite graph GtG_tGt​ with χ(Gt)=t\chi(G_t) = tχ(Gt​)=t and ω(Gt)=2\omega(G_t) = 2ω(Gt​)=2.
  6. Theorem 5.28. Every finite interval graph satisfies χ(G)=ω(G)\chi(G) = \omega(G)χ(G)=ω(G).

The book's proof of Theorem 5.18 uses only the pigeonhole principle, so none of these results lies on its path. The milestones are the chapter's other theorems about the same objects (cycles, degrees, colorings), and they share the goal's definitions. Cayley's formula (Theorem 5.39, nn−2n^{n-2}nn−2 labelled trees on nnn vertices) is already on the platform as GYGraphTheory.cayley_tree_formula and is included as a reference item.

Significance

Dirac's theorem is sharp: the complete bipartite graph Kk,k+1K_{k,k+1}Kk,k+1​ has minimum degree k=⌈n/2⌉−1k = \lceil n/2 \rceil - 1k=⌈n/2⌉−1 and no hamiltonian cycle. It is the starting point for the theory of degree conditions for hamiltonicity (Ore, Pósa, Chvátal) and for its extremal and random-graph analogues. Euler's theorem gives linear-time recognition of eulerian graphs; Theorem 5.21 characterizes bipartite graphs; Proposition 5.25 shows that local sparsity (no triangles) does not bound the chromatic number; Theorem 5.28 is the first step towards perfect graphs.

All of these results have been proved for a long time. What is missing is their formalization. Mathlib provides SimpleGraph, chromaticNumber, cliqueNum, degree-sum identities, and definitions of eulerian and hamiltonian walks. It contains no Dirac theorem, no converse of the eulerian degree condition, and no interval-graph theory. On the platform, FamousTheorems.two_colorable_iff_no_odd_cycle_6b is stated with odd closed walks rather than odd cycles, and triangle_free_chromatic_number asserts only χ≥k\chi \ge kχ≥k rather than χ=t\chi = tχ=t with ω=2\omega = 2ω=2 exactly.

Difficulty

For Dirac's theorem the obvious approach, induction on nnn (deleting a vertex), fails: deleting a vertex lowers degrees, and the hypothesis deg⁡≥⌈n/2⌉\deg \ge \lceil n/2 \rceildeg≥⌈n/2⌉ is not inherited by the smaller graph. The degree condition is global, and so is the argument that uses it. In Lean the argument also has to reverse and splice vertex sequences while keeping distinctness and every adjacency along the way. The sufficiency half of Euler's theorem has to construct a circuit that uses every edge exactly once, not just show that one exists up to parity. In Theorem 5.21 the obstruction is a cycle with distinct vertices; producing one from a failed 2-coloring takes more work than producing an odd closed walk. Proposition 5.25 needs a concrete infinite family of graphs and exact computation of both invariants. The upper bound χ≤t\chi \le tχ≤t is easy, but the lower bound χ≥t\chi \ge tχ≥t is not.

Formalization scope

  • Graphs are SimpleGraph V on a Fintype V. The number of vertices is always Fintype.card V, never a free parameter; in Theorem 5.18 the ceiling ⌈n/2⌉\lceil n/2 \rceil⌈n/2⌉ is (n + 1) / 2 in N\mathbb NN, with Fintype.card V = n.
  • "Hamiltonian" is the book's definition (p. 79), a list of all vertices without repetition with consecutive and first/last entries adjacent. It is not Mathlib's Walk.IsHamiltonianCycle, which needs at least three vertices: under that notion Theorem 5.18 would be false for K2K_2K2​. The goal assumes VVV nonempty, as the book's proof does (nnn positive). Degenerate readings are ruled out: a predicate satisfied by a cycle through only some vertices, or a statement in which nnn is not the number of vertices, would make the goal vacuous or false, and neither is used here.
  • "Eulerian" follows p. 75. Since the book defines it only for graphs without isolated vertices, Theorem 5.13 carries that hypothesis explicitly.
  • Cycles are the book's cycles (distinct vertices, length ≥3\ge 3≥3), not closed walks. Theorem 5.21 is stated with odd cycles.
  • χ\chiχ is Mathlib's chromaticNumber (N∞\mathbb N_\inftyN∞​-valued) and ω\omegaω is cliqueNum; on finite graphs both agree with the book's definitions (pp. 81, 84).
  • Planarity (Section 5.5: Euler's formula 5.32, the bound 3n−63n - 63n−6 in 5.33, Kuratowski 5.34, the Four Color Theorem 5.37) is excluded. The book defines planar drawings through polygonal arcs in R2\mathbb R^2R2 and faces, which is a topological development rather than a chapter mission. Kuratowski's theorem and the Four Color Theorem are not proved in the book.
  • Contributions reusable beyond this mission are welcome: path/cycle manipulation lemmas for list-based sequences, the Euler circuit construction, bipartiteness via distance parity, and the Mycielski construction.

Selected references

  • M. T. Keller and W. T. Trotter, Applied Combinatorics, 2017 Edition, Chapter 5. https://www.appliedcombinatorics.org/
  • G. A. Dirac, Some theorems on abstract graphs, Proc. London Math. Soc. (3) 2 (1952), 69–81. https://doi.org/10.1112/plms/s3-2.1.69
  • L. Euler, Solutio problematis ad geometriam situs pertinentis, Commentarii Academiae Scientiarum Petropolitanae 8 (1741), 128–140 (read 1736). https://scholarlycommons.pacific.edu/euler-works/53/
  • J. Mycielski, Sur le coloriage des graphs, Colloquium Mathematicum 3 (1955), 161–162. https://doi.org/10.4064/cm-3-2-161-162
  • R. M. Karp, Reducibility among combinatorial problems, in Complexity of Computer Computations (1972), 85–103. https://doi.org/10.1007/978-1-4684-2001-2_9
12 thms5 active usersReviewed
🏆Completed
CombinatoricsLinear OptimizationOperations Research+1·Captain: mikedeng1

Theoretical Improvements in Algorithmic Efficiency for Network Flow Problems 3: Capacity Scaling for the Hitchcock ProblemResearch Paper

Motivation

The Hitchcock transportation problem asks how to ship a commodity from mmm supply points to nnn demand points at minimum total cost. It was posed by Hitchcock in 1941 and is one of the founding problems of linear programming and network optimization; it is solved routinely in logistics, and its structure (a bipartite network with supplies, demands and per-unit costs) recurs in assignment, optimal transport and matching.

The classical algorithms for it, the Ford–Fulkerson primal–dual method among them, augment flow one path at a time. With integral data their number of augmentations is bounded only by the total supply ∑iai\sum_i a_i∑i​ai​, which is exponential in the number of binary digits used to write the data. Edmonds and Karp, Theoretical Improvements in Algorithmic Efficiency for Network Flow Problems (J. ACM 19(2), 1972, doi:10.1145/321694.321699), introduced capacity scaling: solve a coarse version of the problem first, then refine one binary digit at a time. Their Theorem 9 (p. 260) bounds the total number of augmentations by a quantity proportional to max⁡(m,n)\max(m,n)max(m,n) times the number of bits of the data, which made the transportation problem, and through standard reductions the minimum-cost flow problem, one of the first network problems with a polynomial-time ("good") algorithm in the sense of Edmonds.

Timeline. Hitchcock (1941) posed the problem; Ford and Fulkerson (1956–1962) gave the primal–dual labeling method and the optimality conditions by node potentials; Edmonds and Karp (1972) gave the scaling method and the bound formalized here. Strongly polynomial algorithms, independent of the size of the numbers, came later (Tardos 1985; Orlin 1988).

Setting

The network of Figure 1 (p. 259) has a source sss, a sink ttt, supply nodes s1,…,sms_1,\dots,s_ms1​,…,sm​ and demand nodes t1,…,tnt_1,\dots,t_nt1​,…,tn​, with m,n≥1m,n\ge 1m,n≥1. Its arcs are (s,si)(s,s_i)(s,si​) with capacity aia_iai​ and cost 000; (si,tj)(s_i,t_j)(si​,tj​) with capacity +∞+\infty+∞ and cost dij≥0d_{ij}\ge 0dij​≥0; (tj,t)(t_j,t)(tj​,t) with capacity bjb_jbj​ and cost 000; and the return arc (t,s)(t,s)(t,s) with capacity +∞+\infty+∞ and cost 000. The supplies aia_iai​ and demands bjb_jbj​ are positive integers with ∑iai=∑jbj=:B\sum_i a_i=\sum_j b_j=:B∑i​ai​=∑j​bj​=:B.

A flow assigns a nonnegative number to every arc, at most the capacity, with inflow equal to outflow at every node. Write f0i=f(s,si)f_{0i}=f(s,s_i)f0i​=f(s,si​), fij=f(si,tj)f_{ij}=f(s_i,t_j)fij​=f(si​,tj​), fj0=f(tj,t)f_{j0}=f(t_j,t)fj0​=f(tj​,t); the value of fff is f(t,s)f(t,s)f(t,s), and a maximum flow is one of largest value. Its cost is ∑i,jdijfij\sum_{i,j} d_{ij} f_{ij}∑i,j​dij​fij​; a flow is extreme if no flow of the same value is cheaper. A flow is pseudo-extreme if there are real ui,vju_i, v_jui​,vj​ with ui−vj+dij≥0u_i-v_j+d_{ij}\ge 0ui​−vj​+dij​≥0 for all i,ji,ji,j and fij=0f_{ij}=0fij​=0 whenever ui−vj+dij>0u_i-v_j+d_{ij}>0ui​−vj​+dij​>0.

An augmenting path relative to fff is a sequence of distinct nodes from sss to ttt in which each step either follows an arc with spare capacity or traverses backwards an arc carrying positive flow; augmenting pushes the minimum spare amount ε\varepsilonε along it and raises f(t,s)f(t,s)f(t,s) by ε\varepsilonε.

For p≥0p\ge 0p≥0, Problem ppp has the same network and costs, with capacities ⌊ai/2p⌋\lfloor a_i/2^p\rfloor⌊ai​/2p⌋ and ⌊bj/2p⌋\lfloor b_j/2^p\rfloor⌊bj​/2p⌋. Choose lll with every ai,bj<2la_i,b_j<2^lai​,bj​<2l. The scaling method solves Problems l−1,l−2,…,0l-1,l-2,\dots,0l−1,l−2,…,0 in turn. Each phase performs augmentations keeping every flow pseudo-extreme, until no augmenting path is left. Problem l−1l-1l−1 starts from the zero flow, and Problem p−1p-1p−1 starts from twice the final flow of Problem ppp.

Formalization targets

Goal: Theorem 9

For every run of the scaling method, the total number ∑p<lKp\sum_{p<l} K_p∑p<l​Kp​ of flow augmentations satisfies

∑p=0l−1Kp  ≤  max⁡(m,n)(2+⌊log⁡2∑i=1maimax⁡(m,n)⌋).\sum_{p=0}^{l-1} K_p \;\le\; \max(m,n)\left(2+\left\lfloor \log_2\frac{\sum_{i=1}^m a_i}{\max(m,n)}\right\rfloor\right).p=0∑l−1​Kp​≤max(m,n)(2+⌊log2​max(m,n)∑i=1m​ai​​⌋).

The bound holds for every lll admissible for the data and every choice of costs, and is stated with the paper's constant exactly.

Milestones

  1. §1.1: augmentation preserves feasibility and raises the value by ε>0\varepsilon>0ε>0; a flow is maximum iff no augmenting path exists.
  2. Theorem 8: a maximum flow is extreme iff there are potentials u0,…,umu_0,\dots,u_mu0​,…,um​, v0,…,vnv_0,\dots,v_nv0​,…,vn​ with (5a)–(5f).
  3. §2.2: a pseudo-extreme maximum flow is extreme.
  4. Lemma 3: if fff is pseudo-extreme in Problem ppp, then 2f2f2f is pseudo-extreme in Problem p−1p-1p−1.
  5. The maximum-flow value of Problem ppp is fp∗=min⁡(∑i⌊ai/2p⌋,∑j⌊bj/2p⌋)f_p^*=\min\big(\sum_i\lfloor a_i/2^p\rfloor,\sum_j\lfloor b_j/2^p\rfloor\big)fp∗​=min(∑i​⌊ai​/2p⌋,∑j​⌊bj​/2p⌋).
  6. Eq. (6): ∑pKp≤f0∗−∑p=1l−1fp∗\sum_p K_p\le f_0^*-\sum_{p=1}^{l-1} f_p^*∑p​Kp​≤f0∗​−∑p=1l−1​fp∗​.
  7. fp∗≥max⁡(0, B/2p−max⁡(m,n))f_p^*\ge\max\big(0,\,B/2^p-\max(m,n)\big)fp∗​≥max(0,B/2p−max(m,n)).

Significance

The result. Theorem 9 shows that scaling reduces the number of augmentations from order BBB to order max⁡(m,n)log⁡2(B/max⁡(m,n))\max(m,n)\log_2(B/\max(m,n))max(m,n)log2​(B/max(m,n)), which is roughly the length of the binary encoding of the data. Combined with the O(mn)O(mn)O(mn) cost of one augmentation it gives a polynomial-time algorithm for the transportation problem; through the reduction of minimum-cost flow to transportation (p. 261) it gives one for minimum-cost flow. Capacity and cost scaling became standard techniques in network optimization and in combinatorial optimization generally. Theorem 8 and the pseudo-extreme criterion are the optimality certificates for transportation, a special case of linear-programming complementary slackness.

Formalizing it. The theorem has been proved since 1972. This mission produces a machine-checked version of the bound, of the exact counting argument (eq. (6)) and of the arithmetic estimate that turns it into the stated constant, together with the potential-based optimality conditions for the transportation network. To our knowledge no machine-checked bound on the number of augmentations of a flow algorithm exists on the platform.

Difficulty

The obvious argument bounds the number of augmentations by the increase in flow value, since each augmentation raises the value by a positive integer. Applied to Problem 0 directly this gives only BBB. The scaling bound needs three things that the naive count does not supply. First, doubling the final flow of Problem ppp must give a feasible, still pseudo-extreme, starting flow for Problem p−1p-1p−1 (Lemma 3). Second, the gap between that start and the optimum of Problem p−1p-1p−1 must be small, which requires the exact maximum-flow value fp∗f_p^*fp∗​ of each scaled problem. Third, the telescoping sum of the gaps must be estimated against log⁡2(B/max⁡(m,n))\log_2(B/\max(m,n))log2​(B/max(m,n)) with the floors handled exactly. Integrality of every intermediate flow is not assumed; it has to be carried along the run from the integral capacities and the zero start.

Formalization scope

Everything lives in the namespace EdmondsKarp.Scaling. Nodes form an inductive type (s, t, src i, dst j). A flow is a structure with components f0, fx, fz, ret for the four arc families, real valued; the infinite capacities are encoded by the absence of an upper bound. IsMaxFlow is a predicate comparing values with every flow, not a supremum. Problem ppp uses natural-number division for ⌊ai/2p⌋\lfloor a_i/2^p\rfloor⌊ai​/2p⌋. Augmenting paths are lists of distinct nodes from s to t whose consecutive pairs have positive residual amount (resCap, valued in WithTop ℝ); they never use the return arc.

A run of the scaling method (IsScalingRun) is a family of phases F p 0,…,F p (K p)F\,p\,0,\dots,F\,p\,(K\,p)Fp0,…,Fp(Kp) for p<lp<lp<l. It starts from 000 in Problem l−1l-1l−1, restarts from 2F p (K p)2F\,p\,(K\,p)2Fp(Kp) in Problem p−1p-1p−1, and advances by one augmentation per step. Every flow is pseudo-extreme, and each phase ends with no augmenting path left. The paper's path-selection rule (minimum weight for the modified reduced costs Δˉ\bar\DeltaΔˉ) is abstracted to this invariant, which the paper states for it, so every run of the paper's method is covered.

The goal compares the count, cast to Z\mathbb{Z}Z, with max⁡(m,n) (2+⌊log⁡2(B/max⁡(m,n))⌋)\max(m,n)\,(2+\lfloor\log_2(B/\max(m,n))\rfloor)max(m,n)(2+⌊log2​(B/max(m,n))⌋), using Real.logb 2 and Int.floor. The constant is the printed one; neither O(⋅)O(\cdot)O(⋅) nor a weaker constant is acceptable. Positivity of all ai,bja_i,b_jai​,bj​ is a hypothesis, because without it the printed bound is false (for m=5m=5m=5, n=1n=1n=1, a=(1,0,0,0,0)a=(1,0,0,0,0)a=(1,0,0,0,0), b=(1)b=(1)b=(1) one augmentation is needed and the bound is negative). A formalization in which runs could be empty or never reach a maximum flow would trivialize the goal; the run predicate forbids this, and it is satisfiable (for instance with m=n=1m=n=1m=n=1, a=b=(1)a=b=(1)a=b=(1), l=1l=1l=1 and one augmentation).

A complete development needs max-flow/min-cut for the bipartite network, integrality of flows along a run, the exact value of fp∗f_p^*fp∗​, the telescoping identity and floor/logarithm estimates. LP duality or complementary slackness is needed only for Theorem 8. The augmenting-path and max-flow lemmas are reusable for other bipartite flow problems. Contributions welcome: proofs of any milestone, and a formalization of the paper's exact Δˉ\bar\DeltaΔˉ path rule showing that it satisfies the invariant.

Selected references

  • J. Edmonds, R. M. Karp, Theoretical Improvements in Algorithmic Efficiency for Network Flow Problems, Journal of the ACM 19(2):248–264, 1972. https://doi.org/10.1145/321694.321699
  • F. L. Hitchcock, The Distribution of a Product from Several Sources to Numerous Localities, Journal of Mathematics and Physics 20:224–230, 1941. https://doi.org/10.1002/sapm1941201224
  • L. R. Ford, D. R. Fulkerson, Flows in Networks, Princeton University Press, 1962. https://doi.org/10.1515/9781400875184
  • É. Tardos, A strongly polynomial minimum cost circulation algorithm, Combinatorica 5:247–255, 1985. https://doi.org/10.1007/BF02579369
  • J. B. Orlin, A faster strongly polynomial minimum cost flow algorithm, Proc. STOC 1988, 377–387. https://doi.org/10.1145/62212.62249
12 thms5 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
🏆Completed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

An Analysis of Several Heuristics for the Traveling Salesman Problem I: Nearest Neighbor Tours Can Be Far from OptimalResearch Paper

Motivation

The traveling salesman problem with the triangle inequality asks for a shortest closed tour through nnn points whose distances form a metric. It is NP-hard, so in practice tours are built by fast construction heuristics, and the natural question is how far such a tour can be from optimal in the worst case. Rosenkrantz, Stearns and Lewis (SIAM J. Comput. 6(3), 1977) gave the first systematic worst-case analysis of the standard heuristics. Their results are reproduced in textbooks on approximation algorithms and combinatorial optimization, and they are the reference point against which later guarantees (Christofides' 3/23/23/2 algorithm, the double-tree 222-approximation) are compared.

The simplest heuristic studied is the nearest neighbor algorithm (Bellmore and Nemhauser, 1968; the "next best method" of Gavett, 1965): from the current node, always move to the closest node not yet visited, and return to the start at the end. The paper shows that this greedy rule is never worse than logarithmic (Theorem 1) and that the logarithm cannot be removed (Theorem 2). This mission is about Theorem 2, the lower bound.

Setting

A traveling salesman graph on nnn nodes is a complete graph with a distance d(a,b)∈Rd(a,b)\in\mathbb Rd(a,b)∈R that is symmetric, d(a,b)=d(b,a)d(a,b)=d(b,a)d(a,b)=d(b,a), nonnegative, d(a,b)≥0d(a,b)\ge 0d(a,b)≥0, and satisfies the triangle inequality d(a,c)≤d(a,b)+d(b,c)d(a,c)\le d(a,b)+d(b,c)d(a,c)≤d(a,b)+d(b,c). A tour lists the nodes in a visiting order τ(0),…,τ(n−1)\tau(0),\dots,\tau(n-1)τ(0),…,τ(n−1) and returns to τ(0)\tau(0)τ(0); its length is the sum of the nnn distances along it. OPTIMAL is the least length of a tour.

The nearest neighbor algorithm starts at an arbitrary node τ(0)\tau(0)τ(0); having reached τ(k)\tau(k)τ(k), it moves to a node τ(k+1)\tau(k+1)τ(k+1) that minimizes d(τ(k),⋅)d(\tau(k),\cdot)d(τ(k),⋅) over the nodes not yet visited, breaking ties arbitrarily; after the last node it returns to τ(0)\tau(0)τ(0). The length of the resulting tour is written NEARNEIBER. Because the start node and the ties are free, one instance has in general several nearest-neighbor tours. A lower bound needs only one of them; an upper bound must hold for all.

The instances of the proof are built from a recursive family of weighted graphs. With li=16(4⋅2i−(−1)i+3)l_i=\frac16(4\cdot 2^i-(-1)^i+3)li​=61​(4⋅2i−(−1)i+3) (so l1,l2,l3,l4=2,3,6,11l_1,l_2,l_3,l_4=2,3,6,11l1​,l2​,l3​,l4​=2,3,6,11), the graph F1F_1F1​ is a triangle with unit weights, and Fi+1F_{i+1}Fi+1​ consists of two copies of FiF_iFi​ joined through one new node by two edges of length 111 and two edges of length lil_ili​. Each FiF_iFi​ has 2i+1−12^{i+1}-12i+1−1 nodes and a path PiP_iPi​ from its start node to its middle node through every node, of length LiL_iLi​ with L1=2L_1=2L1​=2, Li+1=2Li+2liL_{i+1}=2L_i+2l_iLi+1​=2Li​+2li​. The graph GiG_iGi​ adds two closing edges to FiF_iFi​, and Gˉi\bar G_iGˉi​ is the complete graph on the same nodes whose distance is the shortest-path distance of GiG_iGi​.

Formalization targets

Goal: Theorem 2 (p. 566)

For each m>3m>3m>3 there is a traveling salesman graph with n=2m−1n=2^m-1n=2m−1 nodes and a nearest-neighbor tour on it such that

NEARNEIBEROPTIMAL>13lg⁡(n+1)+49.\frac{\mathrm{NEARNEIBER}}{\mathrm{OPTIMAL}}>\frac13\lg(n+1)+\frac49 .OPTIMALNEARNEIBER​>31​lg(n+1)+94​.

The statement is existential in both the instance and the run of the algorithm, exactly as in the paper.

Milestones, in the order the proof uses them

  1. (2.12): the difference equation Li+1=2Li+2liL_{i+1}=2L_i+2l_iLi+1​=2Li​+2li​, L1=2L_1=2L1​=2, has the solution Li=19(6 i 2i+8⋅2i+(−1)i−9)L_i=\frac19(6\,i\,2^i+8\cdot2^i+(-1)^i-9)Li​=91​(6i2i+8⋅2i+(−1)i−9).
  2. Gˉi\bar G_iGˉi​ is a traveling salesman graph: the shortest-path distance of GiG_iGi​ is symmetric, nonnegative and satisfies the triangle inequality.
  3. (2.13)–(2.17): the shortest-path distances in Fi+1F_{i+1}Fi+1​ between the seven named nodes A,…,GA,\dots,GA,…,G of Fig. 1, e.g. AG‾=li+2−2\overline{AG}=l_{i+2}-2AG=li+2​−2.
  4. Property a): every edge of GiG_iGi​ is a shortest path between its endpoints.
  5. Property b): the nearest neighbor algorithm started at the start node of Gˉi\bar G_iGˉi​ can follow PiP_iPi​ and return along the edge of length li−1l_i-1li​−1.
  6. The optimal tour: OPTIMAL(Gˉi)=2i+1−1\mathrm{OPTIMAL}(\bar G_i)=2^{i+1}-1OPTIMAL(Gˉi​)=2i+1−1.
  7. The exact ratio: the tour along PiP_iPi​ has length Li+li−1L_i+l_i-1Li​+li​−1, so its ratio is (Li+li−1)/n(L_i+l_i-1)/n(Li​+li​−1)/n.
  8. The inequality: (Li+li−1)/n>13lg⁡(n+1)+49(L_i+l_i-1)/n>\frac13\lg(n+1)+\frac49(Li​+li​−1)/n>31​lg(n+1)+94​ for i≥3i\ge3i≥3.

The instance for mmm is Gˉm−1\bar G_{m-1}Gˉm−1​.

Significance

Theorem 1 of the same paper shows NEARNEIBER/OPTIMAL≤12⌈lg⁡n⌉+12\mathrm{NEARNEIBER}/\mathrm{OPTIMAL}\le\frac12\lceil\lg n\rceil+\frac12NEARNEIBER/OPTIMAL≤21​⌈lgn⌉+21​ for every nearest-neighbor tour on every traveling salesman graph. Theorem 2 shows that this bound has the right order: no constant-factor guarantee holds for the nearest neighbor rule, and the gap between the two constants (13\frac1331​ against 12\frac1221​) is all that remains. This separates the nearest neighbor rule from the insertion rules analysed later in the same paper, of which nearest and cheapest insertion are within a factor 222 of optimal. It is the standard example of a natural greedy heuristic whose approximation ratio grows with nnn.

The upper bound, Theorem 1, is already on Prove2Me with a machine-checked proof (SupplyChainTheory.nearest_neighbor_bound); its statement notes that the lower-bound instances are not formalized there. This mission supplies them: an explicit recursive family of metric instances, the shortest-path computations that certify it, and the arithmetic of its ratio. The result is proved in the paper; to our knowledge it has not been formalized in any proof assistant. The construction (a recursively defined weighted graph with a closed-form shortest-path table) is also a reusable pattern for other worst-case lower bounds of greedy heuristics.

Difficulty

The arithmetic ((2.12) and the final inequality) is routine. The content is in properties a) and b). A shortest-path distance is an infimum over all walks, and property a) asks that no detour through the recursive structure is shorter than the direct edge, at every level of the recursion. The paper handles this by an induction on (2.13)–(2.17) that tracks only seven nodes per level, and argues that distances inside a copy of FiF_iFi​ are not shortened by embedding it into Fi+1F_{i+1}Fi+1​. Property b) then needs that at each step of PiP_iPi​ the chosen node is at least as close as every unvisited node, including nodes in the other copy and nodes reached through the start or right nodes; ties occur, and the claim is only that some resolution of them follows PiP_iPi​. Checking small cases by computer does not give either property for all iii.

Formalization scope

Nodes of an instance are Fin n, a tour is a permutation of Fin n, the tour length is the sum over consecutive pairs including the closing edge, and OPTIMAL is a minimum over the finite set of permutations. The model is the paper's: symmetric, nonnegative distances with the triangle inequality. The distance structure also carries d(a,a)=0d(a,a)=0d(a,a)=0, a normalization not in the paper; the diagonal never enters a tour length. A nearest-neighbor tour is a permutation in which each step goes to a node at least as close as every unvisited node, from an arbitrary start with arbitrary ties.

Ratios are multiplied out: the goal is (13log⁡2(n+1)+49)⋅OPTIMAL<NEARNEIBER(\frac13\log_2(n+1)+\frac49)\cdot\mathrm{OPTIMAL}<\mathrm{NEARNEIBER}(31​log2​(n+1)+94​)⋅OPTIMAL<NEARNEIBER together with OPTIMAL>0\mathrm{OPTIMAL}>0OPTIMAL>0, the paper's standing assumption (1.1). lg⁡(n+1)\lg(n+1)lg(n+1) is Real.logb 2 of n+1n+1n+1, as printed. Because of the strict inequality and the conjunct OPTIMAL>0\mathrm{OPTIMAL}>0OPTIMAL>0, the all-zero distance does not satisfy the goal, so the statement cannot be met by a degenerate instance.

In the construction the nodes of FiF_iFi​, GiG_iGi​, Gˉi\bar G_iGˉi​ are numbered 0,…,2i+1−20,\dots,2^{i+1}-20,…,2i+1−2 from left to right (start node 000, middle node 2i−12^i-12i−1, right node 2i+1−22^{i+1}-22i+1−2); in Fi+1F_{i+1}Fi+1​ the left copy comes first, then the new node, then the right copy. Graphs are edge lists with real weights and lil_ili​ is defined in R\mathbb RR exactly as in (2.11). The shortest-path distance is the infimum of walk weights over an inductive walk predicate; it would be 000 for two nodes with no connecting walk, a case that does not arise because every GiG_iGi​ and FiF_iFi​ is connected. LiL_iLi​ is defined by its difference equation; its identification with the length of the tour along PiP_iPi​ is milestone 7. All construction statements assume i≥1i\ge1i≥1.

A complete development needs a small library for shortest-path distances of finite weighted edge lists (symmetry, triangle inequality, attainment, behaviour under relabelling and under gluing two graphs at a few nodes); this part is reusable beyond the mission. Contributions welcome: proofs of any milestone, and such general shortest-path lemmas as separate theorems. Theorem 1 is not part of this mission.

Selected references

  • D. J. Rosenkrantz, R. E. Stearns, P. M. Lewis II, An Analysis of Several Heuristics for the Traveling Salesman Problem, SIAM J. Comput. 6(3):563–581, 1977. https://doi.org/10.1137/0206041
  • M. Bellmore, G. L. Nemhauser, The Traveling Salesman Problem: A Survey, Operations Research 16(3):538–558, 1968. https://doi.org/10.1287/opre.16.3.538
  • J. W. Gavett, Three Heuristic Rules for Sequencing Jobs to a Single Production Facility, Management Science 11(8):B166–B176, 1965. https://doi.org/10.1287/mnsc.11.8.B166
  • N. Christofides, Worst-Case Analysis of a New Heuristic for the Travelling Salesman Problem, Report 388, Graduate School of Industrial Administration, Carnegie Mellon University, 1976.
12 thms4 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
🏆Completed
Operations Research·Captain: Shuze Chen

Dynamic Programming and Optimal Control II: Label Correcting MethodsTextbook

Motivation

Label correcting methods are the workhorse family of shortest-path algorithms — Dijkstra's method, Bellman–Ford, SLF/LLL variants and A* all fit the template analyzed in §2.3.1 of Bertsekas, Dynamic Programming and Optimal Control, Vol. I (3rd ed., 2005), where shortest paths appear as the purely deterministic face of dynamic programming. The correctness proof (Prop. 2.3.1) is short on paper but genuinely nondeterministic — any node may be removed from the candidate list, children processed in any order — so a formal proof certifies a whole family of concrete algorithms at once.

Setting

A finite directed graph with arc set A\mathcal{A}A, real arc lengths aija_{ij}aij​, origin sss and destination t≠st \ne st=s (BertsekasSPGraph). Walks are nonempty node lists whose consecutive pairs are arcs (BertsekasIsWalkFrom), with length the sum of arc lengths (BertsekasWalkLength); the shortest distance is the infimum of walk lengths in the extended reals, +∞+\infty+∞ if no walk exists (BertsekasShortestDistance). The standing assumption of §2.3: every cycle has nonnegative length (negative arcs allowed).

The algorithm state (BertsekasLCState) carries labels dj∈R‾d_j \in \overline{\mathbb{R}}dj​∈R, the scalar UPPER, and the candidate list OPEN. Initially ds=0d_s = 0ds​=0, all other labels ∞\infty∞, UPPER =∞= \infty=∞, OPEN ={s}= \{s\}={s}. One iteration (BertsekasLCStep, nondeterministic): remove any iii from OPEN; for each child jjj of iii in any order, if di+aij<min⁡{dj,UPPER}d_i + a_{ij} < \min\{d_j, \text{UPPER}\}di​+aij​<min{dj​,UPPER} set dj:=di+aijd_j := d_i + a_{ij}dj​:=di​+aij​, and put jjj in OPEN if j≠tj \ne tj=t, or update UPPER if j=tj = tj=t. The algorithm terminates when OPEN is empty.

Target

OPEN=∅  ⟹  UPPER=dist⁡(s,t)∈R‾,\text{OPEN} = \varnothing \implies \text{UPPER} = \operatorname{dist}(s, t) \in \overline{\mathbb{R}},OPEN=∅⟹UPPER=dist(s,t)∈R,

for every execution, under the nonnegative arc length assumption of §2.3 (aij≥0a_{ij} \ge 0aij​≥0 for every arc) — BertsekasDP.label_correcting_correctness_of_nonneg_arcs (goal). Milestones: termination — no infinite execution exists, which needs only the weaker nonnegative-cycle assumption (label_correcting_terminates) — and the workhorse invariant that every finite label is the length of an actual walk from sss, which needs neither (label_correcting_invariant).

The nonnegative-arc hypothesis is essential and not a formalization artifact: the algorithm prunes with the test di+aij<min⁡{dj,UPPER}d_i + a_{ij} < \min\{d_j, \mathrm{UPPER}\}di​+aij​<min{dj​,UPPER}, and with a negative arc a longer prefix can still reach ttt more cheaply, so the pruned node is never entered into OPEN. An earlier version of this mission's goal carried only the nonnegative-cycle assumption of §2.1 and was disproved by the counterexample s=0s=0s=0, t=2t=2t=2, a02=1a_{02}=1a02​=1, a01=2a_{01}=2a01​=2, a12=−2a_{12}=-2a12​=−2 (a graph with no cycles at all), where the algorithm terminates with UPPER=1\mathrm{UPPER}=1UPPER=1 while the shortest distance is 000. Exercise 2.7 of the source treats the nonnegative-cycle case, which requires a modified algorithm.

Significance

Prop. 2.3.1 certifies simultaneously breadth-first search, Dijkstra (best-first), depth-first and small-label-first variants — every removal discipline is one refinement of the nondeterministic relation. Formally, the development contributes a reusable small-step framework for label-setting/correcting algorithms on which sharper results (Dijkstra's single-pass property, A* admissibility, §2.3.3) can later be built. The result is classical; the formal content is the induction along the nondeterministic step relation.

Difficulty

Termination is the subtle half: labels do not decrease monotonically along the run in an obvious well-founded way; the book's argument counts the finitely many distinct walk lengths below a bound — this needs the nonnegative-cycle assumption and a careful bound relating labels to simple-path lengths. The invariant proof must thread through the fold over children within a single step.

Formalization scope

Finite node type with decidable equality; arcs as a Finset of ordered pairs; lengths total on V×VV \times VV×V (only arc values matter). The step relation is fully nondeterministic in pivot choice and child order (a permutation quantifier); correctness quantifies over all reachable terminal states — there is no fixed schedule to exploit. Distances live in EReal, so the no-path case is the honest empty infimum, not a sentinel. The trivializing risk of restricting to nonnegative arcs is avoided: only cycles are constrained.

Selected references

  • D. P. Bertsekas, Dynamic Programming and Optimal Control, Vol. I, 3rd ed., Athena Scientific, 2005. (Prop. 2.3.1, §2.3.) http://www.athenasc.com/dpbook.html
  • E. W. Dijkstra, A note on two problems in connexion with graphs, Numer. Math. 1 (1959), 269–271. https://doi.org/10.1007/BF01386390
  • R. Bellman, On a routing problem, Quart. Appl. Math. 16 (1958), 87–90. https://doi.org/10.1090/qam/102435
6 thms4 active users
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Eigenvalues and Expanders: A Regular Bipartite Graph Is a Strong Expander If and Only If λ(G) Is Bounded Away from 0Research Paper

Motivation

Expander graphs are sparse graphs in which every set of vertices has many neighbours. Families of them with bounded degree and expansion bounded away from zero are a basic tool of theoretical computer science: they are the main component of the sorting network of Ajtai, Komlós and Szemerédi (AKS 1983), the building block of superconcentrators and other graphs with strong connectivity properties, and an ingredient of many later constructions in coding theory, derandomization and complexity.

Expansion is hard to certify. Checking that every set of vertices has many neighbours means looking at exponentially many sets, and computing the exact expansion of a graph is coNP-complete. A spectral quantity, by contrast, is computable in polynomial time. N. Alon's paper Eigenvalues and expanders (Combinatorica 6 (1986) 83–96) proves that for regular bipartite graphs the two notions are equivalent: a graph is a strong expander if and only if the second-smallest eigenvalue of its Laplacian is bounded away from 0, with explicit constants in both directions. This is a discrete counterpart of Cheeger's inequality for Riemannian manifolds.

Timeline.

  • 1983: Ajtai, Komlós and Szemerédi use bounded-degree bipartite expanders to build sorting networks of depth O(log⁡n)O(\log n)O(logn).
  • 1984: Tanner (SIAM J. Alg. Disc. Meth. 5) bounds the neighbourhood size of a set in a regular bipartite graph by its second eigenvalue, the direction "eigenvalue gap implies expansion".
  • 1985: Alon and Milman (J. Combin. Theory Ser. B 38) prove isoperimetric inequalities for graphs in terms of λ(G)\lambda(G)λ(G) and introduce enlargers.
  • 1986: Alon proves the converse direction, "expansion implies an eigenvalue gap" (Lemma 2.4 and Theorem 3.4 of the paper).

Setting

All graphs are finite and simple. For a graph G=(V,E)G = (V, E)G=(V,E) and a set X⊆VX \subseteq VX⊆V, N(X)={v∈V:vx∈E for some x∈X}N(X) = \{v \in V : vx \in E \text{ for some } x \in X\}N(X)={v∈V:vx∈E for some x∈X} is the set of neighbours of XXX; it may meet XXX.

The Laplacian of GGG is QG=diag(d(v))v∈V−AGQ_G = \mathrm{diag}(d(v))_{v \in V} - A_GQG​=diag(d(v))v∈V​−AG​, where AGA_GAG​ is the 0–1 adjacency matrix and d(v)d(v)d(v) the degree of vvv. It is symmetric with eigenvalues 0=λ0≤λ1≤⋯≤λn−10 = \lambda_0 \le \lambda_1 \le \dots \le \lambda_{n-1}0=λ0​≤λ1​≤⋯≤λn−1​, counted with multiplicity, and λ(G)=λ1\lambda(G) = \lambda_1λ(G)=λ1​ is its second-smallest eigenvalue. It is positive exactly when GGG is connected.

  • An (n,d,c)(n, d, c)(n,d,c)-magnifier is a graph on nnn vertices with maximal degree ddd in which every X⊆VX \subseteq VX⊆V with ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 satisfies ∣N(X)−X∣≥c∣X∣|N(X) - X| \ge c|X|∣N(X)−X∣≥c∣X∣.
  • An (n,d,ε)(n, d, \varepsilon)(n,d,ε)-enlarger is a graph on nnn vertices with maximal degree ddd and λ(G)≥ε\lambda(G) \ge \varepsilonλ(G)≥ε.
  • A bipartite graph G=(I,O;E)G = (I, O; E)G=(I,O;E) has inputs III, outputs OOO and edges only between III and OOO. It is a strong (n,d,c)(n, d, c)(n,d,c)-expander if ∣I∣=∣O∣=n|I| = |O| = n∣I∣=∣O∣=n, the maximal degree is ddd, and for every X⊆IX \subseteq IX⊆I
∣N(X)∣≥(1+c(1−∣X∣n))∣X∣.|N(X)| \ge \Bigl(1 + c\Bigl(1 - \frac{|X|}{n}\Bigr)\Bigr)|X|.∣N(X)∣≥(1+c(1−n∣X∣​))∣X∣.

Formalization targets

Goal: Theorem 3.4

Let G=(I,O;E)G = (I, O; E)G=(I,O;E) be a ddd-regular bipartite graph with ∣I∣=∣O∣=n|I| = |O| = n∣I∣=∣O∣=n and λ=λ(G)\lambda = \lambda(G)λ=λ(G).

  1. If GGG is a strong (n,d,c)(n, d, c)(n,d,c)-expander then
λ≥c21024+2c2.\lambda \ge \frac{c^2}{1024 + 2c^2}.λ≥1024+2c2c2​.
  1. If λ≥ε\lambda \ge \varepsilonλ≥ε then GGG is a strong (n,d,c)(n, d, c)(n,d,c)-expander with
c=2dε−ε2d2.c = \frac{2d\varepsilon - \varepsilon^2}{d^2}.c=d22dε−ε2​.

Milestones, in the order of the paper

  • Lemma 2.2 (Alon–Milman, already on the platform): for disjoint sets A,BA, BA,B at distance ϱ>1\varrho > 1ϱ>1, b≤(1−a)/(1+(λ/d)aϱ2)b \le (1-a)/(1 + (\lambda/d)a\varrho^2)b≤(1−a)/(1+(λ/d)aϱ2).
  • Corollary 2.3: every (n,d,ε)(n, d, \varepsilon)(n,d,ε)-enlarger is an (n,d,2ε/(d+2ε))(n, d, 2\varepsilon/(d+2\varepsilon))(n,d,2ε/(d+2ε))-magnifier.
  • Eq. (2.1): if fff is an eigenvector of QGQ_GQG​ for λ(G)\lambda(G)λ(G) and ggg its positive part, then ∑uv∈E(g(u)−g(v))2≤λ∑vg2(v)\sum_{uv \in E}(g(u)-g(v))^2 \le \lambda \sum_v g^2(v)∑uv∈E​(g(u)−g(v))2≤λ∑v​g2(v).
  • Lemma 2.4: every (n,d,c)(n, d, c)(n,d,c)-magnifier has λ(G)≥c2/(4+2c2)\lambda(G) \ge c^2/(4 + 2c^2)λ(G)≥c2/(4+2c2).
  • Lemma 3.1: a strong (n,d,c)(n, d, c)(n,d,c)-expander is a (2n,d,c/16)(2n, d, c/16)(2n,d,c/16)-magnifier.
  • Proof of Lemma 3.3, spectrum: the two largest eigenvalues of CTCC^TCCTC, with CCC the I×OI \times OI×O biadjacency matrix, are d2d^2d2 and (d−λ)2(d - \lambda)^2(d−λ)2.
  • Proof of Lemma 3.3, Tanner's bound: ∣N(X)∣≥d2∣X∣/(α(d2−(d−λ)2)+(d−λ)2)|N(X)| \ge d^2|X| / \bigl(\alpha(d^2 - (d-\lambda)^2) + (d-\lambda)^2\bigr)∣N(X)∣≥d2∣X∣/(α(d2−(d−λ)2)+(d−λ)2) with α=∣X∣/n\alpha = |X|/nα=∣X∣/n.
  • Lemma 3.3: a ddd-regular bipartite graph is a strong (n,d,(2dλ−λ2)/d2)(n, d, (2d\lambda - \lambda^2)/d^2)(n,d,(2dλ−λ2)/d2)-expander.

Part (1) of the goal combines Lemmas 3.1 and 2.4; part (2) follows from Lemma 3.3.

Significance

The result. Theorem 3.4 makes expansion of regular bipartite graphs checkable in polynomial time up to a constant-factor loss. A random regular bipartite graph can be generated and its expansion certified by computing one eigenvalue. Lemma 2.4 is one of the first discrete Cheeger inequalities. Together with Corollary 2.3 it shows that magnifiers and enlargers are the same graphs up to the constants, and it underlies the later theory of spectral expanders, including the Alon–Boppana bound and Ramanujan graphs.

The formalization. The results are proved in the paper. The work is to formalize the known proofs. That includes Tanner's eigenvalue bound, which the paper only cites, and a max-flow min-cut argument, for which Mathlib has no general theorem. No machine-checked version of Lemma 2.4, Lemma 3.1, Lemma 3.3 or Theorem 3.4 is known to exist. Lemma 2.2 is already stated and proved on the platform as part of the Alon–Milman mission.

Difficulty

The direction "eigenvalue gap implies expansion" is a variational argument on the spectrum of CTCC^TCCTC. The converse is the hard one. A first attempt bounds λ(G)\lambda(G)λ(G) from below by testing the Rayleigh quotient on indicator vectors of sets. That only gives upper bounds on λ\lambdaλ: any one test vector does. A lower bound has to control every vector orthogonal to the constants at once. The paper first reduces to the positive part of an eigenvector (Eq. (2.1)). It then turns the combinatorial expansion of the graph into an analytic inequality for that function, using a network flow whose existence comes from the max-flow min-cut theorem. Step (ii) of the flow conditions printed on p. 87 is false as stated: the arcs (u,u)(u, u)(u,u) of the network absorb part of the flow. The flow argument has to be repaired before it can be formalized.

Lemma 3.1 has its own obstacle: one-sided expansion of inputs must be converted into expansion of arbitrary vertex sets that mix inputs and outputs. This needs the strong form of expansion; for ordinary expanders the lemma is false.

Formalization scope

Graphs are SimpleGraph V on a Fintype. A bipartite graph lives on the sum type I ⊕ O, and IsIOBipartite forbids edges inside I and inside O. Cardinalities ∣N(X)∣|N(X)|∣N(X)∣ are Set.ncard. "Maximal degree ddd" is read as G.maxDegree ≤ d; every statement is monotone in ddd or fixes ddd by regularity (G.IsRegularOfDegree d). The condition ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 is written 2∣X∣≤n2|X| \le n2∣X∣≤n in N\mathbb{N}N. All constants are real, and every subtraction and division is taken in R\mathbb{R}R.

Reused published items:

  • λ(G)\lambda(G)λ(G) is AlonMilman.Diameter.lambda1, the second-smallest eigenvalue of G.lapMatrix ℝ, which is 000 by convention on fewer than two vertices.
  • N(X)N(X)N(X) is AKSSorting.Core.neighbours.
  • Lemma 2.2 is AlonMilman.Diameter.theorem_2_5, which carries Alon–Milman's standing hypotheses that GGG is connected and n≥2n \ge 2n≥2; outside them the inequality is trivial.

The page omits a few degenerate cases, and the following hypotheses are added for them. Each is necessary, with a counterexample recorded in the item's statement:

  • n≥1n \ge 1n≥1 and c≥0c \ge 0c≥0 in Theorem 3.4 (1);
  • ε>0\varepsilon > 0ε>0 in Theorem 3.4 (2);
  • n≥2n \ge 2n≥2 and c≥0c \ge 0c≥0 in Lemma 2.4;
  • n≥2n \ge 2n≥2 in Lemma 3.1 and in the CTCC^TCCTC statement;
  • d≥1d \ge 1d≥1 in Lemma 3.3;
  • ε≥0\varepsilon \ge 0ε≥0 in Corollary 2.3.

Eq. (2.1) is stated in multiplied form, so no quotient by ∑g2\sum g^2∑g2 appears.

The closing sentences of Theorem 2.5 and Theorem 3.4 ("Thus … one can prove efficiently …") are not formalized. Read as implications between expanders they reduce to monotonicity in ccc, because c′≤cc' \le cc′≤c; their content is algorithmic.

Several encodings would trivialize the mission and are excluded:

  • a λ\lambdaλ other than the published second-smallest Laplacian eigenvalue, in particular one defined as the best constant of a quotient;
  • expansion or magnifier conditions with a negative constant in a hypothesis;
  • ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 with truncating natural-number division.

Contributions are welcome:

  • Tanner's bound in Lean, which is reusable for any regular bipartite graph;
  • a max-flow min-cut theorem for finite networks;
  • the corrected flow lemma behind Eqs. (2.2)–(2.3);
  • the spectral facts about λ(G)\lambda(G)λ(G) for bipartite graphs (λ≤d\lambda \le dλ≤d for n≥2n \ge 2n≥2, and λ=d−σ2(C)\lambda = d - \sigma_2(C)λ=d−σ2​(C)).

Selected references

  • N. Alon, Eigenvalues and expanders, Combinatorica 6 (1986) 83–96. https://doi.org/10.1007/BF02579166
  • N. Alon and V. D. Milman, λ₁, isoperimetric inequalities for graphs, and superconcentrators, J. Combin. Theory Ser. B 38 (1985) 73–88. https://doi.org/10.1016/0095-8956(85)90092-9
  • R. M. Tanner, Explicit concentrators from generalized N-gons, SIAM J. Algebraic Discrete Methods 5 (1984) 287–293. https://doi.org/10.1137/0605030
  • M. Ajtai, J. Komlós and E. Szemerédi, Sorting in c log n parallel steps, Combinatorica 3 (1983) 1–19. https://doi.org/10.1007/BF02579338
16 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
🏆Completed
CombinatoricsLinear OptimizationOperations Research+1·Captain: mikedeng1

Applied Combinatorics VIII: The Max Flow–Min Cut TheoremTextbook

Motivation

Moving as much as possible of something — freight, water, data — from an origin to a destination through connections of limited capacity is one of the basic problems of operations research. Its mathematical form, the maximum flow problem, was posed in the 1950s in work on rail networks and solved independently by Ford and Fulkerson (Maximal flow through a network, Canadian J. Math. 8 (1956)) and by Elias, Feinstein and Shannon (A note on the maximum flow through a network, IRE Trans. Inform. Theory 2 (1956)). The answer, the Max Flow–Min Cut Theorem, is a min–max duality: the largest amount that can be shipped equals the smallest total capacity whose removal disconnects the destination from the origin. It is a standard example of linear-programming duality with a combinatorial proof, and it is the source of Hall's matching theorem, Menger's theorem and Dilworth's theorem via network constructions.

This mission formalizes Chapter 13 of Keller and Trotter's Applied Combinatorics (2017 Edition), together with the two theorems of Chapter 14 that apply it, in the book's own model of a network.

Setting

A network consists of a finite vertex set VVV, a set of directed edges (x,y)(x, y)(x,y), a source SSS and a sink TTT with S≠TS \ne TS=T, and a capacity c(x,y)≥0c(x, y) \ge 0c(x,y)≥0 (a real number) on each edge. The underlying directed graph is an oriented graph: for any two vertices x,yx, yx,y at most one of (x,y)(x, y)(x,y), (y,x)(y, x)(y,x) is an edge. Every edge at SSS points away from SSS and every edge at TTT points into TTT.

A flow is a function ϕ\phiϕ on the edges with 0≤ϕ(x,y)≤c(x,y)0 \le \phi(x, y) \le c(x, y)0≤ϕ(x,y)≤c(x,y), extended by ϕ(x,y)=0\phi(x, y) = 0ϕ(x,y)=0 on pairs that are not edges, satisfying the conservation laws

∑xϕ(S,x)=∑xϕ(x,T),∑xϕ(x,y)=∑xϕ(y,x)(y≠S,T).\sum_x \phi(S, x) = \sum_x \phi(x, T), \qquad \sum_x \phi(x, y) = \sum_x \phi(y, x)\quad (y \ne S, T).x∑​ϕ(S,x)=x∑​ϕ(x,T),x∑​ϕ(x,y)=x∑​ϕ(y,x)(y=S,T).

The value of ϕ\phiϕ is value⁡(ϕ)=∑xϕ(S,x)\operatorname{value}(\phi) = \sum_x \phi(S, x)value(ϕ)=∑x​ϕ(S,x).

A cut is a partition V=L∪UV = L \cup UV=L∪U with S∈LS \in LS∈L, T∈UT \in UT∈U. Its capacity is

c(L,U)=∑x∈L, y∈Uc(x,y),c(L, U) = \sum_{x \in L,\ y \in U} c(x, y),c(L,U)=x∈L, y∈U∑​c(x,y),

summed over the edges directed from LLL to UUU only.

Given a flow ϕ\phiϕ, an edge (x,y)(x, y)(x,y) is used if ϕ(x,y)>0\phi(x, y) > 0ϕ(x,y)>0 and has spare capacity if ϕ(x,y)<c(x,y)\phi(x, y) < c(x, y)ϕ(x,y)<c(x,y). An augmenting path is a sequence P=(x0,…,xm)P = (x_0, \dots, x_m)P=(x0​,…,xm​) of distinct vertices from x0=Sx_0 = Sx0​=S to xm=Tx_m = Txm​=T such that each step either follows an edge (xi−1,xi)(x_{i-1}, x_i)(xi−1​,xi​) with spare capacity (a forward edge) or traverses a used edge (xi,xi−1)(x_i, x_{i-1})(xi​,xi−1​) backwards (a backward edge). Its augmentation amount is δ=min⁡{δ1,δ2}\delta = \min\{\delta_1, \delta_2\}δ=min{δ1​,δ2​}, where δ1\delta_1δ1​ is the least spare capacity of a forward edge and δ2\delta_2δ2​ the least flow on a backward edge (δ=δ1\delta = \delta_1δ=δ1​ when there is no backward edge).

For Chapter 14: in a finite simple graph with bipartition V=V1∪V2V = V_1 \cup V_2V=V1​∪V2​, a matching is a set of edges no two of which share an endpoint; it saturates a vertex that is an endpoint of one of its edges; and N(A)N(A)N(A) is the set of neighbors of the vertices in AAA.

Formalization targets

Goal: the Max Flow–Min Cut Theorem (Theorem 13.10)

For every network there is a real number v0v_0v0​ with

v0=max⁡{value⁡(ϕ):ϕ a flow}=min⁡{c(L,U):V=L∪U a cut},v_0 = \max\{\operatorname{value}(\phi) : \phi \text{ a flow}\} = \min\{c(L, U) : V = L \cup U \text{ a cut}\},v0​=max{value(ϕ):ϕ a flow}=min{c(L,U):V=L∪U a cut},

that is, v0v_0v0​ is attained by some flow and bounds every flow value from above, and v0v_0v0​ is attained by some cut and bounds every cut capacity from below.

Milestones

  • Theorem 13.4. For every flow ϕ\phiϕ and every cut, value⁡(ϕ)≤c(L,U)\operatorname{value}(\phi) \le c(L, U)value(ϕ)≤c(L,U).
  • Proposition 13.7. If PPP is an augmenting path for a flow ϕ\phiϕ of value vvv and δ\deltaδ is its augmentation amount, the function obtained by adding δ\deltaδ on the forward edges of PPP and subtracting δ\deltaδ on its backward edges is a flow of value v+δv + \deltav+δ.
  • Theorem 14.1. If every capacity is an integer, some maximum flow has ϕ(x,y)∈Z\phi(x, y) \in \mathbb Zϕ(x,y)∈Z on every edge.
  • Theorem 14.7 (Hall). In a finite bipartite graph with bipartition V1∪V2V_1 \cup V_2V1​∪V2​ there is a matching saturating every vertex of V1V_1V1​ if and only if ∣N(A)∣≥∣A∣|N(A)| \ge |A|∣N(A)∣≥∣A∣ for every A⊆V1A \subseteq V_1A⊆V1​.

Significance

The result. Theorem 13.10 turns every maximum-flow computation into a certified one: a flow and a cut of equal value prove each other optimal, and Theorem 13.4 shows no certificate can do better. Together with the integrality theorem 14.1 it is the engine behind the combinatorial applications of Chapter 14: maximum matchings in bipartite graphs, Hall's theorem, and the computation of the width of a poset with a minimum chain partition. Beyond the book, the same duality underlies Menger's theorem, König's theorem, the analysis of image segmentation by graph cuts, and the combinatorial theory of totally unimodular linear programs.

Formalizing it. The results are classical and proved. Mathlib has no theory of network flows. The platform has a Max-Flow Min-Cut theorem in the model of Bertsimas and Tsitsiklis (a general digraph on Fin n with capacities in (0,∞](0, \infty](0,∞], value compared in EReal), which does not cover the book's networks with zero capacities and is stated for a different encoding. This mission produces the theory in the book's model: finite oriented networks with real non-negative capacities, flows as functions on vertex pairs, cuts as vertex subsets, and the augmenting-path step that the Ford–Fulkerson labeling algorithm iterates. Hall's theorem is in Mathlib in its indexed-family form; the graph form stated here is new to the platform.

Difficulty

Theorem 13.4 is a finite-sum rearrangement. The difficulty of the goal is the existence of a maximum flow. The textbook argument runs the labeling algorithm until it halts, then reads off a cut from the labeled vertices. With real capacities this algorithm need not halt: with badly chosen augmenting paths and irrational capacities the flow values can converge to a limit strictly below the maximum, so "repeat until no augmenting path exists" does not by itself produce a maximum flow. The existence of an optimal flow is therefore not a by-product of the algorithm's description; it has to be established in its own right before the absence of augmenting paths can be turned into a cut of equal capacity. A formalization that assumes a maximum flow exists proves a strictly weaker statement. Proposition 13.7 is elementary but bookkeeping-heavy: backward edges subtract flow, and conservation must be checked at every interior vertex of the path.

Formalization scope

The vertex set is a type V with [Fintype V] [DecidableEq V]. A network (AppliedComb.Flows.Network) bundles an edge relation adj, the source S and sink T with S ≠ T, and a real capacity function cap, together with the axioms of an oriented graph, the orientation of edges at S and T, and 0 ≤ cap x y on edges. Flows are functions ϕ : V → V → ℝ satisfying IsFlow, which includes ϕ=0\phi = 0ϕ=0 off the edges and keeps the first conservation law as part of the definition, as on the page. The value is ∑xϕ(S,x)\sum_x \phi(S, x)∑x​ϕ(S,x). A cut is its part L : Finset V with S ∈ L, T ∉ L. Augmenting paths are injective maps Fin (m + 1) → V, and δ1,δ2,δ\delta_1, \delta_2, \deltaδ1​,δ2​,δ are computed in WithTop ℝ so that an empty minimum is ⊤\top⊤ and δ=δ1\delta = \delta_1δ=δ1​ when there is no backward edge. Hall's theorem uses Mathlib's SimpleGraph with a given bipartition into two Finsets and matchings as sets of Sym2 V edges.

No explicit constants arise: the chapter has no asymptotic or approximate statements.

The book's sentence of Theorem 13.10 reads "if v0v_0v0​ is the maximum value of a flow and c0c_0c0​ the minimum capacity of a cut, then v0=c0v_0 = c_0v0​=c0​". A formalization that takes v0v_0v0​ and c0c_0c0​ as hypothetical extrema of possibly empty or unattained sets would be trivial or vacuous; the goal here asserts the existence of a maximum flow and a minimum cut at the same number, and the existence of a maximum flow for real capacities is part of what must be proved.

Needed infrastructure: finite-sum manipulation over Finset (reindexing, splitting over L and Lᶜ), existence of maximizers of a linear function over the set of flows, and, for Theorem 14.1, control of integrality. The definitions of networks, flows, cuts and augmenting paths are reusable for Menger's theorem, König's theorem and the chain-partition network of Section 14.3. Contributions of alternative proofs (via linear-programming duality) are welcome.

Selected references

  • M. T. Keller and W. T. Trotter, Applied Combinatorics, 2017 Edition, Chapters 13–14. https://www.appliedcombinatorics.org/book/
  • L. R. Ford and D. R. Fulkerson, Maximal flow through a network, Canadian Journal of Mathematics 8 (1956), 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • P. Elias, A. Feinstein and C. E. Shannon, A note on the maximum flow through a network, IRE Transactions on Information Theory 2 (1956), 117–119. https://doi.org/10.1109/TIT.1956.1056816
  • P. Hall, On representatives of subsets, Journal of the London Mathematical Society 10 (1935), 26–30. https://doi.org/10.1112/jlms/s1-10.37.26
  • U. Zwick, The smallest networks on which the Ford–Fulkerson maximum flow procedure may fail to terminate, Theoretical Computer Science 148 (1995), 165–170. https://doi.org/10.1016/0304-3975(95)00022-O
8 thms3 active usersReviewed
🏆Completed
CombinatoricsProbability·Captain: mikedeng1

Applied Combinatorics VI: Ramsey's Theorem and Erdős's Lower BoundTextbook

Motivation

Ramsey theory studies the principle that complete disorder is impossible: every sufficiently large structure contains a large, perfectly uniform substructure. Its most familiar instance concerns graphs. In any graph on six vertices there are three vertices that are pairwise adjacent or three that are pairwise non-adjacent, and the same phenomenon persists at every scale. The quantity that measures it, the Ramsey number R(m,n)R(m, n)R(m,n), is one of the most studied and least understood functions in combinatorics. Only a handful of values are known exactly (R(3,3)=6R(3,3) = 6R(3,3)=6, R(4,4)=18R(4,4) = 18R(4,4)=18), while R(5,5)R(5,5)R(5,5) is only known to lie between 43 and 49 (Radziszowski, Small Ramsey Numbers).

The subject has a short, well-documented history. F. P. Ramsey proved the general theorem in 1930 as a lemma in decidability (Ramsey 1930). Erdős and Szekeres (1935) gave the upper bound R(m,n)≤(m+n−2m−1)R(m, n) \le \binom{m+n-2}{m-1}R(m,n)≤(m−1m+n−2​) (Erdős–Szekeres 1935). In 1947 Erdős proved an exponential lower bound for the diagonal numbers R(n,n)R(n, n)R(n,n) by counting graphs (Erdős 1947); this argument is now regarded as the origin of the probabilistic method. In 1959 Erdős used the same method to show that graphs of large girth and large chromatic number exist (Erdős 1959). For decades the exponential bases 2\sqrt 22​ and 444 stood essentially unchanged; the upper base was lowered below 444 only in 2023 (Campos–Griffiths–Morris–Sahasrabudhe).

This mission formalizes these results as presented in Chapter 11 of Keller and Trotter, Applied Combinatorics (2017 Edition).

Setting

A graph GGG is a finite simple graph: a finite vertex set with a symmetric, irreflexive adjacency relation (no loops, no multiple edges). A complete subgraph on mmm vertices is a set of mmm pairwise adjacent vertices; an independent set of size nnn is a set of nnn pairwise non-adjacent vertices.

A non-negative integer NNN is a Ramsey bound for (m,n)(m, n)(m,n) if every graph with at least NNN vertices contains a complete subgraph on mmm vertices or an independent set of size nnn. The Ramsey number R(m,n)R(m, n)R(m,n) is the least positive Ramsey bound. In Lean these are AppliedComb.Ramsey.IsRamseyBound m n N and AppliedComb.Ramsey.ramseyNumber m n.

More generally, write [n]={1,…,n}[n] = \{1, \dots, n\}[n]={1,…,n} and C(X,s)C(X, s)C(X,s) for the family of sss-element subsets of XXX. For a string h=(h1,…,hr)h = (h_1, \dots, h_r)h=(h1​,…,hr​), the number R(s:h1,…,hr)R(s : h_1, \dots, h_r)R(s:h1​,…,hr​) is the least positive NNN such that for every n≥Nn \ge Nn≥N and every colouring ϕ:C([n],s)→[r]\phi : C([n], s) \to [r]ϕ:C([n],s)→[r] some colour α\alphaα has a set Hα⊆[n]H_\alpha \subseteq [n]Hα​⊆[n] of size hαh_\alphahα​ all of whose sss-subsets receive colour α\alphaα (hypergraphRamseyNumber s r h).

The girth of a graph is the smallest number of vertices on a cycle, and infinite for a forest; the chromatic number χ(G)\chi(G)χ(G) is the least number of colours in a proper vertex colouring.

Formalization targets

Goal: Erdős's lower bound (Theorem 11.4)

For every positive integer nnn,

R(n,n)  ≥  ne2 2n/2.R(n, n) \;\ge\; \frac{n}{e\sqrt 2}\, 2^{n/2}.R(n,n)≥e2​n​2n/2.

Equivalently, below this threshold there is a graph on each number of vertices with neither a complete subgraph on nnn vertices nor an independent set of size nnn. The statement is for every n≥1n \ge 1n≥1, with no asymptotic slack.

Milestones

  • Lemma 11.1. Every graph with at least six vertices has a complete subgraph on 3 vertices or an independent set of size 3.
  • Theorem 11.2 (Ramsey's Theorem for Graphs). For positive integers m,nm, nm,n the least positive integer R(m,n)R(m, n)R(m,n) exists.
  • Theorem 11.6. For positive integers r,sr, sr,s and h1,…,hr≥sh_1, \dots, h_r \ge sh1​,…,hr​≥s, the least positive integer R(s:h1,…,hr)R(s : h_1, \dots, h_r)R(s:h1​,…,hr​) exists.
  • Theorem 11.7 (Erdős). For all integers g≥3g \ge 3g≥3 and ttt there is a graph with χ(G)>t\chi(G) > tχ(G)>t and girth greater than ggg.

A further item, not a milestone because the book does not number it, records the bound that the proof of Theorem 11.2 establishes: every graph with at least (m+n−2m−1)\binom{m+n-2}{m-1}(m−1m+n−2​) vertices has a complete subgraph on mmm vertices or an independent set of size nnn.

Significance

The goal is the diagonal lower bound that every later improvement is measured against. With the upper bound from the proof of Theorem 11.2 it shows that R(n,n)R(n,n)R(n,n) grows exponentially, with base between 2\sqrt 22​ and 444. No explicit construction is known to give R(n,n)>cnR(n, n) > c^nR(n,n)>cn for any constant c>1c > 1c>1. Theorem 11.7 is the standard example of a statement whose only known proofs for decades were probabilistic, and it shows that chromatic number is not a local property.

On the formal side, Mathlib has cliques, independent sets, girth and chromatic number, but no Ramsey numbers, no Erdős lower bound and no high-girth theorem. The platform already has weaker or differently shaped relatives, all checked for this mission. Erdos1947.ramsey_lower_bound gives a graph on 2⌊k/2⌋2^{\lfloor k/2 \rfloor}2⌊k/2⌋ vertices without monochromatic kkk-sets, a weaker bound than the goal's. BookSixth.high_girth_chromatic is a single-parameter form of Theorem 11.7 in another Lean environment. ramsey_theory_upper_bound is a diagonal 4k4^k4k bound. The goal statement carries the constant 1/(e2)1/(e\sqrt2)1/(e2​) exactly, which requires an explicit, non-asymptotic lower bound for n!n!n! where the book writes "Stirling's approximation".

Difficulty

The obvious route to the goal is to count graphs with a large clique or independent set and compare the result with the total number of graphs. Two steps of that route do not go through as written in the text. First, the book replaces n!n!n! by its Stirling approximation, which is only asymptotic; a statement for every n≥1n \ge 1n≥1 needs an inequality valid for all nnn, and the constant 1/(e2)1/(e\sqrt2)1/(e2​) leaves no room for a cruder estimate such as n!≥(n/e)nn! \ge (n/e)^nn!≥(n/e)n alone. Second, the counting argument yields a graph on each ttt below the threshold, whereas R(n,n)R(n,n)R(n,n) is defined as a least threshold over all graphs with at least that many vertices; the two have to be connected.

Theorem 11.7 needs random graphs with edge probability depending on nnn, a first-moment bound on short cycles and on independent sets, and a deletion step. The book states Theorem 11.6 without proof.

Formalization scope

  • Graphs are Mathlib SimpleGraph V on a finite type V : Type (Theorems 11.1, 11.2, 11.4) or on Fin N (Theorem 11.7). Cliques and independent sets are SimpleGraph.IsNClique and SimpleGraph.IsNIndepSet on a Finset.
  • IsRamseyBound m n N quantifies over every graph with at least NNN vertices, as the book does; ramseyNumber m n is the sInf of the positive Ramsey bounds. Theorem 11.2 is stated as IsLeast {N | 0 < N ∧ IsRamseyBound m n N} (ramseyNumber m n), so its content is the nonemptiness of that set. The same pattern is used for Theorem 11.6.
  • The goal compares real numbers: (n : ℝ) / (Real.exp 1 * Real.sqrt 2) * (2 : ℝ) ^ ((n : ℝ) / 2) ≤ (ramseyNumber n n : ℝ), with a real power. Explicit constant: the book says "use the Stirling approximation … after some algebra"; the statement keeps the book's constant 1/(e2)1/(e\sqrt 2)1/(e2​) and holds for every n≥1n \ge 1n≥1 with no threshold.
  • The bound of the proof of Theorem 11.2 is stated as the Ramsey property at (m+n−2m−1)\binom{m+n-2}{m-1}(m−1m+n−2​), not as an inequality on ramseyNumber, so it cannot hold through an empty defining set.
  • Girth is Mathlib's SimpleGraph.egirth (valued in N∪{∞}\mathbb N \cup \{\infty\}N∪{∞}, ∞\infty∞ for forests), not SimpleGraph.girth, which is 000 on forests. Chromatic number is SimpleGraph.chromaticNumber in N∪{∞}\mathbb N \cup \{\infty\}N∪{∞}. The parameter ttt of Theorem 11.7 is a natural number; negative ttt is trivial.
  • Theorem 11.6 is printed with typos: hi≥sh_i \ge shi​≥s is read for all i=1,…,ri = 1, \dots, ri=1,…,r, the undefined n0n_0n0​ is read as R(s:h1,…,hr)R(s : h_1, \dots, h_r)R(s:h1​,…,hr​), and C([n],s]C([n], s]C([n],s] as C([n],s)C([n], s)C([n],s). Colourings are functions on the subtype of sss-element subsets of Fin n, with colours in Fin r.
  • Trivializing formalizations ruled out. A Ramsey number defined as an arbitrary upper bound, or as a supremum with junk value 000, would make the goal vacuous or false. Here the goal's right-hand side is positive, so it forces the defining set to be nonempty, and every graph is simple on exactly its vertex type, with no loops or multiple edges.
  • Reusable infrastructure: IsRamseyBound/ramseyNumber, the hypergraph version, an all-nnn lower bound for n!n!n!, and counting over the 2(t2)2^{\binom{t}{2}}2(2t​) labelled graphs on ttt vertices. Proofs of any milestone are welcome contributions.

Selected references

  • M. T. Keller and W. T. Trotter, Applied Combinatorics, 2017 Edition, Chapter 11, pp. 229–238. https://www.appliedcombinatorics.org/
  • F. P. Ramsey, On a problem of formal logic, Proc. London Math. Soc. 30 (1930), 264–286. https://doi.org/10.1112/plms/s2-30.1.264
  • P. Erdős and G. Szekeres, A combinatorial problem in geometry, Compositio Math. 2 (1935), 463–470. http://www.numdam.org/item/CM_1935__2__463_0/
  • P. Erdős, Some remarks on the theory of graphs, Bull. Amer. Math. Soc. 53 (1947), 292–294. https://doi.org/10.1090/S0002-9904-1947-08785-1
  • P. Erdős, Graph theory and probability, Canad. J. Math. 11 (1959), 34–38. https://doi.org/10.4153/CJM-1959-003-9
  • S. Radziszowski, Small Ramsey Numbers, Electron. J. Combin. Dynamic Survey DS1. https://doi.org/10.37236/21
  • M. Campos, S. Griffiths, R. Morris and J. Sahasrabudhe, An exponential improvement for diagonal Ramsey, 2023. https://arxiv.org/abs/2303.09521
7 thms3 active usersReviewed
Convex OptimizationLinear algebraOperations Research+1·Captain: mikedeng1

Lifts of Convex Sets and Cone Factorizations III: Stable Set Polytopes Have No Small Semidefinite LiftsResearch Paper

Motivation

Many polytopes of combinatorial optimization have exponentially many facets, yet linear optimization over them is tractable because they are projections of simpler convex sets: affine slices of a nonnegative orthant (linear programming) or of the cone of positive semidefinite matrices (semidefinite programming). The size of such a representation, the number of variables of the extended formulation, is the natural measure of how compactly a polytope can be optimized over. Yannakakis (Expressing combinatorial optimization problems by linear programs, JCSS 1991) characterized polyhedral representations through nonnegative factorizations of the slack matrix. Gouveia, Parrilo and Thomas (arXiv:1111.3164, Mathematics of Operations Research 2013) extended the characterization to lifts into arbitrary closed convex cones, in particular to cones of positive semidefinite matrices.

The stable set polytope of a graph is the standard test case. For a perfect graph on nnn vertices it is a linear image of an affine slice of the cone of (n+1)×(n+1)(n+1)\times(n+1)(n+1)×(n+1) positive semidefinite matrices (Lovász's theta body construction, stated in the paper as Theorem 5.1 with a citation to Lovász & Schrijver, SIAM J. Optim. 1991); this is the reason the maximum weight stable set problem is solvable in polynomial time on perfect graphs. The question addressed by this mission is whether a smaller matrix size could suffice. Theorem 5.2 of Gouveia–Parrilo–Thomas answers it: for every graph on nnn vertices, matrices of size nnn do not suffice.

Setting

Let GGG be a graph with vertex set V={1,…,n}V = \{1,\dots,n\}V={1,…,n}. A set S⊆VS \subseteq VS⊆V is stable if no edge joins two of its elements, and its incidence vector χS∈{0,1}n\chi_S \in \{0,1\}^nχS​∈{0,1}n has (χS)i=1(\chi_S)_i = 1(χS​)i​=1 exactly when i∈Si \in Si∈S. The stable set polytope is

STAB(G)=conv{χS:S stable}⊆Rn.\mathrm{STAB}(G) = \mathrm{conv}\{\chi_S : S \text{ stable}\} \subseteq \mathbb R^n .STAB(G)=conv{χS​:S stable}⊆Rn.

Let S+k\mathcal S^k_+S+k​ be the cone of k×kk \times kk×k real symmetric positive semidefinite matrices, with the trace inner product ⟨A,B⟩=tr(AB)\langle A, B\rangle = \mathrm{tr}(AB)⟨A,B⟩=tr(AB), under which it is self-dual. For a closed convex cone KKK, a set CCC has a KKK-lift if C=π(K∩L)C = \pi(K \cap L)C=π(K∩L) for an affine subspace LLL and a linear map π\piπ; the lift is proper if LLL meets the interior of KKK.

For a polytope PPP with vertices p1,…,pvp_1,\dots,p_vp1​,…,pv​ and facet inequalities h1(x)≥0,…,hf(x)≥0h_1(x) \ge 0, \dots, h_f(x) \ge 0h1​(x)≥0,…,hf​(x)≥0, the slack matrix is the nonnegative v×fv\times fv×f matrix (hj(pi))(h_j(p_i))(hj​(pi​)). A KKK-factorization of a nonnegative matrix MMM assigns ai∈Ka^i \in Kai∈K to each row and bj∈K∗b^j \in K^*bj∈K∗ to each column with ⟨ai,bj⟩=Mij\langle a^i, b^j\rangle = M_{ij}⟨ai,bj⟩=Mij​. In Lean the objects are stab, HasPSDLift, HasConeLift, HasProperConeLift, IsSlackMatrix, HasConeFactorization and HasPSDFactorization in the namespace ConeLifts.StableSet.

Formalization targets

Goal: Theorem 5.2

For every n≥1n \ge 1n≥1 and every graph GGG on nnn vertices,

¬ ∃ L, π:STAB(G)=π(S+n∩L).\neg\ \exists\, L,\ \pi:\quad \mathrm{STAB}(G) = \pi(\mathcal S^n_+ \cap L).¬ ∃L, π:STAB(G)=π(S+n​∩L).

The statement excludes all lifts, proper or not, and holds for every graph, perfect or not.

Milestones

  1. Theorem 3.3 (first sentence). If a full-dimensional polytope PPP with the origin in its interior has a proper KKK-lift, then every slack matrix of PPP admits a KKK-factorization.
  2. Rows of the submatrix. The origin and e1,…,ene_1,\dots,e_ne1​,…,en​ are vertices of STAB(G)\mathrm{STAB}(G)STAB(G).
  3. Columns of the submatrix. For n≥1n \ge 1n≥1, each {x∈STAB(G):xi=0}\{x \in \mathrm{STAB}(G) : x_i = 0\}{x∈STAB(G):xi​=0} is a facet, and some facet does not contain the origin.
  4. The core lemma. For every s∈Rns \in \mathbb R^ns∈Rn the block matrix
S′=(10nsIn)S' = \begin{pmatrix} 1 & 0_n \\ s & I_n\end{pmatrix}S′=(1s​0n​In​​)

has no S+n\mathcal S^n_+S+n​-factorization.

Significance

Theorem 5.2 shows that the semidefinite representation of STAB(G)\mathrm{STAB}(G)STAB(G) for perfect graphs has the smallest possible matrix size: n+1n+1n+1 cannot be lowered to nnn. As Remark 5.3 of the paper notes, the same argument shows that no polytope in Rn\mathbb R^nRn with a vertex at which it locally looks like the nonnegative orthant has an S+n\mathcal S^n_+S+n​-lift. It is also an instance of the factorization method: a statement about all possible semidefinite representations is reduced to a finite obstruction on a small submatrix of the slack matrix.

The theorem is proved in the paper. To the best of our knowledge no machine-checked proof of it, of the factorization theorem for cone lifts, or of any positive semidefinite lower bound for a polytope exists in Mathlib or on this platform. A formalization produces reusable statements about positive semidefinite factorizations, slack matrices and lifts, and a verified instance of the general lower-bound technique.

Difficulty

The step from lifts to factorizations is where the direct argument fails. Theorem 3.3 applies only to proper lifts and only to polytopes with the origin in their interior, while the goal concerns all lifts of a polytope that has the origin as a vertex. Applying Theorem 3.3 to STAB(G)\mathrm{STAB}(G)STAB(G) and an arbitrary lift therefore does not match its hypotheses, and the printed proof does not spell out how the two gaps are closed (see Formalization scope). Theorem 3.3 itself is a consequence of the general factorization theorem of the paper (Theorem 2.4), whose proof rests on conic duality. The core lemma about S′S'S′ is a statement about every family of 2(n+1)2(n+1)2(n+1) positive semidefinite matrices, so it cannot be settled by any finite search.

Formalization scope

Rn\mathbb R^nRn is EuclideanSpace ℝ (Fin n); vertex i+1i+1i+1 of the paper is i : Fin n; graphs are SimpleGraph (Fin n) and stability is SimpleGraph.IsIndepSet. Vertices of a polytope are Set.extremePoints ℝ. S+k\mathcal S^k_+S+k​ is the set of real k×kk\times kk×k matrices satisfying Matrix.PosSemidef (which includes symmetry), and the ambient space of a positive semidefinite lift is all k×kk \times kk×k matrices; this does not change which sets have lifts, because a lift in the symmetric matrices extends linearly and a lift in all matrices restricts to them. Positive semidefinite factorizations require both factor families to be positive semidefinite and use tr(AiBj)\mathrm{tr}(A_iB_j)tr(Ai​Bj​).

Reading decisions: the goal assumes n≥1n \ge 1n≥1, the paper's meaning of "a graph with nnn vertices", since for n=0n = 0n=0 the polytope {0}\{0\}{0} is the image of S+0\mathcal S^0_+S+0​ and the printed statement fails. Milestone 3 also assumes n≥1n \ge 1n≥1, and so does Milestone 1 (Theorem 3.3): in R0\mathbb R^0R0 the point {0}\{0\}{0} has a proper lift to the whole space Rm\mathbb R^mRm, whose dual cone {0}\{0\}{0} cannot factor the slack matrix (1)(1)(1). In Milestone 4 the column ∗n*_n∗n​ is an arbitrary real vector. The slack matrices of Theorem 3.3 are encoded through the identification on p. 9 of the paper: rows are vertices of PPP, columns are extreme points yyy of the polar P∘={y:⟨x,y⟩≤1 ∀x∈P}P^\circ = \{y : \langle x, y \rangle \le 1\ \forall x \in P\}P∘={y:⟨x,y⟩≤1 ∀x∈P}, the canonical entry is 1−⟨p,y⟩1 - \langle p, y\rangle1−⟨p,y⟩, and every slack matrix is the canonical one with positively scaled columns. Facets in Milestone 3 are nonempty proper exposed faces of dimension one less than the polytope.

The goal must not be weakened to proper lifts, and lifts must use equality STAB(G)=π(S+n∩L)\mathrm{STAB}(G) = \pi(\mathcal S^n_+ \cap L)STAB(G)=π(S+n​∩L) with π\piπ linear and LLL affine; with inclusion, or with arbitrary maps, the statement becomes trivial or false. The core lemma is meaningful only with both factor families positive semidefinite; without that requirement S′S'S′ factors trivially.

Beyond the milestones, a complete proof of the goal needs two facts the paper uses without stating them as claims of this proof: (a) an S+n\mathcal S^n_+S+n​-lift that is not proper is a proper lift to a face of S+n\mathcal S^n_+S+n​ (p. 5), every face of S+n\mathcal S^n_+S+n​ is isomorphic to some S+r\mathcal S^r_+S+r​ with r≤nr \le nr≤n (Example 4.2, p. 12), and an S+r\mathcal S^r_+S+r​-factorization yields an S+n\mathcal S^n_+S+n​-factorization; (b) lifts are preserved by affine maps (Proposition 2.9, pp. 6–7), and translating a polytope changes its slack matrices only by positive column scalings, which is how Theorem 3.3 applies to STAB(G)\mathrm{STAB}(G)STAB(G), whose origin is a vertex rather than an interior point. Stating (a) and (b) as separate lemmas is welcome.

Needed infrastructure: positive semidefinite matrices and the trace pairing, the face structure of S+n\mathcal S^n_+S+n​, invariance of lifts under affine maps, and conic duality for Theorem 3.3. All of these are reusable beyond this mission. Contributions welcome: proofs of the milestones, the bridging facts (a) and (b), and alternative routes to the goal.

Selected references

  • J. Gouveia, P. A. Parrilo, R. R. Thomas, Lifts of Convex Sets and Cone Factorizations, Mathematics of Operations Research 38(2):248–264, 2013. arXiv:1111.3164v2. https://arxiv.org/abs/1111.3164
  • M. Yannakakis, Expressing combinatorial optimization problems by linear programs, Journal of Computer and System Sciences 43(3):441–466, 1991. https://doi.org/10.1016/0022-0000(91)90024-Y
  • L. Lovász, A. Schrijver, Cones of matrices and set-functions and 0-1 optimization, SIAM Journal on Optimization 1(2):166–190, 1991. https://doi.org/10.1137/0801013
14 thms3 active usersReviewed
🏆Completed
CombinatoricsGroup TheoryOperations Research+1·Captain: mikedeng1

λ1, Isoperimetric Inequalities for Graphs, and Superconcentrators 2: Cayley Graphs of Finite Quotients of a Property (T) Group Are Linear EnlargersResearch Paper

Motivation

An expander is a sparse graph in which every set of vertices has many neighbours outside itself. Expanders are the building blocks of superconcentrators (sparse directed graphs that route any rrr inputs to any rrr outputs along vertex-disjoint paths), of sorting and switching networks, and of many constructions in complexity theory and coding; the survey of Hoory, Linial and Wigderson (Bull. AMS 2006) describes these uses. Random regular graphs are expanders with high probability, but applications need explicit families with fixed degree and a uniform expansion constant.

Alon and Milman (J. Combin. Theory Ser. B 38 (1985)) replace combinatorial expansion by a spectral quantity, the second-smallest eigenvalue λ1\lambda_1λ1​ of the matrix Q=D−AQ = D - AQ=D−A of a graph, which Fiedler called the algebraic connectivity (Czech. Math. J. 1973). Their Theorem 4.3 shows that a regular graph with λ1\lambda_1λ1​ bounded away from 000 yields an expander. This mission formalizes their Section 4 source of such graphs: Cayley graphs of the finite quotients of a group with Kazhdan's property (T).

Timeline:

  • 1967: Kazhdan introduces property (T) and proves that SL(n,Z)SL(n,\mathbb{Z})SL(n,Z), n≥3n \ge 3n≥3, has it (Funct. Anal. Appl. 1 (1967)).
  • 1973: Margulis uses property (T) to give the first explicit expander family (Probl. Inf. Transm. 9 (1973)).
  • 1981: Gabber and Galil give a variant of Margulis' construction with an explicit expansion constant, proved by Fourier analysis (J. Comput. Syst. Sci. 22 (1981)).
  • 1985: Alon and Milman state the construction in terms of λ1\lambda_1λ1​ (Lemma 4.8, Theorem 4.9) and link it to the concentration property of Section 2 of their paper.

Setting

Let TTT be a finite group. A finite multigraph on TTT is a symmetric matrix M=(Mw,u)M = (M_{w,u})M=(Mw,u​) of nonnegative integers, Mw,uM_{w,u}Mw,u​ being the number of edges joining www and uuu (diagonal entries count loops). It is kkk-regular if every row sums to kkk. Its matrix is Q=diag⁡(d(v))−MQ = \operatorname{diag}(d(v)) - MQ=diag(d(v))−M with d(v)=∑uMv,ud(v) = \sum_u M_{v,u}d(v)=∑u​Mv,u​, a real symmetric positive semidefinite matrix. Its eigenvalues, repeated according to multiplicity, are 0=λ0≤λ1≤⋯≤λ∣T∣−10 = \lambda_0 \le \lambda_1 \le \dots \le \lambda_{|T|-1}0=λ0​≤λ1​≤⋯≤λ∣T∣−1​, and λ1(G)\lambda_1(G)λ1​(G) denotes the second of them.

Let HHH be a group, S⊆HS \subseteq HS⊆H a finite set with S=S−1S = S^{-1}S=S−1, and ϕ:H→T\phi : H \to Tϕ:H→T a homomorphism onto TTT. The Cayley multigraph G(T,ϕ(S))G(T, \phi(S))G(T,ϕ(S)) joins www and uuu by as many edges as there are s∈Ss \in Ss∈S with wu−1=ϕ(s)w u^{-1} = \phi(s)wu−1=ϕ(s). It is ∣S∣|S|∣S∣-regular, and its matrix is Q=∣S∣⋅I−∑s∈Sπ(ϕ(s))Q = |S| \cdot I - \sum_{s\in S} \pi(\phi(s))Q=∣S∣⋅I−∑s∈S​π(ϕ(s)), where π(t)\pi(t)π(t) is the permutation matrix of the left regular representation, (π(t))w,u=1(\pi(t))_{w,u} = 1(π(t))w,u​=1 iff wu−1=tw u^{-1} = twu−1=t.

An (n,k,ε)(n,k,\varepsilon)(n,k,ε)-enlarger (Definition 4.1) is a kkk-regular graph on nnn vertices with λ1≥ε\lambda_1 \ge \varepsilonλ1​≥ε.

A unitary representation π\piπ of HHH in a complex Hilbert space VVV is essentially nontrivial (Definition 4.5) if no nonzero vector is fixed by every π(h)\pi(h)π(h). A discrete group HHH has property (T) (Definition 4.6) if there are ε>0\varepsilon > 0ε>0 and a finite K⊆HK \subseteq HK⊆H such that for every essentially nontrivial unitary representation π\piπ and every unit vector yyy some h∈Kh \in Kh∈K satisfies ∣(π(h)y,y)∣<1−ε|(\pi(h)y, y)| < 1 - \varepsilon∣(π(h)y,y)∣<1−ε.

Formalization targets

Goal: Theorem 4.9

Let HHH have property (T), let SSS be a finite generating set of HHH with S=S−1S = S^{-1}S=S−1, and let ϕi:H→Ti\phi_i : H \to T_iϕi​:H→Ti​ be surjective homomorphisms onto finite groups with ∣Ti∣→∞|T_i| \to \infty∣Ti​∣→∞. Then there is one ε>0\varepsilon > 0ε>0 with

G(Ti,ϕi(S)) is a (∣Ti∣, ∣S∣, ε)-enlarger for every i with ∣Ti∣≥2.G(T_i, \phi_i(S)) \text{ is a } (|T_i|,\ |S|,\ \varepsilon)\text{-enlarger for every } i \text{ with } |T_i| \ge 2 .G(Ti​,ϕi​(S)) is a (∣Ti​∣, ∣S∣, ε)-enlarger for every i with ∣Ti​∣≥2.

The constant is unspecified: the theorem asserts uniformity in iii, not a value.

Milestones, in the order the paper's argument uses them

  1. Lemma 4.7. For any generating set SSS of a property (T) group there is ε>0\varepsilon > 0ε>0 such that every essentially nontrivial unitary representation and every unit vector yyy admit s∈Ss \in Ss∈S with ∣(π(s)y,y)∣<1−ε|(\pi(s)y,y)| < 1-\varepsilon∣(π(s)y,y)∣<1−ε.
  2. Proof of Lemma 4.8, essential nontriviality. For ϕ\phiϕ onto TTT, every nonzero vector of W={v:∑tvt=0}W = \{v : \sum_t v_t = 0\}W={v:∑t​vt​=0} is moved by some π(ϕ(h))\pi(\phi(h))π(ϕ(h)).
  3. Proof of Lemma 4.8, Rayleigh's principle. For the Cayley multigraph, min⁡{(Qy,y):y∈W, ∥y∥=1}=λ1(G)\min\{(Qy,y) : y \in W,\ \|y\| = 1\} = \lambda_1(G)min{(Qy,y):y∈W, ∥y∥=1}=λ1​(G).
  4. Lemma 4.8. With HHH, SSS and ε\varepsilonε as in Lemma 4.7, SSS finite and S=S−1S = S^{-1}S=S−1, and ϕ\phiϕ onto a finite group TTT, the Cayley graph G(T,ϕ(S))G(T,\phi(S))G(T,ϕ(S)) is a (∣T∣,∣S∣,ε)(|T|, |S|, \varepsilon)(∣T∣,∣S∣,ε)-enlarger.

Significance

Theorem 4.9 turns an analytic property of one infinite group into a uniform spectral bound for infinitely many finite graphs of fixed degree. With Theorem 4.3 of the paper it produces explicit families of linear expanders, hence of linear superconcentrators; for H=SL(n,Z)H = SL(n,\mathbb{Z})H=SL(n,Z), n≥3n \ge 3n≥3, and its reductions modulo iii the paper obtains infinitely many explicit families of (n,4,ε)(n, 4, \varepsilon)(n,4,ε)-enlargers. The same mechanism underlies later work on expanders from groups, surveyed in Lubotzky's monograph (Birkhäuser 1994).

The result is proved in the paper, modulo Lemma 4.7, which the paper refers to Margulis for. The formalization adds a checked account of every step, including Lemma 4.7 itself. At the pinned Mathlib revision there is no notion of property (T), of Kazhdan constants, or of the algebraic connectivity of a multigraph, and no machine-checked version of Theorem 4.9 is known on this platform.

Difficulty

The combinatorial and linear-algebra steps are routine; the substance is in two places. First, Lemma 4.7: Definition 4.6 supplies a constant for one finite set KKK, nothing in the definition relates KKK to a given generating set SSS, the paper gives no proof, and the standard references state the textbook form ∥π(h)y−y∥≥ε\|\pi(h)y - y\| \ge \varepsilon∥π(h)y−y∥≥ε rather than the paper's absolute-value form ∣(π(h)y,y)∣<1−ε|(\pi(h)y,y)| < 1-\varepsilon∣(π(h)y,y)∣<1−ε. Second, the uniformity: a spectral gap for each fixed quotient is easy, since a connected graph has λ1>0\lambda_1 > 0λ1​>0, but a bound that does not decay as ∣Ti∣→∞|T_i| \to \infty∣Ti​∣→∞ is exactly what cannot come from any finite computation and must come from property (T) through Lemma 4.8. The Rayleigh quotient of milestone 3 is taken over real vectors, while Lemma 4.7 is stated for complex Hilbert spaces.

Formalization scope

Groups are Lean types with a Group instance; the finite groups TTT carry Fintype and DecidableEq. Graphs are multigraphs given by symmetric matrices Matrix T T ℕ, with loops allowed. This matters: when ϕ\phiϕ identifies two generators or sends one to the identity, the degree is still ∣S∣|S|∣S∣, and a loop contributes 000 to QQQ. Mathlib's SimpleGraph Cayley graph forgets these multiplicities and is not used. λ1\lambda_1λ1​ is the second-smallest eigenvalue with multiplicity of the real symmetric matrix QQQ (via Matrix.IsHermitian.eigenvalues₀). It is defined spectrally, as in the paper, and not as a Rayleigh minimum. It is only meaningful for ∣T∣≥2|T| \ge 2∣T∣≥2, and the goal excludes trivial quotients explicitly. Unitary representations are homomorphisms into the unitary group of bounded operators on a complex Hilbert space in universe Type. Compact subsets of a discrete group are finite sets.

A trivializing formalization is ruled out. Property (T) is not replaced by the hypothesis that the regular representations of the quotients have no almost-invariant vectors, which would make Theorem 4.9 a restatement of its hypothesis. And λ1\lambda_1λ1​ is not defined as the minimum of (Qy,y)(Qy,y)(Qy,y) over zero-sum unit vectors, which would make milestone 3 true by definition.

A complete development needs basic Kazhdan-constant manipulations, Courant–Fischer for real symmetric matrices, and the regular representation of a finite group as a unitary representation. The regularity of Cayley multigraphs and the identity Q=∣S∣I−∑sπ(ϕ(s))Q = |S| I - \sum_s \pi(\phi(s))Q=∣S∣I−∑s​π(ϕ(s)) are short. The spectral and representation-theoretic lemmas are reusable beyond this mission. Contributions of intermediate lemmas, such as the variational characterization of eigenvalues₀ or the invariance of the zero-sum subspace, are welcome.

Selected references

  • N. Alon, V. D. Milman, λ1, Isoperimetric inequalities for graphs, and superconcentrators, J. Combin. Theory Ser. B 38 (1985) 73–88. https://doi.org/10.1016/0095-8956(85)90092-9
  • D. A. Kazhdan, Connection of the dual space of a group with the structure of its closed subgroups, Funct. Anal. Appl. 1 (1967) 63–65. https://doi.org/10.1007/BF01075866
  • G. A. Margulis, Explicit constructions of concentrators, Probl. Inf. Transm. 9 (1973) 325–332. http://mi.mathnet.ru/ppi1162
  • O. Gabber, Z. Galil, Explicit constructions of linear-sized superconcentrators, J. Comput. Syst. Sci. 22 (1981) 407–420. https://doi.org/10.1016/0022-0000(81)90040-4
  • M. Fiedler, Algebraic connectivity of graphs, Czech. Math. J. 23 (1973) 298–305. https://doi.org/10.21136/CMJ.1973.101168
  • A. Lubotzky, Discrete Groups, Expanding Graphs and Invariant Measures, Birkhäuser, 1994. https://doi.org/10.1007/978-3-0346-0332-4
  • B. Bekka, P. de la Harpe, A. Valette, Kazhdan's Property (T), Cambridge University Press, 2008. https://doi.org/10.1017/CBO9780511542749
  • S. Hoory, N. Linial, A. Wigderson, Expander graphs and their applications, Bull. AMS 43 (2006) 439–561. https://doi.org/10.1090/S0273-0979-06-01126-8
12 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations Research·Captain: mikedeng1

λ1, Isoperimetric Inequalities for Graphs, and Superconcentrators 1: A Diameter Bound from λ1Research Paper

Motivation

The eigenvalues of the Laplacian of a graph carry metric information about the graph. The second-smallest one, λ1(G)\lambda_1(G)λ1​(G), was named the algebraic connectivity by Fiedler (Fiedler 1973), who showed it is positive exactly for connected graphs. N. Alon and V. D. Milman (J. Combin. Theory Ser. B 38 (1985) 73–88) showed that a large λ1\lambda_1λ1​ also forces two further properties: small diameter and a concentration of measure phenomenon, in which almost every vertex is close to any set containing half the vertices. They used these facts to build explicit expanders and superconcentrators, which are sparse networks with strong connectivity guarantees used in the theory of computation and in communication network design.

This mission covers Section 2 of that paper, "The Main Tools": the edge-count inequality (Lemma 2.1), the isoperimetric inequalities (Theorems 2.5 and 2.6), and the resulting diameter bound (Theorem 2.7).

Timeline:

  • 1973: Fiedler introduces λ1(G)\lambda_1(G)λ1​(G) as algebraic connectivity and proves λ1≤nn−1min⁡vd(v)\lambda_1 \le \frac{n}{n-1}\min_v d(v)λ1​≤n−1n​minv​d(v).
  • 1985: Alon and Milman prove the isoperimetric and diameter bounds of Section 2.
  • 1986: Alon proves the converse direction, that edge expansion implies a spectral gap (Alon 1986).
  • Later work sharpened the constant in the diameter bound, e.g. Chung 1989.

Setting

Let G=(V,E)G = (V, E)G=(V,E) be a finite, connected, simple graph on n=∣V∣≥2n = |V| \ge 2n=∣V∣≥2 vertices. Write d(v)d(v)d(v) for the degree of a vertex vvv, d=max⁡vd(v)d = \max_v d(v)d=maxv​d(v) for the maximum degree, and AGA_GAG​ for the adjacency matrix. The Laplacian is the V×VV \times VV×V matrix

Q=QG=diag⁡(d(v))v∈V−AG.Q = Q_G = \operatorname{diag}(d(v))_{v \in V} - A_G .Q=QG​=diag(d(v))v∈V​−AG​.

For real functions fff on VVV with scalar product (f,g)=∑vf(v)g(v)(f, g) = \sum_v f(v)g(v)(f,g)=∑v​f(v)g(v), the quadratic form of QQQ is (Qf,f)=∑{u,v}∈E(f(u)−f(v))2≥0(Qf, f) = \sum_{\{u,v\} \in E} (f(u) - f(v))^2 \ge 0(Qf,f)=∑{u,v}∈E​(f(u)−f(v))2≥0. The eigenvalues of QQQ, counted with multiplicity, are real and are written 0=λ0≤λ1≤⋯≤λn−10 = \lambda_0 \le \lambda_1 \le \dots \le \lambda_{n-1}0=λ0​≤λ1​≤⋯≤λn−1​. The algebraic connectivity λ1=λ1(G)\lambda_1 = \lambda_1(G)λ1​=λ1​(G) is the second-smallest of them.

For vertices u,vu, vu,v, dist⁡(u,v)\operatorname{dist}(u, v)dist(u,v) is the number of edges of a shortest path from uuu to vvv. For disjoint vertex sets A,BA, BA,B the paper writes ρ\rhoρ for the distance between them, a=∣A∣/na = |A|/na=∣A∣/n and b=∣B∣/nb = |B|/nb=∣B∣/n for their relative sizes, and EAE_AEA​ (EBE_BEB​) for the set of edges with both endpoints in AAA (in BBB). [x][x][x] denotes the integer part of x≥0x \ge 0x≥0.

Formalization targets

Goal: Theorem 2.7 (p. 79)

dist⁡(u,v)  ≤  2[2d/λ1 log⁡2n]for all u,v∈V.\operatorname{dist}(u, v) \;\le\; 2\left[\sqrt{2d/\lambda_1}\,\log_2 n\right] \qquad\text{for all } u, v \in V.dist(u,v)≤2[2d/λ1​​log2​n]for all u,v∈V.

Milestones, in the order the proof uses them

  1. Section 2, p. 76: 0=λ0<λ10 = \lambda_0 < \lambda_10=λ0​<λ1​ for connected GGG.
  2. Eq. (2.1), Rayleigh's principle: if ∑vf(v)=0\sum_v f(v) = 0∑v​f(v)=0 then (Qf,f)≥λ1∥f∥2(Qf, f) \ge \lambda_1 \|f\|^2(Qf,f)≥λ1​∥f∥2.
  3. Lemma 2.1: for nonempty A,BA, BA,B at distance ρ≥1\rho \ge 1ρ≥1,
λ1n≤1ρ2(1a+1b)(∣E∣−∣EA∣−∣EB∣).\lambda_1 n \le \frac{1}{\rho^2}\Big(\frac1a + \frac1b\Big)\big(|E| - |E_A| - |E_B|\big).λ1​n≤ρ21​(a1​+b1​)(∣E∣−∣EA​∣−∣EB​∣).
  1. Remark 2.3: λ1≤nn−1min⁡vd(v)\lambda_1 \le \frac{n}{n-1}\min_v d(v)λ1​≤n−1n​minv​d(v).
  2. Theorem 2.5: if ρ>1\rho > 1ρ>1 then
b≤1−a1+(λ1/d) aρ2.b \le \frac{1-a}{1 + (\lambda_1/d)\,a\rho^2}.b≤1+(λ1​/d)aρ21−a​.
  1. Theorem 2.6: if every AAA–BBB distance exceeds a real ρ≥1\rho \ge 1ρ≥1, then
b≤(1−a)exp⁡ ⁣(−ln⁡(1+2a)[λ1/(2d) ρ]).b \le (1-a)\exp\!\Big(-\ln(1+2a)\Big[\sqrt{\lambda_1/(2d)}\,\rho\Big]\Big).b≤(1−a)exp(−ln(1+2a)[λ1​/(2d)​ρ]).

Each statement keeps the paper's explicit constants. The goal is the endpoint of this chain and the paper's headline graph-theoretic bound.

Significance

Theorem 2.7 gives, for any family of graphs of bounded maximum degree whose algebraic connectivity stays bounded away from zero, a diameter of order log⁡n\log nlogn. By the paper's Remark 2.8, the 4-regular graphs constructed in its Section 4 show that this order cannot be improved. Theorem 2.6 is a discrete concentration of measure inequality: the proportion of vertices at distance more than ρ\rhoρ from a set of relative size aaa decays exponentially in ρλ1/(2d)\rho\sqrt{\lambda_1/(2d)}ρλ1​/(2d)​. It is the graph analogue of the Gromov–Milman concentration for manifolds, and Section 3 of the paper applies it to cubes and other product graphs. Theorem 2.5 is the input for the construction of expanders from graphs with a spectral gap (Theorem 4.3 of the paper).

All results are proved in the paper, and the formal work here is a machine-checked version of known proofs. As far as could be determined, none of the four inequalities (Lemma 2.1, Theorems 2.5–2.7) has been formalized in Lean or elsewhere. Mathlib has the Laplacian matrix, its positive semidefiniteness, and the relation between its kernel and connected components, but no statement about its second eigenvalue. The spectral facts (milestones 1–2), stated for Mathlib's Matrix.IsHermitian.eigenvalues₀, are reusable for any future work on algebraic connectivity.

Difficulty

The combinatorial steps are short. The work is at the interface between the spectral definition and the quadratic form. Mathlib defines eigenvalues through the spectral theorem for a Hermitian matrix, sorted into a list. Obtaining Rayleigh's principle for the second eigenvalue from that list, with the constant functions as the eigenvector of λ0=0\lambda_0 = 0λ0​=0, takes a Courant–Fischer-type argument over an orthonormal eigenbasis. It does not follow from positive semidefiniteness alone. Strict positivity of λ1\lambda_1λ1​ additionally needs that the kernel of QQQ is one-dimensional for a connected graph.

Theorem 2.6 iterates Theorem 2.5 over a sequence of neighbourhoods {v:dist⁡(v,A)≤jμ}\{v : \operatorname{dist}(v, A) \le j\mu\}{v:dist(v,A)≤jμ} with a real step length μ\muμ, so it needs bookkeeping of integer parts and of real-valued distance thresholds. Theorem 2.7 then combines Theorem 2.6 with Remark 2.3 and needs the estimate 12 2−[log⁡2n]<1/n\tfrac12\, 2^{-[\log_2 n]} < 1/n21​2−[log2​n]<1/n with the integer part kept. Replacing [⋅][\cdot][⋅] by the real number inside it changes the statement.

Formalization scope

  • Graphs are Mathlib SimpleGraph V on a Fintype vertex type with decidable adjacency. Every item assumes G.Connected and 2≤∣V∣2 \le |V|2≤∣V∣ (the goal writes 1<∣V∣1 < |V|1<∣V∣, as the paper does).
  • QQQ is G.lapMatrix ℝ. λ1\lambda_1λ1​ is the mission definition AlonMilman.Diameter.lambda1: the eigenvalue at index n−2n-2n−2 of eigenvalues₀, which lists the eigenvalues in decreasing order. It is 000 by convention when n<2n < 2n<2, a case no theorem uses.
  • λ1\lambda_1λ1​ is defined spectrally. Defining it as the best constant in Eq. (2.1) would make Rayleigh's principle definitional and remove the spectral content of the mission, so that formalization is excluded. Likewise the goal quantifies over all pairs of vertices of a connected graph and does not use SimpleGraph.diam without connectivity, since that is 000 for a disconnected graph.
  • Distances are SimpleGraph.dist (a natural number). "The distance between AAA and BBB is ρ\rhoρ" is encoded as ρ≤dist⁡(u,v)\rho \le \operatorname{dist}(u, v)ρ≤dist(u,v) for all u∈Au \in Au∈A, v∈Bv \in Bv∈B. Because the bounds weaken as ρ\rhoρ decreases, this is equivalent to the paper's exact distance. In Theorem 2.6 ρ\rhoρ is real and the hypothesis is strict.
  • EAE_AEA​ is AlonMilman.Diameter.edgesWithin G A. All counts are cast to R\mathbb RR before subtraction, a=∣A∣/na = |A|/na=∣A∣/n is a real quotient, [x][x][x] is Nat.floor, log⁡2\log_2log2​ is Real.logb 2, and ln⁡\lnln is Real.log.
  • Lemma 2.1 requires A,BA, BA,B nonempty (so a,b>0a, b > 0a,b>0). Theorems 2.5 and 2.6 hold as stated for empty sets and carry no such hypothesis.

Contributions welcome: proofs of any milestone, and in particular general Mathlib-style lemmas for Rayleigh quotients and eigenvalues₀, which have uses beyond this mission.

Selected references

  • N. Alon, V. D. Milman, λ1, isoperimetric inequalities for graphs, and superconcentrators, J. Combin. Theory Ser. B 38 (1985) 73–88. https://doi.org/10.1016/0095-8956(85)90092-9
  • M. Fiedler, Algebraic connectivity of graphs, Czechoslovak Math. J. 23 (1973) 298–305. https://doi.org/10.21136/CMJ.1973.101168
  • N. Alon, Eigenvalues and expanders, Combinatorica 6 (1986) 83–96. https://doi.org/10.1007/BF02579166
  • F. R. K. Chung, Diameters and eigenvalues, J. Amer. Math. Soc. 2 (1989) 187–196. https://doi.org/10.1090/S0894-0347-1989-0965008-X
9 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchTheoretical Computer Science·Captain: mikedeng1

A Linear-Time Algorithm for Finding a Sparse k-Connected Spanning Subgraph of a k-Connected Graph 2: FOREST's Forests E_i Preserve Local Node-Connectivity up to i in a Simple GraphResearch Paper

Motivation

Given a kkk-connected graph, many connectivity algorithms run in time that grows with the number of edges ∣E∣|E|∣E∣. A sparse certificate is a spanning subgraph with only O(k∣V∣)O(k|V|)O(k∣V∣) edges that is still kkk-connected; computing one first and running the expensive algorithm on it replaces ∣E∣|E|∣E∣ by k∣V∣k|V|k∣V∣ in the bound. Finding a kkk-connected spanning subgraph with the minimum number of edges is NP-complete for every fixed k≥2k \ge 2k≥2 (Garey and Johnson, problem GT31), so the question is how cheaply a sparse, not necessarily minimum, certificate can be found.

Nagamochi and Ibaraki (Algorithmica 7 (1992) 583–596) answered this with a single linear-time scanning procedure, FOREST, which partitions the edges into classes E1,E2,…,E∣E∣E_1, E_2, \dots, E_{|E|}E1​,E2​,…,E∣E∣​. They showed that the prefix unions E1∪⋯∪EkE_1 \cup \dots \cup E_kE1​∪⋯∪Ek​ are certificates for edge-connectivity and, for simple graphs, for node-connectivity. The node-connectivity result is the subject of this mission; the edge-connectivity result is the preceding mission of this series.

Timeline:

  • 1980: Galil gives an algorithm testing κ(G)≥k\kappa(G) \ge kκ(G)≥k whose running time depends on ∣E∣|E|∣E∣ (SIAM J. Comput. 9).
  • Before 1992 (as cited on p. 583): Suzuki et al. give O(∣E∣)O(|E|)O(∣E∣)-time algorithms for sparse 2- and 3-node-connected spanning subgraphs; Nishizeki and Poljak find, for general kkk, a kkk-node-connected spanning subgraph with at most k(∣V∣−1)k(|V|-1)k(∣V∣−1) edges in O(∣V∣1/2∣E∣2)O(|V|^{1/2}|E|^2)O(∣V∣1/2∣E∣2) time.
  • 1992: Nagamochi and Ibaraki prove that FOREST, which runs in O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣) time, yields a kkk-node-connected spanning subgraph GkG_kGk​ of every simple kkk-node-connected graph, with ∣E(Gk)∣≤k∣V∣−k(k+1)/2|E(G_k)| \le k|V| - k(k+1)/2∣E(Gk​)∣≤k∣V∣−k(k+1)/2. Their Theorem 3.1 states a stronger, local form.
  • 1993: Cheriyan, Kao and Thurimella isolate "scan-first search" as the general principle behind such certificates (SIAM J. Comput. 22 (1993)).

Setting

A graph G=(V,E)G = (V, E)G=(V,E) has a finite node set VVV with ∣V∣≥2|V| \ge 2∣V∣≥2 and a finite edge set EEE; each edge has an unordered pair of two distinct end nodes. In this mission the graph is simple: no two edges have the same end nodes. For F⊆EF \subseteq EF⊆E, (V,F)(V, F)(V,F) is the spanning subgraph with edge set FFF.

The local node-connectivity κ(x,y;H)\kappa(x, y; H)κ(x,y;H) of nodes x,yx, yx,y in a graph HHH on VVV is ∣V∣−1|V| - 1∣V∣−1 if xxx and yyy are adjacent in HHH, and otherwise the minimum size of a node set W⊆V−{x,y}W \subseteq V - \{x, y\}W⊆V−{x,y} whose deletion leaves no xxx–yyy path. The node connectivity is κ(G)=min⁡x,yκ(x,y;G)\kappa(G) = \min_{x, y} \kappa(x, y; G)κ(G)=minx,y​κ(x,y;G).

Procedure FOREST keeps a label r(v)≥0r(v) \ge 0r(v)≥0 on each node, initially 000. While some node is unscanned, it chooses an unscanned node xxx of largest label; for each unscanned edge e=(x,y)e = (x, y)e=(x,y) it puts eee into the class Er(y)+1E_{r(y)+1}Er(y)+1​, increases r(x)r(x)r(x) by one if r(x)=r(y)r(x) = r(y)r(x)=r(y), and increases r(y)r(y)r(y) by one; then it marks xxx scanned. Ties are broken arbitrarily. The time instants are the states between these elementary operations; Ei∗E^*_iEi∗​ denotes the class iii at an instant, and EiE_iEi​ its final value. Put

Gi=(V, E1∪E2∪⋯∪Ei).G_i = (V,\ E_1 \cup E_2 \cup \dots \cup E_i).Gi​=(V, E1​∪E2​∪⋯∪Ei​).

Formalization targets

Goal: Theorem 3.1

For a simple graph GGG and the classes of any completed run of FOREST, for 1≤i≤∣E∣1 \le i \le |E|1≤i≤∣E∣,

κ(x,y;Gi) ≥ min⁡{κ(x,y;G), i}for any x,y∈V.(3.1)\kappa(x, y; G_i) \ \ge\ \min\{\kappa(x, y; G),\ i\} \qquad \text{for any } x, y \in V. \tag{3.1}κ(x,y;Gi​) ≥ min{κ(x,y;G), i}for any x,y∈V.(3.1)

The statement is local: it holds pair by pair, not only for the global minimum, and for every tie-breaking of the procedure.

Milestones, in the order the proof uses them

  1. Lemma 2.2: at every instant, a node vvv meets EiE_iEi​ exactly for i=1,…,r(v)i = 1, \dots, r(v)i=1,…,r(v).
  2. Lemma 2.4(b): at every instant, a uuu–vvv path in Ej∗E^*_jEj∗​ yields uuu–vvv paths in every Ei∗E^*_iEi∗​, i<ji < ji<j.
  3. In-degree at most one (§2, p. 588): orienting each edge from the earlier-scanned to the later-scanned end, every node has at most one entering arc in each class.
  4. Lemma 3.1: if an xxx–yyy path of Ej∗E^*_jEj∗​ has the form x,u1,…,uk=w,yx, u_1, \dots, u_k = w, yx,u1​,…,uk​=w,y with k=1k = 1k=1 or u1u_1u1​ scanned before www, then any www–xxx and www–yyy paths in Ei∗E^*_iEi∗​ (i<ji < ji<j) share a node other than www.
  5. Lemma 3.2: for a node cut set W={w1,…,wi}W = \{w_1, \dots, w_i\}W={w1​,…,wi​} of Gi+1G_{i+1}Gi+1​ (in scan order) separating a component XXX from the rest YYY, immediately after wtw_twt​ is scanned every XXX–YYY path of Et∗E^*_tEt∗​ passes through wtw_twt​, and Ej∗E^*_jEj∗​ has no XXX–YYY path for t+1≤j≤i+1t + 1 \le j \le i + 1t+1≤j≤i+1.

A companion item states the paper's announcement in §3: GkG_kGk​ is kkk-node-connected for every 1≤k≤κ(G)1 \le k \le \kappa(G)1≤k≤κ(G).

Significance

Theorem 3.1 at i=ki = ki=k shows that the first kkk classes of FOREST form a kkk-node-connected spanning subgraph whenever GGG is, and the edge-count analysis of the companion mission bounds its size by k∣V∣−k(k+1)/2k|V| - k(k+1)/2k∣V∣−k(k+1)/2. Since FOREST runs in linear time, any algorithm testing κ(G)≥k\kappa(G) \ge kκ(G)≥k can be run on GkG_kGk​ instead of GGG; the paper uses this to improve the bound for testing κ(G)≥k\kappa(G) \ge kκ(G)≥k from O(max⁡{k2∣V∣1/2,k∣V∣}∣E∣)O(\max\{k^2|V|^{1/2}, k|V|\}|E|)O(max{k2∣V∣1/2,k∣V∣}∣E∣) to O(max⁡{k3∣V∣3/2,k2∣V∣2})O(\max\{k^3|V|^{3/2}, k^2|V|^2\})O(max{k3∣V∣3/2,k2∣V∣2}), and similar gains for computing the number of node-disjoint paths between two nodes. The local form (3.1) is what makes the sss–ttt applications possible.

The result is proved on paper. As far as a search of the platform shows, neither FOREST nor local node-connectivity has a machine-checked treatment there; Mathlib has no notion of vertex connectivity of a pair of nodes. This mission produces a formal model of FOREST as a nondeterministic transition system and the statements needed to verify the paper's proof step by step.

Difficulty

For edge-connectivity the analogous statement follows from a general principle: any sequence of maximal spanning forests, each taken in what remains of the graph, preserves local edge-connectivity up to its length. The obvious attempt is to prove (3.1) the same way, from the fact that each EiE_iEi​ is a maximal spanning forest of what the earlier classes leave. The paper gives no such argument for node-connectivity: its proof uses the specific scan order of FOREST in an essential way, through the orientation of edges from earlier- to later-scanned nodes and the in-degree bound of that orientation (Lemmas 3.1 and 3.2). The argument tracks, for a hypothetical node cut WWW of size iii in Gi+1G_{i+1}Gi+1​, the classes at the moments the nodes of WWW are scanned, which requires reasoning about intermediate states of the algorithm and about paths in several classes at once. None of this reduces to a static property of the output partition.

Formalization scope

  • Graphs. A node type V and an edge type E, both finite, with ends : E → Sym2 V; loop-freeness is ∀ e, ¬ (ends e).IsDiag and simplicity is Function.Injective ends. The standing assumptions of p. 583 and p. 589 (∣V∣≥2|V| \ge 2∣V∣≥2, no self-loop, a simple graph when node-connectivity is discussed) appear as hypotheses; §2 items are stated for loopless graphs, as on the page.
  • Connectivity. κ(x,y;(V,F))\kappa(x, y; (V, F))κ(x,y;(V,F)) is valued in N∞\mathbb N_\inftyN∞​: ∣V∣−1|V| - 1∣V∣−1 on adjacent pairs, the minimum node cut otherwise, and ⊤\top⊤ when x=yx = yx=y, where (3.1) holds trivially.
  • FOREST. A nondeterministic step relation with three steps (select, scan, finish), a run of length KKK from the initial state, and completion when every node is scanned. Every tie-breaking is allowed, so the theorems quantify over all completed runs. A time instant is a state of a run; "scanned before" compares positions in the run's selection order; "immediately after wtw_twt​ has been scanned" is the state right after the finish step of wtw_twt​.
  • Excluded. The running time "O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣)" and everything in §4 (the connectivity-testing algorithms and their bounds) are not stated: the paper fixes no machine model. The edge bounds on ∣Ei∣|E_i|∣Ei​∣ belong to the edge-connectivity mission.
  • Ruled out. Stating (3.1) for an arbitrary partition into maximal spanning forests, or reading the classes Ej∗E^*_jEj∗​ of Lemmas 3.1–3.2 off the final state, would state a different theorem from the one the paper proves; the classes are those of a run of FOREST at the instant the page specifies.
  • Infrastructure. Reachability avoiding a node set, walks and paths in SimpleGraph, and invariants of the FOREST transition system. The definitions duplicate those of the edge-connectivity mission by design and are candidates for a shared layer. Contributions of general lemmas about the run (label invariants, monotonicity of classes along a run) are welcome.

Selected references

  • H. Nagamochi, T. Ibaraki, A linear-time algorithm for finding a sparse kkk-connected spanning subgraph of a kkk-connected graph, Algorithmica 7 (1992), 583–596. https://doi.org/10.1007/BF01758778
  • Z. Galil, Finding the vertex connectivity of graphs, SIAM J. Comput. 9 (1980), 197–199. https://doi.org/10.1137/0209016
  • J. Cheriyan, M.-Y. Kao, R. Thurimella, Scan-first search and sparse certificates: an improved parallel algorithm for kkk-vertex connectivity, SIAM J. Comput. 22 (1993), 157–174. https://doi.org/10.1137/0222013
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, 1979.
9 thms3 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchTheoretical Computer Science·Captain: mikedeng1

A Linear-Time Algorithm for Finding a Sparse k-Connected Spanning Subgraph of a k-Connected Graph 1: The FOREST Decomposition Preserves Local Edge-ConnectivityResearch Paper

Motivation

Many graph algorithms for connectivity questions run in time proportional to the number of edges. When the question is only whether a graph is kkk-edge-connected, or what its local edge-connectivities are up to a threshold kkk, most edges are irrelevant: a spanning subgraph with O(k∣V∣)O(k|V|)O(k∣V∣) edges already carries the answer. Such a subgraph is called a sparse certificate. Computing one first and running the expensive algorithm on it replaces ∣E∣|E|∣E∣ by k∣V∣k|V|k∣V∣ in the running time of connectivity testing, of Matula-type edge-connectivity algorithms, and of sss–ttt flow computations used as connectivity oracles.

Nagamochi and Ibaraki (Algorithmica 7, 1992) gave a procedure, FOREST, that computes such a certificate for edge-connectivity and, on simple graphs, for node-connectivity, with a single graph search. The same partition of the edges into forests is the engine of their deterministic minimum-cut algorithm for multigraphs (SIAM J. Discrete Math. 5, 1992), which later became the maximum-adjacency ordering of the Stoer–Wagner minimum-cut algorithm (J. ACM 44, 1997).

Timeline:

  • 1927: Menger identifies the minimum number of edges separating two nodes with the maximum number of edge-disjoint paths between them.
  • 1992: Nagamochi and Ibaraki publish Procedure FOREST and prove that its iii-th prefix preserves local edge-connectivity up to iii in multigraphs, and local node-connectivity up to iii in simple graphs.
  • 1993: Cheriyan, Kao and Thurimella (SIAM J. Comput. 22) obtain sparse certificates by scan-first search; Frank, Ibaraki and Nagamochi (J. Graph Theory 17) give a shorter proof for the node-connectivity case.
  • 1994: Nishizeki and Poljak (Discrete Appl. Math. 55) publish the forest-decomposition lemma (Lemma 2.1 below), found independently.

Setting

A graph G=(V,E)G = (V, E)G=(V,E) is finite and undirected, has ∣V∣≥2|V| \ge 2∣V∣≥2 nodes, may have multiple edges (several edges with the same pair of end nodes), and has no self-loop. It is simple if no two edges have the same end nodes. For F⊆EF \subseteq EF⊆E, (V,F)(V, F)(V,F) is the spanning subgraph with edge set FFF. It is a forest if it has no cycle; two parallel edges form a cycle. It is a maximal spanning forest in (V,H)(V, H)(V,H), for F⊆HF \subseteq HF⊆H, if adding any edge of H∖FH \setminus FH∖F to FFF creates a cycle.

The local edge-connectivity λ(x,y;H)\lambda(x, y; H)λ(x,y;H) is the minimum number of edges of HHH whose removal leaves no path from xxx to yyy; parallel edges count separately. It is ∞\infty∞ when x=yx = yx=y.

Procedure FOREST keeps a label r(v)∈Nr(v) \in \mathbb{N}r(v)∈N on every node, initially 000, and classes E1,E2,…,E∣E∣E_1, E_2, \dots, E_{|E|}E1​,E2​,…,E∣E∣​, initially empty. While an unscanned node exists, it picks an unscanned node xxx of largest label; for every unscanned edge eee from xxx to a node yyy, it puts eee into Er(y)+1E_{r(y)+1}Er(y)+1​, increases r(x)r(x)r(x) by one if r(x)=r(y)r(x) = r(y)r(x)=r(y), increases r(y)r(y)r(y) by one, and marks eee scanned; then it marks xxx scanned. Ties among nodes and the order of edges are free. On termination, Gi=(V,E1∪⋯∪Ei)G_i = (V, E_1 \cup \cdots \cup E_i)Gi​=(V,E1​∪⋯∪Ei​).

Formalization targets

Goal: Theorem 2.1 (pp. 588–589), without the running time

For every graph GGG and every completed execution of FOREST on GGG:

  1. every edge lies in exactly one class EiE_iEi​, 1≤i≤∣E∣1 \le i \le |E|1≤i≤∣E∣;
  2. for i=1,…,∣E∣i = 1, \dots, |E|i=1,…,∣E∣,
λ(x,y;Gi)≥min⁡{λ(x,y;G), i}for all x,y∈V;(2.1)\lambda(x, y; G_i) \ge \min\{\lambda(x, y; G),\ i\} \qquad \text{for all } x, y \in V; \tag{2.1}λ(x,y;Gi​)≥min{λ(x,y;G), i}for all x,y∈V;(2.1)
  1. ∣Ei∣≤∣V∣−1|E_i| \le |V| - 1∣Ei​∣≤∣V∣−1 for all iii;
  2. if GGG is simple, ∣Ei∣≤∣V∣−i|E_i| \le |V| - i∣Ei​∣≤∣V∣−i for i≤∣V∣−1i \le |V| - 1i≤∣V∣−1 and Ei=∅E_i = \emptysetEi​=∅ for i≥∣V∣i \ge |V|i≥∣V∣.

Milestones

  • Lemma 2.2 (p. 587): during the execution, a node vvv has incident edges in exactly the classes E1,…,Er(v)E_1, \dots, E_{r(v)}E1​,…,Er(v)​.
  • Lemma 2.3 (p. 587): each (V,Ei)(V, E_i)(V,Ei​) is a forest at every instant.
  • Lemma 2.4 (p. 588): (a) an edge (u,v)(u, v)(u,v) added to EiE_iEi​ has its end nodes joined by a path in Ei−1E_{i-1}Ei−1​; (b) a path in EjE_jEj​ between uuu and vvv yields a path in every EiE_iEi​, i<ji < ji<j.
  • Lemma 2.5 (p. 588): each output (V,Ei)(V, E_i)(V,Ei​) is a maximal spanning forest in G−E1∪⋯∪Ei−1G - E_1 \cup \cdots \cup E_{i-1}G−E1​∪⋯∪Ei−1​.
  • Lemma 2.1 (p. 584): any sequence of successive maximal spanning forests satisfies (2.1).

Companions

  • the sparse certificate (p. 589): if λ(x,y;G)≥k\lambda(x, y; G) \ge kλ(x,y;G)≥k for all x,yx, yx,y, then GkG_kGk​ is kkk-edge-connected with ∣E(Gk)∣≤k(∣V∣−1)|E(G_k)| \le k(|V| - 1)∣E(Gk​)∣≤k(∣V∣−1), and ∣E(Gk)∣≤k∣V∣−k(k+1)/2|E(G_k)| \le k|V| - k(k+1)/2∣E(Gk​)∣≤k∣V∣−k(k+1)/2 for simple GGG;
  • Lemma 2.6 (p. 589): for k≤δ(G)k \le \delta(G)k≤δ(G), GkG_kGk​ has a node of degree exactly kkk.

Significance

(2.1) says that one search produces, for every threshold kkk at once, a subgraph with at most k(∣V∣−1)k(|V|-1)k(∣V∣−1) edges that keeps every local edge-connectivity up to kkk. Any algorithm whose running time grows with ∣E∣|E|∣E∣ can then be run on GkG_kGk​ in place of GGG; §4 of the paper uses this to speed up kkk-connectivity tests and the computation of local connectivities. The same forest partition is the structural fact behind the Nagamochi–Ibaraki and Stoer–Wagner minimum-cut algorithms. Lemma 2.6 shows that the certificate is tight: its edge-connectivity is exactly kkk when λ(G)≥k\lambda(G) \ge kλ(G)≥k.

The results are proved in the paper. No machine-checked proof of them is known. Formalizing them means formalizing a graph search with free tie-breaking as a transition system, reasoning about invariants of all its executions, and proving a cut-counting statement for multigraphs. Mathlib's connectivity notions, such as SimpleGraph.IsEdgeReachable, do not see parallel edges, so the multigraph cut theory here is new.

Difficulty

Lemma 2.1 is a short cut argument once maximality is available. The difficulty is showing that FOREST, which assigns each edge to a class by looking only at the label of one end node, produces maximal forests in the successive residual graphs (Lemma 2.5). The obvious invariant, that the class of an edge is the first forest it does not close a cycle in, is not what line 7 computes. The label r(y)r(y)r(y) records only which classes touch yyy, not which component of each class contains yyy. The paper's argument needs Lemma 2.4(a): at the moment an edge is added to EiE_iEi​ its ends already lie in one tree of Ei−1E_{i-1}Ei−1​. That relies on the choice of the unscanned node of largest label, on the order of lines 8 and 9, and on an argument about the scan order of tree roots.

Formalization scope

  • Graphs. A graph is a finite node type V with ∣V∣≥2|V| \ge 2∣V∣≥2, a finite edge type E, and ends : E → Sym2 V with no diagonal value (no self-loops). Parallel edges are distinct elements of E. Simplicity is injectivity of ends, and edge subsets are Finset E. A forest is an edge set in which every edge is a bridge, a condition that sees parallel edges.
  • Connectivity. λ\lambdaλ is an infimum in ℕ∞ over separating edge sets, ∞\infty∞ at x=yx = yx=y. (2.1) is kept "for all x,yx, yx,y", as printed.
  • FOREST. FOREST is a nondeterministic step relation (select, scan, finish) on explicit states: labels, class index per edge (000 = unscanned), scanned nodes, current node and selection order. The theorems quantify over every run from the initial state, so no tie-breaking rule is fixed. "At some time instant" is a state of the run; "upon completion" is a run whose last state has every node scanned. Lemma 2.2 is stated at every state, not only after a scan block (the other steps change neither labels nor classes). Lemma 2.4(a) assumes i≥2i \ge 2i≥2, since E0E_0E0​ does not exist.
  • Exclusions. Theorem 2.1's clause "is found in O(∣V∣+∣E∣)O(|V| + |E|)O(∣V∣+∣E∣) time", the bucket implementation, and the time bound of the certificate are not formalized: the paper fixes no machine model. The goal consists of the structural conclusions only. The bound printed "if GGG is multiple" is stated for every loopless graph.
  • Non-triviality. The goal is about the classes of a run of FOREST. An arbitrary partition of EEE into maximal spanning forests is Lemma 2.1's hypothesis, not a formalization of Theorem 2.1. A statement in which the classes are unconstrained variables, or in which the run hypotheses cannot be met, would be trivial. A separate sanity file checks that a complete run on the triangle K3K_3K3​ exists and attains ∣E1∣=∣V∣−1|E_1| = |V| - 1∣E1​∣=∣V∣−1, ∣E2∣=∣V∣−2|E_2| = |V| - 2∣E2​∣=∣V∣−2.
  • Welcome contributions. Useful reusable infrastructure includes:
    • a cut and Menger layer for finite multigraphs;
    • forest and bridge lemmas for edge-indexed graphs;
    • invariant-style reasoning over runs.

Selected references

  • H. Nagamochi, T. Ibaraki, A linear-time algorithm for finding a sparse kkk-connected spanning subgraph of a kkk-connected graph, Algorithmica 7 (1992) 583–596. https://doi.org/10.1007/BF01758778
  • H. Nagamochi, T. Ibaraki, Computing edge-connectivity in multigraphs and capacitated graphs, SIAM J. Discrete Math. 5 (1992) 54–66. https://doi.org/10.1137/0405004
  • T. Nishizeki, S. Poljak, kkk-connectivity and decomposition of graphs into forests, Discrete Appl. Math. 55 (1994) 295–301. https://doi.org/10.1016/0166-218X(94)90014-0
  • J. Cheriyan, M.-Y. Kao, R. Thurimella, Scan-first search and sparse certificates: an improved parallel algorithm for kkk-vertex connectivity, SIAM J. Comput. 22 (1993) 157–174. https://doi.org/10.1137/0222013
  • A. Frank, T. Ibaraki, H. Nagamochi, On sparse subgraphs preserving connectivity properties, J. Graph Theory 17 (1993) 275–281. https://doi.org/10.1002/jgt.3190170302
  • M. Stoer, F. Wagner, A simple min-cut algorithm, J. ACM 44 (1997) 585–591. https://doi.org/10.1145/263867.263872
10 thms3 active usersReviewed
🏆Completed
Control TheoryDynamical SystemsLinear algebra·Captain: mikedeng1

On Controllability of Delayed Boolean Control Networks: Trajectory Controllability Avoiding Forbidden Trajectories Iff the Reduced Transition-Count Matrix Is IrreducibleResearch Paper

Motivation

Boolean networks model gene regulatory networks by giving each gene an on/off value and updating all values synchronously by logical rules (Kauffman, 1969). Adding external Boolean inputs gives Boolean control networks (BCNs), in which a designer (a drug, an intervention) chooses the inputs over time. The basic control-theoretic question is controllability: can the inputs drive the network from any configuration to any other? Cheng and Qi (Automatica 2009) answered it for BCNs using the semi-tensor product of matrices, which rewrites a BCN as a linear recursion on canonical basis vectors.

Biological regulation is not instantaneous: transcription and translation introduce delays, so the next value of a gene can depend on several past values. Lu, Zhong, Ho, Tang and Cao (SIAM J. Control Optim. 2016) study delayed BCNs, in which the update reads the last μ\muμ states, and give criteria for two kinds of controllability, including the case where some configurations are dangerous and must be avoided (in biology, states corresponding to disease). The mission formalizes their criteria.

Setting

Write D={0,1}\mathcal D=\{0,1\}D={0,1}. A state is x=(x1,…,xn)∈Dnx=(x_1,\dots,x_n)\in\mathcal D^nx=(x1​,…,xn​)∈Dn and an input value is u∈Dmu\in\mathcal D^mu∈Dm. Fix μ≥1\mu\ge1μ≥1. The delayed BCN (2.2) is

xi(t+1)=fi(u(t),x(t−μ+1),…,x(t)),i=1,…,n,x_i(t+1)=f_i\big(u(t),x(t-\mu+1),\dots,x(t)\big),\qquad i=1,\dots,n,xi​(t+1)=fi​(u(t),x(t−μ+1),…,x(t)),i=1,…,n,

with arbitrary Boolean functions fi:Dm+μn→Df_i:\mathcal D^{m+\mu n}\to\mathcal Dfi​:Dm+μn→D, collected into one map FFF with x(t+1)=F(u(t),X(t))x(t+1)=F(u(t),X(t))x(t+1)=F(u(t),X(t)). A trajectory is the window X(t)=(x(t−μ+1),…,x(t))X(t)=(x(t-\mu+1),\dots,x(t))X(t)=(x(t−μ+1),…,x(t)) of the last μ\muμ states; its last entry is the current state. One step maps X(t)X(t)X(t) under the input u(t)u(t)u(t) to X(t+1)=(x(t−μ+2),…,x(t+1))X(t+1)=(x(t-\mu+2),\dots,x(t+1))X(t+1)=(x(t−μ+2),…,x(t+1)). A control sequence of length kkk is U=(u(0),…,u(k−1))U=(u(0),\dots,u(k-1))U=(u(0),…,u(k−1)), with inputs chosen freely; y(i)=X(i)y(i)=X(i)y(i)=X(i) denotes the trajectory after iii steps from an initial trajectory y(0)y(0)y(0).

The transition-count matrix QQQ has rows and columns indexed by trajectories: Qb,aQ_{b,a}Qb,a​ is the number of input values uuu taking trajectory aaa to trajectory bbb in one step. For a set CtC_tCt​ of forbidden trajectories, QCtQ_{C_t}QCt​​ is QQQ with the rows and columns of CtC_tCt​ replaced by zeros, and QCt\mathbb Q_{C_t}QCt​​ is QQQ with those rows and columns deleted.

The notions of controllability are:

  • Trajectory controllable (Definition 3.1): from every initial trajectory, every trajectory XdX_dXd​ equals X(k)X(k)X(k) for some k≥1k\ge1k≥1 and some control sequence.
  • Trajectory controllable under CtC_tCt​ (Definition 3.11): for all trajectories a,b∉Cta,b\notin C_ta,b∈/Ct​ there are k≥0k\ge0k≥0 and a control sequence steering y(0)=ay(0)=ay(0)=a to y(k)=by(k)=by(k)=b with y(i)∉Cty(i)\notin C_ty(i)∈/Ct​ for i=0,…,ki=0,\dots,ki=0,…,k.
  • State controllable (Definition 4.1): from every initial trajectory, every state equals x(k)x(k)x(k) for some k>0k>0k>0 and some control sequence.

A real square matrix AAA of size N≥2N\ge2N≥2 is reducible (Definition 3.7) if a simultaneous permutation of rows and columns brings it to block upper-triangular form (A11A120A22)\begin{pmatrix}A_{11}&A_{12}\\0&A_{22}\end{pmatrix}(A11​0​A12​A22​​) with square diagonal blocks; it is irreducible otherwise, so every 1×11\times11×1 matrix is irreducible.

N1(k;ya,yb,Ct)\mathbb N_1(k;y_a,y_b,C_t)N1​(k;ya​,yb​,Ct​) counts the control sequences of length kkk steering yay_aya​ to yby_byb​ while avoiding CtC_tCt​; N2(k;a,bs)\mathbb N_2(k;a,b_s)N2​(k;a,bs​) counts those steering the initial trajectory aaa to x(k)=bsx(k)=b_sx(k)=bs​; Ξμp\Xi^{p}_\muΞμp​ is the set of trajectories with current state ppp.

Formalization targets

Goal: Theorem 3.12

the delayed BCN is trajectory controllable under Ct  ⟺  QCt is irreducible,\text{the delayed BCN is trajectory controllable under } C_t\iff \mathbb Q_{C_t}\ \text{is irreducible},the delayed BCN is trajectory controllable under Ct​⟺QCt​​ is irreducible,

for every μ≥1\mu\ge1μ≥1, nnn, mmm, every update map and every forbidden set CtC_tCt​. It is the most general criterion of the paper.

Milestones

  • Proposition 3.5: for k>0k>0k>0, N1(k;ya,yb,Ct)=ybT(QCt)kya\mathbb N_1(k;y_a,y_b,C_t)=y_b^{\mathsf T}(Q_{C_t})^ky_aN1​(k;ya​,yb​,Ct​)=ybT​(QCt​​)kya​.
  • Remark 3: N1(k;ya,yb)=ybTQkya\mathbb N_1(k;y_a,y_b)=y_b^{\mathsf T}Q^ky_aN1​(k;ya​,yb​)=ybT​Qkya​.
  • Theorem 3.10: trajectory controllable   ⟺  \iff⟺ QQQ irreducible.
  • Theorem 4.3: N2(k;a,bs)=∑b∈ΞμbsbTQka\mathbb N_2(k;a,b_s)=\sum_{b\in\Xi^{b_s}_\mu}b^{\mathsf T}Q^kaN2​(k;a,bs​)=∑b∈Ξμbs​​​bTQka.
  • Theorem 5.1: the same count while avoiding forbidden states CsC_sCs​, with QQQ zeroed on the trajectories containing a state of CsC_sCs​.
  • Corollary 5.3: trajectory controllability implies state controllability.

A supporting item, Lemma 3.9 (corrected form), relates Definition 3.7 to positivity of entries of powers for nonnegative matrices of size at least 2.

Significance

The criteria replace a question about control sequences of unbounded length by a finite test on one nonnegative integer matrix of size 2μn2^{\mu n}2μn. The counting identities give more than a yes/no answer: they enumerate the control sequences achieving a transfer, which is what a designer choosing among interventions needs, and they handle forbidden states by passing to forbidden trajectories.

The results are proved in the paper; none of them has a machine-checked proof. Formalizing them yields a reusable Lean model of delayed (and, for μ=1\mu=1μ=1, ordinary) Boolean control networks defined directly through their dynamics, together with a checked bridge from dynamic reachability to the combinatorics of nonnegative matrices. A formal version also pins down points the printed text leaves loose: Lemma 3.9 as printed is false (it characterizes primitive, not irreducible, matrices), and the two different matrices QCtQ_{C_t}QCt​​ (zeroed) and QCt\mathbb Q_{C_t}QCt​​ (deleted) must not be confused.

Difficulty

The central step is the passage between three objects: reachability under the dynamics, positivity of entries of powers of QQQ, and the block-triangular definition of irreducibility. The counting identity for powers of QQQ is a statement about sequences of inputs, not about paths in a graph, so the bijection between control sequences and weighted walks has to be established with the avoidance constraint carried at every intermediate time, including the endpoints. The equivalence between Definition 3.7 and strong connectivity is the classical graph-theoretic characterization, but it has edge cases the obvious argument misses: a 1×11\times11×1 zero matrix is irreducible by Definition 3.7 yet has no positive power, which is exactly why Definition 3.11 allows k=0k=0k=0 while Definition 3.1 does not. A proof that routes through "some power of QCt\mathbb Q_{C_t}QCt​​ is entrywise positive", following the printed Lemma 3.9, proves a false intermediate statement.

Formalization scope

Conventions committed to in Lean:

  • A state is Fin n → Bool (true is the paper's 1); a trajectory is Fin μ → State n with index 0 the oldest state and index μ−1\mu-1μ−1 the current state; μ≥1\mu\ge1μ≥1 is a standing assumption (NeZero μ); n,m≥0n,m\ge0n,m≥0 are arbitrary and no assumption is placed on the update functions.
  • Matrices are indexed by trajectories rather than by the paper's index jjj of δ2μnj\delta^j_{2^{\mu n}}δ2μnj​. The paper's Q=L⋉12mQ=L\ltimes\mathbf 1_{2^m}Q=L⋉12m​ is obtained by the simultaneous relabelling of Lemma 2.6, and every statement (entries of powers, sums over sets of trajectories, Definition 3.7) is invariant under it. QCt\mathbb Q_{C_t}QCt​​ is indexed by the subtype of allowed trajectories.
  • Irreducibility is Definition 3.7 applied to the matrix with entries cast to R\mathbb RR; it is not Mathlib's Matrix.IsIrreducible, which differs on 1×11\times11×1 matrices.
  • Pinned readings: Definition 3.1 and Definition 4.1 use k≥1k\ge1k≥1; Definition 3.11 uses k≥0k\ge0k≥0; Proposition 3.5, Theorem 4.3 and Theorem 5.1 are stated for k>0k>0k>0 (Proposition 3.5 is false at k=0k=0k=0 when ya=yb∈Cty_a=y_b\in C_tya​=yb​∈Ct​); Remark 3 holds for all k≥0k\ge0k≥0. "Avoiding CsC_sCs​" in Theorem 5.1 means that no state x(i)x(i)x(i), i=1−μ,…,ki=1-\mu,\dots,ki=1−μ,…,k, lies in CsC_sCs​, which is the theorem's own middle term N1(k;a,b,ΞCs)\mathbb N_1(k;a,b,\Xi^{C_s})N1​(k;a,b,ΞCs​). Ξμp\Xi^p_\muΞμp​ is defined semantically by eq. (4.3), since the printed index range in (4.4) is a misprint. Lemma 3.9 is included in corrected form with N≥2N\ge2N≥2.

Controllability is defined through the dynamics and control sequences, never as positivity of entries of powers of QQQ; a formalization that defines reachability by (Qk)b,a>0(Q^k)_{b,a}>0(Qk)b,a​>0 would reduce Theorems 3.10 and 3.12 to library facts and is ruled out.

Needed infrastructure: the correspondence between control sequences and products of entries of QQQ (the core of Proposition 3.5), and the equivalence of Definition 3.7 with strong connectivity of the support graph for nonnegative matrices (Mathlib's Matrix.IsIrreducible, Matrix.isIrreducible_iff_exists_pow_pos and Matrix.pow_apply_pos_iff_nonempty_path cover much of the second part for sizes at least 2). Both are reusable beyond this mission, for ordinary BCNs, probabilistic BCNs and finite automata. Contributions of either piece, or of the semi-tensor-product bridge identifying QQQ with L⋉12mL\ltimes\mathbf 1_{2^m}L⋉12m​, are welcome.

Selected references

  • J. Lu, J. Zhong, D. W. C. Ho, Y. Tang, J. Cao, On Controllability of Delayed Boolean Control Networks, SIAM J. Control Optim. 54(2):475–494, 2016. https://doi.org/10.1137/140991820
  • D. Cheng, H. Qi, Controllability and observability of Boolean control networks, Automatica 45(7):1659–1667, 2009. https://doi.org/10.1016/j.automatica.2009.03.006
  • S. A. Kauffman, Metabolic stability and epigenesis in randomly constructed genetic nets, J. Theoret. Biol. 22(3):437–467, 1969. https://doi.org/10.1016/0022-5193(69)90015-0
  • A. Berman, R. J. Plemmons, Nonnegative Matrices in the Mathematical Sciences, SIAM Classics in Applied Mathematics 9, 1994. https://doi.org/10.1137/1.9781611971262
9 thms3 active usersReviewed
CombinatoricsLinear OptimizationOperations Research·Captain: mikedeng1

Cones of Matrices and Set-Functions and 0–1 Optimization II: One Round of N on the Stable Set Polytope Gives Exactly the Odd Hole ConstraintsResearch Paper

Motivation

The stable set problem (vertex packing) asks for a largest set of pairwise non-adjacent nodes of a graph. It is NP-hard, and its polyhedral study, the description of the stable set polytope STAB(G)\mathrm{STAB}(G)STAB(G) by linear inequalities, is one of the most studied topics of polyhedral combinatorics. Classes of valid inequalities (clique, odd hole, odd antihole, wheel constraints) and the graph classes they describe exactly (perfect, ttt-perfect, hhh-perfect graphs) organize much of that literature; see Grötschel, Lovász and Schrijver, Geometric Algorithms and Combinatorial Optimization (Springer, 1988).

Lovász and Schrijver (SIAM J. Optim. 1(2), 1991) introduced a general lift-and-project procedure for 0–1 programs: lift a relaxation KKK into a space of matrices, impose linear conditions that every 0–1 point satisfies, and project back. One round of their operator NNN gives a tighter relaxation N(K)N(K)N(K) that still contains every 0–1 point of KKK; nnn rounds give the 0–1 hull. The procedure is an ancestor of the Sherali–Adams and Lasserre hierarchies, and the stable set problem is its first test case. This mission formalizes the paper's exact description of what one round of NNN does to the fractional stable set polytope: it adds precisely the odd hole constraints.

Setting

Let G=(V,E)G = (V, E)G=(V,E) be a finite graph with no isolated nodes, n=∣V∣n = |V|n=∣V∣. Vectors of RV∪{0}\mathbb{R}^{V \cup \{0\}}RV∪{0} have a distinguished coordinate x0x_0x0​; RV\mathbb{R}^VRV sits inside as the hyperplane H0={x0=1}H_0 = \{x_0 = 1\}H0​={x0​=1}, via x↦(1,x)x \mapsto (1, x)x↦(1,x).

  • FRAC(G)⊆RV\mathrm{FRAC}(G) \subseteq \mathbb{R}^VFRAC(G)⊆RV is the solution set of the nonnegativity constraints xi≥0x_i \ge 0xi​≥0 (i∈Vi \in Vi∈V) and the edge constraints xi+xj≤1x_i + x_j \le 1xi​+xj​≤1 (ij∈Eij \in Eij∈E).
  • FR(G)⊆RV∪{0}\mathrm{FR}(G) \subseteq \mathbb{R}^{V\cup\{0\}}FR(G)⊆RV∪{0} is the cone given by xi≥0x_i \ge 0xi​≥0 and xi+xj≤x0x_i + x_j \le x_0xi​+xj​≤x0​; it is the cone spanned by the vectors (1,x)(1, x)(1,x) with x∈FRAC(G)x \in \mathrm{FRAC}(G)x∈FRAC(G).
  • QQQ is the cone spanned by the 0–1 vectors with x0=1x_0 = 1x0​=1. For a convex cone KKK, its polar cone is K∗={u:uTx≥0 ∀x∈K}K^* = \{u : u^{\mathsf T}x \ge 0 \ \forall x \in K\}K∗={u:uTx≥0 ∀x∈K}.
  • M(K)=M(K,Q)M(K) = M(K, Q)M(K)=M(K,Q) is the set of (n+1)×(n+1)(n+1)\times(n+1)(n+1)×(n+1) matrices Y=(yij)Y = (y_{ij})Y=(yij​) that are symmetric, satisfy yii=y0iy_{ii} = y_{0i}yii​=y0i​ for i∈Vi \in Vi∈V, and satisfy uTYv≥0u^{\mathsf T} Y v \ge 0uTYv≥0 for all u∈K∗u \in K^*u∈K∗, v∈Q∗v \in Q^*v∈Q∗.
  • N(K)={Ye0:Y∈M(K)}N(K) = \{Y e_0 : Y \in M(K)\}N(K)={Ye0​:Y∈M(K)}, and N(G)={x∈RV:(1,x)∈N(FR(G))}N(G) = \{x \in \mathbb{R}^V : (1, x) \in N(\mathrm{FR}(G))\}N(G)={x∈RV:(1,x)∈N(FR(G))}.
  • A set C⊆VC \subseteq VC⊆V is an odd hole if it induces a chordless cycle of odd length ∣C∣≥3|C| \ge 3∣C∣≥3 (triangles included). Its odd hole constraint is ∑i∈Cxi≤12(∣C∣−1)\sum_{i \in C} x_i \le \frac12(|C| - 1)∑i∈C​xi​≤21​(∣C∣−1).

Formalization targets

Goal: Theorem 2.3 (p. 178)

For every finite graph GGG without isolated nodes,

N(G)={x∈RV:xi≥0 (i∈V),  xi+xj≤1 (ij∈E),  ∑i∈Cxi≤12(∣C∣−1) (C an odd hole)}.N(G) = \Big\{x \in \mathbb{R}^V : x_i \ge 0\ (i \in V),\ \ x_i + x_j \le 1\ (ij \in E),\ \ \sum_{i \in C} x_i \le \tfrac12(|C|-1)\ (C \text{ an odd hole})\Big\}.N(G)={x∈RV:xi​≥0 (i∈V),  xi​+xj​≤1 (ij∈E),  i∈C∑​xi​≤21​(∣C∣−1) (C an odd hole)}.

Milestones, in the order the proof uses them

  1. Lemma 1.3 (p. 171): for a convex cone K⊆QK \subseteq QK⊆Q and i∈Vi \in Vi∈V, N(K)⊆(K∩Hi)+(K∩Gi)N(K) \subseteq (K \cap H_i) + (K \cap G_i)N(K)⊆(K∩Hi​)+(K∩Gi​), with Hi={xi=0}H_i = \{x_i = 0\}Hi​={xi​=0}, Gi={xi=x0}G_i = \{x_i = x_0\}Gi​={xi​=x0​}.
  2. Lemma 2.2 (p. 178): if both the deletion and the contraction of some node vvv give inequalities valid for KKK, then aTx≤ba^{\mathsf T}x \le baTx≤b is valid for N(K)N(K)N(K).
  3. Part (1) of the proof of Theorem 2.3 (p. 178): for an odd hole CCC and i∈Ci \in Ci∈C, the deletion and contraction of iii in the odd hole constraint are valid for FRAC(G)\mathrm{FRAC}(G)FRAC(G).
  4. Observation of Section 2.b (p. 177): every Y∈M(FR(G))Y \in M(\mathrm{FR}(G))Y∈M(FR(G)) has yij=0y_{ij} = 0yij​=0 for ij∈Eij \in Eij∈E.
  5. Part (2) of the proof of Theorem 2.3 (p. 178): x∈N(G)x \in N(G)x∈N(G) if and only if some nonnegative symmetric YYY with y00=1y_{00} = 1y00​=1, yi0=yii=xiy_{i0} = y_{ii} = x_iyi0​=yii​=xi​ satisfies xi+xj+xk−1≤yik+yjk≤xkx_i + x_j + x_k - 1 \le y_{ik} + y_{jk} \le x_kxi​+xj​+xk​−1≤yik​+yjk​≤xk​ for all i,j,ki, j, ki,j,k with ij∈Eij \in Eij∈E.
  6. Lemma 2.4 (p. 178): a system a(ij)≤yi+yj≤b(ij)a(ij) \le y_i + y_j \le b(ij)a(ij)≤yi​+yj​≤b(ij), y≥0y \ge 0y≥0, y∣U=0y|_U = 0y∣U​=0 on a graph is infeasible if and only if a walk with a negative alternating sum of one of four types exists.

Significance

Theorem 2.3 gives a complete description of one round of NNN on the stable set problem: the only new constraints are the odd hole constraints. Consequences:

  • For ttt-perfect graphs (those for which nonnegativity, edge and odd hole constraints describe STAB(G)\mathrm{STAB}(G)STAB(G)), N(G)=STAB(G)N(G) = \mathrm{STAB}(G)N(G)=STAB(G).
  • It is the base case for the paper's bounds on the NNN-index of stable set inequalities (Theorem 2.13), and it contrasts with the semidefinite operator N+N_+N+​, which after one round already satisfies clique, odd antihole and wheel constraints.
  • Lemma 2.4 is a combinatorial feasibility criterion for systems with two variables per inequality, useful beyond this paper.

The result has been proved since 1991. At the time of drafting, Prove2Me holds no formalization of it or of any part of the Lovász–Schrijver construction, and Mathlib has none. The mission produces a formal account of the NNN operator on the stable set polytope and a formal proof of the walk criterion for two-variable systems.

Difficulty

The inclusion of N(G)N(G)N(G) in the odd hole system is a short argument once Lemma 1.3 is available. The reverse inclusion is the substance: given xxx satisfying all odd hole constraints, one must exhibit a lifted matrix YYY. A direct appeal to Farkas' lemma yields a certificate with no visible relation to odd cycles; the difficulty is to show that every obstruction to solvability of the matrix system forces a violated odd hole constraint, which is what Lemma 2.4 and the analysis of its four walk types accomplish. Case (d) of that analysis needs the odd hole constraints; the other cases need only the edge constraints. Lemma 2.4 itself is called folklore on the page and is stated without proof there.

A further point: Lemma 2.4 is stated for lower bounds 0≤a0 \le a0≤a, while the lower bounds that arise from the matrix system, xi+xj+xk−1x_i + x_j + x_k - 1xi​+xj​+xk​−1, can be negative.

Formalization scope

  • Coordinates of RV∪{0}\mathbb{R}^{V\cup\{0\}}RV∪{0} are indexed by Option V, with none the coordinate x0x_0x0​. Graphs are Mathlib SimpleGraphs on a finite type VVV with decidable adjacency. Every statement about a graph carries the paper's standing assumption that GGG has no isolated nodes (∀ v, ∃ w, G.Adj v w).
  • MMM is defined by condition (iii), never by its rewritings. Lemma 1.3 and Lemma 2.2 take the cone KKK closed, a hypothesis the paper leaves tacit (its cones are polyhedral); for a non-closed KKK Lemma 1.3 is false. FR(G)\mathrm{FR}(G)FR(G) is polyhedral, so the goal needs no such hypothesis.
  • FR(G)\mathrm{FR}(G)FR(G) is defined by its constraints; this agrees with the cone over FRAC(G)\mathrm{FRAC}(G)FRAC(G) because GGG has no isolated nodes.
  • Lemma 2.2 is stated in cone form: KKK is any closed convex cone inside FR(G)\mathrm{FR}(G)FR(G), and validity is read on the slice x0=1x_0 = 1x0​=1. The paper's extra hypothesis STAB(G)⊆K\mathrm{STAB}(G) \subseteq KSTAB(G)⊆K is dropped, which strengthens the lemma.
  • Deletion and contraction of a node are coefficient vectors on the same graph (coefficients set to 000), not inequalities on the subgraphs G−vG - vG−v and G−Γ(v)−vG - \Gamma(v) - vG−Γ(v)−v.
  • Odd holes are chordless odd cycles including triangles; triangles are needed, as 121\tfrac12\mathbf 121​1 satisfies all other constraints on a triangle.
  • The matrix system of part (2) is stated as an equivalence; the page uses one direction.
  • Lemma 2.4 uses edge values on unordered pairs and strict inequalities, exactly as printed.

A trivializing formalization is ruled out: the goal is the set equality for every graph without isolated nodes, not the existence of a lifted matrix and not a single graph.

Not formalized here: the semidefinite operator N+N_+N+​, the operator N^\hat NN^, algorithmic statements (Theorems 1.6, 2.1, Corollary 2.5), and the set-function results of Section 3.

Reusable beyond this mission: the matrix cone layer (QQQ, MMM, NNN), the stable-set cones, and the two-variable feasibility criterion of Lemma 2.4. Contributions of any of the milestones, and of general facts about polar cones of polyhedral cones in this setting, are welcome.

Selected references

  • L. Lovász and A. Schrijver, Cones of matrices and set-functions and 0–1 optimization, SIAM Journal on Optimization 1(2) (1991) 166–190. https://doi.org/10.1137/0801013
  • M. Grötschel, L. Lovász and A. Schrijver, Geometric Algorithms and Combinatorial Optimization, Springer, 1988. https://doi.org/10.1007/978-3-642-97881-4
  • H. D. Sherali and W. P. Adams, A hierarchy of relaxations between the continuous and convex hull representations for zero-one programming problems, SIAM Journal on Discrete Mathematics 3(3) (1990) 411–430. https://doi.org/10.1137/0403036
13 thms3 active usersReviewed
🏆Completed
CombinatoricsNumber TheoryTheoretical Computer Science·Captain: mikedeng1

Fast Algorithms for Finding Nearest Common Ancestors II: Nearest Common Ancestors in a Complete Binary Tree by Symmetric-Order ArithmeticResearch Paper

Motivation

The nearest common ancestor (nca) problem asks, for a fixed rooted tree and a sequence of vertex pairs (v,w)(v, w)(v,w), for the deepest vertex that is an ancestor of both. It is a basic step in suffix-tree string algorithms and is equivalent to range-minimum queries (Bender, Farach-Colton, 2000). Harel and Tarjan, Fast Algorithms for Finding Nearest Common Ancestors, SIAM J. Comput. 13 (1984) 338–355, gave the first algorithm answering each query on a static tree in constant time on a random-access machine after linear preprocessing.

Their construction reduces the general problem to the case of a complete binary tree, where §3 of the paper shows that nca queries can be answered "by direct calculation" on vertex numbers: multiplication, division, powers of two, the base-two logarithm and bitwise exclusive or. The later simplification of Schieber and Vishkin (1988) is built on the same in-order numbering of a complete binary tree. This mission formalizes that arithmetic core.

Timeline, as reviewed in the paper's §1 (pp. 338–340):

  • 1976: Aho, Hopcroft and Ullman (SIAM J. Comput. 5) give an O(n+mα(m+n,n))O(n + m\alpha(m+n, n))O(n+mα(m+n,n))-time off-line algorithm on a pointer machine, and for static trees a random-access algorithm with O(nlog⁡log⁡n)O(n \log\log n)O(nloglogn) preprocessing and O(log⁡log⁡n)O(\log\log n)O(loglogn) time per query.
  • 1976: van Leeuwen (unpublished report) gives an O(n+mlog⁡log⁡n)O(n + m \log\log n)O(n+mloglogn)-time algorithm for linking roots and static trees that runs on a pointer machine in O(n)O(n)O(n) space.
  • 1980: Harel (Proc. 21st FOCS) gives a preliminary version of the paper's results.
  • 1984: Harel and Tarjan prove that pointer machines need Ω(log⁡log⁡n)\Omega(\log\log n)Ω(loglogn) time per query on static trees (Theorem 1), and give the O(n)O(n)O(n)-preprocessing, O(1)O(1)O(1)-query random-access algorithm whose base case is the subject of this mission.

Setting

Fix d≥0d \ge 0d≥0 and let TTT be the complete binary tree of depth ddd. A vertex is identified with the path from the root to it, a word of at most ddd left or right turns; the root is the empty word and TTT has n=2d+1−1n = 2^{d+1} - 1n=2d+1−1 vertices. Following the paper's Appendix (pp. 354–355):

  • www is an ancestor of vvv (vvv a descendant of www) if the word www is a prefix of the word vvv; every vertex is its own ancestor. vvv and www are unrelated if neither is an ancestor of the other.
  • The depth of vvv is its distance to the root; its height h(v)h(v)h(v) is the length of the longest path from a leaf to vvv, which in TTT is d−depth⁡(v)d - \operatorname{depth}(v)d−depth(v).
  • nca⁡(v,w)\operatorname{nca}(v, w)nca(v,w) is the vertex of greatest depth that is an ancestor of both: the longest common prefix.

The vertices of TTT are numbered from 111 to nnn in symmetric order (in-order): at every vertex, first the left subtree, then the vertex, then the right subtree. sym(v)\mathrm{sym}(v)sym(v) is the number of vvv and sym−1(i)\mathrm{sym}^{-1}(i)sym−1(i) the vertex numbered iii. For d=4d = 4d=4 (Fig. 1 of the paper) the root is 161616, its children 888 and 242424, and the leaves 1,3,5,…,311, 3, 5, \dots, 311,3,5,…,31. i⊕ji \oplus ji⊕j denotes bitwise exclusive or and lg⁡\lglg the base-two logarithm.

Two procedures of §3 use only numbers, heights and ddd:

  • the nca depth algorithm: return d−h(v)d - h(v)d−h(v) if sym(w)∈[sym(v)−2h(v)+1,sym(v)+2h(v)−1]\mathrm{sym}(w) \in [\mathrm{sym}(v) - 2^{h(v)} + 1, \mathrm{sym}(v) + 2^{h(v)} - 1]sym(w)∈[sym(v)−2h(v)+1,sym(v)+2h(v)−1]; else d−h(w)d - h(w)d−h(w) if the same holds with v,wv, wv,w exchanged; else d−⌊lg⁡(sym(v)⊕sym(w))⌋d - \lfloor \lg(\mathrm{sym}(v) \oplus \mathrm{sym}(w)) \rfloord−⌊lg(sym(v)⊕sym(w))⌋;
  • the depth algorithm: given vvv and a depth d2≤depth⁡(v)d_2 \le \operatorname{depth}(v)d2​≤depth(v), with h=d−d2h = d - d_2h=d−d2​, return sym−1(2h+1⌊sym(v)/2h+1⌋+2h)\mathrm{sym}^{-1}\bigl(2^{h+1}\lfloor \mathrm{sym}(v)/2^{h+1}\rfloor + 2^h\bigr)sym−1(2h+1⌊sym(v)/2h+1⌋+2h).

Formalization targets

Goal: the nca algorithm is correct

The algorithm to compute nca⁡(v,w)\operatorname{nca}(v,w)nca(v,w) (p. 342) runs the nca depth algorithm to obtain d0d_0d0​ and then the depth algorithm on (v,d0)(v, d_0)(v,d0​). The goal states that it returns the nearest common ancestor: for all vertices v,wv, wv,w of TTT, with d0d_0d0​ the output of the nca depth algorithm and h=d−d0h = d - d_0h=d−d0​,

sym(nca⁡(v,w))=2h+1⌊sym(v)2h+1⌋+2h.\mathrm{sym}(\operatorname{nca}(v,w)) = 2^{h+1}\left\lfloor \frac{\mathrm{sym}(v)}{2^{h+1}} \right\rfloor + 2^h .sym(nca(v,w))=2h+1⌊2h+1sym(v)​⌋+2h.

Milestones

In the order the paper uses them:

  1. Numbers at height hhh (p. 341): the vertices of height hhh are numbered 2h,3⋅2h,5⋅2h,…2^h, 3\cdot 2^h, 5\cdot 2^h, \dots2h,3⋅2h,5⋅2h,… from left to right.
  2. Lemma 1: h(v)h(v)h(v) is the largest hhh with 2h∣sym(v)2^h \mid \mathrm{sym}(v)2h∣sym(v).
  3. Lemma 2: the descendants of vvv are the vertices numbered in [sym(v)−2h(v)+1,sym(v)+2h(v)−1][\mathrm{sym}(v) - 2^{h(v)} + 1, \mathrm{sym}(v) + 2^{h(v)} - 1][sym(v)−2h(v)+1,sym(v)+2h(v)−1].
  4. Lemma 3: for a height h≥h(v)h \ge h(v)h≥h(v), the height-hhh ancestor of vvv has number 2h+1⌊sym(v)/2h+1⌋+2h2^{h+1}\lfloor \mathrm{sym}(v)/2^{h+1}\rfloor + 2^h2h+1⌊sym(v)/2h+1⌋+2h.
  5. Lemma 4: for unrelated v,wv, wv,w,
h(nca⁡(v,w))=⌊lg⁡(sym(v)⊕sym(w))⌋.h(\operatorname{nca}(v,w)) = \lfloor \lg(\mathrm{sym}(v) \oplus \mathrm{sym}(w)) \rfloor .h(nca(v,w))=⌊lg(sym(v)⊕sym(w))⌋.
  1. The nca depth algorithm returns depth⁡(nca⁡(v,w))\operatorname{depth}(\operatorname{nca}(v,w))depth(nca(v,w)).
  2. The depth algorithm returns the number of the depth-d2d_2d2​ ancestor of vvv.

Two supporting statements pin the definitions to the paper: sym\mathrm{sym}sym is a bijection onto {1,…,2d+1−1}\{1, \dots, 2^{d+1} - 1\}{1,…,2d+1−1}, and the longest common prefix is the deepest common ancestor.

Significance

The constant-time nca computation on complete binary trees is the base case of the whole paper: §§4–5 embed an arbitrary tree into a moderately sized complete binary tree through a compressed tree and a balanced binary tree, and every query ends with the arithmetic of §3. The same idea, that in-order numbers encode ancestry in their low-order bits, underlies the Schieber–Vishkin algorithm. Lemma 1 identifies the height with the 2-adic valuation of the number, and Lemma 4 identifies the nca height with the position of the highest differing bit.

The results are proved in the paper, with the proofs left as "easy to verify". No machine-checked version of this numbering or of these four lemmas is known to exist in Mathlib or on this platform. A formal development supplies proofs of the four lemmas and the two algorithms, and a reusable library connecting in-order ranks of a complete binary tree to binary arithmetic (Nat.log, bitwise xor, 2-adic valuation).

Difficulty

The numbering is defined by a traversal order, while the lemmas speak about divisibility, floor division and exclusive or. The work lies in connecting the rank of a vertex in symmetric order to its closed form (2j+1)⋅2h(v)(2j+1)\cdot 2^{h(v)}(2j+1)⋅2h(v), where jjj is its left-to-right position. That counting argument sums the sizes of the subtrees that precede vvv and is where most of the effort goes. Lemma 4 then needs the observation that two unrelated numbers agree in all bits above the height of their nca and differ in the bit at that height. This is a statement about Nat.testBit of the exclusive or, and it fails for related vertices. The algorithm statements add a case analysis whose first two cases overlap when v=wv = wv=w.

Formalization scope

  • A vertex of the tree of depth ddd is a List Bool of length at most ddd (false = left). Ancestry is the prefix relation, nca⁡\operatorname{nca}nca the longest common prefix, depth the length, and height d−lengthd - \text{length}d−length. None of these structural notions uses the numbering.
  • sym(v)\mathrm{sym}(v)sym(v) is the number of vertices whose in-order sort key is lexicographically at most that of vvv. The key is the path with left ↦0\mapsto 0↦0, right ↦2\mapsto 2↦2, followed by 111. The numbering is not defined by the closed form or by a recursion on numbers: a definition of that kind would make the height-hhh numbering and Lemma 1 immediate and move the content of the mission into an uncheckable definition.
  • ⌊lg⁡x⌋\lfloor \lg x \rfloor⌊lgx⌋ is Nat.log 2 x, which agrees for x≥1x \ge 1x≥1. ⊕\oplus⊕ is ^^^ on N\mathbb NN, and floor division is / on N\mathbb NN.
  • Interval tests a∈[b−c+1,b+c−1]a \in [b - c + 1, b + c - 1]a∈[b−c+1,b+c−1] are written additively as b+1≤a+cb + 1 \le a + cb+1≤a+c and a+1≤b+ca + 1 \le b + ca+1≤b+c. The subtractions d−h(v)d - h(v)d−h(v) and d−d2d - d_2d−d2​ never truncate for heights and depths of vertices.
  • Lemma 3 states explicitly that h≤dh \le dh≤d ("hhh is a height") and that the ancestor exists. The depth algorithm assumes d2≤depth⁡(v)d_2 \le \operatorname{depth}(v)d2​≤depth(v), as printed.
  • sym−1\mathrm{sym}^{-1}sym−1 is not defined as a function. The goal and the depth algorithm state that a vertex has the computed number if and only if it is the nearest common ancestor (respectively the ancestor at depth d2d_2d2​), which says that sym−1\mathrm{sym}^{-1}sym−1 of that number is that vertex.
  • The O(1)O(1)O(1) time bounds are not formalized, since the random-access machine model is out of scope.

Proofs of any milestone are welcome.

Selected references

  • D. Harel and R. E. Tarjan, Fast Algorithms for Finding Nearest Common Ancestors, SIAM J. Comput. 13(2) (1984), 338–355. https://doi.org/10.1137/0213024
  • A. V. Aho, J. E. Hopcroft and J. D. Ullman, On Finding Lowest Common Ancestors in Trees, SIAM J. Comput. 5(1) (1976), 115–132. https://doi.org/10.1137/0205011
  • B. Schieber and U. Vishkin, On Finding Lowest Common Ancestors: Simplification and Parallelization, SIAM J. Comput. 17(6) (1988), 1253–1262. https://doi.org/10.1137/0217079
  • M. A. Bender and M. Farach-Colton, The LCA Problem Revisited, LATIN 2000, LNCS 1776, 88–94. https://doi.org/10.1007/10719839_9
11 thms3 active usersReviewed
🏆Completed
Operations ResearchOptimizationTheoretical Computer Science·Captain: mikedeng1

A New Approach to the Maximum-Flow Problem 2: The Nonsaturating-Push Bound for FIFO Push-RelabelResearch Paper

Motivation

The maximum-flow problem asks how much of a commodity can be sent from a source to a sink through a network whose edges carry capacities. It is a basic model of operations research. Transportation, scheduling, bipartite matching and image segmentation reduce to it, and it is the inner step of many combinatorial algorithms.

Goldberg and Tarjan introduced the push–relabel (preflow) method in A New Approach to the Maximum-Flow Problem (J. ACM 35(4), 1988). Ford–Fulkerson-type algorithms augment along whole source–sink paths. The push–relabel method instead moves excess flow across single edges, guided by integer distance labels on the vertices. Whatever order its local operations are applied in, it is correct and performs O(n2m)O(n^2 m)O(n2m) of them (§3 of the paper). Section 4 shows that one particular order, processing the active vertices first-in, first-out, cuts the dominant term, the number of nonsaturating pushes, to O(n3)O(n^3)O(n3). The method and its FIFO and highest-label variants are the standard practical maximum-flow codes.

Timeline:

  • 1956: Ford and Fulkerson, augmenting paths and max-flow min-cut.
  • 1970–72: Dinic, and Edmonds and Karp, give polynomial augmenting-path bounds.
  • 1974: Karzanov introduces preflows and obtains O(n3)O(n^3)O(n3).
  • 1982: Shiloach and Vishkin give a parallel O(n2log⁡n)O(n^2 \log n)O(n2logn) preflow algorithm with a first-in, first-out flavour.
  • 1988: Goldberg and Tarjan, the generic push–relabel method, the FIFO bound of this mission, and O(nmlog⁡(n2/m))O(nm \log(n^2/m))O(nmlog(n2/m)) with dynamic trees.

Setting

A flow network has a finite vertex set VVV with n=∣V∣n = |V|n=∣V∣, a capacity c(v,w)≥0c(v,w) \ge 0c(v,w)≥0 on every ordered pair, a source sss and a sink t≠st \neq st=s. The edges are the pairs with c(v,w)>0c(v,w) > 0c(v,w)>0, and there are no loops. A preflow is a function fff on vertex pairs with f(v,w)≤c(v,w)f(v,w) \le c(v,w)f(v,w)≤c(v,w) and f(v,w)=−f(w,v)f(v,w) = -f(w,v)f(v,w)=−f(w,v). Its excess e(v)=∑uf(u,v)e(v) = \sum_u f(u,v)e(v)=∑u​f(u,v) must be nonnegative at every v≠sv \neq sv=s. The residual capacity is rf(v,w)=c(v,w)−f(v,w)r_f(v,w) = c(v,w) - f(v,w)rf​(v,w)=c(v,w)−f(v,w). A labeling ddd assigns each vertex a value in N∪{∞}\mathbb{N} \cup \{\infty\}N∪{∞}. A vertex v∉{s,t}v \notin \{s,t\}v∈/{s,t} is active if d(v)<∞d(v) < \inftyd(v)<∞ and e(v)>0e(v) > 0e(v)>0.

The two basic operations (Fig. 1 of the paper) are:

  • push(v,w)(v,w)(v,w), applicable when vvv is active, rf(v,w)>0r_f(v,w) > 0rf​(v,w)>0 and d(v)=d(w)+1d(v) = d(w)+1d(v)=d(w)+1. It sends δ=min⁡(e(v),rf(v,w))\delta = \min(e(v), r_f(v,w))δ=min(e(v),rf​(v,w)) from vvv to www. The push is saturating if rf(v,w)=0r_f(v,w) = 0rf​(v,w)=0 afterwards and nonsaturating otherwise.
  • relabel(v)(v)(v), applicable when vvv is active and d(v)≤d(w)d(v) \le d(w)d(v)≤d(w) for every residual edge (v,w)(v,w)(v,w). It sets d(v)←min⁡{d(w)+1:rf(v,w)>0}d(v) \leftarrow \min\{d(w)+1 : r_f(v,w) > 0\}d(v)←min{d(w)+1:rf​(v,w)>0}.

The algorithm starts by saturating every edge leaving sss, with d(s)=nd(s) = nd(s)=n and d(v)=0d(v) = 0d(v)=0 for v≠sv \neq sv=s.

In the first-in, first-out algorithm (§4), each vertex vvv scans a fixed list L(v)L(v)L(v) of its neighbours through a current edge. The push/relabel(v)(v)(v) operation pushes through the current edge if possible. Otherwise it advances the current edge, or, at the end of the list, returns to the first edge and relabels vvv. Active vertices wait in a queue QQQ, initially {v∈V−{s,t}:c(s,v)>0}\{v \in V - \{s,t\} : c(s,v) > 0\}{v∈V−{s,t}:c(s,v)>0}. The discharge operation removes the front vertex vvv and repeats push/relabel(v)(v)(v) until e(v)=0e(v) = 0e(v)=0 or d(v)d(v)d(v) increases. Every vertex that becomes active meanwhile is appended to QQQ, and vvv is appended too if it is still active. Passes over the queue are defined inductively. Pass 1 consists of the discharges of the initially queued vertices. Pass i+1i+1i+1 consists of the discharges of vertices added during pass iii.

Formalization targets

Goal: Corollary 4.4 (p. 931)

For every network, every edge-list order, every initial queue order, and every run of the FIFO algorithm,

#{nonsaturating pushes}≤4n3.\#\{\text{nonsaturating pushes}\} \le 4n^3 .#{nonsaturating pushes}≤4n3.

The constant is the printed one.

Milestones

  • Lemma 4.1 (p. 929): the push/relabel operation relabels only when relabeling is applicable.
  • Lemma 3.5 (p. 926): from any vertex with positive excess, the source is reachable in the residual graph.
  • Lemma 3.7 (p. 927): at any time, d(v)≤2n−1d(v) \le 2n-1d(v)≤2n−1 for every vertex.
  • Lemma 3.8 (p. 927): at most 2n−12n-12n−1 relabelings per vertex and at most (2n−1)(n−2)<2n2(2n-1)(n-2) < 2n^2(2n−1)(n−2)<2n2 in total.
  • Lemma 4.3 (p. 930): at most 4n24n^24n2 passes over the queue.

Significance

Corollary 4.4 is the combinatorial core of Theorem 4.5, which states that the FIFO algorithm runs in O(n3)O(n^3)O(n3) time. Theorem 4.2 shows that the remaining work of the implementation is O(nm)O(nm)O(nm) plus constant time per nonsaturating push. The bound of Corollary 4.4 is therefore what separates the O(n3)O(n^3)O(n3) FIFO method from the O(n2m)O(n^2 m)O(n2m) bound of the generic method, which matters on dense networks. The same pass-counting argument is reused for the parallel algorithm of §6 and underlies later analyses of highest-label and wave variants.

The results are proved in the paper. Formalizing them adds an analysis of a push–relabel algorithm, which the platform does not yet have. Its existing network-flow material states max-flow min-cut and Ford–Fulkerson termination in an arc-based model with nonnegative flows (the Introduction to Linear Optimization missions). The mission builds a precise operational model of the FIFO implementation, with edge lists, current edges and a queue carrying pass numbers, and states an explicit operation count for it. A companion mission in this series treats the generic algorithm's correctness and its (2n−1)(n−2)+2nm+4n2m(2n-1)(n-2) + 2nm + 4n^2m(2n−1)(n−2)+2nm+4n2m operation bound.

Difficulty

The obvious argument is the potential-function count of §3, over the sum of the labels of active vertices. It yields only 4n2m4n^2 m4n2m and does not use the queue discipline at all. The 4n34n^34n3 bound has to charge nonsaturating pushes to passes over the queue, and then bound the number of passes by the total growth of the labels. Neither step is visible in the generic algorithm, because both depend on the order in which vertices are processed.

Making this rigorous requires invariants of the implementation that the paper uses silently:

  • a vertex is in the queue exactly when it is active, and at most once;
  • pass numbers are nondecreasing along the queue;
  • current edges only move forward between relabelings.

Lemma 4.1 in particular depends on the current-edge scan: an edge passed over earlier stays inadmissible until vvv is relabeled.

Formalization scope

The Lean development works in namespace GoldbergTarjan.FIFO. Vertices form a type V with [Fintype V] [DecidableEq V], and nnn is Fintype.card V. Capacities are c : V → V → ℝ with c ≥ 0 and c v v = 0. Flows are antisymmetric real functions on all ordered pairs, not nonnegative arc flows. Excess is computed from the flow. Labels are in ℕ∞, and the empty minimum in relabel is ⊤.

The state of the algorithm consists of the flow, the labels, the current-edge index cur v into the edge list L v, and the queue Q : List (V × ℕ), each entry tagged with its pass number. Push/relabel (Fig. 3) is a total function, and a discharge (Fig. 4) is a relation carrying the number of push/relabel operations it performs. A run consists of the states S 0, …, S K with S 0 the initial state and consecutive states related by one discharge. The printed variant of Fig. 4, which stops as soon as vvv is relabeled, is the one formalized. Counts are natural numbers over all push/relabel operations of all discharges. The number of passes is the largest pass tag of a discharged entry.

All constants are explicit, exactly as printed:

  • 2n−12n-12n−1 (Lemmas 3.7, 3.8);
  • (2n−1)(n−2)(2n-1)(n-2)(2n−1)(n−2) and 2n22n^22n2 (Lemma 3.8);
  • 4n24n^24n2 (Lemma 4.3);
  • 4n34n^34n3 (Corollary 4.4).

No asymptotic notation is used, and no m≥n−1m \ge n-1m≥n−1 assumption is made.

A model without current edges, where relabeling happens whenever no push applies, would make Lemma 4.1 vacuous and change the algorithm. Pass numbers that are not propagated by the "added during pass iii" rule would make the pass count arbitrary. Both are ruled out by the definitions. A sorry-free check, outside the proposal, exhibits a three-vertex network with two legal discharges, two passes and no nonsaturating push, so the run hypotheses are satisfiable.

Reusable beyond this mission are the network, preflow, push and relabel definitions and Lemma 3.5, which is about an arbitrary preflow. Contributions welcome: invariants of FIFO runs (preflow, valid labeling, queue = active set, cur within bounds), proofs of the milestones, and the reduction of Corollary 4.4 to Lemma 4.3.

Selected references

  • A. V. Goldberg, R. E. Tarjan, A New Approach to the Maximum-Flow Problem, Journal of the ACM 35(4):921–940, 1988. https://doi.org/10.1145/48014.61051
  • A. V. Karzanov, Determining the maximal flow in a network by the method of preflows, Soviet Math. Doklady 15:434–437, 1974.
  • Y. Shiloach, U. Vishkin, An O(n² log n) parallel max-flow algorithm, Journal of Algorithms 3(2):128–146, 1982. https://doi.org/10.1016/0196-6774(82)90013-X
  • L. R. Ford, D. R. Fulkerson, Maximal flow through a network, Canadian Journal of Mathematics 8:399–404, 1956. https://doi.org/10.4153/CJM-1956-045-5
  • J. Edmonds, R. M. Karp, Theoretical improvements in algorithmic efficiency for network flow problems, Journal of the ACM 19(2):248–264, 1972. https://doi.org/10.1145/321694.321699
10 thms3 active usersReviewed
PreviousPage 1 of 4Next

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