Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Combinatorics

630 missions · 352 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

Open278Completed352All630
Graph TheoryLinear algebra·Captain: mikedeng1

Quick Approximation to Matrices and Applications III: The Atoms of a Cut Decomposition with Error εn² Form a 2ε-Pseudo-Regular Partition of a GraphResearch Paper

Motivation

Szemerédi's regularity lemma partitions the vertex set of any graph into a bounded number of parts between which edges are distributed almost randomly. It is a central tool of extremal graph theory, but the number of parts it needs is a tower of exponentials in 1/ϵ1/\epsilon1/ϵ, which makes it unusable in algorithms. In 1996 and 1999 Frieze and Kannan introduced a weaker notion, now called weak regularity: a partition is good if, for every pair of disjoint vertex sets S,TS,TS,T, the number of edges between them is predicted, up to an additive error ϵn2\epsilon n^2ϵn2, by the densities between the parts. This weaker requirement is enough to approximate dense instances of Max-Cut and other problems, and it can be met with only 2O(1/ϵ2)2^{O(1/\epsilon^2)}2O(1/ϵ2) parts.

The paper Frieze and Kannan, Quick approximation to matrices and applications, Combinatorica 19 (1999) derives such partitions from a matrix decomposition. Section 5.1 shows that if the adjacency matrix of a graph is approximated by a sum of a few cut matrices with small error in the cut norm, then the atoms of the sets defining those cut matrices form a pseudo-regular partition. The weak regularity of this paper later became the basis of the cut-distance theory of graph limits (Lovász, Large networks and graph limits, AMS, 2012).

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a simple graph with n=∣V∣n=|V|n=∣V∣ vertices and adjacency matrix A\mathbf AA. For a real V×VV\times VV×V matrix M\mathbf MM and S,T⊆VS,T\subseteq VS,T⊆V write M(S,T)=∑i∈S∑j∈TM(i,j)\mathbf M(S,T)=\sum_{i\in S}\sum_{j\in T}\mathbf M(i,j)M(S,T)=∑i∈S​∑j∈T​M(i,j). A cut matrix CUT(R,C,d)\mathrm{CUT}(R,C,d)CUT(R,C,d) has entry ddd on R×CR\times CR×C and 000 elsewhere. A cut decomposition of width sss is a sum D=D(1)+⋯+D(s)\mathbf D=\mathbf D^{(1)}+\dots+\mathbf D^{(s)}D=D(1)+⋯+D(s) of cut matrices D(t)=CUT(Rt,Ct,dt)\mathbf D^{(t)}=\mathrm{CUT}(R_t,C_t,d_t)D(t)=CUT(Rt​,Ct​,dt​); its error is W=A−D\mathbf W=\mathbf A-\mathbf DW=A−D.

Let P=V1,…,Vk\mathcal P=V_1,\dots,V_kP=V1​,…,Vk​ be a partition of VVV with index set KKK, write Xi=X∩ViX_i=X\cap V_iXi​=X∩Vi​, let e(S,T)e(S,T)e(S,T) be the number of edges between disjoint sets SSS and TTT, and let di,j=e(Vi,Vj)/(∣Vi∣∣Vj∣)d_{i,j}=e(V_i,V_j)/(|V_i||V_j|)di,j​=e(Vi​,Vj​)/(∣Vi​∣∣Vj​∣) be the density between parts. The deviation from regularity is

ΔP(S,T)=e(S,T)−∑i∈K∑j∈Kdi,j∣Si∣∣Tj∣,\Delta_{\mathcal P}(S,T)=e(S,T)-\sum_{i\in K}\sum_{j\in K}d_{i,j}|S_i||T_j|,ΔP​(S,T)=e(S,T)−i∈K∑​j∈K∑​di,j​∣Si​∣∣Tj​∣,

and P\mathcal PP is ϵ\epsilonϵ-pseudo-regular if ∣ΔP(S,T)∣≤ϵn2|\Delta_{\mathcal P}(S,T)|\le\epsilon n^2∣ΔP​(S,T)∣≤ϵn2 for all disjoint S,T⊆VS,T\subseteq VS,T⊆V.

For a partition Q=W1,…,Wq\mathcal Q=W_1,\dots,W_qQ=W1​,…,Wq​, the block average MQ\mathbf M_{\mathcal Q}MQ​ of M\mathbf MM takes the value M(Wi,Wj)/(∣Wi∣∣Wj∣)\mathbf M(W_i,W_j)/(|W_i||W_j|)M(Wi​,Wj​)/(∣Wi​∣∣Wj​∣) on Wi×WjW_i\times W_jWi​×Wj​, and M\mathbf MM is compatible with Q\mathcal QQ if it is constant on every block Wi×WjW_i\times W_jWi​×Wj​. The atom partition of R1,…,Rs,C1,…,CsR_1,\dots,R_s,C_1,\dots,C_sR1​,…,Rs​,C1​,…,Cs​ is the coarsest partition of VVV in which every RtR_tRt​ and CtC_tCt​ is a union of parts; it has at most 4s4^s4s parts.

Formalization targets

Goal: claim (50), p. 204

If the cut decomposition has error at most ϵn2\epsilon n^2ϵn2, that is

∣A(S,T)−D(S,T)∣≤ϵn2for all S,T⊆V,|\mathbf A(S,T)-\mathbf D(S,T)|\le\epsilon n^2\quad\text{for all }S,T\subseteq V,∣A(S,T)−D(S,T)∣≤ϵn2for all S,T⊆V,

then the atom partition V1,…,VkV_1,\dots,V_kV1​,…,Vk​ of the sets Rt,CtR_t,C_tRt​,Ct​ satisfies

∣ΔP(S,T)∣≤2ϵn2for all disjoint S,T⊆V,|\Delta_{\mathcal P}(S,T)|\le 2\epsilon n^2\quad\text{for all disjoint }S,T\subseteq V,∣ΔP​(S,T)∣≤2ϵn2for all disjoint S,T⊆V,

i.e. it is 2ϵ2\epsilon2ϵ-pseudo-regular. The width sss is unconstrained: the statement holds for decompositions of any width, and the width only enters through the bound k≤4sk\le 4^sk≤4s on the number of parts.

Milestones

  1. (47), p. 203. For disjoint S,TS,TS,T: A(S,T)=e(S,T)\mathbf A(S,T)=e(S,T)A(S,T)=e(S,T), AQ(S,T)=∑i,jdi,j∣Si∣∣Tj∣\mathbf A_{\mathcal Q}(S,T)=\sum_{i,j}d_{i,j}|S_i||T_j|AQ​(S,T)=∑i,j​di,j​∣Si​∣∣Tj​∣, hence
A(S,T)−AQ(S,T)=ΔQ(S,T).\mathbf A(S,T)-\mathbf A_{\mathcal Q}(S,T)=\Delta_{\mathcal Q}(S,T).A(S,T)−AQ​(S,T)=ΔQ​(S,T).
  1. D is compatible with P\mathcal PP, p. 204. The sum of the cut matrices is constant on every block of the atom partition.
  2. Lemma 7(a), p. 203. If Q\mathcal QQ refines P\mathcal PP and M\mathbf MM is compatible with P\mathcal PP, then
sup⁡S,T⊆V∣A(S,T)−AQ(S,T)∣≤2sup⁡S,T⊆V∣A(S,T)−M(S,T)∣.\sup_{S,T\subseteq V}|\mathbf A(S,T)-\mathbf A_{\mathcal Q}(S,T)|\le 2\sup_{S,T\subseteq V}|\mathbf A(S,T)-\mathbf M(S,T)|.S,T⊆Vsup​∣A(S,T)−AQ​(S,T)∣≤2S,T⊆Vsup​∣A(S,T)−M(S,T)∣.

Companions

  • Lemma 7(b), p. 204. If M\mathbf MM is compatible with P\mathcal PP, then ∥A−AP∥F≤∥A−M∥F\|\mathbf A-\mathbf A_{\mathcal P}\|_F\le\|\mathbf A-\mathbf M\|_F∥A−AP​∥F​≤∥A−M∥F​.
  • k≤4sk\le 4^sk≤4s, p. 204. The atom partition of 2s2s2s sets has at most 4s4^s4s parts.

Significance

Claim (50) is the bridge from matrix approximation to graph partitions. Combined with Theorem 7 of the same paper (every matrix has a cut decomposition of width at most 1/ϵ21/\epsilon^21/ϵ2 with cut-norm error ϵmn∥A∥F≤ϵn2\epsilon\sqrt{mn}\|\mathbf A\|_F\le\epsilon n^2ϵmn​∥A∥F​≤ϵn2 for a graph), it gives the weak regularity lemma: every graph has a 2ϵ2\epsilon2ϵ-pseudo-regular partition into at most 41/ϵ24^{1/\epsilon^2}41/ϵ2 parts. Pseudo-regularity in turn means that e(S,T)e(S,T)e(S,T) is determined up to 2ϵn22\epsilon n^22ϵn2 by the numbers ∣Si∣,∣Tj∣|S_i|,|T_j|∣Si​∣,∣Tj​∣, which reduces Max-Cut and related problems on dense graphs to a search over a bounded-dimensional space. Lemma 7(a) is of independent use: Section 5.1.1 of the paper applies it again to pass to an equitable refinement.

The results are proved in the paper. Mathlib contains Szemerédi's regularity lemma in its equipartition form, and the platform has a published statement of it; neither the Frieze–Kannan weak regularity statement, the deviation ΔP\Delta_{\mathcal P}ΔP​, nor Lemma 7 has a machine-checked proof known to this mission. The mission produces a formal statement and proof of the partition half of the weak regularity lemma, reusable with any cut decomposition, including the existence result of the companion mission on cut decompositions.

Difficulty

The identity (47) and the compatibility of D\mathbf DD are bookkeeping, but bookkeeping over Mathlib's Finpartition, atomise and the rational-valued edgeDensity, with sums over parts that must be regrouped into sums over vertices. The substance lies in Lemma 7(a). The obvious argument bounds ∣A−AQ∣|\mathbf A-\mathbf A_{\mathcal Q}|∣A−AQ​∣ by ∣A−M∣+∣M−AQ∣|\mathbf A-\mathbf M|+|\mathbf M-\mathbf A_{\mathcal Q}|∣A−M∣+∣M−AQ​∣, and the difficulty is the second term: M−AQ\mathbf M-\mathbf A_{\mathcal Q}M−AQ​ is constant on the blocks of Q\mathcal QQ, but its block sums on arbitrary S,TS,TS,T have to be controlled by block sums of A−M\mathbf A-\mathbf MA−M, on which A\mathbf AA and AQ\mathbf A_{\mathcal Q}AQ​ need not agree. The printed argument handles only disjoint S,TS,TS,T, while the lemma is stated, and needed, for all pairs; a complete proof must cover pairs that overlap.

Formalization scope

All objects are defined in one definitions item in the namespace FriezeKannan.PseudoReg. The vertex set is a Fintype VVV with decidable equality, the graph a SimpleGraph V with decidable adjacency, and its adjacency matrix Mathlib's adjMatrix over R\mathbb RR. A partition of VVV is a Finpartition of Finset.univ, whose parts are the index set KKK; refinement is Mathlib's order Q≤P\mathcal Q\le\mathcal PQ≤P. Densities are SimpleGraph.edgeDensity, cast to R\mathbb RR, for all pairs of parts. For distinct parts this is the paper's d(Vi,Vj)d(V_i,V_j)d(Vi​,Vj​). For a part with itself it is the block average A(Vi,Vi)/∣Vi∣2\mathbf A(V_i,V_i)/|V_i|^2A(Vi​,Vi​)/∣Vi​∣2, not the value e(A,A)/(∣A∣2)e(A,A)/\binom{|A|}{2}e(A,A)/(2∣A∣​) printed in the §5 preamble: with the printed value the identity (47) fails on diagonal blocks, and the proof of (50) uses (47). The atom partition is Mathlib's Finpartition.atomise of the family {Rt}∪{Ct}\{R_t\}\cup\{C_t\}{Rt​}∪{Ct​}. The Frobenius norm is the square root of the sum of squared entries.

The paper takes the cut matrices from its Theorem 2, a randomized algorithm. The goal replaces this by the one property its proof uses: the hypothesis ∣A(S,T)−D(S,T)∣≤ϵn2|\mathbf A(S,T)-\mathbf D(S,T)|\le\epsilon n^2∣A(S,T)−D(S,T)∣≤ϵn2 for all S,T⊆VS,T\subseteq VS,T⊆V (Theorem 2 bounds the cut norm, a maximum over all pairs, by ϵn∥A∥F≤ϵn2\epsilon n\|\mathbf A\|_F\le\epsilon n^2ϵn∥A∥F​≤ϵn2). The paper's symmetrised error Wˉ\bar{\mathbf W}Wˉ is not needed. Lemma 7 is stated for an arbitrary real matrix A\mathbf AA, and the suprema of Lemma 7(a) are written as "every bound ccc on the right-hand family gives the bound 2c2c2c on the left-hand family", over all pairs S,TS,TS,T. No hypothesis is added to any statement.

The partition in the goal is fixed as the atom partition. A formalization that let the partition vary would be trivial: the partition into singletons has Δ≡0\Delta\equiv0Δ≡0.

Needed infrastructure: regrouping sums over a Finpartition, and block sums of matrices constant on the blocks of a partition, on arbitrary and on block-compatible sets. These are reusable for any statement about cut norms of block-constant matrices. Proofs of the milestones, of the companions, and of the goal are all welcome.

Selected references

  • A. Frieze and R. Kannan, Quick approximation to matrices and applications, Combinatorica 19(2):175–220, 1999. https://doi.org/10.1007/s004930050052
  • A. Frieze and R. Kannan, The regularity lemma and approximation schemes for dense problems, Proc. 37th FOCS, 12–20, 1996. https://doi.org/10.1109/SFCS.1996.548459
  • E. Szemerédi, Regular partitions of graphs, Problèmes combinatoires et théorie des graphes, Colloq. Internat. CNRS 260, 399–401, 1978.
  • L. Lovász, Large networks and graph limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
5 thms1 active userReviewed
Graph TheoryOperations Research·Captain: mikedeng1

Edge-Disjoint Spanning Trees of Finite Graphs: A Finite Multigraph Has k Edge-Disjoint Spanning Trees iff Every Partition P of Its Vertices Is Crossed by at Least k(|P| − 1) EdgesResearch Paper

Why edge-disjoint spanning trees

A connected network survives the failure of any single link exactly when it has no bridge, but a stronger and more useful property is to carry several edge-disjoint spanning trees: each tree can broadcast to every node on its own, so kkk such trees give kkk independent routing or broadcast structures, and the network stays connected after any k−1k-1k−1 link failures. The question of when a graph contains kkk edge-disjoint spanning trees was answered in 1961, simultaneously and independently, by W. T. Tutte (On the problem of decomposing a graph into n connected factors) and C. St.J. A. Nash-Williams (Edge-disjoint spanning trees of finite graphs). The answer is a min-max condition over vertex partitions, now called the Tutte–Nash-Williams tree-packing theorem. It is a standard result in combinatorial optimization and the graphic case of Edmonds's matroid base-packing theorem (1965).

Timeline:

  • 1961. Tutte proves the theorem in an equivalent form; Nash-Williams, unaware of Tutte's work, gives a different proof a few months later, the one formalized here.
  • 1964. Nash-Williams proves the companion covering theorem: the edges can be covered by kkk forests iff every vertex set UUU spans at most k(∣U∣−1)k(|U|-1)k(∣U∣−1) edges.
  • 1965. Edmonds derives both results from matroid partition (Minimum partition of a matroid into independent subsets).

Setting

A graph GGG here is a finite unoriented multigraph in which every edge joins two distinct vertices; two vertices may be joined by several edges. Write VVV for its vertex set and EEE for its edge set. A tree on a non-empty vertex set WWW is an edge set TTT, all of whose edges have both ends in WWW, containing no cycle (two parallel edges count as a cycle), such that any two vertices of WWW are joined by a path of TTT-edges. A spanning tree of GGG is a tree on all of VVV using edges of GGG. Spanning trees are edge-disjoint if no two share an edge.

A partition PPP of VVV is a set of non-empty, pairwise disjoint subsets of VVV whose union is VVV; ∣P∣|P|∣P∣ is the number of its members. The crossing edges EP(G)E_P(G)EP​(G) are the edges whose two ends lie in different members of PPP, counted with multiplicity. Throughout, kkk is a fixed positive integer. For X⊆VX \subseteq VX⊆V, EXE_XEX​ is the set of edges with both ends in XXX, eX=∣EX∣e_X = |E_X|eX​=∣EX​∣, and

ΔG(X)=k(∣X∣−1)−eX.\Delta_G(X) = k(|X|-1) - e_X .ΔG​(X)=k(∣X∣−1)−eX​.

Formalization targets

Goal: Theorem 1

G has k edge-disjoint spanning trees  ⟺  ∣EP(G)∣≥k(∣P∣−1) for every partition P of V.(1)G \text{ has } k \text{ edge-disjoint spanning trees} \iff |E_P(G)| \ge k(|P|-1) \text{ for every partition } P \text{ of } V. \tag{1}G has k edge-disjoint spanning trees⟺∣EP​(G)∣≥k(∣P∣−1) for every partition P of V.(1)

Milestones, in the order of the proof

  1. Lemma 1: a tree on WWW has ∣W∣−1|W|-1∣W∣−1 edges. Lemma 2: every connected graph has a spanning tree.
  2. Necessity: kkk edge-disjoint spanning trees imply (1) for every PPP.
  3. (*) On couples [G,g][G,g][G,g] (g≥0g \ge 0g≥0 on vertices, ΔG≥0\Delta_G \ge 0ΔG​≥0 on non-empty sets): critical sets (ΔG=0\Delta_G = 0ΔG​=0) whose intersection is non-empty are closed under ∩\cap∩ and ∪\cup∪. The same holds for crucial sets (Γ=0\Gamma = 0Γ=0, where Γ(X)=ΔG(X)−s+g . Xˉ\Gamma(X) = \Delta_G(X) - s + g\,.\,\bar XΓ(X)=ΔG​(X)−s+g.Xˉ) when the couple is sss-good.
  4. Lemma 3: for s≥1s \ge 1s≥1, an sss-good couple has an (s−1)(s-1)(s−1)-good supercouple with one added edge. Corollary 3A: an sss-good couple has an sss-supercouple, obtained by adding sss edges.
  5. (†) and Lemma 4: a spanning tree of a graph produced by fusions at ξ\xiξ (replacing edges ξη\xi\etaξη, ξζ\xi\zetaξζ by a new edge ηζ\eta\zetaηζ) pulls back to a spanning tree of the original graph.
  6. Lemma 5: a graph with ∣E∣=k(∣V∣−1)|E| = k(|V|-1)∣E∣=k(∣V∣−1) satisfies (1) iff ΔG(X)≥0\Delta_G(X) \ge 0ΔG​(X)≥0 for every non-empty XXX.
  7. Lemma 6: under the inductive hypothesis of the proof, a critical partition other than {V}\{V\}{V} and the partition into singletons yields kkk edge-disjoint spanning trees.

Significance

The theorem gives an exact, checkable certificate in both directions: a family of kkk trees certifies the packing, and a single partition violating (1) certifies that no packing exists. Corollaries include that every 2k2k2k-edge-connected graph has kkk edge-disjoint spanning trees, and that the maximum number of edge-disjoint spanning trees (the strength-type parameter used in network reliability and in Nagamochi–Ibaraki style connectivity algorithms) is the minimum of ⌊∣EP(G)∣/(∣P∣−1)⌋\lfloor |E_P(G)|/(|P|-1) \rfloor⌊∣EP​(G)∣/(∣P∣−1)⌋ over partitions with ∣P∣≥2|P| \ge 2∣P∣≥2. The partition condition is also the prototype of the rank condition in matroid base packing.

The result itself is classical and fully proved. What this mission adds is a machine-checked version in a multigraph setting. Mathlib has trees and spanning trees for simple graphs (SimpleGraph.IsTree, SimpleGraph.Connected.exists_isTree_le), but not edge-disjoint tree packing, and its simple graphs cannot express the parallel edges that the theorem and its proof need. Related items on the platform include the multigraph layer NagamochiIbaraki.EdgeConn.Multigraph (reused here), the Keller–Trotter spanning-forest count for simple graphs, and Edmonds's matroid partition theorem with Nash-Williams's covering corollary. Covering and packing are different statements, so none of these yields Theorem 1 directly. Formalizing Tutte's alternative proof, or deriving Theorem 1 from a formalized matroid union theorem, would also be welcome.

Difficulty

Necessity is a counting argument. The difficulty lies in sufficiency. The obvious induction deletes an edge or contracts a set of vertices, but deleting an edge can destroy (1), and contracting changes the vertex set, so neither reduction preserves the hypothesis on its own. The hardest case is a graph in which only the partition into singletons is critical. There a vertex of degree less than 2k2k2k has to be removed while (1) is preserved on the remaining graph, and the trees found there have to be converted back into trees of the original graph. Both steps need careful bookkeeping of edges incident with that vertex. A simple-graph development cannot carry the argument, because the edges that the reduction adds among the neighbours of the removed vertex may be parallel to existing ones.

Formalization scope

  • Graphs. A vertex type V, an edge type E (both finite, in Type), and ends : E → Sym2 V, the published encoding of NagamochiIbaraki.EdgeConn. Loops are excluded by the hypothesis ∀ e, ¬ (ends e).IsDiag. Parallel edges are allowed, and simplicity is never assumed. A subgraph is a pair (W : Finset V, F : Finset E).
  • Trees. "Tree" means non-empty vertex set, edges inside it, no cycle (IsForest: every edge is a bridge) and connected. ∣T∣=∣W∣−1|T| = |W|-1∣T∣=∣W∣−1 is Lemma 1, not part of the definition. "kkk edge-disjoint spanning trees" is a family Fin k → Finset E of pairwise disjoint spanning trees.
  • Partitions are Mathlib Finpartition (Finset.univ : Finset V). (1) is required for every partition, including the one-part partition and the singletons. All counts involving subtraction (ΔG\Delta_GΔG​, (1), criticality, Γ\GammaΓ) are in Z\mathbb ZZ.
  • Explicit readings. VVV is non-empty in Theorem 1, since the paper's graphs have a vertex and on V=∅V = \emptysetV=∅ condition (1) holds while no tree exists. k≥1k \ge 1k≥1 in every statement. The paper's "X⊂V(G)X \subset V(G)X⊂V(G)" includes X=V(G)X = V(G)X=V(G) and is read as ⊆\subseteq⊆. "Exactly g(ξ)−h(ξ)g(\xi)-h(\xi)g(ξ)−h(ξ) new edges" is written deg⁡new(ξ)+h(ξ)=g(ξ)\deg_{\text{new}}(\xi) + h(\xi) = g(\xi)degnew​(ξ)+h(ξ)=g(ξ), without natural-number subtraction. A supergraph with sss added edges has edge type E ⊕ Fin s, so the added edges always exist; none of them may be a loop. A fusion uses a given fresh edge ρ : E in an ambient edge type, and a sequence of fusions is a list, each valid on the graph produced by the later ones (L=Φ1⋯ΦnGL = \Phi_1\cdots\Phi_n GL=Φ1​⋯Φn​G applies Φn\Phi_nΦn​ first). In (†) and Lemma 4, "is a tree" means "is a spanning tree of GGG", since the vertex set is all of VVV. Lemma 6 carries the paper's inductive hypothesis as an explicit hypothesis: Theorem 1 (admissible implies kkk trees) for all loopless multigraphs with smaller ∣V∣+∣E∣|V|+|E|∣V∣+∣E∣.
  • Ruled out. Defining a spanning tree by an edge count alone, or by acyclicity alone, would make Lemma 1 or the goal trivial; the definitions require both acyclicity and connectivity on the stated vertex set. Taking added edges from a fixed ambient type would make Lemma 3 false. Dropping or globalizing Lemma 6's inductive hypothesis would turn it into the goal or a tautology.
  • Infrastructure. Multigraph trees, the tree edge count and spanning-tree existence, contraction GPG_PGP​ of a partition, and splitting-off at a vertex. All of these are reusable beyond this mission. Contributions of general multigraph lemmas are welcome.

Selected references

  • C. St.J. A. Nash-Williams, Edge-disjoint spanning trees of finite graphs, J. London Math. Soc. 36 (1961), 445–450. https://doi.org/10.1112/jlms/s1-36.1.445
  • W. T. Tutte, On the problem of decomposing a graph into n connected factors, J. London Math. Soc. 36 (1961), 221–230. https://doi.org/10.1112/jlms/s1-36.1.221
  • C. St.J. A. Nash-Williams, Decomposition of finite graphs into forests, J. London Math. Soc. 39 (1964), 12. https://doi.org/10.1112/jlms/s1-39.1.12
  • J. Edmonds, Minimum partition of a matroid into independent subsets, J. Res. Nat. Bur. Standards 69B (1965), 67–72. https://doi.org/10.6028/jres.069B.004
  • C. Berge, Théorie des graphes et ses applications, Dunod, Paris, 1958.
15 thms1 active userReviewed
Graph Theory·Captain: mikedeng1

Graph Minors. X. Obstructions to Tree-Decomposition III: In the θ-Grid, the Separations of Order < θ With Small E(A) Form a Tangle of Order θResearch Paper

Why a grid has a tangle

Robertson and Seymour introduced tangles as a way to describe a region of a graph that cannot be isolated by a separation of low order. A separation divides the edges into two sides, and its order counts the vertices shared by those sides. A tangle makes a consistent choice of the side regarded as small at every order below a threshold. This turns an informal idea of a highly connected region into an object with explicit axioms. In Graph Minors. X, the authors connect tangles with branch-width, tree-width and grid minors; the grid construction is the direct example behind that connection. The paper states the existence of a tangle of order θ\thetaθ in every θ\thetaθ-grid as (7.3) on p. 173. Robertson and Seymour, 1991.

The result is already proved in the paper. The present target is its Lean statement and proof, together with the two numbered results from §7 on which the construction rests. This mission concerns the grid itself. The paper's further assertion that a grid minor induces a tangle in a larger graph uses the separate minor-transfer result (6.1), whose edge identities require additional definitions. Robertson and Seymour, 1991, pp. 170–173.

The grid, its edge boundary and its small sets

Fix an integer θ≥2\theta\ge2θ≥2. The θ\thetaθ-grid GGG has vertices (i,j)(i,j)(i,j) for 1≤i,j≤θ1\le i,j\le\theta1≤i,j≤θ. Two vertices are adjacent exactly when their coordinate differences have absolute values summing to one. For each iii, PiP_iPi​ is the horizontal row with vertices (i,1),…,(i,θ)(i,1),\ldots,(i,\theta)(i,1),…,(i,θ); for each jjj, QjQ_jQj​ is the vertical column. Write E(Pi)E(P_i)E(Pi​) and E(Qj)E(Q_j)E(Qj​) for their edge sets. The Lean development uses the already published definition of this same simple graph, with coordinates shifted to start at zero. The shift changes no adjacency or row and column relation. Robertson and Seymour, 1991, p. 171.

For an edge set X⊆E(G)X\subseteq E(G)X⊆E(G), the edge boundary ∂(X)\partial(X)∂(X) consists of vertices incident both with an edge in XXX and with an edge outside XXX. The source prints “vertices v∈Xv\in Xv∈X” when introducing this set. Since XXX is a set of edges, this is a type slip; the two incidence conditions give the intended reading. An edge set is small if ∣∂(X)∣<θ|\partial(X)|<\theta∣∂(X)∣<θ and it contains no complete row edge set E(Pi)E(P_i)E(Pi​). These are two separate requirements. An edge set can have a small boundary while still containing a whole row. Robertson and Seymour, 1991, pp. 171–172.

A subhypergraph AAA specifies subsets V(A)V(A)V(A) and E(A)E(A)E(A) of the grid's vertices and edges, retaining the incidences of the ambient grid. A separation (A,B)(A,B)(A,B) satisfies A∪B=GA\cup B=GA∪B=G and E(A∩B)=∅E(A\cap B)=\varnothingE(A∩B)=∅; its order is ∣V(A∩B)∣|V(A\cap B)|∣V(A∩B)∣. Although the grid is a simple graph, these are the paper's hypergraph definitions. A tangle T\mathcal TT of order θ\thetaθ contains only separations of order below θ\thetaθ and satisfies three axioms: it chooses an orientation of every such separation; the union of the first sides of any three chosen separations is not GGG; and no chosen first side contains every vertex of GGG. The chosen three may coincide. Robertson and Seymour, 1991, pp. 153–154.

Formalization targets

The first milestone, (7.1), identifies a symmetry between rows and columns under a low boundary:

∣∂(X)∣<θ⟹(∃i, E(Pi)⊆X)⟺(∃j, E(Qj)⊆X).|\partial(X)|<\theta\quad\Longrightarrow\quad \bigl(\exists i,\ E(P_i)\subseteq X\bigr) \Longleftrightarrow \bigl(\exists j,\ E(Q_j)\subseteq X\bigr).∣∂(X)∣<θ⟹(∃i, E(Pi​)⊆X)⟺(∃j, E(Qj​)⊆X).

The second milestone, (7.2), says that three small edge sets cannot cover all the grid's edges:

X1∪X2∪X3=E(G)⟹¬(X1,X2,X3 are all small).X_1\cup X_2\cup X_3=E(G) \quad\Longrightarrow\quad \neg\bigl(X_1,X_2,X_3\text{ are all small}\bigr).X1​∪X2​∪X3​=E(G)⟹¬(X1​,X2​,X3​ are all small).

The goal, (7.3), defines Tθ\mathcal T_\thetaTθ​ rather than taking it as an arbitrary set. It selects precisely the separations (A,B)(A,B)(A,B) of order less than θ\thetaθ for which E(A)E(A)E(A) is small, and asserts

Tθ is a tangle in G of order θ.\mathcal T_\theta\text{ is a tangle in }G\text{ of order }\theta.Tθ​ is a tangle in G of order θ.

These are the numbered §7 results on pp. 171–173, in the order used by the paper. Robertson and Seymour, 1991, pp. 171–173.

What the result supplies

The theorem gives an explicit tangle in a familiar graph for every order θ≥2\theta\ge2θ≥2. A tangle is more than a numerical connectivity bound: it records coherent choices across all separations below its order. In the paper, the grid tangle feeds the link between large grids and tangles of high order. The extension from a grid minor in another graph to a tangle of that graph invokes (6.1), beyond this mission's goal. The converse direction, that large tree-width forces a large grid minor, is (7.4), proved in Graph Minors V (Robertson and Seymour, 1986); the published definitions and statements of that paper on this platform are about grid minors and tree-width, not tangles. Robertson and Seymour, 1991, pp. 170, 173.

Formalizing (7.3) will make the grid tangle available as a checked object alongside the already published grid graph. It will also produce reusable definitions for an edge boundary, small edge sets and the paper's hypergraph separations. The Lean files here are draft statements with open proofs; the published 1986 grid definition is the existing checked ingredient. No machine-checked proof of (7.1), (7.2) or (7.3) is claimed by this proposal.

Where the difficulty lies

The obvious attempt at (7.2) counts boundary vertices column by column. If every column QjQ_jQj​ has all its edges in X1∪X2X_1\cup X_2X1​∪X2​ and meets both sets, each column contributes at least two vertices to ∣∂(X1)∣+∣∂(X2)∣|\partial(X_1)|+|\partial(X_2)|∣∂(X1​)∣+∣∂(X2​)∣, so the two boundaries together have at least 2θ2\theta2θ vertices. This only works when the third set misses the first and last rows; in general the third set can absorb edges of every column, and no single count over the whole grid bounds all three boundaries at once. The statement is a covering property of three sets, not a bound on any one of them, and the paper's argument for it is credited to D. Kleitman and M. Saks. The goal must then meet all three tangle axioms. The page's proof of (7.3) argues the first axiom and says only that the rest follows from (7.2); the third axiom, V(A)≠V(G)V(A)\ne V(G)V(A)=V(G), is part of the definition and must be established, not dropped. Robertson and Seymour, 1991, pp. 172–173.

Formalization scope

The source assumes finite hypergraphs throughout, and §7 fixes θ≥2\theta\ge2θ≥2. Lean uses a finite simple graph on Fin θ × Fin θ, reusing RobertsonSeymour1986.GM5.grid, and views its edge set as the edge type of a local hypergraph. Subhypergraphs hold vertex and edge subsets with inherited incidence. Set.ncard counts boundary vertices and separation order. Rows and columns use all grid edges whose two endpoints have the specified coordinate. This representation keeps every grid edge's identity, even though it is written as a SimpleGraph edge subtype.

IsSmall retains both ∣∂(X)∣<θ|\partial(X)|<\theta∣∂(X)∣<θ and the exclusion of full row edge sets. boundary requires an incident edge on each side of XXX. Both details are necessary for the claimed tangle: an IsSmall without the row clause would call every set with a small boundary small, including sets that contain whole rows, and the resulting set of separations is not a tangle; a boundary defined with "every edge outside XXX" instead of "some edge outside XXX" would be empty and make the size condition vacuous. The formalization keeps the third tangle axiom and does not replace it with an easier condition. The definitions of hypergraph, separation and tangle duplicate the shared encoding used by the other missions for this paper and can later be consolidated. The definitions also include the paper's tangle number and maximum edge size for that shared layer, though the grid goal does not invoke them; on finite inputs their maxima are bounded, and their empty-set value is zero.

The companion missions are Graph Minors. X. Obstructions to Tree-Decomposition I: Tangle Number Equals max(Branch-Width, Largest Edge Size), max(β(G), γ(G)) = θ(G); Graph Minors. X. Obstructions to Tree-Decomposition II: Branch-Width and Tree-Width Satisfy max(β, γ) ≤ ω + 1 ≤ max(⌊3β/2⌋, γ, 1); Graph Minors. X. Obstructions to Tree-Decomposition IV: Mutually Distinguishable Tangles Have a Tree-Decomposition Separating Them by Their Distinctions; and Graph Minors. X. Obstructions to Tree-Decomposition V: A θ-Pervasive Class of Designs Gives a Tree-Decomposition over 𝒮^{3θ−2} ∪ ℛ_{4θ−3}. Those missions address the paper's other main results; this one isolates the grid construction.

Selected references

  • N. Robertson and P. D. Seymour, Graph Minors. X. Obstructions to Tree-Decomposition, Journal of Combinatorial Theory, Series B 52 (1991), 153–190. DOI: 10.1016/0095-8956(91)90061-n.
  • N. Robertson and P. D. Seymour, Graph Minors. V. Excluding a Planar Graph, Journal of Combinatorial Theory, Series B 41 (1986), 92–114. DOI: 10.1016/0095-8956(86)90030-4.
9 thms1 active userReviewed
Graph Theory·Captain: mikedeng1

Graph Minors. X. Obstructions to Tree-Decomposition V: A θ-Pervasive Class of Designs Gives a Tree-Decomposition over 𝒮^{3θ−2} ∪ ℛ_{4θ−3}Research Paper

Motivation

Tangles were introduced by Robertson and Seymour in Graph Minors. X (JCTB 1991) to make precise the idea of a "highly connected region" of a graph. A tangle orients every low-order separation of the graph towards one side, consistently. They are the organising notion of the later papers of the Graph Minors series, which culminate in the graph structure theorem and the proof of Wagner's conjecture.

Section 11 of Graph Minors X, Structure relative to a tangle, asks how local and global structure are related. In the paper's words: suppose that GGG has large tangle number, but relative to each high-order tangle the graph has a structure of a certain kind XXX. What can be inferred about the global structure of GGG? One might guess that GGG has a tree-decomposition into pieces each with structure XXX, but that is false. The theorem of §11 shows that GGG does have a tree-decomposition into pieces which almost have structure XXX. The authors state that they need this for an application in a later paper of the series, cited as [6], Graph minors. XVII. Excluding a non-planar graph (submitted); it appeared as Graph Minors XVI (2003), where the local structure relative to a tangle is a near-embedding in a surface.

Setting

A hypergraph GGG has a finite vertex set V(G)V(G)V(G), a finite edge set E(G)E(G)E(G) and an incidence relation. A subhypergraph G′⊆GG' \subseteq GG′⊆G has V(G′)⊆V(G)V(G') \subseteq V(G)V(G′)⊆V(G) and E(G′)⊆E(G)E(G') \subseteq E(G)E(G′)⊆E(G), with every end of an edge of G′G'G′ in V(G′)V(G')V(G′). A separation of GGG is a pair (A,B)(A, B)(A,B) of subhypergraphs with A∪B=GA \cup B = GA∪B=G and E(A∩B)=∅E(A \cap B) = \emptysetE(A∩B)=∅, and its order is ∣V(A∩B)∣|V(A \cap B)|∣V(A∩B)∣.

For θ≥1\theta \ge 1θ≥1, a tangle of order θ\thetaθ in GGG is a set T\mathcal TT of separations of order <θ< \theta<θ such that (i) for each separation (A,B)(A, B)(A,B) of order <θ< \theta<θ one of (A,B)(A,B)(A,B), (B,A)(B,A)(B,A) is in T\mathcal TT; (ii) no three members (Ai,Bi)(A_i, B_i)(Ai​,Bi​) of T\mathcal TT have A1∪A2∪A3=GA_1 \cup A_2 \cup A_3 = GA1​∪A2​∪A3​=G; (iii) V(A)≠V(G)V(A) \ne V(G)V(A)=V(G) for (A,B)∈T(A, B) \in \mathcal T(A,B)∈T.

A tree-decomposition (T,τ)(T, \tau)(T,τ) of GGG is a tree TTT and subhypergraphs τ(t)\tau(t)τ(t), t∈V(T)t \in V(T)t∈V(T), whose union is GGG, which are pairwise edge-disjoint, and with τ(t)∩τ(t′′)⊆τ(t′)\tau(t) \cap \tau(t'') \subseteq \tau(t')τ(t)∩τ(t′′)⊆τ(t′) whenever t′t't′ lies on the path of TTT between ttt and t′′t''t′′.

A design is a pair (H,M)(H, M)(H,M) of a hypergraph HHH and a set MMM of subsets of V(H)V(H)V(H).

  • The design of a node t0t_0t0​ of (T,τ)(T, \tau)(T,τ) with neighbours t1,…,tkt_1, \dots, t_kt1​,…,tk​ is (τ(t0),{V(τ(t0)∩τ(ti))})(\tau(t_0), \{V(\tau(t_0) \cap \tau(t_i))\})(τ(t0​),{V(τ(t0​)∩τ(ti​))}). A tree-decomposition is over a class S\mathcal SS of designs if the design of every node belongs to S\mathcal SS.
  • (H′,M′)(H', M')(H′,M′) is an nnn-enlargement of (H,M)(H, M)(H,M) if, for some Z⊆V(H′)Z \subseteq V(H')Z⊆V(H′) with ∣Z∣≤n|Z| \le n∣Z∣≤n, H⊆H′H \subseteq H'H⊆H′, V(H′)−V(H)⊆ZV(H') - V(H) \subseteq ZV(H′)−V(H)⊆Z, E(H′)⊆E(H)E(H') \subseteq E(H)E(H′)⊆E(H), and X∩V(H)∈MX \cap V(H) \in MX∩V(H)∈M for every X∈M′X \in M'X∈M′ other than ZZZ. Sn\mathcal S^nSn is the class of nnn-enlargements of members of S\mathcal SS, and Rn\mathcal R_nRn​ is the class of designs with at most nnn vertices.
  • A location in G′G'G′ is a set of separations (Ai,Bi)(A_i, B_i)(Ai​,Bi​) of G′G'G′ with Ai⊆BjA_i \subseteq B_jAi​⊆Bj​ for i≠ji \ne ji=j. Its design is (G′∩B1∩⋯∩Bk,{V(Ai∩Bi)})(G' \cap B_1 \cap \dots \cap B_k, \{V(A_i \cap B_i)\})(G′∩B1​∩⋯∩Bk​,{V(Ai​∩Bi​)}).
  • A class S\mathcal SS is θ\thetaθ-pervasive in GGG if for every subhypergraph G′⊆GG' \subseteq GG′⊆G and every tangle T\mathcal TT in G′G'G′ of order ≥θ\ge \theta≥θ, some location L⊆T\mathcal L \subseteq \mathcal TL⊆T in G′G'G′ has its design in S\mathcal SS.
  • The ZZZ-extension of (H,M)(H, M)(H,M) is (H,M∪{Z})(H, M \cup \{Z\})(H,M∪{Z}).

Formalization targets

Goal: (11.1), p. 185

For θ≥1\theta \ge 1θ≥1 and a class S\mathcal SS that is θ\thetaθ-pervasive in GGG,

G has a tree-decomposition over S3θ−2∪R4θ−3.G \text{ has a tree-decomposition over } \mathcal S^{3\theta-2} \cup \mathcal R_{4\theta-3}.G has a tree-decomposition over S3θ−2∪R4θ−3​.

Milestone: (11.2), pp. 185–186

If S\mathcal SS is θ\thetaθ-pervasive in GGG and ∣Z∣=3θ−2|Z| = 3\theta - 2∣Z∣=3θ−2, then either some separation (A,B)(A, B)(A,B) of order <θ< \theta<θ has

∣(Z∪V(A))∩V(B)∣, ∣(Z∪V(B))∩V(A)∣≤3θ−3,|(Z \cup V(A)) \cap V(B)|,\ |(Z \cup V(B)) \cap V(A)| \le 3\theta - 3,∣(Z∪V(A))∩V(B)∣, ∣(Z∪V(B))∩V(A)∣≤3θ−3,

or some location {(Ai,Bi)}\{(A_i, B_i)\}{(Ai​,Bi​)} in GGG with design in S\mathcal SS has ∣Z∩V(Ai)∣≤∣V(Ai∩Bi)∣<θ|Z \cap V(A_i)| \le |V(A_i \cap B_i)| < \theta∣Z∩V(Ai​)∣≤∣V(Ai​∩Bi​)∣<θ for all iii.

Milestone: (11.3), p. 186

If ∣Z∣≤3θ−2|Z| \le 3\theta - 2∣Z∣≤3θ−2, there is a tree-decomposition of GGG over S′=S3θ−2∪R4θ−3\mathcal S' = \mathcal S^{3\theta-2} \cup \mathcal R_{4\theta-3}S′=S3θ−2∪R4θ−3​ with a node t0t_0t0​ such that Z⊆V(τ(t0))Z \subseteq V(\tau(t_0))Z⊆V(τ(t0​)) and the ZZZ-extension of the design of t0t_0t0​ is in S′\mathcal S'S′. The goal is the case Z=∅Z = \emptysetZ=∅.

Companion: the remark on p. 188

For every θ≥1\theta \ge 1θ≥1, if GGG has no tangle of order θ\thetaθ, then

ω(G)≤4θ−4.\omega(G) \le 4\theta - 4.ω(G)≤4θ−4.

Consequently ω(G)≤4θ(G)\omega(G) \le 4\theta(G)ω(G)≤4θ(G), where θ(G)\theta(G)θ(G) is the tangle number. This is the p. 188 remark following (11.3).

Significance

(11.1) converts local information, stated tangle by tangle, into a global decomposition. Every piece of the decomposition is either small (at most 4θ−34\theta - 34θ−3 vertices) or differs from a structure in S\mathcal SS by at most 3θ−23\theta - 23θ−2 added vertices. This is the kind of local-to-global step the later Graph Minors papers need: there the local structure near each high-order tangle is a near-embedding in a surface, and a result of this shape glues such local structures along a tree. With S=∅\mathcal S = \emptysetS=∅ it gives a tree-width bound ω(G)≤4θ(G)\omega(G) \le 4\theta(G)ω(G)≤4θ(G), a weaker form of the tangle/tree-width duality.

The result has been proved since 1991. To our knowledge no proof assistant formalizes it, and none formalizes the hypergraph tangles or tree-decompositions of this paper. This mission produces machine-checked statements of the paper's notions of design, enlargement, location and pervasiveness, the structure theorem, and its inductive strengthening.

Difficulty

The first idea is to take the pieces in S\mathcal SS themselves as the pieces of the tree-decomposition. That fails, and the paper says so: the pieces can only almost have the structure. Locations for different tangles interact, and the decomposition must absorb up to 3θ−23\theta - 23θ−2 extra vertices into each piece. The induction also has to control a set ZZZ of attachment vertices of bounded size: a piece built at one node must still accept the vertices through which it is glued to the rest of the tree. Keeping ∣Z∣≤3θ−2|Z| \le 3\theta - 2∣Z∣≤3θ−2 throughout the recursion is the delicate point, and it is why the strengthened statement (11.3), not (11.1), is the one that admits induction.

Formalization scope

  • Hypergraphs. A hypergraph is a structure Hypergraph V E with an incidence relation. Finiteness is [Finite V] [Finite E] on every theorem. Subhypergraphs, separations, tangles and tree-decompositions are the series' shared encodings (namespace RobertsonSeymour1991.GM10.Structure). The tree of a tree-decomposition is a SimpleGraph on Fin n satisfying IsTree.
  • Designs. Every design the statements involve has as hypergraph a subhypergraph of the fixed GGG. So a design is a GGG-design G.Sub × Set (Set V), and a class of designs is a set of GGG-designs. Classes need not be closed under isomorphism; the paper never uses that.
  • Relative tangles. Separations, tangles and locations in a subhypergraph G′G'G′ are pairs of subhypergraphs of GGG whose union is G′G'G′. For G′=GG' = GG′=G they agree definitionally with the absolute notions.
  • Pervasiveness. It quantifies over every subhypergraph G′G'G′ and every order θ′≥θ\theta' \ge \thetaθ′≥θ. Restricting to G′=GG' = GG′=G would be a different, stronger theorem. The paper's opening remark in the proof of (11.3), that pervasiveness passes to subhypergraphs, is immediate in this encoding and is not a separate item.
  • Edge cases. The empty set is a location; its design is (G′,∅)(G', \emptyset)(G′,∅).
  • Integer arithmetic. Since θ≥1\theta \ge 1θ≥1, the bounds 3θ−23\theta-23θ−2, 3θ−33\theta-33θ−3 and 4θ−34\theta-34θ−3 are exact natural numbers. Tree-width is an integer, and it is −1-1−1 for an empty vertex set.
  • No trivial instances. IsPervasive is satisfiable: the class of all designs is θ\thetaθ-pervasive, through the empty location. The conclusion is satisfiable too: the one-node decomposition lies over R∣V(G)∣\mathcal R_{|V(G)|}R∣V(G)∣​. So neither the hypothesis nor the conclusion of (11.1) is vacuous.

The other main results of the paper are separate missions of this series:

  • Graph Minors. X. Obstructions to Tree-Decomposition I (tangle number versus branch-width);
  • II (branch-width versus tree-width);
  • III (the grid tangle);
  • IV (the tree-decomposition separating mutually distinguishable tangles).

Contributions are welcome on any item: proofs of (11.2) and (11.3), and of the p. 188 companion. The companion needs two auxiliary facts: a tangle of order ≥θ\ge \theta≥θ in a subhypergraph yields one of order θ\thetaθ in GGG, and tree-width is bounded by bag size.

Selected references

  • N. Robertson, P. D. Seymour, Graph Minors. X. Obstructions to Tree-Decomposition, J. Combin. Theory Ser. B 52 (1991) 153–190. https://doi.org/10.1016/0095-8956(91)90061-n
  • N. Robertson, P. D. Seymour, Graph Minors. XVI. Excluding a non-planar graph, J. Combin. Theory Ser. B 89 (2003) 43–76. https://doi.org/10.1016/S0095-8956(03)00042-X
  • N. Robertson, P. D. Seymour, Graph Minors. V. Excluding a planar graph, J. Combin. Theory Ser. B 41 (1986) 92–114. https://doi.org/10.1016/0095-8956(86)90030-4
7 thms1 active userReviewed
Graph Theory·Captain: mikedeng1

Graph Minors. X. Obstructions to Tree-Decomposition I: Tangle Number Equals max(Branch-Width, Largest Edge Size), max(β(G), γ(G)) = θ(G)Research Paper

Motivation

The Graph Minors series of Robertson and Seymour proves that finite graphs are well-quasi-ordered under the minor relation. A central tool of the series is the duality between two ways of measuring how "tree-like" a graph is: decompositions along a tree, which certify that a graph is thin, and tangles, which certify that it contains a highly connected region that no small separation can split. Graph Minors X (Robertson and Seymour, 1991) introduces tangles and proves the exact minimax relation between the tangle number and branch-width, an invariant within a constant factor of tree-width. Tangles became a standard object of structural graph theory and were later abstracted to connectivity systems, matroids and data clustering; branch-width is the parameter behind dynamic-programming algorithms on graphs and matroids of bounded width.

Setting

A hypergraph GGG consists of a finite vertex set V(G)V(G)V(G), a finite edge set E(G)E(G)E(G) and an incidence relation; an edge may be incident with any number of vertices, its ends, and the size of an edge is its number of ends. Write γ(G)\gamma(G)γ(G) for the maximum size of an edge (γ(G)=0\gamma(G)=0γ(G)=0 if E(G)=∅E(G)=\emptysetE(G)=∅).

A subhypergraph is given by a set of vertices and a set of edges containing all ends of its edges. A separation of GGG is a pair (A,B)(A,B)(A,B) of subhypergraphs with A∪B=GA\cup B=GA∪B=G and no common edge; its order is ∣V(A∩B)∣|V(A\cap B)|∣V(A∩B)∣.

For an integer θ≥1\theta\ge1θ≥1, a tangle of order θ\thetaθ is a set T\mathcal TT of separations of order <θ<\theta<θ such that (i) for every separation (A,B)(A,B)(A,B) of order <θ<\theta<θ, one of (A,B),(B,A)(A,B),(B,A)(A,B),(B,A) lies in T\mathcal TT; (ii) A1∪A2∪A3≠GA_1\cup A_2\cup A_3\ne GA1​∪A2​∪A3​=G for any three members (Ai,Bi)(A_i,B_i)(Ai​,Bi​) of T\mathcal TT; (iii) V(A)≠V(G)V(A)\ne V(G)V(A)=V(G) for every (A,B)∈T(A,B)\in\mathcal T(A,B)∈T. The tangle number θ(G)\theta(G)θ(G) is the largest order of a tangle in GGG, or 000 if there is none.

A ternary tree is a tree whose vertices have valency 111 or 333. A branch-decomposition (T,τ)(T,\tau)(T,τ) of GGG is a ternary tree TTT with a bijection τ\tauτ from its leaves onto E(G)E(G)E(G). Deleting an edge fff of TTT splits the leaves, hence E(G)E(G)E(G), into two parts; the order of fff is the number of vertices of GGG incident with edges on both sides. The width is the largest order of an edge of TTT, and the branch-width β(G)\beta(G)β(G) is the minimum width over all branch-decompositions, with β(G)=0\beta(G)=0β(G)=0 when ∣E(G)∣≤1|E(G)|\le1∣E(G)∣≤1.

The proof goes through an abstract setting. On a finite set EEE, a connectivity function is a symmetric submodular map κ:2E→Z\kappa:2^E\to\mathbb Zκ:2E→Z. A set XXX is efficient if κ(X)≤0\kappa(X)\le0κ(X)≤0. A bias is a set of efficient sets containing one of X,E−XX,E-XX,E−X for each efficient XXX and no three members covering EEE. Its dual object is a tree-labelling: a ternary tree whose incidences carry efficient sets, complementary across each edge, covering EEE at each internal vertex, and absorbed by a given family A\mathcal AA at each leaf.

Formalization targets

Goal: (4.3), p. 165

max⁡(β(G),γ(G))=θ(G)unless γ(G)=0 and V(G)≠∅.\max\bigl(\beta(G),\gamma(G)\bigr)=\theta(G)\qquad\text{unless }\gamma(G)=0\text{ and }V(G)\ne\emptyset.max(β(G),γ(G))=θ(G)unless γ(G)=0 and V(G)=∅.

The exception is genuine: a hypergraph with a vertex and only edges without ends has β=γ=0\beta=\gamma=0β=γ=0 and θ=1\theta=1θ=1.

Milestones

  1. Tangle lemmas of §2: (2.2) closure under (A∪A′,B∩B′)(A\cup A',B\cap B')(A∪A′,B∩B′); (2.3) E(A1∪A2∪A3)≠E(G)E(A_1\cup A_2\cup A_3)\ne E(G)E(A1​∪A2​∪A3​)=E(G) for order ≥2\ge2≥2; (2.4) an edge of size ≥θ\ge\theta≥θ defines a tangle of order θ\thetaθ; (2.5) GGG has a tangle iff V(G)≠∅V(G)\ne\emptysetV(G)=∅; (2.7) the third axiom reduces to separations (Ke,G∖e)(K_e,G\setminus e)(Ke​,G∖e); (2.8) extreme separations have order θ−1\theta-1θ−1, and their large side is well connected.
  2. The abstract minimax of §3: (3.1) a bias extending A\mathcal AA excludes tree-labellings over A\mathcal AA; (3.2) tree-labellings can be made exact on the same tree; (3.3) leaf sets of an exact tree-labelling partition EEE; (3.4) no bias gives an exact tree-labelling; (3.5) the three-way equivalence; (3.6) leaf labels can be taken different from EEE.
  3. The three claims of the proof of (4.3), for γ(G)>0\gamma(G)>0γ(G)>0, k≥γ(G)k\ge\gamma(G)k≥γ(G), κ=κ0−k\kappa=\kappa_0-kκ=κ0​−k and A={{e}}\mathcal A=\{\{e\}\}A={{e}}: biases correspond to tangles of order k+1k+1k+1, exact tree-labellings to β(G)≤k\beta(G)\le kβ(G)≤k, and hence GGG has a tangle of order k+1k+1k+1 iff k<β(G)k<\beta(G)k<β(G).

A companion statement, (4.4), computes θ(Kn)=⌈2n/3⌉\theta(K_n)=\lceil 2n/3\rceilθ(Kn​)=⌈2n/3⌉ and, for n≥3n\ge3n≥3, β(Kn)=⌈2n/3⌉\beta(K_n)=\lceil 2n/3\rceilβ(Kn​)=⌈2n/3⌉.

Significance

The minimax theorem makes tangles the exact obstruction to small branch-width: either a branch-decomposition of width at most kkk exists, or a tangle of order k+1k+1k+1 witnesses that none does. Every later use of tangles as certificates of high connectivity in the Graph Minors series, including the tree-decomposition of a graph into regions indexed by its maximal tangles and the grid theorem, rests on this duality. The abstract lemma (3.5) applies to any connectivity function, so it also gives the duality for matroid branch-width and for other symmetric submodular functions.

The results are proved in the published paper; none has a machine-checked proof that this mission is aware of. The platform has branch-width objects for set functions on subcubic trees, from the Oum–Seymour approximation of clique-width (ApproxCliqueWidth.Certificate.bw_ge_of_wellLinked, bw_le_of_no_wellLinked), which bound branch-width by well-linked sets, a different obstruction. This mission adds a hypergraph tangle library and the exact minimax.

Difficulty

The inequality θ(G)≥γ(G)\theta(G)\ge\gamma(G)θ(G)≥γ(G) and the direction "a tangle of order k+1k+1k+1 excludes width ≤k\le k≤k" are short. The converse, building a branch-decomposition of width ≤k\le k≤k when no tangle of order k+1k+1k+1 exists, is the hard direction. A natural first attempt is to merge small separations greedily into a tree, but the submodular uncrossing step can destroy the leaf conditions already established, and an induction on the edge set does not preserve the tangle hypothesis. A second difficulty is bookkeeping: the trees are ternary, so every construction that cuts, joins or contracts trees must restore valencies 111 and 333, and the passage between leaf labels and edges of GGG needs the labels at the leaves to partition E(G)E(G)E(G).

Formalization scope

The objects live in the namespace RobertsonSeymour1991.GM10.Minimax.

  • Hypergraphs are a structure with an incidence predicate on arbitrary types V, E; every theorem assumes [Finite V] [Finite E], as the paper assumes all hypergraphs finite.
  • Subhypergraphs are determined by their vertex and edge sets.
  • A tangle is a predicate IsTangle G θ 𝒯 with θ≥1\theta\ge1θ≥1 built in, and the three members in the second axiom need not be distinct. The tangle number is a supremum in N\mathbb NN of a set bounded by ∣V(G)∣|V(G)|∣V(G)∣, so it is a maximum.
  • Ternary trees and the trees of branch-decompositions and tree-labellings live on Fin n. A branch-decomposition stores τ\tauτ as an injection of the edges onto the leaves.
  • Branch-width is 000 for ∣E(G)∣≤1|E(G)|\le1∣E(G)∣≤1 and otherwise an infimum over a nonempty set.
  • In §3 subsets are Set E, and E−XE-XE−X is the complement. The incidence (u,e)(u,e)(u,e) with e=uwe=uwe=uw is the adjacent pair (u,w)(u,w)(u,w), and every condition on labels is required on adjacent pairs only.
  • The claims of the proof of (4.3) carry the paper's own context, γ(G)>0\gamma(G)>0γ(G)>0 and k≥γ(G)k\ge\gamma(G)k≥γ(G), as hypotheses.

Junk values would make the goal trivial: a tangle number of 000 from an unbounded set of orders, or a branch-width of 000 from an empty set of decompositions. Neither occurs: every tangle has order at most ∣V(G)∣|V(G)|∣V(G)∣, and a branch-decomposition exists whenever ∣E(G)∣≥2|E(G)|\ge2∣E(G)∣≥2. The goal also keeps the exact exception ¬(γ(G)=0∧V(G)≠∅)\neg(\gamma(G)=0\wedge V(G)\neq\emptyset)¬(γ(G)=0∧V(G)=∅), neither dropped nor strengthened.

A complete development needs finite trees with leaf counting, path and component reasoning in SimpleGraph, and submodular uncrossing arguments. The tangle library of §2 and the abstract minimax of §3 can be reused beyond this mission. Contributions of tree lemmas (ternary trees on Fin n, grafting two trees at a leaf, contracting a valency-2 vertex) are welcome.

The other main results of the paper are separate missions of this series: Graph Minors. X. … II: Branch-Width and Tree-Width Satisfy max(β, γ) ≤ ω + 1 ≤ max(⌊3β/2⌋, γ, 1), III: In the θ-Grid, the Separations of Order < θ With Small E(A) Form a Tangle of Order θ, IV: Mutually Distinguishable Tangles Have a Tree-Decomposition Separating Them by Their Distinctions, and V: A θ-Pervasive Class of Designs Gives a Tree-Decomposition over 𝒮^{3θ−2} ∪ ℛ_{4θ−3}.

Selected references

  • N. Robertson, P. D. Seymour, Graph Minors. X. Obstructions to Tree-Decomposition, J. Combin. Theory Ser. B 52 (1991) 153–190. https://doi.org/10.1016/0095-8956(91)90061-n
  • N. Robertson, P. D. Seymour, Graph Minors. V. Excluding a Planar Graph, J. Combin. Theory Ser. B 41 (1986) 92–114. https://doi.org/10.1016/0095-8956(86)90030-4
  • S. Oum, P. D. Seymour, Approximating clique-width and branch-width, J. Combin. Theory Ser. B 96 (2006) 514–528. https://doi.org/10.1016/j.jctb.2005.10.006
25 thms1 active userReviewed
Graph Theory·Captain: mikedeng1

Graph Minors. X. Obstructions to Tree-Decomposition II: Branch-Width and Tree-Width Satisfy max(β, γ) ≤ ω + 1 ≤ max(⌊3β/2⌋, γ, 1)Research Paper

Motivation

Tree-width measures how closely a graph resembles a tree. It is the central width parameter of the Graph Minors series of N. Robertson and P. D. Seymour, which proves that finite graphs are well-quasi-ordered by the minor relation, and it governs the complexity of many algorithms: problems that are NP-hard in general become solvable in linear time on graphs of bounded tree-width (Courcelle's theorem). Branch-width, introduced in Graph Minors X (Robertson, Seymour 1991), is a second width parameter, defined from ternary trees whose leaves carry the edges of the graph. It is the natural parameter for the paper's main minimax theorem, which equates branch-width with the maximum order of a tangle, the obstruction to small decompositions.

Section 5 of Graph Minors X shows that the two parameters are equivalent, within a factor 3/2. This is what transfers the tangle duality from branch-width to tree-width, and it is why later work (for instance on matroid branch-width and on rank-width) treats branch-width as an interchangeable substitute for tree-width.

Setting

All objects are finite. A hypergraph GGG has a vertex set V(G)V(G)V(G), an edge set E(G)E(G)E(G), and an incidence relation; the ends of an edge are the vertices incident with it, and an edge may have any number of ends. γ(G)\gamma(G)γ(G) is the maximum number of ends of an edge, or 000 if E(G)=∅E(G)=\emptysetE(G)=∅. A subhypergraph consists of subsets of V(G)V(G)V(G) and E(G)E(G)E(G) such that every end of a chosen edge is a chosen vertex; unions, intersections and inclusion of subhypergraphs are taken setwise.

A tree is a connected non-null graph without circuits; its leaves are its vertices of valency at most 111, and it is ternary if every valency is 111 or 333.

A branch-decomposition of GGG is a pair (T,τ)(T,\tau)(T,τ) with TTT a ternary tree and τ\tauτ a bijection from the leaves of TTT onto E(G)E(G)E(G). The order of an edge fff of TTT is the number of vertices vvv of GGG such that vvv is an end both of an edge labelling a leaf in one component of T∖fT\setminus fT∖f and of an edge labelling a leaf in the other. The width of (T,τ)(T,\tau)(T,τ) is the maximum order of an edge of TTT, and the branch-width β(G)\beta(G)β(G) is the minimum width, or 000 if ∣E(G)∣≤1|E(G)|\le 1∣E(G)∣≤1.

A tree-decomposition of GGG is a pair (T,τ)(T,\tau)(T,τ) with TTT a tree and τ(t)\tau(t)τ(t) a subhypergraph of GGG for every t∈V(T)t\in V(T)t∈V(T), such that

  1. the union of all τ(t)\tau(t)τ(t) is GGG;
  2. for distinct t,t′t,t't,t′, the parts τ(t)\tau(t)τ(t) and τ(t′)\tau(t')τ(t′) share no edge;
  3. if t′t't′ lies on the path of TTT between ttt and t′′t''t′′, then τ(t)∩τ(t′′)⊆τ(t′)\tau(t)\cap\tau(t'')\subseteq\tau(t')τ(t)∩τ(t′′)⊆τ(t′).

Its width is max⁡t (∣V(τ(t))∣−1)\max_t\,(|V(\tau(t))|-1)maxt​(∣V(τ(t))∣−1), and the tree-width ω(G)\omega(G)ω(G) is the minimum width; ω(G)=−1\omega(G)=-1ω(G)=−1 exactly when V(G)=∅V(G)=\emptysetV(G)=∅. Each part is a subhypergraph and each edge lies in exactly one part, which differs from the familiar vertex-bag tree-decompositions of graphs, but gives the same tree-width on graphs.

The tangle number θ(G)\theta(G)θ(G) is the maximum order of a tangle in GGG (p. 154), or 000 if there are none.

Formalization targets

Goal: (5.1)

max⁡(β(G),γ(G))  ≤  ω(G)+1  ≤  max⁡(⌊32β(G)⌋, γ(G), 1)for every finite hypergraph G.\max\bigl(\beta(G),\gamma(G)\bigr)\;\le\;\omega(G)+1\;\le\;\max\bigl(\lfloor\tfrac32\beta(G)\rfloor,\ \gamma(G),\ 1\bigr)\qquad\text{for every finite hypergraph } G.max(β(G),γ(G))≤ω(G)+1≤max(⌊23​β(G)⌋, γ(G), 1)for every finite hypergraph G.

Milestones

  • Second inequality of (5.1): ω(G)+1≤max⁡(⌊32β(G)⌋,γ(G),1)\omega(G)+1\le\max(\lfloor\tfrac32\beta(G)\rfloor,\gamma(G),1)ω(G)+1≤max(⌊23​β(G)⌋,γ(G),1).
  • Claims (1), (2), (3) of the proof of (5.1), under the proof's standing assumptions γ(G)>0\gamma(G)>0γ(G)>0, ∣E(G)∣≥2|E(G)|\ge 2∣E(G)∣≥2 and no isolated vertices: there is a tree-decomposition of width ω(G)\omega(G)ω(G) in which (1) every edge eee is the only edge of a leaf part whose vertices are exactly the ends of eee, and parts at vertices of valency ≥2\ge2≥2 have no edges; (2) moreover every leaf part has exactly one edge; (3) moreover every vertex of TTT has valency at most 333. The claims are cumulative.
  • First inequality of (5.1): max⁡(β(G),γ(G))≤ω(G)+1\max(\beta(G),\gamma(G))\le\omega(G)+1max(β(G),γ(G))≤ω(G)+1.

Companion: (5.2)

θ(G)  ≤  ω(G)+1  ≤  32 θ(G).\theta(G)\;\le\;\omega(G)+1\;\le\;\tfrac32\,\theta(G).θ(G)≤ω(G)+1≤23​θ(G).

Significance

The result. (5.1) shows that branch-width and tree-width bound each other up to the factor 3/23/23/2 and the edge-size correction, and the page records that both bounds are attained (KnK_nKn​ with 3∣n3\mid n3∣n for the right-hand side, Kn,nK_{n,n}Kn,n​ minus a perfect matching for the left-hand side). Together with the minimax theorem (4.3) of the same paper, max⁡(β(G),γ(G))=θ(G)\max(\beta(G),\gamma(G))=\theta(G)max(β(G),γ(G))=θ(G) unless γ(G)=0\gamma(G)=0γ(G)=0 and V(G)≠∅V(G)\ne\emptysetV(G)=∅, it yields (5.2): tangles are, up to the factor 3/23/23/2, the exact obstruction to tree-decompositions of small width. That duality is the form in which tree-width enters the later papers of the Graph Minors series.

Formalizing it. The result is proved (1991); no machine-checked proof of it is known to exist. The platform holds tree-width statements for simple graphs with vertex bags (Graph Minors V, RobertsonSeymour1986.GM5.treewidth_le_theta9 and excluding_planar_graph), which are a different theorem about a different object. This mission formalizes the hypergraph tree-decomposition and branch-decomposition of Graph Minors X and the comparison between them. The companion missions of this series are Graph Minors. X. Obstructions to Tree-Decomposition I (the minimax theorem (4.3)), III (the grid tangle (7.3)), IV (the tree-decomposition separating distinguishable tangles (10.3)) and V (the structure theorem (11.1)).

Difficulty

The second inequality converts a branch-decomposition into a tree-decomposition; the difficulty is to obtain the bound with the floor ⌊32β⌋\lfloor\tfrac32\beta\rfloor⌊23​β⌋ and the separate γ(G)\gamma(G)γ(G) term exactly, not up to an additive error.

The first inequality is where the obvious approach stalls: a tree-decomposition of minimum width is not a branch-decomposition. Its tree may have vertices of any valency, its leaves may carry no edge or several, and an edge may sit at an internal vertex. The three claims of the proof normalize an optimal tree-decomposition step by step without increasing its width. Each step modifies a tree on a finite vertex type while preserving three structural conditions (covering, edge-disjointness, the path condition), and the final step, suppressing vertices of valency 222, must be shown to produce a ternary tree with the edge orders bounded by part sizes. Most of the formal work is this tree surgery.

Formalization scope

  • Every theorem quantifies over types V, E with [Finite V] [Finite E] (the paper's "all hypergraphs in this paper are finite", p. 154) and a hypergraph G : Hypergraph V E given by its incidence relation.
  • Trees are SimpleGraph (Fin n) with IsTree (connected, hence non-null, and acyclic). A leaf is a vertex of valency at most 111; valency is (T.neighborSet t).ncard.
  • A branch-decomposition stores τ\tauτ as its inverse, an injective map from edges onto the leaves. β(G)\beta(G)β(G) is set to 000 when ∣E(G)∣≤1|E(G)|\le1∣E(G)∣≤1 and is otherwise the sInf of the widths, which is a genuine minimum because a branch-decomposition exists when ∣E(G)∣≥2|E(G)|\ge2∣E(G)∣≥2.
  • Widths of tree-decompositions are integers, so that ω(G)=−1\omega(G)=-1ω(G)=−1 when V(G)=∅V(G)=\emptysetV(G)=∅. ω(G)\omega(G)ω(G) is the sInf in Z\mathbb ZZ of all www bounding every ∣V(τ(t))∣−1|V(\tau(t))|-1∣V(τ(t))∣−1; the set is nonempty (one-vertex tree with τ=G\tau=Gτ=G) and bounded below by −1-1−1, so the sInf is the minimum width and not a junk value.
  • ⌊32β(G)⌋\lfloor\tfrac32\beta(G)\rfloor⌊23​β(G)⌋ is the natural-number division 3 * β / 2; the ceiling would be a weaker statement. In (5.2), ω(G)+1≤32θ(G)\omega(G)+1\le\tfrac32\theta(G)ω(G)+1≤23​θ(G) is stated as 2(ω(G)+1)≤3θ(G)2(\omega(G)+1)\le3\theta(G)2(ω(G)+1)≤3θ(G), which is equivalent.
  • The claims read "we may assume X" as "some tree-decomposition of width ω(G)\omega(G)ω(G) satisfies X", under the proof's standing assumptions γ(G)>0\gamma(G)>0γ(G)>0, ∣E(G)∣≥2|E(G)|\ge2∣E(G)∣≥2 and no isolated vertices (p. 168).
  • Condition (ii) of a tree-decomposition, edge-disjointness of the parts, is kept as on the page. Replacing it by "every edge lies in some part", or allowing an empty tree, would change the tree-width and trivialize the comparison; both are ruled out.

Infrastructure needed and reusable beyond this mission: tree surgery on SimpleGraph (Fin n) (attaching a leaf, deleting a leaf, splitting a vertex, suppressing valency-2 vertices), the existence of branch-decompositions for ∣E∣≥2|E|\ge2∣E∣≥2, and the computation of sInf-defined widths. Contributions of such lemmas, of proofs of the individual claims, and of small examples (KnK_nKn​, grids) are welcome.

Selected references

  • N. Robertson, P. D. Seymour, Graph Minors. X. Obstructions to Tree-Decomposition, J. Combin. Theory Ser. B 52 (1991) 153–190. https://doi.org/10.1016/0095-8956(91)90061-n
  • N. Robertson, P. D. Seymour, Graph Minors. V. Excluding a Planar Graph, J. Combin. Theory Ser. B 41 (1986) 92–114. https://doi.org/10.1016/0095-8956(86)90030-4
  • B. Courcelle, The monadic second-order logic of graphs. I. Recognizable sets of finite graphs, Information and Computation 85 (1990) 12–75. https://doi.org/10.1016/0890-5401(90)90043-H
10 thms1 active userReviewed
Graph TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

The Geometry of Graphs and Some of Its Algorithmic Applications 1: Every Multicommodity Flow Network with k Source-Sink Pairs Has a Cut S with Cap(S)/Dem(S) ≤ O(log k) · maxflowResearch Paper

Motivation

Routing several commodities through one capacitated network at the same time is a basic problem of network optimization, with applications to communication networks, VLSI layout and divide-and-conquer graph algorithms. Its optimum is the value of a linear program, but the natural certificates of an upper bound are combinatorial: every cut of the network limits how much can be routed across it. For a single commodity the max-flow min-cut theorem says that the best cut certificate is exact. With many commodities it is not, and the size of the gap between the best flow and the best cut is the question of this mission.

  • 1988–1999, Leighton and Rao proved that for uniform demands (one unit between every pair of vertices) the gap is O(log⁡n)O(\log n)O(logn) on an nnn-vertex network, and that the bound is attained on expanders (Leighton–Rao, JACM 1999; FOCS 1988).
  • 1990–1995, Klein, Rao, Agrawal and Ravi extended a polylogarithmic bound to arbitrary demands, O(log⁡Clog⁡D)O(\log C \log D)O(logClogD) in the total capacity CCC and total demand DDD (Combinatorica 15, 1995); Plotkin and Tardos improved it to O(log⁡2k)O(\log^2 k)O(log2k) in the number of commodities kkk (Combinatorica 15, 1995).
  • 1995, Linial, London and Rabinovich (Combinatorica 15, 1995), and independently Aumann and Rabani (SIAM J. Comput. 1998), showed that the gap is bounded by the least distortion of embedding a metric derived from the network into ℓ1\ell_1ℓ1​, and obtained the bound O(log⁡k)O(\log k)O(logk) from Bourgain's embedding theorem. This mission formalizes that theorem, Theorem 4.1 of the Linial–London–Rabinovich paper.

Setting

A multicommodity flow network has a finite vertex set VVV, capacities Ci,j=Cj,i≥0C_{i,j} = C_{j,i} \ge 0Ci,j​=Cj,i​≥0 (zero for non-edges and on the diagonal), and kkk commodities: commodity μ\muμ has a source sμs_\musμ​, a sink tμt_\mutμ​ and a demand Dμ≥0D_\mu \ge 0Dμ​≥0.

A feasible concurrent flow of value λ\lambdaλ sends, for every μ\muμ, λDμ\lambda D_\muλDμ​ units of commodity μ\muμ from sμs_\musμ​ to tμt_\mutμ​: arc flows fμ(i,j)≥0f_\mu(i,j) \ge 0fμ​(i,j)≥0 satisfy conservation of matter, with net outflow λDμ\lambda D_\muλDμ​ at sμs_\musμ​, −λDμ-\lambda D_\mu−λDμ​ at tμt_\mutμ​ and 000 elsewhere, and the total flow of all commodities through each undirected edge {i,j}\{i,j\}{i,j}, in both directions, is at most Ci,jC_{i,j}Ci,j​. The maxflow of the network is the largest such λ\lambdaλ.

For S⊆VS \subseteq VS⊆V, Cap(S)=∑i∈S∑j∉SCi,j\mathrm{Cap}(S) = \sum_{i \in S}\sum_{j \notin S} C_{i,j}Cap(S)=∑i∈S​∑j∈/S​Ci,j​ is the capacity of the edges leaving SSS, and Dem(S)\mathrm{Dem}(S)Dem(S) is the total demand of the pairs that SSS separates, i.e. with ∣S∩{sμ,tμ}∣=1|S \cap \{s_\mu,t_\mu\}| = 1∣S∩{sμ​,tμ​}∣=1. Every flow of value λ\lambdaλ has to push λ Dem(S)\lambda\,\mathrm{Dem}(S)λDem(S) units across the cut, so

maxflow≤Cap(S)Dem(S)for every S with Dem(S)>0.\mathrm{maxflow} \le \frac{\mathrm{Cap}(S)}{\mathrm{Dem}(S)} \quad\text{for every } S \text{ with } \mathrm{Dem}(S) > 0.maxflow≤Dem(S)Cap(S)​for every S with Dem(S)>0.

A semi-metric on VVV is a function d≥0d \ge 0d≥0 with d(i,i)=0d(i,i) = 0d(i,i)=0, symmetry and the triangle inequality; distinct points may be at distance 000, as in the paper.

Formalization targets

Goal: Theorem 4.1

There is an absolute constant C0>0C_0 > 0C0​>0 such that every network in which some commodity has Dμ>0D_\mu > 0Dμ​>0 and sμ≠tμs_\mu \ne t_\musμ​=tμ​ has a set S⊆VS \subseteq VS⊆V with Dem(S)>0\mathrm{Dem}(S) > 0Dem(S)>0 and

Cap(S)Dem(S)≤C0log⁡(max⁡(k,2))⋅maxflow.\frac{\mathrm{Cap}(S)}{\mathrm{Dem}(S)} \le C_0 \log(\max(k,2)) \cdot \mathrm{maxflow}.Dem(S)Cap(S)​≤C0​log(max(k,2))⋅maxflow.

The constant is not fixed; only the shape O(log⁡k)O(\log k)O(logk) in the number of commodities is asserted.

Milestones, in the order the proof uses them

  1. Weak duality (§4, p. 226): λ Dem(S)≤Cap(S)\lambda\,\mathrm{Dem}(S) \le \mathrm{Cap}(S)λDem(S)≤Cap(S) for every feasible λ≥0\lambda \ge 0λ≥0 and every SSS.
  2. LP duality display (p. 227): maxflow=min⁡d∑{i,j}Ci,jdi,j∑μDμdsμ,tμ\mathrm{maxflow} = \min_d \dfrac{\sum_{\{i,j\}} C_{i,j} d_{i,j}}{\sum_\mu D_\mu d_{s_\mu,t_\mu}}maxflow=mind​∑μ​Dμ​dsμ​,tμ​​∑{i,j}​Ci,j​di,j​​ over semi-metrics ddd.
  3. Corollary 3.4, Claim 2 (p. 223): every finite semi-metric space with kkk marked pairs has a non-expanding map into some ℓ1m\ell_1^mℓ1m​ that shrinks each marked pair by at most C1log⁡(max⁡(k,2))C_1 \log(\max(k,2))C1​log(max(k,2)).
  4. Coordinate minimum (p. 227): for points of ℓ1m\ell_1^mℓ1m​, the cost-to-demand ratio is at least its value on the best single coordinate.
  5. 0/1 rounding (pp. 227–228): for symmetric real weights aaa and nonnegative symmetric weights bbb, the minimum of ∑aij∣zi−zj∣/∑bij∣zi−zj∣\sum a_{ij}|z_i - z_j| / \sum b_{ij}|z_i - z_j|∑aij​∣zi​−zj​∣/∑bij​∣zi​−zj​∣ over real vectors zzz is attained at a 0/1 vector.

Significance

The theorem bounds the flow–cut gap: the integrality gap of the linear relaxation of the sparsest cut problem is O(log⁡k)O(\log k)O(logk), so solving one linear program and rounding gives an O(log⁡k)O(\log k)O(logk)-approximation to the sparsest cut. Sparsest-cut approximations are the engine of approximation algorithms for minimum bisection, graph partitioning, crossing number and VLSI layout. The dependence on kkk rather than on the number of vertices matters when few pairs carry demand. The bound is tight up to the constant: on constant-degree expanders with all-pairs unit demands the gap is Ω(log⁡n)\Omega(\log n)Ω(logn) (Leighton–Rao).

The result has been proved since 1995 and is textbook material (Williamson–Shmoys, Ch. 15). It has no machine-checked proof. A formalization adds a verified bridge between concurrent flows, LP duality on networks and ℓ1\ell_1ℓ1​ embeddings of finite metrics; each piece is reusable on its own. Corollary 3.4, Claim 2 in particular is a terminal version of Bourgain's theorem, of which no formal proof exists either.

Difficulty

Two steps carry the weight. The first is the duality between concurrent flows and metrics: Mathlib has no linear programming duality theorem in a form that applies directly, and the dual of the arc-flow program has to be identified with a minimum over semi-metrics, which needs a shortest-path argument on top of LP duality. The second is the embedding: Claim 2 of Corollary 3.4 must achieve distortion O(log⁡k)O(\log k)O(logk) on the kkk marked pairs, independent of ∣X∣|X|∣X∣. A direct application of Bourgain's theorem to the whole space gives O(log⁡∣X∣)O(\log |X|)O(log∣X∣), which is not enough when kkk is much smaller than the number of vertices, and the paper's proof is a one-line reference to an earlier argument.

Formalization scope

  • Vertex sets are finite types V with [Fintype V] [DecidableEq V]; capacities are a function C : V → V → ℝ with symmetry, nonnegativity and zero diagonal as fields of a Network structure. Commodities are indexed by Fin k; coincident or repeated pairs are allowed.
  • Flows are arc flows; the capacity bounds the total over all commodities and both directions. maxflow is the supremum of feasible values; Lean's sSup returns 000 on an unbounded set, which only happens without a commodity of positive demand and distinct endpoints, and the goal assumes one.
  • Sums over edges are over unordered pairs: the page's ∑i≠j\sum_{i \ne j}∑i=j​ in the LP-duality and coordinate displays counts each edge twice, while Cap\mathrm{Cap}Cap counts it once; the statements use 12∑i∑j\tfrac12\sum_i\sum_j21​∑i​∑j​ in both places, which removes the factor-2 slip without changing the proof.
  • Explicit constants. The paper's O(log⁡k)O(\log k)O(logk) in Theorem 4.1 is C0log⁡(max⁡(k,2))C_0 \log(\max(k,2))C0​log(max(k,2)) with C0C_0C0​ quantified before the network; the Ω(1/log⁡k)\Omega(1/\log k)Ω(1/logk) of Corollary 3.4 is 1/(C1log⁡(max⁡(k,2)))1/(C_1\log(\max(k,2)))1/(C1​log(max(k,2))) with C1C_1C1​ quantified before the metric space. The max⁡(k,2)\max(k,2)max(k,2) keeps k=1k = 1k=1 in scope, where log⁡k=0\log k = 0logk=0.
  • Out of scope. The deterministic polynomial-time algorithms of Theorem 4.1 and Corollary 3.4, and the target dimension ℓ1O(n2)\ell_1^{O(n^2)}ℓ1O(n2)​ of the latter, are not stated; the statements assert the existence of the cut and of the embedding.
  • The 0/1 rounding step adds b≥0b \ge 0b≥0, which the page omits but its argument needs.
  • A trivializing formalization is ruled out: the constant C0C_0C0​ is quantified before the network, so it cannot depend on the instance, and maxflow is defined from flows, not as the metric minimum.
  • Contributions welcome: an LP duality theorem for finite-dimensional linear programs in the form needed here, a formal Bourgain-type embedding with terminal distortion, and the max-flow min-cut theorem for undirected capacities.

Selected references

  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15 (1995) 215–245. https://doi.org/10.1007/BF01200757
  • Y. Aumann, Y. Rabani, An O(log k) approximate min-cut max-flow theorem and approximation algorithm, SIAM J. Comput. 27 (1998) 291–301. https://doi.org/10.1137/S0097539794285983
  • T. Leighton, S. Rao, Multicommodity max-flow min-cut theorems and their use in designing approximation algorithms, J. ACM 46 (1999) 787–832. https://doi.org/10.1145/331524.331526
  • J. Bourgain, On Lipschitz embedding of finite metric spaces in Hilbert space, Israel J. Math. 52 (1985) 46–52. https://doi.org/10.1007/BF02776078
  • D. P. Williamson, D. B. Shmoys, The Design of Approximation Algorithms, Cambridge University Press, 2011, Chapter 15. https://doi.org/10.1017/CBO9780511921735
9 thms1 active userReviewed
Operations ResearchProbabilityTheoretical Computer Science·Captain: mikedeng1

Packet Routing and Job-Shop Scheduling in O(Congestion + Dilation) Steps: Edge-Simple Paths with Congestion c and Dilation d Admit an O(c + d)-Step Schedule with Constant-Size QueuesResearch Paper

Motivation

In a store-and-forward network, messages are cut into packets that travel from node to node along wires, one wire per time step, and wait in buffers between moves. Routing such traffic splits into two problems: choosing a path for each packet, and scheduling the packets along their paths, deciding at every step which packets move and which wait. Leighton, Maggs and Rao showed that the second problem always has an essentially optimal solution: once paths are fixed, two simple parameters of the paths determine the routing time up to a constant factor, on every network.

The result separates path selection from timing in the design of routing algorithms for parallel machines, and it reaches beyond networks: a job-shop problem in which every operation takes one unit of time and no job visits a machine twice is the same problem with jobs as packets and machines as edges.

Timeline.

  • 1988: Leighton, Maggs and Rao, extended abstract at FOCS, Universal packet routing algorithms.
  • 1994: the full paper in Combinatorica 14, with the existence theorem of an O(c+d)O(c+d)O(c+d) schedule with constant queues (Theorem 3.4) and a randomized on-line algorithm.
  • 1999: Leighton, Maggs and Richa, Fast algorithms for finding O(congestion + dilation) packet routing schedules, make the construction algorithmic, using the algorithmic Local Lemma of Beck.

Setting

A network is a directed multigraph: a type VVV of nodes, a type EEE of edges, and maps src,tgt:E→V\mathrm{src},\mathrm{tgt}:E\to Vsrc,tgt:E→V. A path is a list of edges e0,…,eℓ−1e_0,\dots,e_{\ell-1}e0​,…,eℓ−1​ with tgt(ek)=src(ek+1)\mathrm{tgt}(e_k)=\mathrm{src}(e_{k+1})tgt(ek​)=src(ek+1​); it is edge-simple if no edge occurs twice. A finite set PPP of packets is given, each with its path.

The dilation ddd is the largest number of edges on a path; the congestion ccc is the largest number of paths through one edge. Since a packet crosses at most one edge per step and an edge carries at most one packet per step, every schedule needs at least max⁡(c,d)\max(c,d)max(c,d) steps.

A schedule assigns to every packet ppp and every index kkk of its path the step τ(p,k)≥1\tau(p,k)\ge1τ(p,k)≥1 at which ppp crosses its kkk-th edge, strictly increasing in kkk. Its length is the last crossing step. A packet waits in its initial queue before its first crossing, in the edge queue at the head of the edge it last crossed between two crossings, and in its final queue after its last crossing. Only edge queues count for queue size: the other two are fixed by the instance. A schedule is valid if at most one packet crosses each edge at each step.

For the proof, the paper measures schedules that are not yet valid through frames: a TTT-frame is a run of TTT consecutive steps, and the relative congestion of a frame is the largest number of packets crossing one edge in it, divided by TTT.

Formalization targets

Goal: Theorem 3.4 (p. 11)

There are absolute constants KKK and QQQ such that every finite set of packets with edge-simple paths of congestion at most ccc and dilation at most ddd, on any network, has a valid schedule with

length≤K (c+d),every edge queue≤Q at every step.\text{length}\le K\,(c+d),\qquad \text{every edge queue}\le Q \text{ at every step}.length≤K(c+d),every edge queue≤Q at every step.

The constants are not fixed: the paper proves existence, and any constants are a valid answer.

Milestones, in the order the proof uses them

  1. Lemma 3.1 (p. 8), the Lovász Local Lemma: events of probability at most ppp with dependence at most b≥1b\ge1b≥1 and 4pb<14pb<14pb<1 all fail together with positive probability.
  2. Lemma 3.2 (p. 8): with congestion and dilation at most ddd, a schedule of length O(d)O(d)O(d) with no waiting in edge queues and at most TTT packets per edge in every frame of size T≥log⁡2dT\ge\log_2 dT≥log2​d.
  3. The recurrences (pp. 12–13): I(1)=log⁡dI^{(1)}=\log dI(1)=logd, I(i+1)=log⁡5I(i)I^{(i+1)}=\log^5 I^{(i)}I(i+1)=log5I(i), r(1)=1r^{(1)}=1r(1)=1, r(i+1)=r(i)(1+κ/log⁡I(i))r^{(i+1)}=r^{(i)}(1+\kappa/\sqrt{\log I^{(i)}})r(i+1)=r(i)(1+κ/logI(i)​) stop at some j=O(log⁡∗d)j=O(\log^* d)j=O(log∗d) with r(j)=O(1)r^{(j)}=O(1)r(j)=O(1).
  4. Lemma 3.5 (p. 14): frame bounds for all sizes from TTT to 2T−12T-12T−1 imply them for all sizes ≥T\ge T≥T.
  5. The final simulation (pp. 12–13): a schedule with relative congestion O(1)O(1)O(1) in frames of constant size, in which every packet waits at most once every k1≥2k_1\ge2k1​≥2 steps, becomes a valid schedule, a constant factor longer, with constant edge queues.

Significance

The theorem shows that congestion and dilation, two quantities read off the paths alone, determine the optimal schedule length up to a constant factor on every network, with buffers of constant size. Path-selection algorithms that minimize c+dc+dc+d therefore yield near-optimal routing, and the same bound holds for unit-time job shops without repeated machines. It is also an often-cited application of the Local Lemma beyond a single round of random choices: the proof applies it O(log⁡∗d)O(\log^* d)O(log∗d) times in succession.

The result has been proved since 1994, and the algorithmic version since 1999. As far as is known, no part of it is machine-checked. The mission produces a formal model of store-and-forward schedules with edge queues and frames, a formal statement of the theorem that excludes the degenerate readings, and formal versions of the steps of its proof.

Difficulty

The naive approach gives each packet a random initial delay and then lets it move without waiting. That gives O(log⁡(Nd))O(\log(Nd))O(log(Nd)) packets per edge per step, and O(c+dlog⁡(Nd))O(c+d\log(Nd))O(c+dlog(Nd)) steps after slowing down. Lemma 3.2 does better only in frames of size log⁡d\log dlogd, not in single steps. Recursing on frames, as in Theorem 3.3, loses a constant factor per level and gives (c+d)2O(log⁡∗(c+d))(c+d)2^{O(\log^* (c+d))}(c+d)2O(log∗(c+d)).

Removing that factor is the central difficulty. Each refinement must keep the relative congestion nearly unchanged, r(i+1)=r(i)(1+O(1)/log⁡I(i))r^{(i+1)}=r^{(i)}(1+O(1)/\sqrt{\log I^{(i)}})r(i+1)=r(i)(1+O(1)/logI(i)​), which requires second-order terms in the tail estimates, delays spread over the block rather than inserted at its start, and careful handling of block boundaries. The constant queue bound needs an invariant: every packet waits at most once every I(i)I^{(i)}I(i) steps. The Local Lemma gives existence only; the construction is non-constructive.

Formalization scope

Conventions committed to in Lean:

  • The network is arbitrary: V E : Type with src tgt : E → V, no finiteness or degree bound. Packets form a Fintype P; paths are List E with matching endpoints (List.IsChain) and List.Nodup.
  • Congestion and dilation are bounds (CongestionLE path c, DilationLE path d), equivalent to exact values because every bound is monotone.
  • A schedule is a Timetable: time p k : ℕ, the step at which packet p crosses its k-th edge, at least 1 and strictly increasing in k. Length at most L means every crossing step is ≤ L. Valid means two different crossings of one edge happen at different steps.
  • The edge-queue size of edge g at the end of step t counts crossings of g at a step ≤ t whose packet crosses its next edge at a step > t. Initial and final queues are not counted.
  • Frames: frameCount τ g t T counts the packets that cross g at a step in [t, t+T). Relative congestion at most r in frames of size ≥ T₀ is frameCount ≤ r·T for all T ≥ 1 with T₀ ≤ T.
  • log is Real.logb 2. Every O(1) and "sufficiently large" is a constant quantified before the instance, and in the goal ∃ K Q comes before the network.
  • Lemma 3.1 is stated on an arbitrary probability space with measurable events. Dependence uses independence from the generated σ-algebra, and the statement adds 1 ≤ b, without which the page's statement is false.

These choices rule out the trivial formalizations: constants chosen after the instance would make K=cdK=cdK=cd suffice; a schedule without edge exclusivity would make the greedy schedule of length ddd a solution; a missing queue bound drops half of the theorem; and a statement without its length bound is solved by sending one packet at a time.

Not stated: Lemmas 3.6–3.10 and the summary of the refinement step (p. 20). They concern the block decomposition and delay-insertion rules of pp. 13–14 and 17–18, which the paper defines only in prose. Contributions are welcome on the Local Lemma (finite or general), the probabilistic estimates of Lemma 3.2, Lemma 3.5, the elementary simulation step, and formal definitions of the block operations from which Lemmas 3.6–3.10 can be stated. The schedule model is reusable for other routing results, such as the on-line algorithm of §2 or the O(c+d)O(c+d)O(c+d) results for leveled networks.

Selected references

  • F. T. Leighton, B. M. Maggs, S. B. Rao, Packet routing and job-shop scheduling in O(congestion + dilation) steps, Combinatorica 14 (1994) 167–186. https://doi.org/10.1007/BF01215349 (this mission cites the authors' manuscript).
  • F. T. Leighton, B. M. Maggs, S. B. Rao, Universal packet routing algorithms, Proc. 29th IEEE FOCS (1988) 256–269.
  • F. T. Leighton, B. M. Maggs, A. W. Richa, Fast algorithms for finding O(congestion + dilation) packet routing schedules, Combinatorica 19 (1999) 375–401. https://doi.org/10.1007/s004930050061
  • P. Erdős, L. Lovász, Problems and results on 3-chromatic hypergraphs and some related questions, in Infinite and Finite Sets, Colloq. Math. Soc. János Bolyai 10 (1975) 609–627.
  • J. Spencer, Ten Lectures on the Probabilistic Method, SIAM (1987), pp. 57–58. https://doi.org/10.1137/1.9780898719918
10 thms1 active userReviewed
Operations ResearchOptimization·Captain: mikedeng1

Optimization and Approximation in Deterministic Sequencing and Scheduling: A Survey 2: The Optimal Preemptive Open Shop Makespan Equals the Largest Machine Load or Job LengthResearch Paper

Motivation

Open shops model production and service systems in which every job must visit every machine, but the order of the visits is free: a car that needs an inspection, a wash and a tyre change, a patient who needs several tests, a student who sits several exams. The survey of Graham, Lawler, Lenstra and Rinnooy Kan (Ann. Discrete Math. 5, 1979) fixed the three-field notation α∣β∣γ\alpha|\beta|\gammaα∣β∣γ that the scheduling literature still uses, and classified the complexity of the problems it can express. Among the polynomially solvable cases, the preemptive open shop with makespan objective, O∣pmtn∣Cmax⁡O|pmtn|C_{\max}O∣pmtn∣Cmax​, is one of the few multi-machine problems whose optimal value has a closed form for any number of machines and jobs.

Timeline.

  • 1976: Gonzalez and Sahni (J. ACM 23) prove that the optimal preemptive open-shop makespan is the largest machine load or job length, and give a polynomial algorithm. In the same paper they solve O2∥Cmax⁡O2\|C_{\max}O2∥Cmax​ in linear time and show O3∥Cmax⁡O3\|C_{\max}O3∥Cmax​ NP-hard.
  • 1978: Lawler and Labetoulle (J. ACM 25) give a linear-programming treatment of preemptive scheduling on unrelated machines and reformulate the open-shop construction in terms of decrementing sets, found by an assignment problem through the Birkhoff–von Neumann theorem.
  • 1979: the survey (§5.2.2, p. 313) presents this construction as the standard argument and records the O(r+min⁡{m4,n4,r2})O(r+\min\{m^4,n^4,r^2\})O(r+min{m4,n4,r2}) bound of Gonzalez (1976), where rrr is the number of nonzero processing times.

Setting

There are mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​ and nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​. Job JjJ_jJj​ consists of operations O1j,…,OmjO_{1j},\dots,O_{mj}O1j​,…,Omj​; operation OijO_{ij}Oij​ must be processed on machine MiM_iMi​ for pij≥0p_{ij}\ge 0pij​≥0 time units. The processing-time matrix is P=(pij)P=(p_{ij})P=(pij​): its rows are machines and its columns are jobs. Every job is available at time 000.

Preemption is allowed: an operation may be interrupted and resumed later. A schedule is a finite list of pieces (i,j,s,e)(i,j,s,e)(i,j,s,e), each meaning that MiM_iMi​ processes JjJ_jJj​ during [s,e)[s,e)[s,e). A schedule is feasible if

  1. every piece satisfies 0≤s≤e0\le s\le e0≤s≤e;
  2. each machine processes at most one job at a time, and each job is processed on at most one machine at a time: two pieces that share a machine or a job do not overlap;
  3. for every pair (i,j)(i,j)(i,j) the pieces of OijO_{ij}Oij​ have total length exactly pijp_{ij}pij​.

The makespan Cmax⁡C_{\max}Cmax​ is the time at which the last piece ends, and Cmax⁡∗C^*_{\max}Cmax∗​ is its minimum over feasible schedules. The load of machine MiM_iMi​ is the row sum ∑jpij\sum_j p_{ij}∑j​pij​, the length of job JjJ_jJj​ is the column sum ∑ipij\sum_i p_{ij}∑i​pij​, and

C=max⁡{max⁡j∑ipij, max⁡i∑jpij}.C=\max\Big\{\max_j \sum_i p_{ij},\ \max_i \sum_j p_{ij}\Big\}.C=max{jmax​i∑​pij​, imax​j∑​pij​}.

A row or column is tight if its sum equals CCC and slack otherwise. A decrementing set is a set SSS of strictly positive entries of PPP with exactly one element in each tight row and each tight column and at most one in each slack row and each slack column.

Formalization targets

Goal: Cmax⁡∗=CC^*_{\max}=CCmax∗​=C

For every T≥0T\ge 0T≥0,

(∃ feasible schedule with Cmax⁡≤T)  ⟺  (∑jpij≤T ∀i  and  ∑ipij≤T ∀j).\big(\exists \text{ feasible schedule with } C_{\max}\le T\big)\iff \Big(\sum_j p_{ij}\le T\ \forall i\ \text{ and }\ \sum_i p_{ij}\le T\ \forall j\Big).(∃ feasible schedule with Cmax​≤T)⟺(j∑​pij​≤T ∀i  and  i∑​pij​≤T ∀j).

This says that the optimal makespan is exactly the largest machine load or job length, and that it is attained.

Milestones (all from §5.2.2, p. 313)

  1. Lower bound Cmax⁡∗≥CC^*_{\max}\ge CCmax∗​≥C.
  2. Existence of a decrementing set for every nonzero nonnegative PPP.
  3. Positive step: for a decrementing set, the largest δ\deltaδ satisfying the constraints (1)–(3) of the survey exists and is positive.
  4. Step property: after replacing each pij∈Sp_{ij}\in Spij​∈S by max⁡{0,pij−δ}\max\{0,p_{ij}-\delta\}max{0,pij​−δ}, the largest line sum is exactly C−δC-\deltaC−δ.
  5. Partial schedule: for each pij∈Sp_{ij}\in Spij​∈S, MiM_iMi​ processes JjJ_jJj​ for min⁡{pij,δ}\min\{p_{ij},\delta\}min{pij​,δ} time units, with no machine or job used twice.
  6. Termination: every run of the procedure reaches P′=(0)P'=(0)P′=(0) within a bounded number of stages.
  7. Joining: the concatenated partial schedules form a feasible schedule with Cmax⁡≤CC_{\max}\le CCmax​≤C.

Significance

The result. The theorem turns an optimization over continuous-time schedules into the computation of m+nm+nm+n sums. It certifies optimality by a counting argument, it is the base case for preemptive open shops with release dates and due dates, and it is used elsewhere in the survey (§4.4.6) to reduce problems on unrelated machines with preemption to open-shop instances. Because a nonnegative matrix whose row and column sums are all equal is a multiple of a doubly stochastic matrix, the theorem is a scheduling form of the Birkhoff–von Neumann decomposition. It also underlies timetabling and edge-colouring results for bipartite multigraphs.

Formalizing it. The theorem has been proved since 1976 and is textbook material. No machine-checked proof is known to exist: Mathlib has the Birkhoff–von Neumann theorem for doubly stochastic matrices but no model of open-shop schedules. A formalization adds a reusable model of preemptive multi-machine schedules with both disjointness requirements, a checked proof of the decrementing-set construction, and a termination argument the survey asserts without proof.

Difficulty

The lower bound is a one-line counting argument. The difficulty is the construction of a schedule of length exactly CCC. Scheduling each machine's operations back to back gives length max⁡i∑jpij\max_i\sum_j p_{ij}maxi​∑j​pij​, but may run one job on two machines at once. Scheduling job by job has the symmetric defect. A greedy list schedule that only respects both constraints can leave machines idle and overshoot CCC. The construction must keep every tight line busy at every moment while never letting a slack line fall behind. The existence of the decrementing set at each stage is the combinatorial core: it is a Hall-type matching condition, not a local choice. Termination is also not automatic, because a careless choice of step length can produce infinitely many shrinking steps.

Formalization scope

All objects live in the namespace SchedSurvey.OPmtn. Machines and jobs are Fin m and Fin n, both 0-based, and processing times and piece endpoints are real numbers; integer data are a special case. A schedule is a List of pieces. Feasibility requires nonnegative start times, disjointness for pieces sharing a machine or a job (touching intervals allowed), and exactly pijp_{ij}pij​ units of processing for every pair (i,j)(i,j)(i,j). Every theorem assumes pij≥0p_{ij}\ge 0pij​≥0.

"CCC is the maximum" is the predicate IsMaxLoad P C: all line sums are at most CCC and one equals CCC. The goal is stated in threshold form and mentions no maximum at all. The display defining CCC on p. 313 prints max⁡i{∑ipij}\max_i\{\sum_i p_{ij}\}maxi​{∑i​pij​} for the second term; the following sentence shows it means the row sums ∑jpij\sum_j p_{ij}∑j​pij​, and the formalization uses those. The hypothesis T≥0T\ge 0T≥0 matters only for P=0P=0P=0, where the empty schedule finishes by every TTT.

A trivializing formalization is ruled out: the goal mentions neither decrementing sets nor δ\deltaδ. Its "if" direction asserts that a schedule exists. Feasibility counts work per pair (machine, job), not per job, and forbids a job from running on two machines at once. Without either requirement the statement would be a different and easier theorem.

A complete development needs: list sums of interval lengths over disjoint intervals, the existence of decrementing sets (via Birkhoff–von Neumann, König's theorem or Hall's theorem on the bipartite graph of positive entries), the step and termination lemmas, and concatenation of schedules. The schedule model and the decrementing-set lemma are reusable for other preemptive shop problems. Contributions of alternative proofs of any milestone, for example a direct Hall-theorem proof of the existence of decrementing sets, are welcome.

Selected references

  • 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
  • T. Gonzalez, S. Sahni, Open shop scheduling to minimize finish time, Journal of the ACM 23 (1976) 665–679. https://doi.org/10.1145/321978.321985
  • E. L. Lawler, J. Labetoulle, On preemptive scheduling of unrelated parallel processors by linear programming, Journal of the ACM 25 (1978) 612–619. https://doi.org/10.1145/322077.322090
9 thms1 active userReviewed
Dynamic ProgrammingGraph TheoryOperations Research·Captain: mikedeng1

The Complexity of Markov Decision Processes IV: In a Deterministic Process a Cheapest T-Arc Walk Is a Walk of Fewer Than n³ Arcs Plus One Repeated CycleResearch Paper

Motivation

Finite-horizon Markov decision processes ask which decisions minimize accumulated cost over a specified number of stages. In a stationary process, the same costs and transitions apply at every stage. This compact input can describe a horizon TTT much larger than the state set, so a method that processes all TTT stages one by one may take time exponential in the length of the written horizon. Papadimitriou and Tsitsiklis identified this issue while comparing sequential and parallel algorithms for decision processes. For the deterministic stationary case, their Theorem 5 states an NC classification, supported by a finite structural description of a cheapest long walk. Papadimitriou and Tsitsiklis (1987)

The structural description is the subject of this mission. It connects a policy optimization problem to bounded-length graph objects and an exact cost identity. This is useful independently of a particular model of parallel computation: it says which shorter walk and cycle costs determine the answer even when the horizon is very large. Papadimitriou and Tsitsiklis (1987), p. 447

Setting

Let SSS be a finite nonempty set of states, with n=∣S∣n=|S|n=∣S∣, and fix an initial state s0s_0s0​. At each state sss, a finite nonempty set DsD_sDs​ contains the available decisions. A decision d∈Dsd\in D_sd∈Ds​ has a certain successor state next⁡(s,d)\operatorname{next}(s,d)next(s,d) and a real cost c(s,d)c(s,d)c(s,d). The model is stationary: these data do not depend on time. The corresponding directed graph has one arc for each decision. It may have loops and parallel arcs, because distinct decisions can share endpoints. This is the deterministic graph translation made in §3 of the paper. Papadimitriou and Tsitsiklis (1987), pp. 444–445

A policy is a function δ(s,t)∈Ds\delta(s,t)\in D_sδ(s,t)∈Ds​ of both state and time. Starting at s0s_0s0​, it produces st+1=next⁡(st,δ(st,t))s_{t+1}=\operatorname{next}(s_t,\delta(s_t,t))st+1​=next(st​,δ(st​,t)). In the convention of the paragraph before Theorem 5, horizon TTT means TTT decisions and therefore a walk with TTT arcs. Its cost is

JT(s0,δ)=∑t=0T−1c(st,δ(st,t)),JT∗(s0)=inf⁡δJT(s0,δ).J_T(s_0,\delta)=\sum_{t=0}^{T-1}c(s_t,\delta(s_t,t)),\qquad J_T^*(s_0)=\inf_\delta J_T(s_0,\delta).JT​(s0​,δ)=t=0∑T−1​c(st​,δ(st​,t)),JT∗​(s0​)=δinf​JT​(s0​,δ).

A walk records the states and the decision arc used at each step; states may repeat. A simple path has no repeated state. A simple cycle has positive length, returns to its starting state, and has no other repeated state. A closed walk may repeat states internally and need not be a simple cycle. The distinction matters because the structural claim uses a simple cycle, whereas the computation in the target identity uses a cheapest closed walk of a chosen length.

Formalization targets

The first milestone identifies JT∗(s0)J_T^*(s_0)JT∗​(s0​) with the minimum cost of a TTT-arc walk from s0s_0s0​. The next milestones record a decomposition into a simple path and simple cycles, the restriction on cycles repeated more than nnn times, and the existence of a cheapest walk with a short prefix and one repeated simple cycle. These are claims in the paragraph preceding Theorem 5, with corrections to the printed exchange and its connectivity justification documented in the formal statements. Papadimitriou and Tsitsiklis (1987), p. 447

For l≤Tl\le Tl≤T and u∈Su\in Su∈S, let wl,uw_{l,u}wl,u​ be the cheapest lll-arc walk from s0s_0s0​ that visits uuu. Let rk,ur_{k,u}rk,u​ be the cheapest closed kkk-arc walk at uuu. Only choices for which these walks exist are included. The mission's goal is

JT∗(s0)=min⁡l≤T, l<n3, u∈S,1≤k≤n, k∣T−l(wl,u+T−lk rk,u),T>n2.J_T^*(s_0)=\min_{\substack{l\le T,\ l<n^3,\ u\in S,\\1\le k\le n,\ k\mid T-l}}\left(w_{l,u}+\frac{T-l}{k}\,r_{k,u}\right),\qquad T>n^2.JT∗​(s0​)=l≤T, l<n3, u∈S,1≤k≤n, k∣T−l​min​(wl,u​+kT−l​rk,u​),T>n2.

The l<n3l<n^3l<n3 bound is the paper's finite compression claim. The divisibility condition makes (T−l)/k(T-l)/k(T−l)/k an integer number of cycle copies. The goal states an exact value, with the inner minima taken over genuine walks and the outer optimum still taken over all policies. Papadimitriou and Tsitsiklis (1987), p. 447

Significance

The identity replaces an optimization over arbitrarily many time-dependent policy decisions with costs of walks having fewer than n3n^3n3 arcs and closed walks of at most nnn arcs. It gives the mathematical statement that must hold before the paper's proposed parallel calculation can be trusted. It also applies when several different decisions have the same endpoints, because each decision remains a separate weighted arc. Papadimitriou and Tsitsiklis (1987), pp. 445, 447

The paper proves the result informally and states the NC classification. The mission seeks machine-checked proofs of the graph and cost statements, including the printed argument's delicate boundary cases. The local proposal contains open Lean theorem statements, not proofs. Reusable outcomes include a decision-walk representation, a cycle decomposition that preserves the multiset of arcs, and finite-horizon policy-to-walk equivalence. A related published result, TropicalLA.exists_repeat, detects repetition for dense max-plus matrices; it does not state this weighted multigraph decomposition or the min-cost identity.

Difficulty

The obstacle is preserving both walk feasibility and the exact number of arcs when repeated cycles are rearranged. A cycle may be the only connection to other parts of a walk, so deleting its last copy can invalidate the remaining sequence. The printed exchange paragraph also reverses the cost-improving direction and appeals to a no-ties perturbation without carrying it through the bounded-prefix conclusion. Those issues affect the claimed n3n^3n3 bound; merely knowing that a long walk repeats a vertex does not establish the target identity. Papadimitriou and Tsitsiklis (1987), p. 447

Formalization scope

Lean represents SSS and each DsD_sDs​ as finite types and uses real costs. The decision sets are nonempty, making policies and arbitrarily long walks available. Policies depend on state and time, as in §2. A walk stores its decisions as well as its states, so loops and parallel arcs are represented. Cycle decompositions preserve endpoints, length, cost, and the multiset of decision arcs. A loop is a simple cycle of length one. Time starts at zero, and an lll-arc walk has states at indices 0,…,l0,\ldots,l0,…,l.

The paper's general finite-horizon definition sums costs from 000 through TTT, while the Theorem 5 paragraph explicitly asks for a TTT-arc walk. This mission follows the latter convention: decisions occur at 0,…,T−10,\ldots,T-10,…,T−1. The goal also imposes l≤Tl\le Tl≤T, which the printed loop leaves implicit, and uses a closed walk for rk,ur_{k,u}rk,u​ because that is what the length-kkk graph calculation computes. It does not assume distinct walk costs. The condition T>n2T>n^2T>n2 is retained from the paper even though related identities may hold in other ranges. Papadimitriou and Tsitsiklis (1987), pp. 444, 447

Theorem 5's NC membership, processor and parallel-time bounds, the membership-in-P discussion, log-space reductions, PSPACE membership, and Corollaries 1–2 are outside this mission. Mathlib does not supply the PRAM and complexity-class framework those claims require. The target is the exact cost identity behind the algorithm. Defining the optimum directly as the right-hand minimum, restricting policies to stationary policies, or allowing an empty policy class would trivialize or change it. Contributions proving the milestones, the goal, or reusable walk and cycle facts are welcome.

Selected references

  • Christos H. Papadimitriou and John N. Tsitsiklis, The Complexity of Markov Decision Processes, Mathematics of Operations Research 12(3), 441–450, 1987. DOI: 10.1287/moor.12.3.441
6 thms1 active userReviewed
Algorithmic Game TheoryOperations ResearchOptimization+1·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments V: Stackelberg Matroid Pricing Is an Assortment Problem Under a Regular Choice ModelResearch Paper

Motivation

In Stackelberg network pricing a leader sets prices on some resources, and a follower then buys the cheapest structure available to them, paying the leader for the priced resources they use. The Stackelberg Minimum Spanning Tree problem, introduced by Cardinal, Demaine, Fiorini, Joret, Langerman, Newman and Weimann (Algorithmica 2011), is the version in which the follower buys a minimum spanning tree. The best known approximation factors for it are those of uniform pricing, which gives every priced edge the same price (Berbeglia–Joret, p. 21).

Independently, revenue management studies assortment optimisation: a seller chooses which products to offer, and customers choose among the offered products according to a discrete choice model. Berbeglia and Joret (arXiv:1606.01371, Algorithmica 2020) analyse the revenue-ordered heuristic (offer every product whose revenue is at least some threshold) under any regular choice model, and prove tight approximation bounds.

§4.6 of that paper shows that the two problems are the same problem: a Stackelberg Matroid pricing instance is an assortment problem under a regular choice model, and uniform pricing is the revenue-ordered heuristic on that model. The bounds of Cardinal et al. on uniform pricing thus become special cases of the general bounds on revenue-ordered assortments. This mission formalizes that correspondence (Theorem 4.16).

Setting

Choice model. Let C\mathcal CC be a finite set of products. A system of choice probabilities assigns to each offered set S⊆CS\subseteq\mathcal CS⊆C and product xxx a probability P(x,S)\mathcal P(x,S)P(x,S), with no-purchase probability P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). It is regular if (i) all probabilities are nonnegative, (ii) P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S, (iii) ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1, and (iv) P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) whenever S⊆S′S\subseteq S'S⊆S′ and x∈S∪{0}x\in S\cup\{0\}x∈S∪{0}. With revenues r:C→R>0r:\mathcal C\to\mathbb R_{>0}r:C→R>0​, rev(S)=∑x∈SP(x,S)r(x)\mathrm{rev}(S)=\sum_{x\in S}\mathcal P(x,S)r(x)rev(S)=∑x∈S​P(x,S)r(x) and OPT=max⁡Srev(S)\mathrm{OPT}=\max_S\mathrm{rev}(S)OPT=maxS​rev(S). The revenue-ordered assortment generated by yyy is {y′∈C:r(y′)≥r(y)}\{y'\in\mathcal C:r(y')\ge r(y)\}{y′∈C:r(y′)≥r(y)}.

Greedy algorithm. Given a family of independent sets, a linear ordering LLL and a set FFF, greedy(F,L)\mathrm{greedy}(F,L)greedy(F,L) scans FFF in the order induced by LLL, starting from ∅\emptyset∅, and adds each element that keeps the current set independent.

Stackelberg Matroid problem. An instance is a matroid M=(E,X)M=(E,\mathcal X)M=(E,X), a bipartition E=R⊔BE=R\sqcup BE=R⊔B into red and blue elements, and red costs c:R→R>0c:R\to\mathbb R_{>0}c:R→R>0​; some base of MMM lies inside RRR. The leader chooses prices p:B→R>0p:B\to\mathbb R_{>0}p:B→R>0​. The customer buys a minimum-weight base by running greedy on R∪BR\cup BR∪B with an ordering L∗L^*L∗ that is non-decreasing in weight (ccc on RRR, ppp on BBB) and puts blue elements first on ties. The leader earns revStack(p,L∗)=∑e∈B∩greedyM(R∪B,L∗)p(e)\mathrm{rev}_{\mathrm{Stack}}(p,L^*)=\sum_{e\in B\cap\mathrm{greedy}_M(R\cup B,L^*)}p(e)revStack​(p,L∗)=∑e∈B∩greedyM​(R∪B,L∗)​p(e).

The assortment instance. Let c1<⋯<ckc_1<\dots<c_kc1​<⋯<ck​ be the distinct red costs. Products are C=B×{c1,…,ck}\mathcal C=B\times\{c_1,\dots,c_k\}C=B×{c1​,…,ck​} with r((e,q))=∣B∣ qr((e,q))=|B|\,qr((e,q))=∣B∣q. The auxiliary matroid M′M'M′ on R∪CR\cup\mathcal CR∪C declares XXX independent when it holds at most one pair (e,q)(e,q)(e,q) per blue eee and (R∩X)∪{e:(e,q)∈X}(R\cap X)\cup\{e:(e,q)\in X\}(R∩X)∪{e:(e,q)∈X} is independent in MMM. With an ordering LLL of R∪CR\cup\mathcal CR∪C that is non-decreasing in cost and puts products before red elements on ties, P((e,q),S)=1/∣B∣\mathcal P((e,q),S)=1/|B|P((e,q),S)=1/∣B∣ if (e,q)∈greedyM′(R∪S,L)(e,q)\in\mathrm{greedy}_{M'}(R\cup S,L)(e,q)∈greedyM′​(R∪S,L) and 000 otherwise.

Formalization targets

Goal: Theorem 4.16 (p. 24)

For the instance above, with B≠∅B\ne\emptysetB=∅ and every admissible LLL:

P is regular,r>0,max⁡S⊆Crev(S)=max⁡p>0, L∗revStack(p,L∗),\mathcal P\ \text{is regular},\quad r>0,\qquad \max_{S\subseteq\mathcal C}\mathrm{rev}(S)=\max_{p>0,\ L^*}\mathrm{rev}_{\mathrm{Stack}}(p,L^*),P is regular,r>0,S⊆Cmax​rev(S)=p>0, L∗max​revStack​(p,L∗),

and uniform pricing corresponds to revenue-ordered assortments: for each cic_ici​ some revenue-ordered assortment earns the revenue of the uniform price cic_ici​, and every revenue-ordered assortment earns the revenue of some uniform price cic_ici​.

Milestones

  1. Lemma 4.14 (p. 23): for F⊆F′F\subseteq F'F⊆F′, ∣greedyM(F′,L)∣≥∣greedyM(F,L)∣|\mathrm{greedy}_M(F',L)|\ge|\mathrm{greedy}_M(F,L)|∣greedyM​(F′,L)∣≥∣greedyM​(F,L)∣ and F∩greedyM(F′,L)⊆greedyM(F,L)F\cap\mathrm{greedy}_M(F',L)\subseteq\mathrm{greedy}_M(F,L)F∩greedyM​(F′,L)⊆greedyM​(F,L).
  2. Property (14) (p. 31): orderings agreeing on a block partition make greedy pick the same number of elements per block.
  3. Lemma 4.15 (p. 23): the Stackelberg revenue does not depend on the compatible ordering.
  4. M′M'M′ is a matroid (p. 32).
  5. Regularity of P\mathcal PP (pp. 32–33).
  6. rev(S)\mathrm{rev}(S)rev(S) equals the Stackelberg revenue of pS(e)=min⁡{q:(e,q)∈S}p_S(e)=\min\{q:(e,q)\in S\}pS​(e)=min{q:(e,q)∈S} (pp. 33–34).
  7. Rounding prices up to the cost levels does not decrease revenue (p. 34).
  8. For prices in {c1,…,ck}∪{+∞}\{c_1,\dots,c_k\}\cup\{+\infty\}{c1​,…,ck​}∪{+∞}, rev(Sp)\mathrm{rev}(S_p)rev(Sp​) equals the Stackelberg revenue of ppp (pp. 34–35).

Significance

Theorem 4.16 transfers every guarantee for revenue-ordered assortments under regular models to uniform pricing in Stackelberg Matroid pricing. Specialising Theorems 3.1, 3.2 and 3.3 of the paper through it gives exactly the three bounds on uniform pricing proved by Cardinal et al. for Stackelberg Minimum Spanning Tree (Theorem 3 of their paper), which those authors showed to be tight (p. 24). It also places uniform pricing inside a general picture: it is a threshold policy for a regular choice model whose purchase probabilities come from a matroid greedy algorithm, and the paper notes that the argument extends to polymatroids.

The result is proved in the paper, with two steps left to the reader (M′M'M′ is a matroid) or justified briefly (rounding prices). No machine-checked proof exists. A formalization also produces a reusable greedy algorithm on finite matroids with its monotonicity properties (Lemma 4.14, (14)), which Mathlib at the pinned revision does not have.

Difficulty

The obvious argument identifies the customer's run of greedy on (M,R∪B)(M,R\cup B)(M,R∪B) with the run of greedy on (M′,R∪S)(M',R\cup S)(M′,R∪S) step by step. That identification fails as stated: the choice probabilities are defined with one fixed ordering LLL of R∪CR\cup\mathcal CR∪C, while the customer's ordering L∗L^*L∗ is any ordering compatible with the prices, and ties among blue elements of equal price, or between a product and a red element of equal cost, may be broken differently. Equality of revenues therefore needs the tie-independence property (14) applied to M′M'M′ as well as to MMM, which in turn needs M′M'M′ to be a matroid. A second gap is axiom (iv) for the no-purchase option: it requires that offering more products never decreases the number of products greedy selects, which is Lemma 4.14 (i)–(ii) for M′M'M′ combined, not a pointwise statement. Finally, the paper's "+∞+\infty+∞" price must be handled with the red base: without a base inside RRR, a blue element priced above all costs can still be bought and the optimum is unbounded.

Formalization scope

  • Elements and orderings. The matroid is a Mathlib Matroid α; RRR, BBB are Finset α with R∪BR\cup BR∪B equal to the ground set; costs and prices are functions α → ℝ, used only on RRR and BBB. A linear ordering is a duplicate-free List covering the set; greedy on FFF scans the whole list and skips elements outside FFF, so FFF is never re-sorted.
  • Greedy is defined for an arbitrary independence predicate on finite sets, so that it applies to M′M'M′ before M′M'M′ is shown to be a matroid.
  • Customer model. Compatibility (15a)–(15b), including blue priority on ties, is part of the definition of an admissible customer ordering. The red base is a field of every instance.
  • Assortment instance. Elements of R∪CR\cup\mathcal CR∪C live in α ⊕ (α × ℝ). The paper's conditions (1)–(2) for M′M'M′ contain two typos (a missing "∈X\in X∈X" and X\mathcal XX for XXX); the intended conditions are formalized. The price +∞+\infty+∞ is the real number 1+∑f∈Rc(f)1+\sum_{f\in R}c(f)1+∑f∈R​c(f), strictly above every red cost.
  • Goal shape. The theorem is stated for the explicit instance of the proof, for every admissible LLL. The optimum of the Stackelberg side is a maximum (IsGreatest) over positive real prices and all compatible customer orderings. Cost levels are quantified as elements of the set of red costs rather than by index.
  • Ruled out. An existential statement ("some regular instance has the same optimum") is met by a single product with revenue OPTStack\mathrm{OPT}_{\mathrm{Stack}}OPTStack​ and is not this theorem; so is a customer model without tie-breaking, a formalization without the red base, or Lemma 4.14 only for F=F′F=F'F=F′.

Contributions welcome: proofs of Lemma 4.14 and (14) for Mathlib matroids (reusable beyond this mission), the matroid property of parallel extensions such as M′M'M′, and the three revenue identities.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019 (Algorithmica, 2020). https://arxiv.org/abs/1606.01371
  • J. Cardinal, E. D. Demaine, S. Fiorini, G. Joret, S. Langerman, I. Newman and O. Weimann, The Stackelberg Minimum Spanning Tree Game, Algorithmica 59(2):129–144, 2011. https://doi.org/10.1007/s00453-009-9299-y
  • A. Schrijver, Combinatorial Optimization: Polyhedra and Efficiency, Vol. B, Algorithms and Combinatorics 24, Springer, 2003 (matroids and the greedy algorithm). https://link.springer.com/book/9783540443896
14 thms1 active userReviewed
Linear OptimizationOperations ResearchTheoretical Computer Science·Captain: mikedeng1

Online Primal-Dual Algorithms for Covering and Packing 3: A Deterministic O(log d log(n/OPT))-Competitive Algorithm for Online Unweighted Set CoverResearch Paper

Motivation

Online set cover is the basic covering problem in which the requests arrive over time. A ground set of elements and a family of sets are known in advance, but which elements must be covered is revealed one element at a time, and each arriving element has to be covered at once by a set chosen irrevocably. The problem models resource placement under unknown demand (facilities, servers, sensors that must serve clients as they appear) and is the prototype for a family of online covering problems.

Alon, Awerbuch, Azar, Buchbinder and Naor (SIAM J. Comput. 2009) gave the first deterministic algorithm, with competitive ratio O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) for nnn elements and mmm sets, and showed that no deterministic algorithm does better than Ω(log⁡mlog⁡n/(log⁡log⁡m+log⁡log⁡n))\Omega(\log m\log n/(\log\log m+\log\log n))Ω(logmlogn/(loglogm+loglogn)) on some instances. Buchbinder and Naor (Math. Oper. Res. 2009) recast the fractional part of that algorithm as an instance of a general online primal-dual scheme for covering and packing linear programs, and turned an offline pessimistic estimator of Srinivasan into an online potential function. The result, in their Section 5.1, is a deterministic algorithm whose ratio O(log⁡dlog⁡(n/OPT))O(\log d\log(n/OPT))O(logdlog(n/OPT)) depends on the maximum element frequency ddd instead of the number of sets mmm, and on the ratio n/OPTn/OPTn/OPT instead of nnn.

Timeline:

  • 2003 (conference), 2009 (journal): Alon et al., deterministic O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) for unweighted online set cover, with a potential ∑j∉Cn2wj\sum_{j\notin C} n^{2w_j}∑j∈/C​n2wj​.
  • 2005 (conference), 2009 (journal): Buchbinder and Naor, the general online fractional covering/packing scheme, and the derandomized rounding of this mission.

Setting

A set-cover instance consists of a finite ground set XXX of nnn elements and a finite family S\mathcal SS of mmm sets. For an element eee, Se\mathcal S_eSe​ is the collection of sets containing eee, and ddd bounds its size: ∣Se∣≤d|\mathcal S_e|\le d∣Se​∣≤d for every eee (the frequency). In the unweighted problem every set costs 111.

Elements arrive in a list σ\sigmaσ. The algorithm maintains:

  • fractional weights w(s)≥0w(s)\ge 0w(s)≥0 for the sets, produced by the paper's Section 3 scheme with {0,1}\{0,1\}{0,1} coefficients: when an element eee arrives that is not yet fractionally covered (∑s∈Sew(s)<1\sum_{s\in\mathcal S_e}w(s)<1∑s∈Se​​w(s)<1), its dual variable y(e)y(e)y(e) is raised to the least value at which
w(s)=max⁡{w(s), 1d(exp⁡(B2c(s)∑k: s∋eky(ek))−1)}(s∋e)w(s)=\max\Big\{w(s),\ \tfrac1d\Big(\exp\Big(\tfrac{B}{2c(s)}\textstyle\sum_{k:\ s\ni e_k}y(e_k)\Big)-1\Big)\Big\}\qquad(s\ni e)w(s)=max{w(s), d1​(exp(2c(s)B​∑k: s∋ek​​y(ek​))−1)}(s∋e)

gives ∑s∈Sew(s)≥1\sum_{s\in\mathcal S_e}w(s)\ge 1∑s∈Se​​w(s)≥1; here B>0B>0B>0 is a parameter;

  • a cover C⊆S\mathcal C\subseteq\mathcal SC⊆S that only grows; CCC is the set of elements covered by C\mathcal CC.

With f(e)=min⁡{1,exp⁡(−α+α∑s∋ew(s))}f(e)=\min\{1,\exp(-\alpha+\alpha\sum_{s\ni e}w(s))\}f(e)=min{1,exp(−α+α∑s∋e​w(s))}, the potential is Φ=Φ1+Φ2\Phi=\Phi_1+\Phi_2Φ=Φ1​+Φ2​ with

Φ1=1−∏e∈X∖C(1−f(e)),Φ2=exp⁡(∑s∈S((ln⁡2)χC(s)−αw(s))−OPT),\Phi_1=1-\prod_{e\in X\setminus C}\big(1-f(e)\big),\qquad \Phi_2=\exp\Big(\sum_{s\in\mathcal S}\big((\ln 2)\chi_{\mathcal C}(s)-\alpha w(s)\big)-OPT\Big),Φ1​=1−e∈X∖C∏​(1−f(e)),Φ2​=exp(s∈S∑​((ln2)χC​(s)−αw(s))−OPT),

where OPTOPTOPT is the optimum number of sets covering the arrived elements, assumed known, r=eln⁡(e/(e−1))r=e\ln(e/(e-1))r=eln(e/(e−1)) and α=max⁡{1,ln⁡(rn/OPT)}\alpha=\max\{1,\ln(rn/OPT)\}α=max{1,ln(rn/OPT)}. The rounding rule: each time the weight of a set sss is augmented, sss is added to C\mathcal CC if this does not increase Φ\PhiΦ.

Formalization targets

Goal: Lemma 5.2

Every arriving element is covered by C\mathcal CC, and at every time

∣C∣ ≤ (2αln⁡(1+d)+1) OPTln⁡2,α=max⁡{1,ln⁡rnOPT}.|\mathcal C|\ \le\ \frac{\big(2\alpha\ln(1+d)+1\big)\,OPT}{\ln 2},\qquad \alpha=\max\Big\{1,\ln\frac{rn}{OPT}\Big\}.∣C∣ ≤ ln2(2αln(1+d)+1)OPT​,α=max{1,lnOPTrn​}.

This is the paper's OPT⋅O(log⁡dlog⁡(n/OPT))OPT\cdot O(\log d\log(n/OPT))OPT⋅O(logdlog(n/OPT)) with the constant its proof gives.

Milestones

  1. Theorem 3.2 (covering half): for every B>0B>0B>0, the {0,1}\{0,1\}{0,1} scheme with frequency bound ℓ\ellℓ yields a fractional cover of cost at most 2ln⁡(1+ℓ)2\ln(1+\ell)2ln(1+ℓ) times that of any fractional cover.
  2. Lemma 5.1 (i): initially Φ≤1\Phi\le 1Φ≤1; Φ>0\Phi>0Φ>0 in every state.
  3. Lemma 5.1 (ii): after the weight of a set is augmented by δ≥0\delta\ge0δ≥0, taking the set or excluding it leaves Φ\PhiΦ no larger than before.

Significance

The bound improves Alon et al.'s O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) whenever sets are many but each element lies in few of them (d≪md\ll md≪m), and whenever the optimum is large compared with nnn. It also shows that the fractional part and the rounding part of an online covering algorithm can be designed separately: any online fractional solution with a competitive guarantee can be rounded deterministically by an online potential function. The same method gives the routing result of Section 5.2 of the paper.

As far as is known, none of these results is machine-checked. Related formal work on the platform covers Alon et al.'s algorithm and its O(log⁡mlog⁡n)O(\log m\log n)O(logmlogn) bound (a different potential and a different fractional update), and the general covering scheme of Buchbinder and Naor's monograph; neither states Lemma 5.1 or Lemma 5.2, nor Theorem 3.2 for the {0,1}\{0,1\}{0,1} scheme with ℓ\ellℓ in place of nnn. A complete development here would give a verified derandomized rounding argument, reusable for other online covering problems.

Difficulty

The difficulty is not in the final inequality, which follows from Φ2≤1\Phi_2\le1Φ2​≤1 in one line, but in keeping Φ≤1\Phi\le1Φ≤1 throughout. The decision to take a set must be made with no knowledge of future elements, and the potential has to account for elements that may never arrive: Φ1\Phi_1Φ1​ ranges over the whole ground set. Lemma 5.1 (ii) asks that, for every current state, one of the two decisions does not increase a non-linear function of all uncovered elements at once and this must hold for every state, not only for the states a particular run reaches. On the fractional side, Theorem 3.2 is only asserted, "along the same lines" as Theorem 3.1, so its constant has to be re-derived with ℓ\ellℓ in place of nnn and with the scheme's continuous increase made discrete.

A tempting shortcut is to bound ∣C∣|\mathcal C|∣C∣ by the number of rounds or by ∑sw(s)\sum_s w(s)∑s​w(s) directly; neither gives a logarithmic factor in n/OPTn/OPTn/OPT, which comes only from the choice of α\alphaα in Φ\PhiΦ.

Formalization scope

The instance is the published OnlinePrimalDual.OnlineSetCover.SetCoverInstance (elements E, set indices T, incidence elemSets, positive costs c), with elementWeight and coveredBy. Unit costs are the hypothesis ∀ s, inst.c s = 1 in Lemma 5.2; Theorem 3.2 is stated for general positive costs. Logarithms are natural (Real.log), because they invert Real.exp.

Committed conventions:

  • The fractional scheme is the discrete form of the continuous increase: in each round y(e)y(e)y(e) is the least t≥0t\ge0t≥0 (an sInf) at which the new constraint holds. Arrival lists may repeat elements; an element already covered changes nothing.
  • The rounding treats each set's increase within a round as one augmentation; the augmented sets are processed one at a time in the order of a list ord containing every set, and a set is added when Φ(w+δs1s,C∪{s})≤Φ(w,C)\Phi(w+\delta_s\mathbf 1_s,\mathcal C\cup\{s\})\le\Phi(w,\mathcal C)Φ(w+δs​1s​,C∪{s})≤Φ(w,C). Lemma 5.2 holds for every order.
  • The algorithm is a function, so every run exists.
  • OPTOPTOPT is a natural number ≥1\ge1≥1 bounding the size of some cover of the arrived elements; d≥1d\ge1d≥1 bounds the frequency of every element. At the true optimum and the maximum frequency this is the paper's statement.
  • Explicit constants replacing O(⋅)O(\cdot)O(⋅): Theorem 3.2's O(log⁡ℓ)O(\log\ell)O(logℓ) is 2ln⁡(1+ℓ)2\ln(1+\ell)2ln(1+ℓ); Lemma 5.2's OPT⋅O(log⁡dlog⁡(n/OPT))OPT\cdot O(\log d\log(n/OPT))OPT⋅O(logdlog(n/OPT)) is (2αln⁡(1+d)+1) OPT/ln⁡2(2\alpha\ln(1+d)+1)\,OPT/\ln2(2αln(1+d)+1)OPT/ln2 with α=max⁡{1,ln⁡(rn/OPT)}\alpha=\max\{1,\ln(rn/OPT)\}α=max{1,ln(rn/OPT)}.
  • The proof of Lemma 5.1 (i) on p. 13 prints ddd where rrr is meant ("exp⁡(−dne−α)\exp(-dne^{-\alpha})exp(−dne−α)", "α≥ln⁡(dn/OPT)\alpha\ge\ln(dn/OPT)α≥ln(dn/OPT)"); the statement uses rrr, and rrr is printed "eln⁡(e/e−1)e\ln(e/e-1)eln(e/e−1)" for eln⁡(e/(e−1))e\ln(e/(e-1))eln(e/(e−1)).

Not in scope: the packing half of Theorem 3.2, Theorem 3.1 for general coefficients, and the doubling wrapper of p. 12 that removes the assumption that OPTOPTOPT is known. The guarantee is for the algorithm run with the stated α\alphaα; a formalization in which Φ\PhiΦ's product ranges only over arrived elements, or in which the weights or the chosen family are free variables constrained by hypotheses instead of being produced by the algorithm (OPTOPTOPT is an input of the algorithm, as the page assumes it known), would be a different and weaker statement and is ruled out.

Contributions welcome: proofs of the three milestones and of the goal, and general lemmas about the fractional round (attainment of the least ttt, monotonicity of the weights) that other online covering missions can reuse.

Selected references

  • N. Buchbinder, J. Naor, Online Primal-Dual Algorithms for Covering and Packing, Mathematics of Operations Research 34(2), 2009. https://doi.org/10.1287/moor.1080.0363
  • N. Alon, B. Awerbuch, Y. Azar, N. Buchbinder, J. Naor, The Online Set Cover Problem, SIAM Journal on Computing 39(2), 2009. https://doi.org/10.1137/060661946
  • A. Srinivasan, Improved approximation guarantees for packing and covering integer programs, SIAM Journal on Computing 29(2), 1999. https://doi.org/10.1137/S0097539796314240
  • N. Buchbinder, J. Naor, The Design of Competitive Online Algorithms via a Primal-Dual Approach, Foundations and Trends in Theoretical Computer Science 3(2–3), 2009. https://doi.org/10.1561/0400000024
10 thms1 active userReviewed
Discrete GeometryGraph Theory·Captain: mikedeng1

Diameter of Polyhedra: Limits of Abstraction 1: The Base Abstraction in Dimension n/4 with n Facets Has Diameter Ω(n²/log n)Research Paper

Motivation

Whether the diameter of a polyhedron is bounded by a polynomial in its number of facets is one of the central open questions of convex geometry, and it bears directly on linear programming: the diameter of the graph of a polyhedron is a lower bound on the number of pivots of any simplex method that moves along edges. Write △u(d,n)\triangle_u(d,n)△u​(d,n) for the largest diameter of a ddd-dimensional polyhedron with nnn facets. The best known bounds are

n−d+⌊d/5⌋  ≤  △u(d,n)  ≤  n1+log⁡d,n - d + \lfloor d/5 \rfloor \;\le\; \triangle_u(d,n) \;\le\; n^{1+\log d},n−d+⌊d/5⌋≤△u​(d,n)≤n1+logd,

due to Klee and Walkup (1967) and Kalai and Kleitman (1992), and the gap between them has stood for decades.

A recurring approach is to isolate a few purely combinatorial properties of polyhedra and to prove bounds for every object having them. Eisenbrand, Hähnle, Razborov and Rothvoß (2010) proposed a base abstraction defined by a single property, shared by all abstractions from which earlier bounds were derived. They showed that the known upper bounds survive in it, and, as their main result, that it admits almost-quadratic diameter. This mission formalizes that lower bound.

Timeline.

  • 1967: Klee and Walkup, lower bound n−d+⌊d/5⌋n-d+\lfloor d/5\rfloorn−d+⌊d/5⌋ for unbounded polyhedra.
  • 1970: Larman, △u(d,n)≤2d−1n\triangle_u(d,n) \le 2^{d-1} n△u​(d,n)≤2d−1n (linear in fixed dimension).
  • 1974: Adler, Dantzig and Murty study abstract polytopes and their diameters.
  • 1992: Kalai and Kleitman, the quasi-polynomial bound n1+log⁡dn^{1+\log d}n1+logd.
  • 2010: Eisenbrand, Hähnle, Razborov and Rothvoß, the base abstraction; both upper bounds hold in it, and its diameter is Ω(n2/log⁡n)\Omega(n^2/\log n)Ω(n2/logn) in dimension n/4n/4n/4.
  • 2012: Santos disproves the Hirsch conjecture for polytopes; the polynomial question remains open.

Setting

Let [n]={1,…,n}[n] = \{1,\dots,n\}[n]={1,…,n} and let ([n]d)\binom{[n]}{d}(d[n]​) be the family of its ddd-element subsets. Consider a graph GGG whose vertex set VVV is a nonempty subfamily of ([n]d)\binom{[n]}{d}(d[n]​). The graph satisfies condition i) if

for each u,v∈Vu, v \in Vu,v∈V there exists a path connecting uuu and vvv whose intermediate vertices all contain u∩vu \cap vu∩v.

The base abstraction Bd,n\mathcal B_{d,n}Bd,n​ is the class of all such graphs; ddd is called the dimension and nnn the number of facets. Write distG(u,v)\mathrm{dist}_G(u,v)distG​(u,v) for the length of a shortest uuu–vvv path, and

D(d,n)=max⁡{distG(u,v):G∈Bd,n, u,v∈V(G)}.D(d,n) = \max\{\mathrm{dist}_G(u,v) : G \in \mathcal B_{d,n},\ u, v \in V(G)\}.D(d,n)=max{distG​(u,v):G∈Bd,n​, u,v∈V(G)}.

The geometry behind the definition: in a non-degenerate polyhedron every vertex lies on exactly ddd facets, so labelling vertices by their facets gives distinct ddd-subsets of [n][n][n]; and two vertices can be joined without leaving the smallest face containing both, i.e. through vertices lying on every facet common to uuu and vvv. Hence △u(d,n)≤D(d,n)\triangle_u(d,n) \le D(d,n)△u​(d,n)≤D(d,n) (the subject of a separate mission of this series).

In the Lean development, EHRR10.LowerBound.CondI V G is condition i), and EHRR10.LowerBound.InB d n V G says that VVV is nonempty, all its members have ddd elements, and condition i) holds.

Formalization targets

Goal: almost-quadratic diameter in dimension n/4n/4n/4

D(n/4, n)=Ω ⁣(n2log⁡n).D(n/4,\, n) = \Omega\!\left(\frac{n^2}{\log n}\right).D(n/4,n)=Ω(lognn2​).

Explicitly, there are c>0c > 0c>0 and NNN such that for every n≥Nn \ge Nn≥N with 4∣n4 \mid n4∣n there is a graph G∈Bn/4, nG \in \mathcal B_{n/4,\,n}G∈Bn/4,n​ and two vertices u,vu, vu,v of GGG with distG(u,v)≥c n2/log⁡n\mathrm{dist}_G(u,v) \ge c\, n^2/\log ndistG​(u,v)≥cn2/logn. The constants ccc and NNN are left unfixed, so the goal survives any improvement in the constant.

There are no milestones: the source is an extended abstract that states the result without a proof and names no intermediate lemma.

Significance

The result. Since △u(d,n)≤D(d,n)\triangle_u(d,n) \le D(d,n)△u​(d,n)≤D(d,n), any proof of a bound on DDD bounds the diameter of all non-degenerate polyhedra. The paper shows that the Kalai–Kleitman and Larman bounds are of this kind. The lower bound shows the limits of the approach: an argument that uses only condition i) cannot prove a linear (Hirsch-type) bound, because D(n/4,n)D(n/4,n)D(n/4,n) grows almost quadratically in nnn. Any proof of a linear bound must use further geometric properties of polyhedra. A polynomial, even quadratic, upper bound for DDD is not excluded by the result.

Formalizing it. The result is proved in the literature (Math. Oper. Res. 35 (2010) 786–794); it has no machine-checked proof. A formal proof requires an explicit combinatorial construction of graphs in the base abstraction with large diameter and a probabilistic existence argument, neither of which is currently on the platform. The platform's Hirsch campaign has the same abstraction in its layer-family form (Hirsch_clf, whose description mentions the Ω(n2/log⁡n)\Omega(n^2/\log n)Ω(n2/logn) families in prose) and the proved upper bound Hirsch.clf_saturated_quadratic_bound, L≤d(n−d)L \le d(n-d)L≤d(n−d), on a subclass that contains the constructions of this paper; neither states the lower bound.

Difficulty

A long path through all of ([n]d)\binom{[n]}{d}(d[n]​) has exponential length but violates condition i): two vertices of the path generally share elements that some intermediate vertex lacks. Condition i) must hold for every pair of vertices simultaneously, which forces the vertices containing any fixed set to induce a connected subgraph. Building a long graph means arranging many ddd-sets so that each of these exponentially many connectivity constraints holds at once while the graph distance between two vertices stays large. The authors announce that their construction relies on disjoint covering designs, and that the existence of such designs with the required parameters is obtained from the Lovász Local Lemma.

Formalization scope

  • [n][n][n] is Fin n; vertices are the ddd-sets themselves, V : Finset (Finset (Fin n)) and G : SimpleGraph V. Distinct vertices carry distinct labels.
  • Condition i) is stated with walks, and requires every vertex of the walk (endpoints included, which is automatic) to contain u∩vu \cap vu∩v. A walk contains a path on a subset of its vertices, so this is the page's condition.
  • "Connected" is encoded as V≠∅V \ne \emptysetV=∅; connectivity follows from condition i).
  • The page writes D(n/4,n)D(n/4,n)D(n/4,n) without rounding; the goal is read literally: 4∣n4 \mid n4∣n and dimension exactly n/4n/4n/4.
  • Ω(⋅)\Omega(\cdot)Ω(⋅) is made explicit as ∃c>0,∃N,∀n≥N\exists c > 0, \exists N, \forall n \ge N∃c>0,∃N,∀n≥N. The logarithm is natural (Real.log); the base only rescales ccc.
  • "D(n/4,n)≥BD(n/4,n) \ge BD(n/4,n)≥B" is stated as the existence of one graph in the class with two vertices at SimpleGraph.dist at least BBB; DDD is not defined as a supremum, so no junk value enters. SimpleGraph.dist returns 000 for unreachable pairs, which could only make the statement harder.
  • Only the lower half of Ω\OmegaΩ is claimed, as on the page.

The goal cannot be trivialized by a non-injective labelling (vertices are the sets), by checking condition i) only on adjacent pairs (it is required for all pairs, with u∩vu\cap vu∩v kept along the whole walk), or by a zero constant (c>0c > 0c>0 is required).

Infrastructure that a solution needs and that is reusable beyond this mission: covering designs and disjoint covering designs over Finset (Fin n), the Lovász Local Lemma (posed on the platform as AppliedComb.ManyFaces.local_lemma_symmetric), and a toolkit for proving condition i) for graphs assembled from blocks. Contributions of any of these as separate theorems are welcome.

Selected references

  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Dagstuhl Seminar Proceedings 10211 (Flexible Network Design), 2010 (the extended abstract this mission cites). http://drops.dagstuhl.de/opus/volltexte/2010/2724
  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoß, Diameter of Polyhedra: Limits of Abstraction, Mathematics of Operations Research 35(4), 786–794, 2010 (full proofs). https://doi.org/10.1287/moor.1100.0470
  • G. Kalai, D. J. Kleitman, A quasi-polynomial bound for the diameter of graphs of polyhedra, Bull. Amer. Math. Soc. 26(2), 315–316, 1992. https://doi.org/10.1090/S0273-0979-1992-00285-9
  • D. G. Larman, Paths on polytopes, Proc. London Math. Soc. 20(3), 161–178, 1970. https://doi.org/10.1112/plms/s3-20.1.161
  • V. Klee, D. W. Walkup, The d-step conjecture for polyhedra of dimension d < 6, Acta Math. 117, 53–78, 1967. https://doi.org/10.1007/BF02395040
  • F. Santos, A counterexample to the Hirsch conjecture, Annals of Mathematics 176(1), 383–412, 2012. https://doi.org/10.4007/annals.2012.176.1.7
  • N. Alon, J. Spencer, The Probabilistic Method, 3rd ed., Wiley, 2008 (Lovász Local Lemma).
2 thms1 active userReviewed
Group TheoryMachine Learning·Captain: mikedeng1

A Characterization of Multiclass Learnability 2: A Concept Class with Natarajan Dimension 1 and Infinite DS DimensionResearch Paper

Motivation

In binary classification the VC dimension decides PAC learnability: a class of {0,1}\{0,1\}{0,1}-valued functions is learnable from finitely many examples exactly when its VC dimension is finite. Multiclass classification, where a predictor outputs one of many labels, arises whenever the label set is large: language models choosing a next token, image recognition over open vocabularies, structured prediction. For finitely many labels the Natarajan dimension plays the role of the VC dimension (Natarajan 1989; Ben-David, Cesa-Bianchi, Haussler and Long 1995). Whether it still characterizes learnability when the label set is infinite stayed open for three decades.

Brukhim, Carmon, Dinur, Moran and Yehudayoff (arXiv:2203.01550, FOCS 2022) settled both directions. Their Theorem A shows that the DS dimension of Daniely and Shalev-Shwartz (COLT 2014, PMLR 35) characterizes multiclass PAC learnability for every label set. Their Theorem 2, the goal of this mission, shows that the Natarajan dimension does not: there is a class whose Natarajan dimension is 111 and whose DS dimension is infinite.

Timeline:

  • 1989: Natarajan introduces his dimension and proves it gives sample-complexity bounds when the label set is finite.
  • 1995: Ben-David, Cesa-Bianchi, Haussler and Long show that, for finite label sets, every "reasonable" extension of the VC dimension characterizes learnability.
  • 2003: Januszkiewicz and Świątkowski construct, for every dimension, finite simplicial complexes without empty squares from coset complexes of finite groups (Comment. Math. Helv. 78(3), 555–583); the multiclass paper uses this construction for its separation.
  • 2014: Daniely and Shalev-Shwartz introduce the DS dimension, prove that finite DS dimension is necessary for learnability, and ask whether it is sufficient.
  • 2022: Brukhim et al. prove that finite DS dimension is sufficient and that the Natarajan dimension fails to characterize learnability for infinite label sets.

Setting

A concept class is a set H⊆YX\mathcal H \subseteq \mathcal Y^{\mathcal X}H⊆YX of functions from a domain X\mathcal XX to a label set Y\mathcal YY, with no finiteness assumption on either. For a sequence S=(x1,…,xn)∈XnS = (x_1, \dots, x_n) \in \mathcal X^nS=(x1​,…,xn​)∈Xn, the projection H∣S⊆Yn\mathcal H|_S \subseteq \mathcal Y^nH∣S​⊆Yn is the set of words (h(x1),…,h(xn))(h(x_1), \dots, h(x_n))(h(x1​),…,h(xn​)), h∈Hh \in \mathcal Hh∈H.

  • SSS is N-shattered if there are f,g:[n]→Yf, g : [n] \to \mathcal Yf,g:[n]→Y with f(i)≠g(i)f(i) \ne g(i)f(i)=g(i) for every iii and H∣S⊇{f(1),g(1)}×⋯×{f(n),g(n)}\mathcal H|_S \supseteq \{f(1), g(1)\} \times \dots \times \{f(n), g(n)\}H∣S​⊇{f(1),g(1)}×⋯×{f(n),g(n)}: the projection contains a copy of the Boolean cube. The Natarajan dimension dN(H)d_N(\mathcal H)dN​(H) is the largest nnn for which some S∈XnS \in \mathcal X^nS∈Xn is N-shattered, or ∞\infty∞.
  • A pseudo-cube of dimension ddd is a non-empty, finite B⊆YdB \subseteq \mathcal Y^dB⊆Yd in which every word hhh has, for every coordinate iii, an iii-neighbour: a word g∈Bg \in Bg∈B with g(i)≠h(i)g(i) \ne h(i)g(i)=h(i) and g(j)=h(j)g(j) = h(j)g(j)=h(j) for j≠ij \ne ij=i. SSS is DS-shattered if H∣S\mathcal H|_SH∣S​ contains an nnn-dimensional pseudo-cube, and the DS dimension dDS(H)d_{DS}(\mathcal H)dDS​(H) is the largest such nnn, or ∞\infty∞.

Every Boolean cube is a pseudo-cube, so dN≤dDSd_N \le d_{DS}dN​≤dDS​. The hexagon {12,32,34,54,56,16}⊆{1,…,6}2\{12, 32, 34, 54, 56, 16\} \subseteq \{1,\dots,6\}^2{12,32,34,54,56,16}⊆{1,…,6}2 is a 2-dimensional pseudo-cube that contains no Boolean square.

The milestones pass through simplicial complexes: downward-closed families of finite sets. A complex is good if it is finite, pure, has a proper coloring rrr of its vertices with dim⁡(C)+1\dim(C)+1dim(C)+1 colors, and satisfies replacement (every vertex of every face can be exchanged for a new vertex). A good complex CCC with coloring rrr defines the class B(C,r)B(C, r)B(C,r) of its top faces, each written as the word listing its vertices by color. A square is a 4-cycle of distinct vertices in the 1-skeleton; it is empty if neither diagonal is an edge. The coset complex CF(H1,…,Hd)C_F(H_1, \dots, H_d)CF​(H1​,…,Hd​) of subgroups of a group FFF has the cosets gHigH_igHi​ as vertices and the sets of cosets with a common point as faces.

Formalization targets

Goal: Theorem 2 (p. 4)

∃ X,Y, H⊆YX:dN(H)=1anddDS(H)=∞.\exists\, \mathcal X, \mathcal Y,\ \mathcal H \subseteq \mathcal Y^{\mathcal X}:\qquad d_N(\mathcal H) = 1 \quad\text{and}\quad d_{DS}(\mathcal H) = \infty.∃X,Y, H⊆YX:dN​(H)=1anddDS​(H)=∞.

Milestones

  1. Theorem 45 (p. 30; Januszkiewicz–Świątkowski): for every d>1d > 1d>1 a finite group FFF and subgroups H1,…,HdH_1, \dots, H_dH1​,…,Hd​ with (⋂j≠iHj)∖Hi≠∅(\bigcap_{j\ne i} H_j) \setminus H_i \ne \emptyset(⋂j=i​Hj​)∖Hi​=∅ for all iii, whose coset complex has no empty squares.
  2. Proposition 46 (p. 31): such a coset complex has dimension d−1d - 1d−1, is good and has no empty squares.
  3. Proposition 42 (p. 28): a ddd-dimensional good complex with a proper coloring rrr yields the (d+1)(d+1)(d+1)-dimensional pseudo-cube B(C,r)B(C, r)B(C,r); conversely every pseudo-cube yields a good complex C(B)C(B)C(B).
  4. Proposition 43 (p. 29): dN(B(C,r))≥2d_N(B(C,r)) \ge 2dN​(B(C,r))≥2 iff CCC has a square v0v1v2v3v_0 v_1 v_2 v_3v0​v1​v2​v3​ with r(v0)=r(v2)r(v_0) = r(v_2)r(v0​)=r(v2​) and r(v1)=r(v3)r(v_1) = r(v_3)r(v1​)=r(v3​).
  5. Corollary 44 (p. 29): a good complex without empty squares gives dN(B(C,r))≤1d_N(B(C, r)) \le 1dN​(B(C,r))≤1 for every proper coloring.
  6. Proof of Theorem 2 (p. 32): for every d≥1d \ge 1d≥1, a ddd-dimensional pseudo-cube with Natarajan dimension exactly 111.

Significance

Theorem 2 shows that the classical generalization of the VC dimension to many labels is the wrong invariant once the label set is infinite: a class can contain no Boolean square at all and still be unlearnable, because it contains pseudo-cubes of every dimension. Combined with the necessity of finite DS dimension, it gives a class that is not PAC learnable although its Natarajan dimension is 111, and it identifies pseudo-cubes, not Boolean cubes, as the relevant combinatorial obstruction. It also links learning theory to a problem studied in geometric group theory, finite "flag-no-square" complexes.

The paper's proof is complete modulo Theorem 45, which it imports from Januszkiewicz–Świątkowski 2003. None of these results is formalized. A formalization would give machine-checked versions of the dictionary between concept classes and properly colored complexes (Propositions 42–44), of the coset-complex translation (Proposition 46), and of the final disjoint-union argument; Theorem 45 itself, which rests on Coxeter-group and topological arguments, is a separate and substantial formalization target.

Difficulty

Infinite complexes that are pure, properly colored, satisfy replacement and have no empty squares are easy to build: grow a tree of faces indefinitely. The definition of a pseudo-cube demands finiteness, and the difficulty is entirely there: one must "fold" such an infinite object into a finite one without creating an empty square. The obvious finite candidate, the group (Z/2)d(\mathbb Z/2)^d(Z/2)d with its coordinate subgroups, produces the Boolean cube, whose complex is full of empty squares. Theorem 45 is the input that resolves this, and it is far beyond the rest of the argument.

Formalization scope

All declarations live in the namespace MulticlassDS.NatGap.

  • Concept classes are Set (X → Y) with arbitrary types; [n][n][n] is Fin n (0-based), and shattering is defined for sequences Fin n → X, as in the paper.
  • Both dimensions are ℕ∞-valued suprema, so "infinite DS dimension" is dsDim H = ⊤. An ℕ-valued supremum would silently return 000 on an unbounded family and would trivialize the goal.
  • The goal requires the Natarajan dimension to be exactly 111; an upper bound alone holds for any class with at most one element.
  • Pseudo-cubes are required to be finite (Definition 5). Without finiteness, the tree classes of Example 8 would already have infinite "DS dimension".
  • Complexes are Set (Finset V). The dimension is the predicate HasDim C d, not a natural-number subtraction, and colors are Fin (d + 1).
  • Replacement is stated with a new vertex u∉fu \notin fu∈/f. The page writes "u≠vu \ne vu=v", but read literally that allows u∈fu \in fu∈f, which makes the condition hold by downward closure and makes Proposition 42 false; the proofs of Propositions 42 and 46 use a new vertex.
  • Coset-complex vertices are left cosets as subsets of the group, not pairs (index, coset).
  • Proposition 46 states dimension d−1d - 1d−1 under d>1d > 1d>1, where the subtraction is exact; the converse of Proposition 42 is indexed by d+1d + 1d+1 and ddd to avoid it.
  • The proof-of-Theorem-2 milestone says "for every ddd"; it is posed for d≥1d \ge 1d≥1, because at d=0d = 0d=0 the only pseudo-cube has Natarajan dimension 000.

Welcome contributions: proofs of Propositions 42–44 and 46 and of the goal from the milestones, which need only finite combinatorics and elementary group theory; and, separately, a formalization of the Januszkiewicz–Świątkowski construction behind Theorem 45. The definitions of pseudo-cubes, the DS dimension and good complexes are reusable by the companion mission on sample compression and by any later work on multiclass learnability.

Selected references

  • N. Brukhim, D. Carmon, I. Dinur, S. Moran, A. Yehudayoff, A Characterization of Multiclass Learnability, arXiv:2203.01550v1, 2022 (FOCS 2022). https://arxiv.org/abs/2203.01550
  • T. Januszkiewicz, J. Świątkowski, Hyperbolic Coxeter groups of large dimension, Comment. Math. Helv. 78(3) (2003), 555–583 (reference [Januszkiewicz and Świątkowski 2003] of arXiv:2203.01550v1, p. 33).
  • A. Daniely, S. Shalev-Shwartz, Optimal learners for multiclass problems, COLT 2014, PMLR 35, 287–316. https://proceedings.mlr.press/v35/
  • B. K. Natarajan, On learning sets and functions, Machine Learning 4 (1989), 67–97. https://doi.org/10.1007/BF00114804
  • S. Ben-David, N. Cesa-Bianchi, D. Haussler, P. M. Long, Characterizations of learnability for classes of {0,…,n}-valued functions, J. Comput. Syst. Sci. 50(1) (1995), 74–86. https://doi.org/10.1006/jcss.1995.1008
10 thms1 active userReviewed
Graph TheoryOperations ResearchOptimization+1·Captain: mikedeng1

Maximal Flow Through a Network II: In an ab-Planar Network Some Chain from Source to Sink Meets Every Cut Exactly OnceResearch Paper

Motivation

The maximum flow problem asks how much of a commodity can be shipped from a source to a sink through a network whose arcs have limited capacities. L. R. Ford, Jr. and D. R. Fulkerson's 1956 paper Maximal Flow Through a Network proved the minimal cut theorem: the largest flow value equals the smallest total capacity of a set of arcs that separates source from sink. That theorem is formalized in the companion mission Maximal Flow Through a Network I.

The second section of the same paper treats a special class of networks, those that remain planar after an arc from source to sink is added. For these networks the paper shows that one particular source–sink chain crosses every minimal separating set exactly once. This structural fact turns the minimal cut theorem into a simple computing procedure: repeatedly push as much flow as possible along such a chain and delete the arcs it saturates. The paper notes that G. Dantzig had conjectured, before the minimal cut theorem was proved, that this procedure yields a maximal flow on planar networks. The same "uppermost path" idea underlies later algorithms for maximum flow in planar graphs with source and sink on a common face (Itai and Shiloach, 1979).

The statement is short and purely combinatorial in its conclusion, but its hypothesis is topological. This mission isolates that theorem.

Setting

A network NNN has a finite set VVV of vertices and a finite set EEE of arcs. Each arc eee joins two distinct end vertices, written tail(e)\mathrm{tail}(e)tail(e) and head(e)\mathrm{head}(e)head(e); arcs carry no direction, and two arcs may join the same pair of vertices. Two distinct vertices are distinguished, the source aaa and the sink bbb, and each arc carries a positive capacity (capacities play no role in the target below).

A chain joining uuu and www is a set CCC of distinct arcs that can be arranged as α1(v0v1),α2(v1v2),…,αk(vk−1vk)\alpha_1(v_0v_1), \alpha_2(v_1v_2), \dots, \alpha_k(v_{k-1}v_k)α1​(v0​v1​),α2​(v1​v2​),…,αk​(vk−1​vk​) with v0=uv_0 = uv0​=u, vk=wv_k = wvk​=w, and the vertices v0,…,vkv_0, \dots, v_kv0​,…,vk​ pairwise distinct; each arc may be traversed in either direction. The empty set is the null chain from uuu to uuu.

A set DDD of arcs is a disconnecting set if every chain joining aaa and bbb contains an arc of DDD. A disconnecting set none of whose proper subsets is disconnecting is a cut.

The network is ab-planar if the graph of NNN, together with one additional arc joining aaa and bbb, can be drawn in the plane without crossings: vertices go to distinct points of R2\mathbb R^2R2; each arc, including the added arc ababab, goes to an injective continuous path between the points of its end vertices; no arc passes through a vertex other than its ends; and two distinct arcs meet only at endpoints of both. In Lean the drawing is the structure ABPlaneDrawing N, and NNN is ab-planar when Nonempty (ABPlaneDrawing N). The section's standing assumption is that no arc of NNN already joins aaa and bbb.

Formalization targets

Goal: Theorem 2 (p. 403)

If NNN is ab-planar, no arc of NNN joins aaa and bbb, and some chain joins aaa and bbb, then

∃ T a chain joining a and b  such that  ∣T∩D∣=1  for every cut D of N.\exists\, T \text{ a chain joining } a \text{ and } b \ \text{ such that }\ |T \cap D| = 1 \ \text{ for every cut } D \text{ of } N.∃T a chain joining a and b  such that  ∣T∩D∣=1  for every cut D of N.

This is FordFulkerson56.Planar.ab_planar_exists_chain_meeting_each_cut_once. "Precisely once" is exact cardinality one, neither "at least once" (true of every chain) nor "at most once".

Milestone: a chain meeting a cut in one prescribed arc (proof of Theorem 2, p. 403)

For every network NNN, every cut DDD and every arc α∈D\alpha \in Dα∈D, there is a chain CCC joining aaa and bbb with C∩D={α}C \cap D = \{\alpha\}C∩D={α}. No planarity is involved; the statement is what the minimality of a cut provides to the proof.

Further item: the Fig. 2 example (p. 403)

In the "gas, water, electricity" graph K3,3K_{3,3}K3,3​ with the arc ababab removed, every chain joining aaa and bbb meets some cut in three arcs. This network is not ab-planar, so the example shows that the planarity hypothesis of Theorem 2 cannot be dropped.

Significance

Theorem 2 and the minimal cut theorem together give the paper's procedure for planar networks: if TTT meets every cut once, then imposing a flow kkk on TTT lowers the value of every cut by exactly kkk, so the minimal cut value, and hence the maximal flow value, drops by kkk. Saturated arcs can then be deleted and the step repeated. Without the "exactly once" property the reduction could overshoot the cut structure, and the greedy step would not be justified. The theorem is also one of the earliest instances of the link between planarity and cut structure that later underlies planar duality arguments for minimum cuts.

The result has been known since 1956 and is not open. No machine-checked version is recorded on the platform, and Mathlib, at the pinned revision, has neither planar graphs nor the Jordan curve theorem. A formal proof would be the first formalized statement about source–sink planar networks in this library, and the counterexample item records, as a checkable fact, that the hypothesis is necessary.

Difficulty

The conclusion is combinatorial while the hypothesis is a drawing in R2\mathbb R^2R2. The paper's proof normalises the drawing (the added arc ababab on the outer boundary, the graph in a vertical strip with aaa on the left line and bbb on the right), selects the "top-most" chain from aaa to bbb, and argues that a chain meeting a cut below the top-most chain must cross another such chain. Each of these steps rests on plane topology: the existence of the outer region, the meaning of "top-most", and the fact that two chains with interleaved endpoints on a boundary must intersect, which is a form of the Jordan curve theorem.

The naive purely combinatorial route fails: the analogous statement for arbitrary networks is false (Fig. 2), so any argument has to use the drawing somewhere. Replacing the drawing by a combinatorial embedding (rotation systems, faces) is possible but then requires proving that the two notions agree, which is again Jordan-curve territory.

Formalization scope

Conventions committed to in the Lean statements:

  • Vertices and arcs are finite types V, E with decidable equality. Arcs are undirected, may be parallel, and have two distinct end vertices. Source and sink are distinct, capacities are positive (structure Network).
  • A chain is a Finset E that is the arc set of some arrangement as a simple path (IsChainWalk, IsChain); the null chain is allowed.
  • IsDisconnecting and IsCut quantify over all chains joining source and sink; a cut is a disconnecting set no proper subset of which is disconnecting.
  • ab-planarity is a plane drawing of the graph with the extra arc indexed by none : Option E, with injective Paths in ℝ × ℝ as arcs.

Hypotheses of the goal: hno_ab, the standing assumption of §2 (no arc joins aaa and bbb, p. 403); hconn, that some chain joins aaa and bbb. The second is not stated in the paper; its proof starts from "the chain joining a and b which is top-most", which presupposes one, and without it the statement is false (if aaa and bbb are disconnected, the empty set is a cut and no chain exists).

The drawing structure is satisfiable (a three-vertex path network has an explicit drawing), so the planarity hypothesis is not vacuous; and it covers the added arc ababab and all crossings, so K3,3K_{3,3}K3,3​ minus ababab is not ab-planar and the goal is not refuted by the paper's own example. A formalization that dropped the arc ababab from the drawing, or quantified over disconnecting sets instead of cuts, would state a false theorem and is ruled out.

A complete development needs basic plane topology for paths in R2\mathbb R^2R2 (a Jordan-curve-type separation lemma for simple closed curves, or an equivalent statement about crossing paths in a strip), together with combinatorial lemmas about chains (concatenation and shortcutting of chains at a common vertex). The topological lemmas are reusable well beyond this mission. Proofs through a combinatorial embedding are welcome, provided the equivalence with ABPlaneDrawing is proved.

Selected references

  • L. R. Ford, Jr. and D. R. Fulkerson, Maximal Flow Through a Network, Canadian Journal of Mathematics 8 (1956), 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • H. Whitney, Non-separable and planar graphs, Transactions of the American Mathematical Society 34 (1932), 339–362. https://doi.org/10.1090/S0002-9947-1932-1501641-2
  • A. Itai and Y. Shiloach, Maximum flow in planar networks, SIAM Journal on Computing 8 (1979), 135–150. https://doi.org/10.1137/0208012
  • H. Whitney, Planar graphs, Fundamenta Mathematicae 21 (1933), 73–84. https://doi.org/10.4064/fm-21-1-73-84
7 thms1 active userReviewed
Graph TheoryOperations Research·Captain: mikedeng1

Graph Minors. V. Excluding a Planar Graph: Bounded Tree-Width without a Planar MinorResearch Paper

Motivation

Tree-width measures how closely a graph resembles a tree. Graphs of bounded tree-width admit dynamic-programming algorithms for problems that are NP-hard in general (colouring, Hamiltonicity, and every property expressible in monadic second-order logic, by Courcelle's theorem), which is why the parameter is central to parameterized complexity and to combinatorial optimization on sparse networks. The question this mission addresses is structural: which excluded substructures force bounded tree-width?

The answer is the Excluded Grid Theorem of Robertson and Seymour: excluding a fixed graph HHH as a minor bounds the tree-width if and only if HHH is planar. The "only if" direction is easy, since grids are planar and have unbounded tree-width. The "if" direction is the content of Graph Minors. V. Excluding a Planar Graph (J. Combin. Theory Ser. B 41 (1986) 92–114). It is a cornerstone of the Graph Minors series that culminates in the Robertson–Seymour theorem (graphs are well-quasi-ordered by the minor relation), and it underlies polynomial-time minor testing for planar HHH (Graph Minors XIII) and the Erdős–Pósa-type results of Sect. 8 of the same paper.

Timeline. Robertson and Seymour, Graph Minors V (1986): tree-width at most an explicit, iterated-exponential function of the grid size. Robertson, Seymour and Thomas, Quickly excluding a planar graph (JCTB 62, 1994): bound 2O(θ5)2^{O(\theta^5)}2O(θ5) for the θ\thetaθ-grid. Chekuri and Chuzhoy (J. ACM 2016): the first polynomial bound. Chuzhoy and Tan (JCTB 2021): O(θ9 polylog θ)O(\theta^9\,\mathrm{polylog}\,\theta)O(θ9polylogθ).

Setting

Graphs are finite. A graph HHH is a minor of GGG if HHH can be obtained by contraction from a subgraph of GGG; equivalently, there are nonempty, pairwise disjoint vertex sets β(w)⊆V(G)\beta(w)\subseteq V(G)β(w)⊆V(G), one per vertex www of HHH, each inducing a connected subgraph, such that every edge ababab of HHH is matched by an edge of GGG between β(a)\beta(a)β(a) and β(b)\beta(b)β(b).

A tree-decomposition of GGG is a tree TTT together with bags Xt⊆V(G)X_t\subseteq V(G)Xt​⊆V(G) (t∈V(T)t\in V(T)t∈V(T)) such that every vertex lies in some bag, both ends of every edge lie in a common bag, and Xt∩Xt′′⊆Xt′X_t\cap X_{t''}\subseteq X_{t'}Xt​∩Xt′′​⊆Xt′​ whenever t′t't′ lies on the path of TTT between ttt and t′′t''t′′. Its width is max⁡t(∣Xt∣−1)\max_t(|X_t|-1)maxt​(∣Xt​∣−1), and the tree-width tw(G)\mathrm{tw}(G)tw(G) is the least width of a tree-decomposition of GGG.

The θ\thetaθ-grid has vertex set {vij:1≤i,j≤θ}\{v_{ij}: 1\le i,j\le\theta\}{vij​:1≤i,j≤θ}, with vijv_{ij}vij​ adjacent to vi′j′v_{i'j'}vi′j′​ exactly when ∣i−i′∣+∣j−j′∣=1|i-i'|+|j-j'|=1∣i−i′∣+∣j−j′∣=1. For even θ≥6\theta\ge 6θ≥6, Fθ\mathcal F_\thetaFθ​ is the class of graphs with no minor isomorphic to the θ\thetaθ-grid. Every planar graph HHH is a minor of some even grid of size at least 6; θ(H)\theta(H)θ(H) denotes the least such size.

The paper fixes explicit parameters. For k≥2k\ge 2k≥2: α(2,n)=n+1\alpha(2,n)=n+1α(2,n)=n+1 and α(k,n)=2nθ4+α(k−1,2nθ4+n+1)\alpha(k,n)=2^{n\theta^4}+\alpha(k-1,2^{n\theta^4}+n+1)α(k,n)=2nθ4+α(k−1,2nθ4+n+1). Then θ1=2α(θ2/2,θ2/2)\theta_1=2\alpha(\theta^2/2,\theta^2/2)θ1​=2α(θ2/2,θ2/2); ϕθ1=θ2/2\phi_{\theta_1}=\theta^2/2ϕθ1​​=θ2/2 and ϕk=ϕk+12ϕk+1θ2\phi_k=\phi_{k+1}2^{\phi_{k+1}\theta^2}ϕk​=ϕk+1​2ϕk+1​θ2; θ2=ϕ0+2ϕ1+⋯+2ϕθ1−1+ϕθ1\theta_2=\phi_0+2\phi_1+\dots+2\phi_{\theta_1-1}+\phi_{\theta_1}θ2​=ϕ0​+2ϕ1​+⋯+2ϕθ1​−1​+ϕθ1​​; θ3=(θ2/2)θ2−1\theta_3=(\theta^2/2)^{\theta_2-1}θ3​=(θ2/2)θ2​−1; θ4=θ2(θ3θ2)+12θ2(θ3θ2/2)\theta_4=\theta_2\binom{\theta_3}{\theta_2}+\tfrac12\theta^2\binom{\theta_3}{\theta^2/2}θ4​=θ2​(θ2​θ3​​)+21​θ2(θ2/2θ3​​); θ5=(θ2/2)θ4−1\theta_5=(\theta^2/2)^{\theta_4-1}θ5​=(θ2/2)θ4​−1; θ6=θ3(θ5θ4)+12θ2(θ5θ2/2)\theta_6=\theta_3\binom{\theta_5}{\theta_4}+\tfrac12\theta^2\binom{\theta_5}{\theta^2/2}θ6​=θ3​(θ4​θ5​​)+21​θ2(θ2/2θ5​​); θ7=α(θ5,θ6)\theta_7=\alpha(\theta_5,\theta_6)θ7​=α(θ5​,θ6​); θ8=3θ5(3θ5−1)/4\theta_8=3\theta_5(3^{\theta_5}-1)/4θ8​=3θ5​(3θ5​−1)/4; θ9=θ7(θ8+1)+1\theta_9=\theta_7(\theta_8+1)+1θ9​=θ7​(θ8​+1)+1.

Two auxiliary structures carry the argument. An (m,n)(m,n)(m,n)-web is a pair of families of paths (A1,…,Am)(A_1,\dots,A_m)(A1​,…,Am​), (B1,…,Bn)(B_1,\dots,B_n)(B1​,…,Bn​), each family vertex-disjoint, every AiA_iAi​ meeting every BjB_jBj​, and all m+nm+nm+n paths pairwise edge-disjoint. An (m,n)(m,n)(m,n)-mesh is the same with arbitrary connected subgraphs in place of paths and without edge-disjointness.

Formalization targets

Goal: (2.1)

For every finite planar graph HHH and every finite graph GGG,

H⪯̸G  ⟹  tw(G)≤θ9(θ(H)).H \not\preceq G \;\Longrightarrow\; \mathrm{tw}(G)\le\theta_9\bigl(\theta(H)\bigr).H⪯G⟹tw(G)≤θ9​(θ(H)).

Principal theorem: (7.3)

For even θ≥6\theta\ge 6θ≥6 and G∈FθG\in\mathcal F_\thetaG∈Fθ​,

tw(G)≤θ9.\mathrm{tw}(G)\le\theta_9.tw(G)≤θ9​.

Intermediate targets

  • Sect. 2: every planar graph is a minor of some even θ\thetaθ-grid, θ≥6\theta\ge 6θ≥6.
  • (3.2): nnn disjoint connected subgraphs meeting each of V1,…,VkV_1,\dots,V_kV1​,…,Vk​, or a hitting set of size <α(k,n)<\alpha(k,n)<α(k,n).
  • (4.1), (4.2), (4.4), (4.5), (4.6): no (θ2,θ2)(\theta_2,\theta_2)(θ2​,θ2​)-web in G∈FθG\in\mathcal F_\thetaG∈Fθ​.
  • (5.1), (5.2), (5.3): no (θ5,θ6)(\theta_5,\theta_6)(θ5​,θ6​)-mesh in G∈FθG\in\mathcal F_\thetaG∈Fθ​.
  • (6.2), (6.3), (6.4): weighted and unweighted balanced-cut lemmas valid for all graphs.
  • (7.1), (7.2): separations of order ≤θ7\le\theta_7≤θ7​ splitting V(G)V(G)V(G), or any X⊆V(G)X\subseteq V(G)X⊆V(G), in ratio 1−θ8−11-\theta_8^{-1}1−θ8−1​.

Significance

The theorem converts a qualitative exclusion (no HHH minor) into a quantitative width bound, and it is the entry point of the structure theory of minor-closed classes: every minor-closed class excluding a planar graph has bounded tree-width, and hence all MSO-definable problems on it are solvable in linear time. It is used in the proof of the graph minor theorem, in minor testing, and in the Erdős–Pósa property for planar minors (Sect. 8 of the paper).

The result is proved and classical; to our knowledge it has not been machine-checked in any proof assistant, and Mathlib has no graph minors, tree-decompositions or Menger's theorem. A formalization produces reusable definitions (branch-set minors, tree-decompositions, separations, grids) and a checked proof of the paper's explicit bound. The goal is stated with the paper's constant θ9\theta_9θ9​, not an optimized one; later improvements are stronger variants, not replacements.

Difficulty

The obvious attempt, building a tree-decomposition greedily from small separations, fails because nothing forces small balanced separations to exist. The whole argument supplies them: a graph without a large grid minor has no large mesh (5.3), and a graph with no large mesh has a balanced separation of bounded order (7.1). The step from no grid to no mesh goes through webs and spiders (Sect. 4) and relies on two results from Graph Minors I ((3.1) and (4.3) of the paper, cited without proof), which themselves depend on Menger's theorem. The balanced-cut lemmas (6.2)–(6.4) rest on Tutte's ordering of 2-connected graphs. None of this infrastructure exists in Mathlib.

Formalization scope

Graphs are Mathlib SimpleGraphs on finite types (Fintype V, DecidableEq V). The paper allows loops and multiple edges; for GGG this changes nothing, since every notion used depends only on adjacency, and for HHH in the goal it specializes the theorem to simple planar graphs. Minors use the branch-set model (IsMinor); planarity (IsPlanar) is the existence of a crossing-free drawing in R2\mathbb R^2R2, mirroring the platform's FourColor.IsPlanar. Tree-width is not defined as an infimum; "tree-width at most www" (TreewidthLE) is the existence of a tree-decomposition with all bags of size ≤w+1\le w+1≤w+1 over a finite tree. θ(H)\theta(H)θ(H) enters the goal as a hypothesis IsLeast {t | Even t ∧ 6 ≤ t ∧ IsMinor H (grid t)} θ, which is satisfiable for every planar HHH by the Sect. 2 milestone, so the goal is not vacuous. Every statement of Sects. 3–7 that mentions θ\thetaθ carries the standing assumption "θ\thetaθ even, θ≥6\theta\ge 6θ≥6" as hypotheses. Rational bounds such as (1−θ8−1)∣V(G)∣(1-\theta_8^{-1})|V(G)|(1−θ8−1​)∣V(G)∣ and 2(3k−1)−1∣V(G)∣2(3^k-1)^{-1}|V(G)|2(3k−1)−1∣V(G)∣ are compared in Q\mathbb QQ.

Two printed statements are corrected. (4.5) is printed for 0≤k<θ20\le k<\theta_20≤k<θ2​ and is stated for 0≤k<θ10\le k<\theta_10≤k<θ1​, the only range on which ϕk+1,ψk+1\phi_{k+1},\psi_{k+1}ϕk+1​,ψk+1​ are defined. (5.1) is false as printed for p=1p=1p=1, q≥1q\ge 1q≥1, so it carries the hypothesis "p=1p=1p=1 implies q=0q=0q=0"; the paper uses it only with p=θ2/2p=\theta^2/2p=θ2/2.

Needed infrastructure, reusable well beyond this mission: Menger's theorem, the Graph Minors I linkage results, Tutte's ordering of 2-connected graphs, and a library of lemmas for minors and tree-decompositions. Proofs of any milestone, of the cited results as separate theorems, and of basic API for the definitions are welcome.

Selected references

  • N. Robertson, P. D. Seymour, Graph Minors. V. Excluding a Planar Graph, J. Combin. Theory Ser. B 41 (1986) 92–114. https://doi.org/10.1016/0095-8956(86)90030-4
  • N. Robertson, P. D. Seymour, Graph Minors. I. Excluding a Forest, J. Combin. Theory Ser. B 35 (1983) 39–61. https://doi.org/10.1016/0095-8956(83)90079-5
  • N. Robertson, P. D. Seymour, R. Thomas, Quickly Excluding a Planar Graph, J. Combin. Theory Ser. B 62 (1994) 323–348. https://doi.org/10.1006/jctb.1994.1073
  • C. Chekuri, J. Chuzhoy, Polynomial Bounds for the Grid-Minor Theorem, J. ACM 63 (2016), Art. 40. https://doi.org/10.1145/2820609
  • J. Chuzhoy, Z. Tan, Towards Tight(er) Bounds for the Excluded Grid Theorem, J. Combin. Theory Ser. B 146 (2021) 219–265. https://doi.org/10.1016/j.jctb.2020.09.010
  • B. Courcelle, The Monadic Second-Order Logic of Graphs. I. Recognizable Sets of Finite Graphs, Information and Computation 85 (1990) 12–75. https://doi.org/10.1016/0890-5401(90)90043-H
28 thms1 active userReviewed
Graph Theory·Captain: Minghui

Formalize the Four Color Theorem in Lean 4Research Paper

Why formalize the Four Color Theorem in Lean 4?

The Four Color Theorem links a short mathematical statement to a large collection of finite checks. For contributors to graph theory and proof assistants, it is a concrete test of whether abstract mathematics, executable verification, and a public theorem statement can share one auditable foundation. The proposed result is a complete Lean 4 formalization, reusing Mathlib and adapting the established proof architecture to Lean.

The mathematical result is established. Robertson, Sanders, Seymour, and Thomas gave a modern computer-assisted proof in 1997. Gonthier subsequently developed a Coq formalization covering both mathematical reasoning and computation, described in his 2005 report and 2008 Notices article. This mission is classified as ResearchPaper because it formalizes these published results. (RSST, Gonthier 2005, Gonthier 2008)

Graphs, drawings, and colors

The root Lean interface uses Mathlib's SimpleGraph: a finite vertex set with a symmetric, irreflexive adjacency relation. Edges in this representation carry no additional identities or multiplicities. The conventional statement also permits parallel edges. Already-proved local infrastructure uses Mathlib's Graph V E to erase parallel edges while preserving a genuine plane drawing and proper vertex colorings; only the actual vertex set must be finite. A proper vertex coloring assigns a color to every vertex and assigns different colors to adjacent vertices. Four available colors means at most four colors are used; there is no requirement that each color occur.

A planar drawing places distinct vertices at distinct points of the real plane and represents each edge by an injective continuous arc joining its endpoints. An arc contains no other vertex, and arcs of different undirected edges meet only at shared endpoints. Reversing the orientation of an edge reverses its parameterization. A graph is planar when such a drawing exists. This is a geometric condition independent of coloring and of the eventual configuration checkers. Gonthier discusses the graph-embedding formulation in Section 2, PDF p. 4 of the 2005 report.

Write VVV for the vertex set, GGG for its adjacency relation, and C={0,1,2,3}C=\{0,1,2,3\}C={0,1,2,3} for the available colors. Empty graphs, isolated vertices, and disconnected graphs are included. The drawing is a witness to a hypothesis; the theorem does not impose coordinates, a prescribed embedding, or a straight-line representation.

Formalization target

For every finite loopless planar graph, establish

∀ G=(V,E),Planar⁡(G)⟹∃c:V→C,∀v,w∈V, {v,w}∈E⟹c(v)≠c(w).\forall\,G=(V,E),\qquad \operatorname{Planar}(G)\Longrightarrow \exists c:V\to C,\quad \forall v,w\in V,\ \{v,w\}\in E\Longrightarrow c(v)\ne c(w).∀G=(V,E),Planar(G)⟹∃c:V→C,∀v,w∈V, {v,w}∈E⟹c(v)=c(w).

The draft Lean target is FourColor.four_color: for every finite vertex type and every SimpleGraph on that type, FourColor.IsPlanar G implies G.Colorable 4. Its graph formulation follows RSST's Section 1, PDF p. 2 of the author-hosted manuscript. The manuscript's PDF page numbers differ from the journal's pagination.

The exact unchanged root signature is:

FourColor.four_color.{u} :
  ∀ (V : Type u) [Finite V] (G : SimpleGraph V),
    FourColor.IsPlanar G → G.Colorable 4

This is an open target signature, not a completed proof. The existing FourColor.FinalAudit.multigraph_four_color_of_necessary_targets conditionally connects the mission obligations to the conventional finite loopless planar multigraph statement. Its drawing has distinct vertex positions and injective continuous arcs for each edge identity; erasing parallel edges is proved to preserve planarity and every proper-coloring constraint. This bridge adds no simplicity, connectedness, triangulation, or nonemptiness assumption to the conventional conclusion. The root target itself remains unchanged.

An equivalent combinatorial formulation is welcome only with the formally verified translations needed to recover this target. If equivalence is claimed, both directions must be proved under precisely stated conventions. In particular, a hypermap coloring theorem alone does not complete this mission.

The original seven structural milestones remain unchanged:

MilestoneDeliverable
1. Graph realizationConstruct a planar plain hypermap whose faces represent exactly the nonisolated graph vertices and whose edge steps encode precisely adjacency.
2. Cubic normalizationConstruct a plain cubic hypermap with six times as many darts, preserving planarity and bridgelessness and transporting a coloring back.
3. Minimal counterexampleChoose a least-dart counterexample within the planar, bridgeless, plain, precubic comparison class.
4. Counterexample structureProve cubicity, connectedness and minimum face arity five as consequences of minimality.
5. Charge conservationProve total face charge 120c120c120c for ccc components under arbitrary rational dart transfers, and a positive-charge face when connected.
6. EliminationComplete the source's reducibility and unavoidability analysis to exclude every minimal counterexample.
7. Hypermap theoremAssemble four-colorability for all planar bridgeless hypermaps.

These correspond to Gonthier's reductions in Section 3, PDF pp. 6–9, and development in Sections 5.1–5.6 of the 2005 report. Exact locations and executable-reference declarations accompany each milestone. Milestones 2, 3 and 6 imply milestone 7; milestones 1 and 7 imply the graph target. Those two conditional implications have been checked in Lean against the exact proposed statements. Milestone 6 now has thirteen concrete supporting targets:

Core targetDeliverable
Catalogue geometryProve geometric admissibility of every one of the fixed 633 configurations.
Catalogue reducibilityKernel-check reducibility of every fixed map and contract using the proved complete checker or a verified refinement.
ReflectionProve that the explicit mirror preserves minimal counterexamples.
Geometric exclusionProve that a C-reducible configuration cannot occur in a minimal counterexample; this includes the Birkhoff and patching arguments.
Presentation soundnessProve the concrete finite presentation checker's generic soundness as one route to coverage.
Seven coverage casesIndependently handle positive hubs of degrees 5, 6, 7, 8, 9, 10 and 11, allowing reflected occurrences and unbounded neighboring arities.
Transfer boundBound each directed transfer by five under the explicit absence-of-configurations hypothesis; proved arithmetic then excludes positive hubs of degree at least 12.

The fixed catalogue and fixed rules are shared by both branches. The local Lean assembly checks that reducibility, geometric exclusion, reflection, the transfer bound, the seven degree cases, and milestones 4–5 imply milestone 6. It then connects to the unchanged hypermap and graph targets. The generic presentation checker offers a sufficient route to the semantic degree targets; its ability to certify every source presentation is not presumed. That route also requires actual presentation certificates and proofs of acceptance for the degree cases. Generic soundness and catalogue geometry alone do not supply those witnesses. Direct proofs of the seven semantic coverage obligations remain valid. Catalogue geometry remains a separate, meaningful finite theorem.

Exact mission obligations

There are 21 open theorem targets: 20 supporting milestones and one root. Every name below has prefix FourColor.. All remain future mission work; the checked conditional assembly is not a proof of any of these targets.

#Exact target nameObligation
1graph_realizationRealize a finite drawn graph by a planar plain hypermap with exact face/adjacency incidence.
2cubic_normalizationConstruct the sixfold plain cubic map with the stated preservation and coloring transport.
3minimal_counterexample_existsSelect a least-dart counterexample in the precubic comparison class.
4minimal_counterexample_structureDerive cubicity, connectedness and face arity at least five from minimality.
5charge_conservationProve total charge and existence of a positive face for a connected host.
6catalogue_embeddableProve the fixed 633 entries satisfy configuration geometry.
7catalogue_reducibility_certificatesEstablish accepted reducibility certificates for every fixed entry.
8mirror_minimal_counterexamplePreserve minimal-counterexample status under the specified mirror.
9reducible_configuration_exclusionExclude a C-reducible occurrence from a minimal counterexample.
10discharge_presentation_soundnessProve generic soundness of the concrete finite presentation checker.
11degree_5_coverageDerive a catalogue occurrence, in either orientation, from a positive degree-5 hub.
12degree_6_coverageEstablish the same semantic occurrence obligation for degree 6.
13degree_7_coverageEstablish the same semantic occurrence obligation for degree 7.
14degree_8_coverageEstablish the same semantic occurrence obligation for degree 8.
15degree_9_coverageEstablish the same semantic occurrence obligation for degree 9.
16degree_10_coverageEstablish the same semantic occurrence obligation for degree 10.
17degree_11_coverageEstablish the same semantic occurrence obligation for degree 11.
18discharge_transfer_boundBound each transfer by five under minimality and absence of catalogue occurrences in either orientation.
19no_minimal_counterexampleAssemble the core argument to exclude all minimal counterexamples.
20hypermap_four_colorAssemble four-colorability of every planar bridgeless hypermap.
21four_colorProve the unchanged finite planar graph root above.

The degree cases quantify over minimal counterexamples with the exact positive score defined by the fixed rules; neighboring degrees have no artificial upper bound. The dependency path uses 16 necessary input targets to derive targets 19, 20 and 21. Targets 6 and 10 support the optional presentation-certificate route. None of targets 19, 20 or 21 is assumed as a shortcut in this assembly.

What completion would provide

The theorem provides a uniform existence guarantee for all finite planar graphs, without a bound on their size. Its Lean development should also make reusable graph embeddings, finite combinatorial maps, coloring transports, and verified finite checkers available to later work.

The existing Rocq development is an executable reference for definitions, dependency structure, and proof behavior. It is not a proof import into Lean. The proposed contribution is a Lean development whose proof objects and computation are justified within Lean's documented foundations. A mechanically translated collection of scripts is not required; contributors should choose abstractions that work well with Mathlib.

Where the difficulty lies

Checking finitely many small graphs does not establish the theorem for arbitrary finite graphs. The development must justify why its finite computational tasks suffice, and must connect their results to the graph statement. The mathematical and computational obligations must meet at explicit, proved interfaces.

There is also a representation gap. The standard target concerns vertex colorings of graphs drawn in the plane, whereas the reference implementation's combinatorial core colors faces of planar bridgeless hypermaps. The latter endpoint is four_color_hypermap in combinatorial4ct.v. Correct treatment of duality, connected components, and isolated vertices is part of the work. Renaming a combinatorial predicate “planar” cannot establish that connection.

Formalization scope and acceptance criteria

Use Lean 4 with a supported, pinned Mathlib revision. The initial draft is checked against Lean v4.30.0 and Mathlib c5ea00351c28e24afc9f0f84379aa41082b1188f. Any later migration must preserve the statements and repeat the relevant checks. Reuse Mathlib's graph and coloring interfaces where appropriate, and develop missing infrastructure as reusable modules. The topological definition must not assume colorability, successful certificate verification, or the conclusion of an intermediate theorem.

The definition layer now provides plane drawings, finite permutation hypermaps, an exact graph/face incidence representation, minimal counterexamples, and rational face charges. Coloring equivalence across a face representation is proved for every positive number of colors, with isolated vertices handled explicitly. Exact Euler equality defines combinatorial planarity; geometric realization is a theorem obligation and cannot be assumed from the name.

The concrete core now defines configurations, ordered boundary traces, contracts, chromograms, Kempe closure, kernel preembeddings and reflected occurrences. The reference reducibility checker has proved soundness and completeness. All 633 maps, rings and contracts are literal data with kernel-validated table/index proofs. The 38 base entries and their 71 symmetrized entries are explicit, and the executable matcher has proved semantic correctness. These infrastructure results do not prove the catalogue's geometric validity, its reducibility, or unavoidability.

The production targets import FourColor.CompactCatalogue.configurations, a compact catalogue containing the same literal maps, rings and contracts. A generic verified decoder constructs validated entries without storing large evaluated proof terms in every import. A separate kernel-checked comparison proves exact ordered equality with the original catalogue, proves that every raw entry is accepted, and proves that none is dropped. The original certificate batches remain available for that audit but are excluded from production target imports. This representation change does not alter any of the 21 mathematical obligations.

Data provenance is pinned to Rocq commit c1d6b1cd5288bea4b067aac13cdde3c18dffe018. Independent fresh-source comparisons check all 633 configurations and every base/expanded rule entry against that snapshot, including actual evaluated Lean payloads. These comparisons establish source correspondence; they are not proofs of reducibility or unavoidability.

The duplicate-data gate rejects unexplained or unintended duplicates and permits verified source-mandated repetition encoding multiplicity or weight. The only repeated rule payload is drule1, deliberately listed twice: the source explains its two-point transfer in discharge.v, lines 17–24 and retains both copies in base_drules (line 145). The Lean rule-match count also counts both copies, preserving weight two. Thus there are 38 base entries with 37 distinct payloads, and 71 expanded entries with 70 distinct payloads. Both copies must remain. No configuration duplicates or other repeated rule payloads are permitted without separately verified source justification.

The reference reducibility evaluator is exponentially expensive. Contributors should implement verified compressed representations or efficient checkers, with acceptance implying the same semantic C-reducibility. Completeness of the reference checker then recovers the stated certificate-existence target. The presentation language is a transparent finite baseline; its generic soundness is an open target, and no completeness claim is made. Source-style quizzes/hubcaps or direct Lean proofs can close the same seven semantic coverage targets. This keeps performance choices separate from the mathematical endpoint.

All computational claims needed by the theorem must be established in Lean. External generators may prepare candidate data or certificates, but their output must be validated by a checker with proved correctness, and the resulting proof must be accepted by the Lean kernel. A recorded successful run of an external program is insufficient. The final theorem and its dependency closure must contain no sorry, admit, unproved custom axiom, or unsupported computational assumption. An axiom audit must identify only Lean/Mathlib's documented foundations. The current 37-declaration conditional/infrastructure audit uses only propext, Classical.choice and Quot.sound; it does not claim an unconditional Four Color proof. Python generators and runtime JSON extraction are outside the mathematical trust boundary and cannot discharge the 21 targets.

Completion requires a reproducible repository in which lake build succeeds from a clean environment and builds the final theorem. Pin toolchain, dependencies, source revisions, and certificate data; document regeneration and verification commands. Maintain a dependency map connecting Lean declarations to exact paper locations and Rocq declarations, and record justified differences in representation. The local multigraph adapter already justifies forgetting parallel-edge multiplicities and transports proper colorings. Any auxiliary restrictions such as connectedness or nonemptiness must be discharged before the root theorem. The current statement and conditional builds verify proposal infrastructure. Success requires proofs of the mission obligations and the final theorem in the reproducible clean build, with no unproved assumptions beyond the documented Lean/Mathlib foundations. All 21 targets remain open at proposal finalization. The submission snapshot records the repository base commit, exact working-file hashes, toolchain, Mathlib revision, payload hashes and audit report; it must not misrepresent uncommitted files as contents of the base commit.

Selected references

  • Georges Gonthier, A Computer-Checked Proof of the Four Colour Theorem, technical report, 2005. Sections 2–5. Microsoft Research PDF.
  • Georges Gonthier, Formal Proof—The Four-Color Theorem, Notices of the AMS 55(11), 1382–1393, 2008. Theorem 1 and the formalization architecture. AMS PDF.
  • Gonthier and Rocq-community contributors, fourcolor, executable formalization, pinned to commit c1d6b1cd5288bea4b067aac13cdde3c18dffe018. Repository.
  • Neil Robertson, Daniel Sanders, Paul Seymour, Robin Thomas, The Four-Colour Theorem, Journal of Combinatorial Theory, Series B 70(1), 2–44, 1997. DOI; author-hosted manuscript.
54 thms1 active userReviewed
Operations ResearchOptimization·Captain: mikedeng1

Branch-and-Price-and-Cut for the Split-Delivery Vehicle Routing Problem with Time Windows: Some Optimal Solution Traverses Each Pair of Reverse Customer Arcs at Most OnceResearch Paper

Motivation

Vehicle routing problems ask for minimum-cost vehicle routes that deliver goods from a depot to a set of customers. In the split-delivery vehicle routing problem with time windows (SDVRPTW) a customer's demand may be served by several vehicles, and each customer must be visited inside a prescribed time window. Allowing split deliveries matters in practice: for the variant without time windows, Dror and Trudeau (1989) showed empirically that splitting can save substantially, and Archetti, Savelsbergh and Speranza (2006) proved the savings can reach 50%.

Exact methods for the SDVRPTW are branch-and-price algorithms, and they rely on structural properties of optimal solutions to prune the search. Desaulniers (Operations Research 58(1), 2010) collects these properties in §2 of his paper and uses the strongest one, Corollary 2, as a family of valid inequalities (constraint (7)) in his branch-and-price-and-cut method.

Timeline (as reported by Desaulniers 2010, §§1–2).

  • 1989–1990. Dror and Trudeau prove, for the SDVRP without time windows and with the triangle inequality, that some optimal solution has no two routes sharing more than one split customer.
  • 2006. Gendreau, Dejax, Feillet and Gueguen observe that the property holds with time windows, derive the arc-based Corollary 1, and remark that elementary routes suffice (Remark 1).
  • 2010. Desaulniers strengthens Corollary 1 to pairs of reverse arcs (Corollary 2) and exploits it as cutting planes.

Setting

An instance has nnn customers N\mathcal NN, a start depot 000 and an end depot n+1n+1n+1 (the same location at the beginning and end of the planning horizon), and the following data: a vehicle capacity Q>0Q > 0Q>0; a demand di>0d_i > 0di​>0 for each customer, which may exceed QQQ; a time window [ev,lv][e_v, l_v][ev​,lv​] for each node, shared by the two depot copies; nonnegative travel times tvwt_{vw}tvw​, which include the service time at vvv; and nonnegative costs cvwc_{vw}cvw​. The arc set A\mathcal AA contains the idle arc (0,n+1)(0, n+1)(0,n+1) and every arc (v,w)(v, w)(v,w), v≠wv \ne wv=w, with ev+tvw≤lwe_v + t_{vw} \le l_wev​+tvw​≤lw​. The triangle inequality tvx≤tvw+twxt_{vx} \le t_{vw} + t_{wx}tvx​≤tvw​+twx​, cvx≤cvw+cwxc_{vx} \le c_{vw} + c_{wx}cvx​≤cvw​+cwx​ is assumed throughout.

A route is a walk 0→v1→⋯→vm→n+10 \to v_1 \to \dots \to v_m \to n+10→v1​→⋯→vm​→n+1 along arcs of A\mathcal AA, customers possibly repeated, with service start times inside the time windows that respect travel times (waiting is allowed), and nonnegative quantities delivered at its visits whose total is at most QQQ. Its cost is the sum of its arc costs. A solution is a finite family of routes, one per vehicle, with no bound on their number; it is feasible if every customer iii receives in total at least did_idi​, and optimal if it is feasible and no feasible solution costs less. For customers i,ji, ji,j, let xijx_{ij}xij​ be the number of times arc (i,j)(i, j)(i,j) is traversed, summed over all routes of a solution, and let A(N)=A∩(N×N)\mathcal A(\mathcal N) = \mathcal A \cap (\mathcal N \times \mathcal N)A(N)=A∩(N×N).

Formalization targets

All four statements assume the triangle inequality and that the instance has a feasible solution, and assert the existence of an optimal solution with a structural property.

Goal: Corollary 2 (p. 181)

∃ optimal solution with xij+xji≤1for all (i,j)∈A(N).\exists \text{ optimal solution with } x_{ij} + x_{ji} \le 1 \quad \text{for all } (i,j) \in \mathcal A(\mathcal N).∃ optimal solution with xij​+xji​≤1for all (i,j)∈A(N).

The two arcs of a pair of reverse customer arcs are used at most once in total. This is the form constraint (7) of the paper gives to the corollary.

Milestones

  1. Remark 1. Some optimal solution has only elementary routes: no route visits a customer twice.
  2. Theorem 1. Some optimal solution has no two distinct routes with two customers in common.
  3. Corollary 1. Some optimal solution has xij≤1x_{ij} \le 1xij​≤1 for all (i,j)∈A(N)(i, j) \in \mathcal A(\mathcal N)(i,j)∈A(N).

Significance

The result. Corollary 2 turns a property of optimal solutions into linear inequalities on arc-flow variables. Adding them to the arc-flow formulation cuts off fractional points of its linear relaxation while keeping an optimal integer solution, which is how Desaulniers uses them. Remark 1 justifies pricing only elementary routes in column generation. Theorem 1 is the combinatorial fact underneath both corollaries and is reused throughout the split-delivery literature.

Formalizing it. The four results are proved in the literature (Dror–Trudeau; Gendreau et al. 2006); Desaulniers states them without proof. No machine-checked version exists. This mission produces a reusable Lean model of the SDVRPTW (instances, arc sets, feasible routes with schedules and delivery patterns, solutions, optimality) and checked proofs of these properties, including the existence of an optimal solution, which the paper takes for granted.

Difficulty

The obvious argument is local: take an optimal solution that violates the property, shift quantities between two routes, remove a visit, and shortcut. Three points make this less routine than it sounds. First, the statements are existential: each exchange must keep the solution optimal and not reintroduce a violation already removed, so one needs a termination measure that decreases under every exchange (Corollary 2 needs elementarity and the Theorem 1 property simultaneously, not two separate optimal solutions). Second, removing a visit is feasible only because the arc set is defined by time windows: the shortcut arc (v,w)(v, w)(v,w) must be shown to exist from the schedule and the triangle inequality on travel times, and the new schedule must be built explicitly. Third, the feasible set is infinite (real quantities, unbounded walks, unbounded number of vehicles), so the existence of an optimal solution is itself a statement to prove, not a hypothesis.

Formalization scope

Namespace SplitDeliveryVRPTW.Known. Nodes are the inductive type Node n (start, cust i for i : Fin n, finish). All quantities, times and costs are real. The arc set is exactly the set the paper defines (read as "if and only if"), with arcs into the start depot and out of the end depot excluded. The triangle inequality for ttt is imposed on pairwise distinct nodes and for ccc on triples of arcs, where the paper's data is defined. A route is a customer list (repetitions allowed) with a time function over path positions and a quantity function over visits. A solution is an indexed family Fin m → Route I, so identical routes may appear twice. Demand satisfaction uses ≥\ge≥, as constraint (2) does. The per-vehicle bound min⁡{di,Q}\min\{d_i, Q\}min{di​,Q} of constraint (14) is omitted because it changes neither the feasible route patterns nor the costs.

Explicit readings of imprecise phrases:

  • "split customer", "in common" (Theorem 1) are undefined in the paper. A customer two distinct routes visit is split, and the routes have it in common. This visit-based reading is at least as strong as a delivery-based one.
  • Corollary 2's wording "at most one arc in set Aij∗\mathcal A^*_{ij}Aij∗​ appears at most once" is read through constraint (7): the total number of traversals of the arcs of Aij∗\mathcal A^*_{ij}Aij∗​ is at most one. It is stated in the equivalent form free of the choice of A∗(N)\mathcal A^*(\mathcal N)A∗(N).
  • "there exists an optimal solution" is conditional on feasibility, which is the hypothesis added.

Ruled-out trivializations: routes are not restricted to elementary walks (that would make Remark 1 definitional), solutions are not sets (that would forbid duplicate routes), optimality compares against solutions with any number of routes and any visit pattern, "split" is never counted over visits by the same route, capacity and time windows are part of route feasibility, and deliveries occur only at visits.

Useful infrastructure, reusable beyond this mission: shortcut lemmas for feasible routes (removing a visit), exchange lemmas between two routes, and existence of an optimum for split-delivery routing. Contributions of any of these as separate lemmas are welcome.

Selected references

  • G. Desaulniers, Branch-and-Price-and-Cut for the Split-Delivery Vehicle Routing Problem with Time Windows, Operations Research 58(1):179–192, 2010. https://doi.org/10.1287/opre.1090.0713
  • M. Dror, P. Trudeau, Savings by split delivery routing, Transportation Science 23(2):141–145, 1989. https://doi.org/10.1287/trsc.23.2.141
  • M. Dror, P. Trudeau, Split delivery routing, Naval Research Logistics 37(3):383–402, 1990. https://doi.org/10.1002/nav.3800370304
  • M. Gendreau, P. Dejax, D. Feillet, C. Gueguen, Vehicle routing with time windows and split deliveries, Technical Report 2006-851, Laboratoire Informatique d'Avignon, 2006.
  • C. Archetti, M. W. P. Savelsbergh, M. G. Speranza, Worst-case analysis for split delivery vehicle routing problems, Transportation Science 40(2):226–234, 2006. https://doi.org/10.1287/trsc.1050.0117
6 thms1 active userReviewed
Graph TheoryOperations ResearchOptimization·Captain: mikedeng1

Scheduling Subject to Resource Constraints: Classification and Complexity III: The Two-Machine Algorithm for Q2 with One Resource and Unit-Time Jobs Is OptimalResearch Paper

Motivation

Many production and computing systems run jobs on parallel machines that also draw on a shared, limited resource: tools, workers, memory, power. Adding such a resource to a scheduling problem can change its complexity entirely. Błażewicz, Lenstra and Rinnooy Kan (DAM 1983) extended the three-field classification α∣β∣γ\alpha\mid\beta\mid\gammaα∣β∣γ of Graham, Lawler, Lenstra and Rinnooy Kan by a resource field resλσρres\lambda\sigma\rhoresλσρ, and determined the complexity of every problem with unit-time jobs on identical or uniform machines under the makespan criterion. Their Fig. 2 separates the maximal polynomially solvable cases from the minimal NP-hard ones.

This mission formalizes the polynomial side. Two identical machines are easy under arbitrary resources (Theorem 1, due to Garey and Johnson, via maximum matching). Three identical machines with one resource are NP-hard in the strong sense (Theorem 4), and so are two uniform machines with unit resources (Theorem 3). What remains for uniform machines is settled by two algorithms: a sorting-and-shifting procedure for two uniform machines with one resource of arbitrary size (Theorem 5), and a bottleneck transportation problem for any number of uniform machines with one resource and 0–1 requirements (Theorem 6). The hardness results are the subject of the companion missions I and II.

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​. Machine MiM_iMi​ has speed qi>0q_i>0qi​>0; every job has unit execution requirement, so it takes time 1/qi1/q_i1/qi​ on MiM_iMi​. Identical machines (PPP) have qi=1q_i=1qi​=1; uniform machines (QQQ) have arbitrary speeds. There are lll resources RhR_hRh​ with positive integer sizes shs_hsh​, and job JjJ_jJj​ needs a nonnegative integer amount rhjr_{hj}rhj​ of RhR_hRh​ throughout its execution. The field resλσρres\lambda\sigma\rhoresλσρ records restrictions: λ\lambdaλ bounds the number of resources, σ\sigmaσ their sizes, ρ\rhoρ the requirements, a dot meaning "part of the input". So res1⋅⋅res1{\cdot}{\cdot}res1⋅⋅ is one resource with arbitrary size and requirements, and res1⋅1res1{\cdot}1res1⋅1 is one resource with requirements in {0,1}\{0,1\}{0,1}.

A schedule gives every job a machine μ(j)\mu(j)μ(j) and a start time Sj≥0S_j\ge0Sj​≥0; the job is executed during [Sj,Cj)[S_j,C_j)[Sj​,Cj​) with Cj=Sj+1/qμ(j)C_j=S_j+1/q_{\mu(j)}Cj​=Sj​+1/qμ(j)​. It is feasible if jobs on the same machine do not overlap and, at every time ttt, the jobs executed at ttt use at most shs_hsh​ of each resource RhR_hRh​. The makespan is Cmax⁡=max⁡jCjC_{\max}=\max_j C_jCmax​=maxj​Cj​. No precedence constraints occur in this mission.

Formalization targets

Goal: Theorem 5, correctness of the algorithm

For Q2∣res1⋅⋅, pj=1∣Cmax⁡Q2\mid res1{\cdot}{\cdot},\,p_j=1\mid C_{\max}Q2∣res1⋅⋅,pj​=1∣Cmax​ with q1≥q2q_1\ge q_2q1​≥q2​: put all jobs on M1M_1M1​ in order of nonincreasing r1jr_{1j}r1j​, then repeatedly move the last job of M1M_1M1​ to the earliest feasible time on M2M_2M2​ after the jobs already there, as long as this strictly reduces Cmax⁡C_{\max}Cmax​. For every order with nonincreasing requirements, the resulting schedule AAA is feasible and

Cmax⁡(A)≤Cmax⁡(σ)for every feasible schedule σ.C_{\max}(A)\le C_{\max}(\sigma)\quad\text{for every feasible schedule }\sigma.Cmax​(A)≤Cmax​(σ)for every feasible schedule σ.

Milestones for the goal

The paper's proof has two steps, both milestones. Call a schedule an (a)–(c) schedule when (a) M1M_1M1​ runs its jobs back to back from time 000 in nonincreasing r1jr_{1j}r1j​, (b) M2M_2M2​ runs its jobs in nondecreasing r1kr_{1k}r1k​, and (c) every requirement on M1M_1M1​ is at least every requirement on M2M_2M2​.

  1. The algorithm's schedule is feasible, is an (a)–(c) schedule, and is best among feasible (a)–(c) schedules.
  2. Every feasible schedule can be transformed into a feasible (a)–(c) schedule with no larger Cmax⁡C_{\max}Cmax​.

Further results

  • Theorem 1. For P2∣res⋅⋅⋅, pj=1∣Cmax⁡P2\mid res{\cdot}{\cdot}{\cdot},\,p_j=1\mid C_{\max}P2∣res⋅⋅⋅,pj​=1∣Cmax​, with GGG the graph joining two jobs when they can run together and SSS a maximum matching of GGG, the optimal makespan is n−∣S∣n-|S|n−∣S∣.
  • Theorem 6. For Q∣res1⋅1, pj=1∣Cmax⁡Q\mid res1{\cdot}1,\,p_j=1\mid C_{\max}Q∣res1⋅1,pj​=1∣Cmax​ with the s1s_1s1​ fastest machines listed first, the optimal makespan equals the optimal value of a bottleneck transportation problem that assigns jobs to slots (machine, position) with cost k/qik/q_ik/qi​, resource jobs only to the s1s_1s1​ fastest machines.

Significance

Theorems 5 and 6 complete the classification of Fig. 2 for uniform machines: every special case of Q∣res⋅⋅⋅, pj=1∣Cmax⁡Q\mid res{\cdot}{\cdot}{\cdot},\,p_j=1\mid C_{\max}Q∣res⋅⋅⋅,pj​=1∣Cmax​ not covered by the hardness theorems has a polynomial algorithm. Theorem 1 is the classical reduction of two-machine resource scheduling to maximum matching, the model case for later work on scheduling with conflict graphs.

The paper proves these results briefly: "clearly" for the first half of Theorem 5, "obviously" for Theorem 1, and a one-paragraph model for Theorem 6. The exchange argument of Theorem 5 is presented "in an informal way" through five steps that pass through fractional, preempted jobs. A machine-checked proof makes these arguments exact on a model with real start times. No formalization of these results is known, and the platform had no statement about resource-constrained scheduling on uniform machines before this mission.

Difficulty

With q1≠q2q_1\ne q_2q1​=q2​ the job boundaries on the two machines are misaligned: a job on M2M_2M2​ overlaps parts of several jobs on M1M_1M1​, so the resource check cannot be done slot by slot, and discrete reasoning on integer time grids does not apply. The exchange argument of Theorem 5 must control the resource usage at every real time while jobs are moved between machines and reordered, and it has to end with a nonpreemptive schedule even though the paper's intermediate steps split jobs. For Theorem 1, the hard direction is the lower bound: a feasible schedule with arbitrary real start times must be converted into a matching, which is a statement about how unit jobs on two machines can overlap. For Theorem 6, one must show that restricting resource jobs to the fastest machines and to back-to-back positions loses nothing.

Formalization scope

  • Model. Jobs, machines and resources are Fin n, Fin m, Fin l (0-based). Speeds are positive reals, sizes positive naturals, requirements naturals. Start times are nonnegative reals, execution intervals are half-open, and the resource constraint is checked at every real time. Cmax⁡=0C_{\max}=0Cmax​=0 for n=0n=0n=0. The model carries a precedence digraph for consistency with the companion missions; every statement here assumes it has no arcs.
  • Implicit hypothesis. Theorems 1 and 5 assume every job fits alone (rhj≤shr_{hj}\le s_hrhj​≤sh​), which the paper leaves unstated; without it no feasible schedule exists.
  • The algorithm is a Lean definition following the page: the order is an argument (any nonincreasing order), "as early as possible" is the earliest start after M2M_2M2​'s last job at which the resource constraint holds throughout, and the loop stops at the first move that does not strictly reduce Cmax⁡C_{\max}Cmax​.
  • Optimality is always stated in full: feasibility plus a lower bound against every feasible schedule. No minimum is written as an infimum of a possibly empty set.
  • Theorem 6 is stated with 0–1 slot assignments, the interpretation the paper gives to xijkx_{ijk}xijk​; the page's constraint ∑k=1m\sum_{k=1}^{m}∑k=1m​ is read as ∑k=1n\sum_{k=1}^{n}∑k=1n​.
  • Not formalized: the running times O(ln2+n5/2)O(ln^2+n^{5/2})O(ln2+n5/2) (Theorem 1), O(nlog⁡n)O(n\log n)O(nlogn) (Theorem 5, including the phrase "This O(n log n) algorithm") and O(n3)O(n^3)O(n3) (Theorem 6), which depend on a machine model the paper does not fix and, for Theorems 1 and 6, on cited matching and transportation algorithms.
  • Ruled out: a formalization of the goal that proves optimality only against (a)–(c) schedules, against schedules with integer start times, or for one fixed tie-breaking order proves less than Theorem 5.

Contributions welcome: lemmas about step functions of resource usage on half-open intervals, a left-shifting lemma for unit-time schedules on two machines, and the exchange steps of Theorem 5 as separate lemmas.

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
  • 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. Even, O. Kariv, An O(n^{2.5}) algorithm for maximum matching in general graphs, Proc. 16th IEEE FOCS (1975) 100–112. https://doi.org/10.1109/SFCS.1975.23
10 thms1 active userReviewed
Dynamical SystemsFormal Verification·Captain: Rizwan G Mir

Bound L4: 33,070,982 <= R Reversible Binary 2D Moore RulesOpen Problem

Bound L4L_4L4​: 33,070,982≤R33,070,982 \le R33,070,982≤R

This mission formalizes the lower bound 33,070,982≤R33,070,982 \le R33,070,982≤R on the number of reversible binary cellular automata on the 3×33 \times 33×3 Moore neighborhood. Extending the conserved-landscape marker families to both centered and off-centered rules.

1 thm1 active userReviewed
Discrete Geometry·Captain: mikedeng1

Convex Polytopes V.5: Effective enumeration of combinatorial typesTextbook

A convex polytope has both geometric coordinates and a finite pattern of faces. Moving its vertices can change distances and angles while preserving which vertices belong to which faces. Enumeration by combinatorial type asks for those incidence patterns, with geometrically different realizations of the same pattern counted together. Grünbaum’s enumeration theorem establishes that the complete collection can be determined algorithmically when the dimension and number of vertices are prescribed Convex Polytopes, §5.5, p.91.

A d-polytope here is a nonempty convex hull of finitely many points in real d-dimensional coordinate space, with affine span equal to the whole space. A face is obtained by maximizing a linear functional over the polytope; the collection of faces also includes the empty face and the polytope itself. Inclusion orders this collection. Two polytopes have the same combinatorial type when their full face collections admit an order isomorphism. This notion retains the incidence structure and forgets metric measurements.

To describe a type with finite data, label the k vertices by the integers from zero through k−1. Record the set of vertex labels belonging to each nonempty proper face. This family of subsets is the polytope’s scheme. A finite list of finite lists of natural numbers encodes such a family. The order of labels and repeated occurrences of a label do not change the represented subset. A scheme is realized only when its recorded subsets are exactly the vertex sets of the nonempty proper faces: it cannot add faces, omit faces, or introduce labels outside the prescribed range.

The target is a single total computable function

E:N×N⟶List⁡(List⁡(List⁡(N))).E : \mathbb N\times\mathbb N\longrightarrow \operatorname{List}(\operatorname{List}(\operatorname{List}(\mathbb N))).E:N×N⟶List(List(List(N))).

For every dimension d and vertex count k, each scheme in E(d,k) must have a realizing d-polytope with exactly k vertices. Every d-polytope with k vertices must realize some scheme in that output. Finally, if polytopes realizing two output positions have isomorphic full face posets, those positions must be equal. These requirements say that the output contains precisely one representative of every combinatorial type. The existential choice of one function precedes both numerical inputs, so the algorithm must work uniformly for all dimensions and vertex counts.

The result supplies an effective finite classification at each prescribed size. It is stronger than the observation that only finitely many incidence families can be written down: those families need not all arise from real convex polytopes. It also addresses termination, since a total algorithm must return its entire finite answer for every input. The theorem is a known mathematical result in the cited textbook. The present formal target asks for a proof of its stated algorithmic conclusion; no machine-checked proof of that conclusion is asserted here.

The principal difficulty is the connection between finite incidence data and geometric realizability. Combinatorial consistency alone does not supply real coordinates for a convex polytope. The source distinguishes this realizability question from the enumeration conclusion and identifies decidability over the real numbers as relevant to it §5.5, p.91. The mission retains the complete enumeration conclusion, including realizability of each answer, coverage of every type, and absence of duplicate types.

The formal representation uses real coordinates without a rationality restriction and permits nonsimplicial polytopes. Vertices are labelled injectively and exhaustively. The number of vertices is expressed as the number of nonempty zero-dimensional exposed faces. Nonemptiness separates the empty face, whose book dimension is −1, from the natural-valued dimension used to count faces. For finite convex hulls the face collections are finite, so this cardinality has its usual meaning.

Natural-number dimensions include dimension zero. The unique point has one vertex and no nonempty proper faces, and therefore uses an empty scheme. This is an explicit extension of the source’s positive-dimensional scheme convention. Unrealizable dimension and vertex-count pairs require an empty output. Finite-polytope faces, their ordered collection, vertex counts, and finite scheme realizations are the concrete objects needed to state the goal; total computability applies to the complete finite output rather than to individual tests alone.

Reference: Branko Grünbaum, Convex Polytopes, second edition, Springer, 2003, §5.5, Theorem 2 (5.5.2), printed p.91, PDF p.117; definitions in §§2.4 and 3.1. Source text.

5 thms1 active userReviewed
Discrete Geometry·Captain: mikedeng1

Convex Polytopes IV: Facet growth of binary polytopes (GRU-M04-BINARY-EXTREMA)Textbook

Binary choices and polyhedral complexity

A vector whose coordinates are all zero or one records a collection of yes-or-no choices. Taking the convex hull of a collection of these vectors gives a 0/1-polytope, also called a binary polytope. Such objects connect discrete choices with geometry: points describe feasible combinations, while supporting inequalities describe restrictions that every feasible combination satisfies. Grünbaum's discussion in §4.9 emphasizes the role of facets of these polytopes in combinatorial optimization and cutting-plane descriptions.

The question here concerns how many facets can occur when the dimension grows. Restricting every generating point to binary coordinates gives a finite, highly structured collection of possible points. Nevertheless, this restriction still allows polytopes with very many distinct boundary faces. The theorem of Bárány and Pór establishes superexponential growth in the largest possible number of facets. The mission concerns the form of this growth asserted in Grünbaum's 2003 notes, printed page 69a, rather than a prescribed numerical constant.

Polytopes, dimension, and facets

Fix a natural number d. The ambient space is the real coordinate space R^d. A convex hull consists of all convex combinations of the generating points, meaning weighted averages with nonnegative weights whose sum is one. A polytope is the convex hull of a finite set of points. It is full-dimensional when its affine span is all of R^d; this rules out a polytope lying in a proper affine subspace.

For a binary polytope, every generating point belongs to {0,1}^d. A subset of this cube is automatically finite. Nonemptiness and full dimension are imposed separately in the target. Consequently, the dimension in the facet bound is the actual dimension of the polytope as well as the dimension of its coordinate space.

An exposed face is the set of all points of the polytope that maximize a given linear functional. A facet is a face of affine dimension d−1. Facets are counted as geometric sets. Two different inequalities that expose the same face do not contribute two facets. Write f_(d−1)(P) for this number.

The growth target

The target is the following existence statement:

∃c>1  ∃D∈N, D≥2,∀d≥D  ∃P⊆Rd,P is a full-dimensional binary polytope,fd−1(P)>cdln⁡d.\exists c>1\;\exists D\in\mathbb N,\ D\ge2,\quad \forall d\ge D\;\exists P\subseteq\mathbb R^d,\quad P\text{ is a full-dimensional binary polytope},\qquad f_{d-1}(P)>c^{d\ln d}.∃c>1∃D∈N, D≥2,∀d≥D∃P⊆Rd,P is a full-dimensional binary polytope,fd−1​(P)>cdlnd.

The constant c and threshold D are chosen before the dimension d. The polytope P may depend on d. The assertion therefore supplies a witness in every sufficiently large dimension. It does not merely assert the existence of one complicated polytope or of witnesses along an unspecified subsequence.

The logarithm is natural. Replacing it by another fixed base greater than one changes the admissible constant c, while preserving the shape of the theorem. The statement leaves both c and D unspecified. The Bárány–Pór paper, Theorem 1.1, supplies the asymptotic context for the book's formulation.

What superexponential growth says

For any fixed c greater than one, the expression c^(d ln d) eventually exceeds A^d for every fixed A greater than one. Thus a single exponential base cannot bound the facet counts of all binary polytopes across dimensions. The binary-coordinate restriction alone does not yield that kind of uniform bound.

The underlying mathematical result is established in the literature. The formalization task is to prove its stated existence conclusion with the precise geometric definitions above. The conclusion concerns actual facets of actual polytopes, so a family of redundant inequalities or a list with repetitions would not satisfy the counting requirement.

Why the existence statement is demanding

A large supply of binary points does not by itself identify the supporting hyperplanes of their convex hull. Counting points and counting facets are different tasks. Moreover, the theorem requires the same growth constant across all sufficiently large dimensions. Verifying individual examples, even examples with many facets, leaves that uniform asymptotic requirement unresolved.

The book describes certain random polytopes as witnesses. Its stated conclusion gives no probability distribution or numerical probability bound. The target records the resulting extremal existence assertion; it imposes no extra probabilistic hypothesis on the witness.

Geometric conventions

The coordinate model is Fin d → ℝ. IsDPolytope requires a nonempty finite convex hull with full affine span. IsZeroOnePolytope specifies the hull of binary-coordinate points. faceCount P k counts nonempty exposed faces with affine dimension k. These concrete notions also apply to other questions about finite-dimensional polytopes.

Face counts use natural cardinality. On the domain of finite polytopes there are only finitely many faces, so this is the ordinary finite count. Requiring D at least two ensures that d−1 represents the facet dimension without a low-dimensional subtraction convention and that the logarithm's argument is positive. Full dimension excludes the whole polytope from the facet count, and nonemptiness excludes the empty face.

Selected references

  • Branko Grünbaum, Convex Polytopes, second edition, Springer, 2003, §4.9, printed p.69a; definitions in §§2.4 and 3.1. Book.
  • Imre Bárány and Attila Pór, On 0-1 Polytopes with Many Facets, Advances in Mathematics 161 (2001), 209–228, Theorem 1.1. Paper.
2 thms1 active userReviewed
Discrete Geometry·Captain: mikedeng1

Convex Polytopes V: Recognition from projections (GRU-M05-PROJECTION-RECOGNITION)Textbook

Recognizing a convex set from its projections

A convex set can be studied through its images in spaces of smaller dimension. Each image records part of the geometry, while losing the information along the directions that are collapsed. A finite convex hull always has finite convex hulls as its affine images. The converse asks whether sufficiently rich projection data can force the original set to have a finite description by points. This mission concerns the recognition theorem attributed to Klee in Grünbaum's Convex Polytopes, §5.1, Theorem 8, printed page 74.

The question concerns all projections into a suitable dimension, rather than an individual view of a set. A single image may conceal directions in which the original set has additional structure. The theorem identifies a condition on the entire family of images that characterizes polytopes among bounded convex sets.

Sets, convex hulls, and affine projections

Fix an integer d≥3d\geq3d≥3. Real coordinate space Rd\mathbb R^dRd consists of vectors with ddd real coordinates. A set K⊆RdK\subseteq\mathbb R^dK⊆Rd is convex if it contains the line segment joining any two of its points. It is bounded if its points remain within some finite distance of the origin. Neither condition requires KKK to fill the ambient space.

The convex hull of a set VVV, written conv⁡(V)\operatorname{conv}(V)conv(V), is the smallest convex set containing VVV. A polytope is the convex hull of a finite set. The generating set need not be a minimal set of vertices. Grünbaum gives this characterization in §3.1, printed page 31. The empty generating set is permitted, as are generating sets contained in a proper affine subspace.

An affine map preserves affine combinations. A surjective affine map f:Rd→Rjf:\mathbb R^d\to\mathbb R^jf:Rd→Rj has a jjj-dimensional target and reaches every point of that target. When j<dj<dj<d, such a map loses dimensions. This represents the singular affine images called projections in §5.1, printed page 71, with coordinates chosen on the target affine space. Translation of the target is allowed.

The recognition target

For every bounded convex set K⊆RdK\subseteq\mathbb R^dK⊆Rd, the goal is the equivalence

K is a polytope⟺∃j∈N,2≤j<d,∀f:Rd↠Rj affine,f(K) is a polytope.K\text{ is a polytope} \quad\Longleftrightarrow\quad \exists j\in\mathbb N,\quad 2\leq j<d,\quad \forall f:\mathbb R^d\twoheadrightarrow\mathbb R^j\text{ affine},\quad f(K)\text{ is a polytope}.K is a polytope⟺∃j∈N,2≤j<d,∀f:Rd↠Rj affine,f(K) is a polytope.

The dimension jjj may depend on KKK, but it is chosen before testing the affine maps. Every surjective affine map with that target dimension is tested. For each such map, the finite set generating the image can be different. No common set of projected vertices or uniform bound on the number of generators is required.

This is one recognition theorem, containing both implications. The selected source grouping does not introduce separate supporting theorem targets.

What the criterion establishes

The theorem characterizes a global finite convex-hull property through lower-dimensional images. In particular, the conclusion concerns the original set itself; it does not merely assert that its closure has a finite generating set. This distinction matters because boundedness and convexity alone do not assert closedness.

The result is a known mathematical theorem in the source. The formalization target is its complete equivalence with the stated domain and quantifier order. A proof of only the preservation of polytopes under affine maps would leave the recognition implication unresolved.

Why the converse requires more than one image

Every tested image can have its own finite generating set. Finiteness of each image does not directly supply a finite generating set that works in the original ambient space. The recognition implication must connect the universal family of images to the geometry of the whole set. Replacing that family by a convenient fixed projection would change the question.

Mathematical conventions

Real ddd-space is represented by functions from Fin d to the real numbers. Boundedness and convexity use the ordinary Mathlib predicates, and being a polytope is expressed directly by existence of a finite set with the specified convex hull. No full-dimensional polytope structure is imposed on KKK.

The hypothesis d≥3d\geq3d≥3 makes explicit the admissible ambient dimension needed for 2≤j<d2\leq j<d2≤j<d. In dimension three, the only available target dimension is two. No assertion of the displayed existential criterion is made in dimensions zero, one, or two. Empty sets, singleton sets, and other lower-dimensional subsets remain within the domain in every admissible ambient dimension.

Affine maps are required to be surjective onto the selected coordinate space. Thus their rank is exactly the selected target dimension, rather than an accidentally smaller rank. Closedness, nonemptiness, rational coordinates, and a predetermined number of vertices are not additional hypotheses. The mathematical ingredients are real convex hulls, finite sets, bounded sets, and affine images.

Selected references

Branko Grünbaum, Convex Polytopes, second edition, Springer, 2003. Theorem 5.1.8, printed page 74 (source.pdf page 100); projection convention, printed page 71 (PDF97); polytope convention, printed page 31 (PDF51). The source pages are included in this package's evidence directory.

1 thm1 active userReviewed
Discrete Geometry·Captain: mikedeng1

Convex Polytopes V: Simplex sections through prescribed points (GRU-M05-PRESCRIBED-SECTIONS)Textbook

A polytope can be represented by cutting a higher-dimensional simplex with an affine flat. The geometry of the cut determines the resulting polytope. Perles's prescribed-point theorem adds a constraint to this representation: the cutting flat must pass through a point chosen in advance inside the simplex. Theorem 5.1.10 in Branko Grünbaum's Convex Polytopes, second edition (2003), shows that this requirement can always be met when the simplex has enough facets. The source is §5.1, printed page 74, with the section convention introduced on printed page 71.

A polytope is the convex hull of finitely many points. In this mission it is nonempty and full-dimensional in its ambient real coordinate space. Thus a d-polytope P in R^d contains enough points to span that space affinely. A face is the set of points on which a supporting linear functional attains its maximum; a facet is a face of dimension d−1. Different supporting functionals can describe the same face, so the facet allowance counts geometric faces rather than inequalities used to describe them.

A k-simplex is the convex hull of k+1 affinely independent points. Write its vertices as v₀,…,vₖ in R^k and its convex hull as T. Affine independence means there is no nontrivial affine relation among these vertices. The simplex therefore has dimension k, and its interior is taken in the whole space R^k. No restriction is placed on its side lengths or angles. A d-flat is a translate of a d-dimensional linear subspace. A section of T by such a flat L is the entire intersection T ∩ L.

The target is the following existence assertion. For every d-polytope P with at most k+1 facets, every k-simplex T in R^k, and every p in the interior of T, there is a d-flat L such that

p∈L,T∩L is affinely equivalent to P.p\in L,\qquad T\cap L\text{ is affinely equivalent to }P.p∈L,T∩L is affinely equivalent to P.

An affine equivalence here preserves affine combinations and is invertible on the affine spans. It can change lengths and angles. In the stated coordinates, it is represented by an injective real affine map A from R^d onto L satisfying

A(P)=T∩L.A(P)=T\cap L.A(P)=T∩L.

All input data, including the interior point p, are universally quantified before the flat and map are chosen. The flat can depend on these data. The source writes the facet allowance as f and the simplex dimension as f−1; the notation here uses f=k+1.

The prescribed-point conclusion gives control beyond the existence of some simplex section representing P. It permits the simplex and an interior point to be fixed while the cutting flat is selected to recover the given polytope. The conclusion concerns its full affine geometry: the image of every point of P lies in the section, and every point of the section belongs to that image. This is the additional representational constraint discussed immediately before Theorem 10 on printed page 74.

The difficulty lies in meeting these requirements simultaneously. A flat through the chosen point can have an intersection of the wrong affine shape. A section with the correct affine shape need not contain the prescribed point. The theorem requires both properties for every allowed simplex and point, while retaining the original facet bound. Neither a special regular simplex nor a selected interior point captures this quantifier structure.

The formal statement uses real coordinate spaces indexed by finite sets, finite convex hulls, affine spans, exposed faces, and affine maps. Nonempty faces are counted by their affine dimension. In dimension zero the facet allowance is automatic: under the convention counting the empty face it contributes at most one facet, which fits the allowance k+1. This avoids interpreting natural subtraction d−1 as a facet dimension at d=0. The zero-dimensional simplex and zero-dimensional polytope remain within the statement.

This is a known geometric theorem whose formal proof remains to be supplied. The target contains one proof placeholder. The definitions specify the underlying geometry concretely. The mission consists of the single prescribed-section theorem; its definitions support that statement, and no separate supporting theorem is included as a milestone. The source is Grünbaum, Convex Polytopes, second edition, Springer, 2003, §5.1, Theorem 10, printed page 74 (source.pdf page 100); see also §5.1, printed page 71 (PDF page 97), for the section convention.

3 thms1 active userReviewed
Discrete Geometry·Captain: mikedeng1

Grünbaum — Gale data and combinatorial type (5.4.5)Textbook

Gale data and combinatorial type

A convex polytope has geometric coordinates and a combinatorial structure: its faces, ordered by inclusion. Different coordinates can describe the same combinatorial type. Gale transforms connect these viewpoints by encoding affine dependencies among the vertices as another point configuration, often in a smaller-dimensional space.

Let P and Q be full-dimensional d-polytopes in real coordinate space, each with n vertices. Choose injective lists V and W containing all their vertices and fix a permutation θ of the n indices. The question is whether this prescribed vertex correspondence extends to an isomorphism between the full face posets of P and Q. The empty face and the whole polytope are included in these posets.

An affine dependency of V is a list of real coefficients a whose sum is zero and whose weighted sum of the vertices is zero. These coefficient lists form a vector space of dimension n−d−1. Choose any basis and put its vectors into the columns of a matrix. The n rows of this matrix form a Gale transform G of V. A Gale transform H of W is constructed in the same way, with its own choice of basis. The rows may repeat and may be zero; there is no general-position assumption.

Theorem 5.4.5 states that the prescribed correspondence extends to an isomorphism of face posets exactly when it preserves every relative-interior test on the Gale configurations. For every subset J of vertex indices, the origin belongs to the relative interior of the convex hull of the rows G(J) if and only if it belongs to the relative interior of the convex hull of H(θ(J)). Relative interior means interior within the affine span of the selected convex hull, not interior in the whole coordinate space.

The equivalence concerns every subset of indices and both directions of implication. It retains the prescribed permutation rather than merely asking whether the polytopes have some combinatorial equivalence. Indexing the Gale points also keeps track of which original vertices they represent when several Gale rows coincide. Taking a set image does not change the convex hull of a chosen subconfiguration.

The empty subset is included: its convex hull has empty relative interior, so both membership tests are false. For a simplex, n=d+1 and the Gale space has dimension zero. Its nonempty row subconfigurations consist of the zero vector, and the relative-interior formulation still applies. Zero-dimensional polytopes are also retained.

This mission is the equivalence criterion itself, with the concrete definitions of a full-dimensional polytope, its face poset, and a Gale transform. The grouping keeps transform construction within the vocabulary of the criterion. The neighboring results on affine and projective equivalence have different conclusions and are outside this goal.

Source: Branko Grünbaum, Convex Polytopes, second edition (2003), §5.4, Theorem 5, printed page 89 (PDF page 115). The affine-dependence construction appears on printed pages 85–86 (PDF pages 111–112).

4 thms1 active userReviewed
PreviousPage 11 of 12Next

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