Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Combinatorics

265 missions · 156 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

Open109Completed156All265
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
Graph Theory·Captain: hao jia

Monochromatic Reachability in Three-Colored Tournaments (OPG-1808)Open Problem

Motivation

Edge-colored tournaments combine a complete orientation with a finite palette. They are a natural setting for comparing local multicolor obstructions with global directed reachability. The question attributed to Sands, Sauer, and Woodrow asks whether three colors force one of two outcomes: a directed triangle whose three arcs all have different colors, or a single vertex that can reach every target along a monochromatic directed path.

The problem was recorded by the Open Problem Garden in 2008. A minimum-counterexample reduction was later restated by Georgakopoulos and Sprüssel in their study of three-colored tournaments. The available project computation excludes counterexamples through eleven vertices, but that package is explicitly candidate_only: it is bounded search evidence, not a proof of the unrestricted theorem.

Setting

A tournament is an orientation of a finite complete simple graph. For each pair of distinct vertices u,vu,vu,v, exactly one of u→vu\to vu→v and v→uv\to uv→u is present. Every directed arc receives one of three labeled colors.

A rainbow directed triangle is a cyclically oriented triangle

a→b→c→aa\to b\to c\to aa→b→c→a

whose three arc colors are pairwise distinct. A transitive three-vertex subtournament is not a directed triangle and is therefore not forbidden merely because its three arcs have different colors.

A vertex sss is a monochromatic source when, for every vertex ttt, there is some color kkk and a directed sss-to-ttt path all of whose arcs have color kkk. The chosen color may depend on ttt; the theorem does not demand one common color for all targets. Length-zero reachability handles t=st=st=s.

The formal domain is nonempty finite tournaments. This nonemptiness convention is stated explicitly because an empty vertex type has neither a rainbow triangle nor a candidate source and would trivialize the negation of the intended question.

Formalization targets

Root theorem

For every nonempty finite tournament TTT with a three-coloring of its arcs,

T has a rainbow directed triangle∨∃s∈V(T) ∀t∈V(T), s reaches t monochromatically.T\text{ has a rainbow directed triangle} \quad\lor\quad \exists s\in V(T)\ \forall t\in V(T),\ \text{$s$ reaches $t$ monochromatically}. T has a rainbow directed triangle∨∃s∈V(T) ∀t∈V(T), s reaches t monochromatically.

No compatibility is required between the colors of paths to different targets, and unused palette colors are permitted.

Finite order milestone

The first milestone freezes the exact bounded claim supported by the replay package:

1≤∣V(T)∣≤11 and no rainbow directed triangle⟹T has a monochromatic source. 1\le |V(T)|\le 11\text{ and no rainbow directed triangle} \quad\Longrightarrow\quad T\text{ has a monochromatic source}.1≤∣V(T)∣≤11 and no rainbow directed triangle⟹T has a monochromatic source.

The statement includes all tournaments and all three-color arc assignments at those orders, not only one symmetry representative. The repository's observations report exhaustive search after a minimum-counterexample reduction, but the Lean theorem remains open until it has an accepted proof.

Significance

The root theorem would turn a local forbidden configuration into a global reachability certificate. Such a result clarifies how orientation and edge color interact: ordinary Gallai decompositions for undirected colored complete graphs cannot be imported unchanged, because the hypothesis forbids only rainbow cyclic triangles and allows rainbow transitive triples.

The formal development creates reusable definitions for colored directed reachability and exposes the direction of every relation. This matters in minimum-counterexample arguments, where an auxiliary arc u→Fvu\to_F vu→F​v may encode that vvv cannot reach uuu; reversing that convention invalidates the cycle reduction. A verified finite milestone would also provide a regression target for SAT, SMT, or exhaustive encodings without elevating their raw output to a universal theorem.

Difficulty

The classical Gallai theorem is not directly applicable. It assumes an undirected complete graph with no rainbow triangle of any orientation, whereas this problem permits a transitive triple with three distinct colors. A proposed partition must therefore control both arc colors and directions between parts.

The minimum-counterexample route yields a useful spanning cycle in an auxiliary nonreachability digraph. It does not itself bound the size of a counterexample. The order-eleven computation terminates because its domain is finite, but no induction from eleven to arbitrary order follows. A proof must add a structural theorem that survives all orientations and allows monochromatic paths of arbitrary length rather than treating reachability bits as independent physical arcs.

Formalization scope

Lean represents the tournament as a binary relation D with looplessness and exactly one orientation on each unordered pair. The coloring is a total function on ordered pairs, but only values on actual arcs are semantically used. Monochromatic reachability is the reflexive transitive closure of arcs of one fixed color. The root and finite theorem quantify over every nonempty finite vertex type.

The finite replay, its solver versions, hashes, and no-witness observations remain external candidate evidence. They do not close the milestone without a checkable certificate or a proof accepted by the platform. Contributions may formalize the minimum-counterexample cycle lemma, build an independently checked finite certificate, isolate a directed decomposition theorem, or prove the root. No contribution may replace a directed rainbow triangle by an undirected one, require the same path color for every target, or assume heredity of failure for arbitrary induced subtournaments.

Selected references

  • Open Problem Garden, Monochromatic reachability versus rainbow triangles, posted 2008. https://www.openproblemgarden.org/op/monochromatic_reachability_vs_rainbow_triangles
  • B. Sands, N. Sauer, and R. Woodrow, On monochromatic paths in edge-coloured digraphs, Journal of Combinatorial Theory, Series B 33 (1982), 271–275.
  • A. Georgakopoulos and P. Sprüssel, On 3-coloured tournaments, 2009. https://arxiv.org/abs/0904.1967
  • A. Trygub, Full Characterization of Color Degree Sequences in Complete Graphs Without Tricolored Triangles, 2023. https://arxiv.org/abs/2304.14579
3 thms1 active userReviewed
Graph TheoryOperations Research·Captain: mikedeng1

On the Graph Structure of Convex Polyhedra in n-Space II: Whitney's Theorem, a Graph Is n-Tuply Connected iff Any Two Points Are Joined by n Disjoint PathsResearch Paper

Motivation

Vertex connectivity measures how robust a network is against the failure of nodes. It can be measured in two ways that look different. One way counts the fewest nodes whose removal disconnects the network. The other counts the routes between two nodes that share no intermediate node. Whitney's theorem (1932) says that the two measures agree for every pair of nodes. It is the vertex form of Menger's theorem, and it underlies reliability analysis of communication and transportation networks, the design of fault-tolerant routing, and much of structural graph theory.

M. L. Balinski's 1961 paper On the graph structure of convex polyhedra in n-space proves that the graph of a bounded full-dimensional polyhedron in nnn-space is nnn-tuply connected (the subject of Mission I of this series). It then invokes Whitney's theorem to conclude that any two vertices of such a polyhedron are joined by nnn disjoint paths. Balinski gives a short new proof of Whitney's theorem through the max-flow min-cut theorem of Ford and Fulkerson and of Dantzig and Fulkerson. That makes the theorem a consequence of linear programming duality. This mission formalizes that part of the paper: the network vocabulary, the max-flow min-cut theorem with capacities on both points and lines, the integrality of maximum flows, and Whitney's theorem itself.

Timeline.

  • 1927: Menger states the disjoint-paths theorem for separating sets.
  • 1932: Whitney proves the characterization of nnn-connected graphs by nnn disjoint paths between every pair of points.
  • 1956: Ford and Fulkerson and Dantzig and Fulkerson prove the max-flow min-cut theorem.
  • 1961: Balinski derives Whitney's theorem from it with a unit-capacity network.

Setting

A graph GGG consists of a finite set VVV of points and a set of lines, each line being a pair of distinct points. A path from psp_sps​ to pkp_kpk​ is a sequence of lines (p1,p2),(p2,p3),…,(pm,pm+1)(p_1,p_2),(p_2,p_3),\dots,(p_m,p_{m+1})(p1​,p2​),(p2​,p3​),…,(pm​,pm+1​) with p1=psp_1 = p_sp1​=ps​, pm+1=pkp_{m+1} = p_kpm+1​=pk​ and m≥1m \ge 1m≥1. Paths are disjoint if they have no point in common except possibly their first and last points.

GGG is nnn-tuply connected if it has at least n+1n+1n+1 points and, for every set XXX of fewer than nnn points, the graph G−XG - XG−X remaining after deleting XXX is connected. GGG has nnn disjoint paths from psp_sps​ to pkp_kpk​ if there are nnn pairwise distinct paths from psp_sps​ to pkp_kpk​, none of which repeats a point, and no two of which share a point other than psp_sps​ and pkp_kpk​.

A network is a connected graph with a capacity c(x)≥0c(x) \ge 0c(x)≥0 on every point and c(e)≥0c(e) \ge 0c(e)≥0 on every line, and with a distinguished source psp_sps​ and sink pkp_kpk​. A flow assigns a number f(C)≥0f(C) \ge 0f(C)≥0 to every path CCC from psp_sps​ to pkp_kpk​, such that for every point xxx and every line eee

∑C∋xf(C)≤c(x),∑C∋ef(C)≤c(e).\sum_{C \ni x} f(C) \le c(x), \qquad \sum_{C \ni e} f(C) \le c(e).C∋x∑​f(C)≤c(x),C∋e∑​f(C)≤c(e).

Its value is val⁡(f)=∑Cf(C)\operatorname{val}(f) = \sum_C f(C)val(f)=∑C​f(C). A disconnecting set is a pair (X,F)(X,F)(X,F) of points and lines that meets every walk from psp_sps​ to pkp_kpk​. Its value is ∑x∈Xc(x)+∑e∈Fc(e)\sum_{x\in X} c(x) + \sum_{e \in F} c(e)∑x∈X​c(x)+∑e∈F​c(e).

The unit network of the proof has capacity 111 on every point except psp_sps​ and pkp_kpk​, and capacity n+1n+1n+1 on every line except the line pspkp_sp_kps​pk​ (if present), which has capacity 111. In Lean these are IsNTuplyConnected, HasNDisjointPaths, IsFlow, flowValue, IsDisconnecting, cutValue, unitCapV and unitCapE, all in the namespace Balinski61.Whitney.

Formalization targets

Goal: Whitney's theorem (p. 434)

For a finite graph GGG with at least two points and any n≥0n \ge 0n≥0:

G is n-tuply connected  ⟺  for all ps≠pk, G has n disjoint paths from ps to pk.G \text{ is } n\text{-tuply connected} \iff \text{for all } p_s \ne p_k,\ G \text{ has } n \text{ disjoint paths from } p_s \text{ to } p_k.G is n-tuply connected⟺for all ps​=pk​, G has n disjoint paths from ps​ to pk​.

Both directions are part of the goal.

Milestones, in the order of the proof

  1. Max-flow min-cut (p. 433). In every network there is a number MMM that is the value of some flow and of some disconnecting set, with every flow of value at most MMM and every disconnecting set of value at least MMM.
  2. Integrality (p. 434). If all capacities are integers, some maximum flow has only integer path flows.
  3. Min-cut in the unit network (p. 434). If GGG is nnn-tuply connected and ps≠pkp_s \ne p_kps​=pk​, every disconnecting set of the unit network has value at least nnn.
  4. Paths from unit flows (p. 434). An integral flow of value at least nnn in the unit network yields nnn disjoint paths from psp_sps​ to pkp_kpk​.
  5. Sufficiency (p. 434). If every pair of distinct points is joined by nnn disjoint paths, GGG is nnn-tuply connected.

Significance

Whitney's theorem turns a statement about all small deletion sets into the existence of explicit, verifiable path systems, and back again. In applications it certifies connectivity by exhibiting paths, and it certifies that connectivity is no larger by exhibiting a separating set. It is the base of the theory of kkk-connected graphs: ear decompositions, the fan lemma, and the structure of minimally kkk-connected graphs all use it. Inside this paper it supplies the COROLLARY that any two vertices of a bounded full-dimensional polyhedron in nnn-space are joined by nnn disjoint edge paths.

None of these results is formalized for vertex connectivity at this Mathlib revision. Mathlib has edge connectivity and connected components but no vertex Menger theorem. The Prove2Me library has max-flow min-cut statements for arc capacities only and integrality results for basic solutions of network LPs, but no flow model with capacities on points. A completed development gives a reusable vertex-capacitated max-flow min-cut theorem for undirected graphs and the first machine-checked Whitney theorem in this library. The results are classical and proved; the remaining work is the formalization.

Difficulty

The sufficiency direction is elementary. The necessity direction needs a global object (a flow, or a family of paths) to exist from purely local hypotheses about deletions. The obvious induction on nnn, which deletes a point and applies the hypothesis to a smaller graph, does not keep the path systems disjoint. In Balinski's route the weight falls on max-flow min-cut and integrality for path flows with capacities on points, neither of which exists in the library. A second difficulty sits in a case the paper's proof skips: a disconnecting set of the unit network may use the line pspkp_sp_kps​pk​ (capacity 111) together with up to n−2n-2n−2 points, and the deletion hypothesis of nnn-tuple connectedness speaks only about points.

Formalization scope

Graphs are Mathlib SimpleGraphs on a Fintype with decidable equality and adjacency. Paths are walks with IsPath. A flow is a real function on the finite type G.Path ps pk of simple paths; "through a point" and "through a line" mean membership in the walk's support and edge list. Capacities are functions V → ℝ and Sym2 V → ℝ. Nonnegativity, connectivity of GGG and ps≠pkp_s \ne p_kps​=pk​ are hypotheses of the network theorems.

The following readings of loose phrases are explicit in the statements:

  • "dropping out n−1n-1n−1 or fewer points" is ∣X∣<n|X| < n∣X∣<n;
  • "nnn disjoint paths" means nnn pairwise distinct simple paths. The printed path syntax permits repeated vertices; in a graph with lines psap_s aps​a, psbp_s bps​b, and pspkp_s p_kps​pk​, the distinct walks ps,a,ps,pkp_s,a,p_s,p_kps​,a,ps​,pk​ and ps,b,ps,pkp_s,b,p_s,p_kps​,b,ps​,pk​ share only their endpoints even though the graph is not 222-tuply connected. The theorem therefore uses its conventional simple-path reading;
  • path flows live on simple paths (merging and shortcutting changes no maximum value);
  • the paper leaves the capacities of psp_sps​ and pkp_kpk​ in the unit network unassigned, and here they are n+1n+1n+1;
  • "the condition is sufficient is obvious" is the full statement that nnn-tuple connectedness follows;
  • the hypothesis ∣V∣≥2|V| \ge 2∣V∣≥2 is added to the goal and to sufficiency, because the paper's "any pair of points" presupposes it and the equivalence fails for a one-point graph.

The max-flow min-cut milestone states that the maximum and the minimum are attained. A statement that only bounds some flow by every cut is satisfied by the zero flow. A connectivity notion without the n+1n+1n+1 point count would make every complete graph nnn-connected for all nnn. Both trivializations are excluded.

Contributions welcome: a vertex-capacitated augmenting-path or LP-duality proof of max-flow min-cut for path flows, integrality by an augmenting-path argument, the unit-network lemmas, and direct combinatorial proofs of Whitney's theorem that bypass flows.

Selected references

  • M. L. Balinski, On the graph structure of convex polyhedra in n-space, Pacific J. Math. 11 (1961), 431–434. https://doi.org/10.2140/pjm.1961.11.431
  • H. Whitney, Congruent graphs and the connectivity of graphs, Amer. J. Math. 54 (1932), 150–168. https://doi.org/10.2307/2371086
  • L. R. Ford, Jr. and D. R. Fulkerson, Maximal flow through a network, Canadian J. Math. 8 (1956), 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • G. B. Dantzig and D. R. Fulkerson, On the max-flow min-cut theorem of networks, in Linear Inequalities and Related Systems, Ann. of Math. Stud. 38, Princeton Univ. Press, 1956, 215–221. https://doi.org/10.1515/9781400881987
  • K. Menger, Zur allgemeinen Kurventheorie, Fund. Math. 10 (1927), 96–115. https://doi.org/10.4064/fm-10-1-96-115
9 thms0 active usersReviewed
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 thms0 active usersReviewed
Machine LearningTheoretical Computer Science·Captain: mikedeng1

A Characterization of Multiclass Learnability 1: Classes of Finite DS Dimension Have n → r Sample Compression Schemes with r Polylogarithmic in nResearch Paper

Motivation

In multiclass classification a learner sees examples (x,y)(x, y)(x,y) with xxx in a domain X\mathcal XX and a label yyy in a set Y\mathcal YY, and must predict labels of new points. When Y\mathcal YY is finite, the Natarajan dimension characterizes PAC learnability, extending the role of the VC dimension in binary classification (Natarajan 1989; Ben-David, Cesa-Bianchi, Haussler, Long 1995). Label sets in practice are often unbounded: structured prediction, ranking, and language modelling all predict from very large or infinite label spaces. For infinite Y\mathcal YY the Natarajan dimension fails to characterize learnability, and the question of which combinatorial parameter does was left open by Daniely and Shalev-Shwartz.

Timeline:

  • 1989–1995. Natarajan, then Ben-David et al. and Haussler–Long: for finite Y\mathcal YY, learnability is equivalent to finite Natarajan dimension, with sample complexity depending on log⁡∣Y∣\log|\mathcal Y|log∣Y∣.
  • 2011–2015. Daniely, Sabato, Ben-David and Shalev-Shwartz show that ERM can fail for multiclass problems with many labels. Daniely and Shalev-Shwartz (COLT 2014) introduce the DS dimension, prove that finite DS dimension is necessary for learnability, and ask whether it is sufficient.
  • 2022. Brukhim, Carmon, Dinur, Moran, Yehudayoff prove sufficiency, so the DS dimension characterizes multiclass PAC learnability, and show that the Natarajan dimension does not.

Setting

A concept class is a set H⊆YX\mathcal H\subseteq\mathcal Y^{\mathcal X}H⊆YX of functions. For a sequence S=(x1,…,xn)S=(x_1,\dots,x_n)S=(x1​,…,xn​) the projection H∣S⊆Yn\mathcal H|_S\subseteq\mathcal Y^nH∣S​⊆Yn is the set of label words (h(x1),…,h(xn))(h(x_1),\dots,h(x_n))(h(x1​),…,h(xn​)), h∈Hh\in\mathcal Hh∈H. A finite non-empty set B⊆YdB\subseteq\mathcal Y^dB⊆Yd is a pseudo-cube if every h∈Bh\in Bh∈B has, in every coordinate iii, a neighbour g∈Bg\in Bg∈B that differs from hhh exactly in coordinate iii. The sequence 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 maximum length of a DS-shattered sequence. The Natarajan dimension dN(H)≤dDS(H)d_N(\mathcal H)\le d_{DS}(\mathcal H)dN​(H)≤dDS​(H) is the same with Boolean cubes ∏i{f(i),g(i)}\prod_i\{f(i),g(i)\}∏i​{f(i),g(i)}, f(i)≠g(i)f(i)\ne g(i)f(i)=g(i), in place of pseudo-cubes.

A sample S∈(X×Y)nS\in(\mathcal X\times\mathcal Y)^nS∈(X×Y)n is H\mathcal HH-realizable if some h∈Hh\in\mathcal Hh∈H is consistent with it. An n→rn\to rn→r sample compression scheme for H\mathcal HH (Littlestone and Warmuth 1986) is a single reconstruction function ρ:(X×Y)r→YX\rho:(\mathcal X\times\mathcal Y)^r\to\mathcal Y^{\mathcal X}ρ:(X×Y)r→YX such that every realizable sample of size nnn contains rrr of its examples S′S'S′ with ρ(S′)\rho(S')ρ(S′) consistent with the whole sample. Logarithms are base 222 throughout.

Formalization targets

Goal: Theorem 36 (p. 22)

For H\mathcal HH with dDS(H)=dDS<∞d_{DS}(\mathcal H)=d_{DS}<\inftydDS​(H)=dDS​<∞ and dN(H)=dNd_N(\mathcal H)=d_NdN​(H)=dN​, and all integers n,t>0n,t>0n,t>0, there is an n→rn\to rn→r sample compression scheme, r≤nr\le nr≤n, with

r≤(dDS+t+1t+1(dDS+t)+103dNlog⁡((dDS+t+1t+1)log⁡(2n)))log⁡(2n).r\le\left(\frac{d_{DS}+t+1}{t+1}(d_{DS}+t)+10^3d_N\log\left(\binom{d_{DS}+t+1}{t+1}\log(2n)\right)\right)\log(2n).r≤(t+1dDS​+t+1​(dDS​+t)+103dN​log((t+1dDS​+t+1​)log(2n)))log(2n).

Milestones

The scheme combines two components, each with its own chain of results:

  • List learning from the DS dimension. Lemma 13 (orientations of out-degree ≤d\le d≤d on Yd+1\mathcal Y^{d+1}Yd+1), Claim 16 (the one-inclusion algorithm is right on some leave-one-out example), Fact 14 (leave-one-out symmetrization), Proposition 32 (a list PAC learner with list size (d+tt)\binom{d+t}{t}(td+t​) and success probability t+1d+t+1\frac{t+1}{d+t+1}d+t+1t+1​), Lemma 39 (an n→r1n\to r_1n→r1​ list compression scheme with r1≤dDS+t+1t+1(dDS+t)log⁡(2n)r_1\le\frac{d_{DS}+t+1}{t+1}(d_{DS}+t)\log(2n)r1​≤t+1dDS​+t+1​(dDS​+t)log(2n) and menu size ≤(dDS+t+1t+1)log⁡(2n)\le\binom{d_{DS}+t+1}{t+1}\log(2n)≤(t+1dDS​+t+1​)log(2n)).
  • Learning from a menu via shifting. Claim 22, Corollary 23, Claim 26, Proposition 27 (avd⁡≤4dE\operatorname{avd}\le4d_Eavd≤4dE​), Corollary 28, Lemma 29 (dE≤5dNlog⁡pd_E\le5d_N\log pdE​≤5dN​logp), Lemma 17 (orientations of out-degree ≤20dNlog⁡p\le20d_N\log p≤20dN​logp on [p]n[p]^n[p]n), Proposition 34 (error ≤20dNlog⁡(p)/n\le20d_N\log(p)/n≤20dN​log(p)/n given a ppp-menu), Lemma 40 (an n→r2n\to r_2n→r2​ compression scheme given a ppp-menu with r2≤103dNlog⁡(p)log⁡(2n)r_2\le10^3d_N\log(p)\log(2n)r2​≤103dN​log(p)log(2n)).

Significance

Theorem 36 is the algorithmic heart of the characterization: by the standard "compression implies generalization" argument it gives PAC learnability of every class of finite DS dimension, with sample complexity O~(dDS3/2/ϵ)\tilde O(d_{DS}^{3/2}/\epsilon)O~(dDS3/2​/ϵ) in the realizable case (t=⌈dDS1/2⌉t=\lceil d_{DS}^{1/2}\rceilt=⌈dDS1/2​⌉), and with the agnostic case following by known reductions. It also exhibits sample compression schemes of size polylogarithmic in nnn for multiclass classes with infinitely many labels, in contrast to the constant-size schemes known for finite VC classes.

The result is proved in the paper; none of it is formalized. The formalization would produce a machine-checked theory of one-inclusion graphs and their orientations, multiclass shifting, the exponential dimension, list learning, and sample compression schemes for arbitrary label sets. These objects recur throughout learning theory (one-inclusion graphs in optimal PAC learning, shifting in VC theory), so the infrastructure is reusable beyond this mission.

Difficulty

The natural first idea, running empirical risk minimization or bounding the Natarajan dimension, fails: classes with Natarajan dimension 111 and infinitely many labels can be unlearnable, and ERM can fail even for learnable classes. The DS dimension gives only a weak guarantee: by Claim 16, among d+1d+1d+1 leave-one-out runs, one is correct. Turning this into a learner requires a list learner whose menus are still of unbounded total size, and then learning with a menu of size ppp, where the obstacle is controlling one-inclusion graph orientations over [p]n[p]^n[p]n by the Natarajan dimension. Multiclass shifting does not preserve the average degree (Example 20), so the binary argument breaks down, and a new potential (avd⁡′\operatorname{avd}'avd′) and a new dimension (dEd_EdE​) are needed. Lemma 13 for infinite classes needs a compactness argument.

Formalization scope

Lean conventions:

  • A class is H : Set (X → Y) with arbitrary types X, Y; sequences and samples are functions on Fin n ([n][n][n] is 0-based).
  • The DS, Natarajan and exponential dimensions are suprema in ℕ∞, so unbounded families give ⊤; hypotheses are written dsDim H = dDS with dDS : ℕ. A pseudo-cube is required to be finite.
  • Logarithms are Real.logb 2. Menu sizes use Set.encard.
  • A compression scheme is a reconstruction function fixed before the sample (∃ ρ, ∀ S, ∃ S'); a subsample may repeat and reorder examples. Theorem 36 states r≤nr\le nr≤n explicitly.
  • Orientations of the one-inclusion graph of V⊆YmV\subseteq\mathcal Y^mV⊆Ym are maps sending a direction iii and a vertex vvv to the head of the edge of direction iii through vvv; the out-degree of vvv counts directions whose head is not vvv.
  • Classes over [p][p][p] use labels Fin p; the shifting condition 1≤g(i)≤∣ef∣1\le g(i)\le|e_f|1≤g(i)≤∣ef​∣ becomes g(i)<∣ef∣g(i)<|e_f|g(i)<∣ef​∣.
  • The one-inclusion algorithm (Algorithms 1 and 3) is parametrized by a permutation-equivariant choice of minimal orientations, the reading under which the paper's leave-one-out proofs are valid; statements about the algorithm hold for every such choice. Its default output on non-realizable input requires a non-empty label set, assumed in Claim 16 and Propositions 32 and 34. Lemma 40 assumes a non-empty label set because it is false for X≠∅=Y\mathcal X\ne\emptyset=\mathcal YX=∅=Y.
  • Distributions are discrete (PMF), and i.i.d. probabilities are sums over Zm\mathcal Z^mZm. The measure-theoretic generality of the paper is not attempted.

A trivial formalization is ruled out by these choices. Placing the reconstruction function after the sample would let it output the consistent hypothesis. A dimension in ℕ defined by sSup would be 000 for infinite dimension. Pseudo-cubes without finiteness would change the dimension (Example 8).

Contributions welcome: proofs of any milestone, in particular the shifting results of §3 (self-contained combinatorics on finite classes), Fact 14 (pure discrete probability), and Lemma 13; general-purpose lemmas about one-inclusion graphs, orientations and sample compression schemes are reusable by other learning-theory missions.

Selected references

  • N. Brukhim, D. Carmon, I. Dinur, S. Moran, A. Yehudayoff, A Characterization of Multiclass Learnability, FOCS 2022; arXiv:2203.01550v1 (2022). https://arxiv.org/abs/2203.01550
  • A. Daniely, S. Shalev-Shwartz, Optimal Learners for Multiclass Problems, COLT 2014. https://arxiv.org/abs/1405.2690
  • N. Littlestone, M. Warmuth, Relating Data Compression and Learnability, unpublished technical report, University of California, Santa Cruz, 1986 (no stable link).
  • D. Haussler, N. Littlestone, M. Warmuth, Predicting {0,1}-Functions on Randomly Drawn Points, Information and Computation 115(2), 1994. https://doi.org/10.1006/inco.1994.1097
  • D. Haussler, P. M. Long, A Generalization of Sauer's Lemma, Journal of Combinatorial Theory, Series A 71(2), 1995. https://doi.org/10.1016/0097-3165(95)90006-3
  • S. Ben-David, N. Cesa-Bianchi, D. Haussler, P. M. Long, Characterizations of Learnability for Classes of {0,…,n}-Valued Functions, JCSS 50(1), 1995. https://doi.org/10.1006/jcss.1995.1008
  • B. K. Natarajan, On Learning Sets and Functions, Machine Learning 4, 1989. https://doi.org/10.1007/BF00114804
20 thms0 active usersReviewed
Operations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 6: The Optimal Value of S_G Lies Between n² − an²(ln 1/a + 2) and n² − an²Research Paper

Motivation

The problem 1 ∣ prec ∣ ∑wjCj1\,|\,\mathrm{prec}\,|\,\sum w_jC_j1∣prec∣∑wj​Cj​ asks for a single-machine sequence of jobs, respecting precedence constraints, that minimizes the weighted sum of completion times. It has been known to be strongly NP-hard since Lawler (1978) and Lenstra and Rinnooy Kan (1978), several different 2-approximation algorithms are known, and closing the approximability gap is listed by Schuurman and Woeginger (1999) as one of ten outstanding open problems in scheduling theory. Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 36(4), 2011) give the first inapproximability result for this problem: under a widely believed complexity assumption it has no polynomial-time approximation scheme (PTAS). The bridge to that result is a quantitative link, Lemma 9.1, between the optimal value of a special bipartite scheduling instance and the maximum edge biclique of a bipartite graph, a problem whose hardness of approximation was established by Ambühl, Mastrolilli and Svensson (FOCS 2007). This mission formalizes that link.

Setting

A schedule of a finite job set is a sequence σ\sigmaσ listing every job once; the machine processes the jobs in that order from time 000 without idle time or pre-emption. Job jjj has a processing time pjp_jpj​ and a weight wjw_jwj​; its completion time CjC_jCj​ is the sum of the processing times of the jobs up to and including jjj, and the value of σ\sigmaσ is val(σ)=∑jwjCj\mathrm{val}(\sigma)=\sum_j w_jC_jval(σ)=∑j​wj​Cj​. Precedence constraints are a relation PPP on jobs: (i,j)∈P(i,j)\in P(i,j)∈P with i≠ji\ne ji=j means job iii must be completed before job jjj starts. A schedule respecting all of them is feasible, and a feasible schedule σ∗\sigma^*σ∗ of least value is optimal.

Let G=(U,V,E)G=(U,V,E)G=(U,V,E) be an nnn-by-nnn bipartite graph: ∣U∣=∣V∣=n|U|=|V|=n∣U∣=∣V∣=n and E⊆U×VE\subseteq U\times VE⊆U×V. An edge biclique is a pair A⊆UA\subseteq UA⊆U, B⊆VB\subseteq VB⊆V with A×B⊆EA\times B\subseteq EA×B⊆E, of value ∣A∣⋅∣B∣|A|\cdot|B|∣A∣⋅∣B∣; the maximum edge biclique problem (Definition 9.1) asks for one of largest value. The bipartite scheduling instance SGS_GSG​ has jobs U∪VU\cup VU∪V and precedence constraints

P=(U×V)∖E,P=(U\times V)\setminus E,P=(U×V)∖E,

so u∈Uu\in Uu∈U must precede v∈Vv\in Vv∈V exactly when (u,v)(u,v)(u,v) is not an edge. Jobs of UUU have p=1p=1p=1, w=0w=0w=0; jobs of VVV have p=0p=0p=0, w=1w=1w=1. Thus val(σ)=∑v∈VCv\mathrm{val}(\sigma)=\sum_{v\in V}C_vval(σ)=∑v∈V​Cv​, where CvC_vCv​ is the number of UUU-jobs scheduled before vvv. For i≥1i\ge1i≥1, σ(i)\sigma(i)σ(i) denotes the number of VVV-jobs scheduled before iii jobs of UUU have been scheduled.

In the Lean development these are weightedCompletion, IsOptimalSchedule, IsEdgeBiclique, maxBicliqueValue, precSG, procSG, weightSG, valSG, IsOptimalSG and vBefore in the namespace SingleMachinePrec.Biclique.

Formalization targets

Goal: Lemma 9.1 (p. 666)

If a maximum edge biclique of GGG has value an2an^2an2 with a∈(0,1]a\in(0,1]a∈(0,1], then SGS_GSG​ has an optimal schedule and every optimal schedule σ∗\sigma^*σ∗ satisfies

n2−an2(ln⁡1a+2)≤val(σ∗)≤n2−an2.n^2-an^2\Bigl(\ln\frac1a+2\Bigr)\le\mathrm{val}(\sigma^*)\le n^2-an^2 .n2−an2(lna1​+2)≤val(σ∗)≤n2−an2.

Milestones (proof of Lemma 9.1, §9, p. 666)

  1. For every edge biclique (A,B)(A,B)(A,B), a schedule in the block order U∖A→B→A→V∖BU\setminus A\to B\to A\to V\setminus BU∖A→B→A→V∖B exists, and every such schedule is feasible with
val(σ)=n2−∣A∣⋅∣B∣.\mathrm{val}(\sigma)=n^2-|A|\cdot|B| .val(σ)=n2−∣A∣⋅∣B∣.
  1. For every schedule, σ(n+1)=n\sigma(n+1)=nσ(n+1)=n and
val(σ)=∑i=1n(σ(i+1)−σ(i))i=n2−∑i=1nσ(i).\mathrm{val}(\sigma)=\sum_{i=1}^n\bigl(\sigma(i+1)-\sigma(i)\bigr)i=n^2-\sum_{i=1}^n\sigma(i).val(σ)=i=1∑n​(σ(i+1)−σ(i))i=n2−i=1∑n​σ(i).
  1. For every feasible schedule and i=1,…,ni=1,\dots,ni=1,…,n,
σ(i)(n−i+1)≤an2,σ(i)≤n.\sigma(i)(n-i+1)\le an^2,\qquad \sigma(i)\le n .σ(i)(n−i+1)≤an2,σ(i)≤n.

Significance

Lemma 9.1 shows that the optimal value of SGS_GSG​ determines the maximum edge biclique of GGG up to a factor of order ln⁡(1/a)\ln(1/a)ln(1/a) in the "area above the work line" n2−val(σ∗)n^2-\mathrm{val}(\sigma^*)n2−val(σ∗). Combined with the hardness of approximating maximum edge biclique (Theorem 9.1, cited from Ambühl, Mastrolilli and Svensson 2007) it yields Theorem 9.2: 1 ∣ prec ∣ ∑wjCj1\,|\,\mathrm{prec}\,|\,\sum w_jC_j1∣prec∣∑wj​Cj​ has no PTAS unless SAT can be decided by a probabilistic algorithm in time 2Nϵ2^{N^\epsilon}2Nϵ for every ϵ>0\epsilon>0ϵ>0. It also makes precise the two-dimensional Gantt chart picture of Eastman, Even and Isaacs (1964) and of Goemans and Williamson (2000), in which every point on the work line of a schedule defines an edge biclique.

The lemma is proved in the paper; it is not formalized anywhere to our knowledge. A formal proof certifies the combinatorial core of the no-PTAS result independently of the complexity-theoretic layer, and its definitions (the bipartite instance SGS_GSG​, edge bicliques, the profile σ(i)\sigma(i)σ(i)) are reusable for the gap inequality behind Theorem 9.2.

Difficulty

The upper bound is a direct computation on one explicit schedule. The lower bound is a statement about every feasible schedule, of which there are exponentially many, and it must hold with the explicit constant 222 and the factor ln⁡(1/a)\ln(1/a)ln(1/a) for every a∈(0,1]a\in(0,1]a∈(0,1]. The printed argument splits the sum at i=(1−a)ni=(1-a)ni=(1−a)n and uses ⌊an⌋\lfloor an\rfloor⌊an⌋, treating ananan as an integer; for general aaa (for example n=3n=3n=3, value 222, an=2/3an=2/3an=2/3) the split point is not an integer, so the printed estimate does not apply verbatim and the constant 222 has to be re-checked for non-integral ananan. On the formal side, the value identity requires relating completion times in a list to counting UUU-jobs before each VVV-job, with ties among zero-length jobs.

Formalization scope

Jobs are the disjoint union U ⊕ V of two finite types with Fintype.card U = Fintype.card V = n; EEE is a relation U → V → Prop. A schedule is a duplicate-free list containing every job; feasibility is the published LawlerPrec.MinMax.IsFeasible and completion times are the published MooreLateJobs.Shared.completionTime (time 000 start, no idle time). Processing times and weights are reals, here in {0,1}\{0,1\}{0,1}. The maximum edge biclique value is the maximum of ∣A∣⋅∣B∣|A|\cdot|B|∣A∣⋅∣B∣ over all edge bicliques, the empty ones included, so the hypothesis a>0a>0a>0 means E≠∅E\ne\emptysetE=∅. The logarithm is natural (Real.log).

Conventions and readings committed to:

  • The goal is stated for every optimal schedule, and the existence of an optimal schedule is a separate conclusion, so the bounds cannot hold vacuously. Proving the bounds for one particular schedule, or for an optimal value defined as an infimum that could be a junk default, would not be this lemma.
  • No integrality hypothesis on ananan is added.
  • Milestones 2 and 3 are stated for every schedule (respectively every feasible schedule), not only for σ∗\sigma^*σ∗; milestone 1 states the value of the block-order schedule as an equality, where the paper writes "≤⋯=\le\cdots=≤⋯=".
  • The paper's P=(U×V)∖EP=(U\times V)\setminus EP=(U×V)∖E is irreflexive; feasibility only constrains distinct jobs, so it agrees with the reflexive partial order of §1.

Not formalized: Theorem 9.1 (cited hardness of maximum edge biclique) and Theorem 9.2 (no PTAS under a complexity assumption); no polynomial-time or complexity-theoretic statement appears in the mission. Contributions welcome: proofs of the three milestones and of the goal; Mathlib's bounds on harmonic numbers (Mathlib/NumberTheory/Harmonic/Bounds.lean) are the relevant library.

Selected references

  • C. Ambühl, M. Mastrolilli, N. Mutsanas, O. Svensson, On the Approximability of Single-Machine Scheduling with Precedence Constraints, Mathematics of Operations Research 36(4):653–669, 2011. https://doi.org/10.1287/moor.1110.0512
  • C. Ambühl, M. Mastrolilli, O. Svensson, Inapproximability results for sparsest cut, optimal linear arrangement, and precedence constraint scheduling, Proc. 48th IEEE FOCS, 329–337, 2007 (reference [4] of the paper).
  • W. L. Eastman, S. Even, I. M. Isaacs, Bounds for the optimal scheduling of n jobs on m processors, Management Science 11(2):268–279, 1964 (reference [11]).
  • M. X. Goemans, D. P. Williamson, Two-dimensional Gantt charts and a scheduling algorithm of Lawler, SIAM J. Discrete Math. 13(3):281–294, 2000 (reference [15]).
  • P. Schuurman, G. J. Woeginger, Polynomial time approximation algorithms for machine scheduling: ten open problems, J. Scheduling 2(5):203–213, 1999 (reference [36]).
8 thms0 active usersReviewed
Graph TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 1: A k-Fold Realizer of Size t Yields a Vertex Cover of Expected Weight at Most (2 − 2/(t/k)) Times OptimalResearch Paper

Motivation

Single-machine scheduling with precedence constraints, written 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ in the notation of Graham et al., asks for an order of nnn weighted jobs on one machine that respects a given partial order and minimizes the weighted sum of completion times. The problem is strongly NP-hard (Lawler 1978; Lenstra and Rinnooy Kan 1978), and closing its approximability gap is listed by Schuurman and Woeginger among ten outstanding open problems in scheduling theory. Several 2-approximation algorithms are known (Schulz 1996; Hall et al. 1997; Chudak and Hochbaum 1999; Chekuri and Motwani 1999; Margot et al. 2003).

A line of work by Chudak and Hochbaum, Correa and Schulz, and Ambühl and Mastrolilli showed that the problem is a special case of minimum weighted vertex cover in a graph built from the precedence order. Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 36(4), 2011) observed that this graph is the graph of incomparable pairs of dimension theory, and used that identification to obtain (2−2/f)(2-2/f)(2−2/f)-approximations for orders of fractional dimension at most fff. This mission formalizes that framework: the identification of the two graphs and the rounding guarantee of the paper's Theorem 5.1.

Setting

An instance SSS consists of a finite set NNN of jobs, a partial order PPP on NNN (reflexive, antisymmetric, transitive; (i,j)∈P(i,j)\in P(i,j)∈P with i≠ji\ne ji=j means job iii finishes before job jjj starts), processing times pj≥0p_j\ge 0pj​≥0 and weights wj≥0w_j\ge 0wj​≥0.

Two jobs x,yx,yx,y are incomparable, x∥yx\parallel yx∥y, when neither (x,y)(x,y)(x,y) nor (y,x)(y,x)(y,x) lies in PPP. The set inc⁡(P)\operatorname{inc}(P)inc(P) of incomparable pairs consists of ordered pairs and is closed under swapping. A linear extension of PPP is a linear order L⊇PL\supseteq PL⊇P on NNN; it reverses (x,y)∈inc⁡(P)(x,y)\in\operatorname{inc}(P)(x,y)∈inc(P) when y<xy<xy<x in LLL. A nonempty multiset L1,…,LtL_1,\dots,L_tL1​,…,Lt​ of linear extensions is a k:tk:tk:t-realizer if every incomparable pair is reversed by at least kkk of them. The fractional dimension fdim⁡(P)\operatorname{fdim}(P)fdim(P) is the least ratio t/kt/kt/k over all k:tk:tk:t-realizers.

The vertex cover graph GPSG^S_PGPS​ has the incomparable pairs as nodes; nodes (i,j)(i,j)(i,j) and (k,ℓ)(k,\ell)(k,ℓ) are adjacent if j=kj=kj=k and i=ℓi=\elli=ℓ, or j=kj=kj=k and (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P, or (i,ℓ),(k,j)∈P(i,\ell),(k,j)\in P(i,ℓ),(k,j)∈P (symmetrically closed). Node (i,j)(i,j)(i,j) has weight w(i,j)=piwjw_{(i,j)}=p_iw_jw(i,j)​=pi​wj​, and w(C)=∑u∈Cwuw(C)=\sum_{u\in C}w_uw(C)=∑u∈C​wu​. OPT\mathrm{OPT}OPT is the minimum weight of a vertex cover of GPSG^S_PGPS​. The LP relaxation [CS-LP] asks for x∈[0,1]inc⁡(P)x\in[0,1]^{\operatorname{inc}(P)}x∈[0,1]inc(P) with xu+xv≥1x_u+x_v\ge1xu​+xv​≥1 on every edge, minimizing ∑uwuxu\sum_u w_ux_u∑u​wu​xu​. For a solution xxx write Va={u:xu=a}V_a=\{u: x_u=a\}Va​={u:xu​=a}, and for a linear extension LLL let I1/2(L)I_{1/2}(L)I1/2​(L) be the pairs of V1/2V_{1/2}V1/2​ reversed in LLL.

The graph of incomparable pairs GPG_PGP​ (Felsner and Trotter 2000) also has the incomparable pairs as vertices; two of them are adjacent when the pair of them is a minimal set of incomparable pairs that no linear extension reverses entirely.

Formalization targets

Goal: Theorem 5.1 (p. 658)

For an instance SSS, a k:tk:tk:t-realizer L1,…,LtL_1,\dots,L_tL1​,…,Lt​ of PPP, and a half-integral optimal solution xxx of [CS-LP], put Ci=V1∪(V1/2∖I1/2(Li))C_i=V_1\cup\bigl(V_{1/2}\setminus I_{1/2}(L_i)\bigr)Ci​=V1​∪(V1/2​∖I1/2​(Li​)). Then every CiC_iCi​ is a vertex cover of GPSG^S_PGPS​, and

1t∑i=1tw(Ci)  ≤  (2−2t/k)OPT.\frac1t\sum_{i=1}^t w(C_i)\;\le\;\Bigl(2-\frac{2}{t/k}\Bigr)\mathrm{OPT}.t1​i=1∑t​w(Ci​)≤(2−t/k2​)OPT.

Milestones

  • Proposition 3.2 (p. 657): GPS=GPG^S_P=G_PGPS​=GP​.
  • Footnote 4 (p. 659): the pairs reversed by a linear extension are independent in GPSG^S_PGPS​.
  • Eq. (4): 1t∣{i:Li reverses u}∣≥k/t\frac1t|\{i: L_i\text{ reverses }u\}|\ge k/tt1​∣{i:Li​ reverses u}∣≥k/t for every incomparable pair uuu.
  • Eq. (5): 1t∑iw(I1/2(Li))≥kt w(V1/2)\frac1t\sum_i w(I_{1/2}(L_i))\ge \frac kt\,w(V_{1/2})t1​∑i​w(I1/2​(Li​))≥tk​w(V1/2​).
  • Hochbaum's observation (§5, p. 659): for half-integral feasible xxx, V1∪CV_1\cup CV1​∪C covers GPSG^S_PGPS​ whenever CCC covers GPS[V1/2]G^S_P[V_{1/2}]GPS​[V1/2​].
  • Eqs. (6)–(8): 1t∑iw(Ci)≤w(V1)+(1−kt)w(V1/2)≤2(1−kt)(w(V1)+12w(V1/2))≤(2−2t/k)OPT\frac1t\sum_iw(C_i)\le w(V_1)+(1-\frac kt)w(V_{1/2})\le 2(1-\frac kt)(w(V_1)+\frac12w(V_{1/2}))\le(2-\frac2{t/k})\mathrm{OPT}t1​∑i​w(Ci​)≤w(V1​)+(1−tk​)w(V1/2​)≤2(1−tk​)(w(V1​)+21​w(V1/2​))≤(2−t/k2​)OPT when PPP is not a linear order.

Significance

The result. Combined with the cited Theorem 2.1 (Ambühl–Mastrolilli 2009; Correa–Schulz 2005), which turns an α\alphaα-approximate vertex cover of GPSG^S_PGPS​ into an α\alphaα-approximate schedule, Theorem 5.1 gives a (2−2/f)(2-2/f)(2−2/f)-approximation for 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ whenever the precedence order has an efficiently samplable realizer with t/k≤ft/k\le ft/k≤f. The paper applies it to interval orders (3/23/23/2), convex bipartite orders and semiorders (4/34/34/3), orders of bounded degree and orders of interval dimension two; for the earlier special classes it matched or improved the best known ratios, and for the last two it gave the first results. Proposition 3.2 makes the dimension theory of posets (realizers, critical pairs, fractional dimension) directly available to the vertex cover approach.

Formalizing it. The results are proved in the paper; to the best of current knowledge none of them has a machine-checked proof. The mission produces a Lean development of incomparable pairs, linear extensions, kkk-fold realizers and the hypergraph of incomparable pairs, which other dimension-theory missions can reuse, and a verified LP-rounding argument for half-integral vertex cover solutions under a distribution of independent sets.

Difficulty

Most of the rounding argument is arithmetic over finite sums. The central nontrivial step is the inclusion GP⊆GPSG_P\subseteq G^S_PGP​⊆GPS​ in Proposition 3.2: for two incomparable pairs that are not adjacent under the three-case rule, one must construct a single linear extension reversing both. This needs an extension of PPP by two new comparabilities whose transitive closure is still antisymmetric, followed by Szpilrajn's theorem; checking that every potential cycle is excluded by the three cases is the actual content. The opposite inclusion, and footnote 4, follow from transitivity and antisymmetry of linear orders. A second point of care is the inequality t≥2kt\ge 2kt≥2k used in step (7): it is not part of the definition of a realizer, and follows from each linear extension reversing exactly one of (x,y)(x,y)(x,y) and (y,x)(y,x)(y,x).

Formalization scope

  • Jobs form a finite type N; the precedence order is a relation P : N → N → Prop with IsPartialOrder, carried by the structure Instance. Processing times and weights are nonnegative reals.
  • inc⁡(P)\operatorname{inc}(P)inc(P) is the subtype IncPair P of N × N; a linear extension is a relation with IsLinearOrder containing P; a k:tk:tk:t-realizer is a family Fin t → LinearExtension P with t>0t>0t>0. Reversal of (x,y)(x,y)(x,y) means y<xy<xy<x in LLL throughout; the page's "y>xy>xy>x" in Eq. (4) and "Prob[j>i]\mathrm{Prob}[j>i]Prob[j>i]" in Eq. (5) are the same family of inequalities because inc⁡(P)\operatorname{inc}(P)inc(P) is symmetric.
  • GPSG^S_PGPS​ is the symmetric closure of the printed three-case rule on distinct nodes. GPG_PGP​ is defined through linear extensions and hyperedge minimality, never through the three-case rule, so Proposition 3.2 is a genuine statement and not a definitional equality.
  • [CS-LP] drops the constant term ∑jpjwj+∑(i,j)∈Ppiwj\sum_jp_jw_j+\sum_{(i,j)\in P}p_iw_j∑j​pj​wj​+∑(i,j)∈P​pi​wj​ of [CS-IP], which does not affect optimality. OPT\mathrm{OPT}OPT is the minimum weight of a vertex cover of GPSG^S_PGPS​, taken over a finite nonempty family.
  • Not formalized: "efficiently samplable", "polynomial time" and "randomized algorithm". The expectation over a uniformly sampled LiL_iLi​ is stated as the average 1t∑i=1t\frac1t\sum_{i=1}^tt1​∑i=1t​, which is equivalent to it and stronger than the existence of one good index. The existence of a half-integral optimal [CS-LP] solution (Nemhauser–Trotter, cited) is a hypothesis on xxx. The conversion of a vertex cover into a schedule (Theorem 2.1, cited) is not formalized; the goal is stated for vertex covers of GPSG^S_PGPS​.
  • Constants: 2−2/(t/k)2-2/(t/k)2−2/(t/k) in real arithmetic; it equals 2−2k/t2-2k/t2−2k/t, and equals 222 when k=0k=0k=0.
  • The paper assumes fdim⁡(P)≥2\operatorname{fdim}(P)\ge2fdim(P)≥2, i.e. PPP is not a linear order. Eqs. (6)–(8) carry that hypothesis as the paper does; the goal omits it because for a linear order both sides are 000.
  • Conclusion (a), that each CiC_iCi​ is a vertex cover, is part of the goal and is not assumed.

Contributions welcome: proofs of the milestones in any order, and a reusable Szpilrajn-style lemma producing a linear extension that reverses a prescribed set of compatible incomparable pairs.

Selected references

  • C. Ambühl, M. Mastrolilli, N. Mutsanas, O. Svensson, On the approximability of single-machine scheduling with precedence constraints, Math. Oper. Res. 36(4):653–669, 2011. https://doi.org/10.1287/moor.1110.0512
  • C. Ambühl, M. Mastrolilli, Single machine precedence constrained scheduling is a vertex cover problem, Algorithmica 53(4), 2009 (reference [2] of the paper).
  • J. R. Correa, A. S. Schulz, Single-machine scheduling with precedence constraints, Math. Oper. Res. 30(4):1005–1021, 2005. https://doi.org/10.1287/moor.1050.0158
  • G. R. Brightwell, E. R. Scheinerman, Fractional dimension of partial orders, Order 9(2):139–158, 1992 (reference [7]).
  • S. Felsner, W. T. Trotter, Dimension, graph and hypergraph coloring, Order 17(2):167–177, 2000 (reference [13]).
  • D. S. Hochbaum, Efficient bounds for the stable set, vertex cover and set packing problems, Discrete Appl. Math. 6(3):243–254, 1983 (reference [20]).
  • G. L. Nemhauser, L. E. Trotter, Vertex packings: structural properties and algorithms, Math. Programming 8(1):232–248, 1975 (reference [29]).
10 thms0 active usersReviewed
Graph TheoryOperations ResearchTheoretical Computer Science·Captain: mikedeng1

On the Approximability of Single-Machine Scheduling with Precedence Constraints 4: Vertex Cover in Connected Graphs of Degree ≤ 3 Reduces to Weighted Vertex Cover for Interval-Order InstancesResearch Paper

Motivation

In the single-machine scheduling problem 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​, a set NNN of nnn jobs, each with a processing time pj≥0p_j\ge 0pj​≥0 and a weight wj≥0w_j\ge 0wj​≥0, is processed on one machine without interruption, subject to precedence constraints given by a partial order PPP on NNN. The aim is to minimize the weighted sum of completion times ∑jwjCj\sum_j w_jC_j∑j​wj​Cj​. The problem is strongly NP-hard for general precedence constraints (Lawler 1978; Lenstra and Rinnooy Kan 1978), and its approximability was a recurring open question in scheduling theory (Schuurman and Woeginger 1999).

A line of work by Chudak and Hochbaum, Correa and Schulz, and Ambühl and Mastrolilli showed that the problem is a special case of minimum weighted vertex cover in a graph GPSG^S_PGPS​ built from the instance. Many problems on partial orders become polynomial when the order is an interval order, so it is natural to ask whether this one does too. Section 7 of Ambühl, Mastrolilli, Mutsanas and Svensson (Math. Oper. Res. 2011) answers no: the problem stays NP-hard on interval orders. The proof is a reduction from vertex cover in connected graphs of maximum degree 3. This mission formalizes the correctness of that reduction.

Setting

A poset P=(N,P)P=(N,P)P=(N,P) is read as a reflexive relation: (x,y)∈P(x,y)\in P(x,y)∈P means x≤yx\le yx≤y. Jobs x,yx,yx,y are incomparable if neither (x,y)(x,y)(x,y) nor (y,x)(y,x)(y,x) is in PPP, and inc⁡(P)\operatorname{inc}(P)inc(P) is the set of ordered incomparable pairs. The vertex cover graph GPSG^S_PGPS​ has one node (i,j)(i,j)(i,j) for each incomparable pair, weighted piwjp_iw_jpi​wj​. Two nodes (i,j)(i,j)(i,j) and (k,ℓ)(k,\ell)(k,ℓ) are adjacent if j=kj=kj=k and i=ℓi=\elli=ℓ, or j=kj=kj=k and (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P, or (i,ℓ)∈P(i,\ell)\in P(i,ℓ)∈P and (k,j)∈P(k,j)\in P(k,j)∈P. Write w(CI)w(C_I)w(CI​) for the minimum weight of a vertex cover of GISG^S_IGIS​.

A poset is an interval order if each element xxx can be assigned a closed real interval [ax,bx][a_x,b_x][ax​,bx​] such that x<yx<yx<y if and only if bx<ayb_x<a_ybx​<ay​.

The reduction starts from a graph G=(V,E)G=(V,E)G=(V,E) with vertices v1,…,vNv_1,\dots,v_Nv1​,…,vN​ and a spanning tree T=(V,ET)T=(V,E_T)T=(V,ET​) rooted at v1v_1v1​, numbered so that each parent comes before its children. The paper uses a breadth-first search tree.

  • Stage 1. The graph G′G'G′ is built from TTT. Each viv_ivi​ gets a pendant path vi−u1i−u2iv_i - u^i_1 - u^i_2vi​−u1i​−u2i​. Each non-tree edge {vi,vj}∈E∖ET\{v_i,v_j\}\in E\setminus E_T{vi​,vj​}∈E∖ET​ with i<ji<ji<j gets the path vi−e1ij−e2ij−u2jv_i - e^{ij}_1 - e^{ij}_2 - u^j_2vi​−e1ij​−e2ij​−u2j​. The non-tree edges themselves are not edges of G′G'G′.
  • Stage 2. The scheduling instance SSS has jobs s0s_0s0​, s1,…,sNs_1,\dots,s_Ns1​,…,sN​, m1,…,mNm_1,\dots,m_Nm1​,…,mN​, e1,…,eNe_1,\dots,e_Ne1​,…,eN​, and bijb_{ij}bij​ for each non-tree edge. Their intervals, processing times and weights are given in a table on p. 662. For example, sjs_jsj​ has interval [i,j][i,j][i,j], processing time 1/kj1/k^j1/kj and weight kik^iki, where viv_ivi​ is the parent of vjv_jvj​. The precedence constraints III are the interval order of these intervals. With nnn the number of jobs, the parameter is k=n2+1k=n^2+1k=n2+1.
  • The set DDD. It is {(s0,s1)}∪{(si,sj):vi parent of vj}∪{(si,mi),(mi,ei)}∪{(si,bij),(bij,mj)}\{(s_0,s_1)\}\cup\{(s_i,s_j): v_i \text{ parent of } v_j\}\cup\{(s_i,m_i),(m_i,e_i)\}\cup\{(s_i,b_{ij}),(b_{ij},m_j)\}{(s0​,s1​)}∪{(si​,sj​):vi​ parent of vj​}∪{(si​,mi​),(mi​,ei​)}∪{(si​,bij​),(bij​,mj​)}. The graph GI′G'_IGI′​ is the subgraph of GISG^S_IGIS​ induced by DDD.

Formalization targets

Goal: Theorem 7.1 (p. 661)

For every connected graph GGG of maximum degree at most 333, every parent-first spanning tree TTT and every m∈Nm\in\mathbb Nm∈N, the precedence constraints III of SSS form an interval order, and

G has a vertex cover of size≤m  ⟺  ⌊w(CI)⌋≤m+∣V∣+∣E∖ET∣.G \text{ has a vertex cover of size} \le m \iff \lfloor w(C_I)\rfloor \le m + |V| + |E\setminus E_T|.G has a vertex cover of size≤m⟺⌊w(CI​)⌋≤m+∣V∣+∣E∖ET​∣.

Milestones, in the order the proof uses them

  • Claim 1 (p. 662): τ(G′)=τ(G)+∣V∣+∣E∖ET∣\tau(G') = \tau(G)+|V|+|E\setminus E_T|τ(G′)=τ(G)+∣V∣+∣E∖ET​∣, where τ\tauτ is the vertex cover number.
  • Remark 7.1 (p. 662): for jobs with intervals [a,b][a,b][a,b] and [c,d][c,d][c,d] and a≤da\le da≤d, pi≤1/k⌈b⌉p_i\le 1/k^{\lceil b\rceil}pi​≤1/k⌈b⌉ and wj≤k⌈c⌉w_j\le k^{\lceil c\rceil}wj​≤k⌈c⌉. On incomparable pairs piwj∈{1}∪[0,1/k]p_iw_j\in\{1\}\cup[0,1/k]pi​wj​∈{1}∪[0,1/k]. Moreover, piwj≥kp_iw_j\ge kpi​wj​≥k forces b<cb<cb<c, and piwj=1p_iw_j=1pi​wj​=1 forces ⌈b⌉=⌈c⌉\lceil b\rceil=\lceil c\rceil⌈b⌉=⌈c⌉.
  • Claim 2 (p. 663): an incomparable pair (i,j)(i,j)(i,j) has piwj=1p_iw_j=1pi​wj​=1 if it is in DDD, and piwj≤1/kp_iw_j\le 1/kpi​wj​≤1/k otherwise.
  • Claim 3 (p. 663): GI′≅G′G'_I\cong G'GI′​≅G′.
  • §7, p. 664: with k=n2+1k=n^2+1k=n2+1, ∑(i,j)∈inc⁡(I)∖Dpiwj<1\sum_{(i,j)\in\operatorname{inc}(I)\setminus D}p_iw_j<1∑(i,j)∈inc(I)∖D​pi​wj​<1, and hence w(CI′)=⌊w(CI)⌋w(C'_I)=\lfloor w(C_I)\rfloorw(CI′​)=⌊w(CI​)⌋.

Significance

The result. Interval orders are a standard tractable class: several scheduling and order-theoretic problems that are hard in general become polynomial on them (Papadimitriou and Yannakakis 1979). Theorem 7.1 puts 1∣prec∣∑wjCj1|\mathrm{prec}|\sum w_jC_j1∣prec∣∑wj​Cj​ outside this pattern. Section 6 of the same paper shows that the problem nonetheless has a 3/23/23/2-approximation on interval orders, so hardness and approximability are separated on this class. The paper also remarks that the proof makes weighted vertex cover NP-hard to approximate within some factor r>1r>1r>1 on the graphs GISG^S_IGIS​ arising from interval orders.

Formalizing it. The theorem is proved in the paper. Nothing in this mission is open, and none of it has been machine-checked before. The work splits into the following parts:

  • a gadget argument on unweighted vertex covers (Claim 1, after Alimonti and Kann);
  • an exact case analysis of incomparable pairs in a concrete interval order (Remark 7.1, Claim 2);
  • a graph isomorphism (Claim 3);
  • a rounding argument that links weighted and unweighted optima.

The definitions of GPSG^S_PGPS​ and of minimum-weight vertex covers are shared with the other missions of this series.

Difficulty

The construction is explicit, and each step is elementary. The work is in the bookkeeping. Claim 2 requires classifying every incomparable pair of jobs, including pairs of different kinds such as (bij,sℓ)(b_{ij}, s_\ell)(bij​,sℓ​), by comparing ceilings of interval endpoints. Half-integer endpoints (mim_imi​, bijb_{ij}bij​) are exactly what separates weight-one pairs from comparable ones. Claim 3 requires checking adjacency in GISG^S_IGIS​ for all pairs of nodes of DDD in both directions. The paper writes out two cases in each direction and calls the rest similar.

Claim 1 has a direction that is not simply local. A vertex cover of G′G'G′ that misses both endpoints of a non-tree edge has to be repaired by swapping gadget vertices, and the repair must be repeated without increasing the size.

A natural first idea is to treat the light nodes (weight at most 1/k1/k1/k) as negligible one at a time. This does not suffice: the argument needs their total weight to stay below 111, which is what forces kkk to grow with n2n^2n2.

Formalization scope

  • Graph and tree. GGG is a SimpleGraph (Fin N); vertex vi+1v_{i+1}vi+1​ is i, and the root is index 0. The tree is a TreeLayout: a parent function returning none exactly at the root, with each parent of smaller index and adjacent in GGG. The statements hold for every such layout. This is stronger than the paper's breadth-first tree, and the proof uses only "parent before child".
  • Hypotheses of the goal. Connectivity and the degree bound ((G.neighborSet v).ncard ≤ 3) are kept as in the paper. They matter only for the NP-completeness of the source problem.
  • Jobs. The jobs form an inductive type with one constructor per row of the table. Their order is a PartialOrder instance: x≤yx\le yx≤y iff x=yx=yx=y or bx<ayb_x<a_ybx​<ay​. Processing times and weights are real numbers. Section 1 of the paper asks for nonnegative integers, but the instance uses 1/kj1/k^j1/kj and the formalization follows the instance as printed.
  • Constants. The constants are explicit: k=n2+1k=n^2+1k=n2+1 with nnn the cardinality of the job type, and c=∣V∣+∣E∖ET∣c=|V|+|E\setminus E_T|c=∣V∣+∣E∖ET​∣. Remark 7.1 and Claim 2 are stated for every real k>1k>1k>1.
  • Optimum values. w(CI)w(C_I)w(CI​) is a minimum over the finite family of vertex covers. Unweighted cover numbers are Mathlib's SimpleGraph.vertexCoverNum. The floor is Nat.floor, which agrees with the integer floor because w(CI)≥0w(C_I)\ge 0w(CI​)≥0.
  • Not formalized. The goal's wording ("NP-hard") is not formalized. Neither are the NP-completeness of degree-3 vertex cover (Garey, Johnson and Stockmeyer), the polynomial size of the construction, or Theorem 2.1 (cited), which turns a vertex cover of GISG^S_IGIS​ into a schedule. What is stated is the correctness of the reduction: the instance has interval-order constraints, and its optimum decides the vertex cover question.
  • No trivialization. The instance SSS is built from GGG and TTT exactly as in the table. The goal quantifies over all graphs and layouts, never over an instance SSS assumed to have the properties.
  • Contributions. Contributions are welcome on any milestone. Claims 1 and 3 are independent of the weights, and Claim 2 is independent of the graph theory.

Selected references

  • C. Ambühl, M. Mastrolilli, N. Mutsanas, O. Svensson, On the Approximability of Single-Machine Scheduling with Precedence Constraints, Math. Oper. Res. 36(4):653–669, 2011. https://doi.org/10.1287/moor.1110.0512
  • C. Ambühl, M. Mastrolilli, Single machine precedence constrained scheduling is a vertex cover problem, Algorithmica 53(4):488–503, 2009. https://doi.org/10.1007/s00453-008-9251-1
  • J. R. Correa, A. S. Schulz, Single machine scheduling with precedence constraints, Math. Oper. Res. 30(4):1005–1021, 2005. https://doi.org/10.1287/moor.1050.0158
  • P. Alimonti, V. Kann, Some APX-completeness results for cubic graphs, Theoret. Comput. Sci. 237(1–2):123–134, 2000. https://doi.org/10.1016/S0304-3975(98)00158-3
  • M. R. Garey, D. S. Johnson, L. Stockmeyer, Some simplified NP-complete graph problems, Theoret. Comput. Sci. 1(3):237–267, 1976. https://doi.org/10.1016/0304-3975(76)90059-1
  • C. H. Papadimitriou, M. Yannakakis, Scheduling interval-ordered tasks, SIAM J. Comput. 8(3):405–409, 1979. https://doi.org/10.1137/0208031
10 thms0 active usersReviewed
Linear algebraOperations Research·Captain: mikedeng1

On the Abstract Properties of Linear Dependence 4: Orthogonal Subspaces Have Dual MatroidsResearch Paper

Motivation

Hassler Whitney introduced matroids in On the Abstract Properties of Linear Dependence (Amer. J. Math. 57, 1935, doi:10.2307/2371182) to capture what linear dependence among the columns of a matrix and the circuit structure of a graph have in common. One of the central constructions of the paper is duality. For graphs, duality exists only for planar graphs (Whitney, Non-separable and planar graphs, Trans. AMS 34, 1932); for matroids Whitney showed that every matroid has a dual, and that duality has a direct linear-algebra model: the matroid of a subspace of Rn\mathbb R^nRn and the matroid of its orthogonal complement are duals.

Matroid duality underlies much of combinatorial optimization: the duality between cycles and cuts in graphs, the relation between a linear code and its dual code, the exchange of rank and corank in matroid intersection and partition, and the treatment of network flows as duals of potential problems. Whitney's §§11–13 are where this structure is first defined and where its linear-algebra meaning (Theorem 28) is established.

Setting

A matroid MMM is a finite set of elements e1,…,ene_1,\dots,e_ne1​,…,en​ with a rank function rrr on its subsets; equivalently, a family of independent sets, or of bases (maximal independent sets). For a subset NNN write ρ(N)\rho(N)ρ(N) for its number of elements and

n(N)=ρ(N)−r(N)n(N) = \rho(N) - r(N)n(N)=ρ(N)−r(N)

for its nullity; r(M)r(M)r(M) and n(M)n(M)n(M) are the rank and nullity of the whole set of elements.

Duality (§11). Let MMM and M′M'M′ be matroids and σ\sigmaσ a one-to-one correspondence between their elements. M′M'M′ is a dual of MMM (via σ\sigmaσ) if, for every subset NNN of MMM, with N′N'N′ the complement in M′M'M′ of σ(N)\sigma(N)σ(N),

r(N′)=r(M′)−n(N).(11.1)r(N') = r(M') - n(N). \tag{11.1}r(N′)=r(M′)−n(N).(11.1)

The matroid of a subspace (§12). Let EnE_nEn​ be nnn-dimensional Euclidean space with coordinates x1,…,xnx_1,\dots,x_nx1​,…,xn​, and let HHH be a hyperplane through the origin, which in Whitney's usage is a linear subspace of any dimension. For a set SSS of coordinates, project HHH onto the coordinate subspace ES′E'_SES′​ spanned by the axes xix_ixi​, i∈Si\in Si∈S. The matroid associated with HHH has elements e1,…,ene_1,\dots,e_ne1​,…,en​, one per coordinate, and gives the set {ei:i∈S}\{e_i: i\in S\}{ei​:i∈S} the rank

rH(S)=dim⁡πS(H).r_H(S) = \dim \pi_S(H).rH​(S)=dimπS​(H).

If HHH is the row space of a matrix M\mathbf MM, then rH(S)r_H(S)rH​(S) is the rank of the columns of M\mathbf MM indexed by SSS, so the associated matroid is the column matroid of M\mathbf MM.

In the Lean development these objects are WhitneyMatroid.Components.nullity (shared), IsDualVia M M' σ, IsDual M M', subspaceRank H S and IsAssociated M H.

Formalization targets

Goal: Theorem 28

Let HHH be a subspace of EnE_nEn​ and H′=H⊥H' = H^\perpH′=H⊥ its orthogonal complement. If MMM and M′M'M′ are the matroids associated with HHH and H′H'H′, then MMM and M′M'M′ are duals under the correspondence of equal coordinates:

∀N⊆{e1,…,en}:rM′(N‾)=r(M′)−nM(N).\forall N\subseteq\{e_1,\dots,e_n\}:\qquad r_{M'}(\overline N) = r(M') - n_M(N).∀N⊆{e1​,…,en​}:rM′​(N)=r(M′)−nM​(N).

Milestones

  1. Theorem 27. For every subspace HHH there is exactly one matroid associated with HHH.
  2. Theorem 6 (already proved on the platform). All bases have the same number of elements.
  3. Theorem 7. BBB is a base iff r(B)=r(M)r(B)=r(M)r(B)=r(M) and n(B)=0n(B)=0n(B)=0.
  4. Theorem 8. If BBB is a base and NNN is independent, then N∪N′N\cup N'N∪N′ is a base for some N′⊆BN'\subseteq BN′⊆B.
  5. Theorem 20. If M′M'M′ is a dual of MMM, then r(M′)=n(M)r(M') = n(M)r(M′)=n(M) and n(M′)=r(M)n(M') = r(M)n(M′)=r(M).
  6. Theorem 23. M′M'M′ is a dual of MMM via σ\sigmaσ iff, for every BBB, BBB is a base of MMM exactly when the complement of σ(B)\sigma(B)σ(B) is a base of M′M'M′.
  7. Theorem 21. Duality is symmetric.
  8. Theorem 22. Every matroid has a dual.

Significance

The result. Theorem 28 identifies abstract duality with orthogonal complementation. With Theorem 23 it says that the column matroid of a real matrix whose rows span HHH and the column matroid of a matrix whose rows span H⊥H^\perpH⊥ have complementary bases. This is the basis of the standard representation of the dual of a represented matroid ([Ir∣A][I_r\mid A][Ir​∣A] and [−AT∣In−r][-A^{\mathsf T}\mid I_{n-r}][−AT∣In−r​]), of the cycle/cocycle duality of graphs viewed through incidence matrices, and of the fact that a matroid representable over a field has a dual representable over the same field. Theorems 20–23 are the basic facts every later treatment of matroid duality starts from: duality exchanges rank and nullity, is symmetric, always exists, and is characterized by base complements.

Formalizing it. All of these results are classical and proved in Whitney's paper; what this mission adds is a machine-checked development of Whitney's own definition of duality, the rank identity (11.1), rather than the base-complement definition used in modern libraries. Mathlib defines the dual matroid M∗M^*M∗ by base complements and proves that it is a matroid; it does not, at the pinned revision, contain the rank formula (11.1) for M∗M^*M∗, the matroid of a real subspace, or Theorem 28. Theorem 6 is already proved on the platform and enters as a reference. Theorems 7 and 8 are close to Mathlib lemmas and are footholds rather than new content.

Difficulty

The difficulty of the goal is not in the combinatorics but in the passage between the two descriptions of a subspace. The natural first idea is to compare the ranks rH(S)r_H(S)rH​(S) and rH⊥(S‾)r_{H^\perp}(\overline S)rH⊥​(S) directly by counting dimensions of HHH and H⊥H^\perpH⊥. This does not close: rH(S)r_H(S)rH​(S) is the dimension of a projection of HHH onto coordinates, while dim⁡H⊥=n−dim⁡H\dim H^\perp = n - \dim HdimH⊥=n−dimH only controls H⊥H^\perpH⊥ as a whole, and the rank of a complementary coordinate set in M′M'M′ is not determined by the dimensions of HHH and H⊥H^\perpH⊥ alone. Some relation between coordinate projections of HHH and the structure of H⊥H^\perpH⊥ on the complementary coordinates has to be established. Theorem 27 is itself nontrivial in Lean: the rank function rHr_HrH​ has to be shown to satisfy the matroid axioms, which needs a careful treatment of coordinate projections and of submodularity of dimension. On the abstract side, Theorem 23 needs both directions of the passage between the rank identity and base complements, including the counting step r(M)+r(M′)=ρ(M)r(M) + r(M') = \rho(M)r(M)+r(M′)=ρ(M).

Formalization scope

  • Matroids are Mathlib Matroid α on a finite type α whose ground set is all of α (Whitney's matroid is its finite set of elements). This finiteness and the ground set convention are part of the definitions; all theorems are stated under them.
  • Ranks and nullities. Ranks are Mathlib's eRk, finite on a finite type and converted to integers; the nullity n(N)=ρ(N)−r(N)n(N)=\rho(N)-r(N)n(N)=ρ(N)−r(N) and the identity (11.1) are computed in Z\mathbb ZZ. No truncated subtraction appears.
  • Duality is Whitney's rank identity (11.1) for every subset, for a fixed bijection σ : α ≃ β (IsDualVia) or some bijection (IsDual). It is deliberately not defined as Mathlib's M✶: with that definition Theorem 23 would be nearly definitional and Theorem 28 would lose the rank identity. A formalization that replaces (11.1) by the base-complement condition, or states Theorem 28 only for bases, is not this mission's goal.
  • Euclidean space is EuclideanSpace ℝ (Fin n) with its standard inner product; the orthogonal hyperplane is Hᗮ. The field is R\mathbb RR, as in the paper; the analogous statement over other fields with the dot-product annihilator is a generalization and not part of this mission.
  • The associated matroid is a predicate (IsAssociated M H): ground set Fin n, and the rank of every set SSS of coordinates equals the finrank of the image of HHH under the coordinate projection onto SSS (not of H∩ES′H\cap E'_SH∩ES′​, a different set function). Existence and uniqueness is Theorem 27.
  • Implicit hypotheses made explicit: the ground sets are finite and equal to the whole element type; Whitney's "dimension rrr" and "dimension n−rn-rn−r" in Theorem 28 are consequences of H′=H⊥H'=H^\perpH′=H⊥ and are not hypotheses. The goal covers n=0n=0n=0, H={0}H=\{0\}H={0} and H=EnH=E_nH=En​.
  • Infrastructure needed and reusable: the matroid of a real subspace (equivalently the column matroid of a real matrix), with its rank function in terms of coordinate projections; the rank formula of the dual matroid; the dimension identity relating projections of HHH and sections of H⊥H^\perpH⊥. All are reusable for the other missions of this series (the Fano matroid and binary matroids) and for any later work on represented matroids. Proofs of the milestones, alternative proofs of Theorem 28 through Mathlib's M✶, and supporting lemmas on coordinate projections are welcome.

Selected references

  • H. Whitney, On the Abstract Properties of Linear Dependence, American Journal of Mathematics 57 (1935), 509–533. https://doi.org/10.2307/2371182
  • 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
  • J. Oxley, Matroid Theory, 2nd ed., Oxford Graduate Texts in Mathematics 21, Oxford University Press, 2011, Chapter 2 (duality). https://doi.org/10.1093/acprof:oso/9780198566946.001.0001
  • Mathlib, Mathlib.Combinatorics.Matroid (Dual, Rank). https://github.com/leanprover-community/mathlib4
12 thms0 active usersReviewed
Discrete GeometryLinear Optimization·Captain: mikedeng1

A Counterexample to the Hirsch Conjecture I: A 43-Dimensional Polytope with 86 Facets Whose Diameter Exceeds 43Research Paper

Motivation

The Hirsch conjecture, stated by Warren M. Hirsch in a 1957 letter to George Dantzig, asserts that the combinatorial diameter of a ddd-dimensional polytope with nnn facets is at most n−dn-dn−d. The diameter is the largest number of edge steps needed to walk between two vertices of the polytope. The question matters to linear optimization: every edge-following pivot rule of the simplex method, starting at a vertex of the feasible region, needs at least as many pivots as the graph distance to an optimal vertex. A polynomial bound on polytope diameters is a necessary condition for a strongly polynomial simplex method, which is still unknown.

Timeline of the conjecture before the counterexample:

  • 1967, Klee and Walkup: the conjecture fails for unbounded polyhedra, and for bounded polytopes it is equivalent to the ddd-step conjecture (the case n=2dn=2dn=2d). They also proved it for n−d≤5n-d\le 5n−d≤5.
  • 1992, Kalai and Kleitman: a quasi-polynomial upper bound nlog⁡2d+2n^{\log_2 d+2}nlog2​d+2 on the diameter.
  • 2010, Santos: a counterexample in dimension 43 with 86 facets (arXiv:1006.2814, published in Ann. of Math. 176 (2012) 383–412). It is the subject of this mission.
  • 2012, Matschke, Santos and Weibel: smaller 5-dimensional spindles, giving counterexamples in dimension 20 (arXiv:1202.4701).

The polynomial Hirsch conjecture, which asks only for a bound polynomial in nnn and ddd, remains open.

Setting

A polytope in Rd\mathbb R^dRd is written as Hpoly(a,b)={x:⟨ai,x⟩≤bi, i=1,…,n}\mathrm{Hpoly}(a,b)=\{x:\langle a_i,x\rangle\le b_i,\ i=1,\dots,n\}Hpoly(a,b)={x:⟨ai​,x⟩≤bi​, i=1,…,n}. A facet presentation is one in which Hpoly(a,b)\mathrm{Hpoly}(a,b)Hpoly(a,b) has nonempty interior and no inequality is redundant. For a bounded polytope it then has dimension ddd, and its nnn inequalities correspond one to one to its facets. A step between two vertices crosses an edge (a one-dimensional face). The polytope is non-Hirsch when its diameter exceeds n−dn-dn−d.

A spindle (Santos, Definition 1.4) is a polytope with two distinguished vertices u,vu,vu,v such that every facet contains exactly one of them. Its length is the graph distance from uuu to vvv. The polar objects are prismatoids: polytopes with two parallel facets Q+,Q−Q^+,Q^-Q+,Q− containing all vertices, whose width is the dual-graph distance between Q+Q^+Q+ and Q−Q^-Q−.

Santos's explicit object is the 5-prismatoid QQQ whose vertices are the 48 rows of Table 1 of the paper: 1+,…,24+1^+,\dots,24^+1+,…,24+ with x5=1x_5=1x5​=1 and 1−,…,24−1^-,\dots,24^-1−,…,24− with x5=−1x_5=-1x5​=−1. Its polar QΔ={x∈R5:⟨r,x⟩≤1 for each row r}Q^\Delta=\{x\in\mathbb R^5:\langle r,x\rangle\le1\ \text{for each row } r\}QΔ={x∈R5:⟨r,x⟩≤1 for each row r} is a spindle with apices e5e_5e5​ and −e5-e_5−e5​. Table 2 of the paper lists the 322 facets of QQQ, which are the 322 vertices of QΔQ^\DeltaQΔ: ±e5\pm e_5±e5​ and the points 1c0(±c1,±c2,±c3,±c4,c5)\frac1{c_0}(\pm c_1,\pm c_2,\pm c_3,\pm c_4,c_5)c0​1​(±c1​,±c2​,±c3​,±c4​,c5​) of twenty types B,B′,…,K,K′B,B',\dots,K,K'B,B′,…,K,K′ with sixteen sign patterns each. The symmetry group Σ\SigmaΣ of QQQ has order 64; its index-two subgroup Σ+\Sigma^+Σ+ preserves Q+Q^+Q+ and Q−Q^-Q−.

Formalization targets

Goal: Corollary 1.7

∃ a:Fin 86→R43, b∈R86:Hpoly(a,b) nonempty, bounded, a facet presentation, and diam⁡Hpoly(a,b)>43.\exists\, a:\mathrm{Fin}\,86\to\mathbb R^{43},\ b\in\mathbb R^{86}:\quad \mathrm{Hpoly}(a,b)\ \text{nonempty, bounded, a facet presentation, and}\ \operatorname{diam}\mathrm{Hpoly}(a,b)>43.∃a:Fin86→R43, b∈R86:Hpoly(a,b) nonempty, bounded, a facet presentation, and diamHpoly(a,b)>43.

The goal fixes the dimension (43) and the number of facets (86), as the paper does. It does not assert the exact diameter (44, remarked on p. 24).

Milestones

  1. §3, first bullet (polar form): QΔQ^\DeltaQΔ is a bounded 48-facet spindle with apices ±e5\pm e_5±e5​, the rows tight at e5e_5e5​ being 1+,…,24+1^+,\dots,24^+1+,…,24+.
  2. Theorem 4.1(1) (polar form): the 322 points of Table 2 are distinct vertices of QΔQ^\DeltaQΔ; the letter classes are Σ+\Sigma^+Σ+-orbits and the six Σ\SigmaΣ-orbits are A ∪ L, B ∪ K, C ∪ J, D ∪ I, E ∪ H, F ∪ G.
  3. Theorem 4.1(2) (polar form): QΔQ^\DeltaQΔ has no other vertices.
  4. Theorem 3.1 (polar form): the distance from e5e_5e5​ to −e5-e_5−e5​ in the graph of QΔQ^\DeltaQΔ is exactly 6.
  5. Theorem 1.6: a 5-spindle with 48 facets, 322 vertices and length exactly six exists.
  6. Inductive step of Theorem 2.6 (polar form): a ddd-spindle with n>2dn>2dn>2d facets and length at least lll yields a (d+1)(d+1)(d+1)-spindle with n+1n+1n+1 facets and length at least l+1l+1l+1.
  7. Theorem 1.5 (strong ddd-step theorem): a ddd-spindle with nnn facets and length at least lll yields an (n−d)(n-d)(n−d)-spindle with 2n−2d2n-2d2n−2d facets and length at least l+n−2dl+n-2dl+n−2d, non-Hirsch when l>dl>dl>d.

Significance

Corollary 1.7 settles the Hirsch conjecture in the negative. Theorem 1.5 is the reusable part: it turns any spindle whose length exceeds its dimension into a counterexample. That made the later search for small spindles (Matschke–Santos–Weibel) a finite, low-dimensional problem. Theorem 1.6 is the concrete certificate. Its claims about Tables 1 and 2 are finite statements about explicit rational data.

The result is proved, and several pieces are already formalized on Prove2Me. Hirsch.santos_counterexample (Proved) states that some polytope violates the bound n−dn-dn−d, with no dimension or facet count. Hirsch.strong_dstep_spindle (Proved) is the "in particular" clause of Theorem 1.5, counting inequalities rather than facets. Lemmas 2.2 and 2.4 of the paper are Proved in polar form (Hirsch.spindle_inward_row_push_graph, Hirsch.spindle_wedge_vertex_graph_projection), as is the bound n≥2dn\ge2dn≥2d for spindles (Hirsch.spindle_n_ge_two_d). This mission adds the dimension-43, 86-facet statement with facets counted exactly, the general length bound l+n−2dl+n-2dl+n−2d with the spindle structure preserved, and a machine-checkable account of Santos's specific spindle.

Difficulty

The obvious first idea for the 5-dimensional milestones is brute computation. This is feasible in principle but heavy in Lean. Showing that 322 given points are extreme requires exhibiting, for each, five linearly independent tight rows among 48. Showing that no other vertex exists requires a complete enumeration argument rather than a membership check. Exact distance 6 requires the full adjacency structure between the 322 vertices, not only a path.

For the general theorems, the difficulty is that the construction is a one-point suspension followed by a perturbation. Edges and facets of the new polytope must be controlled, with the face lattice changing only in a prescribed way. Diameter is not monotone under arbitrary perturbations, and in particular a walk can become shorter. Keeping the facet count exact (not merely the inequality count) adds an irredundancy argument at each step.

Formalization scope

Polytopes are H-presentations Hirsch.Hpoly a b with a : Fin n → EuclideanSpace ℝ (Fin d), from the platform definition Hirsch_model. "A ddd-polytope with nnn facets" is boundedness plus IsFacetPresentation: nonempty interior and every row irredundant. Without this predicate, padding with redundant rows would change the bound n−dn-dn−d, and the statement would collapse to the platform's dimension-free Hirsch.santos_counterexample. Diameter and length use the platform's padded walks (Hirsch.DiamLE, Hirsch.Reach). "Length at least LLL" means no walk of fewer than LLL steps, and "length exactly LLL" adds a walk of LLL steps. The junk-valued Hirsch.gdist is never used.

All prismatoid statements of the paper (§3, Theorems 2.6, 3.1, 4.1) are stated for the polar spindle, which the paper licenses (pp. 6, 8, 9). Vertices of QQQ become rows, facets of QQQ become vertices, the width becomes the apex distance, and Q±Q^\pmQ± become ±e5\pm e_5±e5​. Natural-number expressions n−dn-dn−d, 2n−2d2n-2d2n−2d, l+(n−2d)l+(n-2d)l+(n−2d) are faithful because spindles have n≥2dn\ge2dn≥2d. Theorem 1.5's "length lll" is read as "length at least lll", which is equivalent by monotonicity. Tables 1 and 2 are written entry by entry in Lean in the page's order.

Useful contributions include proofs of the finite Table 1/Table 2 facts by certified computation, reusable lemmas on one-point suspensions and on facet presentations under perturbation, and a proof of Corollary 1.7 from milestones 5 and 7.

Selected references

  • F. Santos, A counterexample to the Hirsch conjecture, Ann. of Math. 176 (2012) 383–412; cited version arXiv:1006.2814v3. https://arxiv.org/abs/1006.2814 , https://doi.org/10.4007/annals.2012.176.1.7
  • V. Klee and D. W. Walkup, The d-step conjecture for polyhedra of dimension d < 6, Acta Math. 117 (1967) 53–78. https://doi.org/10.1007/BF02395040
  • G. Kalai and D. J. Kleitman, A quasi-polynomial bound for the diameter of graphs of polyhedra, Bull. Amer. Math. Soc. 26 (1992) 315–316. https://arxiv.org/abs/math/9204233
  • B. Matschke, F. Santos and C. Weibel, The width of five-dimensional prismatoids, Proc. London Math. Soc. 110 (2015) 647–672. https://arxiv.org/abs/1202.4701
15 thms0 active usersReviewed
Operations ResearchOptimization·Captain: mikedeng1

Assortment Optimization under Variants of the Nested Logit Model 4: With Dissimilarity Parameters at Most One, the Knapsack-Relaxation and Singleton LP Optimum Scaled by 2 Is Feasible for the Full LPResearch Paper

Motivation

Assortment optimization asks which set of products a firm should offer when customers choose among the offered products according to a discrete choice model; it underlies shelf-space planning in retail and fare-class control in airline revenue management (Talluri and van Ryzin, 2004). Under the nested logit model products are grouped into nests, and a customer first picks a nest and then a product inside it. Davis, Gallego and Topaloglu (Operations Research, 2014; DOI 10.1287/opre.2014.1256) map out how hard this problem is across variants of the model.

When every nest dissimilarity parameter is at most one and a customer who chose a nest always buys there, offering the top-revenue products of each nest is optimal (Theorem 4 of the paper, the subject of an earlier mission of this series). Once a customer may leave a nest without buying — a partially-captured nest — that structure breaks and the problem becomes NP-hard (Theorem 8). This mission targets the paper's response: a small, explicitly constructed family of candidate assortments per nest from which a linear program recovers a solution within a factor of two of optimal.

Setting

There are mmm nests MMM and, in each nest, nnn products N={1,…,n}N = \{1, \dots, n\}N={1,…,n}. Product jjj of nest iii has a revenue rij≥0r_{ij} \ge 0rij​≥0 and a preference weight vij>0v_{ij} > 0vij​>0, with ri1≥⋯≥rinr_{i1} \ge \dots \ge r_{in}ri1​≥⋯≥rin​. Nest iii has a dissimilarity parameter γi>0\gamma_i > 0γi​>0 and a within-nest no-purchase weight vi0≥0v_{i0} \ge 0vi0​≥0; v0≥0v_0 \ge 0v0​≥0 is the weight of choosing no nest. For an assortment Si⊆NS_i \subseteq NSi​⊆N,

Vi(Si)=vi0+∑j∈Sivij,Ri(Si)=∑j∈SirijvijVi(Si),Ri(∅)=0,V_i(S_i) = v_{i0} + \sum_{j \in S_i} v_{ij}, \qquad R_i(S_i) = \frac{\sum_{j \in S_i} r_{ij} v_{ij}}{V_i(S_i)},\quad R_i(\emptyset)=0,Vi​(Si​)=vi0​+j∈Si​∑​vij​,Ri​(Si​)=Vi​(Si​)∑j∈Si​​rij​vij​​,Ri​(∅)=0,

and the expected revenue of (S1,…,Sm)(S_1, \dots, S_m)(S1​,…,Sm​) is Π=∑iVi(Si)γiRi(Si)/(v0+∑iVi(Si)γi)\Pi = \sum_i V_i(S_i)^{\gamma_i} R_i(S_i) / (v_0 + \sum_i V_i(S_i)^{\gamma_i})Π=∑i​Vi​(Si​)γi​Ri​(Si​)/(v0​+∑i​Vi​(Si​)γi​). The optimal expected revenue Z∗Z^*Z∗ is the optimal value of the linear program

(3)min⁡ xs.t.v0x≥∑i∈Myi,yi≥Vi(Si)γi(Ri(Si)−x)  ∀Si⊆N, i∈M,\text{(3)}\quad \min\ x \quad\text{s.t.}\quad v_0 x \ge \sum_{i \in M} y_i,\qquad y_i \ge V_i(S_i)^{\gamma_i}\big(R_i(S_i) - x\big)\ \ \forall S_i \subseteq N,\ i \in M,(3)min xs.t.v0​x≥i∈M∑​yi​,yi​≥Vi​(Si​)γi​(Ri​(Si​)−x)  ∀Si​⊆N, i∈M,

and problem (4) is the same program with the second family of constraints imposed only for a chosen collection of candidate assortments in each nest.

Throughout, γi≤1\gamma_i \le 1γi​≤1 for every nest and the vi0v_{i0}vi0​ are arbitrary. For a capacity ϵi≥0\epsilon_i \ge 0ϵi​≥0, the knapsack value Ki(ϵi)K_i(\epsilon_i)Ki​(ϵi​) is the largest ∑j∈Srijvij\sum_{j \in S} r_{ij} v_{ij}∑j∈S​rij​vij​ over assortments SSS with ∑j∈Svij≤ϵi\sum_{j \in S} v_{ij} \le \epsilon_i∑j∈S​vij​≤ϵi​ (display (9)). Its continuous relaxation (11) allows fractional zij∈[0,1(vij≤ϵi)]z_{ij} \in [0, \mathbf 1(v_{ij} \le \epsilon_i)]zij​∈[0,1(vij​≤ϵi​)] under the same capacity. The greedy solution z^i(ϵi)\hat z_i(\epsilon_i)z^i​(ϵi​) of (11) fills the capacity with the products of weight at most ϵi\epsilon_iϵi​ in revenue order, each fully while it fits and the next one fractionally, and

S^i(ϵi)={j∈N:z^ij(ϵi)=1}.\hat S_i(\epsilon_i) = \{ j \in N : \hat z_{ij}(\epsilon_i) = 1 \}.S^i​(ϵi​)={j∈N:z^ij​(ϵi​)=1}.

Problem (10) replaces the per-assortment constraints of (3) by yi≥max⁡ϵi≥0(vi0+ϵi)γi[Ki(ϵi)/(vi0+ϵi)−x]y_i \ge \max_{\epsilon_i \ge 0} (v_{i0}+\epsilon_i)^{\gamma_i}[K_i(\epsilon_i)/(v_{i0}+\epsilon_i) - x]yi​≥maxϵi​≥0​(vi0​+ϵi​)γi​[Ki​(ϵi​)/(vi0​+ϵi​)−x].

Formalization targets

Goal: Theorem 10 (p. 24)

Let (x^,y^)(\hat x, \hat y)(x^,y^​) be an optimal solution of (4) with candidate collections {S^i(ϵi):ϵi∈[0,∞]}∪{{j}:j∈N}\{\hat S_i(\epsilon_i) : \epsilon_i \in [0,\infty]\} \cup \{\{j\} : j \in N\}{S^i​(ϵi​):ϵi​∈[0,∞]}∪{{j}:j∈N}. Then

(2x^, 2y^)  is feasible for (3).(2\hat x,\ 2\hat y) \ \text{ is feasible for (3).}(2x^, 2y^​)  is feasible for (3).

Milestones

  1. Per-nest identity (proof of Lemma 9, p. 23). For x≥0x \ge 0x≥0, max⁡SiVi(Si)γi(Ri(Si)−x)=max⁡ϵi≥0(vi0+ϵi)γi[Ki(ϵi)/(vi0+ϵi)−x]\max_{S_i} V_i(S_i)^{\gamma_i}(R_i(S_i) - x) = \max_{\epsilon_i \ge 0}(v_{i0}+\epsilon_i)^{\gamma_i}[K_i(\epsilon_i)/(v_{i0}+\epsilon_i) - x]maxSi​​Vi​(Si​)γi​(Ri​(Si​)−x)=maxϵi​≥0​(vi0​+ϵi​)γi​[Ki​(ϵi​)/(vi0​+ϵi​)−x].
  2. Lemma 9 (p. 23). Problems (3) and (10) have the same optimal solutions.
  3. Relaxation (p. 23). Every feasible point of (9) is feasible for (11), so K^i(ϵi)≥Ki(ϵi)\hat K_i(\epsilon_i) \ge K_i(\epsilon_i)K^i​(ϵi​)≥Ki​(ϵi​).
  4. Greedy solution (pp. 23–24). z^i(ϵi)\hat z_i(\epsilon_i)z^i​(ϵi​) is optimal for (11) and has at most one fractional component.
  5. Sign (A.3, p. 45). x^≥0\hat x \ge 0x^≥0.
  6. Inequalities (28) and (29) (A.3, pp. 45–46). In both cases — z^i(ϵ)\hat z_i(\epsilon)z^i​(ϵ) with and without a fractional component — 2y^i≥(vi0+ϵ)γi[Ki(ϵ)/(vi0+ϵ)−2x^]2\hat y_i \ge (v_{i0}+\epsilon)^{\gamma_i}[K_i(\epsilon)/(v_{i0}+\epsilon) - 2\hat x]2y^​i​≥(vi0​+ϵ)γi​[Ki​(ϵ)/(vi0​+ϵ)−2x^].

Two further statements accompany the goal: the factor-two revenue guarantee obtained from Theorem 10 and Theorem 1 of the paper, and the fact that every S^i(ϵi)\hat S_i(\epsilon_i)S^i​(ϵi​) is one of the at most 1+n21 + n^21+n2 assortments NijkN^k_{ij}Nijk​, the first jjj products by revenue among the kkk lightest.

Significance

Theorem 10 turns an NP-hard assortment problem into a linear program with 1+m1 + m1+m variables and 1+m(1+n+n2)1 + m(1 + n + n^2)1+m(1+n+n2) constraints whose solution is within a factor of two of optimal. The construction is explicit: the candidates are defined by a greedy rule, not by an optimization oracle. The same template, a restricted linear program whose doubled optimum is feasible for the full one, is reused in §6 of the paper for the most general instances, and Lemma 9's knapsack reformulation is the link to the classical approximation theory of knapsack problems (Williamson and Shmoys, 2011).

The theorem is proved in the paper. No machine-checked proof of it, of Lemma 9, or of greedy optimality for the continuous knapsack with an eligibility bound exists on the platform. Formalizing it yields a checked factor-two guarantee and a reusable fractional-knapsack development.

Difficulty

The obvious argument would compare the restricted program (4) with (3) constraint by constraint. That fails: (3) has one constraint per subset of products, and most subsets are not candidates. The comparison has to pass through the knapsack reformulation (10), which requires showing that a maximum over all subsets equals a maximum over a one-dimensional capacity parameter, using γi≤1\gamma_i \le 1γi​≤1 and x≥0x \ge 0x≥0 in an essential way. The second obstacle is that the greedy assortment S^i(ϵi)\hat S_i(\epsilon_i)S^i​(ϵi​) keeps only the fully taken products, so its value can fall short of the continuous knapsack value, and no single candidate assortment need attain the knapsack bound. With dissimilarity parameters above one the monotonicity behind the reformulation is lost, and §6 of the paper needs a different factor.

Formalization scope

Products are Fin n (indices 0,…,n−10, \dots, n-10,…,n−1), nests a finite type, and every quantity is real. Powers are Real.rpow; x/0=0x/0 = 0x/0=0, which gives Ri(∅)=0R_i(\emptyset) = 0Ri​(∅)=0. An optimal solution of a linear program is a feasible pair whose xxx is minimal among feasible pairs. The constraint "yi≥max⁡ϵi≥0(… )y_i \ge \max_{\epsilon_i \ge 0}(\dots)yi​≥maxϵi​≥0​(…)" of (10) is stated in constraint form, for every ϵi≥0\epsilon_i \ge 0ϵi​≥0, so no real supremum is taken. Ki(ϵ)K_i(\epsilon)Ki​(ϵ) is defined for ϵ≥0\epsilon \ge 0ϵ≥0 only; its placeholder value for ϵ<0\epsilon < 0ϵ<0 is never used. Ties in revenue (and, for NijkN^k_{ij}Nijk​, in weight) are broken by index. The candidate collection is taken over real ϵi≥0\epsilon_i \ge 0ϵi​≥0; ϵi=∞\epsilon_i = \inftyϵi​=∞ adds nothing, since every capacity of at least ∑jvij\sum_j v_{ij}∑j​vij​ already gives S^i=N\hat S_i = NS^i​=N.

Standing assumptions and added hypotheses: γi≤1\gamma_i \le 1γi​≤1 for every nest (the section's assumption) on the goal and on every model milestone; the pins vij>0v_{ij} > 0vij​>0, rij≥0r_{ij} \ge 0rij​≥0, γi>0\gamma_i > 0γi​>0 and the revenue ordering, shared by the series; n≥1n \ge 1n≥1 on Lemma 9, on x^≥0\hat x \ge 0x^≥0 and on (28)/(29), the paper's nonempty NNN; and v0>0v_0 > 0v0​>0 on the factor-two revenue guarantee, where Theorem 1 of the paper fails without it.

The greedy assortments S^i(ϵi)\hat S_i(\epsilon_i)S^i​(ϵi​) are defined explicitly. Quantifying over arbitrary optimal solutions of (11) instead would change the candidate collection and is not the paper's theorem. The goal states feasibility for the full program (3) and does not mention knapsack values, the greedy solution or the case split. A formalization that weakens the conclusion to feasibility for (10), or that drops the singletons from the candidate collection, is not a solution.

Needed infrastructure: fractional knapsack optimality of the greedy rule with an eligibility bound, monotonicity of t↦tγt \mapsto t^{\gamma}t↦tγ and t↦tγ−1t \mapsto t^{\gamma - 1}t↦tγ−1 for γ≤1\gamma \le 1γ≤1, and finite maximization over subsets. The fractional-knapsack lemmas are reusable beyond this mission. Proofs of any milestone, and alternative decompositions of the goal, are welcome.

Selected references

  • J. M. Davis, G. Gallego, H. Topaloglu, Assortment Optimization under Variants of the Nested Logit Model, Operations Research 62(2), 2014 (revised manuscript of June 18, 2013, cited here). DOI 10.1287/opre.2014.1256
  • K. T. Talluri, G. J. van Ryzin, Revenue Management Under a General Discrete Choice Model of Consumer Behavior, Management Science 50(1), 15–33, 2004. DOI 10.1287/mnsc.1030.0147
  • D. P. Williamson, D. B. Shmoys, The Design of Approximation Algorithms, Cambridge University Press, 2011. DOI 10.1017/CBO9780511921735
11 thms0 active usersReviewed
Graph TheoryLinear algebraTheoretical Computer Science·Captain: mikedeng1

Explicit Expanders of Every Degree and Size 3: Deleting Far-Apart Tree-Like Vertices of a Near-Ramanujan Graph and Matching Their Neighbours Keeps λ ≤ 2√(d−1) + εResearch Paper

Motivation

Sparse graphs with small nontrivial eigenvalues, expanders, are basic objects in combinatorics and theoretical computer science. They are used in error-correcting codes, derandomization, sorting networks, and the analysis of random walks. The Alon–Boppana bound says that a ddd-regular graph on nnn vertices has a nontrivial eigenvalue of absolute value at least 2d−1−o(1)2\sqrt{d-1}-o(1)2d−1​−o(1) (Alon 1986; Nilli 1991). Graphs that reach 2d−12\sqrt{d-1}2d−1​ are Ramanujan graphs.

The classical explicit Ramanujan graphs of Lubotzky, Phillips and Sarnak (1988) and Margulis exist only for degrees d=p+1d = p+1d=p+1 with ppp prime, and only for very sparse sequences of vertex counts. Constructions with λ≤2d−1+ε\lambda\le 2\sqrt{d-1}+\varepsilonλ≤2d−1​+ε for every degree came from Mohanty, O'Donnell and Paredes (STOC 2020, arXiv:1909.06988), but their graphs also do not have every number of vertices. Alon (arXiv:2003.11673, Combinatorica 41, 2021) asked for near-Ramanujan graphs of every degree and every large size. Theorem 1.3 of that paper answers this up to ε\varepsilonε: for every ddd, every ε>0\varepsilon>0ε>0 and every large nnn with ndndnd even there is an explicit (n,d,λ)(n,d,\lambda)(n,d,λ)-graph with λ≤2d−1+ε\lambda\le 2\sqrt{d-1}+\varepsilonλ≤2d−1​+ε.

Setting

A (n,d,λ)(n,d,\lambda)(n,d,λ)-graph is a ddd-regular simple graph on nnn vertices in which every nontrivial eigenvalue of the adjacency matrix AAA has absolute value at most λ\lambdaλ. The trivial eigenvalue is ddd, with the constant eigenvector 1\mathbf 11. Equivalently, every eigenvalue μ\muμ of AAA with an eigenvector f≠0f\ne0f=0, ∑vf(v)=0\sum_v f(v)=0∑v​f(v)=0, satisfies ∣μ∣≤λ|\mu|\le\lambda∣μ∣≤λ.

Distances dist⁡(v,w)\operatorname{dist}(v,w)dist(v,w) are graph distances, and they are ∞\infty∞ between components. The kkk-neighbourhood of a vertex vvv is B(v,k)={w:dist⁡(v,w)≤k}B(v,k)=\{w:\operatorname{dist}(v,w)\le k\}B(v,k)={w:dist(v,w)≤k}. The kkk-neighbourhood of an edge uvuvuv is B(u,k)∪B(v,k)B(u,k)\cup B(v,k)B(u,k)∪B(v,k), and NiN_iNi​ is the set of vertices at distance exactly iii from {u,v}\{u,v\}{u,v}. A set contains no cycle if the subgraph induced on it is a forest. A ball contains at most one cycle if its induced subgraph has at most as many edges as vertices.

The construction starts from a ddd-regular graph HHH on a vertex set VVV and a set U⊆VU\subseteq VU⊆V. Write N(U)N(U)N(U) for the set of neighbours of UUU, and let mmm be a perfect matching on N(U)N(U)N(U). Then H′H'H′ is the subgraph induced on V∖UV\setminus UV∖U, MMM is the graph of matching edges {x,m(x)}\{x,m(x)\}{x,m(x)}, and G=H′∪MG=H'\cup MG=H′∪M.

Formalization targets

Goal: Theorem 1.3, relative to the input graph

Let d≥3d\ge3d≥3, ε>0\varepsilon>0ε>0, r=⌈2/ε⌉r=\lceil 2/\varepsilon\rceilr=⌈2/ε⌉. Suppose HHH is an (N,d,2d−1+ε/2)(N,d,2\sqrt{d-1}+\varepsilon/2)(N,d,2d−1​+ε/2)-graph in which the (2r+4)(2r+4)(2r+4)-neighbourhood of every vertex contains at most one cycle, and r≤log⁡d−1Nr\le\log_{d-1}Nr≤logd−1​N. Then for every uuu with ududud even and u≤N/(2d2r+3)u\le N/(2d^{2r+3})u≤N/(2d2r+3),

∃ G on N−u vertices:G is an (N−u, d, 2d−1+ε)-graph.\exists\, G \text{ on } N-u \text{ vertices}:\quad G \text{ is an } \bigl(N-u,\ d,\ 2\sqrt{d-1}+\varepsilon\bigr)\text{-graph}.∃G on N−u vertices:G is an (N−u, d, 2d−1​+ε)-graph.

The hypotheses on HHH are what Theorem 3.3 (Mohanty–O'Donnell–Paredes) supplies, and that theorem is not formalized.

Milestones

  • Lemma 3.1 (p. 10). A ddd-regular graph whose (2r+4)(2r+4)(2r+4)-balls contain at most one cycle has a set UUU with ∣U∣≥n/(2d2r+3)|U|\ge n/(2d^{2r+3})∣U∣≥n/(2d2r+3), cycle-free (r+1)(r+1)(r+1)-balls, and pairwise distances ≥2r+3\ge 2r+3≥2r+3.
  • Lemma 3.2 (p. 11). If the rrr-neighbourhood of an edge uvuvuv contains no cycle and Af=μfA f=\mu fAf=μf with μ≥2d−1\mu\ge2\sqrt{d-1}μ≥2d−1​, then
∑w∈Nif2(w) ≥ ∑w∈Ni−1f2(w),1≤i≤r.\sum_{w\in N_i}f^2(w)\ \ge\ \sum_{w\in N_{i-1}}f^2(w),\qquad 1\le i\le r .w∈Ni​∑​f2(w) ≥ w∈Ni−1​∑​f2(w),1≤i≤r.
  • The variational characterization of nontrivial eigenvalues (§2.4, p. 8).
  • In G=H′∪MG=H'\cup MG=H′∪M: GGG is ddd-regular on ∣V∣−∣U∣|V|-|U|∣V∣−∣U∣ vertices and AG=AH′+AMA_G=A_{H'}+A_MAG​=AH′​+AM​. Matching edges have cycle-free (r−1)(r-1)(r−1)-neighbourhoods and pairwise disjoint rrr-neighbourhoods.
  • Inequalities (9), (10), (11) (p. 13), and the spectral step: for every admissible UUU and mmm, GGG is an (N−∣U∣,d,2d−1+ε)(N-|U|,d,2\sqrt{d-1}+\varepsilon)(N−∣U∣,d,2d−1​+ε)-graph.

Significance

Theorem 1.3 shows that the size restrictions of algebraic Ramanujan constructions cost nothing spectrally: up to an arbitrarily small ε\varepsilonε, the Alon–Boppana bound is attained by explicit graphs on every admissible vertex count. The deletion method is local. It turns any near-Ramanujan graph whose short cycles are sparse into graphs of all nearby sizes, so it applies to future constructions as well. Lemma 3.2 is a self-contained delocalization statement in the tradition of Kahale 1995: eigenvectors of eigenvalues at least 2d−12\sqrt{d-1}2d−1​ in absolute value cannot concentrate near tree-like edges.

The result is proved on paper. To our knowledge none of it is formalized; Mathlib has adjacency matrices, extended graph distance and acyclicity, but no theory of expanders. A complete development would give machine-checked versions of a delocalization lemma, of the greedy selection of far-apart vertices away from short cycles, and of the variational eigenvalue bound for induced subgraphs. It would also check two points the paper passes over. The proof of Theorem 1.3 treats only positive eigenvalues λ≥2d−1\lambda\ge2\sqrt{d-1}λ≥2d−1​. And its claim that the rrr-neighbourhood of a matching edge is cycle-free fails when two deleted vertices are at distance exactly 2r+32r+32r+3. This mission states the corrected forms (see Formalization scope).

Difficulty

The spectral bound for GGG does not follow from interlacing alone. Deleting vertices is harmless, since by (9) the quadratic form of H′H'H′ is controlled by HHH. But the added matching contributes up to ∑x∈N(U)f(x)2\sum_{x\in N(U)}f(x)^2∑x∈N(U)​f(x)2 to ftAGff^tA_GfftAG​f, which can be as large as ∥f∥2\|f\|^2∥f∥2 for an eigenvector concentrated on N(U)N(U)N(U). The obvious estimate therefore gives only λ≤2d−1+1+ε/2\lambda\le 2\sqrt{d-1}+1+\varepsilon/2λ≤2d−1​+1+ε/2. Closing the gap requires showing that an eigenvector of a large eigenvalue spreads its mass over the rrr layers around each matching edge (Lemma 3.2). That in turn needs those neighbourhoods to be trees in GGG and pairwise disjoint, which is where Lemma 3.1's choice of UUU is used. The combinatorial part, tracking distances and cycles in GGG when GGG mixes edges of HHH with matching edges, is the main formalization burden.

Formalization scope

  • Representation. Vertex sets are finite types. Graphs are Mathlib SimpleGraphs with real adjacency matrices adjMatrix ℝ. The (n, d, λ) predicate requires IsRegularOfDegree d, symmetry, row sums ddd, and ∣μ∣≤λ|\mu|\le\lambda∣μ∣≤λ for every eigenpair (μ,f)(\mu,f)(μ,f) with f≠0f\ne0f=0, ∑f=0\sum f=0∑f=0. Distances use the extended SimpleGraph.edist, never dist (which is 000 across components). Cycle conditions are on induced subgraphs, and "at most one cycle" on a ball is ∣E∣≤∣V∣|E|\le|V|∣E∣≤∣V∣. Deleted vertices are a Finset U; the new graph lives on the subtype {v // v ∉ U}. The matching is a fixed-point-free involution of N(U)N(U)N(U).
  • Explicit quantities replacing the paper's asymptotics. The paper writes "sufficiently large nnn" and u=o(n)u=o(n)u=o(n). The goal instead takes any u≤N/(2d2r+3)u\le N/(2d^{2r+3})u≤N/(2d2r+3) (the size Lemma 3.1 guarantees) with ududud even, plus Lemma 3.1's side condition r≤log⁡d−1Nr\le\log_{d-1}Nr≤logd−1​N. The equality r=⌈2/ε⌉r=\lceil2/\varepsilon\rceilr=⌈2/ε⌉ is used as ⌈2/ε⌉+∈N\lceil2/\varepsilon\rceil_+\in\mathbb N⌈2/ε⌉+​∈N. The paper's "every degree ddd" becomes d≥3d\ge3d≥3, the range of its proof.
  • Corrections. (11) and Lemma 3.2's companion are stated for ∣μ∣≥2d−1|\mu|\ge2\sqrt{d-1}∣μ∣≥2d−1​, both signs. The matching-edge note is stated for the (r−1)(r-1)(r−1)-neighbourhood, which still yields the factor 1/r1/r1/r in (11). Lemma 3.2 itself is stated as printed.
  • Out of scope. Theorem 3.3 ([18]) is a cited input: its graph is the hypothesis HHH. All claims of explicitness and polynomial running time are out of scope, as is §4's remark on applying the method to LPS graphs directly.
  • No trivialization. The input hypotheses are exactly Theorem 3.3's conclusions plus Lemma 3.1's side condition, and they are met by high-girth Ramanujan graphs. No hypothesis mentions the spectrum or Rayleigh quotients of the constructed graph, and the goal's graph must be ddd-regular on exactly N−uN-uN−u vertices.
  • Reusable infrastructure. Welcome contributions include the variational characterization for symmetric matrices with constant row sums, a forest edge-count lemma for balls, BFS-layer structure of cycle-free balls in regular graphs, and the edge-disjoint decomposition AG=AH′+AMA_{G}=A_{H'}+A_MAG​=AH′​+AM​. Each is useful beyond this mission.

Selected references

  • N. Alon, Explicit expanders of every degree and size, arXiv:2003.11673v1, 2020; Combinatorica 41 (2021). https://arxiv.org/abs/2003.11673
  • S. Mohanty, R. O'Donnell, P. Paredes, Explicit near-Ramanujan graphs of every degree, STOC 2020. https://arxiv.org/abs/1909.06988
  • A. Lubotzky, R. Phillips, P. Sarnak, Ramanujan graphs, Combinatorica 8 (1988). https://doi.org/10.1007/BF02126799
  • N. Alon, Eigenvalues and expanders, Combinatorica 6 (1986). https://doi.org/10.1007/BF02579166
  • A. Nilli, On the second eigenvalue of a graph, Discrete Mathematics 91 (1991). https://doi.org/10.1016/0012-365X(91)90112-F
  • N. Kahale, Eigenvalues and expansion of regular graphs, J. ACM 42 (1995). https://doi.org/10.1145/210118.210136
13 thms0 active usersReviewed
Graph Theory·Captain: mikedeng1

The Strong Perfect Graph Theorem VI: A Berge Graph Containing a Long Odd Prism Admits a Proper 2-Join, a Balanced Skew Partition or a Proper Homogeneous PairResearch Paper

Motivation

A graph is perfect if every induced subgraph has chromatic number equal to its clique number. Berge conjectured in 1961 that a graph is perfect exactly when it contains no odd hole and no odd antihole. Chudnovsky, Robertson, Seymour and Thomas proved this strong perfect graph theorem in 2006 (Ann. of Math. 164 (2006), 51–229). Perfect graphs matter beyond graph theory: for them the stable set polytope is described by clique inequalities, so maximum weight stable set and colouring problems become polynomially solvable linear programs (Grötschel, Lovász and Schrijver).

The proof reduces the theorem to a decomposition statement (1.3 of the paper): every Berge graph is basic or admits one of a few decompositions. That statement is in turn proved in twelve steps, listed as 1.8.1–1.8.12. Each step excludes one kind of configuration from a Berge graph without the decompositions. This mission is step 1.8.5, restated as 13.4: it deals with Berge graphs that contain a long odd prism but no appearance of K4K_4K4​. It is the only step of the proof that needs proper 2-joins in the complement and proper homogeneous pairs.

Setting

All graphs are finite and simple; G‾\overline{G}G is the complement of GGG. A path is an induced subgraph which is a path, and its length is its number of edges. An antipath is a path of G‾\overline{G}G. A hole is an induced cycle of length at least 444, and an antihole is the complement of a hole of G‾\overline{G}G. GGG is Berge if every hole and antihole of GGG has even length.

A prism consists of two disjoint triangles {a1,a2,a3}\{a_1,a_2,a_3\}{a1​,a2​,a3​}, {b1,b2,b3}\{b_1,b_2,b_3\}{b1​,b2​,b3​} and three paths PiP_iPi​ from aia_iai​ to bib_ibi​, such that the only edges between different paths are the triangle edges. It is long if some PiP_iPi​ has length >1>1>1, even if all three lengths are even, and odd otherwise.

A subdivision HHH of a graph JJJ replaces every edge of JJJ by a track, these tracks being disjoint except at their ends. JJJ appears in GGG if L(H)L(H)L(H) is isomorphic to an induced subgraph of GGG for some bipartite subdivision HHH of JJJ, where LLL is the line graph.

The decompositions in the conclusion (pp. 53–54):

  • A proper 2-join is a partition (X1,X2)(X_1,X_2)(X1​,X2​) of V(G)V(G)V(G), with disjoint nonempty Ai,Bi⊆XiA_i,B_i\subseteq X_iAi​,Bi​⊆Xi​, such that the only edges between X1X_1X1​ and X2X_2X2​ are all edges between A1A_1A1​ and A2A_2A2​ and all edges between B1B_1B1​ and B2B_2B2​. Every component of G∣XiG|X_iG∣Xi​ meets AiA_iAi​ and BiB_iBi​. If G∣XiG|X_iG∣Xi​ is a path between single vertices AiA_iAi​ and BiB_iBi​, it has odd length ≥3\ge 3≥3.
  • A skew partition is a partition (A,B)(A,B)(A,B) of V(G)V(G)V(G) with AAA not connected and BBB not anticonnected. It is balanced if no odd path joins two nonadjacent vertices of BBB through AAA, and no odd antipath joins two adjacent vertices of AAA through BBB.
  • A proper homogeneous pair is a pair (A,B)(A,B)(A,B) of disjoint nonempty sets such that every other vertex is complete or anticomplete to AAA, and complete or anticomplete to BBB. All four combinations must occur.

The intermediate objects come from Sections 11–13 of the paper:

  • A strip S=(A,C,B)S=(A,C,B)S=(A,C,B): every vertex of V(S)=A∪B∪CV(S)=A\cup B\cup CV(S)=A∪B∪C lies on a rung, a path from AAA to BBB with interior in CCC.
  • A step: two disjoint rungs joined exactly by an edge at each end.
  • A step-connected strip: steps cover V(S)V(S)V(S) and connect AAA and BBB.
  • Left-stars and right-stars: vertices complete to AAA (resp. BBB) and anticomplete to the rest of V(S)V(S)V(S).
  • A banister: a path from a left-star to a right-star whose interior sees nothing of V(S)V(S)V(S).
  • A staircase K=(S,a0-R0-b0)K=(S,a_0\text{-}R_0\text{-}b_0)K=(S,a0​-R0​-b0​): a step-connected strip with a banister of length ≥3\ge 3≥3. A staircase can be maximal or strongly maximal.
  • Three kinds of breaker: sets around a strip or staircase whose presence forces a balanced skew partition.

Formalization targets

Goal: 13.4

Let GGG be Berge with no appearance of K4K_4K4​ in GGG or in G‾\overline{G}G, and suppose GGG contains a long odd prism as an induced subgraph. Then

G or G‾ admits a proper 2-join, or G admits a balanced skew partition, or G admits a proper homogeneous pair.G \text{ or } \overline{G} \text{ admits a proper 2-join, or } G \text{ admits a balanced skew partition, or } G \text{ admits a proper homogeneous pair.}G or G admits a proper 2-join, or G admits a balanced skew partition, or G admits a proper homogeneous pair.

Milestones

In the order of the paper's argument:

  • 11.3: in a Berge graph with no even prism, every rung of a step-connected strip and every banister has odd length.
  • 11.4: under no appearance of K4K_4K4​ and no even prism, no anticonnected set QQQ has the six properties listed in the statement.
  • 11.5: a 1-breaker forces a balanced skew partition.
  • 12.1: relative to a maximal staircase, every outside vertex is of exactly one of three types (minor; major; a star with a neighbour on R0R_0R0​).
  • 12.3: a connected set containing a left-star and attaching to B∪CB\cup CB∪C contains a major vertex or a banister.
  • 12.4: a 2-breaker forces a balanced skew partition.
  • 13.3: a 3-breaker forces a balanced skew partition.

Significance

13.4 removes long prisms from the analysis. Combined with 10.6 (the even prism), it shows that a recalcitrant graph contains no long prism in GGG or G‾\overline{G}G, which places it in the class F5\mathcal F_5F5​ (p. 154). The later steps (double diamonds, odd wheels, pseudowheels, wheels) all assume this. The step-connected strip and staircase method developed here is also the paper's model for growing a maximal structure and then classifying how the rest of the graph attaches to it.

The theorem has been proved since 2006; no machine-checked proof of it or of any of its steps is known. The proof of 13.4 also cites these results of the same paper, which are posed in other missions of this series:

  • 10.6 (the even-prism step), posed in mission V;
  • 7.2 (equal parity of the paths of a prism), posed in mission V;
  • 2.1, 2.4, 2.6, 2.7, 4.2, 4.3, 4.5, 4.6 (the Roussel–Rubio lemma and the skew-partition toolkit), posed in mission II.

Difficulty

The difficulty is the volume of case analysis behind every statement. The paper does not prove the exact analogue of the even-prism result 10.6. It warns (p. 127) that it does not know whether that analogue holds, and adds the two extra outcomes instead. The obvious first idea is to take a long odd prism and analyse attachments to it as in Section 10. The paper does not get the result that way (p. 127): it replaces two of the three paths by a maximal step-connected strip. Controlling how every remaining vertex or connected set attaches to such a strip is what the breaker results do. Maximality is also delicate. "Strongly maximal" refers to staircases of the complement, so GGG and G‾\overline{G}G must be handled in one framework.

Formalization scope

Graphs are SimpleGraph V on a Fintype vertex type with decidable equality. G‾\overline{G}G is Gᶜ.

  • Paths, holes, rungs, banisters. Paths are lists of distinct vertices, adjacent exactly when consecutive (induced). Holes are lists adjacent exactly when cyclically consecutive. Rungs are listed from their end in AAA to their end in BBB, and banisters from the left-star to the right-star.
  • Vertex sets. Connectedness of a vertex set is reachability inside G.induce X, so ∅\emptyset∅ is connected. Anticonnectedness is the same notion in Gᶜ.
  • Prisms. A prism is three paths whose cross adjacencies are exactly the two triangles. Each path has length at least 111, so the triangles are disjoint.
  • Appearances. A subdivision of JJJ is an injection of V(J)V(J)V(J) together with one track per edge. The tracks are internally disjoint and cover every vertex and edge. An appearance is a graph embedding of L(H)L(H)L(H) into GGG, and embeddings reflect adjacency.
  • Staircases and breakers. A staircase is a quadruple (A,C,B,R0)(A,C,B,R_0)(A,C,B,R0​). Maximality quantifies over all staircases of GGG, strong maximality also over staircases of G‾\overline{G}G. The breakers are predicates on these data.

Every hypothesis of the paper is kept. "No appearance of K4K_4K4​" is in GGG and in G‾\overline{G}G for the goal, and in GGG for the milestones. The goal also keeps "Berge", "no even prism" and "no 1-/2-breaker" where the page has them. 12.1's "exactly one" is an exclusive disjunction. 11.4 is a non-existence statement.

A formalization with non-induced paths, with the complement dropped from "one of G,G‾G,\overline{G}G,G", or with a strip, staircase or breaker predicate that no graph satisfies would make the targets empty or false. The definitions here are checked against a concrete graph: the 8-vertex prism with path lengths 1,1,31,1,31,1,3 is a long odd prism, and it carries a staircase. Contributions are welcome: proofs of the milestones, and reusable material on induced paths, holes and line graphs of subdivisions.

Selected references

  • M. Chudnovsky, N. Robertson, P. Seymour, R. Thomas, The strong perfect graph theorem, Annals of Mathematics 164 (2006), 51–229. https://doi.org/10.4007/annals.2006.164.51
  • C. Berge, Färbung von Graphen, deren sämtliche bzw. deren ungerade Kreise starr sind, Wiss. Z. Martin-Luther-Univ. Halle-Wittenberg 10 (1961), 114–115.
  • V. Chvátal, Star-cutsets and perfect graphs, J. Combin. Theory Ser. B 39 (1985), 189–199. https://doi.org/10.1016/0095-8956(85)90049-8
  • V. Chvátal, N. Sbihi, Bull-free Berge graphs are perfect, Graphs and Combinatorics 3 (1987), 127–139. https://doi.org/10.1007/BF01788536
  • G. Cornuéjols, W. H. Cunningham, Compositions for perfect graphs, Discrete Mathematics 55 (1985), 245–254. https://doi.org/10.1016/0012-365X(85)90051-7
  • M. Grötschel, L. Lovász, A. Schrijver, Geometric Algorithms and Combinatorial Optimization, Springer, 1988. https://doi.org/10.1007/978-3-642-97881-4
21 thms0 active usersReviewed
Graph TheoryLinear algebraTheoretical Computer Science·Captain: mikedeng1

Explicit Expanders of Every Degree and Size 2: Attaching New Vertices to a (p+1)-Regular Ramanujan Graph and Adding Loops Keeps Every Nontrivial Eigenvalue at Most √(2(p+1)) + √p + o(1)Research Paper

Motivation

Sparse graphs whose adjacency spectrum is concentrated near zero, expanders, are used throughout theoretical computer science: in error-correcting codes, derandomization, sorting and routing networks, and the construction of pseudorandom objects (Hoory, Linial and Wigderson, survey). The best possible spectral expansion for a ddd-regular graph is governed by the Alon–Boppana bound 2d−12\sqrt{d-1}2d−1​, and graphs attaining it, Ramanujan graphs, were constructed explicitly by Lubotzky, Phillips and Sarnak (LPS 1988) and by Margulis. These constructions exist only for special degrees (d=p+1d = p+1d=p+1 with ppp prime) and special numbers of vertices (orders of PSL(2,Fq)PSL(2,\mathbb F_q)PSL(2,Fq​) or SL(2,Fq)SL(2,\mathbb F_q)SL(2,Fq​)). Applications often need a graph of a prescribed size nnn.

N. Alon's paper Explicit expanders of every degree and size (arXiv:2003.11673v1; Combinatorica 41, 2021) shows how to obtain explicit near-Ramanujan graphs on exactly nnn vertices. This mission formalizes the spectral core of its Theorem 1.2: a Ramanujan graph on mmm vertices can be enlarged to n=m+rn = m + rn=m+r vertices, with degree raised by one, while the nontrivial eigenvalues stay within a constant factor of optimal.

Setting

Let VVV be a finite set of m≥1m \ge 1m≥1 vertices. A (n,d,λ)(n,d,\lambda)(n,d,λ)-graph is a ddd-regular graph on nnn vertices whose adjacency matrix AAA satisfies ∣μ∣≤λ|\mu| \le \lambda∣μ∣≤λ for every nontrivial eigenvalue μ\muμ, that is, every eigenvalue other than the top eigenvalue ddd of the constant vector 1\mathbf 11. For a symmetric AAA with A1=d 1A\mathbf 1 = d\,\mathbf 1A1=d1, the nontrivial eigenvalues are those with an eigenvector f≠0f \ne 0f=0 satisfying ∑vf(v)=0\sum_v f(v) = 0∑v​f(v)=0. Graphs may carry loops, at most one per vertex, and a loop adds one to the degree: it is a diagonal entry 111 of AAA.

Fix an integer p≥0p \ge 0p≥0 and let HHH be an (m,p+1,2p)(m, p+1, 2\sqrt p)(m,p+1,2p​)-graph on VVV, a (p+1)(p+1)(p+1)-regular Ramanujan graph. Let R={u1,…,ur}R = \{u_1, \dots, u_r\}R={u1​,…,ur​} be rrr new vertices and let W1,…,Wr⊆VW_1, \dots, W_r \subseteq VW1​,…,Wr​⊆V be pairwise disjoint sets of p+2p+2p+2 vertices each. Put W=⋃iWiW = \bigcup_i W_iW=⋃i​Wi​ and L=V∖WL = V \setminus WL=V∖W. The graph GGG on U=V∪RU = V \cup RU=V∪R is obtained from HHH by joining each uiu_iui​ to every vertex of WiW_iWi​ and adding one loop at each vertex of LLL. Its adjacency matrix is

AG=AH+AR+AL,A_G = A_H + A_R + A_L,AG​=AH​+AR​+AL​,

where AHA_HAH​ is the adjacency matrix of HHH (zero on RRR), ARA_RAR​ that of the stars joining uiu_iui​ to WiW_iWi​, and ALA_LAL​ the diagonal matrix of the loops. Every vertex of GGG has degree p+2p+2p+2.

Formalization targets

Goal: Theorem 1.2, spectral core

AG is an (m+r,  p+2,  2(p+1)+p+(p+1) rm) matrix.A_G \text{ is an } \Big(m+r,\; p+2,\; \sqrt{2(p+1)} + \sqrt p + \frac{(p+1)\,r}{m}\Big)\text{ matrix.}AG​ is an (m+r,p+2,2(p+1)​+p​+m(p+1)r​) matrix.

The paper states λ≤2(d−1)+d−1+o(1)\lambda \le \sqrt{2(d-1)} + \sqrt{d-1} + o(1)λ≤2(d−1)​+d−1​+o(1) for d=p+2d = p+2d=p+2; its proof gives 2(p+1)+p+o(1)\sqrt{2(p+1)} + \sqrt p + o(1)2(p+1)​+p​+o(1), which is stronger, and the error term it produces is (p+1)r/m(p+1)r/m(p+1)r/m. The goal is parametrised by HHH, rrr and the sets WiW_iWi​, so it does not depend on how mmm and rrr are chosen.

Milestones

  1. The variational characterization of the nontrivial eigenvalues: for a symmetric matrix with constant row sums and λ≥0\lambda \ge 0λ≥0, ∣μ∣≤λ|\mu| \le \lambda∣μ∣≤λ for every nontrivial eigenvalue if and only if ∣ftAf∣≤λ∥f∥2|f^tAf| \le \lambda\|f\|^2∣ftAf∣≤λ∥f∥2 whenever ∑f=0\sum f = 0∑f=0.
  2. The Cauchy–Schwarz display: ∑Uf=0\sum_U f = 0∑U​f=0 implies ∣∑Vf∣2=∣∑Rf∣2≤∣R∣∑Rf2|\sum_V f|^2 = |\sum_R f|^2 \le |R| \sum_R f^2∣∑V​f∣2=∣∑R​f∣2≤∣R∣∑R​f2.
  3. Inequality (3): ∣ftAHf∣≤b2(p+1)+c2 2p|f^tA_Hf| \le b^2(p+1) + c^2\, 2\sqrt p∣ftAH​f∣≤b2(p+1)+c22p​ with b2=(∑Vf)2/mb^2 = (\sum_V f)^2/mb2=(∑V​f)2/m and c2=∑Vf2−b2c^2 = \sum_V f^2 - b^2c2=∑V​f2−b2.
  4. Display (4): ftALf=∑v∈Lf2(v)f^tA_Lf = \sum_{v\in L} f^2(v)ftAL​f=∑v∈L​f2(v).
  5. Inequality (5): ∣ftARf∣≤p+2x∑Rf2+x∑Wf2|f^tA_Rf| \le \frac{p+2}{x}\sum_R f^2 + x\sum_W f^2∣ftAR​f∣≤xp+2​∑R​f2+x∑W​f2 for every x>0x > 0x>0.
  6. Inequality (6): for ∑Uf=0\sum_U f = 0∑U​f=0 and x>0x > 0x>0,
∣ftAGf∣≤(2p+1)∑Lf2+(2p+x)∑Wf2+p+2x∑Rf2+(p+1)rm∑Rf2.|f^tA_Gf| \le (2\sqrt p+1)\sum_L f^2 + (2\sqrt p+x)\sum_W f^2 + \frac{p+2}{x}\sum_R f^2 + (p+1)\frac rm \sum_R f^2.∣ftAG​f∣≤(2p​+1)L∑​f2+(2p​+x)W∑​f2+xp+2​R∑​f2+(p+1)mr​R∑​f2.

Significance

With HHH the Lubotzky–Phillips–Sarnak graph on m=∣SL(2,Fq)∣m = |SL(2,\mathbb F_q)|m=∣SL(2,Fq​)∣ vertices for the largest suitable prime qqq with m≤nm \le nm≤n, and r=n−mr = n - mr=n−m, the distribution of primes in arithmetic progressions gives r=o(m)r = o(m)r=o(m), and the goal yields an explicit (n,p+2,λ)(n, p+2, \lambda)(n,p+2,λ)-graph with λ≤(1+2)d−1+o(1)\lambda \le (1+\sqrt2)\sqrt{d-1} + o(1)λ≤(1+2​)d−1​+o(1) for every sufficiently large nnn. This is within a factor of about 1.211.211.21 of the Ramanujan bound 2d−12\sqrt{d-1}2d−1​, for every number of vertices, by an elementary modification of an existing graph. The statement is useful independently of LPS: any Ramanujan graph, or any graph with a bound on its nontrivial eigenvalues, can be padded to a nearby size in the same way.

The result is proved in the paper. No formalization of it, of the (n,d,λ)(n,d,\lambda)(n,d,λ) notion, or of the variational characterization of nontrivial eigenvalues for regular graphs exists on the platform. The mission produces a checked version of the spectral argument, and the variational characterization (milestone 1) is a general fact about symmetric matrices with constant row sums that applies to any spectral expander argument.

Difficulty

The vertices of WWW and LLL lie in the old graph HHH, whose spectrum is controlled, but the new vertices of RRR are not; and a vector orthogonal to 1\mathbf 11 on UUU need not be orthogonal to the constant vector on VVV. Bounding ftAGff^tA_GfftAG​f by applying the Ramanujan bound for HHH to fff restricted to VVV therefore fails: the restriction has a component along the trivial eigenvector of HHH, whose eigenvalue p+1p+1p+1 is large. The argument must show that this component is small, of order r/mr/mr/m, and must balance the star edges between RRR and WWW against the loops on LLL so that every vertex class gets the same coefficient. The naive bound ∣ftARf∣≤∥AR∥ ∥f∥2=p+2 ∥f∥2|f^tA_Rf| \le \|A_R\|\,\|f\|^2 = \sqrt{p+2}\,\|f\|^2∣ftAR​f∣≤∥AR​∥∥f∥2=p+2​∥f∥2 added to 2p2\sqrt p2p​ for HHH and 111 for LLL gives a constant larger than 2(p+1)+p\sqrt{2(p+1)}+\sqrt p2(p+1)​+p​; the stated constant needs the weighted estimate.

On the Lean side, milestone 1 concerns the spectrum of a symmetric matrix on the invariant subspace 1⊥\mathbf 1^\perp1⊥, while Mathlib states the spectral theorem for the whole space.

Formalization scope

  • Vertices of GGG are the disjoint union V⊕Fin rV \oplus \mathrm{Fin}\, rV⊕Finr. GGG is represented by its real adjacency matrix, since it has loops; HHH is a Mathlib SimpleGraph with adjMatrix.
  • The (n,d,λ)(n,d,\lambda)(n,d,λ) predicate is stated for matrices: ∣V∣=n|V| = n∣V∣=n, symmetry, A1=d 1A\mathbf 1 = d\,\mathbf 1A1=d1, and ∣μ∣≤λ|\mu| \le \lambda∣μ∣≤λ for every eigenpair (μ,f)(\mu, f)(μ,f) with f≠0f \ne 0f=0 and ∑f=0\sum f = 0∑f=0. For simple graphs, ddd-regularity is added.
  • The paper's o(1)o(1)o(1) terms are replaced by the explicit quantities its proof produces: (p+1)r/m(p+1)r/m(p+1)r/m in the goal, and (p+1)rm∑Rf2(p+1)\frac rm\sum_R f^2(p+1)mr​∑R​f2 in (6). Inequality (3) is stated with the corrected relation b2+c2=∑Vf2b^2 + c^2 = \sum_V f^2b2+c2=∑V​f2; the paper's "b2+c2=1b^2 + c^2 = 1b2+c2=1" holds only for unit restrictions.
  • The bound uses p=d−2\sqrt p = \sqrt{d-2}p​=d−2​, as in the proof and the abstract, which implies the printed d−1\sqrt{d-1}d−1​.
  • ppp is any natural number. The hypothesis "ppp prime, p≡1(mod4)p \equiv 1 \pmod 4p≡1(mod4)" serves only to obtain HHH from LPS, and HHH is a hypothesis here. The sets WiW_iWi​ are arbitrary pairwise disjoint sets of size p+2p+2p+2, not the paper's consecutive blocks of a numbering of SL(2,Fq)SL(2,\mathbb F_q)SL(2,Fq​).
  • Out of scope: the existence of the prime qqq and the estimate n−m=o(m)n - m = o(m)n−m=o(m); the numbering of SL(2,Fq)SL(2,\mathbb F_q)SL(2,Fq​); the "strongly explicit" and polynomial-time claims; the LPS construction (Theorem 2.1, cited); the variant that replaces loops by a matching for even nnn.
  • The goal cannot be satisfied trivially: dropping the condition ∑f=0\sum f = 0∑f=0 makes it false, since p+2p+2p+2 is always an eigenvalue, and for p≥2p \ge 2p≥2 and small r/mr/mr/m the bound is below p+2p + 2p+2 (for p=5p = 5p=5 it is about 5.70+6r/m5.70 + 6r/m5.70+6r/m). At p=1p = 1p=1 the bound 3+2r/m3 + 2r/m3+2r/m is at least the degree 333, so that case holds trivially; it is the paper's statement there as well.
  • Contributions welcome: proofs of each milestone, especially the variational characterization, which is reusable for mission 3 of this series and for any regular-graph spectral argument.

Selected references

  • N. Alon, Explicit expanders of every degree and size, arXiv:2003.11673v1, 2020; Combinatorica 41 (2021). https://arxiv.org/abs/2003.11673
  • A. Lubotzky, R. Phillips, P. Sarnak, Ramanujan graphs, Combinatorica 8 (1988) 261–277. https://doi.org/10.1007/BF02126799
  • S. Hoory, N. Linial, A. Wigderson, Expander graphs and their applications, Bull. AMS 43 (2006) 439–561. https://doi.org/10.1090/S0273-0979-06-01126-8
9 thms0 active usersReviewed
Number Theory·Captain: mikedeng1

Explicit Expanders of Every Degree and Size 1: For Distinct Primes q₁, q₂, Every Large n Has an LPS Vertex Count Q(q₁, q₂, s, t) Between n and n + o(n)Research Paper

Motivation

An (n,d,λ)(n,d,\lambda)(n,d,λ)-graph is a ddd-regular graph on nnn vertices in which every eigenvalue of the adjacency matrix other than the top eigenvalue ddd has absolute value at most λ\lambdaλ. Graphs of this kind with λ\lambdaλ small compared with ddd are expanders. They are used in derandomization, error-correcting codes, sorting networks and many other constructions in theoretical computer science. A Ramanujan graph achieves λ≤2d−1\lambda\le 2\sqrt{d-1}λ≤2d−1​, which is asymptotically optimal.

The classical explicit Ramanujan graphs of Lubotzky, Phillips and Sarnak (LPS, 1988) exist only for special vertex counts, such as q(q2−1)/2q(q^2-1)/2q(q2−1)/2 for a prime qqq or the size of a quaternion group modulo mmm. N. Alon, in Explicit expanders of every degree and size (arXiv:2003.11673, 2020; Combinatorica 41, 2021), asks for explicit (n,d,λ)(n,d,\lambda)(n,d,λ)-graphs with λ≤(2+o(1))d\lambda\le(2+o(1))\sqrt dλ≤(2+o(1))d​ for every degree ddd and every number of vertices nnn. His Proposition 1.1 builds such graphs out of LPS graphs, and one ingredient is purely arithmetic. If the available vertex counts of an LPS family are dense enough, so that for every large nnn one is within a factor 1+o(1)1+o(1)1+o(1) of nnn, then a few vertices can be added or removed while keeping the spectral bound.

This mission formalizes that ingredient, Lemma 2.2 of the paper (p. 7). It concerns the vertex counts of the LPS graphs H(p,q1sq2t)H(p,q_1^sq_2^t)H(p,q1s​q2t​) for two fixed primes q1,q2q_1,q_2q1​,q2​ and all exponents s,t≥1s,t\ge1s,t≥1. With q1,q2q_1,q_2q1​,q2​ fixed, the construction needs no large primes, which is why Proposition 1.1 is strongly explicit for every fixed degree.

Setting

Fix natural numbers q1,q2q_1,q_2q1​,q2​. For natural numbers s,ts,ts,t define the LPS vertex count

Q(q1,q2,s,t)=q13(s−1) q23(t−1)⋅q1(q1−1)(q1+1)2⋅q2(q2−1)(q2+1)2.Q(q_1,q_2,s,t)=q_1^{3(s-1)}\,q_2^{3(t-1)}\cdot\frac{q_1(q_1-1)(q_1+1)}{2}\cdot\frac{q_2(q_2-1)(q_2+1)}{2}.Q(q1​,q2​,s,t)=q13(s−1)​q23(t−1)​⋅2q1​(q1​−1)(q1​+1)​⋅2q2​(q2​−1)(q2​+1)​.

When q1,q2q_1,q_2q1​,q2​ are distinct primes congruent to 111 modulo 4p4p4p and s,t≥1s,t\ge1s,t≥1, this is the number of vertices of the LPS Cayley graph H(p,q1sq2t)H(p,q_1^sq_2^t)H(p,q1s​q2t​) (§2.3, p. 6). Lemma 2.2 itself only requires q1,q2q_1,q_2q1​,q2​ to be distinct primes. Both fractions are integers, because q(q−1)(q+1)q(q-1)(q+1)q(q−1)(q+1) is a product of three consecutive integers.

The proof uses the real number α=log⁡q1/log⁡q2\alpha=\log q_1/\log q_2α=logq1​/logq2​, and the fractional part {x}=x−⌊x⌋∈[0,1)\{x\}=x-\lfloor x\rfloor\in[0,1){x}=x−⌊x⌋∈[0,1), written x mod 1x \bmod 1xmod1 in the paper.

Formalization targets

Goal: Lemma 2.2

For distinct primes q1,q2q_1,q_2q1​,q2​ there is a function g:N→Rg:\mathbb N\to\mathbb Rg:N→R with g(n)=o(n)g(n)=o(n)g(n)=o(n) such that for all sufficiently large nnn there are integers s,t≥1s,t\ge1s,t≥1 with

n≤Q(q1,q2,s,t)≤n+g(n).n\le Q(q_1,q_2,s,t)\le n+g(n).n≤Q(q1​,q2​,s,t)≤n+g(n).

Equivalently, the ratio between consecutive elements of {Q(q1,q2,s,t):s,t≥1}\{Q(q_1,q_2,s,t):s,t\ge1\}{Q(q1​,q2​,s,t):s,t≥1} tends to 111. The statement fixes no rate for ggg, matching the paper's o(n)o(n)o(n).

Milestones, in the order the proof of Lemma 2.2 uses them (p. 7)

  1. For distinct primes q1,q2q_1,q_2q1​,q2​, α=log⁡q1/log⁡q2\alpha=\log q_1/\log q_2α=logq1​/logq2​ is irrational.
  2. For irrational α\alphaα and every δ>0\delta>0δ>0 there is k1≥1k_1\ge1k1​≥1 with 0<{k1α}<δ0<\{k_1\alpha\}<\delta0<{k1​α}<δ.
  3. For distinct primes q1,q2q_1,q_2q1​,q2​ and every μ>0\mu>0μ>0 there are k1≥1k_1\ge1k1​≥1, k2≥0k_2\ge0k2​≥0 with
1≤q1k1q2k2≤1+μ.1\le \frac{q_1^{k_1}}{q_2^{k_2}}\le 1+\mu .1≤q2k2​​q1k1​​​≤1+μ.
  1. For distinct primes, μ>0\mu>0μ>0 and k1≥1k_1\ge1k1​≥1, if 1≤q1k1/q2k2≤1+μ1\le q_1^{k_1}/q_2^{k_2}\le1+\mu1≤q1k1​​/q2k2​​≤1+μ, then for s>k1s>k_1s>k1​ and t≥1t\ge1t≥1
1≤Q(q1,q2,s,t)Q(q1,q2,s−k1,t+k2)≤(1+μ)3.1\le\frac{Q(q_1,q_2,s,t)}{Q(q_1,q_2,s-k_1,t+k_2)}\le(1+\mu)^3 .1≤Q(q1​,q2​,s−k1​,t+k2​)Q(q1​,q2​,s,t)​≤(1+μ)3.

Significance

The result. Lemma 2.2 is the step of Proposition 1.1 that turns a family of Ramanujan graphs with sparse vertex counts into a family whose vertex counts approximate every large nnn up to a factor 1+o(1)1+o(1)1+o(1). The deviation from nnn can then be absorbed by the general packing argument of §2.1 of the paper. Without it, the construction would have to search for large primes depending on nnn, and the result would be explicit but not strongly explicit.

The formalization. The lemma is proved in the paper, in one paragraph. To our knowledge no machine-checked proof exists, and the platform has no statement of it. The only related platform item is the irrationality of log⁡2/log⁡3\log 2/\log 3log2/log3, the case q1=2q_1=2q1​=2, q2=3q_2=3q2​=3 of milestone 1. Formalizing the lemma requires a quantitative inhomogeneous step that the paper leaves implicit ("implying the desired result"). Its last sentence also contains a misprint that the formalization corrects (see below). The Diophantine milestones 1–3 are reusable for any argument about the multiplicative density of {q1aq2b}\{q_1^a q_2^b\}{q1a​q2b​}, for instance the ratio of consecutive elements of {2a3b}\{2^a3^b\}{2a3b}.

Difficulty

Each milestone is short. The difficulty is in making the paper's last sentence ("implying the desired result") into a proof. Taking sss or ttt large separately does not work: changing sss or ttt by one multiplies QQQ by q13q_1^3q13​ or q23q_2^3q23​, a fixed factor larger than 111, so the values obtained by varying one exponent leave gaps of a constant ratio. The bound must hold for every large nnn, not just along a subsequence, and the exponents must stay positive throughout. Milestone 2 is the classical fact that the multiples of an irrational number are dense modulo 111. It needs a pigeonhole argument, not just the irrationality.

Formalization scope

  • Representation. QQQ is a natural-number-valued Lean definition given by the explicit formula above, not the cardinality of a quaternion group. The graph-count interpretation needs the section's congruence conditions on p,q1,q2p,q_1,q_2p,q1​,q2​; Lemma 2.2 is an arithmetic statement for all distinct primes q1,q2q_1,q_2q1​,q2​. Every use of QQQ in the theorems has positive s,ts,ts,t, so natural-number subtraction s−1s-1s−1 is exact, and the division by 222 is exact since q(q−1)(q+1)q(q-1)(q+1)q(q−1)(q+1) is even. Ratios and the bound n+g(n)n+g(n)n+g(n) are computed in R\mathbb RR after casting. The logarithm is Real.log, and the fractional part is Int.fract.
  • o(n). The paper writes n≤Q≤n+o(n)n\le Q\le n+o(n)n≤Q≤n+o(n) for every large nnn. The goal states it literally: ∃g, g=o(n)\exists g,\ g=o(n)∃g, g=o(n) (Mathlib IsLittleO along atTop) and, eventually in nnn, ∃s,t≥1\exists s,t\ge1∃s,t≥1 with n≤Q≤n+g(n)n\le Q\le n+g(n)n≤Q≤n+g(n). This is equivalent to the form "for every μ>0\mu>0μ>0, every large nnn has s,t≥1s,t\ge1s,t≥1 with n≤Q≤(1+μ)nn\le Q\le(1+\mu)nn≤Q≤(1+μ)n". The paper's proof yields the second form, with the factor (1+μ)3(1+\mu)^3(1+μ)3 for arbitrary μ\muμ.
  • Misprint. The paper's last sentence compares Q(q1,q2,s,t)Q(q_1,q_2,s,t)Q(q1​,q2​,s,t) with Q(q1,q2,s−k1,t−k2)Q(q_1,q_2,s-k_1,t-k_2)Q(q1​,q2​,s−k1​,t−k2​) for s,t≥max⁡{k1,k2}s,t\ge\max\{k_1,k_2\}s,t≥max{k1​,k2​}. As printed the ratio is q13k1q23k2q_1^{3k_1}q_2^{3k_2}q13k1​​q23k2​​, which is not close to 111. Milestone 4 uses the intended pair (s−k1,t+k2)(s-k_1,t+k_2)(s−k1​,t+k2​), with s>k1s>k_1s>k1​ so that s−k1≥1s-k_1\ge1s−k1​≥1.
  • Ruling out a trivial reading. The lower bound n≤Q(q1,q2,s,t)n\le Q(q_1,q_2,s,t)n≤Q(q1​,q2​,s,t) alone holds for every nnn by taking sss large. The content of the goal is the upper bound with a sublinear excess, and s,ts,ts,t must be positive. In milestone 3 the condition k1≥1k_1\ge1k1​≥1 excludes the trivial witness k1=k2=0k_1=k_2=0k1​=k2​=0.
  • Out of scope. The derivation of QQQ as the vertex count of H(p,q1sq2t)H(p,q_1^sq_2^t)H(p,q1s​q2t​) (Hensel's lemma and the Chinese remainder theorem, pp. 6–7), Theorem 2.1 (the LPS graphs are Ramanujan, cited from Lubotzky–Phillips–Sarnak), Proposition 1.1, the §2.1 packing argument, and every running-time ("explicit", "strongly explicit") claim.
  • Infrastructure. Only Mathlib is needed: unique factorization for milestone 1, Int.fract and a pigeonhole or Dirichlet-approximation argument for milestone 2, and real exponentiation and asymptotics for the goal. Contributions of alternative proofs of milestone 2, for instance via Mathlib's Dirichlet approximation theorem, are welcome.

Selected references

  • N. Alon, Explicit expanders of every degree and size, arXiv:2003.11673v1, 2020; Combinatorica 41 (2021). https://arxiv.org/abs/2003.11673 , https://doi.org/10.1007/s00493-020-4429-x
  • A. Lubotzky, R. Phillips, P. Sarnak, Ramanujan graphs, Combinatorica 8 (1988) 261–277. https://doi.org/10.1007/BF02126799
  • S. Hoory, N. Linial, A. Wigderson, Expander graphs and their applications, Bull. AMS 43 (2006) 439–561. https://doi.org/10.1090/S0273-0979-06-01126-8
6 thms0 active usersReviewed
Linear algebraOperations Research·Captain: mikedeng1

On the Abstract Properties of Linear Dependence 5: The Seven-Element Fano Matroid Corresponds to No Real MatrixResearch Paper

Motivation

Whitney's 1935 paper On the Abstract Properties of Linear Dependence introduced matroids: finite sets equipped with a rank function, or equivalently a family of independent sets, obeying a few postulates abstracted from the linear dependence of the columns of a matrix. The obvious first question about such an abstraction is whether it is genuinely more general than its model, that is, whether there are matroids that do not arise from any matrix. Section 16 of the paper answers it with a seven-element example, now called the Fano matroid F7F_7F7​, and proves that no real matrix corresponds to it.

The question has had a long life. Representability of matroids over a given field is a central theme of matroid theory: Tutte (1958) characterized the matroids representable over the field with two elements by a single excluded minor, the four-point line U2,4U_{2,4}U2,4​, and the regular matroids by three excluded minors, U2,4U_{2,4}U2,4​, F7F_7F7​ and its dual; and Seymour's decomposition of regular matroids (1980) rests on the same objects. Whitney's §16 is the starting point of this line: the first proof that the abstract postulates admit matroids outside linear algebra over R\mathbb RR.

Timeline:

  • 1935. Whitney defines matroids, the circuit matrix of a matrix, and proves (§16) that the seven-element matroid M′M'M′ corresponds to no real matrix; in a footnote he credits Saunders MacLane with finding that M′M'M′ corresponds to no matrix and identifying it with a finite projective geometry. On p. 533 he exhibits a matrix of integers mod 2 for M′M'M′.
  • 1958. Tutte characterizes binary and regular matroids by excluded minors; F7F_7F7​ appears as an excluded minor for regularity (Tutte 1958).

Setting

Let M=(aij)\mathbf M=(a_{ij})M=(aij​) be an m×nm\times nm×n matrix with columns C1,…,CnC_1,\dots,C_nC1​,…,Cn​. For a set NNN of columns, let r(N)r(N)r(N) be the rank of the submatrix formed by those columns. Regarding the columns as abstract elements gives a matroid MMM on {C1,…,Cn}\{C_1,\dots,C_n\}{C1​,…,Cn​} with rank function rrr: the matroid of M\mathbf MM. A matroid corresponds to M\mathbf MM if it is the matroid of M\mathbf MM, with elements matched to columns.

A circuit of a matroid is a minimal dependent set. For a circuit P={i1,…,ip}P=\{i_1,\dots,i_p\}P={i1​,…,ip​} of the matroid of M\mathbf MM, there are numbers b1,…,bnb_1,\dots,b_nb1​,…,bn​ with ∑jaijbj=0\sum_j a_{ij}b_j=0∑j​aij​bj​=0 for every row iii, and bj≠0b_j\neq 0bj​=0 exactly for j∈Pj\in Pj∈P; the set of such vectors is written Zi1⋯ipZ_{i_1\cdots i_p}Zi1​⋯ip​​ when only the support condition is meant. Stacking one such row per circuit gives the circuit matrix M′\mathbf M'M′ of M\mathbf MM, determined up to nonzero factors on its rows.

A fundamental set of circuits of a matroid MMM with nullity n(M)=ρ(M)−r(M)n(M)=\rho(M)-r(M)n(M)=ρ(M)−r(M) (ρ\rhoρ the number of elements) is a family of circuits P1,…,PqP_1,\dots,P_qP1​,…,Pq​ with q=n(M)q=n(M)q=n(M) such that the elements can be ordered e1,…,ene_1,\dots,e_ne1​,…,en​ with en−q+i∈Pie_{n-q+i}\in P_ien−q+i​∈Pi​ and en−q+j∉Pie_{n-q+j}\notin P_ien−q+j​∈/Pi​ for j>ij>ij>i; it is strict if en−q+j∉Pie_{n-q+j}\notin P_ien−q+j​∈/Pi​ for every j≠ij\neq ij=i.

The matroid M′M'M′ of §16 has elements 1,…,71,\dots,71,…,7; its bases (maximal independent sets) are all three-element sets except

124,135,167,236,257,347,456.(16.1)124,\quad 135,\quad 167,\quad 236,\quad 257,\quad 347,\quad 456. \qquad (16.1)124,135,167,236,257,347,456.(16.1)

Formalization targets

Goal: §16, pp. 529–530

∃ M′and∀m ∀ M∈Rm×7: M′ is not the matroid of M.\exists\,M' \quad\text{and}\quad \forall m\ \forall\,\mathbf M\in\mathbb R^{m\times 7}:\ M' \text{ is not the matroid of } \mathbf M .∃M′and∀m ∀M∈Rm×7: M′ is not the matroid of M.

The number of rows is arbitrary; the existence clause makes the non-existence statement non-vacuous.

Milestones

  1. §12. Every real matrix has a matroid: the ranks of column submatrices satisfy the rank postulates.
  2. §14, (14.1). Every real matrix has a circuit matrix.
  3. Theorem 29. The rows of a fundamental set of circuits form a base for the rows of the circuit matrix, so r(M′)=q=n(M)r(\mathbf M')=q=n(\mathbf M)r(M′)=q=n(M).
  4. Lemma 10. The support of a vector in the row space HHH of a circuit matrix is a union of circuits.
  5. Lemma 11. Two vectors of HHH with the same circuit as support are proportional.
  6. Theorem 32. For a circuit matrix normalised along a strict fundamental set, a minor DDD vanishes iff an associated q×qq\times qq×q minor D′D'D′ vanishes, iff some circuit avoids a prescribed set of columns.
  7. §16, rank of M′M'M′. The rank of a kkk-set is kkk for k≤2k\le 2k≤2, 333 for k≥4k\ge 4k≥4, and for k=3k=3k=3 it is 222 on (16.1) and 333 otherwise.
  8. p. 533. M′M'M′ is the matroid of an explicit 3×73\times 73×7 matrix of integers mod 2.

Significance

The result. The theorem separates the abstract notion of matroid from linear dependence over R\mathbb RR: some matroids are not real-representable. It also exhibits that representability depends on the field, because the same matroid is the matroid of a matrix over the integers mod 2 (milestone 8). Everything later written about representability over particular fields, excluded-minor characterizations, and the gap between abstract and linear matroids starts from this distinction. Theorem 32 is of independent interest: it translates statements about circuits of a represented matroid into the vanishing of minors of a normalised circuit matrix.

Formalizing it. The result is classical and its proof is short on paper, but it is not formalized in Mathlib, which has matroids (Matroid, circuits, ranks) but no column matroid of a matrix with a rank-of-submatrix characterization, no circuit matrix, and no Fano matroid. The mission produces those objects and the bridge lemmas (Theorem 29, Lemmas 10–11, Theorem 32) that connect matroid circuits with linear algebra of the circuit matrix. No machine-checked proof of the non-representability of the Fano matroid over R\mathbb RR in Lean is known to the curators.

Difficulty

The obvious attempt is a direct search: suppose a real m×7m\times 7m×7 matrix has M′M'M′ as its matroid and derive a contradiction from the seven dependent triples. This does not work as stated. Each rank condition is a determinantal (nonlinear) condition on the entries, the number of rows mmm is unbounded, and a representation is determined only up to row operations and column scalings, so there is no finite case check and no single linear computation that settles the question. The contradiction has to come from an argument that is invariant under these symmetries, and the milestones (circuit vectors determined up to scaling, fundamental sets spanning, circuits detected by minors) are what such an argument needs to be stated in. The field also matters: the argument must use that 2≠02\neq 02=0 in R\mathbb RR, since over a field of characteristic 2 the statement is false (milestone 8).

Formalization scope

  • Elements and matrices. Matroids are Mathlib Matroids whose ground set is the whole (finite) type. The Fano matroid lives on Fin 7, Whitney's element kkk being k - 1; the seven triples are written out literally. Matrices are Matrix (Fin m) ι K; "the matroid of M\mathbf MM" means: ground set everything, and the rank M.eRk N of every finite set NNN of columns equals Matrix.rank of the column submatrix.
  • Field. The goal and Lemmas 10–11, Theorems 29 and 32 are stated over R\mathbb RR, as in the paper; the predicate "matroid of a matrix" is stated over any field so that the mod-2 milestone uses the same notion.
  • Circuit matrix. Rows are determined up to nonzero factors, so "circuit matrix" is a predicate on a matrix together with a bijection between its rows and the circuits; every theorem holds for every such choice.
  • Nullity and indices. q=n(M)q=n(M)q=n(M) is written q+r(M)=ρ(M)q+r(M)=\rho(M)q+r(M)=ρ(M) in extended naturals, with no truncated subtraction. In Theorem 32, n=p+qn=p+qn=p+q, the complement of i1,…,isi_1,\dots,i_si1​,…,is​ is given as an order embedding of Fin t with s+t=qs+t=qs+t=q, and determinants are of square submatrices in the paper's row and column order.
  • Ruling out trivial readings. The goal includes the existence of M′M'M′; without it "every matroid with these bases has no real matrix" could hold vacuously. The goal quantifies over every number of rows; fixing m=3m=3m=3 would be a weaker statement.

Reusable beyond this mission: the matroid of a matrix over a field, the circuit matrix, fundamental sets of circuits, and the Fano matroid. Contributions welcome: proofs of the milestones, and a proof of the goal by any route, including one that does not go through Theorem 32.

Selected references

  • H. Whitney, On the Abstract Properties of Linear Dependence, American Journal of Mathematics 57 (1935), 509–533. https://doi.org/10.2307/2371182
  • W. T. Tutte, A homotopy theorem for matroids, I, II, Transactions of the American Mathematical Society 88 (1958), 144–174. https://doi.org/10.2307/1993244
  • J. Oxley, Matroid Theory, 2nd ed., Oxford University Press, 2011. https://doi.org/10.1093/acprof:oso/9780198566946.001.0001
  • O. Veblen and J. W. Young, Projective Geometry, Vol. I, Ginn, 1910 (cited by Whitney for the finite projective geometry).
13 thms0 active usersReviewed
Linear algebraOperations Research·Captain: mikedeng1

On the Abstract Properties of Linear Dependence 6: Every Matroid Satisfying (C*) Is Represented by a Matrix of Integers Mod 2Research Paper

Motivation

Whitney's 1935 paper introduced matroids as an abstraction of linear dependence among the columns of a matrix. Most of the paper works over the real numbers; its appendix asks which matroids arise from matrices of integers mod 2, that is, matrices with entries 0 and 1 in which rank and dependence are computed over the two-element field. These are today's binary matroids. They include the cycle matroids of graphs (Whitney closes the paper by noting that graphs correspond to mod-2 matrices with exactly two ones in each column) and they are the setting of several later structure theorems: Tutte's excluded-minor characterization of binary matroids (Tutte 1958), Seymour's decomposition of regular matroids (Seymour 1980) and Seymour's theory of binary clutters and max-flow min-cut (Seymour 1977), which underlies parts of combinatorial optimization.

Whitney's answer is an intrinsic postulate, (C*), on the circuits of the matroid, stated without reference to any matrix, and a constructive representation theorem (Theorem 37): a matroid satisfying (C*) is the matroid of a mod-2 matrix, and the matrix is unique once the columns of one base are fixed.

Setting

A matroid MMM on elements e1,…,ene_1, \dots, e_ne1​,…,en​ is given by its independent sets; its circuits are its minimal dependent sets, its rank r(M)r(M)r(M) is the size of a base, and its nullity is n(M)=n−r(M)n(M) = n - r(M)n(M)=n−r(M). Here MMM is a Mathlib Matroid (Fin n) whose ground set is all of Fin n.

Subsets of the elements are added mod 2: a sum of finitely many sets is the set of elements lying in an odd number of them (for two sets, the symmetric difference). A cycle is a sum mod 2 of circuits; the empty sum is the null cycle ∅\emptyset∅. A set is a true sum of sets that have no common elements and whose union it is. Postulate (C*) requires that each cycle be a true sum of circuits.

With n=r+qn = r + qn=r+q, a family P1,…,PqP_1, \dots, P_qP1​,…,Pq​ is a strict fundamental set of circuits with respect to en−q+1,…,ene_{n-q+1}, \dots, e_nen−q+1​,…,en​ if q=n(M)q = n(M)q=n(M), each PiP_iPi​ is a circuit, and PiP_iPi​ contains en−q+ie_{n-q+i}en−q+i​ but no other en−q+je_{n-q+j}en−q+j​.

For a matrix M\mathbf MM over the integers mod 2 with columns C1,…,CnC_1, \dots, C_nC1​,…,Cn​, columns are independent (mod 2) if no non-null subset of them sums to the zero column. The matroid corresponding to M\mathbf MM has the column indices as elements and these independent sets.

Formalization targets

Goal: Theorem 37 (p. 533)

Let MMM satisfy (C*), with elements e1,…,ene_1, \dots, e_ne1​,…,en​ and base {e1,…,en−q}\{e_1, \dots, e_{n-q}\}{e1​,…,en−q​}. For every matrix M1\mathbf M_1M1​ mod 2 (any number of rows) whose n−qn - qn−q columns are independent mod 2,

∃! M=(M1∣Cn−q+1⋯Cn)  whose corresponding matroid is M.\exists!\ \mathbf M = (\mathbf M_1 \mid C_{n-q+1} \cdots C_n) \ \text{ whose corresponding matroid is } M.∃! M=(M1​∣Cn−q+1​⋯Cn​)  whose corresponding matroid is M.

Milestones

  1. Theorem 9 (p. 517): if e1,…,en−qe_1, \dots, e_{n-q}e1​,…,en−q​ is a base, there is a unique strict fundamental set of circuits with respect to en−q+1,…,ene_{n-q+1}, \dots, e_nen−q+1​,…,en​.
  2. Appendix, p. 531: (C*) implies the circuit postulate (C₂), for any family of sets.
  3. Theorem 33: under (C*), the circuits are exactly the minimal non-null cycles.
  4. Theorem 34: under (C*), the cycles are exactly the 2q2^q2q sums mod 2 of a strict fundamental set.
  5. Theorem 35: two (C*)-matroids with a common strict fundamental set have the same circuits.
  6. Theorem 36: any P1,…,PqP_1, \dots, P_qP1​,…,Pq​ with en−q+i∈Pi⊆{e1,…,en−q,en−q+i}e_{n-q+i} \in P_i \subseteq \{e_1, \dots, e_{n-q}, e_{n-q+i}\}en−q+i​∈Pi​⊆{e1​,…,en−q​,en−q+i​} is the strict fundamental set of exactly one (C*)-matroid.
  7. Appendix, p. 532: the matroid of a matrix mod 2 exists, satisfies (C*), and its cycles are the supports of the mod-2 dependencies among the columns.

Milestone 7 and the goal together characterize binary matroids as the matroids satisfying (C*).

Significance

The result. Theorem 37 and the p. 532 claim give an intrinsic, matrix-free description of the matroids representable over the two-element field, and Theorem 36 parametrizes all of them by qqq arbitrary subsets of a base. Uniqueness in Theorem 37 says that a binary representation is determined by the columns of one base; in modern terms, binary matroids are uniquely representable over GF(2) up to row operations. Every later theory of binary matroids, including graphic and cographic matroids, Tutte's excluded-minor theorem and Seymour's decomposition, starts from this equivalence.

Formalizing it. The results are proved in the paper and in textbooks (e.g. Oxley, Matroid Theory, Ch. 9) but, at the Mathlib revision used here, there is no notion of a matroid represented by a matrix over a field, and no binary-matroid theory. On Prove2Me, the existing binary objects (SeymourMFMC.Binary.*) are binary clutters defined through blockers, not matroids represented by mod-2 matrices. This mission produces the representation predicate for mod-2 matrices, the cycle space of a matroid, and the equivalence between (C*) and binary representability.

Difficulty

Writing down candidate columns is not the hard part; showing that the matroid of the completed matrix is MMM itself, and not merely a matroid sharing some of its circuits, is. Whitney's example at the end of §9 exhibits two different matroids with a common strict fundamental set, so agreement on fundamental circuits does not by itself identify a matroid; any argument must use (C*) on both the given matroid and the matroid of the matrix. A naive comparison of independent sets column by column does not close this gap. Uniqueness likewise depends on the independence mod 2 of the prescribed columns: without it, different completions can give the same matroid.

Formalization scope

  • Matroids are Mathlib Matroid (Fin (r + q)) with ground set Set.univ; Whitney's eke_kek​ is k - 1, his e1,…,en−qe_1, \dots, e_{n-q}e1​,…,en−q​ is the range of Fin.castAdd q, and en−q+ie_{n-q+i}en−q+i​ is Fin.natAdd r (i - 1). Writing n=r+qn = r + qn=r+q removes natural-number subtraction; qqq is not a free parameter, since {e1,…,er}\{e_1, \dots, e_r\}{e1​,…,er​} is required to be a base.
  • Sums mod 2 count parity of membership (sumMod2); cycles are sums over finite sets of circuits; true sums are unions over finite pairwise-disjoint sets of circuits; (C*) is SatisfiesCStar on the circuit family {C | M.IsCircuit C}. These definitions take the circuit family as a parameter, so that the (C₂) milestone is posed for an arbitrary family of sets, as Whitney poses it.
  • A strict fundamental set includes the nullity condition r(M)+q=ρ(M)r(M) + q = \rho(M)r(M)+q=ρ(M), stated in N∞\mathbb N_\inftyN∞​.
  • Matrices are Matrix (Fin m) (Fin n) (ZMod 2) with any mmm; independence mod 2 of columns is LinearIndepOn (ZMod 2) of the columns (the rows of the transpose). IsMatroidOf M A compares all independent sets, not only bases.
  • Ruled out: the goal is not satisfied by any statement that compares only the bases of one size, by an existence-only statement without uniqueness, or by real (instead of mod-2) independence.
  • Tacit hypotheses made explicit: the matroid's ground set is exactly e1,…,ene_1, \dots, e_ne1​,…,en​ (ρ(M)=n\rho(M) = nρ(M)=n); the elements and matroids are finite.

Contributions welcome: the general fact that the matroid of a vector family over a field exists (a reusable Matroid.ofFun-style construction over any field), the cycle-space lemmas, and proofs of the milestones in any order.

Selected references

  • H. Whitney, On the Abstract Properties of Linear Dependence, American Journal of Mathematics 57 (1935), 509–533. https://doi.org/10.2307/2371182
  • W. T. Tutte, A homotopy theorem for matroids, I, II, Transactions of the AMS 88 (1958), 144–174. https://doi.org/10.2307/1993244
  • P. D. Seymour, The matroids with the max-flow min-cut property, Journal of Combinatorial Theory Ser. B 23 (1977), 189–222. https://doi.org/10.1016/0095-8956(77)90031-4
  • P. D. Seymour, Decomposition of regular matroids, Journal of Combinatorial Theory Ser. B 28 (1980), 305–359. https://doi.org/10.1016/0095-8956(80)90075-1
  • J. Oxley, Matroid Theory, 2nd ed., Oxford University Press, 2011. https://doi.org/10.1093/acprof:oso/9780198566946.001.0001
11 thms0 active usersReviewed
Operations Research·Captain: mikedeng1

On the Abstract Properties of Linear Dependence 3: Two Elements Share a Component Iff Some Circuit Contains BothResearch Paper

Motivation

Hassler Whitney's 1935 paper On the Abstract Properties of Linear Dependence introduced matroids: finite sets of elements carrying an abstract rank function that behaves like the rank of a set of vectors. Part II of the paper opens with the decomposition of a matroid into components. The question it answers is basic to every later use of matroids: when does a matroid split into independent pieces, and how can the pieces be recognized?

For the matroid of a graph (elements = edges, rank = number of vertices minus number of connected pieces spanned) the components are the 2-connected blocks of the graph, and Whitney's theorem recovers the classical fact that two edges lie in a common block exactly when they lie on a common cycle. Whitney had studied separability of graphs in Non-separable and planar graphs (1932), and footnote 11 of the 1935 paper points out that the theorem identifies König's "Glieder" of a graph with components. Matroid connectivity built on this notion runs through later structure theory: Tutte's higher connectivity, Seymour's decomposition of regular matroids, and the matroid minors project all start from the separation of a matroid into components.

Setting

A matroid MMM on a finite ground set EEE is given here by Mathlib's Matroid structure, with rank function r(X)r(X)r(X) for X⊆EX\subseteq EX⊆E (Mathlib's M.eRk X) and circuits, the minimal dependent sets (M.IsCircuit). Whitney treats every subset X⊆EX\subseteq EX⊆E as a matroid in its own right, a submatroid, with the rank function of MMM restricted to subsets of XXX. For sets he writes M1+M2M_1+M_2M1​+M2​ for the union, ρ(N)\rho(N)ρ(N) for the number of elements of NNN, and

n(N)=ρ(N)−r(N)n(N) = \rho(N) - r(N)n(N)=ρ(N)−r(N)

for the nullity of NNN.

Rank is subadditive: r(X1+X2)≤r(X1)+r(X2)r(X_1+X_2)\le r(X_1)+r(X_2)r(X1​+X2​)≤r(X1​)+r(X2​). A submatroid XXX is separable if it can be divided into two disjoint groups X1,X2X_1, X_2X1​,X2​, each containing at least one element, with

r(X)=r(X1)+r(X2),r(X) = r(X_1) + r(X_2),r(X)=r(X1​)+r(X2​),

and non-separable otherwise. Every single element is non-separable. A component of MMM is a maximal non-separable part of MMM: a nonempty non-separable set K⊆EK\subseteq EK⊆E contained in no strictly larger non-separable subset of EEE.

Formalization targets

Goal: Theorem 19

For two distinct elements e1≠e2e_1\neq e_2e1​=e2​ of EEE,

(∃K component of M: e1,e2∈K)  ⟺  (∃P circuit of M: e1,e2∈P).\bigl(\exists K \text{ component of } M:\ e_1, e_2\in K\bigr) \iff \bigl(\exists P \text{ circuit of } M:\ e_1, e_2\in P\bigr).(∃K component of M: e1​,e2​∈K)⟺(∃P circuit of M: e1​,e2​∈P).

Components are defined by the rank function, circuits by dependence; the goal asserts that the two descriptions agree.

Milestones (§10, in the paper's order)

  • Theorem 11. If r(M1+M2)=r(M1)+r(M2)r(M_1+M_2)=r(M_1)+r(M_2)r(M1​+M2​)=r(M1​)+r(M2​), M1′⊆M1M_1'\subseteq M_1M1′​⊆M1​ and M2′⊆M2M_2'\subseteq M_2M2′​⊆M2​, then r(M1′+M2′)=r(M1′)+r(M2′)r(M_1'+M_2')=r(M_1')+r(M_2')r(M1′​+M2′​)=r(M1′​)+r(M2′​).
  • Theorem 12. Under the same rank additivity, a non-separable M′⊆M1+M2M'\subseteq M_1+M_2M′⊆M1​+M2​ lies in M1M_1M1​ or in M2M_2M2​.
  • Theorem 13. Two non-separable sets with a common element have a non-separable union.
  • Theorem 14. Distinct components are disjoint.
  • Theorem 15. The components cover EEE, and no other family of components does.
  • Theorem 16. A set is non-separable of nullity 111 if and only if it is a circuit.
  • Lemma 9. If M1+M2M_1+M_2M1​+M2​ is non-separable, with M1,M2M_1, M_2M1​,M2​ nonempty and disjoint, some circuit inside M1+M2M_1+M_2M1​+M2​ meets both.
  • Theorem 17. A non-separable set of nullity n>0n>0n>0 is built from a circuit by n−1n-1n−1 steps, each adding a set of elements that forms a circuit with elements already present, through non-separable sets of nullity 1,2,…,n1,2,\dots,n1,2,…,n.
  • Theorem 18. For distinct nonempty non-separable M1,…,MpM_1,\dots,M_pM1​,…,Mp​ covering EEE, the following are equivalent: they are the components; they are pairwise disjoint and no circuit meets two of them; r(E)=∑ir(Mi)r(E)=\sum_i r(M_i)r(E)=∑i​r(Mi​).

Significance

The result. Theorem 19 makes the component decomposition computable from circuits alone and shows that "lying on a common circuit" is an equivalence relation on distinct elements, a fact that is not evident from the circuit axioms. Theorem 18 adds that the decomposition is the unique one with additive rank. Together they are the starting point of matroid connectivity: the direct-sum decomposition of a matroid, the reduction of many matroid problems (representability, duality of components, Whitney's own Theorems 24–26 on duals of components) to the connected case, and the higher-connectivity theory that followed.

Formalizing it. The results are classical and proved in the paper; nothing here is open. To our knowledge Mathlib at the pinned revision has no notion of matroid connectivity or components, so this mission produces the first machine-checked development of Whitney's §10: the rank-based definition of separability, the disjoint decomposition into components, the circuit characterization, and the ear-type construction of non-separable matroids (Theorem 17). These are reusable for any later formalization of matroid connectivity, including Whitney's results on duals of components.

Difficulty

The two directions of Theorem 19 rest on different machinery. That two elements on a common circuit lie in one component follows from the rank theory (Theorems 13 and 16). The converse is the substantial direction: a component is defined by the failure of rank additivity, which only says that every division of the component is crossed by some circuit (Lemma 9). It does not directly give one circuit through two prescribed elements. Combining circuits that cross different divisions into a single circuit through both e1e_1e1​ and e2e_2e2​ requires the circuit elimination property together with a minimality argument over subsets of the component; the naive attempt of chaining overlapping circuits from e1e_1e1​ to e2e_2e2​ gives a connected chain of circuits, not one circuit.

Formalization scope

  • Representation. A matroid is Mathlib's Matroid α with [M.Finite]; Whitney's matroids are finite. A submatroid is a subset X⊆X\subseteqX⊆ M.E with the rank M.eRk restricted to its subsets; results that Whitney states for "a matroid M=M1+M2M = M_1 + M_2M=M1​+M2​" are stated for subsets of an ambient finite matroid, which is the same statement applied to the submatroid M1+M2M_1+M_2M1​+M2​.
  • Ranks are Mathlib's ℕ∞-valued M.eRk, finite on a finite matroid, so (10.1) is an equation of natural numbers. Nullity is computed in Z\mathbb ZZ as the number of elements minus the rank.
  • Definitions. IsSeparable M X requires two nonempty, disjoint groups with union XXX and additive rank; without nonemptiness every set would be separable. IsNonSeparable M X adds X⊆X\subseteqX⊆ M.E. IsComponent M K requires KKK nonempty, non-separable, and maximal; nonemptiness excludes the empty set, which is vacuously non-separable.
  • Tacit hypotheses made explicit. In Theorem 19 the two elements are distinct: for e1=e2e_1=e_2e1​=e2​ a coloop is its own component and lies on no circuit. In Theorem 18 the sets M1,…,MpM_1,\dots,M_pM1​,…,Mp​ are distinct and nonempty: a loop listed twice would satisfy (3) but not (2), and an empty set would satisfy (2) and (3) but not (1). In Theorems 11 and 12 the two parts need not be disjoint, as Whitney's use of M1+M2M_1+M_2M1​+M2​ for overlapping sets in Theorem 13 indicates; the statements hold in that generality.
  • Ruled out. Components must not be defined as the classes of the relation "lie on a common circuit": that would make the goal a tautology. Here they are the rank-defined maximal non-separable sets of §10, and circuits are Mathlib's Matroid.IsCircuit.
  • Infrastructure. Solvers will need submodularity of M.eRk and circuit elimination (both in Mathlib), the relation between circuits of M ↾ X and circuits of M inside XXX (Matroid.restrict_isCircuit_iff), and finiteness arguments for maximal non-separable sets. Proofs of the milestones, alternative proofs of the goal, and lemmas relating components to Mathlib's direct sums of matroids are all welcome.

Selected references

  • H. Whitney, On the Abstract Properties of Linear Dependence, American Journal of Mathematics 57 (1935), 509–533. https://doi.org/10.2307/2371182
  • 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
  • D. König, Acta Litterarum ac Scientiarum Szeged, vol. 6, pp. 155–179, as cited by Whitney in footnote 11 (p. 159 for the notion of "Glied").
  • J. Oxley, Matroid Theory, 2nd ed., Oxford University Press, 2011, Chapter 4 (connectivity).
13 thms0 active usersReviewed
Graph Theory·Captain: mikedeng1

The Strong Perfect Graph Theorem IV: A Berge Graph Whose Appearances of K4 Are All Degenerate Is Double Split, Decomposes, or Has No Appearance of K4Research Paper

Perfect graphs and the decomposition of Berge graphs

A graph is perfect if every induced subgraph has chromatic number equal to its clique number. Perfect graphs are the graphs for which colouring and clique problems behave as linear programs do: the stable-set polytope of a perfect graph is described by its clique inequalities, so maximum weight stable sets and minimum colourings can be computed in polynomial time (Grötschel, Lovász & Schrijver 1988). In 1961 Berge conjectured that a graph is perfect exactly when it has no odd hole and no odd antihole. Chudnovsky, Robertson, Seymour and Thomas proved this, the strong perfect graph theorem, in Ann. of Math. 164 (2006).

The proof is a decomposition theorem: every Berge graph is basic or admits one of a few decompositions. The paper reaches it through twelve steps, 1.8.1–1.8.12 (p. 59), each handling graphs that contain a certain configuration. This mission poses step 1.8.3, Theorem 9.6. It handles Berge graphs that contain the line graph of a bipartite subdivision of K4K_4K4​, all such line graphs being degenerate.

Setting

All graphs are finite and simple. G‾\overline{G}G denotes the complement of GGG. A hole is an induced cycle of length at least 444, an antihole is a hole of G‾\overline{G}G, and GGG is Berge if all its holes and antiholes have even length. A path is always an induced path, and an antipath is a path of G‾\overline{G}G. The length of either is its number of edges.

Line graphs and appearances. The line graph L(H)L(H)L(H) has vertex set E(H)E(H)E(H), two edges adjacent when they share an end. HHH is a subdivision of JJJ if it arises from JJJ by replacing every edge by a track (a path, not necessarily induced), these tracks disjoint except for their ends. JJJ appears in GGG if, for some bipartite subdivision HHH of JJJ, L(H)L(H)L(H) is isomorphic to an induced subgraph of GGG; L(H)L(H)L(H) is then an appearance of JJJ. For J=K4J = K_4J=K4​ the appearance is degenerate if some 4-cycle of HHH contains the four vertices of degree three. A K4K_4K4​-enlargement is a 3-connected graph with a proper subgraph isomorphic to a subdivision of K4K_4K4​. An appearance L(H)L(H)L(H) is overshadowed if some branch of HHH of odd length ≥3\ge 3≥3, with ends b1,b2b_1, b_2b1​,b2​, has a vertex of GGG nonadjacent to at most one edge at b1b_1b1​ and at most one edge at b2b_2b2​.

Knots and striations. A knot (P1,P2,Q1,Q2)(P_1, P_2, Q_1, Q_2)(P1​,P2​,Q1​,Q2​) is formed by two paths PiP_iPi​ with ends ai,bia_i, b_iai​,bi​ and two antipaths QjQ_jQj​ with ends xj,yjx_j, y_jxj​,yj​. They are pairwise disjoint and of length ≥1\ge 1≥1, P1P_1P1​ is anticomplete to P2P_2P2​, Q1Q_1Q1​ is complete to Q2Q_2Q2​, and the ends are joined in a prescribed twisted pattern (pp. 107–108). A degenerate appearance of K4K_4K4​ is a knot. A strip (A,C,B)(A, C, B)(A,C,B) is a family of paths ("rungs") from AAA to BBB through CCC; an antistrip is a strip of G‾\overline{G}G. A striation LLL is made of m≥2m \ge 2m≥2 strips and n≥2n \ge 2n≥2 antistrips. All rungs and antirungs are odd, the strips are pairwise anticomplete, the antistrips pairwise complete, and every strip is parallel or co-parallel to every antistrip, with enough "twists" between them (p. 112). A striation is maximal if no striation has a strictly larger vertex set. The paper defines when a set of vertices is local for a knot or striation and when it resolves one.

Outcomes. A double split graph has its vertices partitioned into {ai},{bi}\{a_i\}, \{b_i\}{ai​},{bi​} (m≥2m \ge 2m≥2) and {cj},{dj}\{c_j\}, \{d_j\}{cj​},{dj​} (n≥2n \ge 2n≥2). Each aibia_ib_iai​bi​ is an edge and each cjdjc_jd_jcj​dj​ a nonedge, distinct pairs {ai,bi}\{a_i,b_i\}{ai​,bi​} are anticomplete and distinct pairs {cj,dj}\{c_j,d_j\}{cj​,dj​} complete to each other, and every {ai,bi}\{a_i,b_i\}{ai​,bi​} and {cj,dj}\{c_j,d_j\}{cj​,dj​} are joined by exactly two disjoint edges. A skew partition (A,B)(A, B)(A,B) of V(G)V(G)V(G) has G∣AG|AG∣A disconnected and G‾∣B\overline{G}|BG∣B disconnected. It is balanced if no odd path joins nonadjacent vertices of BBB through AAA and no odd antipath joins adjacent vertices of AAA through BBB. A proper 2-join is a partition (X1,X2)(X_1, X_2)(X1​,X2​) of V(G)V(G)V(G) whose only cross edges are complete joins A1A_1A1​–A2A_2A2​ and B1B_1B1​–B2B_2B2​, with the side conditions of p. 53.

Formalization targets

Goal: Theorem 9.6 (p. 116)

Let GGG be Berge, with every appearance of K4K_4K4​ in GGG and in G‾\overline{G}G degenerate and no induced subgraph of GGG isomorphic to L(K3,3)L(K_{3,3})L(K3,3​). Then

G is double split ∨ G admits a balanced skew partition ∨ G or G‾ admits a proper 2-join ∨ K4 appears in neither G nor G‾.G \text{ is double split} \ \lor\ G \text{ admits a balanced skew partition} \ \lor\ G \text{ or } \overline{G} \text{ admits a proper 2-join} \ \lor\ K_4 \text{ appears in neither } G \text{ nor } \overline{G}.G is double split ∨ G admits a balanced skew partition ∨ G or G admits a proper 2-join ∨ K4​ appears in neither G nor G.

Milestones

  1. 9.1 (p. 108): in a knot of a Berge graph all four paths and antipaths are odd, and either both paths or both antipaths have length one.
  2. 9.3 (p. 109): a connected set FFF whose attachments to a knot are not local either contains a vertex whose neighbourhood resolves the knot, or attaches in one of three special ways ("up to symmetry").
  3. 9.4 (p. 112): the neighbourhood in V(L)V(L)V(L) of a vertex outside a maximal striation LLL is local or resolves LLL.
  4. 9.5 (p. 113): if every vertex of a connected set FFF outside V(L)V(L)V(L) has a local neighbourhood, then the attachments of FFF in V(L)V(L)V(L) are local.

9.3–9.5 assume that no K4K_4K4​-enlargement appears in GGG or G‾\overline{G}G and that no appearance of K4K_4K4​ in GGG or G‾\overline{G}G is overshadowed. An optional, non-milestone item poses 9.7 (p. 118): a Berge graph with an appearance of K4K_4K4​ is a line graph or the complement of one, a double split graph, or admits a proper 2-join (in GGG or G‾\overline{G}G) or a balanced skew partition.

Significance

9.6 is the step of the proof that produces double split graphs, one of the five basic classes. In the main argument it follows step 1.8.1 (5.1, nondegenerate appearances of K4K_4K4​) and is combined with it in 9.7. Through 9.7 it gives the first half of 13.5: every recalcitrant graph belongs to the class F5\mathcal{F}_5F5​ and so contains no appearance of K4K_4K4​ in GGG or G‾\overline{G}G.

The theorem has been proved since 2006. We know of no machine-checked proof of the strong perfect graph theorem or of any of its steps. Mathlib has line graphs and graph embeddings but no subdivisions, appearances or decompositions of Berge graphs. This mission produces a faithful Lean statement of step 1.8.3 and of the four lemmas its proof rests on, together with Lean definitions of knots, strips and striations.

Difficulty

The obvious approach would take a degenerate appearance of K4K_4K4​ and study how each remaining vertex attaches to it, as §§5–6 do for nondegenerate appearances. This does not close. A degenerate appearance can be read as a line graph or as its complement, so the analysis of a vertex in GGG and in G‾\overline{G}G has to be run at once. Single vertices can also be absorbed into larger structures that the line-graph analysis does not see. The proof grows the appearance to a maximal striation and classifies attachments to it, and the hard steps are 9.4 and 9.5. Their proofs need 9.3 in every case, and they use maximality to refute configurations that would let the striation grow.

Formalization scope

Graphs are SimpleGraph V on a Fintype with decidable equality. Every object is defined as on the page, in the namespace StrongPerfectGraph.DoubleSplit.

  • Paths and antipaths are induced and given as vertex lists. A list fixes the labelling of the ends: ai,bia_i, b_iai​,bi​ are the first and last vertices of PiP_iPi​, and xj,yjx_j, y_jxj​,yj​ those of QjQ_jQj​.
  • The empty set is connected (p. 54).
  • Line graphs are Mathlib's SimpleGraph.lineGraph; "isomorphic to an induced subgraph" is an induced embedding ↪g. Subdivisions HHH and enlargements J′J'J′ range over graphs on Fin k.
  • "Every appearance of K4K_4K4​ is degenerate" is the absence of a nondegenerate appearance. The hypothesis on L(K3,3)L(K_{3,3})L(K3,3​) concerns GGG only. Every other hypothesis concerns both GGG and G‾\overline{G}G, as on the page.
  • A striation is a structure with mmm strips and nnn antistrips indexed by Fin m, Fin n.
  • "Up to symmetry" in 9.3 is the paper's exchange of P1,P2P_1, P_2P1​,P2​ and Q1,Q2Q_1, Q_2Q1​,Q2​ with the ends renamed so that the result is again a knot. Both compatible renamings are allowed, and so is their composite, the reversal of all four.

Dropping the bars of the complement would make the theorem false or vacuous. So would reading "path" as a non-induced path, or encoding a decomposition so that it always exists. The statements use the complement explicitly, induced paths throughout, and the full definitions of p. 53–54.

The proof of 9.6 also uses results of the same paper that are posed in other missions of this series. These are 2.1, 2.2, 4.1 and 4.2, posed in mission II (skew partitions), and 5.3, 5.8, 6.1 and 7.5, posed in mission III (line graphs). They are not posed again here. The bridge 9.2 between knots and line graphs, whose proof the paper omits as obvious, is not posed either. Proofs of 9.1, of 9.3, and reusable lemmas about knots and striations are welcome.

Selected references

  • M. Chudnovsky, N. Robertson, P. Seymour, R. Thomas, The strong perfect graph theorem, Annals of Mathematics 164 (2006), 51–229. https://doi.org/10.4007/annals.2006.164.51
  • C. Berge, Färbung von Graphen, deren sämtliche bzw. deren ungerade Kreise starr sind, Wiss. Z. Martin-Luther-Univ. Halle-Wittenberg Math.-Natur. Reihe 10 (1961), 114.
  • M. Grötschel, L. Lovász, A. Schrijver, Geometric Algorithms and Combinatorial Optimization, Springer, 1988. https://doi.org/10.1007/978-3-642-97881-4
19 thms0 active usersReviewed
Operations Research·Captain: mikedeng1

On the Abstract Properties of Linear Dependence 1: The Rank Postulates and the Independence Postulates Are EquivalentResearch Paper

Motivation

In 1935 Hassler Whitney asked which properties of linear dependence among the columns of a matrix can be stated without reference to the matrix at all. His answer, On the Abstract Properties of Linear Dependence (American Journal of Mathematics 57, 1935), introduced the matroid: a finite set together with a rank function obeying three short postulates. The notion now underlies combinatorial optimization (the greedy algorithm is optimal exactly on matroids, and matroid intersection and partition generalize bipartite matching and arborescence packing), graph theory (graphic and cographic matroids), coding theory and the study of linear representations over finite fields.

A defining feature of the subject is that the same structure can be axiomatized in several apparently unrelated ways: by rank, by independent sets, by bases, by circuits. Each axiom system is convenient for different arguments, and passing between them, a so-called cryptomorphism, is routine in practice. Whitney's paper is where these equivalences first appear. Part I, §§2–4 and §6 (pp. 510–514), derives the basic properties of rank from the rank postulates, deduces from them the postulates for independent sets, and shows that the two systems are equivalent. This mission formalizes that first equivalence.

Setting

Let MMM be a finite set of elements e1,…,ene_1, \dots, e_ne1​,…,en​. Following Whitney, write N+eN + eN+e for N∪{e}N \cup \{e\}N∪{e}, M1+M2M_1 + M_2M1​+M2​ for the union and M1M2M_1 M_2M1​M2​ for the intersection of subsets; ρ(N)\rho(N)ρ(N) is the number of elements of NNN.

A rank system is a function rrr on the subsets of MMM satisfying

  • (R₁) r(∅)=0r(\emptyset) = 0r(∅)=0;
  • (R₂) for every subset NNN and element e∉Ne \notin Ne∈/N, r(N+e)=r(N)r(N + e) = r(N)r(N+e)=r(N) or r(N+e)=r(N)+1r(N + e) = r(N) + 1r(N+e)=r(N)+1;
  • (R₃) for every subset NNN and elements e1,e2∉Ne_1, e_2 \notin Ne1​,e2​∈/N, if r(N+e1)=r(N+e2)=r(N)r(N + e_1) = r(N + e_2) = r(N)r(N+e1​)=r(N+e2​)=r(N) then r(N+e1+e2)=r(N)r(N + e_1 + e_2) = r(N)r(N+e1​+e2​)=r(N).

The nullity of NNN is n(N)=ρ(N)−r(N)n(N) = \rho(N) - r(N)n(N)=ρ(N)−r(N), and NNN is independent when n(N)=0n(N) = 0n(N)=0. The increment of (3.1) is Δ(M′,N)=r(M′+N)−r(M′)\Delta(M', N) = r(M' + N) - r(M')Δ(M′,N)=r(M′+N)−r(M′), written Δ(M′,e)\Delta(M', e)Δ(M′,e) when N={e}N = \{e\}N={e}.

An independence system is a predicate "independent" on the subsets of MMM satisfying

  • (I₁) any subset of an independent set is independent;
  • (I₂) if NNN and N′N'N′ are independent and N′N'N′ has exactly one element more than NNN, then N+e′N + e'N+e′ is independent for some e′∈N′e' \in N'e′∈N′ with e′∉Ne' \notin Ne′∈/N.

From an independence system one recovers a rank by letting r(N)r(N)r(N) be the number of elements in a largest independent subset of NNN. In the Lean development these objects are IsRankSystem, nullity, Delta, indepOfRank, IsIndepSystem and rankOfIndep, all in the namespace WhitneyMatroid.RankIndep.

Formalization targets

Goal: (R) and (I) are equivalent (§6, p. 514)

  1. If rrr satisfies (R₁)–(R₃), then {N:ρ(N)=r(N)}\{N : \rho(N) = r(N)\}{N:ρ(N)=r(N)} satisfies (I₁), (I₂), contains ∅\emptyset∅, and
r(N)=max⁡{ρ(I):I⊆N, ρ(I)=r(I)}for every N.r(N) = \max\{\rho(I) : I \subseteq N,\ \rho(I) = r(I)\} \quad \text{for every } N.r(N)=max{ρ(I):I⊆N, ρ(I)=r(I)}for every N.
  1. If "independent" satisfies (I₁), (I₂) and ∅\emptyset∅ is independent, then r(N)=max⁡{ρ(I):I⊆N independent}r(N) = \max\{\rho(I) : I \subseteq N \text{ independent}\}r(N)=max{ρ(I):I⊆N independent} satisfies (R₁)–(R₃), and NNN is independent if and only if ρ(N)=r(N)\rho(N) = r(N)ρ(N)=r(N).

Both translations and both round trips are part of the goal: Whitney's conclusion is not only that each system implies the other but that "the definitions of the rank and the independence or dependence of any subset of MMM agree under the two systems".

Milestones

  • Lemma 1 (p. 510): r(N)≥0r(N) \ge 0r(N)≥0, n(N)≥0n(N) \ge 0n(N)≥0, and N⊆M′N \subseteq M'N⊆M′ implies r(N)≤r(M′)r(N) \le r(M')r(N)≤r(M′), n(N)≤n(M′)n(N) \le n(M')n(N)≤n(M′).
  • Lemma 2 (p. 510): any subset of an independent set is independent, which is (I₁).
  • Lemma 3 (p. 511): Δ(M+e2,e1)≤Δ(M,e1)\Delta(M + e_2, e_1) \le \Delta(M, e_1)Δ(M+e2​,e1​)≤Δ(M,e1​).
  • Lemma 4 (p. 511): Δ(M+N,e)≤Δ(M,e)\Delta(M + N, e) \le \Delta(M, e)Δ(M+N,e)≤Δ(M,e).
  • Theorem 3 (p. 511): Δ(M+N2,N1)≤Δ(M,N1)\Delta(M + N_2, N_1) \le \Delta(M, N_1)Δ(M+N2​,N1​)≤Δ(M,N1​); equivalently
r(M+N1+N2)≤r(M+N1)+r(M+N2)−r(M),r(M1+M2)≤r(M1)+r(M2)−r(M1M2).r(M + N_1 + N_2) \le r(M + N_1) + r(M + N_2) - r(M), \qquad r(M_1 + M_2) \le r(M_1) + r(M_2) - r(M_1 M_2).r(M+N1​+N2​)≤r(M+N1​)+r(M+N2​)−r(M),r(M1​+M2​)≤r(M1​)+r(M2​)−r(M1​M2​).
  • §4 (pp. 511–512): the independent sets of a rank system satisfy (I₂).

Significance

The equivalence makes the rank function and the family of independent sets two descriptions of one object. Every later result of Whitney's paper, and of matroid theory generally, moves between them without comment: the circuit postulates of §5 and §8, the base postulates of §7 and the duality of §§11–13 are all phrased through rank or independence as convenient. Theorem 3 is the submodularity of rank, the property that connects matroids to submodular function minimization and polymatroids; here it is derived from the purely local postulates (R₁)–(R₃), which constrain the rank only under the addition of one or two elements.

On the formal side, Mathlib defines Matroid through independent sets (with constructors from other axiom systems) and proves submodularity of its rank; the platform has submodularity for Mathlib matroids (FamousTheorems.matroid_rank_submodular_7a). Neither starts from Whitney's local rank postulates. What this mission adds is a machine-checked derivation of the global properties of rank from (R₁)–(R₃) and of Whitney's original equivalence, stated for his own postulates, so that the later missions of this series, which work from the same postulates, rest on a verified foundation. The result itself has been settled since 1935; the open work is the formal proof.

Difficulty

The postulates (R₂) and (R₃) are local: they speak about adding at most two elements to a set. Monotonicity and the bound r(N)≤ρ(N)r(N) \le \rho(N)r(N)≤ρ(N) follow by adding elements one at a time, but the submodular inequality relates arbitrary sets, and nothing in (R₃) mentions more than two new elements. The gap between the local and the global statement is the substance of Lemmas 3, 4 and Theorem 3, and the deduction of (I₂) depends on it.

In the converse direction the rank is defined as a maximum over independent subsets, while (I₂) only augments a set from an independent set with exactly one element more; (R₃) for the derived rank is a statement about three sets that are not given in that form. The round trips are where the two halves meet, and each depends on the global properties of the first half rather than on the postulates alone.

Formalization scope

Elements form a type α with [Fintype α] [DecidableEq α]; subsets are Finset α, and the matroid MMM is the whole type. Ranks, nullities and increments take values in ℤ, so differences never truncate; Whitney allows any number, but (R₁) and (R₂) force nonnegative integers. Postulate (I₂) is stated in Whitney's form with N'.card = N.card + 1, not the general augmentation for ρ(N)<ρ(N′)\rho(N) < \rho(N')ρ(N)<ρ(N′). Postulates (R₂), (R₃) keep their hypotheses e,e1,e2∉Ne, e_1, e_2 \notin Ne,e1​,e2​∈/N. Lemmas 3, 4 and Theorem 3 are stated for arbitrary subsets M,N,N1,N2M, N, N_1, N_2M,N,N1​,N2​; Lemma 1's monotonicity for arbitrary N⊆M′N \subseteq M'N⊆M′, the form in which the paper uses it. rankOfIndep is the supremum of cardinalities over the independent members of the powerset.

The paper takes for granted that the empty set is independent in system (I). Without that hypothesis, the predicate declaring nothing independent satisfies (I₁) and (I₂) vacuously, the supremum defining the rank returns 000, and the round trip fails; the goal therefore assumes ∅\emptyset∅ independent in part 2 and proves it in part 1. A formalization in terms of Mathlib's Matroid would make the goal a restatement of library facts, since Mathlib's matroids are independence systems by construction; the goal is deliberately about the postulates as predicates on functions and on families of sets.

A complete development needs only finite set combinatorics (Finset.card, induction on finite sets, Finset.sup). The derived lemmas (monotonicity, submodularity, (I₂)) are reusable for any later work from Whitney's rank postulates, including the circuit-postulate equivalence of the companion mission. A bridge from rank systems to Mathlib's Matroid (via IndepMatroid.ofFinset) would be a welcome addition but is not part of the goal. Proofs of the milestones in any order are welcome; Theorem 3 and the §4 deduction are the natural first targets.

Selected references

  • H. Whitney, On the Abstract Properties of Linear Dependence, American Journal of Mathematics 57 (1935), no. 3, 509–533. https://doi.org/10.2307/2371182
  • J. Oxley, Matroid Theory, 2nd ed., Oxford Graduate Texts in Mathematics 21, Oxford University Press, 2011. https://doi.org/10.1093/acprof:oso/9780198566946.001.0001
  • J. Kung (ed.), A Source Book in Matroid Theory, Birkhäuser, 1986. https://doi.org/10.1007/978-1-4684-9199-9
  • Mathlib, Mathlib.Data.Matroid (matroids via independent sets; IndepMatroid.ofFinset). https://leanprover-community.github.io/mathlib4_docs/Mathlib/Data/Matroid/Basic.html
8 thms0 active usersReviewed
Graph Theory·Captain: mikedeng1

The Strong Perfect Graph Theorem V: A Berge Graph with No Nondegenerate Appearance of K4 Containing an Even Prism Is a Nine-Vertex Even Prism or DecomposesResearch Paper

Motivation

A perfect graph is a finite simple graph in which every induced subgraph has chromatic number equal to its largest clique size. This equality gives a structural reason that a clique lower bound on the number of colours is attainable for every induced part of the graph. Berge proposed a forbidden-subgraph description of perfect graphs: a graph should be perfect exactly when it has neither an odd hole nor an odd antihole. Chudnovsky, Robertson, Seymour and Thomas proved that statement in The strong perfect graph theorem, Annals of Mathematics 164 (2006), Theorem 1.2. Their proof separates possible configurations in a Berge graph and shows that each either belongs to a controlled class or admits a decomposition.

This mission isolates the even-prism step, Theorem 10.6 of that paper. A prism is one of the configurations that can appear in a Berge graph even though odd holes and antiholes do not. The step matters because its conclusion leaves only a specific nine-vertex graph or one of two decompositions that the larger proof handles elsewhere. It is the result identified as step 1.8.4 in the authors’ outline (paper, pp. 59 and 124).

Setting

All graphs here are finite and simple. A path means an induced path; a single vertex is allowed as a path of length zero. A hole is an induced cycle with at least four vertices, and an antihole is a hole in the complementary graph. A graph GGG is Berge if every hole of GGG and of its complement G‾\overline GG has even length.

A prism has two disjoint triangles, A={a1,a2,a3}A=\{a_1,a_2,a_3\}A={a1​,a2​,a3​} and B={b1,b2,b3}B=\{b_1,b_2,b_3\}B={b1​,b2​,b3​}, joined by three pairwise vertex-disjoint induced paths RiR_iRi​ from aia_iai​ to bib_ibi​. Between distinct paths, the only edges are those in AAA and those in BBB. The prism is even when all three RiR_iRi​ have even length. “GGG contains an even prism” means that such paths exist as an induced configuration in GGG; GGG may have other vertices. “GGG is an even prism” means the paths cover every vertex of GGG (paper, pp. 93 and 119).

An appearance of K4K_4K4​ is an induced copy in GGG of the line graph of a bipartite subdivision HHH of the four-vertex complete graph. It is nondegenerate if no four-cycle of HHH contains all four branch vertices. Only appearances in GGG are excluded in this mission; appearances in G‾\overline GG are allowed by the hypotheses (paper, pp. 72, 74–75).

A proper 2-join partitions the vertices into two sides with specified, nonempty attachment sets. The cross edges are exactly the two complete attachment pairs; every component of either side meets both of its attachment sets. If a side is itself a path between singleton attachment sets, that path has odd length at least three. A balanced skew partition divides the vertices into AAA and BBB so that AAA is disconnected, BBB is disconnected in the complement, and two path parity conditions hold: no odd path crosses AAA between nonadjacent vertices of BBB, and no odd antipath crosses BBB between adjacent vertices of AAA (paper, pp. 53–54).

Formalization targets

Prism lemmas

The numbered milestones are Theorems 7.2–7.4 and 10.5. They assert common parity of the three prism paths, common neighbours of an anticonnected set at both end triangles, preservation of two neighbours under replacement of one even prism path, and the balanced skew partition forced by a major vertex. A vertex is major when it is adjacent to at least two vertices of each end triangle. These milestones match the paper’s statements on pp. 93 and 123 (paper).

Goal: Theorem 10.6

For a Berge graph GGG with no nondegenerate appearance of K4K_4K4​ in GGG,

G contains an even prism⟹(G is an even prism and ∣V(G)∣=9)  ∨  G admits a proper 2-join  ∨  G admits a balanced skew partition.G\text{ contains an even prism} \quad\Longrightarrow\quad \bigl(G\text{ is an even prism and }|V(G)|=9\bigr) \;\lor\; G\text{ admits a proper 2-join} \;\lor\; G\text{ admits a balanced skew partition}.G contains an even prism⟹(G is an even prism and ∣V(G)∣=9)∨G admits a proper 2-join∨G admits a balanced skew partition.

The first case describes the entire graph, not just an induced nine-vertex subgraph. The statement fixes all outcomes exactly as in Theorem 10.6, p. 124.

Significance

Theorem 10.6 removes even prisms from the unresolved part of the strong perfect graph theorem’s structural argument. If the graph is larger than the exceptional prism and has no nondegenerate K4K_4K4​ appearance, the theorem supplies a proper 2-join or a balanced skew partition. Subsequent results can work with those decompositions instead of treating arbitrary attachments to a prism (paper, §10 and the outline at 1.8.4).

The mathematical result is proved in the 2006 paper. The work here is to give its graph objects and statements machine-checkable meanings, then formalize the known proof. This proposal contains open Lean theorem statements and sorry-free definitions; the proof obligations remain for solvers. The definitions of induced paths, holes, subdivisions, and decomposition predicates can also support other steps of this paper. The Roussel–Rubio lemma and the balanced-skew-partition results of §§2–4 are proved in the same paper and posed in mission II of this series. The prism-attachment result 10.4 is posed in mission VI, where it is used most directly.

Difficulty

An outside connected set can attach to several parts of a prism without containing a single major vertex. Its attachments need not be local to one path or one triangle, so checking vertices one at a time does not decide which decomposition exists. The paper’s §10 distinguishes several attachment patterns; the evenness of the paths and the exclusion of a nondegenerate K4K_4K4​ appearance restrict them, but do not themselves give a 2-join or skew partition by a one-line parity argument. The larger proof must also account for attachments throughout the graph while preserving the full definitions of both decomposition outcomes (paper, pp. 119–127).

Formalization scope

Lean represents a graph as SimpleGraph V with a finite vertex type. The prism is three lists of vertices, each an induced path; the lists are disjoint, and the cross-edge condition admits exactly the two end triangles. Reversing a list changes its orientation but not the underlying graph configuration. A hole uses a cyclic list with the closing edge; Berge checks holes in both GGG and G‾\overline GG. Connectivity is reachability in an induced graph, so the empty vertex set is connected as the paper says. A K4K_4K4​ subdivision uses six tracks on a finite carrier, with all its vertices and edges accounted for; the appearance uses an induced graph embedding of its line graph. The exception checks that the prism covers the whole graph and that ∣V(G)∣=9|V(G)|=9∣V(G)∣=9.

These encodings require genuinely induced paths and the paper’s nondegenerate appearance condition. Dropping either would change the theorem. The proper 2-join includes component reachability and the odd-path special case; the balanced skew partition includes both path parity clauses. The nine-vertex graph formed by two triangles and three two-edge paths has a separate sorry-free Lean witness, so the exceptional outcome is nonvacuous. Useful contributions include the numbered prism lemmas, attachment analysis for 10.6, and reusable results about finite induced paths and graph subdivisions.

Selected references

  • Maria Chudnovsky, Neil Robertson, Paul Seymour and Robin Thomas, The strong perfect graph theorem, Annals of Mathematics 164 (2006), 51–229. DOI: 10.4007/annals.2006.164.51.
12 thms0 active usersReviewed
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 thms0 active usersReviewed
PreviousPage 4 of 5Next

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