Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Combinatorics

627 missions · 347 completed

The mathematics of finite and discrete structures — counting the arrangements of a set, deciding when a configuration meeting prescribed constraints can exist, and characterizing the patterns such structures are forced to contain. It encompasses enumerative and extremal combinatorics, graph theory, design theory, and additive combinatorics, with deep ties to algebra, probability, and computer science.

Missions

Open280Completed347All627
🏆Completed
Graph TheoryNumber Theory·Captain: Xiang Huang

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
Complexity TheoryGraph 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
🏆Completed
Convex OptimizationOperations ResearchOptimization·Captain: mikedeng1

Convexity and Steinitz's Exchange Property III: Fenchel-Type Min-Max Duality with Primal and Dual Integrality for M-Concave and M-Convex FunctionsResearch Paper

Motivation

Several classical min-max theorems of combinatorial optimization say that a discrete maximization problem and a continuous minimization problem have the same optimal value, and that both have integral optimal solutions when the data are integral. Edmonds' polymatroid intersection theorem (1970), Fujishige's Fenchel-type duality for submodular functions (1984), Frank's discrete separation theorem for a submodular/supermodular pair (1982), and the potential characterizations of weighted matroid intersection (Frank's weight splitting theorem, 1981; Iri and Tomizawa's criterion for the assignment problem, 1976) are instances. Murota's paper Convexity and Steinitz's exchange property, 1996 places all of them under one theorem: a Fenchel-type min-max formula for a pair of an M-concave and an M-convex function, with integrality on both sides.

Timeline:

  • 1970: Edmonds proves the polymatroid intersection theorem.
  • 1982: Frank proves the discrete separation theorem for submodular/supermodular set functions, with integrality.
  • 1984: Fujishige proves a Fenchel-type min-max theorem for submodular functions.
  • 1976–1981: Iri and Tomizawa characterize optimality for independent assignment by potentials; Frank proves the weight splitting theorem for weighted matroid intersection (1981).
  • Early 1990s: Dress and Wenzel introduce valuated matroids.
  • 1995–1996: Murota proves the valuated matroid intersection theorem (SIAM J. Discrete Math. 9, 1996) and the M-concave intersection theorem (Bonn report, 1995), and in the present paper the Fenchel-type duality (Theorem 6.4).
  • Later: the result becomes the central duality theorem of discrete convex analysis (Murota, Discrete Convex Analysis, SIAM, 2003).

Setting

Let VVV be a finite nonempty set. For u∈Vu\in Vu∈V, χu∈ZV\chi_u\in\mathbb Z^Vχu​∈ZV is its characteristic vector; for x∈RVx\in\mathbb R^Vx∈RV, supp⁡±(x)\operatorname{supp}^{\pm}(x)supp±(x) are the sets of coordinates where xxx is positive or negative, x(X)=∑v∈Xx(v)x(X)=\sum_{v\in X}x(v)x(X)=∑v∈X​x(v), and ⟨p,x⟩=∑vp(v)x(v)\langle p,x\rangle=\sum_v p(v)x(v)⟨p,x⟩=∑v​p(v)x(v).

A finite integral base set is a finite nonempty B⊆ZVB\subseteq\mathbb Z^VB⊆ZV such that for x,y∈Bx,y\in Bx,y∈B and u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) some v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) has x−χu+χv∈Bx-\chi_u+\chi_v\in Bx−χu​+χv​∈B. These are exactly the integer points of integral base polytopes of submodular systems. B‾\overline BB is the convex hull of BBB.

A function ω:B→R\omega:B\to\mathbb Rω:B→R has the exchange property (EXC), and is called M-concave, if for x,y∈Bx,y\in Bx,y∈B and u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) some v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) has x−χu+χv, y+χu−χv∈Bx-\chi_u+\chi_v,\ y+\chi_u-\chi_v\in Bx−χu​+χv​, y+χu​−χv​∈B and

ω(x)+ω(y)≤ω(x−χu+χv)+ω(y+χu−χv).\omega(x)+\omega(y)\le\omega(x-\chi_u+\chi_v)+\omega(y+\chi_u-\chi_v).ω(x)+ω(y)≤ω(x−χu​+χv​)+ω(y+χu​−χv​).

A function ζ\zetaζ is M-convex when −ζ-\zeta−ζ is M-concave.

For ω:B1→R\omega:B_1\to\mathbb Rω:B1​→R and ζ:B2→R\zeta:B_2\to\mathbb Rζ:B2​→R the concave conjugate and convex conjugate are

ω∘(p)=min⁡x∈B1(⟨p,x⟩−ω(x)),ζ∙(p)=max⁡x∈B2(⟨p,x⟩−ζ(x)),\omega^\circ(p)=\min_{x\in B_1}\big(\langle p,x\rangle-\omega(x)\big),\qquad \zeta^\bullet(p)=\max_{x\in B_2}\big(\langle p,x\rangle-\zeta(x)\big),ω∘(p)=x∈B1​min​(⟨p,x⟩−ω(x)),ζ∙(p)=x∈B2​max​(⟨p,x⟩−ζ(x)),

and the concave closure and convex closure are ω^(b)=inf⁡p(⟨p,b⟩−ω∘(p))\hat\omega(b)=\inf_p(\langle p,b\rangle-\omega^\circ(p))ω^(b)=infp​(⟨p,b⟩−ω∘(p)) and ζˇ(b)=sup⁡p(⟨p,b⟩−ζ∙(p))\check\zeta(b)=\sup_p(\langle p,b\rangle-\zeta^\bullet(p))ζˇ​(b)=supp​(⟨p,b⟩−ζ∙(p)); they are finite exactly on B1‾\overline{B_1}B1​​ and B2‾\overline{B_2}B2​​.

The primal problem maximizes ω(x)−ζ(x)\omega(x)-\zeta(x)ω(x)−ζ(x) over x∈B1∩B2x\in B_1\cap B_2x∈B1​∩B2​; the relaxed primal problem maximizes ω^(b)−ζˇ(b)\hat\omega(b)-\check\zeta(b)ω^(b)−ζˇ​(b) over b∈B1‾∩B2‾b\in\overline{B_1}\cap\overline{B_2}b∈B1​​∩B2​​; the dual problem minimizes ζ∙(p)−ω∘(p)\zeta^\bullet(p)-\omega^\circ(p)ζ∙(p)−ω∘(p) over p∈RVp\in\mathbb R^Vp∈RV. A maximum over an empty family is −∞-\infty−∞.

Formalization targets

Goal: Theorem 6.4

If ω\omegaω and −ζ-\zeta−ζ satisfy (EXC), then

max⁡x∈B1∩B2(ω(x)−ζ(x))=max⁡b∈B1‾∩B2‾(ω^(b)−ζˇ(b))=inf⁡p∈RV(ζ∙(p)−ω∘(p)),\max_{x\in B_1\cap B_2}\big(\omega(x)-\zeta(x)\big)=\max_{b\in\overline{B_1}\cap\overline{B_2}}\big(\hat\omega(b)-\check\zeta(b)\big)=\inf_{p\in\mathbb R^V}\big(\zeta^\bullet(p)-\omega^\circ(p)\big),x∈B1​∩B2​max​(ω(x)−ζ(x))=b∈B1​​∩B2​​max​(ω^(b)−ζˇ​(b))=p∈RVinf​(ζ∙(p)−ω∘(p)),

with (P1) a finite dual infimum forces B1∩B2≠∅B_1\cap B_2\neq\emptysetB1​∩B2​=∅, and (P2) if B1∩B2≠∅B_1\cap B_2\neq\emptysetB1​∩B2​=∅ all values are finite and equal and the infimum is attained. If ω,ζ\omega,\zetaω,ζ are integer-valued, the infimum may be taken over p∈ZVp\in\mathbb Z^Vp∈ZV and is attained there when finite.

Milestones

  1. Lemma 6.3 (weak duality): for arbitrary ω,ζ\omega,\zetaω,ζ on finite nonempty sets, primal ≤\le≤ relaxed === dual (the Fenchel identity (6.5)).
  2. Lemma 6.1: (−f)∘(p)=−f∙(−p)(-f)^\circ(p)=-f^\bullet(-p)(−f)∘(p)=−f∙(−p) and (−f)∧=−fˇ(-f)^\wedge=-\check f(−f)∧=−fˇ​ on B‾\overline BB.
  3. Lemma 4.5: an M-concave ω\omegaω satisfies ω^=ω\hat\omega=\omegaω^=ω on BBB.
  4. Theorem 2.1: (B1) is equivalent to being the integer points of an integral submodular (or supermodular) base polytope, with the describing functions max⁡x∈Bx(X)\max_{x\in B}x(X)maxx∈B​x(X) and min⁡x∈Bx(X)\min_{x\in B}x(X)minx∈B​x(X).
  5. Theorem 6.5 (Frank's discrete separation theorem, cited in the paper).
  6. Lemma 6.7: four equivalent forms of boundedness of the dual problem.
  7. Theorem 6.6 (the M-concave intersection theorem, cited in the paper): optimality of x∗x^*x∗ for ω1+ω2\omega_1+\omega_2ω1​+ω2​ is equivalent to a potential p∗p^*p∗ with x∗x^*x∗ maximizing both ω1[−p∗]\omega_1[-p^*]ω1​[−p∗] and ω2[p∗]\omega_2[p^*]ω2​[p∗], integral when the data are.

Significance

The formula gives, in one statement, the integrality of an optimal solution of the relaxed primal problem (the essential content of the first half, as the paper observes on p. 296) and of the dual problem. The paper presents it as a unification of two groups of theorems: Edmonds' polymatroid intersection theorem, Fujishige's Fenchel-type duality and Frank's discrete separation theorem on one side, and Iri and Tomizawa's potential characterization for independent assignment with its extensions by Fujishige and Frank (weight splitting) on the other. In the paper it yields the primal and dual separation theorems (Theorems 6.8, 6.9) and the convolution results (Theorems 6.10, 6.11), and it is the prototype of the Fenchel-type duality of discrete convex analysis.

All results here are proved in the literature; none is known to be formalized. Mathlib has no submodular base polytopes, no matroid intersection theorem and no discrete convex analysis. A formal proof of Theorem 6.4 would also require formal proofs of the two cited results, Frank's discrete separation theorem and the M-concave intersection theorem, which the paper uses without proof.

Difficulty

Lemma 6.3 is polyhedral convex duality and holds for any functions. The content is equality with the integral problem: the relaxed maximum over the polytope B1‾∩B2‾\overline{B_1}\cap\overline{B_2}B1​​∩B2​​ must be attained at an integer point. For general finite sets it is not, and the intersection of two integral polytopes generally has fractional vertices. Both the integrality of B1‾∩B2‾\overline{B_1}\cap\overline{B_2}B1​​∩B2​​ (Edmonds) and the existence of an integral optimal potential depend on the exchange structure; a direct argument from the definitions of conjugates does not see it. The dual integrality claim, that ppp can be taken integral, is again specific to (EXC) and fails for general concave extensions.

Formalization scope

Lean conventions, all in namespace SteinitzExchange.Duality:

  • VVV is a type with [Fintype V] [DecidableEq V] [Nonempty V]; integer vectors are V → ℤ, real vectors V → ℝ; finite sets of integer vectors are Finset (V → ℤ).
  • A function on BBB is a total (V → ℤ) → ℝ used only at points of BBB. M-convexity of ζ\zetaζ is (EXC) for fun x => -ζ x; ω\omegaω lives on B1B_1B1​ and ζ\zetaζ on B2B_2B2​, which are distinct sets in general.
  • Conjugates are real-valued min/max over the finite set. The closures are real ⨅/⨆ over p∈RVp\in\mathbb R^Vp∈RV and are evaluated only on the convex hulls, where they equal the paper's values; off the hulls they carry a junk value instead of ∓∞\mp\infty∓∞, which no statement uses.
  • The three optimal values are in EReal, as suprema and infima of coerced reals, so no ∞−∞\infty-\infty∞−∞ occurs. EReal's supremum of the empty family is −∞-\infty−∞, the paper's convention. The dual infimum is never a real ⨅ (which would return 000 when unbounded and make (P1) meaningless).
  • Every "max" of the page includes attainment: (P2) asserts points xxx, bbb, ppp at which the three values are achieved; the integral dual infimum is attained when it is not −∞-\infty−∞.
  • "Integer-valued" means ω(x)∈Z\omega(x)\in\mathbb Zω(x)∈Z on B1B_1B1​ and ζ(x)∈Z\zeta(x)\in\mathbb Zζ(x)∈Z on B2B_2B2​; integral potentials and separating vectors are V → ℤ.
  • Theorem 2.1's "∀X⊂V\forall X\subset V∀X⊂V" is read as all X⊆VX\subseteq VX⊆V.

Formalizations that would trivialize the goal are excluded: an unrestricted real infimum for the dual, a convex closure built from ζ∘\zeta^\circζ∘ instead of ζ∙\zeta^\bulletζ∙, a single base set for both functions, and a relaxed maximum taken over all of RV\mathbb R^VRV instead of B1‾∩B2‾\overline{B_1}\cap\overline{B_2}B1​​∩B2​​.

Needed infrastructure: finite convex hulls and polyhedral Fenchel duality, submodular base polytopes and their integrality, Frank's separation theorem, and the valuated intersection theorem. The submodular-system layer (Theorem 2.1, Theorem 6.5) is reusable beyond this mission. Contributions to any milestone, including proofs of the two cited theorems, are welcome.

Selected references

  • K. Murota, Convexity and Steinitz's exchange property, Advances in Mathematics 124 (1996) 272–311. https://doi.org/10.1006/aima.1996.0084
  • K. Murota, Valuated matroid intersection I: optimality criteria, SIAM J. Discrete Math. 9 (1996) 545–561.
  • K. Murota, Submodular flow problem with a nonseparable cost function, Report 95843-OR, Forschungsinstitut für Diskrete Mathematik, Universität Bonn, 1995 (source of Theorem 6.6).
  • A. Frank, An algorithm for submodular functions on graphs, Annals of Discrete Mathematics 16 (1982) 97–120 (source of Theorem 6.5).
  • A. Frank, A weighted matroid intersection algorithm, J. Algorithms 2 (1981) 328–336.
  • J. Edmonds, Submodular functions, matroids and certain polyhedra, in: Combinatorial Structures and Their Applications, Gordon and Breach, New York, 1970, 69–87.
  • S. Fujishige, Theory of submodular programs: a Fenchel-type min-max theorem and subgradients of submodular functions, Mathematical Programming 29 (1984) 142–155.
  • M. Iri and N. Tomizawa, An algorithm for finding an optimal "independent assignment", J. Oper. Res. Soc. Japan 19 (1976) 32–57.
  • K. Murota, Discrete Convex Analysis, SIAM, 2003. https://doi.org/10.1137/1.9780898718508
18 thms7 active usersReviewed
🏆Completed
Convex OptimizationOperations ResearchOptimization·Captain: mikedeng1

Convexity and Steinitz's Exchange Property II: The Local Supermodularity Theorem for the Concave ConjugateResearch Paper

Motivation

Matroids and their integral generalizations, integral base polytopes, are the combinatorial structures on which the greedy algorithm is exact. Edmonds' theory relates them to submodular and supermodular set functions: a polytope is a base polytope exactly when its support function, restricted to 0/10/10/1 vectors, is supermodular and the greedy formula evaluates it everywhere. Dress and Wenzel's valuated matroids (1990) and Murota's M-concave functions carry the exchange axiom from sets to functions on sets. This paper (Adv. Math. 124, 1996) sets up the resulting theory of discrete concave functions on base sets, later developed into discrete convex analysis (Murota, Discrete Convex Analysis, SIAM 2003).

The question behind this mission is how the set-level correspondence between exchange and supermodularity extends to functions. Section 5 of the paper answers it with the Local Supermodularity Theorem: the exchange property of a function is a supermodularity property of its concave conjugate, holding locally at every point.

Setting

Let VVV be a finite nonempty set, n=∣V∣n=|V|n=∣V∣. For u∈Vu\in Vu∈V let χu∈ZV\chi_u\in\mathbb Z^Vχu​∈ZV be the unit vector, for X⊆VX\subseteq VX⊆V let χX\chi_XχX​ be its characteristic vector, x(X)=∑v∈Xx(v)x(X)=\sum_{v\in X}x(v)x(X)=∑v∈X​x(v), and ⟨p,x⟩=∑vp(v)x(v)\langle p,x\rangle=\sum_v p(v)x(v)⟨p,x⟩=∑v​p(v)x(v). For a finite B⊆ZVB\subseteq\mathbb Z^VB⊆ZV, B‾\overline BB is its convex hull.

A finite integral base set is a finite nonempty B⊆ZVB\subseteq\mathbb Z^VB⊆ZV such that

(B1)x,y∈B, u∈supp⁡+(x−y) ⇒ ∃v∈supp⁡−(x−y): x−χu+χv∈B.\text{(B1)}\quad x,y\in B,\ u\in\operatorname{supp}^+(x-y)\ \Rightarrow\ \exists v\in\operatorname{supp}^-(x-y):\ x-\chi_u+\chi_v\in B.(B1)x,y∈B, u∈supp+(x−y) ⇒ ∃v∈supp−(x−y): x−χu​+χv​∈B.

A function ω:B→R\omega:B\to\mathbb Rω:B→R satisfies the exchange property (EXC) (is M-concave) if for x,y∈Bx,y\in Bx,y∈B and u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) there is v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) with x−χu+χv, y+χu−χv∈Bx-\chi_u+\chi_v,\ y+\chi_u-\chi_v\in Bx−χu​+χv​, y+χu​−χv​∈B and ω(x)+ω(y)≤ω(x−χu+χv)+ω(y+χu−χv)\omega(x)+\omega(y)\le\omega(x-\chi_u+\chi_v)+\omega(y+\chi_u-\chi_v)ω(x)+ω(y)≤ω(x−χu​+χv​)+ω(y+χu​−χv​). Write ω[p](x)=ω(x)+⟨p,x⟩\omega[p](x)=\omega(x)+\langle p,x\rangleω[p](x)=ω(x)+⟨p,x⟩ and argmax⁡(g)\operatorname{argmax}(g)argmax(g) for the maximizers of ggg on BBB.

The support function of BBB is ψ∘(p)=min⁡{⟨p,x⟩∣x∈B}\psi^\circ(p)=\min\{\langle p,x\rangle\mid x\in B\}ψ∘(p)=min{⟨p,x⟩∣x∈B}. A positively homogeneous h:RV→Rh:\mathbb R^V\to\mathbb Rh:RV→R is "matroidal" if

  • (C1) X↦h(χX)X\mapsto h(\chi_X)X↦h(χX​) is supermodular, and
  • (C2) h(p)=∑j=1n(pj−pj+1) h(χVj)h(p)=\sum_{j=1}^n(p_j-p_{j+1})\,h(\chi_{V_j})h(p)=∑j=1n​(pj​−pj+1​)h(χVj​​) whenever V={v1,…,vn}V=\{v_1,\dots,v_n\}V={v1​,…,vn​} with p(v1)≥⋯≥p(vn)p(v_1)\ge\dots\ge p(v_n)p(v1​)≥⋯≥p(vn​), pj=p(vj)p_j=p(v_j)pj​=p(vj​), Vj={v1,…,vj}V_j=\{v_1,\dots,v_j\}Vj​={v1​,…,vj​}, pn+1=0p_{n+1}=0pn+1​=0.

The concave conjugate is ω∘(p)=min⁡{⟨p,x⟩−ω(x)∣x∈B}\omega^\circ(p)=\min\{\langle p,x\rangle-\omega(x)\mid x\in B\}ω∘(p)=min{⟨p,x⟩−ω(x)∣x∈B}, the concave closure is ω^(b)=inf⁡p{⟨p,b⟩−ω∘(p)}\hat\omega(b)=\inf_p\{\langle p,b\rangle-\omega^\circ(p)\}ω^(b)=infp​{⟨p,b⟩−ω∘(p)}, the subdifferential is ∂ω∘(p0)={b∣ω∘(p)−ω∘(p0)≤⟨p−p0,b⟩ ∀p}\partial\omega^\circ(p_0)=\{b\mid\omega^\circ(p)-\omega^\circ(p_0)\le\langle p-p_0,b\rangle\ \forall p\}∂ω∘(p0​)={b∣ω∘(p)−ω∘(p0​)≤⟨p−p0​,b⟩ ∀p}, and the localization of ω∘\omega^\circω∘ at p0p_0p0​ is L^(ω∘,p0)(p)=inf⁡{⟨p,b⟩∣b∈∂ω∘(p0)}\hat L(\omega^\circ,p_0)(p)=\inf\{\langle p,b\rangle\mid b\in\partial\omega^\circ(p_0)\}L^(ω∘,p0​)(p)=inf{⟨p,b⟩∣b∈∂ω∘(p0​)}.

Formalization targets

Goal: the Local Supermodularity Theorem (Theorem 5.3, corrected)

For ω\omegaω on a finite integral base set BBB,

ω satisfies (EXC)  ⟺  (ω=ω^ on B) and (L^(ω∘,p0) is "matroidal" for every p0∈RV).\omega\ \text{satisfies (EXC)}\iff\Big(\omega=\hat\omega\ \text{on}\ B\Big)\ \text{and}\ \Big(\hat L(\omega^\circ,p_0)\ \text{is "matroidal" for every}\ p_0\in\mathbb R^V\Big).ω satisfies (EXC)⟺(ω=ω^ on B) and (L^(ω∘,p0​) is "matroidal" for every p0​∈RV).

The printed Theorem 5.3 has only the second condition on the right. Its "only if" direction holds as printed; its "if" direction is false without the first condition, and a separate item of the mission states the counterexample: B={(2,0),(1,1),(0,2)}B=\{(2,0),(1,1),(0,2)\}B={(2,0),(1,1),(0,2)}, ω=(0,−10,0)\omega=(0,-10,0)ω=(0,−10,0).

Milestones

  1. Theorem 2.1: (B1) is equivalent to BBB being the integer points of an integral submodular (equivalently, supermodular) system, whose defining function is determined by BBB.
  2. Theorem 5.1: if B=ZV∩B‾B=\mathbb Z^V\cap\overline BB=ZV∩B, then BBB satisfies (B1) iff ψ∘\psi^\circψ∘ is "matroidal".
  3. Lemma 5.2: sums of "matroidal" functions are "matroidal".
  4. Theorem 4.4: (EXC) holds iff every argmax⁡(ω[p])\operatorname{argmax}(\omega[p])argmax(ω[p]) satisfies (B1).
  5. Eq. (5.12): L^(ω∘,p0)(p)=min⁡{⟨p,x⟩∣x∈argmax⁡(ω[−p0])}\hat L(\omega^\circ,p_0)(p)=\min\{\langle p,x\rangle\mid x\in\operatorname{argmax}(\omega[-p_0])\}L^(ω∘,p0​)(p)=min{⟨p,x⟩∣x∈argmax(ω[−p0​])}.

Significance

The result. Theorem 5.3 is the function-level version of Theorem 5.1. Condition (C1) is a supermodularity condition, so the theorem expresses (EXC) as "a collection of local supermodularity" properties of ω∘\omega^\circω∘, in the same way that (B1) corresponds to supermodularity of a support function. In the paper this characterization of the conjugate side underlies the Fenchel-type duality of Section 6, and more generally the conjugacy between M-concave and L-convex functions in discrete convex analysis.

The formalization. No part of this theory is formalized in Lean or on this platform: base sets, (EXC), "matroidal" functions and concave conjugates of functions on base sets are all new. The mission also corrects the published statement: the reduction from localizations to base sets needs every integer point of conv⁡(argmax⁡ ω[−p0])\operatorname{conv}(\operatorname{argmax}\,\omega[-p_0])conv(argmaxω[−p0​]) to be a maximizer, and the concave-closure condition supplies this. A machine-checked proof would settle both the corrected theorem and the counterexample. Theorem 2.1 and Lemma 5.2 are classical but have no formal proof either.

Difficulty

ω∘\omega^\circω∘ depends only on the concave closure ω^\hat\omegaω^, so any characterization of (EXC) through ω∘\omega^\circω∘ alone cannot see values of ω\omegaω below ω^\hat\omegaω^. That is why the goal needs the extra clause. The "only if" direction needs the full theory of Section 4: M-concave functions coincide with their concave closure, and all their maximizer sets are base sets. Theorem 5.1 needs the greedy algorithm on integral base polytopes, together with the fact that the base polytope of an integral supermodular function has integral vertices. Theorem 2.1 is the folklore statement that polyhedral and exchange descriptions agree, and the paper does not prove it. Eq. (5.12) needs the subdifferential of a finite minimum of affine functions to be the convex hull of the active gradients, stated globally rather than only near p0p_0p0​.

Formalization scope

Integer vectors are V → ℤ, real vectors V → ℝ, with [Fintype V] [DecidableEq V] [Nonempty V]. A finite subset of ZV\mathbb Z^VZV is a Finset (V → ℤ). A function on BBB is a total function (V → ℤ) → ℝ whose values off BBB are never used. The mission commits to the following readings:

  • Minima. ψ∘\psi^\circψ∘, ω∘\omega^\circω∘ are real infima over the finite set BBB, hence minima for nonempty BBB (every statement has BBB nonempty). ω^\hat\omegaω^ is a real infimum used only at points of BBB, where it is bounded below.
  • Localization. L^\hat LL^ is a real sInf over the subdifferential, defined by (5.8)–(5.9) exactly, not by the formula (5.12). Eq. (5.12) is stated with IsLeast, so it asserts attainment, not just the value.
  • (C2). It is required for every bijection Fin n ≃ V along which ppp is non-increasing. This is equivalent to "for some" such indexing. "Matroidal" includes positive homogeneity but not concavity.
  • Theorem 2.1. The page's "∀X⊂V\forall X\subset V∀X⊂V" is read as all X⊆VX\subseteq VX⊆V. The set functions are integer-valued, and "Moreover" is read strongly: every fff (resp. ggg) as in (b) (resp. (c)) equals the displayed max (resp. min).
  • Theorem 4.4. "argmax⁡(ω[p])‾\overline{\operatorname{argmax}(\omega[p])}argmax(ω[p])​ is an integral base polytope" is read as "argmax⁡(ω[p])\operatorname{argmax}(\omega[p])argmax(ω[p]) satisfies (B1)", following Lemma 4.3 and the proof of Theorem 5.3. The literal convex-hull reading makes the "if" direction false (same counterexample).
  • Theorem 5.1 keeps the page's hypothesis B=ZV∩B‾B=\mathbb Z^V\cap\overline BB=ZV∩B.

Trivializing formalizations are ruled out. Defining L^\hat LL^ by (5.12) would reduce the goal to Theorems 4.4 and 5.1. A "matroidal" without (C2) would be satisfied by support functions of non-base sets. An ω∘\omega^\circω∘ taken as a supremum would reverse the sign conventions.

The development needs: the greedy algorithm and integrality for integral base polytopes, supergradients of polyhedral concave functions, and the Section 4 results of the paper (concave closure of M-concave functions, Lemma 4.3). The base-set and "matroidal" layers can be reused beyond this mission. Proofs of any milestone, of the counterexample, and a proof of the "only if" direction on its own are all welcome.

Selected references

  • K. Murota, Convexity and Steinitz's exchange property, Advances in Mathematics 124 (1996), 272–311. https://doi.org/10.1006/aima.1996.0084
  • A. W. M. Dress, W. Wenzel, Valuated matroids: a new look at the greedy algorithm, Applied Mathematics Letters 3 (1990), 33–35.
  • S. Fujishige, Submodular Functions and Optimization, 2nd ed., Annals of Discrete Mathematics 58, Elsevier, 2005.
  • L. Lovász, Submodular functions and convexity, in Mathematical Programming: The State of the Art, Springer, 1983, 235–257. https://doi.org/10.1007/978-3-642-68874-4_10
  • K. Murota, Discrete Convex Analysis, SIAM, 2003. https://doi.org/10.1137/1.9780898718508
15 thms7 active usersReviewed
🏆Completed
Convex OptimizationOperations ResearchOptimization·Captain: mikedeng1

Convexity and Steinitz's Exchange Property I: The Extension Theorem — M-Concavity Is Concave Extendability with Integral Base Polytope MaximizersResearch Paper

Motivation

Linear optimization over the bases of a matroid, over the integer points of a polymatroid, or over the flows of a network is well understood: the greedy algorithm is exact, the feasible sets are the integer points of polytopes described by submodular functions, and min-max theorems of Edmonds and Frank hold with integrality. Nonlinear objectives on the same sets are much less uniform. Valuated matroids (Dress and Wenzel, 1990; see Murota 2003) showed that a quantitative form of the Steinitz exchange axiom is exactly what keeps the greedy algorithm exact for a nonlinear weight. Kazuo Murota's paper Convexity and Steinitz's Exchange Property (Adv. Math. 124 (1996) 272–311) extends this exchange axiom from matroid bases to the integer points of arbitrary integral base polytopes, names the resulting functions M-concave, and proves that they are the discrete counterpart of concave functions. The paper is the starting point of discrete convex analysis (Murota, Discrete Convex Analysis, SIAM 2003), which is now used in auction theory (gross-substitutes valuations are M♮-concave), inventory and resource allocation, and combinatorial optimization.

This mission covers the first of the paper's three characterizations of M-concavity: the Extension Theorem (Theorem 4.6).

Setting

Let VVV be a finite nonempty set. For u∈Vu\in Vu∈V, χu∈ZV\chi_u\in\mathbb Z^Vχu​∈ZV is the characteristic vector of uuu. For x∈RVx\in\mathbb R^Vx∈RV, supp⁡+(x)={v∣x(v)>0}\operatorname{supp}^+(x)=\{v\mid x(v)>0\}supp+(x)={v∣x(v)>0}, supp⁡−(x)={v∣x(v)<0}\operatorname{supp}^-(x)=\{v\mid x(v)<0\}supp−(x)={v∣x(v)<0}, ∥x∥=∑v∣x(v)∣\|x\|=\sum_v|x(v)|∥x∥=∑v​∣x(v)∣, and ⟨p,x⟩=∑vp(v)x(v)\langle p,x\rangle=\sum_v p(v)x(v)⟨p,x⟩=∑v​p(v)x(v).

A finite integral base set is a finite nonempty B⊆ZVB\subseteq\mathbb Z^VB⊆ZV such that for x,y∈Bx,y\in Bx,y∈B and u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) there is v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) with x−χu+χv∈Bx-\chi_u+\chi_v\in Bx−χu​+χv​∈B (axiom (B1)). Examples are the incidence vectors of the bases of a matroid. Its convex hull B‾\overline BB is an integral base polytope; in general, an integral base polytope is the convex hull of some finite integral base set.

A function ω:B→R\omega:B\to\mathbb Rω:B→R satisfies the exchange property (EXC), and is called M-concave, if for all x,y∈Bx,y\in Bx,y∈B and u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) there is v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) with x−χu+χv∈Bx-\chi_u+\chi_v\in Bx−χu​+χv​∈B, y+χu−χv∈By+\chi_u-\chi_v\in By+χu​−χv​∈B and

ω(x)+ω(y)≤ω(x−χu+χv)+ω(y+χu−χv).\omega(x)+\omega(y)\le\omega(x-\chi_u+\chi_v)+\omega(y+\chi_u-\chi_v).ω(x)+ω(y)≤ω(x−χu​+χv​)+ω(y+χu​−χv​).

The local exchange property (EXCloc_{\mathrm{loc}}loc​) asks only, for x,y∈Bx,y\in Bx,y∈B with ∥x−y∥=4\|x-y\|=4∥x−y∥=4, for some u∈supp⁡+(x−y)u\in\operatorname{supp}^+(x-y)u∈supp+(x−y) and some v∈supp⁡−(x−y)v\in\operatorname{supp}^-(x-y)v∈supp−(x−y) with the same conclusion.

For p∈RVp\in\mathbb R^Vp∈RV, ω[p](x)=ω(x)+⟨p,x⟩\omega[p](x)=\omega(x)+\langle p,x\rangleω[p](x)=ω(x)+⟨p,x⟩, and argmax⁡(ω)={x∈B∣ω(x)≥ω(y) ∀y∈B}\operatorname{argmax}(\omega)=\{x\in B\mid\omega(x)\ge\omega(y)\ \forall y\in B\}argmax(ω)={x∈B∣ω(x)≥ω(y) ∀y∈B}. For any g:B→Rg:B\to\mathbb Rg:B→R, the concave conjugate is g∘(p)=min⁡x∈B(⟨p,x⟩−g(x))g^\circ(p)=\min_{x\in B}(\langle p,x\rangle-g(x))g∘(p)=minx∈B​(⟨p,x⟩−g(x)) and the concave closure is g^(b)=inf⁡p∈RV(⟨p,b⟩−g∘(p))\hat g(b)=\inf_{p\in\mathbb R^V}(\langle p,b\rangle-g^\circ(p))g^​(b)=infp∈RV​(⟨p,b⟩−g∘(p)), a concave function that is finite exactly on B‾\overline BB. A function ωˉ:B‾→R\bar\omega:\overline B\to\mathbb Rωˉ:B→R extends ω\omegaω if ωˉ=ω\bar\omega=\omegaωˉ=ω on BBB.

Formalization targets

Goal: the Extension Theorem (Theorem 4.6)

For a finite integral base set BBB and ω:B→R\omega:B\to\mathbb Rω:B→R,

ω satisfies (EXC)  ⟺  ∃ ωˉ:B‾→R concave, ωˉ∣B=ω, ∀p: argmax⁡B‾(ωˉ[p]) is an integral base polytope.\omega\ \text{satisfies (EXC)}\iff\exists\,\bar\omega:\overline B\to\mathbb R\ \text{concave},\ \bar\omega|_B=\omega,\ \forall p:\ \operatorname{argmax}_{\overline B}(\bar\omega[p])\ \text{is an integral base polytope}.ω satisfies (EXC)⟺∃ωˉ:B→R concave, ωˉ∣B​=ω, ∀p: argmaxB​(ωˉ[p]) is an integral base polytope.

Milestones

  • Lemma 3.2 (p. 282): under (EXCloc_{\mathrm{loc}}loc​), for y=x−χu0−χu1+χv0+χv1∈By=x-\chi_{u_0}-\chi_{u_1}+\chi_{v_0}+\chi_{v_1}\in By=x−χu0​​−χu1​​+χv0​​+χv1​​∈B, ωp(y)−ωp(x)≤max⁡(π00+π11,π01+π10)\omega_p(y)-\omega_p(x)\le\max(\pi_{00}+\pi_{11},\pi_{01}+\pi_{10})ωp​(y)−ωp​(x)≤max(π00​+π11​,π01​+π10​) with πij=ωp(x−χui+χvj)−ωp(x)\pi_{ij}=\omega_p(x-\chi_{u_i}+\chi_{v_j})-\omega_p(x)πij​=ωp​(x−χui​​+χvj​​)−ωp​(x) (−∞-\infty−∞ off BBB).
  • Theorem 3.1 (p. 282): (EXC)   ⟺  \iff⟺ (EXCloc_{\mathrm{loc}}loc​).
  • Theorem 2.2 (p. 280): (EXC) for ω\omegaω implies (EXC) for every ω[p]\omega[p]ω[p].
  • Lemma 4.3 (p. 285): under (EXC), argmax⁡(ω)\operatorname{argmax}(\omega)argmax(ω) is an integral base set.
  • Lemma 4.1 (p. 285): g^≥g\hat g\ge gg^​≥g on BBB; max⁡B‾g^=max⁡Bg\max_{\overline B}\hat g=\max_B gmaxB​g^​=maxB​g; argmax⁡(g^)=argmax⁡(g)‾\operatorname{argmax}(\hat g)=\overline{\operatorname{argmax}(g)}argmax(g^​)=argmax(g)​.
  • Lemma 4.2 (p. 285): (g[p0])∘(p)=g∘(p−p0)(g[p_0])^\circ(p)=g^\circ(p-p_0)(g[p0​])∘(p)=g∘(p−p0​) and (g[p0])∧=g^+⟨p0,⋅⟩(g[p_0])^\wedge=\hat g+\langle p_0,\cdot\rangle(g[p0​])∧=g^​+⟨p0​,⋅⟩ on B‾\overline BB.
  • Theorem 4.4 (p. 286): (EXC)   ⟺  \iff⟺ argmax⁡(ω[p])\operatorname{argmax}(\omega[p])argmax(ω[p]) is an integral base set for every ppp.
  • Lemma 4.5 (p. 288): under (EXC), ω^=ω\hat\omega=\omegaω^=ω on BBB.

Significance

The Extension Theorem identifies a combinatorial axiom with a convex-analytic property: M-concave functions are exactly the restrictions to lattice points of concave functions on the base polytope whose linear perturbations are all maximized on integral base polytopes. It is the reason the M-concave class supports a convex-analysis-style theory at all: local optimality implies global optimality, maximizers of linear perturbations are well behaved, and conjugacy (the paper's Theorems 5.3 and 6.4, separate missions of this series) can be developed. Theorem 3.1 on its own is widely used to verify M-concavity in applications, since it reduces the exchange axiom to pairs at distance four.

All results here were proved in 1996. To our knowledge none of them has a machine-checked proof; Mathlib has matroids on sets but no integral base sets in ZV\mathbb Z^VZV, no M-concave functions and no concave closure of a function on a finite set. A formal proof of this chain would be a first formal development of discrete convex analysis.

Difficulty

The equivalence of (EXC) with its local version (Theorem 3.1) is not a routine induction on ∥x−y∥\|x-y\|∥x−y∥: the exchange inequality for a far pair does not follow from the inequalities along a path of distance-4 pairs, because the exchange partner vvv must be chosen consistently with the prescribed uuu. For the "if" direction of Theorem 4.4, knowing that every maximizer set is an integral base set says nothing directly about the values of ω\omegaω at non-maximizing points; turning this global information on maximizers into the local inequality (EXCloc_{\mathrm{loc}}loc​) requires a supporting hyperplane of the concave closure at a well-chosen point and the integrality of the intersection of an integral base polytope with a box (a cited result on submodular systems). Theorem 4.6 then needs the concave closure to agree with ω\omegaω on BBB (Lemma 4.5), which fails for general ω\omegaω.

Formalization scope

All declarations live in the namespace SteinitzExchange.Extension. The ground set is a type V with [Fintype V] [DecidableEq V] [Nonempty V]; integer vectors are V → ℤ, real vectors V → ℝ, and toReal embeds the former into the latter. BBB is a Finset (V → ℤ). A function ω:B→R\omega:B\to\mathbb Rω:B→R is a total (V → ℤ) → ℝ whose values are only ever read at points required to be in BBB. B‾\overline BB is Mathlib's convexHull ℝ of the image of BBB. Pinned readings:

  1. Integral base polytope means the convex hull of a finite nonempty set satisfying (B1) (by the paper's Theorem 2.1 this is its meaning), not "a polytope with integer vertices".
  2. The concave closure is a real infimum; it is the paper's value on B‾\overline BB and a junk value 000 off B‾\overline BB (the paper's −∞-\infty−∞), so every statement uses it only on B‾\overline BB. argmax⁡(g^)\operatorname{argmax}(\hat g)argmax(g^​) and argmax⁡(ωˉ[p])\operatorname{argmax}(\bar\omega[p])argmax(ωˉ[p]) range over B‾\overline BB only; the concave conjugate is a minimum over the nonempty finite BBB.
  3. Lemma 3.2's maximum with −∞-\infty−∞ entries is stated as: for one of the two pairings both exchanged points lie in BBB and the bound holds for that pairing.
  4. Theorem 4.4 and Lemma 4.3 conclude that argmax⁡(ω[p])\operatorname{argmax}(\omega[p])argmax(ω[p]) is itself an integral base set. Read literally ("its convex hull is an integral base polytope") the "if" direction of Theorem 4.4 is false: B={(2,0),(1,1),(0,2)}B=\{(2,0),(1,1),(0,2)\}B={(2,0),(1,1),(0,2)} with ω=(0,−1,0)\omega=(0,-1,0)ω=(0,−1,0) is a counterexample. The paper's proof, its gloss in Lemma 4.3 and its use on p. 292 all take the integral-base-set reading. Theorem 4.6 needs no such adjustment and is stated as printed.
  5. Theorem 2.2 carries the standing assumption of §2.3 that ω\omegaω satisfies (EXC).

Trivializing formalizations are ruled out: the extension ωˉ\bar\omegaωˉ must agree with ω\omegaω on BBB and be concave on B‾\overline BB, the argmax is over B‾\overline BB and not over RV\mathbb R^VRV, and an integral base polytope is never empty.

A complete development needs basic facts on integral base sets (the equivalence of (B1) with the simultaneous exchange (B2), B=ZV∩B‾B=\mathbb Z^V\cap\overline BB=ZV∩B, and the paper's cited Theorem 2.1 relating them to submodular functions), the representation (4.3) of the concave closure as a maximum of convex combinations, and supporting hyperplanes of polyhedral concave functions. The layer of integral base sets and M-concave functions is reusable for the two other missions of this series and for any later formalization of discrete convex analysis; contributions of general-purpose lemmas about it are welcome.

Selected references

  • K. Murota, Convexity and Steinitz's exchange property, Advances in Mathematics 124 (1996) 272–311. https://doi.org/10.1006/aima.1996.0084
  • K. Murota, Discrete Convex Analysis, SIAM Monographs on Discrete Mathematics and Applications, 2003. https://doi.org/10.1137/1.9780898718508
21 thms7 active usersReviewed
🏆Completed
Number Theory·Captain: aarontcao

The Komlos-Sulyok-Szemeredi bound: every finite set of reals has a Sidon subset of size c sqrt nResearch Paper

Call a set of reals a Sidon set when all its pairwise sums are distinct: if a+b=c+da + b = c + da+b=c+d with all four in the set, then {a,b}={c,d}\{a,b\} = \{c,d\}{a,b}={c,d}.

The goal. There is an absolute constant c>0c > 0c>0 such that every finite set XXX of positive reals contains a Sidon subset SSS with ∣S∣≥c∣X∣|S| \ge c\sqrt{|X|}∣S∣≥c∣X∣​.

This is the lower bound half of Erdos problem 530, which Riddell posed and which asks for the order of the largest guaranteed Sidon subset. That problem is open: it asks whether the guarantee is asymptotically N1/2N^{1/2}N1/2, and the constant is not known. What is settled is the order, by Komlos, Sulyok, and Szemeredi, Linear problems in combinatorial number theory, Acta Math. Acad. Sci. Hungar. 26 (1975) 113-121, as a case of a general theorem about linear equations. Erdos had previously observed the cube-root lower bound and the matching (1+o(1))N1/2(1+o(1))N^{1/2}(1+o(1))N1/2 upper bound from A={1,…,N}A = \{1, \dots, N\}A={1,…,N}. A second and much shorter proof is in Bailleul and Riblet, arXiv:2605.03181.

The exponent is the whole problem

A one-paragraph argument gives ∣S∣≥c∣X∣1/3|S| \ge c|X|^{1/3}∣S∣≥c∣X∣1/3: take a Sidon subset SSS of maximum size, and note that every xxx outside it satisfies x=c+d−bx = c + d - bx=c+d−b or x=(c+d)/2x = (c+d)/2x=(c+d)/2 for elements of SSS, so ∣X∣≤3∣S∣3|X| \le 3|S|^3∣X∣≤3∣S∣3.

That cube root is not a weak first attempt, it is the ceiling for any argument that only counts. An arithmetic progression of length nnn has additive energy of order n3n^3n3, so a probabilistic argument cannot see the difference between it and a generic set. Getting from 1/31/31/3 to 1/21/21/2 requires using the structure of the set, and that is what both published proofs do.

The idea both proofs share

Compress, then pigeonhole against a known Sidon set.

An arbitrary finite set of reals has no arithmetic to work with, so first move it into Z\mathbb{Z}Z: a finite set spans a finite dimensional Q\mathbb{Q}Q-vector space, and a generic rational functional separates its points while preserving every relation a+b=c+da + b = c + da+b=c+d. Then squeeze the resulting integers into an interval of length comparable to their number, keeping a constant fraction of them and keeping the property that a Sidon subset of the image lifts to one of the original. Finally intersect with a translate of the Erdos-Turan Sidon set, which has about N\sqrt{N}N​ elements inside {0,…,N−1}\{0, \dots, N-1\}{0,…,N−1}. A set of size Θ(n)\Theta(n)Θ(n) inside [1,n][1,n][1,n] meets some translate of a Sidon set of size n\sqrt{n}n​ in order n\sqrt{n}n​ points, and that intersection is Sidon.

The two proofs differ only in the compression step, and the mission carries both.

The two routes

The 1975 route compresses in four lemmas driven by a remainder map: choose a modulus qqq dividing no difference, dilate so that the remainders are small, and observe that a small remainder map preserves a+b=c+da + b = c + da+b=c+d. Finding the modulus needs a prime counting bound.

The 2026 route replaces all four with one averaging lemma over a real rotation parameter θ\thetaθ, keeping the elements whose fractional part of amθam\thetaamθ is below 1/21/21/2, where no carry occurs. No prime counting appears anywhere.

Notes on the formalization

Every item is stated in Mathlib primitives alone, so the mission needs no definition items. The Sidon condition, the Erdos-Turan construction, and the reduction relation of the 1975 route are all written out at each use.

The published 2026 proof finishes with Singer's 1938 covering of Z/(q2+q+1)Z\mathbb{Z}/(q^2+q+1)\mathbb{Z}Z/(q2+q+1)Z by q+1q+1q+1 Sidon sets. Mathlib has no perfect difference sets, so the mission uses averaging over translates instead. It does the same job at the same order and gives a worse constant, which costs nothing because the goal asserts only that some c>0c > 0c>0 exists.

Two lemmas of the 1975 paper are deliberately absent. A local formalization of Lemma 2 and Lemma 6 turned out to be false as stated, machine-checked in both cases, so neither is offered here as a milestone. Those are errors in that rendering rather than in the paper, and the 2026 route reaches the goal without either. Lemma 1' is absent for the same practical reason: the 2026 route does not need it.

17 thms7 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations Research·Captain: mikedeng1

Cores of Convex Games: The Core of a Convex Game Is Its Unique von Neumann-Morgenstern Stable SetResearch Paper

Motivation

A cooperative game with transferable utility assigns to every coalition of players the total payoff the coalition can secure on its own. Two solution concepts for such games go back to the foundations of game theory: the core, the set of payoff divisions no coalition can improve upon, and the stable set (von Neumann–Morgenstern solution), a set of divisions that is internally consistent and externally absorbing under the relation of domination. For general games the two concepts behave badly: the core may be empty, stable sets may fail to exist (Lucas 1968), and when they exist there are usually many of them.

Lloyd Shapley's paper Cores of Convex Games (Int. J. Game Theory 1, 1971) isolates a class of games, the convex games (supermodular characteristic functions), on which all of this becomes well behaved. Convex games arise in cost allocation, in bankruptcy and airport problems, in scheduling and sequencing games, and in any setting with increasing returns to cooperation; the supermodular functions behind them are the same objects studied as polymatroid rank functions in combinatorial optimization (Edmonds 1970). For such games the paper shows that the core is nonempty, that its faces fit together in a rigid combinatorial pattern, that its vertices are exactly the marginal-contribution vectors, and that the core is the unique stable set.

Setting

Let N={1,…,n}N=\{1,\dots,n\}N={1,…,n} be a finite set of players. A game is a function vvv from subsets of NNN to the reals with v(∅)=0v(\emptyset)=0v(∅)=0. It is convex if

v(S)+v(T)≤v(S∪T)+v(S∩T)for all S,T⊆N.v(S)+v(T)\le v(S\cup T)+v(S\cap T)\qquad\text{for all } S,T\subseteq N.v(S)+v(T)≤v(S∪T)+v(S∩T)for all S,T⊆N.

A payoff vector is a∈RNa\in\mathbb R^Na∈RN, and a(S)=∑i∈Saia(S)=\sum_{i\in S}a_ia(S)=∑i∈S​ai​. It is feasible if a(N)≤v(N)a(N)\le v(N)a(N)≤v(N). The core CCC is the set of feasible aaa with a(S)≥v(S)a(S)\ge v(S)a(S)≥v(S) for every S⊆NS\subseteq NS⊆N; in particular a(N)=v(N)a(N)=v(N)a(N)=v(N) on CCC.

For a nonempty coalition SSS, the face CSC_SCS​ is the set of core points with a(S)=v(S)a(S)=v(S)a(S)=v(S); by convention C∅=CC_\emptyset=CC∅​=C, and CN=CC_N=CCN​=C. The family {CS}\{C_S\}{CS​} is the core configuration. It is complete if no CSC_SCS​ is empty, and regular if CN≠∅C_N\ne\emptysetCN​=∅ and

CS∩CT⊆CS∪T∩CS∩Tfor all S,T⊆N.C_S\cap C_T\subseteq C_{S\cup T}\cap C_{S\cap T}\qquad\text{for all } S,T\subseteq N.CS​∩CT​⊆CS∪T​∩CS∩T​for all S,T⊆N.

For an ordering ω\omegaω of the players, Sω,kS_{\omega,k}Sω,k​ is the set of the first kkk players, and the marginal vector aωa^\omegaaω pays each player iii its marginal contribution v(Sω,ω(i))−v(Sω,ω(i)−1)v(S_{\omega,\omega(i)})-v(S_{\omega,\omega(i)-1})v(Sω,ω(i)​)−v(Sω,ω(i)−1​).

A payoff vector bbb is dominated by aaa if some nonempty coalition SSS has a(S)≤v(S)a(S)\le v(S)a(S)≤v(S) and ai>bia_i>b_iai​>bi​ for all i∈Si\in Si∈S. A set VVV of feasible vectors is stable if every feasible vector is either a member of VVV or dominated by a member of VVV, but not both.

Formalization targets

Goal: Theorem 8

C is stable, and every stable set V equals C(v convex).C \text{ is stable, and every stable set } V \text{ equals } C \qquad (v \text{ convex}).C is stable, and every stable set V equals C(v convex).

The goal contains both halves of the page's statement: stability of the core, and uniqueness ("the unique von Neumann–Morgenstern solution").

Milestones, in the order the argument uses them

  • Lemma 1 (p. 18) and Lemma 2 (p. 19): for a regular configuration, a point on two nested faces CS∩CTC_S\cap C_TCS​∩CT​ with ∣T∖S∣≥2|T\setminus S|\ge2∣T∖S∣≥2 can be moved to a face CQC_QCQ​ of an intermediate coalition, and a point of CSC_SCS​ to CS∩CS∪{j}C_S\cap C_{S\cup\{j\}}CS​∩CS∪{j}​, keeping its coordinates on SSS.
  • Theorem 2 (p. 18): in a regular configuration CS1∩⋯∩CSm≠∅C_{S_1}\cap\cdots\cap C_{S_m}\ne\emptysetCS1​​∩⋯∩CSm​​=∅ for every strictly increasing chain S1⊂⋯⊂SmS_1\subset\cdots\subset S_mS1​⊂⋯⊂Sm​; in particular a regular configuration is complete.
  • Theorem 4 (p. 21): the core of a convex game is nonempty.
  • Theorem 5 (p. 22): a game is convex if and only if its core configuration is regular.
  • Two claims of §4.3 (p. 24): every stable set contains the core, and no stable set properly includes another.
  • The claim that opens the proof of Theorem 8 (p. 24): in a convex game every feasible vector outside the core is dominated by a core point.

The mission also states Theorem 3 (p. 19), the vertices of a regular core are exactly the marginal vectors aωa^\omegaaω, as a further item that is not on the path to the goal.

Significance

The result. Theorem 8 gives, for a natural and widely occurring class of games, a complete answer to the existence and uniqueness questions for von Neumann–Morgenstern solutions, which are open or negative in general. Theorems 3 and 5 describe the core of a convex game explicitly as the polytope spanned by the n!n!n! marginal vectors, the combinatorial description that underlies later work on the Shapley value, the Weber set, and the polymatroid greedy algorithm. Theorem 5 is the geometric characterization of supermodularity through the face structure of the core.

Formalizing it. All results in this mission are proved in the paper; none has a machine-checked proof on the platform. Theorem 4 is already stated on the platform (as part of a statement that also puts every marginal vector and the Shapley value in the core) and enters the mission as an existing item. The remaining work is a formal development of face configurations of the core, of stable sets and domination, and of the passage from supermodularity to the geometry of the core. The definitions of stable set and domination are general and reusable for any transferable-utility game.

Difficulty

The internal half of stability is immediate from the definitions: a core point cannot be dominated by any vector satisfying a coalition constraint a(S)≤v(S)a(S)\le v(S)a(S)≤v(S). Uniqueness also follows from two short observations. The substance is external stability: every feasible vector outside the core must be dominated by a core point, and the dominating vector has to be produced explicitly. The obvious attempt, raising the payoffs of one violated coalition and leaving the other coordinates of bbb unchanged, does not in general produce a core point, and nothing in the definition of the core alone controls how the core meets the hyperplane of a given coalition; that control is what the face theory of §3 is about. For non-convex games the external half genuinely fails, so no argument that ignores convexity can succeed.

Formalization scope

Players are Fin n (a relabelling of the paper's finite set NNN), a game is f : Finset (Fin n) → ℝ, payoff vectors are Fin n → ℝ, and a(S)a(S)a(S) is ∑ i ∈ S, a i. The existing platform definitions Supermodularity.Cooperative.IsConvexGame (v(∅)=0v(\emptyset)=0v(∅)=0 plus supermodularity on all subsets), Core, InitialCoalition and GreedyPayoff (the marginal vectors, orderings being permutations of Fin n) are reused; the reused Theorem 4 statement is Supermodularity.Cooperative.convex_game_core_and_shapley.

Conventions committed to:

  • Wherever the page says "a game", the hypothesis is exactly v(∅)=0v(\emptyset)=0v(∅)=0; convexity is IsConvexGame.
  • Faces satisfy C∅=CC_\emptyset=CC∅​=C literally: the tightness condition is imposed only for nonempty SSS.
  • Regularity includes CN≠∅C_N\ne\emptysetCN​=∅, as on the page.
  • Lemmas 1–2 and Theorems 2–3 assume a regular configuration, not convexity, as on the page.
  • S⊂⊂TS\subset\subset TS⊂⊂T is S⊊TS\subsetneq TS⊊T with ∣T∣−∣S∣≥2|T|-|S|\ge2∣T∣−∣S∣≥2; Lemma 1's two preassigned elements are distinct.
  • An increasing sequence of m≥1m\ge1m≥1 coalitions is a strictly monotone map from Fin (m + 1).
  • "Vertex" is Set.extremePoints ℝ.
  • Domination requires a nonempty coalition and strict coordinate inequalities; stable sets consist of feasible vectors and the "either … or …, but not both" condition ranges over feasible vectors, following the page rather than the classical imputation-based variant.

A formalization in which the dominating coalition may be empty, in which regularity omits CN≠∅C_N\ne\emptysetCN​=∅, or in which the goal asserts stability without uniqueness does not state the paper's theorem and is ruled out.

Welcome contributions: proofs of the milestones in any order, general lemmas about faces of polytopes cut out by set-function inequalities, and reusable API for domination and stable sets.

Selected references

  • L. S. Shapley, Cores of Convex Games, International Journal of Game Theory 1 (1971), 11–26. https://doi.org/10.1007/BF01753431
  • J. von Neumann and O. Morgenstern, Theory of Games and Economic Behavior, Princeton University Press, 1944.
  • J. Edmonds, Submodular functions, matroids, and certain polyhedra, in Combinatorial Structures and Their Applications, Gordon and Breach, 1970, 69–87. https://doi.org/10.1007/3-540-36478-1_2
  • W. F. Lucas, A game with no solution, Bulletin of the American Mathematical Society 74 (1968), 237–239. https://doi.org/10.1090/S0002-9904-1968-12039-2
  • D. M. Topkis, Supermodularity and Complementarity, Princeton University Press, 1998, §5.2.
20 thms6 active usersReviewed
🏆Completed
Number Theory·Captain: aarontcao

Long-Wagner Conjecture 5.1: cube-free subsets of Z/2^nZ have density at most 5/8Open Problem

Call A⊆Z/2nZA \subseteq \mathbb{Z}/2^n\mathbb{Z}A⊆Z/2nZ cube-free if no triple x,y,zx, y, zx,y,z has all seven of xxx, yyy, zzz, x+yx+yx+y, y+zy+zy+z, z+xz+xz+x, x+y+zx+y+zx+y+z inside AAA. The triple is unconstrained, so a degenerate one counts. Write f(n)f(n)f(n) for the largest size of a cube-free subset.

The conjecture. f(n)≤582nf(n) \le \frac{5}{8} 2^nf(n)≤85​2n for every nnn.

This is Conjecture 5.1 of Jason Long and Adam Zsolt Wagner, The largest projective cube-free subsets of Z2n\mathbb{Z}_{2^n}Z2n​, arXiv:1810.01225. It has been open since October 2018, and a 2026 journal paper still names it as conjectured: Yuchen Meng, On Cube-Free Problems, Electron. J. Combin. 33(1) (2026) #P1.16.

The constant is attained

The bound is sharp, and the extremal set is explicit: A={v:v mod 8∈{1,3,4,5,7}}A = \{v : v \bmod 8 \in \{1,3,4,5,7\}\}A={v:vmod8∈{1,3,4,5,7}}, the odd residues together with those congruent to 4 mod 8. Its size is 2n−1+2n−3=582n2^{n-1} + 2^{n-3} = \frac{5}{8} 2^n2n−1+2n−3=85​2n. In the layer language of Long and Wagner this is C3=L1∪L3C_3 = L_1 \cup L_3C3​=L1​∪L3​.

What is known

The conjecture holds for unions of layers. That is Long-Wagner Theorem 1.10 at d=3d = 3d=3, and it is the largest class on which the conjectured constant is proved.

For arbitrary sets the best published unconditional bound is f(n)<232nf(n) < \frac{2}{3} 2^nf(n)<32​2n. Meng calls this bound "quite trivial" and gives it in one paragraph for every cyclic group, so it should not be read as progress toward 5/85/85/8. The residual gap is exactly 23−58=124\frac{2}{3} - \frac{5}{8} = \frac{1}{24}32​−85​=241​, that is 2n/242^n/242n/24 elements.

Small values are f(1)=1f(1) = 1f(1)=1, f(2)=2f(2) = 2f(2)=2, f(3)=5f(3) = 5f(3)=5, f(4)=10f(4) = 10f(4)=10, f(5)=20f(5) = 20f(5)=20, f(6)=40f(6) = 40f(6)=40, f(7)=80f(7) = 80f(7)=80, matching 2n−1+2n−32^{n-1} + 2^{n-3}2n−1+2n−3 from n=3n = 3n=3 on.

State those values honestly. They come from solver searches, Gurobi in Long and Wagner for n≤7n \le 7n≤7 and an independent SAT reproduction. The SAT half that matters, the unsatisfiability of "a cube-free set of size 81 exists at n=7n = 7n=7", is a solver claim with no proof certificate checked and no kernel check behind it. The witness half is verified: a set of exactly 80 elements was produced and re-checked cube-free. So f(7)≥80f(7) \ge 80f(7)≥80 is solid and f(7)≤80f(7) \le 80f(7)≤80 is not certified. Nothing in this mission rests on either.

What the items are

The goal item is the conjecture itself, for n≥4n \ge 4n≥4, and it is open. Every other item is a milestone that is proved mathematics, and the two closed instances n=4n = 4n=4 and n=5n = 5n=5 are stated separately because they are the only cases of the goal that a proof assistant has actually settled here.

The chain runs: the base case mod 8 by exhaustion, monotonicity under subsets, the bridge between the membership form and the Finset form of the forbidden configuration, sharpness, the odd-residue tight case, the layer-union theorem, the two-thirds bound, and then n=4n = 4n=4 and n=5n = 5n=5.

Notes on the formalization

Six definitions live in one definition item, Def_Z2nCubeFreeLayers: HasCube, CubeFree, config, ConfigFree, layerIdx and IsLayerUnion. layerIdx is written through the 2-adic valuation rather than through a congruence, because the congruence form leaves 000 in no layer at all and needs the last layer special-cased.

CubeFree and ConfigFree are two encodings of the same condition and they are not definitionally equal, because config collapses duplicates on a degenerate triple. Their equivalence is a milestone rather than an assumption.

18 thms6 active usersReviewed
🏆Completed
Graph TheoryTheoretical 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
Graph Theory·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
Discrete Geometry·Captain: mysticflounder

Erdős Problems 97 and 96: Convex Point Sets and Unit DistancesOpen Problem

Closed — negative resolution of Erdős Problems 96 and 97

Adam McKenna closed this mission on 13 September 2026 following Unit distances in convex polygons, by Liam Kruer, Jensen Kohlmeyer, and Liam Price. Their construction gives strictly convex point sets with Ω(n log log n) unit-distance pairs and arbitrarily large minimum unit-distance degree, answering both questions and the general fixed-k version of Problem 97 negatively.

Paper and complete Lean source. All credit for the counterexample and its formalization belongs to those authors. Adam McKenna prepared the Prove2Me adapters.

Do not start further proof attempts or solver runs for the affirmative conjectures. Existing statements, conditional lemmas, partial proofs, and milestones remain as historical work. The owner has authorized closure assuming the external result is correct; individual theorem pages report Prove2Me verification status.


Historical mission description

Motivation

The mission is to prove the combined open goal

Problem 97  ∧  Problem 96\text{Problem 97} \;\land\; \text{Problem 96}Problem 97∧Problem 96

for finite point sets in strictly convex position in the Euclidean plane.

Why Problems 97 and 96 belong together

Problem 97 gives the local step needed for Problem 96. Assume Problem 97. Every nonempty convex-independent finite set then has a vertex with at most three neighbors at each positive radius, in particular at radius 111. Delete that vertex and preserve convex independence. Apply the same step to every subset created by deletion until no points remain. Charge each unordered unit-distance pair to the first endpoint deleted. Each deleted vertex receives at most three charges, so an nnn-point set determines at most 3n3n3n unordered unit-distance pairs. This gives the Problem 96 bound and therefore O(n)O(n)O(n). The package uses this one-way dependency; it does not seek a reverse implication.

Setting

Let A⊂R2A\subset\mathbb R^2A⊂R2 be finite. Strict convex position means that every point of AAA is an extreme point of the convex hull of AAA. For p∈Ap\in Ap∈A, the pinned multiplicity at radius r>0r>0r>0 counts points q∈Aq\in Aq∈A with ∥p−q∥=r\lVert p-q\rVert=r∥p−q∥=r. Problem 97 asks for a point where no radius has four such other points. Problem 96 counts unordered pairs at distance 111, then takes the supremum over convex-independent nnn-point sets.

The historical progression is part of the setting. Erdős’s 1946 paper posed an earlier three-neighbor version. His 1987 account reports Danzer’s convex nonagon in which every vertex has three equidistant witnesses, and asks about four witnesses. Fishburn and Reeds’s 1992 work gives a 20-vertex convex configuration with the same unit distance at every vertex, placing the local question beside the unit-distance problem.

Target

The Problem 97 target is the canonical statement that every nonempty finite convex-independent AAA has no four-equidistant-point property:

∀A,A≠∅  →  ConvexIndep⁡(A)  →  ¬HasNEquidistantProperty⁡(4,A).\forall A,\quad A\ne\varnothing\;\to\;\operatorname{ConvexIndep}(A) \;\to\;\neg\operatorname{HasNEquidistantProperty}(4,A).∀A,A=∅→ConvexIndep(A)→¬HasNEquidistantProperty(4,A).

The Problem 96 target is the canonical asymptotic statement

Uc(n)=O(n),U_c(n)=O(n),Uc​(n)=O(n),

where Uc(n)U_c(n)Uc​(n) is the supremum of the unordered unit-distance counts determined by convex-independent nnn-point sets. The bound is asymptotic; the Problem 97 route would give the stronger explicit bound Uc(n)≤3nU_c(n)\le3nUc​(n)≤3n for every natural number nnn.

Significance

The package records a formal proof route joining a pinned geometric obstruction to a global extremal bound. A successful Problem 97 proof would immediately settle Problem 96 with the explicit constant 333, while preserving the combinatorial meaning of the count. It also separates the historical three-neighbor constructions from the still-open four-neighbor assertion.

Difficulty

The source proof reduces Problem 97 to strong induction on ∣A∣|A|∣A∣. Its counting engine follows Dumitrescu's 2006 isosceles-count method, with the cap-witness refinements used in the source attributed to Nivasch--Pach--Pinchasi--Zerbib (2013). This engine forces every counterexample to have at least nine points; a finite geometric analysis excludes exactly nine points; and the remaining step must produce a removable vertex for every larger minimal counterexample. The removable-vertex statement carries the induction hypothesis that every strictly smaller nonempty convex 4-equidistant set is contradictory. That large-cardinality geometric step remains open, so both headline targets remain open. Finite computational certificates can support local cases but do not replace the universal geometric statement.

Counterexample routes

Problem 97 is open, so the mission also records the parallel negative route. The source formalization calls a nonempty convex-independent finite set with the four-equidistant property a Problem97.IsCounterexample. Constructing one such set would refute Problem 97 and therefore refute the mission's affirmative conjunction, regardless of whether Problem 96 remains true. The counterexample milestone keeps this resolution path visible beside the nonexistence proof. A successful witness must use exact coordinates or exact algebraic data from which Lean verifies both strict convex position and the four-equidistant property; a numerical approximation or a realizable incidence pattern alone is insufficient.

Problem 96 has its own negative route. Because its claim is asymptotic, one finite convex configuration cannot refute it. A counterexample must instead give convex-independent point sets at arbitrarily large cardinalities whose unit-distance counts exceed every proposed linear constant. The mission tracks this superlinear-family statement separately, together with a reduction from it to the exact negation of Problem 96. This keeps both possible outcomes visible: a direct or Problem-97-derived linear upper bound, and an explicit family proving that no such bound exists.

Formalization scope

The canonical source is pinned at commit 757d852766f377f7c1a0ffeeef6d3526bc0cb7a4. It contains the formal source statements for Problem 97 and Problem 96. The source repository reports closed proofs of the conditional bridge to the 3n3n3n bound (conditional three-times bound), the ∣A∣≥9|A|\ge9∣A∣≥9 counting milestone (nine-point counting bound), and the exact nine-point exclusion (exact nine-point exclusion theorem). The remaining large-cardinality milestone is the removable-vertex step, with its minimality hypothesis retained. The current platform mission contains accepted transfers of the counting argument, the conditional bridge, and the exact nine-point exclusion, while the removable-vertex step remains open. Its definitions make convex independence and the positive-radius condition explicit; no theorem is assumed inside a definition. Singletons and two-point sets are included in Problem 97, while Problem 96's counting definitions also include the empty set. The source repository uses Lean v4.27.0; these mission statements target the platform's v4.33.1. Source-proof transfer and revalidation remain separate work. The Lean declarations and proofs are this project's own formalization. The Dumitrescu and Nivasch--Pach--Pinchasi--Zerbib citations record mathematical provenance; they do not indicate that a paper proof was imported or machine-checked directly.

These source results establish the intended dependency graph: the P97 universal root feeds low-unit-degree extraction, strong induction, and then the P96 supremum bound. The platform mission records those contracts and milestones; it does not claim to have transplanted their proof bodies. The milestones include the two canonical roots, their conditional bridge, the |A| ≥ 9 count, the n = 9 exclusion, the |A| > 9 removable-vertex step, the documented Danzer nine-point three-neighbor example, the parallel goal of constructing a Problem 97 counterexample, and the superlinear-family route to a counterexample to Problem 96.

References

  • Erdős, On Sets of Distances of n Points (1946), DOI.
  • Erdős, Some Combinatorial and Metric Problems in Geometry (1987), scan.
  • Fishburn–Reeds, Unit Distances Between Vertices of a Convex Polygon (1992), publisher record.
  • Dumitrescu, On Distinct Distances from a Vertex of a Convex Polygon (2006), Springer record; provenance for the source counting method.
  • Nivasch–Pach–Pinchasi–Zerbib, The Number of Distinct Distances from a Vertex of a Convex Polygon (2013), arXiv:1207.1266; provenance for the cap-witness refinements used by the source formalization.
80 thms6 active usersReviewed
🏆Completed
Number Theory·Captain: ShouqiaoWang

Erdős Problem 390: Exact Second-Order AsymptoticResearch Paper

Determine the exact second-order term in the least possible largest factor in a factorization of n!n!n! into distinct integers exceeding nnn, with the proposed rational constant 4029639598/259700381854029639598/259700381854029639598/25970038185.

97 thms6 active usersReviewed
🏆Completed
Graph TheoryTheoretical 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 thms5 active usersReviewed
🏆Completed
Graph Theory·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
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

Competitive Paging Algorithms IV: An Algorithm Competitive against Several Others Exists iff the Reciprocal Ratios Sum to at Most 1Research Paper

Motivation

Paging is the problem of managing a fast memory that holds kkk pages out of nnn: when a requested page is not in fast memory (a page fault), some resident page must be evicted, and the cost of an algorithm is its number of faults. Practitioners have many eviction rules. Least-recently-used (LRU) performs well on real workloads but can be kkk times worse than the optimal off-line schedule; the randomized marking algorithm of the same paper is 2Hk2H_k2Hk​-competitive and so has better worst-case guarantees. Fiat, Karp, Luby, McGeoch, Sleator and Young asked in 1991 whether one on-line algorithm can combine the advantages of several given ones, and answered the question exactly: the attainable combinations of ratios are characterized by one inequality (arXiv:cs/0205038, §6).

The question of combining on-line algorithms has since become a theme of its own: combining heuristics with worst-case-safe algorithms, and, more recently, combining machine-learned predictions with robust fallbacks, both ask for the same kind of guarantee against several reference algorithms at once.

Setting

A type (k,n)(k,n)(k,n) consists of kkk servers and a finite set MMM of nnn vertices with the uniform metric: two distinct vertices are at distance 111. This is paging: vertices are pages, the vertices covered by servers are the pages in fast memory, and a server move is a page fault.

A deterministic on-line algorithm AAA of type (k,n)(k,n)(k,n) has an initial configuration of its kkk servers and, after each request r∈Mr\in Mr∈M, moves servers so that some server covers rrr; its configuration after a request sequence depends only on that sequence. Its cost CA(σ)C_A(\sigma)CA​(σ) on a request sequence σ\sigmaσ is the total distance its servers travel, i.e. the number of server moves.

For algorithms AAA and BBB of the same type and a constant ccc, AAA is ccc-competitive against BBB if there is a constant aaa such that for every request sequence σ\sigmaσ

CA(σ)≤c⋅CB(σ)+a.C_A(\sigma)\le c\cdot C_B(\sigma)+a .CA​(σ)≤c⋅CB​(σ)+a.

A sequence c∗=(c(1),…,c(m))c^*=(c(1),\dots,c(m))c∗=(c(1),…,c(m)) of positive reals is realizable if for every type (k,n)(k,n)(k,n) and every mmm deterministic on-line algorithms B(1),…,B(m)B(1),\dots,B(m)B(1),…,B(m) of that type there is a deterministic on-line algorithm AAA of the same type that is c(i)c(i)c(i)-competitive against B(i)B(i)B(i) for every iii.

Formalization targets

Goal: Theorem 6

For m≥1m\ge1m≥1 and positive reals c(1),…,c(m)c(1),\dots,c(m)c(1),…,c(m),

c∗ is realizable  ⟺  ∑1≤i≤m1c(i)≤1.c^*\ \text{is realizable}\iff \sum_{1\le i\le m}\frac1{c(i)}\le 1 .c∗ is realizable⟺1≤i≤m∑​c(i)1​≤1.

Milestones

In the order of the paper's proof:

  1. Punishments are paid for. If AAA punishes BBB at a time step (an AAA-interval on a vertex vvv ends at that step and contains the end of a BBB-interval on vvv that began no later), then BBB has moved a server; the number of such steps is at most CB(σ)C_B(\sigma)CB​(σ).
  2. A fault leaves room to punish. If ∣SA∣=k|S_A|=k∣SA​∣=k, ∣SB∣≤k|S_B|\le k∣SB​∣≤k, x∈SBx\in S_Bx∈SB​ and x∉SAx\notin S_Ax∈/SA​, then some u∈SAu\in S_Au∈SA​ is not in SBS_BSB​.
  3. The greedy quota claim. If ∑i1/c(i)≤1\sum_i 1/c(i)\le 1∑i​1/c(i)≤1 and each unit of cost punishes the B(i)B(i)B(i) minimizing c(i)(PUN(i)+1)c(i)(\mathrm{PUN}(i)+1)c(i)(PUN(i)+1) (other algorithms may be punished incidentally), then after cost rrr every B(i)B(i)B(i) has been punished at least ⌊r/c(i)⌋\lfloor r/c(i)\rfloor⌊r/c(i)⌋ times.
  4. Shuttle algorithms. With 2m−12m-12m−1 servers on 2m2m2m vertices there are mmm algorithms, each keeping all vertices outside its own pair covered, no two of which move at the same step; in particular their total cost on any σ\sigmaσ is at most ∣σ∣|\sigma|∣σ∣.
  5. A forcing adversary. With 2m−12m-12m−1 servers on 2m2m2m vertices every algorithm can be forced to move at each of NNN steps, so CA(τ(N))≥NC_A(\tau(N))\ge NCA​(τ(N))≥N.

Significance

The result. Theorem 6 is an exact characterization, not a bound: the region of simultaneously attainable ratios against arbitrary deterministic paging algorithms is {c:∑1/c(i)≤1}\{c:\sum 1/c(i)\le 1\}{c:∑1/c(i)≤1}. For example, any two paging algorithms can be combined into one that is 222-competitive against each, and no better symmetric pair is possible in general. Combined with Theorem 7 of the same paper (not part of this mission), the same region is attainable against randomized algorithms, which is how LRU's practical behaviour and the marking algorithm's 2Hk2H_k2Hk​ worst-case guarantee can be obtained within constant factors by one algorithm.

Formalizing it. The theorem has been proved since 1991; no machine-checked proof is known. A formal proof produces a reusable notion of competitiveness of one on-line algorithm against another, built on the published KServer_model definitions, and a formal account of the scheduling fact at the core of the sufficiency proof.

Difficulty

Sufficiency looks like an averaging argument, but the combined algorithm cannot simulate the B(i)B(i)B(i) and follow one of them: switching between their configurations costs up to kkk per switch, which no additive constant absorbs. The accounting has to charge each of AAA's faults to a specific move of a specific B(i)B(i)B(i), and the charge must be injective; the paper's claim that CB(σ)C_B(\sigma)CB​(σ) is at least the number of punishments is where this happens, and it depends on how server intervals are matched. The allocation of faults to algorithms is then a deadline-scheduling problem whose feasibility is exactly ∑1/c(i)≤1\sum 1/c(i)\le 1∑1/c(i)≤1, and the floor functions make the counting delicate at the boundary. The paper's own definition of punishment only counts intervals that start with a move, so the first kkk faults of AAA (servers on their initial vertices) need separate treatment; they are absorbed by the additive constant.

Necessity needs the right family of hard instances: the mmm algorithms must never move at the same step, which pins the type to (2m−1,2m)(2m-1,2m)(2m−1,2m).

Formalization scope

The Lean development works in the namespace CompetitivePaging.Combining and imports the published KServer_model definitions: KServer.OnlineAlgorithm k M (a configuration map from request prefixes to Fin k → M with a serving condition) and OnlineAlgorithm.cost. Committed conventions:

  • a type (k,n)(k,n)(k,n) is any k : ℕ and any finite M : Type with a metric in which distinct points are at distance 111; realizability quantifies over all of them, never over one fixed type;
  • servers are labelled; each algorithm has its own initial configuration, and the additive constant aaa is chosen before the request sequence;
  • c(i)>0c(i)>0c(i)>0 and m≥1m\ge1m≥1 are hypotheses of the goal, as in the paper; without positivity, 1/0=01/0=01/0=0 in Lean would make a zero ratio free;
  • time ttt is the step processing the ttt-th request; the paper's PUN\mathrm{PUN}PUN counts time steps.

Trivializing encodings are ruled out: realizability is not stated for a single fixed type, the metric is not the metric of Fin n, and the competitive constant is not allowed to depend on the request sequence.

A complete proof needs the construction of the punishing algorithm as a KServer.OnlineAlgorithm (a lazy, injective algorithm whose moves depend on the prefix and on the B(i)B(i)B(i)'s configurations), the injective charging argument, the scheduling lemma, and the explicit shuttle algorithms. The scheduling lemma and the charging lemma are independent of paging and reusable. Proofs of any milestone are welcome, as are alternative statements of the sufficiency construction.

Selected references

  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive Paging Algorithms, J. Algorithms 12(4):685–699, 1991. doi:10.1016/0196-6774(91)90041-V; preprint arXiv:cs/0205038.
  • D. D. Sleator, R. E. Tarjan, Amortized efficiency of list update and paging rules, Comm. ACM 28(2):202–208, 1985. doi:10.1145/2786.2793
  • M. S. Manasse, L. A. McGeoch, D. D. Sleator, Competitive algorithms for server problems, J. Algorithms 11(2):208–230, 1990. doi:10.1016/0196-6774(90)90003-W
9 thms5 active usersReviewed
🏆Completed
Graph TheoryLinear 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
🏆Completed
Complexity TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

Scheduling Subject to Resource Constraints: Classification and Complexity I: Unit-Time Chains on Two Identical Machines with One Unit Resource Are Strongly NP-hardResearch Paper

Resource constraints and the easy/hard borderline in scheduling

Machine scheduling asks how to assign jobs to machines over time so that a criterion such as the makespan Cmax⁡C_{\max}Cmax​, the time at which the last job completes, is as small as possible. In practice jobs also compete for scarce resources beyond the machines themselves: tools, operators, memory, power. Błażewicz, Lenstra and Rinnooy Kan (DAM 1983) extended the standard three-field classification α∣β∣γ\alpha\mid\beta\mid\gammaα∣β∣γ of Graham, Lawler, Lenstra and Rinnooy Kan (1979) by a resource field resλσρres\lambda\sigma\rhoresλσρ. They then settled the complexity of every problem with parallel identical or uniform machines, unit-time jobs, precedence constraints and the Cmax⁡C_{\max}Cmax​ criterion. Their Fig. 2 separates the maximal polynomially solvable problems from the minimal NP-hard ones, and it has been the reference map for resource-constrained scheduling since.

Brief timeline of the problems involved:

  • 1975. Garey and Johnson (SIAM J. Comput. 4) show that P2∣res⋯ ,pj=1∣Cmax⁡P2\mid res\cdots, p_j=1\mid C_{\max}P2∣res⋯,pj​=1∣Cmax​ is solvable in polynomial time via matchings, and that P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ and P2∣res1⋅⋅,tree,pj=1∣Cmax⁡P2\mid res1\cdot\cdot, tree, p_j=1\mid C_{\max}P2∣res1⋅⋅,tree,pj​=1∣Cmax​ are NP-hard in the strong sense, by reduction from 3-PARTITION.
  • 1976. Ullman (Complexity of sequencing problems, in Coffman, ed., Computer & Job/Shop Scheduling Theory, Wiley) gives strong NP-hardness of P2∣res111,prec,pj=1∣Cmax⁡P2\mid res111, prec, p_j=1\mid C_{\max}P2∣res111,prec,pj​=1∣Cmax​ under arbitrary precedence constraints.
  • 1983. Błażewicz, Lenstra and Rinnooy Kan prove Theorem 7: chains suffice. Two identical machines, one resource of size one, requirements in {0,1}\{0,1\}{0,1} and chain-like precedence already give a strongly NP-hard problem. The result dominates both earlier two-machine results.

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​. Every job has processing time 111 on every machine, each machine handles at most one job at a time, and jobs are not preempted. There are lll resources; resource RhR_hRh​ has a positive integer size shs_hsh​, the amount available at any time, and job JjJ_jJj​ has a nonnegative integer requirement rhjr_{hj}rhj​, the amount it holds throughout its execution. A directed acyclic graph HHH on the jobs gives the precedence constraints: if HHH has a path from jjj to kkk (Jj→JkJ_j\to J_kJj​→Jk​), then JjJ_jJj​ must complete before JkJ_kJk​ starts. The precedence is chain-like when every vertex of HHH has indegree and outdegree at most one.

A schedule gives each job a machine and a real start time SjS_jSj​; the job occupies [Sj,Sj+1)[S_j, S_j+1)[Sj​,Sj​+1) and completes at Cj=Sj+1C_j = S_j+1Cj​=Sj​+1. It is feasible if jobs on one machine do not overlap, precedence is respected, and at every time ttt the jobs running at ttt require at most shs_hsh​ of each resource RhR_hRh​. The makespan is Cmax⁡=max⁡jCjC_{\max} = \max_j C_jCmax​=maxj​Cj​.

The problem P2∣res111,chain,pj=1∣Cmax⁡P2\mid res111, chain, p_j=1\mid C_{\max}P2∣res111,chain,pj​=1∣Cmax​ restricts this to m=2m=2m=2, one resource (λ=1\lambda=1λ=1) of size 111 (σ=1\sigma=1σ=1), every requirement at most 111 (ρ=1\rho=1ρ=1), and chain-like precedence. The problem P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ has m=3m=3m=3, one resource of arbitrary size and requirements, and no precedence.

3-PARTITION: given ttt, a positive integer bbb and positive integers a1,…,a3ta_1,\dots,a_{3t}a1​,…,a3t​ with ∑jaj=tb\sum_j a_j = tb∑j​aj​=tb and 14b<aj<12b\tfrac14 b<a_j<\tfrac12 b41​b<aj​<21​b, can {1,…,3t}\{1,\dots,3t\}{1,…,3t} be split into ttt disjoint 3-element sets SiS_iSi​ with ∑j∈Siaj=b\sum_{j\in S_i}a_j=b∑j∈Si​​aj​=b?

A problem is NP-hard in the strong sense if it remains NP-hard when every number of the instance is written in unary.

Formalization targets

Goal: Theorem 7

3-PARTITION is NP-hard in the strong sense  ⟹  P2∣res111, chain, pj=1∣Cmax⁡ is NP-hard in the strong sense.\text{3-PARTITION is NP-hard in the strong sense} \;\Longrightarrow\; P2\mid res111,\ chain,\ p_j=1\mid C_{\max}\ \text{is NP-hard in the strong sense.}3-PARTITION is NP-hard in the strong sense⟹P2∣res111, chain, pj​=1∣Cmax​ is NP-hard in the strong sense.

The hypothesis is Garey and Johnson's theorem on 3-PARTITION, which the paper cites and does not prove. The conclusion concerns the decision version: given an instance and y∈Ny\in\mathbb Ny∈N, is there a feasible schedule with Cmax⁡≤yC_{\max}\le yCmax​≤y?

Milestones

  1. Proof of Theorem 4, the saturation equivalence. For positive bbb, aja_jaj​ with ∑jaj=tb\sum_j a_j=tb∑j​aj​=tb, the P3∣res1⋅⋅P3\mid res1\cdot\cdotP3∣res1⋅⋅ instance with 3t3t3t unit jobs, resource size bbb and requirements aja_jaj​ has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t iff the 3-PARTITION instance has a solution.
  2. Theorem 4 (Garey and Johnson). Under the same hypothesis as the goal, P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ is NP-hard in the strong sense.
  3. Proof of Theorem 7, "if". A 3-PARTITION solution yields a feasible schedule of the constructed two-machine instance with Cmax⁡=2tbC_{\max}=2tbCmax​=2tb.
  4. Proof of Theorem 7, "only if". A feasible schedule of the constructed instance with Cmax⁡≤2tbC_{\max}\le 2tbCmax​≤2tb yields a 3-PARTITION solution.

Significance

Theorem 7 is the sharpest hardness result of the paper's classification. Without resources, two-machine unit-time scheduling with arbitrary precedence is polynomial (Coffman and Graham, Acta Inform. 1972); without precedence, it is polynomial under arbitrary resources (Theorem 1 of the paper). The theorem shows that combining the weakest nontrivial versions of both constraints, chains and one unit resource, already crosses the borderline. The paper's §4.1 extends the same reduction to the ∑Cj\sum C_j∑Cj​ and Lmax⁡L_{\max}Lmax​ criteria.

On the formal side, the mission provides a machine-checked model of resource-constrained scheduling with real start times, a definition of NP-hardness in the strong sense on top of the platform's Turing-machine formalization of P\mathrm PP and NP\mathrm{NP}NP, and 3-PARTITION as a reusable source problem. As far as the platform's corpus shows, none of Theorems 4 and 7, 3-PARTITION, or strong NP-hardness has been formalized before. Both theorems are proved in the literature; what remains is to formalize the reductions and their polynomial running time.

Difficulty

The combinatorial heart is the "only if" direction: a schedule of length 2tb2tb2tb must be shown to be rigid. Start times are arbitrary reals, so the first obstacle is to show that both machines are busy throughout [0,2tb)[0,2tb)[0,2tb), that the chain LLL forces unit spacing, and that the primed jobs of the chains Kj′K'_jKj′​ can only run in the intervals the chain LLL leaves free of the resource. Only after this is established can the index sets SiS_iSi​ be read off. Arguing on integer time slots from the start is not enough: the model allows fractional start times, and ruling them out is part of the proof.

The second obstacle is the complexity layer. NP-hardness is stated with respect to polynomial-time many-one reductions computed by one-tape Turing machines. The reduction from 3-PARTITION therefore has to be implemented and its running time bounded on unary codes. The constructed instance has 4tb4tb4tb jobs, which is polynomial in the unary length of the 3-PARTITION instance; this is exactly why the reduction proves hardness in the strong sense.

Formalization scope

  • Model. Jobs are Fin n and machines Fin m, 0-based. Only identical machines with unit processing times are modelled. Start times are real, execution intervals are half-open, and the resource constraint is imposed at every real time. Precedence is the transitive closure of the arc list of HHH. Cmax⁡=0C_{\max}=0Cmax​=0 for an empty instance.
  • Decision version. Thresholds yyy are natural numbers; this narrower class makes the hardness statement stronger.
  • Encoding. An instance is described by its list of numbers (n,m,ln,m,ln,m,l, the sizes, the requirements row by row, the number of arcs and the arcs, then yyy). The unary language is the set of unary codes of yes-instances over a two-letter alphabet. No pairing function is used. The class conditions (two machines, one unit resource, requirements at most one, chain-like acyclic HHH) are part of the yes-predicate.
  • Strong sense. Strong NP-hardness is NP-hardness of the unary language. This is equivalent to Garey and Johnson's definition, which bounds the largest number by a polynomial in the instance length.
  • Cited hypothesis. The goal and Theorem 4 assume strong NP-hardness of 3-PARTITION (with 14b<aj<12b\tfrac14 b<a_j<\tfrac12 b41​b<aj​<21​b) and nothing else. Stating the goal as a bare reduction between the two languages, or adding P≠NP\mathrm P\ne\mathrm{NP}P=NP, would not be Theorem 7.
  • Constructions. The two scheduling instances built from a 3-PARTITION instance are explicit definitions following the page, not arbitrary instances with a property.
  • Reuse. The scheduling model and the strong-NP-hardness layer are shared with the other missions of this series; 3-PARTITION serves any strong NP-hardness proof by number partitioning.

Welcome contributions: proofs of the four milestones; a formalized polynomial-time implementation of the reduction on unary codes; general lemmas about composing polynomial-time reductions on the one-tape machine model.

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 J. Comput. 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, 1979.
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Ann. Discrete Math. 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • J. D. Ullman, Complexity of sequencing problems, in: E. G. Coffman, Jr., ed., Computer & Job/Shop Scheduling Theory, Wiley, 1976, 139–164.
  • E. G. Coffman, Jr., R. L. Graham, Optimal scheduling for two-processor systems, Acta Informatica 1 (1972) 200–213. https://doi.org/10.1007/BF00288685
  • S. Cook, The P versus NP problem, Clay Mathematics Institute official problem description.
11 thms5 active usersReviewed
🏆Completed
Complexity Theory·Captain: marwahaha

Block Sensitivity vs. Sensitivity: Reach 2.05Open Problem

Establish a separation exponent of at least 2.052.052.05 between block sensitivity and sensitivity of total Boolean functions.

The target is an infinite family fk:{0,1}nk→{0,1}f_k:\{0,1\}^{n_k}\to\{0,1\}fk​:{0,1}nk​→{0,1} such that

s(fk)→∞,bs(fk)≥s(fk)2.05−o(1)as k→∞.s(f_k)\to\infty,\qquad bs(f_k)\ge s(f_k)^{2.05-o(1)}\quad\text{as }k\to\infty.s(fk​)→∞,bs(fk​)≥s(fk​)2.05−o(1)as k→∞.

Here o(1)o(1)o(1) tends to zero with kkk. Sensitivity is the maximum number of individual coordinate flips that change the output at a common input; block sensitivity allows pairwise disjoint nonempty sets of coordinates instead. Both maxima range over all inputs.

The Lean goal uses the campaign’s shared definition: for every ε>0\varepsilon>0ε>0 and every integer S≥2S\ge2S≥2, there must be a total Boolean function fff with s(f)≥Ss(f)\ge Ss(f)≥S and bs(f)≥s(f)2.05−εbs(f)\ge s(f)^{2.05-\varepsilon}bs(f)≥s(f)2.05−ε. The decimal 2.052.052.05 means exactly 41/2041/2041/20.

One sufficient route is a finite function fff, an input xxx, and integers A,B>1A,B>1A,B>1 satisfying

s(f)≤A,bs(f,x)≥B,A2.05≤B.s(f)\le A,\qquad bs(f,x)\ge B,\qquad A^{2.05}\le B.s(f)≤A,bs(f,x)≥B,A2.05≤B.

The proved theorem SensitivityCampaign.finite_certificate then gives the asymptotic conclusion. Direct family constructions are equally welcome.

The starting milestone is the campaign’s proved exponent 2.02052.02052.0205. The target 2.052.052.05 is a challenge to prove, not an established campaign record. Further background and construction parameters appear in the campaign proposal.

6 thms4 active usersReviewed
🏆Completed
Complexity Theory·Captain: marwahaha

Block Sensitivity vs. Sensitivity: Reach 2.007Open Problem

Establish the separation exponent α≥2.007\alpha\ge2.007α≥2.007 for total Boolean functions, improving on the depth-20 specialization of OpenAI's construction.

Precisely, for every ε>0\varepsilon>0ε>0 and every integer S≥2S\ge2S≥2, construct a function fff with s(f)≥Ss(f)\ge Ss(f)≥S and bs(f)≥s(f)2.007−εbs(f)\ge s(f)^{2.007-\varepsilon}bs(f)≥s(f)2.007−ε.

A finite certificate is sufficient: supply a function and an input xxx with s(f)≤As(f)\le As(f)≤A, bs(f,x)≥Bbs(f,x)\ge Bbs(f,x)≥B, A,B>1A,B>1A,B>1, and A2.007≤BA^{2.007}\le BA2.007≤B, then apply the shared certificate theorem. Sharper estimates, improved parameters, and different constructions are all eligible; direct family proofs are also welcome.

This is an open research target. The preceding benchmarks are 222 classically and approximately 2.0064708282.0064708282.006470828 from the depth-20 seed. Numerical searches guide discovery but do not replace certified bounds. Intermediate improvements are welcome as separate campaign entries.

38 thms4 active usersReviewed
🏆Completed
Complexity Theory·Captain: marwahaha

Block Sensitivity vs. Sensitivity: Depth-20 ImprovementResearch Paper

Certify α≥2.006470828\alpha\ge2.006470828α≥2.006470828 by specializing OpenAI's general seed theorem to d=20d=20d=20 and M=921M=9^{21}M=921.

The construction gives

s(f)≤A=44M21,bs(f,0)≥B=M2(2M2+1)20.s(f)\le A=44M^{21},\qquad bs(f,0)\ge B=M^2(2M^2+1)^{20}.s(f)≤A=44M21,bs(f,0)≥B=M2(2M2+1)20.

Here bs(f,0)bs(f,0)bs(f,0) is block sensitivity at the all-zero input. The ratio log⁡B/log⁡A\log B/\log AlogB/logA is approximately 2.006470828227379…2.006470828227379\ldots2.006470828227379…. Prove the exact inequality A2.006470828≤BA^{2.006470828}\le BA2.006470828≤B, then reuse the foundation's finite-certificate theorem.

The goal uses the campaign's shared definitions and requires arbitrarily large sensitivity. This is a stronger specialization of the existing construction, not a claim of a new construction or global parameter optimality.

39 thms4 active usersReviewed
🏆Completed
Linear OptimizationOptimal TransportOptimization·Captain: mikedeng1

Computational Optimal Transport V: A Feasible Pair of Kantorovich Potentials Is Optimal or Strictly Improves Along a Direction (1_S, −1_S′)Textbook

Motivation

The discrete optimal transport problem between two histograms is a linear program, and much of the algorithmic theory of optimal transport is the theory of solving that linear program well. Chapter 3 of Peyré and Cuturi's Computational Optimal Transport (Foundations and Trends in Machine Learning, 2019) surveys the classical combinatorial solvers: the network simplex, dual ascent methods, and the auction algorithm. Dual ascent methods work entirely with the dual variables, the Kantorovich potentials (f,g)(f, g)(f,g): they keep a feasible pair of potentials and improve it step by step until it is optimal. The specialisation of this idea to assignment problems is the Hungarian algorithm of Kuhn (1955), and its general form is the primal-dual method for network flow problems presented in Bertsimas and Tsitsiklis, Introduction to Linear Optimization (1997, §7.7), on which §3.6 of the book is modelled.

The mathematical engine of every dual ascent method is one alternative: a feasible pair of potentials is either already optimal, or it can be improved along a direction of a very special, combinatorial form. This mission formalizes that alternative (Proposition 3.6 of the book) together with the results it rests on.

Setting

Fix integers n,m≥1n, m \ge 1n,m≥1, a cost matrix C∈Rn×mC \in \mathbb{R}^{n\times m}C∈Rn×m, and two histograms a∈Σna \in \Sigma_na∈Σn​, b∈Σmb \in \Sigma_mb∈Σm​, where Σn={a∈R+n:∑iai=1}\Sigma_n = \{a \in \mathbb{R}^n_+ : \sum_i a_i = 1\}Σn​={a∈R+n​:∑i​ai​=1} is the probability simplex. Write [ ⁣[n] ⁣]={1,…,n}[\![n]\!] = \{1, \dots, n\}[[n]]={1,…,n} for the row indices and, following the book, [ ⁣[m] ⁣]′={1′,…,m′}[\![m]\!]' = \{1', \dots, m'\}[[m]]′={1′,…,m′} for the column indices.

The couplings between aaa and bbb form the transportation polytope

U(a,b)={P∈R+n×m:P1m=a, PT1n=b},U(a,b) = \{P \in \mathbb{R}^{n\times m}_+ : P\mathbb{1}_m = a,\ P^{\mathsf T}\mathbb{1}_n = b\},U(a,b)={P∈R+n×m​:P1m​=a, PT1n​=b},

and the primal problem is LC(a,b)=min⁡P∈U(a,b)⟨C,P⟩\mathrm{L}_C(a,b) = \min_{P\in U(a,b)} \langle C, P\rangleLC​(a,b)=minP∈U(a,b)​⟨C,P⟩ with ⟨C,P⟩=∑i,jCi,jPi,j\langle C, P\rangle = \sum_{i,j} C_{i,j}P_{i,j}⟨C,P⟩=∑i,j​Ci,j​Pi,j​. A pair of vectors (f,g)∈Rn×Rm(f, g) \in \mathbb{R}^n \times \mathbb{R}^m(f,g)∈Rn×Rm is dual feasible, written (f,g)∈R(C)(f,g) \in R(C)(f,g)∈R(C), when fi+gj≤Ci,jf_i + g_j \le C_{i,j}fi​+gj​≤Ci,j​ for all (i,j)(i, j)(i,j). The dual problem (3.4) is

LC(a,b)=max⁡(f,g)∈R(C)⟨f,a⟩+⟨g,b⟩.\mathrm{L}_C(a,b) = \max_{(f,g)\in R(C)} \langle f, a\rangle + \langle g, b\rangle .LC​(a,b)=(f,g)∈R(C)max​⟨f,a⟩+⟨g,b⟩.

For a feasible pair (f,g)(f, g)(f,g), a pair of indices (i,j′)(i, j')(i,j′) is balanced if fi+gj=Ci,jf_i + g_j = C_{i,j}fi​+gj​=Ci,j​ and inactive if fi+gj<Ci,jf_i + g_j < C_{i,j}fi​+gj​<Ci,j​. A matrix PPP and a pair (f,g)(f, g)(f,g) are complementary if Pi,j>0P_{i,j} > 0Pi,j​>0 implies Ci,j=fi+gjC_{i,j} = f_i + g_jCi,j​=fi​+gj​, that is, PPP is supported on balanced pairs. For S⊂[ ⁣[n] ⁣]S \subset [\![n]\!]S⊂[[n]] the vector 1S∈Rn\mathbb{1}_S \in \mathbb{R}^n1S​∈Rn has ones at the indices in SSS and zeros elsewhere; likewise 1S′∈Rm\mathbb{1}_{S'} \in \mathbb{R}^m1S′​∈Rm for S′⊂[ ⁣[m] ⁣]′S' \subset [\![m]\!]'S′⊂[[m]]′.

Formalization targets

Goal: Proposition 3.6 (p. 416)

For a∈Σna \in \Sigma_na∈Σn​, b∈Σmb \in \Sigma_mb∈Σm​ and (f,g)∈R(C)(f,g) \in R(C)(f,g)∈R(C): either (f,g)(f, g)(f,g) is optimal for (3.4), or there exist S⊂[ ⁣[n] ⁣]S \subset [\![n]\!]S⊂[[n]], S′⊂[ ⁣[m] ⁣]′S' \subset [\![m]\!]'S′⊂[[m]]′ and ε0>0\varepsilon_0 > 0ε0​>0 such that for every 0<ε≤ε00 < \varepsilon \le \varepsilon_00<ε≤ε0​,

(f~,g~)=(f,g)+ε(1S,−1S′)∈R(C)and⟨f~,a⟩+⟨g~,b⟩>⟨f,a⟩+⟨g,b⟩.(\tilde f, \tilde g) = (f, g) + \varepsilon(\mathbb{1}_S, -\mathbb{1}_{S'}) \in R(C) \quad\text{and}\quad \langle \tilde f, a\rangle + \langle \tilde g, b\rangle > \langle f, a\rangle + \langle g, b\rangle .(f~​,g~​)=(f,g)+ε(1S​,−1S′​)∈R(C)and⟨f~​,a⟩+⟨g~​,b⟩>⟨f,a⟩+⟨g,b⟩.

Milestones

  1. Proposition 3.3 (p. 405). If P∈U(a,b)P \in U(a,b)P∈U(a,b) and (f,g)∈R(C)(f,g) \in R(C)(f,g)∈R(C) are complementary, then PPP is primal optimal and (f,g)(f,g)(f,g) is dual optimal.
  2. Proposition 3.5 (p. 415). If every balanced pair (i,j′)(i, j')(i,j′) with i∈Si \in Si∈S has j′∈S′j' \in S'j′∈S′, then (f,g)+ε(1S,−1S′)∈R(C)(f,g) + \varepsilon(\mathbb{1}_S, -\mathbb{1}_{S'}) \in R(C)(f,g)+ε(1S​,−1S′​)∈R(C) for all sufficiently small ε>0\varepsilon > 0ε>0.
  3. Objective change (proof of Proposition 3.6, pp. 416–417). The step changes the dual objective by exactly ε(1STa−1S′Tb)\varepsilon(\mathbb{1}_S^{\mathsf T}a - \mathbb{1}_{S'}^{\mathsf T}b)ε(1ST​a−1S′T​b).
  4. Labeled sets (proof of Proposition 3.6, pp. 416–417). If no coupling in U(a,b)U(a,b)U(a,b) is complementary to (f,g)(f,g)(f,g), then there are S,S′S, S'S,S′ with every balanced pair leaving SSS landing in S′S'S′, and
1STa−1S′Tb>0.\mathbb{1}_S^{\mathsf T}a - \mathbb{1}_{S'}^{\mathsf T}b > 0 .1ST​a−1S′T​b>0.

Significance

The result. Proposition 3.6 says that non-optimality of a feasible dual pair is always witnessed by a direction with entries in {0,±1}\{0, \pm 1\}{0,±1}, determined by two index sets, which keeps the pair feasible for a positive step and strictly improves the objective. This is what makes dual ascent a finite combinatorial method rather than a generic linear-programming iteration: the search for an ascent direction reduces to a maximum-flow computation on the bipartite graph of balanced pairs. Combined with the step length of Proposition 3.5, it is the primal-dual method, which reduces to the Hungarian algorithm on assignment problems. Milestone 1 is the optimality certificate of complementary slackness, used throughout the book's chapter 3.

Formalizing it. These results are classical and proved in the book and in Bertsimas and Tsitsiklis; no machine-checked proof of them is known to exist. The mission produces a formal statement and proof of the optimality-or-ascent alternative for the discrete Kantorovich dual, together with the complementary-slackness certificate in the transport setting. Milestone 4 is a Hall-type statement (a supply–demand theorem on a bipartite graph) whose formal proof is reusable for other matching and transportation results.

Difficulty

Milestones 1–3 are short computations. The content of the goal is in milestone 4. The obvious argument is linear-programming duality: if (f,g)(f,g)(f,g) is not optimal, some feasible direction improves the objective. But LP duality only produces an arbitrary real direction; the claim that a direction of the form (1S,−1S′)(\mathbb{1}_S, -\mathbb{1}_{S'})(1S​,−1S′​) suffices is a combinatorial statement. The book obtains S,S′S, S'S,S′ as the labeled nodes of a maximal flow (Ford–Fulkerson) and argues through max-flow/min-cut; the flow bookkeeping printed on p. 417 is garbled, so the step from "no complementary coupling exists" to "1STa>1S′Tb\mathbb{1}_S^{\mathsf T}a > \mathbb{1}_{S'}^{\mathsf T}b1ST​a>1S′T​b for a set closed under balanced edges" has to be supplied carefully. Separately, connecting "not optimal" to "no complementary coupling exists" needs the converse direction of complementary slackness or strong duality for the transport problem.

Formalization scope

Row and column indices are Fin n and Fin m (0-based); the primed column set [ ⁣[m] ⁣]′[\![m]\!]'[[m]]′ is just Fin m. Cost matrices and couplings are Matrix (Fin n) (Fin m) ℝ; histograms and potentials are real functions on Fin n, Fin m; Σn\Sigma_nΣn​ is Mathlib's stdSimplex ℝ (Fin n). Index sets S,S′S, S'S,S′ are Finsets and 1S\mathbb{1}_S1S​ is indicatorVec S. "Optimal for Problem (3.4)" is the predicate IsDualOptimal: feasible, and no feasible pair has a larger objective (no real infimum or supremum is used, so no junk value can enter). The book's "for a small enough ε>0\varepsilon > 0ε>0" is stated as "for every ε∈(0,ε0]\varepsilon \in (0, \varepsilon_0]ε∈(0,ε0​]", which is equivalent because R(C)R(C)R(C) is convex and the objective is linear.

Standing assumptions, all from the book: (f,g)(f, g)(f,g) is dual feasible (p. 415, "In what follows, (f,g)(f,g)(f,g) is a feasible dual pair in R(C)R(C)R(C)"); a,ba, ba,b are histograms in the simplex (used for the goal and milestone 4; without equal masses the dual is unbounded and optimality fails for every pair). Milestone 4 replaces the book's flow hypothesis "the throughput is strictly smaller than 1" by the equivalent flow-free statement "no coupling in U(a,b)U(a,b)U(a,b) is complementary to (f,g)(f,g)(f,g)"; flows, capacities and the labeling algorithm appear in no statement.

The goal admits a trivializing formalization: allowing an arbitrary direction (u,v)(u, v)(u,v) instead of (1S,−1S′)(\mathbb{1}_S, -\mathbb{1}_{S'})(1S​,−1S′​) turns Proposition 3.6 into plain non-optimality of a linear program. The statement here fixes the direction to (1S,−1S′)(\mathbb{1}_S, -\mathbb{1}_{S'})(1S​,−1S′​) exactly as the book does.

A complete development needs finite LP duality or complementary slackness for the transportation problem, and a max-flow/min-cut or Hall-type theorem on bipartite graphs with vertex capacities. Both are reusable well beyond this mission. Contributions welcome: proofs of the milestones, a self-contained proof of milestone 4 by induction or by a max-flow argument, and alternative proofs of the goal through LP duality plus a vertex argument.

Selected references

  • G. Peyré, M. Cuturi, Computational Optimal Transport, Foundations and Trends in Machine Learning 11(5–6):355–607, 2019. https://doi.org/10.1561/2200000073 (§3.1–3.3, §3.6, pp. 400–417)
  • D. Bertsimas, J. N. Tsitsiklis, Introduction to Linear Optimization, Athena Scientific, 1997, §7.7 (the primal-dual method) and pp. 305–308 (Ford–Fulkerson and the labeling algorithm).
  • H. W. Kuhn, The Hungarian method for the assignment problem, Naval Research Logistics Quarterly 2(1–2):83–97, 1955. https://doi.org/10.1002/nav.3800020109
  • 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
9 thms4 active usersReviewed
🏆Completed
Group TheoryNumber Theory·Captain: mikedeng1

New Proofs of Plünnecke-type Estimates for Product Sets in Groups 1: In an Abelian Group, |A + B| ≤ α|A| Gives One Nonempty X ⊆ A with |X + hB| ≤ αʰ|X| for Every hResearch Paper

Motivation

A basic question of additive combinatorics asks how a set grows under repeated addition. If a finite set AAA in an abelian group satisfies ∣A+B∣≤α∣A∣|A + B| \le \alpha |A|∣A+B∣≤α∣A∣ for some finite set BBB, how large can the iterated sumsets hB=B+⋯+BhB = B + \cdots + BhB=B+⋯+B be? Plünnecke's inequality gives the answer used throughout the subject: some nonempty part of AAA grows by at most a factor α\alphaα per added copy of BBB. Combined with Ruzsa's triangle inequality it becomes the Plünnecke–Ruzsa inequality ∣mB−nB∣≤αm+n∣A∣|mB - nB| \le \alpha^{m+n}|A|∣mB−nB∣≤αm+n∣A∣, which underlies Freĭman-type structure theorems, sum–product estimates and the polynomial Freĭman–Ruzsa programme.

Timeline.

  • 1970: H. Plünnecke proves the inequality for iterated sumsets in abelian groups with a graph-theoretic method (commutative layered graphs and Menger's theorem).
  • 1989: I. Z. Ruzsa simplifies Plünnecke's argument and combines it with the triangle inequality to bound ∣mB−nB∣|mB - nB|∣mB−nB∣.
  • 2011: G. Petridis gives a short elementary proof, without graphs (arXiv:1101.3507, Combinatorica 32, 2012). His argument also shows that one subset X⊆AX \subseteq AX⊆A works for every number of summands hhh, and it extends to triple products CXBCXBCXB in non-abelian groups.

Setting

Let GGG be a group with finite subsets AAA, BBB, CCC. The product set is AB={ab:a∈A, b∈B}AB = \{ab : a \in A,\ b \in B\}AB={ab:a∈A, b∈B}, and ∣⋅∣|\cdot|∣⋅∣ denotes cardinality. In an abelian group, written additively, the sumset is A+B={a+b:a∈A, b∈B}A + B = \{a + b : a \in A,\ b \in B\}A+B={a+b:a∈A, b∈B}, and for h≥0h \ge 0h≥0 the iterated sumset hBhBhB is B+⋯+BB + \cdots + BB+⋯+B with hhh summands, 0B={0}0B = \{0\}0B={0}. It is not the dilate {hb:b∈B}\{hb : b \in B\}{hb:b∈B}.

A nonempty subset X⊆AX \subseteq AX⊆A has minimal growth under BBB if

∣XB∣∣X∣≤∣ZB∣∣Z∣for every nonempty Z⊆A.\frac{|XB|}{|X|} \le \frac{|ZB|}{|Z|} \quad\text{for every nonempty } Z \subseteq A.∣X∣∣XB∣​≤∣Z∣∣ZB∣​for every nonempty Z⊆A.

Such an XXX exists whenever AAA is nonempty, since AAA has finitely many subsets. In Lean, GGG is a Group or AddCommGroup with decidable equality, the sets are Finset G with the pointwise operations of Mathlib (A * B, A + B, h • B), and cardinalities are #A cast to ℝ.

Formalization targets

Goal: Theorem 3.1

For finite AAA, BBB in an abelian group with AAA nonempty and ∣A+B∣≤α∣A∣|A + B| \le \alpha|A|∣A+B∣≤α∣A∣, there is a nonempty X⊆AX \subseteq AX⊆A with

∣X+hB∣≤αh∣X∣for every h≥0.|X + hB| \le \alpha^h |X| \qquad \text{for every } h \ge 0.∣X+hB∣≤αh∣X∣for every h≥0.

The point is the order of quantifiers: one XXX serves all hhh.

Milestones

  1. Proposition 2.1 (p. 5). If ∣XB∣/∣X∣≤∣ZB∣/∣Z∣|XB|/|X| \le |ZB|/|Z|∣XB∣/∣X∣≤∣ZB∣/∣Z∣ for every Z⊆XZ \subseteq XZ⊆X, then ∣CXB∣≤∣XB∣∣X∣ ∣CX∣|CXB| \le \frac{|XB|}{|X|}\,|CX|∣CXB∣≤∣X∣∣XB∣​∣CX∣ for every finite CCC, in any group. This is already on the platform, proved, as FamousTheorems.pluennecke_petridis (additive notation, cross-multiplied), and enters the mission as a reference.
  2. Theorem 1.5 (p. 3). In any group, ∣AB∣≤α∣A∣|AB| \le \alpha|A|∣AB∣≤α∣A∣ with AAA nonempty gives a nonempty X⊆AX \subseteq AX⊆A with
∣CXB∣≤α∣CX∣for every finite C.|CXB| \le \alpha |CX| \quad\text{for every finite } C.∣CXB∣≤α∣CX∣for every finite C.
  1. Inductive step of Theorem 3.1 (p. 7). For X⊆AX \subseteq AX⊆A of minimal growth and ∣A+B∣≤α∣A∣|A+B| \le \alpha|A|∣A+B∣≤α∣A∣:
∣X+hB∣≤α ∣X+(h−1)B∣(h≥1).|X + hB| \le \alpha\,|X + (h-1)B| \quad (h \ge 1).∣X+hB∣≤α∣X+(h−1)B∣(h≥1).

Significance

Theorem 3.1 is the form of Plünnecke's inequality used in applications. Two consequences follow directly. Taking A=BA = BA=B gives ∣hB∣≤αh∣B∣|hB| \le \alpha^h|B|∣hB∣≤αh∣B∣ whenever ∣2B∣≤α∣B∣|2B| \le \alpha|B|∣2B∣≤α∣B∣. Combining one XXX for all hhh with Ruzsa's triangle inequality gives ∣mB−nB∣≤αm+n∣A∣|mB - nB| \le \alpha^{m+n}|A|∣mB−nB∣≤αm+n∣A∣ (Theorem 1.2 of the paper). With a separate XXX for each hhh, the second step would need the same set for mmm and nnn, and that is exactly what Theorem 3.1 provides. Theorem 1.5 is the non-abelian ingredient: missions 2 and 3 of this series apply it to obtain explicit Plünnecke-type bounds for ∣Bh∣|B^h|∣Bh∣ and ∣SBh∣|SB^h|∣SBh∣ in arbitrary groups.

Formalizing it. The results are proved and classical. Mathlib already contains Proposition 2.1 (Finset.pluennecke_petridis_inequality_mul) and the Plünnecke–Ruzsa inequality Finset.pluennecke_ruzsa_inequality_nsmul_sub_nsmul_add. Its proof runs through a private lemma that bounds ∣X+hB∣|X + hB|∣X+hB∣ for one fixed minimiser, but no public statement in Mathlib or on the platform asserts Theorem 3.1 or Theorem 1.5 in the form "there exists a nonempty XXX for all hhh" (respectively "for all CCC"). This mission adds those statements as reusable, citable nodes in the paper's numbering.

Difficulty

A first attempt iterates the hypothesis directly: ∣A+2B∣≤α∣A+B∣|A + 2B| \le \alpha|A + B|∣A+2B∣≤α∣A+B∣ would follow if A+BA + BA+B again grew by a factor α\alphaα under BBB, but nothing controls the growth of A+BA + BA+B. The ratio ∣ZB∣/∣Z∣|ZB|/|Z|∣ZB∣/∣Z∣ is not monotone in ZZZ, and passing to subsets of AAA can only help if the subset is chosen correctly. Choosing a different good subset for each hhh gives Theorem 1.1 but loses the uniformity. The work is to show that one fixed subset controls all the products CXBCXBCXB simultaneously, and in the non-abelian Theorem 1.5 to keep the order of the factors CCC, XXX, BBB throughout.

Formalization scope

  • Theorem 3.1 and the inductive step are stated in an AddCommGroup with DecidableEq. Theorem 1.5 is stated in a Group, with no commutativity, as on the page. The factor order C⋅X⋅BC \cdot X \cdot BC⋅X⋅B is preserved.
  • hBhBhB is h • B (the hhh-fold sumset; 0 • B = {0}), and "for all hhh" ranges over all natural numbers, including the true case h=0h = 0h=0.
  • α\alphaα is an arbitrary real number; cardinalities are cast to ℝ. No hypothesis α≥1\alpha \ge 1α≥1 is added: it follows from the hypotheses once AAA and BBB are nonempty (with BBB empty they force only α≥0\alpha \ge 0α≥0, and every conclusion is still true).
  • Nonempty witness. The page says "there exists X⊆AX \subseteq AX⊆A". Since ∅+hB=∅\emptyset + hB = \emptyset∅+hB=∅, the empty set would satisfy the printed statements trivially, and that formalization is ruled out: the goal and Theorem 1.5 assume AAA nonempty and assert XXX nonempty, matching the paper's choice of XXX as a minimiser of a ratio defined on nonempty sets.
  • Minimal growth in the step milestone is written cross-multiplied, ∣X+B∣ ∣Z∣≤∣Z+B∣ ∣X∣|X+B|\,|Z| \le |Z+B|\,|X|∣X+B∣∣Z∣≤∣Z+B∣∣X∣, over nonempty Z⊆AZ \subseteq AZ⊆A; the step is stated for h≥1h \ge 1h≥1 with natural-number subtraction guarded by that hypothesis.
  • Ruzsa's covering lemma and Ruzsa's triangle inequality, quoted in §4, are in Mathlib (Finset.ruzsa_covering_mul, Finset.ruzsa_triangle_inequality_mul_mulInv_mul) and are not items of this mission.
  • No new definitions are needed. Welcome contributions: proofs of the three draft statements, and, separately, the deduction of Theorem 1.2 from Theorem 3.1.

Selected references

  • G. Petridis, New proofs of Plünnecke-type estimates for product sets in groups, arXiv:1101.3507v3, 2011; Combinatorica 32 (2012). https://arxiv.org/abs/1101.3507 , https://doi.org/10.1007/s00493-012-2818-5
  • H. Plünnecke, Eine zahlentheoretische Anwendung der Graphentheorie, J. Reine Angew. Math. 243 (1970) 171–183. https://doi.org/10.1515/crll.1970.243.171
  • I. Z. Ruzsa, An application of graph theory to additive number theory, Scientia, Ser. A 3 (1989) 97–109 (no online version).
  • T. Tao, V. Vu, Additive Combinatorics, Cambridge University Press, 2006. https://doi.org/10.1017/CBO9780511755149
4 thms4 active usersReviewed
🏆Completed
Algorithmic Game TheoryMechanism Design·Captain: mikedeng1

Strategy-Proofness and Arrow's Conditions 1: A Strict Voting Procedure with at Least Three Possible Outcomes Is Strategy-Proof Iff It Is DictatorialResearch Paper

Motivation

A committee that chooses one alternative by voting would like a voting procedure under which no member can gain by misreporting their preferences. Such a procedure is called strategy-proof. Gibbard (Econometrica 41, 1973) and Satterthwaite (J. Econ. Theory 10, 1975) proved independently that, once at least three outcomes are possible, the only strategy-proof procedures are dictatorial. The result is the starting point of the theory of strategy-proof mechanisms: it is why mechanism design restricts preference domains (single-peaked preferences, quasi-linear utilities) or allows randomization.

This mission formalizes Theorem 1 of Satterthwaite's paper, in the version circulated as Northwestern University CMS-EMS Discussion Paper No. 122 (revised December 12, 1974). Its formulation has one feature that the textbook version lacks: the hypothesis is on the range of the procedure, not on the whole alternative set, and the dictatorship is accordingly relative to that range.

Setting

A committee has a finite set InI_nIn​ of n≥1n\ge 1n≥1 individuals and a finite set SmS_mSm​ of mmm alternatives. In a strict committee every individual's ballot BiB_iBi​ is a strong order on SmS_mSm​, a ranking without ties; x Bˉi yx\,\bar B_i\,yxBˉi​y means that BiB_iBi​ ranks xxx strictly above yyy. The set of strong orders is ρm\rho_mρm​, and a ballot set is B=(B1,…,Bn)∈ρmnB=(B_1,\dots,B_n)\in\rho^n_mB=(B1​,…,Bn​)∈ρmn​.

A strict voting procedure vnmv^{nm}vnm assigns to each ballot set a single alternative vnm(B)∈Smv^{nm}(B)\in S_mvnm(B)∈Sm​. Its range TpT_pTp​ is the set of alternatives it selects at some ballot set, and p=∣Tp∣p=|T_p|p=∣Tp​∣, so 1≤p≤m1\le p\le m1≤p≤m.

Individual iii can manipulate vnmv^{nm}vnm at BBB if some ballot Bi′B_i'Bi′​ gives

vnm(B1,…,Bi′,…,Bn) Bˉi vnm(B1,…,Bi,…,Bn),v^{nm}(B_1,\dots,B_i',\dots,B_n)\ \bar B_i\ v^{nm}(B_1,\dots,B_i,\dots,B_n),vnm(B1​,…,Bi′​,…,Bn​) Bˉi​ vnm(B1​,…,Bi​,…,Bn​),

an outcome iii strictly prefers by the standard of the original ballot. The procedure is strategy-proof if no individual can manipulate it at any ballot set.

The procedure is dictatorial if some individual iii exists such that, at every ballot set BBB, vnm(B)v^{nm}(B)vnm(B) is the element of the range TpT_pTp​ that BiB_iBi​ ranks highest. It is fully dictatorial if in addition Tp=SmT_p=S_mTp​=Sm​ and partially dictatorial if Tp⊊SmT_p\subsetneq S_mTp​⊊Sm​.

The proof uses three further notions. Condition U requires that if every BiB_iBi​ ranks the same x∈Tpx\in T_px∈Tp​ above all other elements of TpT_pTp​, then vnm(B)=xv^{nm}(B)=xvnm(B)=x. The procedure is weak alternative-excluding if Tp⊊SmT_p\subsetneq S_mTp​⊊Sm​, and strong alternative-excluding if in addition it satisfies Condition U. Finally θT(Bi)\theta_T(B_i)θT​(Bi​) is the restriction of BiB_iBi​ to a subset TTT.

Formalization targets

Goal: Theorem 1 (p. 11)

For a strict committee with n≥1n\ge 1n≥1 and m≥p≥3m\ge p\ge 3m≥p≥3,

vnm is strategy-proof  ⟺  vnm is dictatorial.v^{nm}\ \text{is strategy-proof}\iff v^{nm}\ \text{is dictatorial}.vnm is strategy-proof⟺vnm is dictatorial.

Milestones

  1. Lemma 1 (p. 13): a strategy-proof strict procedure satisfies Condition U.
  2. Lemma 2 (p. 15): with one individual and three alternatives, a strategy-proof procedure is fully dictatorial or strong alternative-excluding.
  3. Lemma 3 (pp. 16–17): for n+1n+1n+1 individuals and three alternatives, individuals 1,…,n1,\dots,n1,…,n cannot manipulate vn+1,3v^{n+1,3}vn+1,3 if and only if each of the procedures B↦vn+1,3(B,Bn+1)B\mapsto v^{n+1,3}(B,B_{n+1})B↦vn+1,3(B,Bn+1​), one for each ballot of individual n+1n+1n+1, is strategy-proof.
  4. Lemma 4 (p. 22): the induction step from nnn to n+1n+1n+1 individuals for three alternatives.
  5. Lemma 5 (p. 25): for n≥1n\ge 1n≥1, m≥3m\ge 3m≥3, a strategy-proof procedure is fully dictatorial or strong alternative-excluding.
  6. Lemma 6 (p. 26): for n≥2n\ge 2n≥2, if θT(Ci)=θT(Di)\theta_{T}(C_i)=\theta_{T}(D_i)θT​(Ci​)=θT​(Di​) for all iii, with T=TpT=T_pT=Tp​, then vnm(C)=vnm(D)v^{nm}(C)=v^{nm}(D)vnm(C)=vnm(D).
  7. Proof of Theorem 1, (27)–(28) (p. 27): in the case of a strategy-proof, strong alternative-excluding vnmv^{nm}vnm with m>p≥3m>p\ge 3m>p≥3, it induces a strategy-proof procedure vnpv^{np}vnp on the alternatives of its range with vnp[θT(B1),…,θT(Bn)]=vnm(B)v^{np}[\theta_T(B_1),\dots,\theta_T(B_n)]=v^{nm}(B)vnp[θT​(B1​),…,θT​(Bn​)]=vnm(B).

Significance

Theorem 1 characterizes the strategy-proof strict voting procedures completely: given p≥3p\ge 3p≥3, they are exactly the procedures that let one individual pick, from a fixed menu TpT_pTp​, their favourite element. The classical statement assumes the procedure can select every alternative; Theorem 1 shows that excluding alternatives from the menu does not help, since a procedure with any menu of at least three alternatives is still governed by one individual. The paper's correspondence theorem with Arrow's social welfare functions (its Theorem 2) and its extension to ballots with indifference (Theorem 1') build on this statement.

The full-range case is on the platform already: AGT.gibbard_satterthwaite (Theorem 9.8 of Nisan–Roughgarden–Tardos–Vazirani, Algorithmic Game Theory) is proved and states that an incentive compatible social choice function onto at least three alternatives has a dictator; it is included here as a reference item and is exactly the case Tp=SmT_p=S_mTp​=Sm​ of Lemma 5. The open item MechanismDesign.SocialChoice.gibbard_satterthwaite (Börgers, Proposition 8.1) poses the full-range equivalence. Theorem 1 contains both as the special case Tp=SmT_p=S_mTp​=Sm​. What this mission adds is the range-restricted statement and the paper's own route to it: the induction on the number of individuals for three alternatives (Lemmas 2–4) and the reduction of a procedure with a proper range to its range (Lemma 6 and (27)–(28)). The paper only sketches Lemma 4 and defers the induction on mmm to Satterthwaite's dissertation, so the formal development fills these gaps.

Difficulty

The "if" direction is immediate. For the "only if" direction the natural first idea is to apply the full-range theorem to the procedure viewed as a map into TpT_pTp​. This fails as stated: the procedure is still defined on rankings of all mmm alternatives, and nothing yet says that its choice ignores how individuals rank alternatives outside TpT_pTp​. That is the content of Lemma 6, which needs strategy-proofness and Condition U; only after it does the procedure factor through rankings of TpT_pTp​, and only then is the factor itself strategy-proof. The full-range case itself is the classical theorem, whose known proofs (via Arrow's theorem, or by pivotal-voter arguments) are combinatorial and long.

Formalization scope

  • Individuals and alternatives are finite types ι and A; nnn = Fintype.card ι, mmm = Fintype.card A, and the n+1n+1n+1 individuals of Lemmas 3–4 are Option ι, with none the individual n+1n+1n+1.
  • Strict ballots and strategy-proofness come from the published definition agt_social (namespace AGT): a ballot is a relation r with IsStrictTotalOrder A r, r x y is the paper's x Bˉi yx\,\bar B_i\,yxBˉi​y, and strategy-proofness is AGT.IncentiveCompatible, which is the paper's definition (1) on ρmn\rho^n_mρmn​.
  • The range is taken over strict ballot sets only, and ppp is its Set.ncard. Theorem 1 assumes n≥1n\ge 1n≥1 and p≥3p\ge 3p≥3; m≥pm\ge pm≥p holds automatically, so m≥3m\ge 3m≥3 is not a separate hypothesis.
  • Dictatorial reads the paper's dictator function fTif^i_TfTi​ as "the BiB_iBi​-maximal element of TpT_pTp​", which for strong orders is unique.
  • θT(Ci)=θT(Di)\theta_T(C_i)=\theta_T(D_i)θT​(Ci​)=θT​(Di​) is encoded as agreement of the two strict relations on T×TT\times TT×T, and the procedure vnpv^{np}vnp of (28) lives on the subtype of the range.
  • Lemmas keep their printed hypotheses (m≥3m\ge 3m≥3, p≥1p\ge 1p≥1, Lemma 6's n≥2n\ge 2n≥2), even where the conclusion would hold without them. The range restriction milestone keeps the proof's case assumptions of strong alternative exclusion and m>p≥3m>p\ge 3m>p≥3.

Two encodings would trivialize the goal and are excluded: a range taken over all relation profiles (junk values would make p≥3p\ge 3p≥3 meaningless), and a dictatorship that only constrains ballot sets at which the dictator's overall favourite lies in TpT_pTp​, or that asks for the dictator's favourite in all of SmS_mSm​ (false whenever Tp⊊SmT_p\subsetneq S_mTp​⊊Sm​). The bound p≥3p\ge 3p≥3 is needed: a procedure with two possible outcomes, such as majority vote between two fixed alternatives, is strategy-proof and not dictatorial, and an imposed procedure (p=1p=1p=1) is dictatorial only vacuously.

Contributions are welcome at every milestone, and the restriction-to-range step (Lemma 6 and (27)–(28)) is reusable for any statement that transfers results from full-range to range-restricted procedures.

Selected references

  • M. A. Satterthwaite, Strategy-proofness and Arrow's Conditions: Existence and Correspondence Theorems for Voting Procedures and Social Welfare Functions, Northwestern University CMS-EMS Discussion Paper No. 122, rev. Dec. 12, 1974; published in J. Econ. Theory 10(2) (1975) 187–217. https://doi.org/10.1016/0022-0531(75)90050-2
  • A. Gibbard, Manipulation of Voting Schemes: A General Result, Econometrica 41(4) (1973) 587–601. https://doi.org/10.2307/1914083
  • N. Nisan, Introduction to Mechanism Design (for Computer Scientists), Chapter 9 of N. Nisan, T. Roughgarden, É. Tardos, V. V. Vazirani (eds.), Algorithmic Game Theory, Cambridge University Press, 2007. https://doi.org/10.1017/CBO9780511800481
  • T. Börgers, An Introduction to the Theory of Mechanism Design, Oxford University Press, 2015, Chapter 8. https://doi.org/10.1093/acprof:oso/9780199734023.001.0001
11 thms4 active usersReviewed
🏆Completed
Machine LearningOperations ResearchTheoretical Computer Science·Captain: mikedeng1

The Weighted Majority Algorithm 7: On a Class Shattered by Its Domain, Every Deterministic Learner's Mistake Bounds Satisfy Σ 2^(−M_i) < 2Research Paper

Motivation

An online learner predicts a binary label for each instance before seeing that instance's correct label. Its guarantee is usually measured by its worst possible number of mistakes against a target function. Such guarantees can differ greatly across target functions: an algorithm might learn one target quickly while taking much longer on another. Littlestone and Warmuth's analysis of the Weighted Majority Algorithm asks which collections of individual mistake guarantees can coexist for one deterministic learner. In Section 7 of their 1994 paper, they isolate a condition under which these guarantees obey a strict numerical constraint. The condition concerns how the target functions vary together over the instance domain, so the result applies without a probability distribution on instances or labels.

The theorem also supports the paper's comparison between arbitrary deterministic learners and Weighted Majority. Given individual mistake bounds, its inequality is the numerical condition that makes a suitable initial weighting of the target functions possible in Theorem 7.1 and Corollary 7.1 of the same paper. The present mission centers on the inequality itself, which is independent of a particular implementation of Weighted Majority.

Setting

Let XXX be an instance domain, possibly finite or infinite, and let φ1,…,φn:X→{0,1}\varphi_1,\ldots,\varphi_n:X\to\{0,1\}φ1​,…,φn​:X→{0,1} be a finite family of target functions. A trial presents an instance x∈Xx\in Xx∈X to a learner, asks it to predict either 000 or 111, and then reveals a label. A finite sequence is consistent with φi\varphi_iφi​ when each revealed label equals φi(x)\varphi_i(x)φi​(x) at that trial's instance. The learner may use the current instance and the complete history of earlier labeled trials when predicting; it does not see the current label first. No sampling or independence assumption is made.

The family is shattered by its domain when every pattern b=(b1,…,bn)∈{0,1}nb=(b_1,\ldots,b_n)\in\{0,1\}^nb=(b1​,…,bn​)∈{0,1}n appears as the vector of function values at some instance:

∀b∈{0,1}n∃x∈X∀i∈{1,…,n}, φi(x)=bi.\forall b\in\{0,1\}^n\quad\exists x\in X\quad\forall i\in\{1,\ldots,n\},\ \varphi_i(x)=b_i.∀b∈{0,1}n∃x∈X∀i∈{1,…,n}, φi​(x)=bi​.

This is a dual notion of shattering. The familiar VC condition asks whether a class of functions realizes every labeling of chosen instances. Here the domain realizes every joint value pattern of a chosen family of functions. The two conditions should not be interchanged.

Let AAA be any deterministic online learner. Write MiM_iMi​ for the supremum of the number of mistakes made by AAA over every finite sequence consistent with φi\varphi_iφi​. Since mistake counts are natural numbers, a finite supremum is an attained natural-number maximum. The theorem assumes each MiM_iMi​ is finite. In the Lean representation, this is a separate finiteness condition on each target's extended-natural mistake bound.

Formalization targets

Strict constraint on individual mistake bounds

The goal is Theorem 7.2 of Littlestone and Warmuth, pp. 244–245. For every domain, every finite family shattered by that domain, and every deterministic online learner with finite individual mistake bounds,

∑i=1n2−Mi<2.\sum_{i=1}^{n}2^{-M_i}<2.i=1∑n​2−Mi​<2.

The inequality is strict. Replacing it by a weak inequality would lose the stated result and the threshold used in the paper's later comparison with Weighted Majority. The theorem quantifies over the learner itself; it does not select a specially behaved learning rule.

Arithmetic milestones

Lemma 7.1 concerns positive real numbers with integral adjacent quotients. When a positive target sss is no larger than their sum and s/r1s/r_1s/r1​ is integral, some prefix sums exactly to sss. Lemma 7.2 concerns positive powers of two with integer exponents. If a power 2ℓ2^\ell2ℓ lies at least as high as each individual weight and at most as high as the total weight, some subset sums exactly to 2ℓ2^\ell2ℓ. Both statements and their exact source wording are recorded as the mission's milestones; they are the numbered arithmetic results used in Theorem 7.2.

Significance

The theorem limits the set of claims that one deterministic learner can make simultaneously. Small MiM_iMi​ contribute large terms to the sum, so a learner cannot have arbitrarily small worst-case mistake counts for all members of a domain-shattered family. This is a worst-case statement about every finite consistent trial sequence, rather than an average-case statement about a chosen data distribution. The paper's Corollary 7.1 uses the resulting inequality to select initial weights for a Weighted Majority comparison, but that later algorithmic statement is outside this mission's goal.

Formalizing this result adds a reusable dual-shattering predicate and an exact constraint on per-target online mistake bounds to the existing online-learning definitions. The paper proves the theorem, while the mission's Lean goal and its two milestone lemmas are open proof obligations. The separate published online-learning module provides the definitions of deterministic predictors, trial histories, mistake counts, and extended-natural worst-case bounds; it does not already prove this theorem.

Difficulty

For a single target, an arbitrary learner may have a small finite mistake bound. The difficulty is relating those bounds across different targets without placing extra restrictions on the learner. Domain shattering gives the adversary control over the simultaneous predictions of all functions at an instance, but the labels of one trial sequence must still remain consistent with a target. A simple count of targets or of possible patterns does not directly produce the strict weighted sum, particularly when the different MiM_iMi​ are unequal. The two arithmetic lemmas address exact sums of structured weights; the main theorem must also connect such weights to the learner's worst-case mistakes.

Formalization scope

The function family is indexed by Fin n, with zero-based Lean indices corresponding to the paper's 1,…,n1,\ldots,n1,…,n. The domain XXX is an arbitrary type and may be infinite. Labels and predictions are Bool. ShatteredByDomain states that every Boolean pattern is realized by some domain point. The published UnderstandingML_Online definition gives OnlineAlg X Bool, finite trial histories, mistakes, and mistakeBound. The latter is valued in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞} and ranges over all finite sequences labeled by the target. Each MiM_iMi​ is converted to a natural number only after assuming its bound is finite. The powers 2−Mi2^{-M_i}2−Mi​ and their finite sum are real numbers. The n=0n=0n=0 conclusion is the empty sum 0<20<20<2; the shattering predicate itself then requires an instance to realize the unique empty pattern.

The two milestones use real-valued weights. Lemma 7.1 assumes n≥1n\ge1n≥1 because the source refers to r1r_1r1​. Lemma 7.2 represents the source's integral base-two logarithms by equations ri=2kir_i=2^{k_i}ri​=2ki​ with integer kik_iki​, which also enforce positivity. Its maximum hypothesis is written pointwise, so an empty family is harmless: the other inequality 2ℓ≤02^\ell\le02ℓ≤0 cannot hold. A coordinate family on X={0,1}nX=\{0,1\}^nX={0,1}n is verified locally to satisfy the shattering predicate, ruling out a vacuous encoding of that central assumption. Contributions proving either arithmetic milestone or the main theorem within these exact definitions are in scope.

Selected references

  • Nick Littlestone and Manfred K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 212–261, 1994. DOI: 10.1006/inco.1994.1009.
  • Shai Shalev-Shwartz and Shai Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 21. DOI: 10.1017/CBO9781107298019.
5 thms4 active usersReviewed
🏆Completed
Graph TheoryLinear OptimizationOperations Research·Captain: mikedeng1

Matching, Euler Tours and the Chinese Postman I: The Convex Hull of Nonnegative Integer Parity Solutions Is the Polyhedron of Odd-Set InequalitiesResearch Paper

Motivation

The Chinese postman problem, posed by Meigu Guan (Kwan Mei-Ko) in 1962, asks for a shortest closed walk in a graph that traverses every edge at least once: the route of a postman who must cover every street and return to the post office. It is the basic model of arc routing (street sweeping, snow ploughing, meter reading, inspection of networks), and it is one of the first combinatorial optimization problems shown to be solvable in polynomial time by matching methods.

Edmonds and Johnson, Matching, Euler tours and the Chinese postman (Mathematical Programming 5, 1973), reduce the postman problem to a parity problem: choose nonnegative integers xex_exe​, the number of extra traversals of each edge, so that every node becomes even, at minimum total length. Their §3 gives a complete linear description of this parity problem. The description is the first instance of what is now called the TTT-join polyhedron (with TTT the set of odd-degree nodes), a standard object of polyhedral combinatorics and a frequent example of total dual integrality.

Timeline:

  • 1962: Guan states the postman problem and an optimality criterion for it.
  • 1965: Edmonds proves the matching polyhedron theorem and gives the blossom algorithm for weighted matching (Edmonds 1965).
  • 1973: Edmonds and Johnson specialise both to the postman problem, giving the polyhedron in the edge variables alone (this mission) and a blossom algorithm that solves it (Edmonds–Johnson 1973).

Setting

A graph GGG has a finite set NNN of nodes and a finite set EEE of edges; each edge meets two different nodes, and several edges may meet the same pair of nodes. The degree of a node nnn is the number of edges meeting it, ∑e∈Eane\sum_{e\in E} a_{ne}∑e∈E​ane​, where ane=1a_{ne}=1ane​=1 if eee meets nnn and 000 otherwise. Set bn=1b_n=1bn​=1 when the degree of nnn is odd (an odd node) and bn=0b_n=0bn​=0 otherwise.

An edge eee meets a set S⊆NS\subseteq NS⊆N when exactly one of its two ends lies in SSS. A set S⊆NS\subseteq NS⊆N is an odd set when it contains an odd number of odd nodes, and any number of even nodes.

The parity points X⊆REX\subseteq\mathbb R^EX⊆RE are the vectors with

xe∈Z,xe≥0(e∈E),∑e∈Eanexe≡bn(mod2)(n∈N),x_e\in\mathbb Z,\quad x_e\ge 0\quad(e\in E),\qquad \sum_{e\in E}a_{ne}x_e\equiv b_n \pmod 2\quad(n\in N),xe​∈Z,xe​≥0(e∈E),e∈E∑​ane​xe​≡bn​(mod2)(n∈N),

the conditions (3.1), (3.2), (3.6) of the paper. Adding xex_exe​ copies of each edge eee to GGG makes every degree even, and an Euler tour of the enlarged graph is a postman tour of GGG. The postman polyhedron is

P={x∈RE: xe≥0 (e∈E),  ∑{xe:e meets S}≥1 for every odd set S},P=\Big\{x\in\mathbb R^E:\ x_e\ge 0\ (e\in E),\ \ \sum\{x_e : e\text{ meets }S\}\ge 1\ \text{for every odd set } S\Big\},P={x∈RE: xe​≥0 (e∈E),  ∑{xe​:e meets S}≥1 for every odd set S},

the conditions (3.2) and (3.5). For lengths c∈REc\in\mathbb R^Ec∈RE the objective is z=∑ecexez=\sum_e c_e x_ez=∑e​ce​xe​. The dual linear program has a variable ySy_SyS​ for each odd set, constraints yS≥0y_S\ge 0yS​≥0 and ∑{yS:e meets S}≤ce\sum\{y_S: e\text{ meets }S\}\le c_e∑{yS​:e meets S}≤ce​, and objective v=∑SySv=\sum_S y_Sv=∑S​yS​.

Formalization targets

Goal: the polyhedron theorem (p. 93)

conv⁡X=P.\operatorname{conv} X = P .convX=P.

Both inclusions are claimed. The equality is between subsets of RE\mathbb R^ERE, so it also asserts that the convex hull of XXX is closed.

Milestones

In the paper's order of argument:

  1. (3.13), p. 95: an integer solution of (3.6) has ∑{xe:e meets S}≡1(mod2)\sum\{x_e : e\text{ meets }S\}\equiv 1 \pmod 2∑{xe​:e meets S}≡1(mod2) for every odd set SSS.
  2. p. 95: conv⁡X⊆P\operatorname{conv}X\subseteq PconvX⊆P.
  3. (3.10), p. 94: weak duality, v≤zv\le zv≤z for x∈Px\in Px∈P and yyy dual feasible.
  4. (3.11)–(3.12), p. 95: complementary slackness makes xxx and yyy optimal.
  5. p. 95: for c≥0c\ge 0c≥0 there is an integral parity point xxx and a dual feasible yyy satisfying (3.11) and (3.12). The paper proves this with the algorithm of §4.
  6. p. 96: if some ce<0c_e<0ce​<0, zzz is unbounded below on PPP and on XXX.
  7. (3.19), p. 96: for c≥0c\ge 0c≥0, the minimum of zzz over PPP is attained at a parity point.

There are two companion statements. The first (pp. 93–94) says that in ∑eanexe−2wn=bn\sum_e a_{ne}x_e-2w_n=b_n∑e​ane​xe​−2wn​=bn​ with x≥0x\ge 0x≥0, integrality of wnw_nwn​ already forces wn≥0w_n\ge 0wn​≥0. The second (p. 91) is the polyhedron theorem in the variables (x,w)(x,w)(x,w): the convex hull of the integer solutions of (3.1)–(3.3) is the set of solutions of (3.2), (3.2′), (3.3), (3.5).

Significance

The theorem turns the parity problem, an integer program, into a linear program over an explicitly described polyhedron. Optimal postman tours therefore come with dual certificates: a packing of odd cuts whose total value matches the tour's extra length. The same description underlies the theory of TTT-joins and TTT-cuts (Schrijver 2003, Chapter 29).

No machine-checked proof of the TTT-join polyhedron theorem exists, on Prove2Me or, to our knowledge, in Mathlib. This mission states it in the multigraph setting of the paper. Milestones 1–4 and 6 are elementary and can serve as warm-ups. Milestone 5, total dual integrality of the odd-cut system, is the substantive step. The paper obtains it from a blossom algorithm; the milestone is stated as an existence result, so any proof of it closes it.

Difficulty

The inclusion conv⁡X⊆P\operatorname{conv}X\subseteq PconvX⊆P is a parity count. The reverse inclusion is the theorem. Two first ideas fail. Total unimodularity does not apply, because the odd-cut constraint matrix is not totally unimodular and PPP has many constraints whose tightness interacts. Reading the result off the matching polyhedron theorem needs the auxiliary loops wnw_nwn​ and an elimination step, and the 1965 matching polyhedron describes perfect matchings of a simple graph, not this parity system on a multigraph. Some argument for integrality of the minimum for every c≥0c\ge 0c≥0 is needed, and the paper's is an algorithm with a dual certificate. A further point is that PPP and XXX are unbounded: equality of hulls is not just agreement of minima over bounded objectives, and the case of objectives with a negative coefficient has to be handled separately (milestone 6).

Formalization scope

A graph is the local ChinesePostman.Polyhedron.Graph V E: a map ends : E → Sym2 V with a proof that no edge is a loop, over finite types V and E. Mathlib's SimpleGraph is not used, because the paper allows parallel edges. Degrees, odd nodes and bnb_nbn​ are computed from the graph; the set of odd nodes is never a free parameter. "eee meets SSS" means exactly one end in SSS. Odd sets range over all subsets S⊆NS\subseteq NS⊆N with an odd number of odd nodes, with no nonemptiness or properness condition. Parity points are real images of integer vectors E → ℤ, so integrality is part of XXX. Dual variables are functions Finset V → ℝ whose values matter only on odd sets: the sign condition and every sum are restricted to the odd sets. Lengths ccc are arbitrary reals except where the page assumes ce≥0c_e\ge 0ce​≥0 (milestones 5 and 7).

Trivializing readings are ruled out: (3.5) taken over nonempty proper sets only, "meets" read as "has at least one end in", the odd nodes taken as an arbitrary set TTT, or XXX allowed to contain non-integral points would each give a different, and in some cases false or trivial, statement.

Every definition the statements need is in one file, ChinesePostman.Polyhedron.Setting. A complete development needs finite-dimensional convexity (Mathlib's convexHull) and linear programming duality or an integrality argument. A Lean proof of total dual integrality for odd-cut systems would be reusable for matching and TTT-join problems beyond this paper. Proofs of any milestone, and alternative proofs of milestone 5, are welcome.

Selected references

  • J. Edmonds and E. L. Johnson, Matching, Euler tours and the Chinese postman, Mathematical Programming 5 (1973) 88–124. https://doi.org/10.1007/BF01580113
  • J. Edmonds, Maximum matching and a polyhedron with 0,1-vertices, Journal of Research of the National Bureau of Standards 69B (1965) 125–130. https://doi.org/10.6028/jres.069B.013
  • M. Guan (Kwan Mei-Ko), Graphic programming using odd or even points, Chinese Mathematics 1 (1962) 273–277.
  • A. Schrijver, Combinatorial Optimization: Polyhedra and Efficiency, Springer, 2003, Chapter 29 (T-joins and T-cuts). https://doi.org/10.1007/978-3-540-44389-6
11 thms4 active usersReviewed
PreviousPage 1 of 14Next

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