Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Discrete Geometry

156 missions · 79 completed

Missions

Open77Completed79All156
Captain: xuanji

Circle packing in a square: exact constantsTextbook

Motivation

Packing congruent circles into a square is a classical problem in discrete geometry: for each natural number nnn, choose a common radius as large as possible while keeping all disks inside the square and preventing overlap. Every exact value requires two logically distinct achievements: an explicit configuration attaining the proposed radius and a proof that no configuration can do better.

The character of those proofs changes sharply with nnn. The first cases admit short geometric arguments; later cases use contact-graph analysis, specialized case divisions, or computer-assisted global optimization with interval arithmetic. Formalizing the resulting constants therefore provides a growing benchmark for extremal geometry, real algebra, finite configurations, and verified computation in Lean.

This is an open-ended formalization mission. It begins with the exact constants currently represented by theorem-backed milestones, but it is not restricted to a fixed terminal value of nnn. Further milestones may be added whenever an exact packing value and its rigorous optimality argument are identified and stated precisely enough for formalization.

Setting

A point is a pair of real coordinates. For points p=(x,y)p=(x,y)p=(x,y) and q=(x′,y′)q=(x',y')q=(x′,y′), squared Euclidean distance is

sqDist⁡(p,q)=(x−x′)2+(y−y′)2.\operatorname{sqDist}(p,q)=(x-x')^2+(y-y')^2.sqDist(p,q)=(x−x′)2+(y−y′)2.

For a real radius rrr, a point lies in the inner square when both coordinates belong to the closed interval [r,1−r][r,1-r][r,1−r]. This is exactly the coordinate condition saying that a closed disk of radius rrr, centered at that point, is contained in the unit square.

The predicate Packable⁡(n,r)\operatorname{Packable}(n,r)Packable(n,r) requires 0≤r≤120\le r\le \tfrac120≤r≤21​ and a family of nnn centers in the inner square such that the squared distance between every two distinctly indexed centers is at least (2r)2(2r)^2(2r)2. Equality is allowed, so tangent disks are admitted. Radius zero is also admitted.

Define

rn=sup⁡{r∈R:Packable⁡(n,r)}r_n=\sup\{r\in\mathbb R:\operatorname{Packable}(n,r)\}rn​=sup{r∈R:Packable(n,r)}

and define the optimal covered-area fraction by

cn=nπrn2.c_n=n\pi r_n^2.cn​=nπrn2​.

It is often convenient to use the equivalent point-separation constant dnd_ndn​, the greatest possible minimum pairwise distance among nnn points in the unit square. The conversion is

rn=dn2(1+dn),cn=nπ(dn2(1+dn))2.r_n=\frac{d_n}{2(1+d_n)}, \qquad c_n=n\pi\left(\frac{d_n}{2(1+d_n)}\right)^2.rn​=2(1+dn​)dn​​,cn​=nπ(2(1+dn​)dn​​)2.

The Lean definitions use a supremum rather than assuming in advance that an optimal packing is attained.

Current exact-value milestones

The mission currently contains theorem-backed milestones for the following values:

nnnExact separation or area valueProof character in the supplied notes
222d2=2d_2=\sqrt2d2​=2​, hence c2=π(3−22)c_2=\pi(3-2\sqrt2)c2​=π(3−22​)diagonal bound
333d3=6−2d_3=\sqrt6-\sqrt2d3​=6​−2​minimum enclosing square of a triangle
444d4=1d_4=1d4​=1, hence c4=π/4c_4=\pi/4c4​=π/4convex hull and perimeter
555d5=1/2d_5=1/\sqrt2d5​=1/2​four-cell pigeonhole argument
666d6=13/6d_6=\sqrt{13}/6d6​=13​/6case-specific geometric proof
777d7=4−23d_7=4-2\sqrt3d7​=4−23​hand proof and later computer verification
888d8=2−3d_8=\sqrt{2-\sqrt3}d8​=2−3​​case-specific geometric proof
999d9=1/2d_9=1/2d9​=1/2, hence c9=π/4c_9=\pi/4c9​=π/4classical geometric proof
161616d16=1/3d_{16}=1/3d16​=1/3, hence c16=π/4c_{16}=\pi/4c16​=π/4theoretical grid-optimality proof
252525d25=1/4d_{25}=1/4d25​=1/4, hence c25=π/4c_{25}=\pi/4c25​=π/4theoretical grid-optimality proof
363636d36=1/5d_{36}=1/5d36​=1/5, hence c36=π/4c_{36}=\pi/4c36​=π/4theoretical grid-optimality proof

For rows stated using dnd_ndn​, the corresponding milestone for cnc_ncn​ uses the conversion formula above. The equalities are claims about the supremum-defined packing constants, not merely about the displayed candidate configurations.

An extensible mission

The milestone list is intended to grow. The supplied survey notes classify n=2,…,33n=2,\ldots,33n=2,…,33 and n=36n=36n=36 as rigorously solved in the cited literature, while distinguishing n=34n=34n=34 and n=35n=35n=35 as not rigorously closed in the cited 2021 account. Many of the computer-assisted cases do not have a simple radical expression in the supplied notes. Before such a case is linked to a Lean theorem, its primary source must provide a precise candidate value, algebraic characterization, certified enclosure, or optimal-configuration certificate that can be stated faithfully.

A new milestone should identify:

  1. the precise value or exact characterization being formalized;
  2. an attaining configuration or a certified existence argument;
  3. a universal upper bound or global-optimality certificate;
  4. the primary source and exact theorem, equation, or certificate location;
  5. any trusted computational artifact and the arithmetic guarantees it requires.

Numerical evidence and strong bounds are valuable, but they must be labeled as bounds rather than exact-value milestones. Conversely, newly published exact results for larger nnn may be added without changing the underlying definitions.

Proof obligations

Every exact-value milestone must connect the proposed value to Packable, r_n, and c_n. Constructing a configuration establishes only a lower bound. An upper-bound argument without attainability also does not establish equality. A complete proof must bridge both directions through the supremum definition.

The proof methods may include:

  • elementary diameter, pigeonhole, convexity, or enclosing-shape arguments;
  • normalization between disk centers and point-separation configurations;
  • contact-graph and boundary-constraint analysis;
  • finite case decompositions;
  • interval arithmetic and formally checked branch-and-bound certificates;
  • exact algebraic identities needed to convert dnd_ndn​ into rnr_nrn​ and cnc_ncn​.

Shortcuts that redefine rnr_nrn​, dnd_ndn​, or cnc_ncn​ to equal a desired answer are excluded. The constants must remain consequences of the common geometric model.

Mission structure

The root theorem CirclePackingConstants.c_all is the conjunction of the eleven exact-value milestones currently in the mission, covering n=2,3,4,5,6,7,8,9,16,25,36n=2,3,4,5,6,7,8,9,16,25,36n=2,3,4,5,6,7,8,9,16,25,36. Its proof sketch reduces the root directly to those milestone theorems, so the mission remains open until every current exact value is proved.

The milestone theorems remain separately reusable and independently auditable. When further exact values are added, a successor aggregate theorem can extend the conjunction and become the new root without replacing the shared definitions or invalidating earlier results.

This structure allows elementary cases, historical hand proofs, and computer-assisted certificates to progress independently while remaining part of one cumulative library of exact circle-packing constants.

Formalization scope

The Lean model uses ℝ × ℝ for points and an explicit coordinate formula for squared Euclidean distance. Disk containment is represented by inclusive coordinate inequalities. Nonoverlap is represented by a weak squared-distance inequality, so tangency is permitted. The indexing type is Fin n, and the definitions apply to every natural number, including zero.

The definition bundle contains only Point, sqDist, InInnerSquare, Packable, r_n, and c_n. Solvers may introduce normalization maps, separation bounds, explicit configurations, supremum lemmas, contact structures, certificate checkers, and radical or polynomial identities as auxiliary declarations.

Selected references

  • User-supplied notes, Circles in squares: constants, proofs, and what is actually known, supplied September 12, 2026. The notes summarize the exact small-nnn formulas, grid cases, historical proof taxonomy, and computer-assisted frontier used to organize this mission.
  • J. Schaer and A. Meir, “On a geometric extremum problem,” Canadian Mathematical Bulletin 8 (1965), 21–27.
  • J. Schaer, “The densest packing of nine circles in a square,” Canadian Mathematical Bulletin 8 (1965), 273–277.
  • B. L. Schwartz, “Separating points in a square,” Journal of Recreational Mathematics 3 (1970), 195–204.
  • J. B. M. Melissen, “Densest packing of six equal circles in a square,” Elemente der Mathematik 49 (1994), 27–31.
  • M. C. Markot, “Improved interval methods for solving circle packing problems in the unit square,” Journal of Global Optimization 81 (2021), 773–803.
  • Erich Friedman, Circles in Squares, Erich's Packing Center, for background tables and diagrams of candidate packings.
88 thms8 active users
CombinatoricsLinear 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 thms4 active usersReviewed
Computational GeometryGraph Theory·Captain: hao jia

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

Motivation

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

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

Setting

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

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

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

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

Formalization targets

The two-part OPG record

The source records both

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

and

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

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

Published first part

A milestone formalizes the stronger published statement

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

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

Universal bound

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

Significance

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

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

Difficulty

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

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

Formalization scope

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

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

Selected references

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

An Improved Kalai-Kleitman Bound for the Diameter of a Polyhedron: Every d-Polyhedron with n Facets, Bounded or Not, Has Graph Diameter at Most (n-d)^(log2 d)Research Paper

Motivation

The simplex method for linear programming moves from vertex to vertex along the edges of a polyhedron. A lower bound on the number of pivots of any edge-following method is the graph diameter of the feasible polyhedron: the largest number of edges needed to connect two of its vertices. Whether this diameter is bounded by a polynomial in the dimension ddd and the number nnn of facets is the polynomial Hirsch conjecture, open since the 1950s. The original Hirsch conjecture, n−dn-dn−d, was refuted for unbounded polyhedra by Klee and Walkup (1967) and for bounded ones by Santos (2012) (doi:10.4007/annals.2012.176.1.7).

Timeline.

  • 1964–1966. Klee proves that a three-dimensional polyhedron with nnn facets, bounded or not, has diameter at most n−3n-3n−3 (doi:10.4153/CJM-1964-058-3, doi:10.1137/0113061, doi:10.2140/pjm.1966.17.249).
  • 1967. Klee and Walkup show that the Hirsch bound fails for unbounded polyhedra in dimension four (doi:10.1007/BF02395039).
  • 1992. Kalai and Kleitman prove the first subexponential bound, nlog⁡2d+2n^{\log_2 d+2}nlog2​d+2, for all polyhedra (doi:10.1090/S0273-0979-1992-00285-9, arXiv:math/9204233).
  • 2010. Eisenbrand, Hähnle, Razborov and Rothvoss show that nlog⁡2d+1n^{\log_2 d+1}nlog2​d+1 holds for a broad combinatorial abstraction of polyhedra (doi:10.1287/moor.1100.0470).
  • 2014. Todd sharpens the Kalai–Kleitman bound to (n−d)log⁡2d(n-d)^{\log_2 d}(n−d)log2​d for all polyhedra, bounded or not (arXiv:1402.3579; SIAM J. Discrete Math. 28 (2014), doi:10.1137/140962310). This is the result of this mission.

Setting

A polyhedron P⊆RdP\subseteq\mathbb R^dP⊆Rd is the set of solutions of finitely many linear inequalities,

P=Hpoly⁡(a,b)={x∈Rd: ⟨ai,x⟩≤bi, i=1,…,n}.P=\operatorname{Hpoly}(a,b)=\{x\in\mathbb R^d:\ \langle a_i,x\rangle\le b_i,\ i=1,\dots,n\}.P=Hpoly(a,b)={x∈Rd: ⟨ai​,x⟩≤bi​, i=1,…,n}.

It need not be bounded. A vertex of PPP is a point vvv for which some half-space HHH satisfies P∩H={v}P\cap H=\{v\}P∩H={v}, equivalently an extreme point of PPP. Two vertices v≠wv\neq wv=w are adjacent when P∩H=[v,w]P\cap H=[v,w]P∩H=[v,w] for some half-space HHH; the segment [v,w][v,w][v,w] is then an edge. A path of length kkk is a sequence of vertices v=v0,…,vk=wv=v_0,\dots,v_k=wv=v0​,…,vk​=w with consecutive vertices adjacent, ρP(v,w)\rho_P(v,w)ρP​(v,w) is the least such kkk, and the diameter is δ(P)=max⁡ρP(v,w)\delta(P)=\max\rho_P(v,w)δ(P)=maxρP​(v,w) over pairs of vertices. A (d,n)(d,n)(d,n)-polyhedron is a ddd-dimensional polyhedron with nnn facets, and

Δ(d,n)=max⁡{δ(P): P a (d,n)-polyhedron},\Delta(d,n)=\max\{\delta(P):\ P \text{ a } (d,n)\text{-polyhedron}\},Δ(d,n)=max{δ(P): P a (d,n)-polyhedron},

the maximum taken over bounded and unbounded polyhedra alike. All logarithms are to base 2.

Formalization targets

Goal: Theorem 1

Δ(d,n)≤(n−d)log⁡d(1≤d≤n),Δ(1,1)=0.\Delta(d,n)\le (n-d)^{\log d}\quad(1\le d\le n),\qquad \Delta(1,1)=0.Δ(d,n)≤(n−d)logd(1≤d≤n),Δ(1,1)=0.

The goal fixes no constant beyond the paper's own exponent and leaves no case of ddd, nnn or boundedness open.

Milestones, in the order the proof uses them

  1. The identity (n−d)log⁡d=dlog⁡(n−d)(n-d)^{\log d}=d^{\log(n-d)}(n−d)logd=dlog(n−d) (p. 2).
  2. Two vertices on a common facet are at distance at most Δ(d−1,n−1)\Delta(d-1,n-1)Δ(d−1,n−1) (p. 2, in the proof of Lemma 1).
  3. Lemma 1 (Kalai–Kleitman): for 2≤d≤⌊n/2⌋2\le d\le\lfloor n/2\rfloor2≤d≤⌊n/2⌋, Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,⌊n/2⌋)+2\Delta(d,n)\le\Delta(d-1,n-1)+2\Delta(d,\lfloor n/2\rfloor)+2Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,⌊n/2⌋)+2.
  4. Δ(d,d)=0\Delta(d,d)=0Δ(d,d)=0.
  5. Δ(1,n)≤1\Delta(1,n)\le 1Δ(1,n)≤1.
  6. Δ(2,n)≤n−2\Delta(2,n)\le n-2Δ(2,n)≤n−2.
  7. Δ(3,n)≤n−3\Delta(3,n)\le n-3Δ(3,n)≤n−3 (Klee), for unbounded polyhedra too.
  8. If n<2dn<2dn<2d, then Δ(d,n)≤Δ(d−1,n−1)\Delta(d,n)\le\Delta(d-1,n-1)Δ(d,n)≤Δ(d−1,n−1).
  9. In the ten cases d=4, 8≤n≤11d=4,\ 8\le n\le 11d=4, 8≤n≤11; d=5, 10≤n≤12d=5,\ 10\le n\le 12d=5, 10≤n≤12; d=6, 12≤n≤13d=6,\ 12\le n\le 13d=6, 12≤n≤13; d=7, n=14d=7,\ n=14d=7, n=14, the bound ⌊(n−d)log⁡d⌋\lfloor (n-d)^{\log d}\rfloor⌊(n−d)logd⌋ satisfies the recursion of Lemma 1.
  10. For d≥4d\ge4d≥4, n−d≥8n-d\ge8n−d≥8, n≥2dn\ge 2dn≥2d: (d−1)log⁡(n−d)+2(n/2−d)log⁡d+2≤dlog⁡(n−d)(d-1)^{\log(n-d)}+2(n/2-d)^{\log d}+2\le d^{\log(n-d)}(d−1)log(n−d)+2(n/2−d)logd+2≤dlog(n−d). For n>2dn>2dn>2d this is the displayed inequality; at n=2dn=2dn=2d the second power is zero, consistent with Δ(d,d)=0\Delta(d,d)=0Δ(d,d)=0.

Significance

The result. The bound (n−d)log⁡d(n-d)^{\log d}(n−d)logd is never larger than the Kalai–Kleitman bound. It equals the true value n−dn-dn−d in dimensions 1 and 2, and it depends on nnn only through the Hirsch quantity n−dn-dn−d. Since (n−d)log⁡d=dlog⁡(n−d)(n-d)^{\log d}=d^{\log(n-d)}(n−d)logd=dlog(n−d), it is also symmetric under linear-programming duality, which exchanges the roles of ddd and n−dn-dn−d. For fixed n−dn-dn−d it is polynomial in ddd. It applies to the unbounded feasible regions that arise in practice in linear programming.

Formalizing it. The theorem has been proved since 2014. The platform has a proof of the bounded case only: Hirsch.todd_bound assumes d≥3d\ge3d≥3, a nonempty interior and Bornology.IsBounded. This mission's goal extends it to all polyhedra, bounded or not, and to d=1,2d=1,2d=1,2. That is the statement the paper proves, and the one an analysis of the simplex method on unbounded feasible regions needs. The unbounded form of Lemma 1 and the common-facet bound are open on the platform as items of the Kalai–Kleitman series and are reused here. Klee's unbounded three-dimensional bound has no formal proof on the platform: the published Hirsch.klee_three_dimensional_bound and Hirsch.dimension_three_bound cover bounded polyhedra only. The bounded recursion Hirsch.kalai_kleitman_recursion and the bounded facet reduction Hirsch.facet_reduction are proved and may serve as models.

Difficulty

Two steps are hard. The first is the unbounded case. Statements proved for polytopes do not transfer directly, because an unbounded polyhedron may have no vertices or a single one, unbounded edges (rays) that are not edges of the graph, and redundant or implicit-equality rows. Each geometric milestone must hold for these. The second is Klee's three-dimensional bound, which the paper cites without proof. Assuming boundedness to apply the published polytope versions is ruled out: Klee and Walkup showed that unbounded and bounded diameters differ from dimension four on, so the unbounded case carries content of its own. The arithmetic milestones are elementary but sensitive to edge values of log⁡\loglog and the real power.

Formalization scope

All statements use the published definitions Hirsch.Hpoly, Hirsch.Adj and Hirsch.DiamLE (Hirsch_model). Vertices are Set.extremePoints ℝ. DiamLE P B says that any two extreme points are joined by a walk of exactly BBB steps, each step staying put or crossing an edge. So the bound includes connectivity of the vertex–edge graph.

  • Inequality reading. "Δ(d,n)≤B\Delta(d,n)\le BΔ(d,n)≤B" is stated as: for every a:Fin n→Rda:\mathrm{Fin}\,n\to\mathbb R^da:Finn→Rd and bbb, DiamLE (Hpoly a b) B. This ranges over all polyhedra given by nnn inequalities in Rd\mathbb R^dRd, including lower-dimensional, empty and redundant ones. It is equivalent to the paper's statement: a lower-dimensional polyhedron of dimension d′d'd′ with n′n'n′ facets has n′−d′≤n−dn'-d'\le n-dn′−d′≤n−d.
  • No boundedness, nonemptiness or full-dimensionality is assumed anywhere.
  • The logarithm is Real.logb 2 and the power is Real.rpow. A distance is at most (n−d)log⁡d(n-d)^{\log d}(n−d)logd exactly when it is at most the floor ⌊(n−d)log⁡d⌋\lfloor (n-d)^{\log d}\rfloor⌊(n−d)logd⌋, abbreviated toddBound d n in the milestone on the ten cases. The clause Δ(1,1)=0\Delta(1,1)=0Δ(1,1)=0 is a separate conjunct, because 00=10^0=100=1 in Real.rpow.
  • Disclosed additions: nonzero rows in the n<2dn<2dn<2d step (as in the referenced Lemma 1 and the d=2d=2d=2 bound). In the general step, n≥2dn\ge2dn≥2d: the paper secures it by treating n<2dn<2dn<2d first. The second power is written with base n/2−dn/2-dn/2−d; it equals the printed expression when n>2dn>2dn>2d and takes the intended value zero when n=2dn=2dn=2d, where the printed logarithm is undefined.
  • Natural-number subtraction appears in n−3n-3n−3 and n−2n-2n−2 (both harmless, since for small nnn there are no vertices) and in Rd−1\mathbb R^{d-1}Rd−1 for facets.

Trivializing encodings are excluded: no Bornology.IsBounded hypothesis; no definition of Δ\DeltaΔ by sSup over N\mathbb NN, which would be 000 on unbounded families; no use of Hirsch.gdist, which is 000 when no walk exists; and no dropping of the Δ(1,1)=0\Delta(1,1)=0Δ(1,1)=0 clause.

Contributions are welcome on every milestone. A complete development needs facial structure of unbounded polyhedra (facets as polyhedra one dimension lower, tight rows at vertices), the relaxation argument of Lemma 1, Klee's three-dimensional argument, and real-power and logarithm estimates. The facet machinery and Klee's bound are reusable for any diameter result on unbounded polyhedra.

Selected references

  • M. J. Todd, An improved Kalai–Kleitman bound for the diameter of a polyhedron, SIAM J. Discrete Math. 28 (2014); arXiv:1402.3579v2. arXiv:1402.3579, doi:10.1137/140962310
  • 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. doi:10.1090/S0273-0979-1992-00285-9
  • V. Klee, Diameters of polyhedral graphs, Canad. J. Math. 16 (1964), 602–614. doi:10.4153/CJM-1964-058-3
  • V. Klee, Paths on polyhedra I, J. SIAM 13 (1965), 946–956. doi:10.1137/0113061
  • V. Klee, Paths on polyhedra II, Pacific J. Math. 17 (1966), 249–262. doi:10.2140/pjm.1966.17.249
  • V. Klee and D. W. Walkup, The d-step conjecture for polyhedra of dimension d < 6, Acta Math. 117 (1967), 53–78. doi:10.1007/BF02395039
  • F. Eisenbrand, N. Hähnle, A. Razborov, T. Rothvoss, Diameter of polyhedra: limits of abstraction, Math. Oper. Res. 35 (2010), 786–794. doi:10.1287/moor.1100.0470
  • F. Santos, A counterexample to the Hirsch conjecture, Ann. of Math. 176 (2012), 383–412. doi:10.4007/annals.2012.176.1.7
13 thms3 active usersReviewed
Mathematical Physics·Captain: Lucas

Loop Quantum Gravity (Rovelli 2011): Minkowski's theorem and the area gapResearch Paper

Motivation

Loop quantum gravity (LQG) is a proposal for a quantum theory of the gravitational field in which space itself is built from discrete quanta. Its kinematics rests on two mathematical facts that C. Rovelli's review Loop quantum gravity: the first 25 years (Class. Quantum Grav. 28 (2011) 153002, arXiv:1012.4707) singles out. First, the spin–geometry theorem of Penrose, in the extended form based on a theorem of Minkowski (1897): the angular-momentum vectors attached to the links around a node of a spin network, once they satisfy the closure (gauge-invariance) relation, are exactly the data of a flat polyhedron, so every quantum of space carries a geometry (§2.1). Second, the discreteness of geometric spectra: areas take values 8πγℏG∑njn(jn+1)8\pi\gamma\hbar G\sum_n\sqrt{j_n(j_n+1)}8πγℏG∑n​jn​(jn​+1)​ with half-integer spins jnj_njn​, with a strictly positive area gap (§2.1, §2.3, §4.3), which the review calls "the most solid result in loop gravity".

Timeline. Minkowski proved the existence and uniqueness of a convex polyhedron with prescribed face normals and areas in 1897. Penrose introduced spin networks and the spin–geometry theorem around 1971. The discreteness of area and volume in LQG was found by Rovelli and Smolin (1995), with the full area spectrum computed by Ashtekar–Lewandowski and others in 1995–1997.

Setting

Work in Euclidean three-space E3=R3\mathbb E^3=\mathbb R^3E3=R3 with its standard inner product ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩. Given finitely many vectors u1,…,uv∈E3u_1,\dots,u_v\in\mathbb E^3u1​,…,uv​∈E3 and real numbers h1,…,hvh_1,\dots,h_vh1​,…,hv​, the half-space polyhedron is

P(u,h)={x∈E3: ⟨ul,x⟩≤hl for all l},P(u,h)=\{x\in\mathbb E^3:\ \langle u_l,x\rangle\le h_l\ \text{for all } l\},P(u,h)={x∈E3: ⟨ul​,x⟩≤hl​ for all l},

and its lll-th face is Fl=P(u,h)∩{x:⟨ul,x⟩=hl}F_l=P(u,h)\cap\{x:\langle u_l,x\rangle=h_l\}Fl​=P(u,h)∩{x:⟨ul​,x⟩=hl​}. The area of a set lying in a plane orthogonal to a nonzero vector uuu is its two-dimensional Lebesgue measure, computed as the measure of its orthogonal projection onto the plane u⊥u^\perpu⊥. For a nonzero vector nnn, n^=n/∥n∥\hat n=n/\|n\|n^=n/∥n∥ is its direction.

On the quantum side, a node of a spin network with vvv links carries vectors L⃗1,…,L⃗v\vec L_1,\dots,\vec L_vL1​,…,Lv​ (the left-invariant derivative operators), subject to the closure relation ∑lL⃗l=0\sum_l \vec L_l=0∑l​Ll​=0 (eq. (6)). The review's classical content is: the L⃗l\vec L_lLl​ behave as the area-weighted normals of a polyhedron, whose face areas are ∣L⃗l∣|\vec L_l|∣Ll​∣ and whose volume (for v=4v=4v=4) is given by the triple product (eq. (10)). Spins j∈{0,12,1,32,… }j\in\{0,\tfrac12,1,\tfrac32,\dots\}j∈{0,21​,1,23​,…} label links; the quantity j(j+1)\sqrt{j(j+1)}j(j+1)​ is the square root of the SU(2)SU(2)SU(2) Casimir. The area spectrum (eq. (26)) is the set of numbers 8πγℏG∑njn(jn+1)8\pi\gamma\hbar G\sum_n\sqrt{j_n(j_n+1)}8πγℏG∑n​jn​(jn​+1)​ over finite families of spins, with γ\gammaγ the Barbero–Immirzi parameter, ℏ\hbarℏ the reduced Planck constant and GGG Newton's constant.

Formalization targets

Goal: Minkowski's theorem in the form of the spin–geometry theorem

For nonzero n1,…,nv∈E3n_1,\dots,n_v\in\mathbb E^3n1​,…,nv​∈E3 with pairwise distinct directions n^l\hat n_ln^l​, spanning E3\mathbb E^3E3, and with ∑lnl=0\sum_l n_l=0∑l​nl​=0:

∃ h: P(n^,h) bounded, with nonempty interior, and area⁡(Fl)=∥nl∥  ∀l,\exists\,h:\ P(\hat n,h)\ \text{bounded, with nonempty interior, and}\ \operatorname{area}(F_l)=\|n_l\|\ \ \forall l,∃h: P(n^,h) bounded, with nonempty interior, and area(Fl​)=∥nl​∥  ∀l,

and any two such polyhedra are translates of each other.

Milestones (in the order of the review)

  1. Eq. (6), closure: for a bounded polyhedron with nonempty interior and distinct unit normals ulu_lul​, ∑larea⁡(Fl) ul=0\sum_l \operatorname{area}(F_l)\,u_l=0∑l​area(Fl​)ul​=0.
  2. Eq. (10), tetrahedron volume: with face area vectors L1,L2,L3L_1,L_2,L_3L1​,L2​,L3​ of three faces through a vertex, V2=29 ∣L1⋅(L2×L3)∣V^2=\tfrac29\,|L_1\cdot(L_2\times L_3)|V2=92​∣L1​⋅(L2​×L3​)∣.
  3. Eq. (13), area gap: min⁡{j(j+1):j∈{12,1,32,… }}=32\min\{\sqrt{j(j+1)}: j\in\{\tfrac12,1,\tfrac32,\dots\}\}=\tfrac{\sqrt3}{2}min{j(j+1)​:j∈{21​,1,23​,…}}=23​​.
  4. Eq. (18), area gap in physical units: the least positive element of the area spectrum is 8πℏGγ 328\pi\hbar G\gamma\,\tfrac{\sqrt3}{2}8πℏGγ23​​.
  5. Eq. (26), discreteness: for γ,ℏ,G>0\gamma,\hbar,G>0γ,ℏ,G>0 the area spectrum has only finitely many elements below any bound.

Significance

The goal theorem is the precise mathematical content behind the statement that "the algebraic structure of the momentum operators determines the existence of a metric at each node and therefore equips each quantum of space with a geometry" (§2.1). Its necessity half (milestone 1) is the classical counterpart of the gauge-invariance constraint (6); its sufficiency and uniqueness halves are what allow LQG states to be read as quantum polyhedra and underlie the "twisted geometry" and polyhedral-volume constructions. The spectral milestones encode the review's discrete-geometry claims in the elementary, operator-free form in which they are stated in eqs. (13), (18) and (26).

All statements here are known theorems; none is open. A search of Mathlib (at the revision used to check these statements) found no theory of convex polytopes, facet areas, the Brunn–Minkowski inequality, or Minkowski's theorem for polyhedra.

Difficulty

The necessity direction (closure) is a divergence-theorem or projection argument, but requires relating the facets of an H-described polytope to its boundary. The existence direction of Minkowski's theorem is a variational argument (maximize volume subject to a linear constraint on support numbers) whose difficulty lies in showing that the optimum has all prescribed faces with positive area; uniqueness is classically deduced from the equality case of the Brunn–Minkowski (or Minkowski mixed-volume) inequality. The spectral milestones are elementary.

Formalization scope

Lean represents E3\mathbb E^3E3 as EuclideanSpace ℝ (Fin 3). A polyhedron is given by a finite family of normals and support numbers (halfspacePolyhedron), faces by polyhedronFace, and face area by planarArea, the Haar (Lebesgue) measure of the orthogonal projection onto the plane orthogonal to the normal, valued in [0,∞][0,\infty][0,∞]; areas are compared as extended reals so that an infinite area cannot be confused with 000. Unit normals are obtained as n/∥n∥n/\|n\|n/∥n∥ (unitDir). The cross product is Mathlib's crossProduct transported to EuclideanSpace. Spins are encoded as k/2k/2k/2 with k∈Nk\in\mathbb Nk∈N, and the area spectrum sums over finite multisets of spins.

The goal excludes trivializing readings: the spanning and distinct-direction hypotheses are the classical non-degeneracy conditions of Minkowski's theorem (without them the existence claim is false), and the conclusion demands a bounded polyhedron with nonempty interior, so empty or flat polyhedra do not count. The review's eq. (10) omits numerical constants; the milestone uses the exact Euclidean constant 2/92/92/9. The operator-theoretic statements of the review (spectra of area and volume operators on L2(SU(2)L/SU(2)N)L^2(SU(2)^L/SU(2)^N)L2(SU(2)L/SU(2)N)) are out of scope; only their stated numerical content is formalized.

Reusable infrastructure: H-polytopes in R3\mathbb R^3R3, facet areas, the closure identity, and the Brunn–Minkowski equality case. Contributions of general convex-geometry lemmas are welcome.

Selected references

  • C. Rovelli, Loop quantum gravity: the first 25 years, Class. Quantum Grav. 28 (2011) 153002. https://doi.org/10.1088/0264-9381/28/15/153002 (arXiv:1012.4707)
  • H. Minkowski, Allgemeine Lehrsätze über die konvexen Polyeder, Nachr. Ges. Wiss. Göttingen (1897) 198–219.
  • E. Bianchi, P. Doná, S. Speziale, Polyhedra in loop quantum gravity, Phys. Rev. D 83 (2011) 044035. https://arxiv.org/abs/1009.3402
  • A. D. Alexandrov, Convex Polyhedra, Springer (2005).
7 thms3 active usersReviewed
CombinatoricsLinear Optimization·Captain: mikedeng1

The d-Step Conjecture for Polyhedra of Dimension d < 6 1: The Bounded 5-Step Conjecture Holds, Every 5-Polytope with 10 Facets Has Diameter at Most 5, i.e. Δb(5, 10) = 5Research Paper

Motivation

The simplex method and every other edge-following linear programming algorithm moves from vertex to vertex of the feasible region along its edges. The number of steps such an algorithm needs in the worst case is at least the combinatorial diameter of the feasible polyhedron: the largest distance, counted in edges, between two of its vertices. Klee and Walkup introduced the two functions that measure this worst case, Δ(d,n)\Delta(d,n)Δ(d,n) for polyhedra and Δb(d,n)\Delta_b(d,n)Δb​(d,n) for bounded polyhedra of dimension ddd with nnn facets, and wrote that Δ(d,n)\Delta(d,n)Δ(d,n) "represents, in a sense, the number of iterations required to solve the 'worst' linear program of nnn inequalities in ddd variables using the 'best' edge-following algorithm" (Klee and Walkup 1967, p. 54).

The Hirsch conjecture, reported by Dantzig in 1963, asserted Δ(d,n)≤n−d\Delta(d,n)\le n-dΔ(d,n)≤n−d. Its diagonal cases Δ(d,2d)≤d\Delta(d,2d)\le dΔ(d,2d)≤d and Δb(d,2d)≤d\Delta_b(d,2d)\le dΔb​(d,2d)≤d are the general and the bounded ddd-step conjectures.

Timeline.

  • 1957–1963: W. M. Hirsch's conjecture, reported in Dantzig's Linear Programming and Extensions (1963, pp. 160, 168), with the claim that it holds for m≤4m\le 4m≤4.
  • Before 1967: Klee determines Δ(3,n)=n−3\Delta(3,n)=n-3Δ(3,n)=n−3 and Δb(3,n)=⌊2n/3⌋−1\Delta_b(3,n)=\lfloor 2n/3\rfloor-1Δb​(3,n)=⌊2n/3⌋−1 and settles Δb(d,n)\Delta_b(d,n)Δb​(d,n) for n≤d+4n\le d+4n≤d+4 (Canad. J. Math. 16 and the papers cited as [9], [12], [13] by Klee and Walkup, p. 55).
  • 1967: Klee and Walkup (Acta Math. 117) prove Δb(5,10)=5\Delta_b(5,10)=5Δb​(5,10)=5, show that Δb(d,2d)≤d\Delta_b(d,2d)\le dΔb​(d,2d)≤d for all ddd is equivalent to the bounded Hirsch conjecture, and disprove the unbounded conjecture with Δ(4,8)=5\Delta(4,8)=5Δ(4,8)=5.
  • 2012: Santos disproves the bounded Hirsch conjecture with a 43-dimensional polytope with 86 facets (Ann. Math. 176). The bounded ddd-step conjecture therefore fails for some ddd, and the smallest failing ddd is unknown; it is at least 6.

Setting

A polyhedron in Rd\mathbb R^dRd is a nonempty set P={x:⟨ai,x⟩≤bi, i=1,…,n}P=\{x:\langle a_i,x\rangle\le b_i,\ i=1,\dots,n\}P={x:⟨ai​,x⟩≤bi​, i=1,…,n}; a polytope is a bounded polyhedron. PPP is a ddd-polyhedron with nnn facets when it has nonempty interior and none of the nnn inequalities is redundant; then the rows are in bijection with the facets, and row iii is tight at vvv exactly when vvv lies on the iii-th facet. The vertices of PPP are its extreme points. Two vertices u≠vu\ne vu=v span an edge when the segment [u,v][u,v][u,v] is an extreme subset of PPP. A path of length lll is a sequence of lll edges [x0,x1],…,[xl−1,xl][x_0,x_1],\dots,[x_{l-1},x_l][x0​,x1​],…,[xl−1​,xl​], and δP(x,y)\delta_P(x,y)δP​(x,y) is the least length of a path from xxx to yyy. Then

Δb(d,n)=max⁡{δP(x,y):P a d-polytope with n facets, x,y vertices of P}.\Delta_b(d,n)=\max\{\delta_P(x,y): P\text{ a }d\text{-polytope with }n\text{ facets},\ x,y\text{ vertices of }P\}.Δb​(d,n)=max{δP​(x,y):P a d-polytope with n facets, x,y vertices of P}.

PPP is simple when every vertex lies on exactly ddd facets. A face is PPP itself or the intersection of PPP with a supporting hyperplane, never ∅\emptyset∅; its dimension is that of its affine hull. A (k1,…,kr)(k_1,\dots,k_r)(k1​,…,kr​)-path from xxx to yyy is a sequence of faces K1,…,KrK_1,\dots,K_rK1​,…,Kr​ with x∈K1x\in K_1x∈K1​, y∈Kry\in K_ry∈Kr​, dim⁡Ki=ki\dim K_i=k_idimKi​=ki​ and Ki∩Ki+1≠∅K_i\cap K_{i+1}\ne\emptysetKi​∩Ki+1​=∅. A ddd-dimensional Dantzig figure (P,x,y)(P,x,y)(P,x,y) is a ddd-polyhedron PPP with 2d2d2d facets and two vertices x,yx,yx,y such that exactly ddd facets are incident to xxx and the other ddd to yyy. A ddd-polytope QQQ has property A when it has between 2d2d2d and 2d+22d+22d+2 facets and the following holds: if its facets are split into two classes X,Y\mathfrak X,\mathfrak YX,Y of at most d+1d+1d+1 facets each, and both sets XXX, YYY of vertices all of whose facets lie in X\mathfrak XX, respectively Y\mathfrak YY, are nonempty, then some vertex of XXX is joined to some vertex of YYY by a path of length at most ddd.

Formalization targets

Goal: 4.3 THEOREM

Δb(5,10)=5.\Delta_b(5,10)=5.Δb​(5,10)=5.

In two halves: every 555-polytope with 101010 facets satisfies δP(x,y)≤5\delta_P(x,y)\le 5δP​(x,y)≤5 for all vertices x,yx,yx,y, and some 555-polytope with 101010 facets has two vertices with δP(x,y)>4\delta_P(x,y)>4δP​(x,y)>4.

Milestones, in the order the proof uses them

  1. 2.8 THEOREM (bounded part): Δb(d,n)\Delta_b(d,n)Δb​(d,n) is attained by two vertices of a simple ddd-polytope with nnn facets, and for n≥2dn\ge 2dn≥2d by two vertices on no common facet. For n=2dn=2dn=2d this is a bounded simple Dantzig figure.
  2. 3.4 PROPOSITION (a): every simple Dantzig figure admits a (k1,1,k3)(k_1,1,k_3)(k1​,1,k3​)-path for positive k1+1+k3=dk_1+1+k_3=dk1​+1+k3​=d.
  3. 3.1 PROPOSITION: (a) facial paths in all ddd-dimensional simple Dantzig figures, (b) facet-avoiding facial paths in simple ddd-polytopes with 2d+12d+12d+1 facets, and (c) facet-avoiding (1,k1,…,kr)(1,k_1,\dots,k_r)(1,k1​,…,kr​)-paths in bounded (d+1)(d+1)(d+1)-dimensional figures are equivalent.
  4. 3.4 PROPOSITION (b): every bounded simple Dantzig figure admits a (1,k2,1,k4)(1,k_2,1,k_4)(1,k2​,1,k4​)-path for positive 1+k2+1+k4=d1+k_2+1+k_4=d1+k2​+1+k4​=d.
  5. Remark after 3.4: adjacent terms of a facial path in a simple figure merge into their sum.
  6. Remark after 3.4, p. 68: bounded simple Dantzig figures of dimension d≥2d\ge 2d≥2 admit (1,d−2,1)(1,d-2,1)(1,d−2,1)-paths.
  7. 4.1 PROPOSITION: if the middle face QQQ of a (1,d−2,1)(1,d-2,1)(1,d−2,1)-path ([x,x′],Q,[y′,y])([x,x'],Q,[y',y])([x,x′],Q,[y′,y]) has property A, then δP(x,y)≤d\delta_P(x,y)\le dδP​(x,y)≤d.
  8. 4.2 LEMMA: every simple 333-polytope with 666, 777 or 888 facets has property A.

Significance

The result. By the paper's reductions, Δb(d,2d)≤d\Delta_b(d,2d)\le dΔb​(d,2d)≤d for a given ddd is equivalent to the bounded Hirsch bound Δb(d′,d′+d)≤d\Delta_b(d',d'+d)\le dΔb​(d′,d′+d)≤d for every d′d'd′. So 4.3 gives Δb(d,n)≤n−d\Delta_b(d,n)\le n-dΔb​(d,n)≤n−d for all n≤d+5n\le d+5n≤d+5, which is the column n−d=5n-d=5n−d=5 of the paper's Fig. 2. Since Santos's counterexample, the smallest ddd with Δb(d,2d)>d\Delta_b(d,2d)>dΔb​(d,2d)>d is a well-defined open number. This theorem was the first to show it is at least 666; later computer-assisted work treated further small cases, including the bounded 666-step case (Bremner and Schewe 2011). Santos's construction starts from exactly the reduction to Dantzig figures (the "strong ddd-step theorem" for spindles) that this paper uses.

Formalizing it. The result has been proved since 1967; no machine-checked proof of Δb(5,10)=5\Delta_b(5,10)=5Δb​(5,10)=5 exists. The platform has the Hirsch bound in dimension ≤3\le 3≤3 (Hirsch.dimension_three_bound, Hirsch.klee_three_dimensional_bound), the Klee–Walkup wedge reduction at the level of inequality descriptions (Hirsch.balanced_hpoly_transfer), products (Hirsch.diamLE_prod) and Santos's counterexample. The mission adds the facial-path calculus of §3 and the property-A argument of §4. Both are reusable in any later formal treatment of small ddd-step cases.

Difficulty

The obvious approach, a direct case analysis of 555-polytopes with 101010 facets, is out of reach: there are far too many combinatorial types, and nothing bounds the search. The paper first reduces to simple polytopes, and the reduction itself (2.8) rests on a perturbation lemma (2.6) whose proof is deferred to a separate paper by Walkup: simplicity must be gained without shrinking any distance. The facial-path results then pass between bounded and unbounded figures through truncation and the wedge (3.1), and a statement about all unbounded ddd-dimensional figures is what yields one about bounded (d+1)(d+1)(d+1)-dimensional ones. Finally the (d−2)(d-2)(d−2)-face QQQ must be shown to have property A. That property fails for some simple polytopes of classes (4,10)(4,10)(4,10) and (5,11)(5,11)(5,11) (p. 69), so 4.2 depends on the dimension being 333 and on the facet count being at most 888. Its proof uses Euler's relation and Balinski's theorem that the graph of a 333-polytope is 333-connected (posed on the platform as Balinski61.Connectivity.disjoint_paths_between_vertices).

Formalization scope

Polyhedra are Hirsch.Hpoly a b with a : Fin n → EuclideanSpace ℝ (Fin d); edges and walks are Hirsch.Adj and Hirsch.Reach from the published Hirsch_model and Hirsch_walk. The conventions are:

  • "ddd-polyhedron with nnn facets" is a facet presentation: nonempty interior and no redundant row. The rows are exactly the facets; a statement over arbitrary inequality systems would count redundant rows as facets and is not the paper's.
  • δP(x,y)≤L\delta_P(x,y)\le LδP​(x,y)≤L is Hirsch.Reach P L x y, a walk of LLL steps that may stand still. Δb(5,10)=5\Delta_b(5,10)=5Δb​(5,10)=5 is an upper bound for every polytope plus a witness, never a supremum of reals. A version proving only the upper bound, or only the lower bound, is not the goal.
  • Simple means pointed with exactly ddd tight rows at every vertex. On a facet presentation this is equivalent to the paper's "exactly ddd edges at each vertex".
  • Faces are nonempty sets {x∈P:⟨ai,x⟩=bi, i∈I}\{x\in P:\langle a_i,x\rangle=b_i,\ i\in I\}{x∈P:⟨ai​,x⟩=bi​, i∈I}, and dimension is the dimension of the affine span. A Dantzig figure has Fin (2 * d) rows.
  • Property A is defined for a bounded face QQQ of dimension mmm: its facets are the (m−1)(m-1)(m−1)-faces of PPP inside QQQ, and the path runs in QQQ. This keeps 4.1, where QQQ is a (d−2)(d-2)(d−2)-face of a polytope in Rd\mathbb R^dRd, and 4.2, where Q=PQ=PQ=P, under one notion. Facet counts use Set.encard.
  • Explicit readings: 2.8's "Δb\Delta_bΔb​ can be realized as the distance … of a simple polytope" is stated polytope by polytope (δP(u,v)≤δP′(x,y)\delta_P(u,v)\le\delta_{P'}(x,y)δP​(u,v)≤δP′​(x,y)). Only the Δb\Delta_bΔb​ sentences of 2.8 are stated. 3.4's joint hypothesis on k1,…,k4k_1,\dots,k_4k1​,…,k4​ is split by part. The merging remark carries the needed hypothesis ki+ki+1≤dk_i+k_{i+1}\le dki​+ki+1​≤d. The printed "3.5 implies" on p. 67 is read as 3.4.

The paper cites two inputs without proving them: Lemma 2.6 (Walkup's simple perturbation) and Balinski's theorem. A complete development needs face lattices of H-polyhedra, wedges and truncations, and the planar-graph arguments of 4.2. The face and facial-path layer is reusable for the companion mission on Δ(4,8)=5\Delta(4,8)=5Δ(4,8)=5 and for other diameter problems. Contributions are welcome at every milestone, including proofs of 2.6.

Selected references

  • V. Klee and D. W. Walkup, The d-step conjecture for polyhedra of dimension d < 6, Acta Mathematica 117 (1967), 53–78. https://doi.org/10.1007/BF02395040
  • G. B. Dantzig, Linear Programming and Extensions, Princeton University Press, 1963.
  • V. Klee, Diameters of polyhedral graphs, Canadian Journal of Mathematics 16 (1964), 602–614. https://doi.org/10.4153/CJM-1964-058-5
  • M. L. Balinski, On the graph structure of convex polyhedra in n-space, Pacific Journal of Mathematics 11 (1961), 431–434. https://doi.org/10.2140/pjm.1961.11.431
  • F. Santos, A counterexample to the Hirsch conjecture, Annals of Mathematics 176 (2012), 383–412. https://doi.org/10.4007/annals.2012.176.1.7
  • D. Bremner and L. Schewe, Edge-graph diameter bounds for convex polytopes with few facets, Experimental Mathematics 20 (2011), 229–237. https://doi.org/10.1080/10586458.2011.564966
12 thms3 active usersReviewed
CombinatoricsLinear Optimization·Captain: mikedeng1

A Quasi-Polynomial Bound for the Diameter of Graphs of Polyhedra: Every Polyhedron Given by n Inequalities in ℝ^d, Bounded or Not, Has Graph Diameter at Most n^(log₂ d + 2)Research Paper

Motivation

The simplex method solves a linear program by walking along the edges of the feasible polyhedron, from vertex to vertex. However cleverly a pivot rule chooses its steps, it cannot reach the optimum in fewer steps than the graph distance from the starting vertex. The largest possible graph diameter of a polyhedron therefore bounds from below the worst-case number of pivots of every pivot rule, and whether that diameter is polynomial in the size of the problem is one of the basic open questions of linear optimization (Klee–Kleinschmidt 1987).

Timeline.

  • 1957: Hirsch conjectures that a ddd-dimensional polyhedron with nnn facets has graph diameter at most n−dn-dn−d (reported in Dantzig, Linear Programming and Extensions, 1963).
  • 1967: Klee and Walkup show the conjecture false for unbounded polyhedra, with the lower bound n−d+[d/5]n-d+[d/5]n−d+[d/5] for n≥2dn\ge2dn≥2d (Acta Math. 1967).
  • 1967–1970: Barnette proves Δ(d,n)≤n3d−3\Delta(d,n)\le n3^{d-3}Δ(d,n)≤n3d−3, and Larman improves it to n2d−3n2^{d-3}n2d−3 (Proc. LMS 1970); both are linear in nnn but exponential in ddd.
  • 1990–1992: Kalai proves a subexponential bound, and then Kalai and Kleitman prove the quasi-polynomial bound nlog⁡d+2n^{\log d+2}nlogd+2 that is the subject of this mission (Bull. AMS 1992, arXiv:math/9204233).
  • 2012: Santos disproves the Hirsch conjecture for bounded polytopes (Ann. Math. 2012). Whether a polynomial bound holds remains open; the best general upper bounds known are refinements of the exponent of the Kalai–Kleitman bound.

Setting

A polyhedron in Rd\mathbb R^dRd is the solution set of finitely many linear inequalities,

P={x∈Rd: ⟨ai,x⟩≤bi, i=1,…,n},ai∈Rd, bi∈R.P=\{x\in\mathbb R^d:\ \langle a_i,x\rangle\le b_i,\ i=1,\dots,n\},\qquad a_i\in\mathbb R^d,\ b_i\in\mathbb R .P={x∈Rd: ⟨ai​,x⟩≤bi​, i=1,…,n},ai​∈Rd, bi​∈R.

It may be bounded or unbounded, and empty. Its vertices are its extreme points. Two distinct vertices u,vu,vu,v are adjacent when the segment [u,v][u,v][u,v] is an edge of PPP, a one-dimensional face; unbounded one-dimensional faces (rays) are not edges. The graph G(P)G(P)G(P) consists of these vertices and edges, and its diameter δ(P)\delta(P)δ(P) is the largest graph distance between two vertices. Kalai and Kleitman write Δ(d,n)\Delta(d,n)Δ(d,n) for the largest diameter of the graph of a ddd-dimensional polyhedron with nnn facets.

For a set FFF of row indices, the relaxation QFQ_FQF​ is the polyhedron obtained from PPP by ignoring the inequalities with index outside FFF. A walk of BBB steps is a sequence w0,…,wBw_0,\dots,w_Bw0​,…,wB​ of points in which each step stays put or crosses an edge; "diameter at most BBB" means that any two vertices are joined by such a walk, which includes connectivity of the graph.

Formalization targets

Goal: Theorem 1

Δ(d,n)≤nlog⁡d+2.\Delta(d,n)\le n^{\log d+2}.Δ(d,n)≤nlogd+2.

Stated for polyhedra given by inequalities: for all d,nd,nd,n and all a∈(Rd)na\in(\mathbb R^d)^na∈(Rd)n, b∈Rnb\in\mathbb R^nb∈Rn, the graph of PPP is connected and has diameter at most ⌊nlog⁡2d+2⌋\lfloor n^{\log_2 d+2}\rfloor⌊nlog2​d+2⌋. No nonemptiness, boundedness or full-dimensionality is assumed.

Milestones, in the order the proof uses them

  1. If a vertex ω\omegaω of PPP satisfies with equality only rows in FFF, then ω\omegaω is a vertex of QFQ_FQF​ (published).
  2. If an edge of QFQ_FQF​ starts at a point of PPP and is not an edge of PPP, then it meets the hyperplane of an ignored row at a vertex of PPP adjacent to its start, or at its start (published).
  3. The claim kv≤Δ(d,[n/2])k_v\le\Delta(d,[n/2])kv​≤Δ(d,[n/2]): if every vertex within distance kkk of vvv lies only on rows of FFF, and G(QF)G(Q_F)G(QF​) has diameter at most BBB, then every vertex within distance kkk of vvv is within distance BBB of vvv.
  4. The facet step: any two vertices of a face {x∈P:⟨ai,x⟩=bi}\{x\in P:\langle a_i,x\rangle=b_i\}{x∈P:⟨ai​,x⟩=bi​} with ai≠0a_i\ne0ai​=0 are joined, inside that face, by a walk of at most Δ(d−1,n−1)\Delta(d-1,n-1)Δ(d−1,n−1) edges.
  5. The recursion Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,[n/2])+2\Delta(d,n)\le\Delta(d-1,n-1)+2\Delta(d,[n/2])+2Δ(d,n)≤Δ(d−1,n−1)+2Δ(d,[n/2])+2, for polyhedra with nonzero rows.
  6. The arithmetic: N(d,n)=⌊nlog⁡2d+2⌋N(d,n)=\lfloor n^{\log_2 d+2}\rfloorN(d,n)=⌊nlog2​d+2⌋ satisfies N(d−1,n−1)+2N(d,⌊n/2⌋)+2≤N(d,n)N(d-1,n-1)+2N(d,\lfloor n/2\rfloor)+2\le N(d,n)N(d−1,n−1)+2N(d,⌊n/2⌋)+2≤N(d,n) for d≥2d\ge2d≥2, n≥d+1n\ge d+1n≥d+1.

Significance

The bound is polynomial in nnn for each fixed ddd, and nO(log⁡d)n^{O(\log d)}nO(logd) overall. Before it, every bound was exponential in ddd or subexponential in nnn. It remains, up to improvements of the exponent, the best known upper bound on the diameter of polyhedra, and thus the best known limit on how badly an optimal pivot rule must behave. Its proof is also the template for the later abstractions of the problem (connected layer families, subset partition graphs) in which the polynomial Hirsch question is studied.

The bounded case is already formalized on the platform: Hirsch.kalai_kleitman_bound states the theorem for nonempty bounded polytopes and is Proved, together with the bounded forms of the recursion (Hirsch.kalai_kleitman_recursion), of the facet step (Hirsch.facet_reduction, Hirsch.facet_walk) and of the zero-row reduction (Hirsch.diamLE_of_nonzero_rows), and the one-step relaxation lemma Hirsch.relaxation_no_new_neighbours. Kalai and Kleitman state Theorem 1 for all polyhedra: their introduction contrasts unbounded polyhedra (where Hirsch fails) with bounded ones. This mission formalizes the theorem as the paper states it, including unbounded polyhedra, where the bounded proofs do not apply as written. The combinatorial abstraction of the same bound (EHRR10.KalaiKleitman.diameter_le_kalai_kleitman, for the Eisenbrand–Hähnle–Razborov–Rothvoß model) is a different object and is not reused.

Difficulty

The bounded argument relies on compactness: faces are polytopes, every relaxation used in the induction is again a polytope, and the facet of a polytope is a polytope of one dimension less. The relaxation QFQ_FQF​ contains PPP and is typically unbounded even when PPP is bounded, and for an unbounded polyhedron a face on a facet hyperplane can be unbounded or have no vertices at all. Each step must therefore hold for general polyhedra: that a shortest walk in G(QF)G(Q_F)G(QF​) between two vertices of the ball uses only edges of PPP, and that the graph of a face, identified with a polyhedron in Rd−1\mathbb R^{d-1}Rd−1, sits inside G(P)G(P)G(P). The "maximal positive number kvk_vkv​" of the paper need not exist, and the proof must handle the two cases (no radius qualifies; every radius qualifies) separately. Zero rows are allowed in the goal and must be removed before the recursion applies, since a zero row with bi=0b_i=0bi​=0 is tight everywhere.

Formalization scope

Lean represents Rd\mathbb R^dRd as EuclideanSpace ℝ (Fin d) and the polyhedron as Hirsch.Hpoly a b with a : Fin n → ℝ^d, b : Fin n → ℝ, from the published definition Hirsch_model, which also supplies adjacency Hirsch.Adj (the segment is an extreme subset) and the diameter predicate Hirsch.DiamLE. The mission commits to these explicit readings:

  • "Δ(d,n)≤N\Delta(d,n)\le NΔ(d,n)≤N" is ∀ a b, DiamLE (Hpoly a b) N over nnn inequalities in Rd\mathbb R^dRd. This is equivalent to the facet formulation, since a ddd-dimensional polyhedron is cut out by its facet inequalities and a polyhedron given by nnn inequalities in Rd\mathbb R^dRd has dimension at most ddd and at most nnn facets.
  • "log" is log⁡2\log_2log2​ (Real.logb 2), the power is the real power, and the integer bound is its floor ⌊·⌋₊, which is equivalent because distances are natural numbers.
  • "[n/2][n/2][n/2]" is ⌊n/2⌋\lfloor n/2\rfloor⌊n/2⌋ (ℕ division); "at most n/2n/2n/2 facets" is a row set of cardinality at most ⌊n/2⌋\lfloor n/2\rfloor⌊n/2⌋.
  • An ignored row is replaced by the vacuous inequality ⟨0,x⟩≤1\langle0,x\rangle\le1⟨0,x⟩≤1, which keeps the index type Fin n.
  • Walks have exactly BBB steps with stationary steps allowed, so "BBB steps" means "at most BBB edges".
  • The paper's kvk_vkv​ is not defined; milestone 3 quantifies over every radius with the ball property.
  • Milestones 4 and 5 assume nonzero rows; the goal does not.

A statement with a boundedness or nonemptiness hypothesis, with the graph diameter measured by a distance function whose value on unreachable pairs is 000, or with the bound read with Nat.log or the natural logarithm, is not the target. Contributions welcome: proofs of the general-polyhedron facet step and of milestone 3, a general zero-row reduction, the arithmetic of milestone 6, and the base cases d≤1d\le1d≤1 and n≤dn\le dn≤d.

Selected references

  • 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://doi.org/10.1090/S0273-0979-1992-00285-9 (arXiv:math/9204233)
  • 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
  • V. Klee and P. Kleinschmidt, The d-step conjecture and its relatives, Math. Oper. Res. 12 (1987), 718–755. https://doi.org/10.1287/moor.12.4.718
  • D. G. Larman, Paths on polytopes, Proc. London Math. Soc. (3) 20 (1970), 161–178. https://doi.org/10.1112/plms/s3-20.1.161
  • F. Santos, A counterexample to the Hirsch conjecture, Ann. of Math. 176 (2012), 383–412. https://doi.org/10.4007/annals.2012.176.1.7
9 thms3 active usersReviewed
Linear OptimizationOperations ResearchProbability·Captain: mikedeng1

A Friendly Smoothed Analysis of the Simplex Method: Gaussian-Perturbed Shadows Have at Most 3 + 128πe²d²√(log n)σ⁻²(1 + 4σ√(d log n))(1 + 16σ√(log n)) Expected EdgesResearch Paper

Why the simplex method needs smoothed analysis

The simplex method solves linear programs quickly in practice, yet for most pivot rules there are inputs on which it takes exponentially many steps (Klee and Minty, 1972). Smoothed analysis, introduced by Spielman and Teng (JACM 2004), explains the gap: it measures the expected running time when every constraint row of a worst-case input is perturbed by a small random amount of size σ\sigmaσ. For the shadow vertex pivot rule, the running time is controlled by a purely geometric quantity, the number of edges of a random polygon, and the smoothed complexity of the method is essentially an upper bound on that number.

Timeline (bounds on the expected shadow size of smoothed unit LPs, as surveyed in §1.1 of the paper). Borgwardt (1987) proved Θ(d1.5log⁡n)\Theta(d^{1.5}\sqrt{\log n})Θ(d1.5logn​) for rows drawn from a centered Gaussian, asymptotically as n→∞n\to\inftyn→∞; this is the case aˉi=0\bar a_i=0aˉi​=0 of the model below. Spielman and Teng (2004) gave the first smoothed bound, O(d3nσ−6+d6nlog⁡3n)O(d^3n\sigma^{-6}+d^6n\log^3 n)O(d3nσ−6+d6nlog3n). Deshpande and Spielman (FOCS 2005) improved the dependence on σ\sigmaσ to O(dn2log⁡n σ−2+d2n2log⁡2n)O(dn^2\log n\,\sigma^{-2}+d^2n^2\log^2 n)O(dn2lognσ−2+d2n2log2n). Vershynin (SIAM J. Comput. 2009) reduced the dependence on nnn to polylogarithmic, O(d3σ−4+d5log⁡2n)O(d^3\sigma^{-4}+d^5\log^2 n)O(d3σ−4+d5log2n). Dadush and Huiberts (arXiv:1711.05667, STOC 2018, SIAM J. Comput. 2020) proved O(d2log⁡n σ−2+d2.5log⁡n σ−1+d2.5log⁡1.5n)O(d^2\sqrt{\log n}\,\sigma^{-2}+d^{2.5}\log n\,\sigma^{-1}+d^{2.5}\log^{1.5}n)O(d2logn​σ−2+d2.5lognσ−1+d2.5log1.5n) for Gaussian perturbations by a simpler argument, which also gives bounds for Laplace perturbations.

Setting

Fix dimensions n≥d≥3n\ge d\ge3n≥d≥3 and a shadow plane, a two-dimensional linear subspace W⊆RdW\subseteq\mathbb R^dW⊆Rd, chosen independently of the randomness. The constraint vectors a1,…,an∈Rda_1,\dots,a_n\in\mathbb R^da1​,…,an​∈Rd are independent, and ai∼Nd(aˉi,σ)a_i\sim N_d(\bar a_i,\sigma)ai​∼Nd​(aˉi​,σ) is Gaussian with center aˉi\bar a_iaˉi​, ∥aˉi∥≤1\|\bar a_i\|\le1∥aˉi​∥≤1, and standard deviation σ>0\sigma>0σ>0 in every coordinate, i.e. with density (2πσ2)−d/2e−∥x−aˉi∥2/(2σ2)(2\pi\sigma^2)^{-d/2}e^{-\|x-\bar a_i\|^2/(2\sigma^2)}(2πσ2)−d/2e−∥x−aˉi​∥2/(2σ2).

Let Q=conv⁡(a1,…,an)Q=\operatorname{conv}(a_1,\dots,a_n)Q=conv(a1​,…,an​). The shadow polygon is Q∩WQ\cap WQ∩W. An edge of a convex set KKK is a segment [u,v][u,v][u,v], u≠vu\ne vu=v, which is an extreme subset of KKK; ∣edges⁡(K)∣|\operatorname{edges}(K)|∣edges(K)∣ is the number of edges. The number of pivots of the shadow vertex method on the unit LP {x:Ax≤1}\{x:Ax\le\mathbf1\}{x:Ax≤1} between two objectives spanning WWW is at most ∣edges⁡(Q∩W)∣|\operatorname{edges}(Q\cap W)|∣edges(Q∩W)∣ (the paper's Lemma 11 and Theorem 12).

The proof passes through general row distributions with density μ\muμ and mean yyy, described by four parameters: the log-Lipschitz constant LLL (μ(x)≤eL∥x−x′∥μ(x′)\mu(x)\le e^{L\|x-x'\|}\mu(x')μ(x)≤eL∥x−x′∥μ(x′)); the line variance τ2\tau^2τ2 (the least variance of μ\muμ restricted to a line); the nnn-th deviation rnr_nrn​ (the least rrr with ∫r∞Pr⁡[∣(X−y)Tθ∣≥t] dt≤r/n\int_r^\infty\Pr[|(X-y)^\mathsf T\theta|\ge t]\,dt\le r/n∫r∞​Pr[∣(X−y)Tθ∣≥t]dt≤r/n for all unit θ\thetaθ); and the cutoff radius Rn,dR_{n,d}Rn,d​ (the least RRR with Pr⁡[∥X−y∥≥R]≤1/(d(nd))\Pr[\|X-y\|\ge R]\le 1/(d\binom nd)Pr[∥X−y∥≥R]≤1/(d(dn​))). The Gaussian is not log-Lipschitz; the Laplace–Gaussian distribution LGd(aˉ,σ,r)LG_d(\bar a,\sigma,r)LGd​(aˉ,σ,r), whose density equals the Gaussian one inside the ball of radius rσr\sigmarσ around aˉ\bar aaˉ and decays like e−∥x−aˉ∥r/σe^{-\|x-\bar a\|r/\sigma}e−∥x−aˉ∥r/σ outside, replaces it.

Formalization targets

Goal: Theorem 13 with explicit constants

E[∣edges⁡(conv⁡(a1,…,an)∩W)∣]≤3+128πe2d2log⁡nσ2(1+4σdlog⁡n)(1+16σlog⁡n).\mathbb E\big[|\operatorname{edges}(\operatorname{conv}(a_1,\dots,a_n)\cap W)|\big]\le 3+\frac{128\pi e^2 d^2\sqrt{\log n}}{\sigma^2}\big(1+4\sigma\sqrt{d\log n}\big)\big(1+16\sigma\sqrt{\log n}\big).E[∣edges(conv(a1​,…,an​)∩W)∣]≤3+σ2128πe2d2logn​​(1+4σdlogn​)(1+16σlogn​).

The paper states O(d2log⁡n σ−2+d2.5log⁡n σ−1+d2.5log⁡1.5n)O(d^2\sqrt{\log n}\,\sigma^{-2}+d^{2.5}\log n\,\sigma^{-1}+d^{2.5}\log^{1.5}n)O(d2logn​σ−2+d2.5lognσ−1+d2.5log1.5n); the constant above is the one its proof produces (Theorem 22's proof, Lemma 45, Lemma 46).

Intermediate targets (milestones)

  • the parametrized shadow bound (Theorem 22), E∣edges⁡∣≤2+8πe2d1.5(L/τ)(1+Rn,d)(1+4rn)\mathbb E|\operatorname{edges}|\le 2+8\pi e^2d^{1.5}(L/\tau)(1+R_{n,d})(1+4r_n)E∣edges∣≤2+8πe2d1.5(L/τ)(1+Rn,d​)(1+4rn​);
  • its ingredients: the edge-counting reduction (Lemma 25), the perimeter bound 2π(1+4rn)2\pi(1+4r_n)2π(1+4rn​) (Lemmas 19, 26), the conditioning on bounded diameter (Lemma 29), the chord combination calculus (Lemmas 35, 36), and the two expected-length bounds (Lemmas 39, 40, 41), together with the trade-off LR(1/2)≥d/3LR(1/2)\ge d/3LR(1/2)≥d/3 (Lemma 21);
  • the Gaussian specialisation: Laplace–Gaussian tails (Lemma 44), its parameters (Lemma 45), and the comparison with the Gaussian (Lemma 46).

Significance

The bound gives the best known smoothed complexity of a simplex method, with only log⁡n\sqrt{\log n}logn​ dependence on the number of constraints. Section 4 of the paper turns it into a complete two-phase shadow vertex algorithm with O(d2log⁡n σ−2+d3log⁡1.5n)O(d^2\sqrt{\log n}\,\sigma^{-2}+d^3\log^{1.5}n)O(d2logn​σ−2+d3log1.5n) expected pivots. The parametrized Theorem 22 applies to any perturbation family with bounded parameters, so it also yields shadow bounds for Laplace and other noise models.

The result is proved; it has not been machine-checked. A formalization adds a verified chain from four distributional parameters to a polygon edge count. The intermediate statements have independent value: an expected-perimeter bound for random polygons, the deterministic chord-length identity for simplices, and tail estimates for the Laplace–Gaussian distribution.

Difficulty

The expected number of edges is a sum, over all (nd)\binom nd(dn​) index sets, of the probability that their simplex meets WWW in an edge; there are far too many terms for per-set estimates to give a bound polylogarithmic in nnn. The paper's route (inspired by Kelner and Spielman) instead compares the expected perimeter with the expected length of an edge, and the hard part is a lower bound on the expected length of an edge conditioned on it appearing. That conditioning is on a probability-zero configuration (the simplex lying in a fixed hyperplane), and its analysis needs Blaschke's change of variables, the shape decomposition of the simplex, and log-Lipschitz shifts of densities. On the Gaussian side, the density is not log-Lipschitz, which forces the Laplace–Gaussian detour and a comparison of edge counts under two different laws.

Formalization scope

Rd\mathbb R^dRd is EuclideanSpace ℝ (Fin d), [n][n][n] is Fin n, log⁡\loglog is the natural logarithm. The Gaussian is the published SmoothedSimplex.Shadow.gaussian, the edge relation is the published Hirsch.Adj, and the rows' joint law is Measure.pi. Edges are segments counted once; the perimeter is the sum of edge lengths. Expectations are lower Lebesgue integrals in [0,∞][0,\infty][0,∞], and every edge-count bound also asserts almost-everywhere measurability of the count, so the integral is a genuine expectation. Conditional expectations are stated multiplied out (no division by a probability), and the parameters of Definitions 15–18 are upper-bound predicates (lower bound for the line variance), which is how every statement of the paper uses them.

Hypotheses added relative to the page, each implicit there: σ>0\sigma>0σ>0; L,τ>0L,\tau>0L,τ>0 and R,r≥0R,r\ge0R,r≥0 in Theorem 22; positivity of log-Lipschitz densities; the standing n≥d≥3n\ge d\ge3n≥d≥3 of §3 in Lemma 44. Lemmas 39 and 41 are stated on the explicit conditional densities that the paper derives in Lemma 37, instead of through regular conditional distributions; Lemma 39 takes the conclusion d/3≤LRd/3\le LRd/3≤LR of Lemma 21 as a hypothesis, as its proof uses it. Two printed slips are corrected: Lemma 46's formula is stated for conv⁡(a1,…,an)∩W\operatorname{conv}(a_1,\dots,a_n)\cap Wconv(a1​,…,an​)∩W, as its sentence and proof have it, and Lemma 35 refers to Definition 33 (not 34).

The goal mentions only Gaussian rows, WWW, and the edge count of conv⁡(a)∩W\operatorname{conv}(a)\cap Wconv(a)∩W. It cannot be trivialized by an unsatisfiable parameter hypothesis or a bound on edge lengths, since it has neither; and it counts edges of the polygon, not of the polytope conv⁡(a)\operatorname{conv}(a)conv(a).

A complete development needs: Gaussian and Laplace–Gaussian tail bounds; measurability of the edge count of a random polygon; the perimeter monotonicity of convex sets; Blaschke's change of variables (Theorem 7 of the paper, quoted from Blaschke 1935), which Mathlib does not have; and the regular conditioning of Lemmas 32, 37 and 38. The Blaschke formula, the perimeter facts for convex polygons, and the tail bounds are reusable well beyond this mission. Contributions to any milestone, or to these supporting results as separate theorems, are welcome.

Selected references

  • D. Dadush, S. Huiberts, A Friendly Smoothed Analysis of the Simplex Method, SIAM J. Comput. 49(5), 2020 (STOC 2018); preprint arXiv:1711.05667v4, 2019. https://arxiv.org/abs/1711.05667
  • D. A. Spielman, S.-H. Teng, Smoothed analysis of algorithms: Why the simplex algorithm usually takes polynomial time, J. ACM 51(3), 2004. https://doi.org/10.1145/990308.990310
  • R. Vershynin, Beyond Hirsch Conjecture: Walks on Random Polytopes and Smoothed Complexity of the Simplex Method, SIAM J. Comput. 39(2), 2009. https://doi.org/10.1137/070683386
  • A. Deshpande, D. A. Spielman, Improved smoothed analysis of the shadow vertex simplex method, FOCS 2005 (cited as [DS05] in the paper).
  • V. Klee, G. J. Minty, How good is the simplex algorithm?, in Inequalities III, Academic Press, 1972 (cited as [KM70] in the paper).
20 thms3 active usersReviewed
Convex OptimizationOperations ResearchOptimization·Captain: Shuze Chen

Discrete Convex Analysis XXVIII: The Conjugacy TheoremTextbook

Motivation

Chapter 8 is where discrete convex analysis explains why it needed two separate notions — M-convexity (exchangeability) and L-convexity (submodularity) — rather than one. The answer is conjugacy: under the classical Legendre-Fenchel transform, the two classes turn out to be exactly dual to each other, the discrete analogue of the fact that convex analysis's transform is self-dual within a single class of convex functions. Mission 10-conjugacy-i proved the integer-lattice version of this fact (Theorem 8.12) but explicitly deferred the polyhedral version — Theorem 8.4, the chapter's own headline "Conjugacy theorem" — noting it needed a real-variable M-/L-convex-function layer the series had not yet built. That layer now exists, built across missions 23-24-ch06*-mconvexfunctions and 26-27-ch07*-lconvexfunctions. This mission proves Theorem 8.4 and its companions: the polar-cone correspondence it induces, its nonpolyhedral generalization, the separation and Fenchel-duality theorems for M♮-/L♮-convex functions, and the basic theory of M2-convex functions (sums of M-convex functions), which the Edmonds intersection theorem's own combinatorics is built from.

Setting

Fix a finite ground set VVV. For f:RV→R∪{+∞}f : \mathbb R^V \to \mathbb R \cup \{+\infty\}f:RV→R∪{+∞}, the Legendre-Fenchel transform is f∙(p)=sup⁡x[⟨p,x⟩−f(x)]f^\bullet(p) = \sup_x [\langle p,x\rangle - f(x)]f∙(p)=supx​[⟨p,x⟩−f(x)]. A polyhedral convex function fff is M-convex (f∈M[R→R]f \in M[\mathbb R \to \mathbb R]f∈M[R→R]) if it satisfies (M-EXC[R]); ggg is L-convex (g∈L[R→R]g \in L[\mathbb R \to \mathbb R]g∈L[R→R]) if it satisfies (SBF[R]) and (TRF[R]). A concave function hhh is always represented via h2=−hh_2 = -hh2​=−h, an ordinary convex function, so every "f≥hf \ge hf≥h" hypothesis is restated as "f+h2≥0f + h_2 \ge 0f+h2​≥0" — an equivalent formulation avoiding any need to represent −∞-\infty−∞ in the codomain. A polyhedral cone's polar is C∘={y:⟨y,x⟩≤0 ∀x∈C}C^\circ = \{y : \langle y,x\rangle \le 0\ \forall x \in C\}C∘={y:⟨y,x⟩≤0 ∀x∈C}. A function is M2-convex if it is the sum of two M-convex functions.

Formalization targets

Goal: the conjugacy theorem (Theorem 8.4)

The classes of polyhedral M-convex functions and polyhedral L-convex functions are in one-to-one correspondence under the Legendre-Fenchel transform: f∈M⇒f∙∈Lf \in M \Rightarrow f^\bullet \in Lf∈M⇒f∙∈L, g∈L⇒g∙∈Mg \in L \Rightarrow g^\bullet \in Mg∈L⇒g∙∈M, and the transform is an involution (f∙∙=ff^{\bullet\bullet}=ff∙∙=f, g∙∙=gg^{\bullet\bullet}=gg∙∙=g) on each class, with the identical statement for the M♮^\natural♮/L♮^\natural♮ variants. This is the theorem mission 10-conjugacy-i deferred, citing exactly the missing infrastructure this series has since built.

Supporting structural targets

Twelve further results build the surrounding theory. Proposition 8.2 gives the easy two-variable case of the general submodularity-preservation fact (Theorem 8.1, already a milestone of mission 10-conjugacy-i); Proposition 8.3 is the technical minimizer-difference lemma the goal's harder direction is built from. Theorem 8.5 derives the M-convex/L-convex cone polarity from the goal, and Theorem 8.6 extends the correspondence beyond the polyhedral case to general closed proper convex functions. Proposition 8.14 and Theorems 8.15-8.16 build the separation theory for M♮-/L♮-convex and concave function pairs, with integral witnesses when the functions are integer valued; Theorem 8.21 (parts 1-2) derives the Fenchel-type strong-duality equality these separation theorems make possible. Propositions 8.29-8.30 and Theorem 8.31 (plus Theorem 8.32, found by direct reading immediately after 8.31) build the basic theory of M2-convex functions: their domains and minimizer sets are M2-convex, they are integrally convex, and their global optimality reduces to a finite local check.

Significance

The goal is the theorem that retroactively explains this entire series' two-track structure: missions 20-25 (M-convex sets and functions) and 08/21/26-28 (L-convex sets and functions) are not two independent theories that happen to share techniques — they are conjugate images of each other, so every theorem proved on one side has a dual counterpart automatically available on the other via Theorem 8.4. This is made concrete immediately: Theorem 8.5's cone polarity and the diagram the book draws connecting M0[R]M_0[\mathbb R]M0​[R], 0L[R→R]0L[\mathbb R\to\mathbb R]0L[R→R], and submodular set functions S[R]S[\mathbb R]S[R] (already correspondences this series proved independently, in missions 24-ch06d-mconvexfunctions and 28-ch07d-lconvexfunctions) are shown to be facets of one single conjugacy fact rather than three separate coincidences. The separation and Fenchel duality theorems (8.15, 8.16, 8.21) are the discrete analogues of the two theorems every convex optimization course opens with, and the book is explicit that they are not corollaries of the classical versions plus convex extensibility — they carry genuinely combinatorial content, specializing to Frank's discrete separation theorem and Edmonds's intersection theorem as examples the book itself gives.

None of these results are open — they are Murota's account of the duality at the heart of discrete convex analysis, the reason the theory needed two dual notions rather than one. What this mission contributes is a faithful, machine-checked formal statement of each, completing a theorem mission 10-conjugacy-i explicitly left for a future session once the necessary polyhedral apparatus existed, and including one result (Theorem 8.32) the platform's own automated extractor missed; no comparable formalization exists on the platform (see Formalization scope).

Difficulty

The naive approach to the goal's harder direction (L⇒M) would try to verify the exchange inequality for g∙g^\bulletg∙ directly from the definition of the transform; the book's actual proof instead identifies the exchange inequality with a statement about weighted minimizers of ggg itself via Proposition 8.3 (the minimizer-difference bound), converting a claim about the conjugate function into a claim about ggg's own combinatorial structure — a genuine change of perspective, not a direct calculation. Proposition 8.3's own proof is the hardest single argument in this block: it derives the minimizer-difference bound by a contradiction argument that constructs an explicit pair of "worse" minimizers via a join/meet perturbation and derives a strict inequality from Theorem 7.29's translation inequality — a multi-step combinatorial argument with no direct shortcut.

Formalization scope

Ground-set elements are a Fintype V with DecidableEq; convex functions are WithTop ℝ valued throughout (never EReal, except for the Legendre-Fenchel transform itself, whose defining supremum/infimum can genuinely be infinite). All thirteen numbered results found in this chunk's page range — the twelve in BRIEF.md's own table plus Theorem 8.32 — are placed, with one documented scope reduction: Theorem 8.21 states only its real-attainment parts (1)-(2), not the integer-attainment refinement of parts (3)-(4), which needs a separate argument no other result in this chunk requires — see HARD.md. Concave functions hhh are always represented via h2=−hh_2 = -hh2​=−h and every inequality f≥hf \ge hf≥h restated as f+h2≥0f + h_2 \ge 0f+h2​≥0, avoiding WithTop ℝ negation entirely. This chunk's own BRIEF.md inherited the chapters-4-7 page-offset boilerplate (printed = PDF −-− 19); chapter 8 uses offset 18, confirmed against the PDF's own footers — every citation here uses the corrected offset. This mission's base vocabulary is redeclared from missions 10-conjugacy-i, 20-ch04b-mconvexsets, 21-ch05b-lconvexsets, 23-24-ch06*-mconvexfunctions, and 26-27-ch07*-lconvexfunctions rather than imported, since sibling drafts in this series cannot yet reference one another. Contributions completing any of the thirteen sorrys are welcome; the goal and Proposition 8.3 carry the most independent proof content.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003. DOI: 10.1137/1.9780898718508.
  • K. Murota and A. Shioura, "M-convex function on generalized polymatroid," Mathematics of Operations Research, 24 (1999), pp. 95-105 [152] (the polyhedral M-/L-convex conjugacy theory this mission's real-variable results are drawn from).
75 thms3 active usersReviewed
Number TheoryOperations Research·Captain: mikedeng1

Minkowski's Convex Body Theorem and Integer Programming: Lattice-Free Convex Bodies Meet Few Translates of an Integral SubspaceResearch Paper

Motivation

Integer programming asks whether a system of linear inequalities Ax≤bAx\le bAx≤b has a solution x∈Znx\in\mathbb Z^nx∈Zn. In fixed dimension nnn it is solvable in polynomial time: Lenstra (1983) proved this by showing that a convex body without integer points is "flat" in some integral direction, so that the search splits into few lower-dimensional subproblems. Kannan's 1987 paper in Mathematics of Operations Research sharpened this approach. It computes a Korkine–Zolotarev ("reduced") basis of a lattice, solves the shortest and closest vector problems exactly in nO(n)n^{O(n)}nO(n) operations, and runs integer programming in O(n9n/2s)O(n^{9n/2}s)O(n9n/2s) arithmetic operations. Underneath the algorithm sits a purely geometric statement, Theorem (5.5): a lattice-free convex body meets only boundedly many integer translates of some integral subspace.

Timeline.

  • Korkine and Zolotareff (1873): the reduced bases used here.
  • Minkowski (1896): a symmetric convex body of volume greater than 2n2^n2n contains a nonzero integer point.
  • Khinchine (1948): lattice-free convex bodies have lattice width bounded by a function of nnn alone (the flatness theorem).
  • Lenstra (1983): integer programming in fixed dimension is polynomial, via a flat direction.
  • Kannan (1987, this paper): Theorem (5.5), with subspaces VVV of any dimension between 111 and n−1n-1n−1 and an explicit bound n2(n−dim⁡V)n^{2(n-\dim V)}n2(n−dimV).
  • Kannan and Lovász (1988), Banaszczyk et al. (1999), and later work: polynomial bounds on the flatness constant.

Setting

Rn\mathcal R^nRn is Euclidean space with dot product (a,b)(a,b)(a,b) and length ∣a∣|a|∣a∣, and Zn\mathbb Z^nZn is the set of integer vectors. For linearly independent b1,…,bm∈Rkb_1,\dots,b_m\in\mathcal R^kb1​,…,bm​∈Rk, the lattice L(b1,…,bm)L(b_1,\dots,b_m)L(b1​,…,bm​) is the set of integer combinations ∑jλjbj\sum_j\lambda_jb_j∑j​λj​bj​, λj∈Z\lambda_j\in\mathbb Zλj​∈Z, and b1,…,bmb_1,\dots,b_mb1​,…,bm​ is a basis. Gram–Schmidt orthogonalisation gives b1∗,…,bm∗b_1^*,\dots,b_m^*b1∗​,…,bm∗​ and unit vectors uj=bj∗/∣bj∗∣u_j=b_j^*/|b_j^*|uj​=bj∗​/∣bj∗​∣, and bi(j)=(bi,uj)b_i(j)=(b_i,u_j)bi​(j)=(bi​,uj​), so bi=∑jbi(j)ujb_i=\sum_jb_i(j)u_jbi​=∑j​bi​(j)uj​ and bj(j)=∣bj∗∣b_j(j)=|b_j^*|bj​(j)=∣bj∗​∣. The determinant is d(L)=∏j∣bj∗∣d(L)=\prod_j|b_j^*|d(L)=∏j​∣bj∗​∣. Λ1(L)\Lambda_1(L)Λ1​(L) is the length of a shortest nonzero vector of LLL. The projected lattice Lj(b1,…,bm)L_j(b_1,\dots,b_m)Lj​(b1​,…,bm​) is the image of LLL under orthogonal projection onto the complement of span⁡(b1,…,bj−1)\operatorname{span}(b_1,\dots,b_{j-1})span(b1​,…,bj−1​). A basis is reduced (Definition 2.6) if bj(j)=Λ1(Lj)b_j(j)=\Lambda_1(L_j)bj​(j)=Λ1​(Lj​) for every jjj and ∣bi(j)∣≤bj(j)/2|b_i(j)|\le b_j(j)/2∣bi​(j)∣≤bj​(j)/2 for i>ji>ji>j.

A convex body is a convex set of positive volume, which for a convex set means nonempty interior. A subspace VVV has a basis of integer vectors if it is the real span of integer vectors. Its integer translates are the sets z+Vz+Vz+V with z∈Znz\in\mathbb Z^nz∈Zn.

In Lean, Rk\mathcal R^kRk is EuclideanSpace ℝ (Fin k), a basis is b : Fin m → EuclideanSpace ℝ (Fin k), the lattice is lattice b = Submodule.span ℤ (Set.range b), ∣bj∗∣|b_j^*|∣bj∗​∣ is gsLen b j, bi(j)b_i(j)bi​(j) is gsCoeff b i j, d(L)d(L)d(L) is latticeDet b, Λ1\Lambda_1Λ1​ is lambdaOne, Lj+1L_{j+1}Lj+1​ is projLattice b j, and a reduced basis is IsReduced b.

Formalization targets

Goal: Theorem (5.5), corrected reading

For n≥2n\ge2n≥2 and every bounded convex set K⊆RnK\subseteq\mathcal R^nK⊆Rn with nonempty interior and K∩Zn=∅K\cap\mathbb Z^n=\emptysetK∩Zn=∅ there is a subspace VVV spanned by integer vectors with 1≤dim⁡V≤n−11\le\dim V\le n-11≤dimV≤n−1 and

#{ z+V:z∈Zn, (z+V)∩K≠∅ } ≤ n2(n−dim⁡V).\#\{\,z+V : z\in\mathbb Z^n,\ (z+V)\cap K\ne\emptyset\,\}\ \le\ n^{2(n-\dim V)}.#{z+V:z∈Zn, (z+V)∩K=∅} ≤ n2(n−dimV).

The printed theorem allows "an iii dimensional space VVV" with 1≤i≤n1\le i\le n1≤i≤n and bound n2(n−i+1)n^{2(n-i+1)}n2(n−i+1). Taken literally that is trivial (V=RnV=\mathcal R^nV=Rn, one translate). The proof on the same page takes V=span⁡(b1,…,bi−1)V=\operatorname{span}(b_1,\dots,b_{i-1})V=span(b1​,…,bi−1​), of dimension i−1i-1i−1, and remarks that this "ensures that the subspace VVV is always of dimension at least 1". The goal states that reading.

Milestones

  1. Theorem (1.11), Minkowski's convex body theorem (referenced from the platform, in Mathlib's general form).
  2. Theorem (1.12): every mmm-dimensional lattice has a nonzero vector with ∣v∣≤m d(L)1/m|v|\le\sqrt m\,d(L)^{1/m}∣v∣≤m​d(L)1/m.
  3. Proposition 1.9: a primitive lattice vector belongs to some basis.
  4. Proposition 2.16, existence form: every lattice has a reduced basis.
  5. Proposition 4.2: for any b0b_0b0​ with projection bˉ0\bar b_0bˉ0​ onto the span, some b∈Lb\in Lb∈L has ∣b−bˉ0∣≤12(∑jbj(j)2)1/2≤m2max⁡jbj(j)|b-\bar b_0|\le\frac12(\sum_jb_j(j)^2)^{1/2}\le\frac{\sqrt m}2\max_jb_j(j)∣b−bˉ0​∣≤21​(∑j​bj​(j)2)1/2≤2m​​maxj​bj​(j).
  6. Proposition 4.3: for a reduced basis and iii maximising bi(i)b_i(i)bi​(i), the tail (λi,…,λm)(\lambda_i,\dots,\lambda_m)(λi​,…,λm​) of every closest lattice point to b0b_0b0​ lies in an explicit set of at most mm−i+1m^{m-i+1}mm−i+1 integer vectors.

Significance

Theorem (5.5) is a structural form of the flatness theorem. For dim⁡V=n−1\dim V=n-1dimV=n−1 it says that a lattice-free convex body meets fewer than n2n^2n2 consecutive integer hyperplanes of some integral direction. For smaller dim⁡V\dim VdimV it gives a finer decomposition of Zn\mathbb Z^nZn into translates, each a lower-dimensional integer program. This is the recursion behind fixed-dimension integer programming, and statements of this form are used in lattice-point enumeration, in the geometry of numbers (covering minima), and in cutting-plane theory (lattice-free bodies define split and intersection cuts). Propositions 4.2 and 4.3 are the correctness core of exact closest-vector enumeration.

The results are proved in the literature, though Theorem (5.5) is proved "albeit sketchily" in the paper itself. As far as is known, none of them is formalized: Mathlib has Minkowski's convex body theorem and the ZLattice API, but not Gram–Schmidt lattice invariants, Korkine–Zolotarev bases, Hermite-type bounds, nearest-plane rounding, or any flatness theorem. This mission produces the first machine-checked versions. It also corrects three statements that are wrong as printed (below), so the formal statements are the ones that can be relied on.

Difficulty

The naive route to (5.5) is to take a flat direction directly: bound the lattice width of KKK and count hyperplanes. That needs a flatness theorem with an explicit bound below n2n^2n2, which is itself the hard part. The paper's argument instead needs John's theorem (every convex body lies between an ellipsoid and its nnn-fold dilation), a reduced basis of the transformed lattice, and a counting argument across projected lattices that combines Minkowski's bound on each LiL_iLi​ with the covering estimate of Proposition 4.2. None of John's theorem, reduced bases or the projected-lattice counting is in Mathlib.

Proposition 4.3 is also delicate as printed: the per-coordinate count on p. 24 undercounts the integers in a closed interval, so the printed arithmetic cannot be transcribed as it stands. Proposition 2.16 in the paper is the correctness of the algorithm SHORTEST. Here only the existence of a reduced basis is needed, which requires attainment of Λ1\Lambda_1Λ1​ on every projected lattice and a lifting argument (Proposition 1.9).

Formalization scope

Conventions: indices are 0-based (Fin m), so the paper's LjL_jLj​ is projLattice b (j-1) and its bound nn−i+1n^{n-i+1}nn−i+1 is m ^ (m - i). Gram–Schmidt is Mathlib's unnormalised gramSchmidt. Lattices are Submodule ℤs of a real Euclidean space generated by a linearly independent family, and m≤km\le km≤k is allowed, because (1.12) and 4.2 are applied to projected lattices. The goal counts translates as sets with Set.encard, so the bound includes finiteness. KKK is assumed convex, bounded and with nonempty interior, but not closed.

Three printed statements are corrected, and the corrections are recorded in each item's Formalization Note.

  • (1.12)'s constant 12n\frac12\sqrt n21​n​ is false for n≤7n\le7n≤7 (for example L=ZL=\mathbb ZL=Z, or the hexagonal lattice) and is replaced by n\sqrt nn​, the constant the paper's own later proofs use.
  • Proposition 4.2's second sentence is stated for bˉ0\bar b_0bˉ0​ instead of b0b_0b0​.
  • Proposition 4.3 fails at n=1n=1n=1 and is stated for m≥2m\ge2m≥2 with the proof's explicit candidate set TTT, since an existential TTT is satisfied by the set of tails of closest points and says nothing.

Trivializing formalizations are ruled out: the goal forbids dim⁡V=n\dim V=ndimV=n, which gives one translate, and dim⁡V=0\dim V=0dimV=0, where no translate meets KKK. It requires nonempty interior (the empty set would satisfy everything) and counts with encard (an infinite count cannot become 000).

Out of scope: the paper's algorithms (SHORTEST, SELECT-BASIS, ENUMERATE, CLP, CLP′, ILP) and their operation and bit counts (Theorems 2.17, 3.9, 4.5, 5.4), because Mathlib has no cost model. Also out of scope is §6 (NP-completeness of the L2L_2L2​ closest vector problem and Cook reductions), because Mathlib has no complexity classes. The definitions of this mission (lattice, gsLen, gsCoeff, latticeDet, lambdaOne, projLattice, IsReduced) are reusable for any later work on lattice reduction. Contributions are welcome at every level: John's theorem, Hermite-type bounds, Korkine–Zolotarev existence, and the counting lemmas.

Selected references

  • R. Kannan, Minkowski's Convex Body Theorem and Integer Programming, Mathematics of Operations Research 12(3):415–440, 1987. https://doi.org/10.1287/moor.12.3.415
  • H. W. Lenstra Jr., Integer programming with a fixed number of variables, Mathematics of Operations Research 8(4):538–548, 1983. https://doi.org/10.1287/moor.8.4.538
  • R. Kannan, L. Lovász, Covering minima and lattice-point-free convex bodies, Annals of Mathematics 128(3):577–602, 1988. https://doi.org/10.2307/1971436
  • A. K. Lenstra, H. W. Lenstra Jr., L. Lovász, Factoring polynomials with rational coefficients, Mathematische Annalen 261:515–534, 1982. https://doi.org/10.1007/BF01457454
  • F. John, Extremum problems with inequalities as subsidiary conditions, Studies and Essays presented to R. Courant, 1948, 187–204.
9 thms3 active usersReviewed
Captain: Lucas

The moving sofa problem: optimality of Gerver's sofaResearch Paper

Motivation

The moving sofa problem, posed by Leo Moser in 1966, asks for the planar shape of largest area that can be moved around a right-angled corner in a hallway of unit width — one of the most famous unsolved problems of everyday geometry. Joseph Gerver exhibited in 1992 a sofa of area 2.2195... built from five boundary curves solving a system of balancing ordinary differential equations, and conjectured it to be the largest possible; Dan Romik later derived the curves' exact formulas, and Kallus and Romik proved the best upper bound before the source. The source resolves the problem: Baek proves that Gerver's sofa is optimal, so the sofa constant is Gerver's area 2.2195..., and the sixty-year-old problem is closed.

The mission formalizes the proof of Theorem 1.1.1 of the source — Gerver's sofa attains the maximum area among all moving sofas — following the paper's own attack path: reduce to monotone sofas and their supporting rotation paths, express the sofa area through the cap and niche of the path, control the rotation angle of large sofas, take Hausdorff limits of maximum polygon sofas to obtain balanced maxima, establish the injectivity of the supporting path of a maximizer, and dominate every injectivity-constrained path by a Gerver sofa through the concavity of the area upper bound Q.

The fragment and its conventions

The arena is the L-shaped hallway L, the union of the horizontal side (-∞, 1] × [0, 1] and the vertical side [0, 1] × (-∞, 1]. A shape is a nonempty compact subset of the plane — the source's standing convention, adopted here — and a moving sofa is a nonempty, compact, connected shape that admits a movement around the corner. A movement is encoded by two continuous paths: a rotation angle path Θ and a translation path W, producing the rigid motion Φ_t(p) = rot(Θ t) p + W t of the unit interval, with Θ 0 = 0, the initial image inside the horizontal side, every intermediate image inside the hallway, the final image inside the vertical side, and total clockwise rotation -Θ 1 = ω, the movement's rotation angle. The area of a planar set is its two-dimensional Lebesgue measure as a real number; the maximum sofa area α_max is the supremum of the areas of all moving sofas.

The monotone arena replaces sofas by their supporting rotation paths. A rotation path x with angle ω ∈ (0, π/2] determines strips, quarter-planes Q_t⁺ \ Q_t⁻ at the outer and inner corners of the translated hallways, the parallelogram P_ω, the fan F_ω, the cap (Equation 1.3), the niche (Equation 1.4), and the monotone sofa of the path — the common part of the strips and the hallways (Equation 1.2). A supporting path pins down the canonical inner-corner path through contact: both outer walls of every supporting hallway touch the sofa. Balanced maximum sofas are stated at the level of sofas: a moving sofa that is the Hausdorff limit of maximum polygon sofas whose uniform angle sets run over a strictly increasing sequence of powers of two, with Hausdorff convergence phrased through Euclidean ε-neighborhoods rather than Mathlib's metric thickening API. The Gerver sofa is characterized by the paper's own structural properties: a moving sofa equal to the monotone sofa of its rotation path, whose supporting-hallway structure is carried by the five curves of Definition 8.4.3 with contact on the outer and inner walls and satisfying the balancing ODEs of Theorem 8.4.2 on the five intervals cut by angles 0 < φ < θ < π/4, derivatives taken one-sided within each interval. The source defers the curves' explicit formulas to Gerver's and Romik's papers; the existential form — some sofa satisfying this structure — is used for the goal and the dominating bound, since that is what the proof of Theorem 1.1.1 actually produces.

Six conventions deserve disclosure. First, the platform publishes definitions as modules Definitions.Def_<definition_name>, so the three definition bundles — the hallway arena, the monotone arena, and the Gerver sofa structure with Mamikon's swept region — are published as platform definitions (Gerver_arena, Gerver_monotone, Gerver_sofa) and attached to the proposal as reference items; the theorem drafts import them. Second, the injectivity condition abridges the source's Definition 6.1.2 by dropping its auxiliary clause on the boundary curve σ_K, keeping the continuous differentiability and the two strict derivative inequalities that the endgame consumes. Third, balancedness is stated at the sofa level rather than through the paper's rotation-path data, the sofa-level form being what the compactness argument yields. Fourth, the generalized Mamikon theorem is stated for continuously differentiable tangent segment families with continuous lengths and pairwise non-overlap on the open parameter rectangle, over turning intervals of length below π — exactly the smooth, non-self-overlapping families to which the endgame applies it. Fifth, the movement's initial term is a pure translation (Θ 0 = 0), matching the source's Definition 1.1.2 with its footnote. Sixth, all statements are developed against the platform's Mathlib build (0df444a), with inner products, rotations, scalings, and neighborhoods given by small self-contained definitions so that every draft compiles without instance-hunting.

Formalization targets

The goal is Theorem 1.1.1: some Gerver sofa has area at least that of every moving sofa. The supporting chain, in the paper's order, each step a milestone except the last:

  • Monotone reduction: every moving sofa is contained in a monotone moving sofa of at least its area — the intersection of the supporting hallways.
  • Cap-niche identity: a moving monotone sofa is its cap minus its niche.
  • Area functional: the sofa area is the cap area minus the niche area whenever the niche lies in the cap.
  • Angle control: a moving sofa of area at least 2.2 admits a movement of rotation angle in [arccos(1/2.2), π/2].
  • Balanced maxima: for each rotation angle, a balanced maximum sofa exists as the Hausdorff limit of maximum polygon sofas.
  • Right-angle reduction: a balanced maximum sofa of large area has a rotated copy admitting a movement of angle π/2.
  • Injectivity: the supporting rotation path of a maximum-area balanced sofa of angle π/2 satisfies the injectivity condition.
  • Mamikon: the area swept by a smooth non-overlapping tangent segment family is half the integral of the squared lengths.
  • Dominating bound: some Gerver sofa bounds the area of every monotone sofa of an injectivity-constrained supporting path of angle π/2.

The goal follows by assembling the chain: every moving sofa sits under a monotone one, maximizers are balanced maxima admitting the full right angle, their supporting paths are injectivity-constrained, and the dominating bound caps them by a Gerver sofa.

Significance

The moving sofa problem is a benchmark of planar geometry that has resisted exact solution since 1966; the source closes it. The proof is not a computational verification but a structural argument — compactness produces extremal objects, contact geometry pins their supporting paths, injectivity rules out self-overlap, and a concavity argument of a kind familiar from calculus of variations delivers the dominating bound. Each layer is a theorem of lasting value independently of the goal: the cap-niche decomposition, the existence of balanced maxima, the injectivity of maximizer paths, and the generalized Mamikon theorem all predate or outlive the final step. A Lean formalization makes each layer independently checkable, and the mission's chain structure lets solvers claim credit for any layer while the community assembles the whole.

Difficulty

The formalization spans real analysis, topology, convex geometry, and measure theory, at a level well beyond elementary competition material but below the frontier of ongoing formalization projects. The monotone reduction and the cap-niche identity require careful handling of intersections of translated half-planes and the additivity of Lebesgue measure; the balanced-maximum existence requires sequential compactness in the Hausdorff topology; injectivity requires differentiability arguments on paths defined by contact conditions; the endgame requires the concavity of the functional Q on triples of convex bodies. The definitions are engineered to be self-contained — plain real arithmetic for inner products, rotations, and neighborhoods — so that the labor concentrates on the geometry rather than on API friction. Individual milestones are of moderate difficulty; the goal requires assembling all of them, which is the source's own achievement.

Formalization scope

Every item quantifies only over the objects it needs — hallways, shapes, rotation paths, angle sets, segment families; no item carries a vacuous hypothesis: moving sofas, monotone sofas, polygon intersections, and Gerver sofas all exist, and the inequalities are in the intended directions. The fragment excludes the explicit trigonometric formulas of Gerver's curves (the goal needs only their structural characterization), the source's sections on the ambidextrous sofa and non-convex generalizations, the computer-search evidence for the sofa constant's decimals, and any statement about shapes that are not moving sofas. Nothing in the mission asserts the uniqueness of the Gerver sofa or the value of the sofa constant as a decimal; it formalizes what the source proves — optimality.

Selected references

  • J. Baek, Optimality of Gerver's Sofa, arXiv:2411.19826 (2024) — the mission's source.
  • L. Moser, Problem 66-11, Moving Furniture Through a Hallway, SIAM Review 8 (1966), p. 381 — the problem.
  • J. L. Gerver, On Moving a Sofa around a Corner, Geometriae Dedicata 42 (1992), pp. 267–283 — the sofa and the balancing ODEs.
  • D. Romik, Differential Equations and Exact Solutions in the Moving Sofa Problem, Experimental Mathematics 27 (2018), pp. 316–330 — the curves' exact formulas.
  • Y. Kallus and D. Romik, Improved Upper Bounds in the Moving Sofa Problem, Advances in Mathematics 340 (2018), pp. 960–982 — the preceding upper bound.
13 thms2 active usersReviewed
CombinatoricsGraph Theory·Captain: mikedeng1

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

Motivation

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

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

Setting

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

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

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

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

Formalization targets

Goal: Theorem 4.5 (p. 3:15)

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

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

Milestones, in the order the paper uses them

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

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

Significance

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

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

Difficulty

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

Formalization scope

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

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

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

Selected references

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

Slope-field perturbations of the two-cylinder coveringOpen Problem

Motivation

A slope-field perturbation of a two-cylinder cover trades a reduction in projected area against the extra overlap needed to maintain coverage. The slope-field manuscript studies when the net cost becomes strictly smaller. The companion radial-sweep manuscript establishes finite triangular approximation with arbitrary tagged equal subdivisions. Both manuscripts map to the same published Lean goal, so this draft covers their shared target. The source manuscript presents the surrounding research claim.

Setting

A cylinder is a measurable finite-area planar base plus a line in its perpendicular unit direction. Ruled sets vary the segment direction over compact planar labels. Radial sweeps use p₀+ρ(a+tb)+s(e+w(t)), with 0≤ρ≤R(t), |s|≤L, and w′(t)=λ(t)(a+tb).

Formalization target

Prove the seven-part fourPaperMain conjunction: the explicit angular-area construction; C¹ hyperbolic approximation; C¹ square-zero-derivative approximation; regular-tetrahedron counterexamples with parallelogram bases and with bounded open bases; aligned radial triangular approximation; and actual-sweep projection bounds. In particular, the tetrahedron clauses require

∑i∣Bi∣<12Amin⁡(K),∑i∣Bi∣∣πui⊥K∣<12.\sum_i|B_i|<\frac12 A_{\min}(K),\qquad\sum_i\frac{|B_i|}{|\pi_{u_i^\perp}K|}<\frac12.i∑​∣Bi​∣<21​Amin​(K),i∑​∣πui⊥​​K∣∣Bi​∣​<21​.

The selected formal target is OAI.CylinderCovering.fourPaperMain.

Significance and status

The selected published goal is a joint endpoint shared by four related manuscripts, and therefore includes more than the slope-field application alone. Its finite covers retain both cost inequalities. The published target is currently Open; the manuscript's argument and a checked Lean proof are separate deliverables.

Difficulty

A local area decrease cannot guarantee a global cover. All regularity, exact coverage and error bounds in the separate conjuncts must hold in their stated regimes.

Formalization scope

Hyperbolic and nilpotent clauses use their C¹ definitions; radial data require a nonempty parameter interval, positive segment half-length, continuous nonnegative cap and alignment derivative. All tagged equal subdivisions are covered. The radial clause requires triangular-base area to converge to weightedArea(D), while actual sweep projections have eventual cost at most weightedArea(D)+ε and the corresponding limsup bound. The separate RuledApproximation.fullMain is also retained as a related reference but is not mislabeled as one of these seven conjuncts.

Selected references

  • OpenAI, Finite triangular approximation of radial sweeps, preprint, 2026. Pinned manuscript.
  • OpenAI, Slope-field perturbations of the two-cylinder covering, preprint, 2026. Pinned manuscript.
  • OpenAI, accompanying formal statements, commit adc7f1241b42. Goal source.
30 thms2 active usersReviewed
CombinatoricsGraph Theory·Captain: wurtle

The Euclidean plane is not five-colorableResearch Paper

Motivation

The Hadwiger–Nelson problem asks for the chromatic number of the plane χ(R2)\chi(\mathbb R^2)χ(R2): the least number of colours needed to colour every point of the Euclidean plane so that any two points at distance exactly 111 receive different colours. Equivalently, it is the chromatic number of the infinite graph whose vertices are the points of the plane and whose edges are the unit-distance pairs. Since 1950 the answer has been known to lie between 444 and 777, and it is one of the best known open problems in combinatorial geometry. Its difficulty is that the colour classes may be arbitrary sets: no measurability or regularity is assumed, so analytic tools that work for measurable colourings do not apply directly.

This mission asks for a formal proof that five colours do not suffice, as claimed in an OpenAI preprint dated September 23, 2026 (source), so that 6≤χ(R2)≤76\le\chi(\mathbb R^2)\le76≤χ(R2)≤7. The preprint has not been peer reviewed and its claims have not been independently verified; on the platform the Lean goal is open.

Background

  • 1950 — Nelson poses the problem and observes χ≥4\chi\ge4χ≥4; Isbell finds a hexagonal 777-colouring (recounted by Soifer, Mathematics Competitions 2003).
  • 1951 — de Bruijn and Erdős: an infinite graph is kkk-colourable iff every finite subgraph is (Indag. Math. 1951).
  • 1961 — Hadwiger publishes the problem with bounds 444 and 777 (Elem. Math. 1961); the Moser spindle gives a 7-vertex obstruction to three colours (Canad. Math. Bull. 1961).
  • 1973–2005 — Woodall studies region colourings (JCTA 1973); Townsend proves map-type colourings need six colours (JCTA 1981; Geombinatorics 2005).
  • 1981 — Falconer proves that measurable colourings need at least five colours (JCTA 1981).
  • 2018–2020 — de Grey constructs a finite unit-distance graph that is not 444-colourable, so χ≥5\chi\ge5χ≥5 (arXiv:1804.02385); Heule and Parts shrink it (arXiv:1805.12181, arXiv:2010.12665); Exoo–Ismailescu give another proof (DCG 2020).
  • 2025 — Sokolov and Voronov prove seven colours are needed for polygonal map-type colourings (arXiv:2502.01958).
  • September 2026 — The OpenAI preprint claims χ(R2)≥6\chi(\mathbb R^2)\ge6χ(R2)≥6 for arbitrary colourings (Theorem 1.1, p. 2).

Setting

A proper kkk-colouring of the plane is any function c:R2→{1,…,k}c:\mathbb R^2\to\{1,\dots,k\}c:R2→{1,…,k} such that c(x)≠c(y)c(x)\ne c(y)c(x)=c(y) whenever ∥x−y∥=1\|x-y\|=1∥x−y∥=1, with ∥⋅∥\|\cdot\|∥⋅∥ the Euclidean norm. No measurability, continuity or regularity of the colour classes c−1(i)c^{-1}(i)c−1(i) is assumed. χ(R2)\chi(\mathbb R^2)χ(R2) is the least kkk admitting a proper kkk-colouring.

Formalization targets

Milestone: the upper bound χ(R2)≤7\chi(\mathbb R^2)\le7χ(R2)≤7 (Theorem 1.1, p. 2)

There is a proper 777-colouring of the plane: colour the Voronoi hexagons of a triangular lattice (circumradius r=2/5r=2/5r=2/5) by the seven cosets of an index-seven similar sublattice, assigning boundary points to any incident hexagon. This classical construction is reproved in the deduction of Theorem 1.1 (pp. 2–3), including all boundary points.

Goal: no proper five-colouring (Theorem 1.1, p. 2)

There is no function c:R2→{1,…,5}c:\mathbb R^2\to\{1,\dots,5\}c:R2→{1,…,5} with c(x)≠c(y)c(x)\ne c(y)c(x)=c(y) whenever ∥x−y∥=1\|x-y\|=1∥x−y∥=1. Together with the milestone this gives 6≤χ(R2)≤76\le\chi(\mathbb R^2)\le76≤χ(R2)≤7.

Significance

The result itself. It raises the lower bound for the chromatic number of the plane from five (de Grey, 2018) to six, leaving only 666 and 777. Unlike de Grey's bound it is not witnessed by a finite graph; the proof transfers the problem to measurable colourings for every number of colours (Theorem 1.3, p. 2) and then excludes weak measurable five-colourings (Theorem 1.4, p. 2). The transfer theorem is of independent interest: in ZFC, a proper kkk-colouring exists iff a weak measurable kkk-colouring does. The preprint also derives a positive lower bound on the invariant-mean frequency of monochromatic unit pairs for every five-colouring (Corollary 1.5, p. 5).

Formalizing it. The statements are elementary, but the proof uses ergodic theory (a Furstenberg–Zimmer compact-extension tower and rigidity of invariant measures on a character group, Theorem 2.3, p. 6), spectral theory, density points, Fourier decay of circle measure, and planar topology. A formal proof would remove any doubt about a long and varied argument. The finite-graph bound χ≥5\chi\ge5χ≥5 has been checked by SAT solvers, but no formal proof of a bound that is not witnessed by a finite graph is known.

Difficulty

By de Bruijn–Erdős, χ≥6\chi\ge6χ≥6 is equivalent to the existence of a finite unit-distance graph that is not 555-colourable, but no such graph is known, and SAT-based searches have not produced one. Measurable methods (Falconer's density-point argument) cannot be applied to arbitrary colour classes, and a measurable colouring can have boundaries far too irregular for map-type interface arguments such as Townsend's or Sokolov–Voronov's. The preprint must therefore (i) produce measurable data from an arbitrary colouring without losing the unit-distance exclusions, via invariant averaging over algebraic rotations and translations and the rigidity Theorem 2.3, and (ii) construct the connected interfaces needed for a geometric contradiction from measure estimates alone (Sections 5–7), before excluding label cycles of length 333, 444 and 555 (Section 8).

Formalization scope

  • The goal works in ℂ with ‖p - q‖ = 1; a colouring is any function ℂ → Fin 5. ProperColoring 5 c requires distinct colours at every unit-distance pair. No measurability is assumed, so the statement is about arbitrary colourings, as in the source.
  • The milestone works in EuclideanSpace ℝ (Fin 2) with dist x y = 1 and asks for some c : Plane → Fin 7; the two planes are isometric.
  • Both statements are in ZFC-style classical Lean; no choice-free reading is intended. Neither is vacuous.

Welcome contributions: the de Bruijn–Erdős compactness theorem, the Moser spindle, Lebesgue density points in R2\mathbb R^2R2, and the measurable/unrestricted transfer of Theorem 1.3.

Selected references

  • OpenAI, The Euclidean plane is not five-colorable, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-Euclidean-plane-is-not-five-colorable-September-23-2026/paper.pdf
  • H. Hadwiger, Ungelöste Probleme Nr. 40, Elem. Math. (1961).
  • L. Moser, W. Moser, Solution to Problem 10, Canad. Math. Bull. (1961). https://doi.org/10.1017/S0008439500025765
  • N. G. de Bruijn, P. Erdős, A colour problem for infinite graphs and a problem in the theory of relations, Indag. Math. (1951). https://users.renyi.hu/~p_erdos/1951-01.pdf
  • K. J. Falconer, The realization of distances in measurable subsets covering Rn\mathbb R^nRn, J. Combin. Theory Ser. A (1981). https://doi.org/10.1016/0097-3165(81)90014-5
  • D. R. Woodall, Distances realized by sets covering the plane, J. Combin. Theory Ser. A (1973). https://doi.org/10.1016/0097-3165(73)90020-4
  • S. P. Townsend, Every 5-coloured map in the plane contains a monochrome unit, J. Combin. Theory Ser. A (1981). https://doi.org/10.1016/0097-3165(81)90046-7
  • A. D. N. J. de Grey, The chromatic number of the plane is at least 5, Geombinatorics (2018). https://arxiv.org/abs/1804.02385v3
  • M. J. H. Heule, Computing small unit-distance graphs with chromatic number 5, Geombinatorics (2018). https://arxiv.org/abs/1805.12181
  • G. Exoo, D. Ismailescu, The chromatic number of the plane is at least 5: a new proof, Discrete Comput. Geom. (2020). https://doi.org/10.1007/s00454-019-00058-1
  • G. Sokolov, V. Voronov, On the chromatic number of the plane for map-type colorings, arXiv:2502.01958 (2025). https://arxiv.org/abs/2502.01958v1
  • A. Soifer, The 50th anniversary of one problem: the chromatic number of the plane & its relatives, Mathematics Competitions (2003). https://www.wfnmc.org/Journal%202003%201.pdf
10 thms2 active usersReviewed
Harmonic AnalysisProbability·Captain: mikedeng1

On Lattices, Learning with Errors, Random Linear Codes, and Cryptography 2: For r ≥ √2·η_ε(L), ε ≤ 1/10, a Discrete Gaussian Sample Lies Outside Any Proper Subspace with Probability ≥ 1/10Research Paper

Motivation

Regev's paper (J. ACM 56(6), Article 34, 2009) introduces the learning with errors (LWE) problem and shows that solving it is at least as hard, for quantum algorithms, as approximating standard worst-case lattice problems such as the shortest independent vectors problem (SIVP) and GapSVP. LWE is the basis of most lattice-based public-key cryptography, including the post-quantum schemes standardized since.

The reduction runs through a sampling problem: produce samples of the discrete Gaussian distribution DL,rD_{L,r}DL,r​ on a lattice LLL. The last step converts such a sampler into a solver for SIVP, which asks for nnn short linearly independent lattice vectors. That conversion rests on one geometric fact about DL,rD_{L,r}DL,r​, Lemma 3.15 of the paper (p. 34:31): when the width rrr is at least 2\sqrt22​ times the smoothing parameter of LLL, a sample is not trapped in any proper subspace. The lemma first appeared, with only a proof sketch, in the preliminary version of Micciancio and Regev's work on worst-case to average-case reductions (SIAM J. Comput. 2007); Regev's article gives the full proof. This mission formalizes that lemma. The paper's main theorem, a quantum reduction (Theorem 3.1), is outside this series.

Setting

Work in Rn\mathbb{R}^nRn, n≥1n \ge 1n≥1, with the Euclidean inner product ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩ and norm ∥⋅∥\|\cdot\|∥⋅∥.

  • A lattice L⊂RnL \subset \mathbb{R}^nL⊂Rn is the set of integer combinations of nnn linearly independent vectors. Its dual is L∗={y∈Rn:⟨x,y⟩∈Z ∀x∈L}L^* = \{y \in \mathbb{R}^n : \langle x, y\rangle \in \mathbb{Z}\ \forall x \in L\}L∗={y∈Rn:⟨x,y⟩∈Z ∀x∈L}, and det⁡(L)\det(L)det(L) is the volume of a fundamental domain; det⁡(L∗)=1/det⁡(L)\det(L^*) = 1/\det(L)det(L∗)=1/det(L).
  • The Gaussian function of width s>0s > 0s>0 is ρs(x)=exp⁡(−π∥x∥2/s2)\rho_s(x) = \exp(-\pi\|x\|^2/s^2)ρs​(x)=exp(−π∥x∥2/s2) (Eq. (4)); for a countable set AAA, ρs(A)=∑x∈Aρs(x)\rho_s(A) = \sum_{x \in A}\rho_s(x)ρs​(A)=∑x∈A​ρs​(x).
  • The discrete Gaussian DL,rD_{L,r}DL,r​ is the probability distribution on LLL with mass DL,r(x)=ρr(x)/ρr(L)D_{L,r}(x) = \rho_r(x)/\rho_r(L)DL,r​(x)=ρr​(x)/ρr​(L) (Eq. (6)).
  • For ϵ>0\epsilon > 0ϵ>0, the smoothing parameter ηϵ(L)\eta_\epsilon(L)ηϵ​(L) is the smallest s>0s > 0s>0 with ρ1/s(L∗∖{0})≤ϵ\rho_{1/s}(L^* \setminus \{0\}) \le \epsilonρ1/s​(L∗∖{0})≤ϵ (Definition 2.10, p. 34:19).
  • A subspace HHH is a real linear subspace of Rn\mathbb{R}^nRn; "dimension at most n−1n-1n−1" means HHH is proper.

Formalization targets

Goal: Lemma 3.15 (p. 34:31)

For an nnn-dimensional lattice LLL, 0<ϵ≤1100 < \epsilon \le \tfrac1{10}0<ϵ≤101​, r≥2 ηϵ(L)r \ge \sqrt2\,\eta_\epsilon(L)r≥2​ηϵ​(L) and every subspace HHH with dim⁡H≤n−1\dim H \le n - 1dimH≤n−1,

Pr⁡x∼DL,r[x∉H]=∑x∈L∖Hρr(x)ρr(L)  ≥  110.\Pr_{x \sim D_{L,r}}[x \notin H] = \sum_{x \in L \setminus H} \frac{\rho_r(x)}{\rho_r(L)} \;\ge\; \frac{1}{10}.x∼DL,r​Pr​[x∈/H]=x∈L∖H∑​ρr​(L)ρr​(x)​≥101​.

The constants 2\sqrt22​, 110\tfrac1{10}101​ and 110\tfrac1{10}101​ are the printed ones.

Milestones, in attack order

  1. Lemma 2.14 (Poisson summation, p. 34:20), already proved on the platform: f(L)=det⁡(L∗)f^(L∗)f(L) = \det(L^*)\hat f(L^*)f(L)=det(L∗)f^​(L∗), with f^(w)=∫f(x)e−2πi⟨x,w⟩dx\hat f(w) = \int f(x)e^{-2\pi i\langle x, w\rangle}dxf^​(w)=∫f(x)e−2πi⟨x,w⟩dx.
  2. The remark after Definition 2.10 (p. 34:19): s↦ρ1/s(L∗∖{0})s \mapsto \rho_{1/s}(L^* \setminus \{0\})s↦ρ1/s​(L∗∖{0}) is continuous and strictly decreasing on (0,∞)(0,\infty)(0,∞), tends to ∞\infty∞ at 000 and to 000 at ∞\infty∞, and ϵ↦ηϵ(L)\epsilon \mapsto \eta_\epsilon(L)ϵ↦ηϵ​(L) is its inverse.
  3. The expectation chain (proof of Lemma 3.15, p. 34:31): for a unit vector www,
Expx∼DL,r[e−π⟨w,x⟩2/r2]=det⁡(L∗)rn2 ρr(L)∑y∈L∗e−πr2(∥y∥2−⟨w,y⟩2/2)≤det⁡(L∗)rn2 ρr(L)ρ2/r(L∗)≤det⁡(L∗)rn2 ρr(L)(1+ϵ).\mathop{\mathrm{Exp}}_{x\sim D_{L,r}}\bigl[e^{-\pi\langle w,x\rangle^2/r^2}\bigr] = \frac{\det(L^*)r^n}{\sqrt2\,\rho_r(L)}\sum_{y\in L^*}e^{-\pi r^2(\|y\|^2 - \langle w,y\rangle^2/2)} \le \frac{\det(L^*)r^n}{\sqrt2\,\rho_r(L)}\rho_{\sqrt2/r}(L^*) \le \frac{\det(L^*)r^n}{\sqrt2\,\rho_r(L)}(1+\epsilon).Expx∼DL,r​​[e−π⟨w,x⟩2/r2]=2​ρr​(L)det(L∗)rn​y∈L∗∑​e−πr2(∥y∥2−⟨w,y⟩2/2)≤2​ρr​(L)det(L∗)rn​ρ2​/r​(L∗)≤2​ρr​(L)det(L∗)rn​(1+ϵ).
  1. The normalization bound (p. 34:31): ρr(L)=det⁡(L∗)rnρ1/r(L∗)≥det⁡(L∗)rn\rho_r(L) = \det(L^*)r^n\rho_{1/r}(L^*) \ge \det(L^*)r^nρr​(L)=det(L∗)rnρ1/r​(L∗)≥det(L∗)rn for every r>0r > 0r>0.
  2. The conclusion of the proof (p. 34:31): the expectation is at most 12(1+ϵ)<0.9\tfrac1{\sqrt2}(1+\epsilon) < 0.92​1​(1+ϵ)<0.9.

Significance

The result. Lemma 3.15 gives Corollary 3.16: n2n^2n2 independent samples of DL,rD_{L,r}DL,r​ contain nnn linearly independent vectors except with exponentially small probability. Since samples of DL,rD_{L,r}DL,r​ have norm at most about n r\sqrt n\, rn​r with high probability, a sampler at width r≈2 ηϵ(L)r \approx \sqrt2\,\eta_\epsilon(L)r≈2​ηϵ​(L) therefore yields short independent vectors, which is how Regev's reduction reaches SIVP (§3.3). The same property is used in later work wherever discrete Gaussian samples must generate a full-rank set of lattice vectors.

Formalizing it. The lemma is proved in the paper; it has no machine-checked proof that we know of. The work this mission asks for is the formalization of the known proof: the Gaussian-sum analysis of a discrete Gaussian on a general lattice in Mathlib's ZLattice framework, built on the platform's proved Poisson summation formula for lattices. The remark after Definition 2.10 (well-definedness of the smoothing parameter) and the normalization bound ρr(L)≥det⁡(L∗)rn\rho_r(L) \ge \det(L^*)r^nρr​(L)≥det(L∗)rn are used across the whole paper (Lemmas 3.2, 3.12, Claims 3.8, 3.9) and are reusable beyond this mission.

Difficulty

The obvious argument fails at the first step: there is no direct way to bound the mass of DL,rD_{L,r}DL,r​ on a hyperplane, because a lattice can place many short vectors in one hyperplane and the normalization ρr(L)\rho_r(L)ρr​(L) is not known in closed form. The paper's route replaces the indicator of HHH by a Gaussian in the direction orthogonal to HHH and moves to the dual lattice, where the smoothing hypothesis controls the sum. Formally this requires Poisson summation for an anisotropic Gaussian (a Gaussian with a non-scalar covariance), i.e. the Fourier transform of x↦exp⁡(−π(∥x∥2+⟨w,x⟩2)/r2)x \mapsto \exp(-\pi(\|x\|^2 + \langle w,x\rangle^2)/r^2)x↦exp(−π(∥x∥2+⟨w,x⟩2)/r2), which is not a stated Mathlib lemma, together with summability of Gaussian sums over lattices and the identification of the dual lattice used by the Poisson formula with the one in the paper. The smoothing remark needs continuity and strict monotonicity of a lattice theta series and its limits at 000 and ∞\infty∞.

Formalization scope

  • Rn\mathbb{R}^nRn is EuclideanSpace ℝ (Fin n) with Lebesgue measure. A lattice is L : Submodule ℤ (EuclideanSpace ℝ (Fin n)) with [DiscreteTopology L] [IsZLattice ℝ L] (discrete, full rank), the same shape as the referenced Poisson summation item. L∗L^*L∗ is LinearMap.BilinForm.dualSubmodule (innerₗ _) L; det⁡(L∗)\det(L^*)det(L∗) is written (ZLattice.covolume L)−1(\texttt{ZLattice.covolume } L)^{-1}(ZLattice.covolume L)−1.
  • All Gaussian sums are real tsums over the lattice or its dual; they are summable, so the values are the true sums. DL,rD_{L,r}DL,r​ is a mass function on the points of LLL, and the probability of x∉Hx \notin Hx∈/H is the sum of the masses of the lattice points outside HHH.
  • ηϵ(L)\eta_\epsilon(L)ηϵ​(L) is the real infimum of {s>0:ρ1/s(L∗∖{0})≤ϵ}\{s > 0 : \rho_{1/s}(L^* \setminus \{0\}) \le \epsilon\}{s>0:ρ1/s​(L∗∖{0})≤ϵ}; that it is attained is milestone 2, not a definitional shortcut.
  • Added hypotheses, all disclosed in the items: n≥1n \ge 1n≥1 in the goal (at n=0n = 0n=0 the natural-number subtraction n−1n - 1n−1 is 000, the only subspace of dimension 000 is the whole space R0\mathbb{R}^0R0, and the probability is 000); ϵ>0\epsilon > 0ϵ>0 (Definition 2.10 is for ϵ>0\epsilon > 0ϵ>0; for ϵ≤0\epsilon \le 0ϵ≤0 the Lean infimum is over the empty set, equals 000, and the hypothesis r≥2 ηϵ(L)r \ge \sqrt2\,\eta_\epsilon(L)r≥2​ηϵ​(L) would be vacuous); r>0r > 0r>0 (a Gaussian width; for n≥1n \ge 1n≥1 it follows from the other hypotheses).
  • Rotation. The proof's "assume without loss of generality that (1,0,…,0)(1,0,\dots,0)(1,0,…,0) is orthogonal to HHH" is formalized coordinate-free: milestones 3 and 5 are stated for an arbitrary unit vector www, with x1x_1x1​ replaced by ⟨w,x⟩\langle w, x\rangle⟨w,x⟩.
  • Lemma 2.14 is printed "for any function fff", the technical conditions omitted (p. 34:20). The referenced platform item states it for Schwartz functions fff, in translated form; every function the proof applies it to is a Gaussian, hence Schwartz.
  • Not a trivializing formalization: the probability is not stated through a Bochner integral or an unconstrained tsum, the smoothing hypothesis is not vacuous (ϵ>0\epsilon > 0ϵ>0), and the goal is the probability bound, not the expectation bound of milestone 3.
  • Contributions welcome: Fourier transforms of anisotropic Gaussians on EuclideanSpace, summability and monotonicity of lattice theta series, and the smoothing remark, each reusable in the sibling missions of this series.

Selected references

  • O. Regev, On Lattices, Learning with Errors, Random Linear Codes, and Cryptography, Journal of the ACM 56(6), Article 34, 2009. https://doi.org/10.1145/1568318.1568324
  • D. Micciancio and O. Regev, Worst-Case to Average-Case Reductions Based on Gaussian Measures, SIAM Journal on Computing 37(1), 2007. https://doi.org/10.1137/S0097539705447360
9 thms2 active usersReviewed
Linear OptimizationNumber TheoryOperations Research·Captain: mikedeng1

Integer Programming with a Fixed Number of Variables 1: If B(p,r) ⊂ τK ⊂ B(p,R) with R/r ≤ c₁ and τK Misses the Lattice, Fewer Than 1 + c₁c₂√n Lattice Hyperplanes H + kbₙ Meet B(p,R)Research Paper

Motivation

The integer linear programming feasibility problem asks, for an integer m×nm\times nm×n matrix AAA and b∈Zmb\in\mathbb Z^mb∈Zm, whether some x∈Znx\in\mathbb Z^nx∈Zn satisfies Ax≤bAx\le bAx≤b. It is NP-complete, so no algorithm polynomial in the length of the data is expected. H. W. Lenstra, Jr. (Math. Oper. Res. 8 (1983) 538–548) showed that the picture changes when the number of variables nnn is fixed: the problem is then solvable in polynomial time. The proof is a short argument from the geometry of numbers, and it became the starting point of lattice methods in integer programming, including Kannan's improved algorithm and the later flatness-based branching schemes.

Lenstra's introduction (p. 538) describes the idea as transforming the problem into an equivalent one in which "either the existence of a vector x∈Znx \in \mathbb Z^nx∈Zn satisfying Ax≤bAx \le bAx≤b is obvious; or it is known that the last coordinate of any such xxx belongs to an interval whose length is bounded by a constant only depending on nnn." This mission formalizes the geometric statement behind that dichotomy, §1 of the paper.

Timeline. For n=2n = 2n=2, Hirschberg and Wong and Kannan gave polynomial algorithms in special cases, and Scarf treated the case n=2n=2n=2 completely; the case of general fixed nnn was conjectured by Hirschberg–Wong and Scarf (as recounted on p. 538). Lenstra's paper (received 1981, published 1983) settled the conjecture. Its lattice reduction step uses the algorithm of Lenstra, Lenstra and Lovász (1982). Kannan (1987) later reduced the dependence on nnn to nO(n)n^{O(n)}nO(n).

Setting

Work in Rn\mathbb R^nRn with the Euclidean length ∣⋅∣|\cdot|∣⋅∣, and write B(p,z)={x∈Rn:∣x−p∣≤z}B(p,z)=\{x\in\mathbb R^n : |x-p|\le z\}B(p,z)={x∈Rn:∣x−p∣≤z} for the closed ball with centre ppp and radius z>0z>0z>0.

A lattice is L=∑i=1nZbiL=\sum_{i=1}^n \mathbb Z b_iL=∑i=1n​Zbi​ for a basis b1,…,bnb_1,\dots,b_nb1​,…,bn​ of Rn\mathbb R^nRn. Its determinant d(L)=∣det⁡(b1,…,bn)∣d(L)=|\det(b_1,\dots,b_n)|d(L)=∣det(b1​,…,bn​)∣ (the bib_ibi​ as columns) is the volume of the parallelepiped ∑i[0,1) bi\sum_i[0,1)\,b_i∑i​[0,1)bi​ and does not depend on the basis. Hadamard's inequality (6) says d(L)≤∏i∣bi∣d(L)\le\prod_i|b_i|d(L)≤∏i​∣bi​∣. A basis is reduced in Lenstra's sense if, for a constant c2c_2c2​ depending only on nnn,

∏i=1n∣bi∣≤c2⋅d(L).(7)\prod_{i=1}^n |b_i| \le c_2\cdot d(L). \tag{7}i=1∏n​∣bi​∣≤c2​⋅d(L).(7)

Singling out bnb_nbn​, let H=∑i=1n−1RbiH=\sum_{i=1}^{n-1}\mathbb R b_iH=∑i=1n−1​Rbi​ be the hyperplane spanned by the other vectors, L′=∑i=1n−1ZbiL'=\sum_{i=1}^{n-1}\mathbb Z b_iL′=∑i=1n−1​Zbi​ the lattice inside it, and hhh the distance of bnb_nbn​ to HHH. Then L⊂⋃k∈Z(H+kbn)L\subset\bigcup_{k\in\mathbb Z}(H+kb_n)L⊂⋃k∈Z​(H+kbn​): the lattice lies on a family of parallel hyperplanes H+kbnH + kb_nH+kbn​ at successive distance hhh.

In the paper, K={x:Ax≤b}K=\{x : Ax\le b\}K={x:Ax≤b} is bounded with positive volume, and a linear map τ\tauτ makes it "round":

B(p,r)⊂τK⊂B(p,R)(3),R/r≤c1(4),B(p,r)\subset\tau K\subset B(p,R)\quad(3),\qquad R/r\le c_1\quad(4),B(p,r)⊂τK⊂B(p,R)(3),R/r≤c1​(4),

with c1c_1c1​ depending only on nnn. Then K∩Zn=∅K\cap\mathbb Z^n=\emptysetK∩Zn=∅ if and only if τK∩L=∅\tau K\cap L=\emptysetτK∩L=∅ for L=τZnL=\tau\mathbb Z^nL=τZn, a lattice with basis bi=τ(ei)b_i=\tau(e_i)bi​=τ(ei​).

Formalization targets

Goal (§1, p. 541)

Let b1,…,bnb_1,\dots,b_nb1​,…,bn​ be a basis with (7), numbered so that ∣bn∣=max⁡i∣bi∣|b_n|=\max_i|b_i|∣bn​∣=maxi​∣bi​∣, and let τK\tau KτK satisfy (3) and (4). Then

τK∩L≠∅ort−1<c1c2n,\tau K\cap L\neq\emptyset\quad\text{or}\quad t-1<c_1c_2\sqrt n,τK∩L=∅ort−1<c1​c2​n​,

where ttt is the number of hyperplanes H+kbnH+kb_nH+kbn​, k∈Zk\in\mathbb Zk∈Z, that meet B(p,R)B(p,R)B(p,R) (finiteness included).

Milestones, in the order the paper uses them

  1. LEMMA (8), p. 540: every xxx has y∈Ly\in Ly∈L with ∣x−y∣2≤14(∣b1∣2+⋯+∣bn∣2)|x-y|^2\le\frac14(|b_1|^2+\cdots+|b_n|^2)∣x−y∣2≤41​(∣b1​∣2+⋯+∣bn​∣2).
  2. (10), p. 540: if ∣bn∣|b_n|∣bn​∣ is maximal, ∣x−y∣≤12n ∣bn∣|x-y|\le\frac12\sqrt n\,|b_n|∣x−y∣≤21​n​∣bn​∣.
  3. (11), p. 540: d(L)=h⋅d(L′)d(L)=h\cdot d(L')d(L)=h⋅d(L′).
  4. (12), pp. 540–541: under (7), c2−1∣bn∣≤h≤∣bn∣c_2^{-1}|b_n|\le h\le|b_n|c2−1​∣bn​∣≤h≤∣bn​∣.
  5. p. 541: if B(p,r)⊂τKB(p,r)\subset\tau KB(p,r)⊂τK and τK∩L=∅\tau K\cap L=\emptysetτK∩L=∅, then r<12n ∣bn∣r<\frac12\sqrt n\,|b_n|r<21​n​∣bn​∣.
  6. p. 541: if precisely ttt hyperplanes H+kbnH+kb_nH+kbn​ meet B(p,R)B(p,R)B(p,R), then t−1≤2R/ht-1\le 2R/ht−1≤2R/h.

The goal leaves c1c_1c1​ and c2c_2c2​ free, so it remains valid for any reduction algorithm and any rounding procedure; the companion mission (Integer Programming with a Fixed Number of Variables 2) supplies c1=2n3/2c_1 = 2n^{3/2}c1​=2n3/2, and the LLL algorithm supplies c2=2n(n−1)/4c_2=2^{n(n-1)/4}c2​=2n(n−1)/4.

Significance

The result. The goal is the branching step of Lenstra's algorithm. It says that a lattice-point-free convex body that is well rounded is "flat" in the lattice direction dual to HHH: only boundedly many lattice hyperplanes meet it. The search for a lattice point in dimension nnn therefore reduces to fewer than 1+c1c2n1+c_1c_2\sqrt n1+c1​c2​n​ searches in dimension n−1n-1n−1, which with nnn fixed gives polynomial running time. Sharper forms of the same dichotomy, such as Kannan's bounds, are the basis of later lattice-based integer programming algorithms.

Formalizing it. The results are classical and proved on paper; to our knowledge none of the statements of §1 is machine-checked. Hadamard's inequality (6) is in Mathlib (Orientation.abs_volumeForm_apply_le). The published platform theorem KannanLattice.Core.exists_lattice_point_near_projection (Kannan 1987, Proposition 4.2) proves a stronger form of (8) and (10), with Gram–Schmidt lengths in place of ∣bi∣|b_i|∣bi​∣; for L=ZnL=\mathbb Z^nL=Zn, (10) is the published QFS.exists_lattice_mem_closedBall. What remains is the "base times height" identity (11) relating a full-rank determinant to the distance of a vector from a hyperplane, the spacing bound (12), the count of hyperplanes meeting a ball, and their assembly. This is one of the two geometric ingredients of Lenstra's polynomial algorithm; the algorithm itself and its running time are not formalized.

Difficulty

The argument is short on paper, and each step is "clear" there. The work lies in connecting three descriptions of the same geometry that Mathlib keeps apart: the determinant of a coordinate matrix (d(L)d(L)d(L)), the metric distance of a point to a subspace (hhh), and the (n−1)(n-1)(n−1)-dimensional volume of L′L'L′. Identity (11) needs exactly this bridge, and (12) needs Hadamard's inequality for the lower-dimensional lattice L′L'L′ inside a hyperplane of Rn\mathbb R^nRn, not for a full-rank matrix. Bounding the number of hyperplanes meeting a ball needs the normal direction of HHH and the fact that the hyperplanes are hhh apart along it. A tempting shortcut is to define hhh as the length of the last Gram–Schmidt vector and d(L)d(L)d(L) as the product of Gram–Schmidt lengths: that turns (11) into a tautology and is ruled out below.

Formalization scope

  • Rn\mathbb R^nRn is EuclideanSpace ℝ (Fin n), so norms are Euclidean and balls round; B(p,z)B(p,z)B(p,z) is Metric.closedBall p z with z>0z>0z>0 as a hypothesis where the paper writes a ball.
  • Where bnb_nbn​ is singled out, n=k+1n=k+1n=k+1, the basis is b : Fin (k+1) → ℝⁿ, bnb_nbn​ = b (Fin.last k) and b1,…,bn−1b_1,\dots,b_{n-1}b1​,…,bn−1​ = b ∘ Fin.castSucc. "Basis for LLL" is linear independence of the nnn vectors over R\mathbb RR, and LLL is the published KannanLattice.Core.lattice b, the Z\mathbb ZZ-span.
  • d(L)d(L)d(L) is the paper's definition, ∣det⁡∣|\det|∣det∣ of the column matrix (latDet). hhh is Metric.infDist from bnb_nbn​ to the span HHH. d(L′)d(L')d(L′), for which the paper writes no formula, is the published KannanLattice.Core.latticeDet (the product of Gram–Schmidt lengths of b1,…,bn−1b_1,\dots,b_{n-1}b1​,…,bn−1​, the (n−1)(n-1)(n−1)-volume of their parallelepiped).
  • The hyperplanes H+kbnH+kb_nH+kbn​ meeting B(p,R)B(p,R)B(p,R) are indexed by the set of integers hitIndices b p R; ttt is its cardinality, its finiteness is part of every conclusion (the cardinality of an infinite set is 000 in Lean), and t−1t-1t−1 is computed in R\mathbb RR.
  • The constants "only depending on nnn", c1c_1c1​ and c2c_2c2​, are arbitrary real numbers. The reducedness (7) is a hypothesis, and so is the maximality ∣bi∣≤∣bn∣|b_i|\le|b_n|∣bi​∣≤∣bn​∣, which enters (10), the radius bound and the goal only.
  • τK\tau KτK is an arbitrary set X⊆RnX\subseteq\mathbb R^nX⊆Rn satisfying (3); §1 never uses that KKK is a polyhedron or the map τ\tauτ, and every lattice is τZn\tau\mathbb Z^nτZn for τ(ei)=bi\tau(e_i)=b_iτ(ei​)=bi​. The paper's "for some p∈τKp\in\tau Kp∈τK" follows from (3).
  • The identity (11) is stated without (7), which the paper assumes in the surrounding sentence but does not use.
  • Ruled out: defining d(L)d(L)d(L) or hhh through Gram–Schmidt lengths (which makes (11) bookkeeping), dropping finiteness of the hyperplane count, and replacing the dichotomy "some lattice point lies in τK\tau KτK" by a statement about one computed lattice point.

Welcome contributions: the bridge between Matrix.det of a basis and products of Gram–Schmidt lengths (reusable for any lattice development), Hadamard's inequality for a family of fewer than nnn vectors, and the count of lattice hyperplanes meeting a ball.

Selected references

  • H. W. Lenstra, Jr., Integer Programming with a Fixed Number of Variables, Mathematics of Operations Research 8(4) (1983) 538–548. https://doi.org/10.1287/moor.8.4.538
  • A. K. Lenstra, H. W. Lenstra, Jr., L. Lovász, Factoring polynomials with rational coefficients, Mathematische Annalen 261 (1982) 515–534. https://doi.org/10.1007/BF01457454
  • R. Kannan, Minkowski's convex body theorem and integer programming, Mathematics of Operations Research 12(3) (1987) 415–440. https://doi.org/10.1287/moor.12.3.415
10 thms2 active usersReviewed
Linear OptimizationOperations Research·Captain: mikedeng1

On Linear Characterizations of Combinatorial Optimization Problems III: Zero-One and Integer-Programming-Type Problems Have Small Facial DescriptionsResearch Paper

Motivation

Many discrete optimization problems can be described by linear inequalities, even when their feasible solutions are integer vectors. A linear description lets one study the geometry of the feasible set and connect combinatorial methods with linear programming. The number of inequalities may be large, so a more basic question is whether each inequality needs coefficients of manageable size. Karp and Papadimitriou separated this coefficient-size issue from the number and computational recognizability of the inequalities in their study of facial descriptions Karp and Papadimitriou, MIT/LCS/TM-154, pp. 3–7.

Their Lemma 1 identifies two common families for which small coefficients suffice: zero-one problems and problems given by an integer linear system. This result establishes a geometric premise used by the paper's later complexity discussion. It says nothing by itself about finding the inequalities efficiently or recognizing whether a proposed inequality belongs to a description. Those computational conditions belong to the paper's subsequent theorems Karp and Papadimitriou, pp. 5–8.

Setting

A combinatorial optimization problem here consists of a set LLL of binary strings, a dimension n(z)n(z)n(z) for each z∈Lz\in Lz∈L, and a set S(z)S(z)S(z) of feasible nonnegative integer vectors of that dimension. The geometric object is the rational convex hull CH(S(z))⊆Qn(z)\mathrm{CH}(S(z))\subseteq\mathbb Q^{n(z)}CH(S(z))⊆Qn(z). The paper calls the rationals RRR in its notation. The values of nnn and SSS outside LLL carry no meaning. A feasible set may be empty; its convex hull is then empty as well Karp and Papadimitriou, Definition 1 and footnote, p. 3.

A facial description FFF is a collection of triples ⟨z,f,g⟩\langle z,f,g\rangle⟨z,f,g⟩, where z∈Lz\in Lz∈L, f∈Zn(z)f\in\mathbb Z^{n(z)}f∈Zn(z), and g∈Zg\in\mathbb Zg∈Z. For every rational vector xxx and every z∈Lz\in Lz∈L, membership in CH(S(z))\mathrm{CH}(S(z))CH(S(z)) must be equivalent to satisfying all inequalities f⋅x≤gf\cdot x\le gf⋅x≤g indexed by triples of FFF with first component zzz. A facial description is small when one polynomial in ∣z∣+n(z)|z|+n(z)∣z∣+n(z) bounds the binary size of every coefficient of every triple in the collection. The collection itself need not be small Karp and Papadimitriou, pp. 4–5.

A zero-one problem has only vectors with entries zero or one in each S(z)S(z)S(z). An integer-programming-type problem gives an integral matrix A(z)A(z)A(z) and right-hand side b(z)b(z)b(z) and takes S(z)S(z)S(z) to be all nonnegative integer vectors satisfying A(z)x≤b(z)A(z)x\le b(z)A(z)x≤b(z). The dimensions of these data depend on zzz. Their total binary coefficient size is controlled by the length of the encoded input Karp and Papadimitriou, pp. 5–6.

Formalization targets

Lemma 1

The mission goal is the complete statement for both classes:

∀C,(ZeroOne⁡(C)∨IPType⁡(C))⟹∃F,FacialDescription⁡(C,F)∧Small⁡(C,F).\forall C,\quad \bigl(\operatorname{ZeroOne}(C)\lor\operatorname{IPType}(C)\bigr) \Longrightarrow \exists F,\quad \operatorname{FacialDescription}(C,F)\land\operatorname{Small}(C,F).∀C,(ZeroOne(C)∨IPType(C))⟹∃F,FacialDescription(C,F)∧Small(C,F).

The milestone list records three claims used in the paper's proof. In a zero-one hull, every vertex is a zero-one vector and there are no nonzero rays. For integer hulls, the cited vertex-size result gives a uniform polynomial bound on the coordinates of vertices. Finally, the Cramer's-rule calculation bounds an integral normal fff and intercept ggg by (2nx)n(2nx)^n(2nx)n when the geometric data have coordinates bounded by xxx Karp and Papadimitriou, pp. 5–7.

Significance

Lemma 1 permits descriptions of the whole hull, including lower-dimensional hulls that require equations and unbounded integer hulls that have recession directions. It gives a coefficient bound for every input with one polynomial, rather than a polynomial chosen separately for each instance. This is the premise needed when the paper treats a facial description as a language of short encodable inequalities. The later assertion that an NP-recognizable small facial description places the decision problem in co-NP has a distinct hypothesis and is handled by the first mission of this series Karp and Papadimitriou, Theorem 1, pp. 7–8.

The mathematical result is known; the mission asks for a Lean proof of its precise rational, input-indexed version. The local declarations are proof obligations and do not yet carry verified proofs. Existing platform work includes a rational integer-hull definition reused here, a related open theorem of Cook–Gerards–Schrijver–Tardos that bounds normal coefficients under different hypotheses, and proved real-carrier versions of Meyer's polyhedrality and finite-generation results. The related open theorem does not supply Lemma 1: it has neither this zero-one case nor the bound on the intercept ggg. The real-carrier theorems supply neither the input-size bound nor the paper's rational convention.

Difficulty

The zero-one case has bounded coordinates, but a complete description must also handle an empty hull and a hull contained in a proper affine subspace. Bounding only the facet inequalities of a full-dimensional nonempty polytope does not describe those cases. The integer-programming case adds unbounded directions, and the paper's printed identification of extreme rays with rows of AAA is false. A faithful treatment must bound suitable recession directions and then obtain a uniform coefficient bound for the inequalities describing the hull Karp and Papadimitriou, pp. 5–7.

Formalization scope

The Lean model uses List Bool for inputs, Fin (n z) → ℤ for feasible vectors, and Fin (n z) → ℚ for the convex hull. The published CookSensitivity.ChvatalRank.polyhedron and integerHull definitions provide the general rational integer-programming substrate. This mission adds only the paper's input-indexed problem, facial description, smallness condition, and two named problem classes. The three language-recognition requirements in the paper's Definition 1 are omitted because Lemma 1 never uses them; the resulting theorem applies to every such geometric problem, including every c.o.p. in the paper's narrower definition.

The polynomial bound is represented by (∣z∣+n(z))k+k(|z|+n(z))^k+k(∣z∣+n(z))k+k with a single natural exponent kkk chosen before every input and inequality. For an integer-programming-type input, the sum s(A,b)=∑a⌈log⁡2(1+∣a∣)⌉s(A,b)=\sum_a\lceil\log_2(1+|a|)\rceils(A,b)=∑a​⌈log2​(1+∣a∣)⌉ over all entries of AAA and bbb must be at most ∣z∣|z|∣z∣. This makes the paper's phrase “zzz specifies AAA and bbb” an explicit binary-size condition. The cited vertex milestone bounds coordinates in terms of sss alone, as the page states: a zero column gives a line direction rather than a vertex. Its statement follows the paper's unrestricted Ax≤bAx\le bAx≤b formulation; the IP-type goal additionally imposes x≥0x\ge0x≥0. The paper calls bbb an nnn-vector, but it has one entry per row of AAA and is an mmm-vector here. The Cramer's-rule milestone allows rows of size 2x2x2x, as vertex differences can reach that size.

Empty feasible sets, dimension zero, lower-dimensional hulls, and unbounded IP hulls remain in the goal. A description restricted to full-dimensional facets, or a smallness exponent chosen separately for each input, would not meet it. A complete development needs rational convex hulls, integer-hull geometry, finite descriptions of rational polyhedra, bounds on vertices and recession generators, and finite-dimensional determinant estimates. The rational integer-hull definitions and the coefficient estimates are reusable beyond this mission; contributions that establish the three listed milestones or the remaining geometric steps are in scope.

Selected references

  • Richard M. Karp and Christos H. Papadimitriou, On Linear Characterizations of Combinatorial Optimization Problems, MIT/LCS/TM-154, February 1980, pp. 3–7. Report scan. Later published in SIAM Journal on Computing 11 (1982), 620–632, DOI.
6 thms2 active usersReviewed
Linear OptimizationOperations ResearchOptimization+1·Captain: mikedeng1

Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time 1: The Expected Shadow of a Gaussian-Perturbed Polytope Has Polynomially Many VerticesResearch Paper

Why the shadow of a perturbed polytope matters

The simplex method solves linear programs very fast in practice, yet for most pivot rules there are inputs on which it takes exponentially many steps (Klee and Minty, 1972, for Dantzig's rule; Goldfarb, 1983, for the shadow-vertex rule). Average-case analyses (Borgwardt, 1980s; Smale, 1983) explained good behaviour on random inputs, but random inputs look nothing like real ones. Spielman and Teng introduced smoothed analysis to close this gap: the input is chosen by an adversary and then perturbed by a small Gaussian, and the running time is measured in expectation over the perturbation. They proved that the shadow-vertex simplex method has smoothed complexity polynomial in the number of constraints nnn, the dimension ddd and 1/σ1/\sigma1/σ (Spielman–Teng, J. ACM 2004; this mission follows the preprint arXiv:cs/0111050v7). The work received the Gödel Prize (2008) and the Fulkerson Prize (2009).

Timeline. Borgwardt (1977–1987) bounded the expected number of shadow-vertex pivots for rotationally symmetric random data. Spielman and Teng (2001, STOC; journal 2004) proved the first smoothed bound, with a shadow bound of order nd3/σ6nd^3/\sigma^6nd3/σ6 — the theorem of this mission. Deshpande and Spielman (FOCS 2005) improved the shadow bound, Vershynin (2009) reduced the dependence on nnn to polylogarithmic, and Dadush and Huiberts (STOC 2018) obtained O(d2log⁡n σ−2)O(d^2\sqrt{\log n}\,\sigma^{-2})O(d2logn​σ−2) for small σ\sigmaσ.

Setting

Fix d≥3d\ge3d≥3 and n>dn>dn>d. The data are vectors a1,…,an∈Rda_1,\dots,a_n\in\mathbb R^da1​,…,an​∈Rd, the constraint vectors of the linear program max⁡⟨z∣x⟩\max\langle z|x\ranglemax⟨z∣x⟩ subject to ⟨ai∣x⟩≤1\langle a_i|x\rangle\le1⟨ai​∣x⟩≤1 for all iii. Each aia_iai​ is a Gaussian of standard deviation σ\sigmaσ centered at a point aˉi\bar a_iaˉi​ with ∥aˉi∥≤1\|\bar a_i\|\le1∥aˉi​∥≤1: it has density

μi(a)=(12π σ)de−∥a−aˉi∥2/2σ2,\mu_i(a)=\Big(\tfrac{1}{\sqrt{2\pi}\,\sigma}\Big)^d e^{-\|a-\bar a_i\|^2/2\sigma^2},μi​(a)=(2π​σ1​)de−∥a−aˉi​∥2/2σ2,

and the aia_iai​ are independent (joint density ∏iμi(ai)\prod_i\mu_i(a_i)∏i​μi​(ai​)).

For a direction q∈Rdq\in\mathbb R^dq∈Rd, optSimpq(a1,…,an)\mathrm{optSimp}_q(a_1,\dots,a_n)optSimpq​(a1​,…,an​) is the set of index sets I⊆{1,…,n}I\subseteq\{1,\dots,n\}I⊆{1,…,n} with ∣I∣=d|I|=d∣I∣=d such that (ai)i∈I(a_i)_{i\in I}(ai​)i∈I​ is linearly independent, the simplex △(AI)=ConvHull(ai:i∈I)\triangle(A_I)=\mathrm{ConvHull}(a_i:i\in I)△(AI​)=ConvHull(ai​:i∈I) is a facet of ConvHull(0,a1,…,an)\mathrm{ConvHull}(0,a_1,\dots,a_n)ConvHull(0,a1​,…,an​), and qqq lies in the cone {∑i∈Iαiai:αi≥0}\{\sum_{i\in I}\alpha_ia_i:\alpha_i\ge0\}{∑i∈I​αi​ai​:αi​≥0}. In polar terms, III is the set of tight constraints at the vertex of the feasible polyhedron that maximizes ⟨q∣x⟩\langle q|x\rangle⟨q∣x⟩.

For linearly independent t,zt,zt,z, the shadow Shadowt,z(a1,…,an)\mathrm{Shadow}_{t,z}(a_1,\dots,a_n)Shadowt,z​(a1​,…,an​) is the set of index sets III that belong to optSimpq\mathrm{optSimp}_qoptSimpq​ for some nonzero q∈Span(t,z)q\in\mathrm{Span}(t,z)q∈Span(t,z). Its size is the number of vertices of the projection of the feasible polyhedron onto the plane Span(t,z)\mathrm{Span}(t,z)Span(t,z); the shadow-vertex method walks along this polygon, one pivot per vertex. Finally

D(n,d,σ)=58,888,678 nd3min⁡(σ, 1/(3dln⁡n))6.\mathcal D(n,d,\sigma)=\frac{58{,}888{,}678\,nd^3}{\min\big(\sigma,\,1/(3\sqrt{d\ln n})\big)^6}.D(n,d,σ)=min(σ,1/(3dlnn​))658,888,678nd3​.

Formalization targets

Goal: Theorem 4.0.1 (Shadow Size)

Ea1,…,an[ ∣Shadowt,z(a1,…,an)∣ ]≤D(n,d,σ)\mathbb E_{a_1,\dots,a_n}\big[\,|\mathrm{Shadow}_{t,z}(a_1,\dots,a_n)|\,\big]\le\mathcal D(n,d,\sigma)Ea1​,…,an​​[∣Shadowt,z​(a1​,…,an​)∣]≤D(n,d,σ)

for every d≥3d\ge3d≥3, n>dn>dn>d, every pair of linearly independent t,zt,zt,z, every σ>0\sigma>0σ>0 and all centers of norm at most 111.

Milestones

The milestones follow the paper's proof, leaves first.

  • Probability tools: the chi-square bound (Corollary 2.4.6), the combination lemma (Lemma 2.3.5), almost polynomial densities (Lemma 2.3.7), and comparing Gaussian tails (Lemma 2.4.11).
  • Reduction: the measure of the event P={∥ai∥≤2 ∀i}P=\{\|a_i\|\le2\ \forall i\}P={∥ai​∥≤2 ∀i} (Proposition 4.0.5), and the discretization of the shadow into mmm equally spaced directions (Lemma 4.0.6).
  • Angle bound: the probability, conditioned on PPP, that the ray through a fixed unit vector qqq passes within angle ε\varepsilonε of the boundary of its optimal facet is O(nd3ε/σ6)O(nd^3\varepsilon/\sigma^6)O(nd3ε/σ6) (Lemma 4.0.7, from Lemma 4.0.11).
  • Distance and incidence: in Blaschke coordinates ai=Rωbi+sqa_i=R_\omega b_i+sqai​=Rω​bi​+sq, a deterministic split (Lemma 4.0.12), a distance bound (Lemmas 4.1.1–4.1.3) and an angle-of-incidence bound (Lemmas 4.2.1–4.2.3).

Significance

The result. Theorem 4.0.1 is the geometric heart of the smoothed analysis of the simplex method. Section 4.3 of the paper extends it to arbitrary centers, covariances and right-hand sides, and Section 5 combines these extensions with a two-phase method to show that the simplex method has polynomial smoothed complexity. The same shadow bound underlies later analyses of the simplex method, of perturbed polytopes' diameters, and of condition numbers of random linear programs.

Formalizing it. The theorem has been proved, and improved constants are known, but none of this is machine-checked. A formal proof would verify a long and delicate argument: a change of variables of integral geometry (Blaschke's formula), several conditional-density estimates, and explicit constants in the millions. The mission also produces reusable statements about Gaussian vectors and convex hulls of random points.

Difficulty

The obvious approach is to count, for each candidate facet III, the probability that III appears in the shadow; there are (nd)\binom nd(dn​) candidates, so a union bound is exponential in ddd. The paper avoids this by discretizing the angle of qqq (Lemma 4.0.6) and bounding, for each fixed direction, the probability that the optimal facet changes within a small angular step. That needs a lower bound on the angle between qqq and the boundary of its optimal facet, conditioned on the facet being optimal. The conditioning changes the distribution of a1,…,ada_1,\dots,a_da1​,…,ad​, so the bound cannot come from the Gaussian density alone. The proof changes variables to the facet's normal ω\omegaω, offset sss and in-plane coordinates bib_ibi​ (Corollary 2.5.3), whose Jacobian contributes the factors ⟨ω∣q⟩\langle\omega|q\rangle⟨ω∣q⟩ and Vol(△(b))\mathrm{Vol}(\triangle(b))Vol(△(b)). It then shows that both the distance of the origin to a face of the in-plane simplex and the angle of incidence ⟨ω∣q⟩\langle\omega|q\rangle⟨ω∣q⟩ are unlikely to be small. Measure-theoretic bookkeeping is as hard as the geometry: densities known only up to normalization, conditioning on events of positive measure, and the measure-zero degeneracies the paper sets aside.

Formalization scope

Points live in EuclideanSpace ℝ (Fin d). Constraint vectors are indexed by Fin n (0-based), so the paper's {1,…,d}\{1,\dots,d\}{1,…,d} is {i:i<d}\{i:i<d\}{i:i<d}. The Gaussian of standard deviation σ\sigmaσ centered at ccc is Lebesgue measure with the density above, and the joint law is the product measure. Lemma 4.0.6 also uses Mathlib's multivariateGaussian with a positive definite covariance. Expectations of shadow sizes are lower Lebesgue integrals of [0,∞][0,\infty][0,∞]-valued counts, and their measurability is part of each conclusion. "Density proportional to ν\nuν" and conditional probabilities are stated cross-multiplied, ∫Eν≤bound⋅∫ν\int_{E}\nu\le\text{bound}\cdot\int\nu∫E​ν≤bound⋅∫ν, so no 0/00/00/0 appears.

The shadow is the set of index sets III, and the direction q=0q=0q=0 is excluded. Including it would add every facet of ConvHull(0,a1,…,an)\mathrm{ConvHull}(0,a_1,\dots,a_n)ConvHull(0,a1​,…,an​) to the shadow, since 000 lies in every cone, and make the goal false. ang(q,∅)=∞\mathrm{ang}(q,\emptyset)=\inftyang(q,∅)=∞ is represented exactly in [0,∞][0,\infty][0,∞], never by a real infimum. Where the paper omits a hypothesis it uses, it is added and recorded in the item: the standing assumptions d≥3d\ge3d≥3, n>dn>dn>d and σ≤1/(3dln⁡n)\sigma\le1/(3\sqrt{d\ln n})σ≤1/(3dlnn​) (Lemma 4.2.3 is false without a bound on σ\sigmaσ), unit length of the reference vector qqq, s≥0s\ge0s≥0, and ε>0\varepsilon>0ε>0 for strict inequalities. Lemma 2.3.7 is stated with ≤\le≤ rather than the page's <<<, which fails in an edge case.

Infrastructure a complete development needs: Gaussian tail and chi-square estimates; faces and facets of convex hulls; the Blaschke change of variables and the latitude–longitude change of variables on the sphere (not in Mathlib); surface measure on Sd−1S^{d-1}Sd−1 (Mathlib's Measure.toSphere); and the disintegration of the joint law used in the combination lemma. The Gaussian estimates, the combination lemma and the Blaschke formula are useful beyond this mission. Proofs of any milestone, and of supporting lemmas such as the change-of-variables formulas, are welcome.

Selected references

  • D. A. Spielman, S.-H. Teng, Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time, arXiv:cs/0111050v7, 2003. https://arxiv.org/abs/cs/0111050v7
  • D. A. Spielman, S.-H. Teng, Smoothed analysis of algorithms: Why the simplex algorithm usually takes polynomial time, J. ACM 51(3):385–463, 2004. https://doi.org/10.1145/990308.990310
  • K. H. Borgwardt, The Simplex Method: A Probabilistic Analysis, Springer, 1987.
  • V. Klee, G. J. Minty, How good is the simplex algorithm?, in Inequalities III, Academic Press, 1972, 159–175.
  • A. Deshpande, D. A. Spielman, Improved smoothed analysis of the shadow vertex simplex method, FOCS 2005, 387–396.
  • R. Vershynin, Beyond Hirsch conjecture: walks on random polytopes and smoothed complexity of the simplex method, SIAM J. Comput. 39(2):646–678, 2009. https://doi.org/10.1137/070683386
  • D. Dadush, S. Huiberts, A friendly smoothed analysis of the simplex method, STOC 2018; arXiv:1711.05667. https://arxiv.org/abs/1711.05667
29 thms2 active usersReviewed
Convex OptimizationOperations ResearchOptimization·Captain: Shuze Chen

Discrete Convex Analysis XVII: Fenchel Duality and Linear-Programming IntegralityTextbook

Motivation

Duality is the organizing principle of convex optimization: a minimization problem's optimal value equals a maximization problem's optimal value, and this coincidence, rather than being a lucky accident, follows from a separating-hyperplane argument that applies whenever the two problems' feasible regions are shaped compatibly enough. Werner Fenchel formalized this in the 1950s for pairs of convex and concave functions related by the Legendre-Fenchel transform, and the resulting Fenchel duality theorem specializes, for linear objectives over polyhedral feasible regions, to linear programming duality — the fact, central to the entire theory of combinatorial optimization, that a linear program's optimal value can always be certified from above and below by a pair of primal and dual feasible solutions. Murota's Discrete Convex Analysis (SIAM, 2003) collects this classical machinery, together with the integrality theory that lets it produce combinatorial (integer-valued) certificates rather than merely real ones, as the technical foundation the rest of the book builds its discrete theory on top of.

Setting

For f:Rn→R∪{+∞}f : \mathbb R^n \to \mathbb R \cup \{+\infty\}f:Rn→R∪{+∞}, the epigraph is epi⁡f={(x,Y):Y≥f(x)}\operatorname{epi} f = \{(x,Y) : Y \ge f(x)\}epif={(x,Y):Y≥f(x)}, and fff is convex iff epi⁡f\operatorname{epi} fepif is a convex set; fff is proper if additionally its effective domain dom⁡f={x:f(x)<+∞}\operatorname{dom} f = \{x : f(x) < +\infty\}domf={x:f(x)<+∞} is nonempty, and closed if epi⁡f\operatorname{epi} fepif is topologically closed. A function h:Rn→R∪{−∞}h : \mathbb R^n \to \mathbb R \cup \{-\infty\}h:Rn→R∪{−∞} is concave, proper, closed analogously via its hypograph. The convex conjugate is f∙(p)=sup⁡x{⟨p,x⟩−f(x)}f^\bullet(p) = \sup_x\{\langle p,x\rangle - f(x)\}f∙(p)=supx​{⟨p,x⟩−f(x)}, and the concave conjugate h∘(p)=inf⁡x{⟨p,x⟩−h(x)}h^\circ(p) = \inf_x\{\langle p,x\rangle - h(x)\}h∘(p)=infx​{⟨p,x⟩−h(x)}. The relative interior ri⁡S\operatorname{ri} SriS of a set SSS is the interior of SSS relative to its affine hull. A function is polyhedral if its epigraph (or hypograph) is a finite intersection of half-spaces. Given an m×nm \times nm×n matrix AAA, b∈Rmb \in \mathbb R^mb∈Rm, c∈Rnc \in \mathbb R^nc∈Rn, the primal and dual linear programs are min⁡{c⊤x:Ax=b, x≥0}\min\{c^\top x : Ax=b,\ x\ge0\}min{c⊤x:Ax=b, x≥0} and max⁡{b⊤y:A⊤y≤c}\max\{b^\top y : A^\top y \le c\}max{b⊤y:A⊤y≤c}, with feasible regions PPP, DDD. A matrix is totally unimodular if every square submatrix has determinant 000, 111, or −1-1−1. A discrete set S⊆ZnS \subseteq \mathbb Z^nS⊆Zn is hole free if S=Sˉ∩ZnS = \bar S \cap \mathbb Z^nS=Sˉ∩Zn, where Sˉ\bar SSˉ is the convex hull of SSS's real embedding; the discrete Minkowski sum is S1+S2={x1+x2:x1∈S1,x2∈S2}S_1+S_2 = \{x_1+x_2 : x_1\in S_1, x_2\in S_2\}S1​+S2​={x1​+x2​:x1​∈S1​,x2​∈S2​}.

Formalization targets

Goal (Theorem 3.6, Fenchel duality). For proper convex fff and proper concave hhh satisfying at least one of four alternative conditions — a relative-interior condition on dom⁡f∩dom⁡h\operatorname{dom} f \cap \operatorname{dom} hdomf∩domh, a polyhedrality condition on the same, or the analogous pair of conditions on dom⁡f∙∩dom⁡h∘\operatorname{dom} f^\bullet \cap \operatorname{dom} h^\circdomf∙∩domh∘ together with closedness of fff, hhh —

inf⁡x{f(x)−h(x)}=sup⁡p{h∘(p)−f∙(p)},\inf_x\{f(x)-h(x)\} = \sup_p\{h^\circ(p)-f^\bullet(p)\},xinf​{f(x)−h(x)}=psup​{h∘(p)−f∙(p)},

with the extremum on the appropriate side attained whenever the common value is finite. This is the mission's capstone: the four alternative hypotheses make it the most broadly applicable statement of the four convex-duality results in this mission, each of the other three being either a special case in substance (Theorem 3.5, separation, which 3.6 is proved from) or a literal specialization to linear data (Theorem 3.10, LP duality).

Supporting milestones. Theorem 3.2 (biconjugation: f∙f^\bulletf∙ is always closed proper convex, and g∙∙=gg^{\bullet\bullet}=gg∙∙=g for closed proper convex ggg); Theorem 3.5 (the separation theorem for convex/concave functions, under two of Theorem 3.6's four hypotheses); Theorem 3.9 (the Farkas lemma, equality form); Theorem 3.10 (LP duality: weak duality, strong duality with attainment, and complementary slackness); Theorem 3.13 (total unimodularity of the constraint matrix guarantees an integral optimal solution whenever an optimal solution exists); Proposition 3.14 (an explicit potential function certifying a minimum-weight bipartite perfect matching, via the totally unimodular incidence-matrix LP); and Proposition 3.16 (for a translation-invariant family of hole-free discrete sets, the property that discrete disjointness implies closure disjointness is equivalent to the discrete Minkowski sum matching the integer points of the closures' Minkowski sum).

Significance

Fenchel duality is the single result from which the separation theorem, LP duality, and (via the totally-unimodular incidence matrix of a bipartite graph) the combinatorial duality underlying weighted bipartite matching all descend, in one unbroken chain of specialization; formalizing this chain in one mission exhibits that structure directly, rather than treating each result as an independent fact. Proposition 3.16 plays a different role: it is the chapter's warning that naive discrete analogues of convexity (hole-freeness) do not automatically inherit convexity's good closure properties under Minkowski sums, which is exactly the gap the book's later M-convexity and L-convexity machinery is built to close — this mission's Proposition 3.16 is therefore the motivating negative result for the rest of the book's positive theory, not a loose end. So far as a platform search shows, no existing formalization matches this chunk's specific combination of extended-valued (possibly ±∞\pm\infty±∞) functions, the four-alternative Fenchel duality hypothesis, or the bipartite-matching-via-total-unimodularity argument; the one related platform result (VectorSpaceOpt.fenchel_duality, from Luenberger) is for real-valued functions on general normed spaces under a single relative-interior-and-solidness hypothesis, a different generality from the extended-valued, four-hypothesis statement here.

Difficulty

The naive approach to Theorem 3.6 tries to prove the duality gap is zero directly from the definitions of the two conjugates, which only gives the easy inequality inf⁡≥sup⁡\inf \ge \supinf≥sup (a one-line computation, shown in the book's own proof in three lines); the substantive content is the reverse inequality, and it genuinely fails without a constraint-qualification hypothesis like (a1)-(b2) — Example 3.8 in the book exhibits a convex/concave pair with inf⁡=0≠−1=sup⁡\inf = 0 \ne -1 = \supinf=0=−1=sup when none of the four conditions hold. The book's actual route reduces Theorem 3.6 to the separation theorem (Theorem 3.5) applied to fff shifted down by the (assumed finite) infimum, which produces the separating affine function directly; this is why Theorem 3.5, although logically a special case in spirit, earns its own milestone rather than being subsumed silently.

Formalization scope

All convex and concave functions are represented uniformly as (V → ℝ) → EReal-valued (Fintype V), rather than mixing WithTop ℝ for convex and WithBot ℝ for concave functions, so that Theorem 3.2's biconjugate — whose properness is a conclusion, not an assumption — has a well-defined codomain without extra casts. Convexity is defined via the epigraph being a convex subset of the ordinary real vector space (V→R)×R(V\to\mathbb R)\times\mathbb R(V→R)×R (Mathlib's Convex ℝ), following the book's own equivalent characterization, rather than unfolding the direct inequality definition, which would require a extended-arithmetic scalar-multiplication convention (0\cdot(+\infty)=0) that Mathlib does not provide for EReal. The relative interior is defined directly from the book's own metric-ball-intersected-with-affine-hull description, since Mathlib has no relative-interior primitive at the pinned revision. Polyhedra are finite intersections of explicit half-spaces. A bipartite perfect matching is represented as a bijection between the two vertex sides restricted to the edge set — a faithful, not narrower, representation since every perfect matching between equal-size parts arises this way. The formalization does not trivialize: Theorem 3.6's four hypotheses are carried in full (not reduced to the easiest single case), and no result is stated only for finite-valued (never ±∞\pm\infty±∞) functions, which would discard the entire point of the extended-value convex-analysis framework this chapter sets up for the rest of the book. Infrastructure needed beyond Mathlib's Convex, Matrix, and EReal API: all epigraph/hypograph, conjugate, relative-interior, and polyhedral apparatus is defined fresh in DiscreteConvex.IntegralConvexityB; a contribution proving any of the seven milestones independently, or supplying Mathlib-quality relative-interior lemmas, would be a natural entry point.

Selected references

  • K. Murota, Discrete Convex Analysis, SIAM, 2003, DOI 10.1137/1.9780898718508, Chapter 3.
  • R. T. Rockafellar, Convex Analysis, Princeton University Press, 1970.
  • A. Schrijver, Theory of Linear and Integer Programming, Wiley, 1986.
37 thms2 active usersReviewed
CombinatoricsLinear OptimizationOperations Research·Captain: mikedeng1

On Sub-determinants and the Diameter of Polyhedra: A Polynomial Diameter Bound in the Largest SubdeterminantResearch Paper

Motivation

The combinatorial diameter of a polyhedron is the largest distance, in its vertex-edge graph, between two vertices. It is a lower bound on the number of pivots any edge-following method such as the simplex method needs in the worst case, which is why the polynomial Hirsch conjecture — the diameter of P={x∈Rn:Ax≤b}P = \{x \in \mathbb{R}^n : Ax \le b\}P={x∈Rn:Ax≤b} is bounded by a polynomial in mmm and nnn — is a central open question of linear optimization and discrete geometry. The best general upper bound is quasi-polynomial, m1+log⁡nm^{1+\log n}m1+logn (Kalai–Kleitman 1992); the original Hirsch bound m−nm - nm−n is false for polytopes (Santos 2012).

A different line of work bounds the diameter by the arithmetic of the constraint matrix instead of its size. For an integer matrix AAA let Δ\DeltaΔ be the largest absolute value of a sub-determinant of AAA. Dyer and Frieze (1994) showed that for totally unimodular AAA (Δ=1\Delta = 1Δ=1) the diameter is polynomial, O(m16n3(log⁡mn)3)O(m^{16} n^3 (\log mn)^3)O(m16n3(logmn)3). Bonifas, Di Summa, Eisenbrand, Hähnle and Niemeier (SoCG 2012; Discrete Comput Geom 52, 2014) improved and generalized this to O(Δ2n4log⁡nΔ)O(\Delta^2 n^4 \log n\Delta)O(Δ2n4lognΔ) for all polyhedra and O(Δ2n3.5log⁡nΔ)O(\Delta^2 n^{3.5} \log n\Delta)O(Δ2n3.5lognΔ) for polytopes, bounds that do not depend on the number mmm of inequalities. This mission formalizes the polytope case.

Setting

Let A∈Zm×nA \in \mathbb{Z}^{m\times n}A∈Zm×n with rows a1,…,ama_1,\dots,a_ma1​,…,am​, let b∈Rmb \in \mathbb{R}^mb∈Rm, and let P={x∈Rn:Ax≤b}P = \{x \in \mathbb{R}^n : Ax \le b\}P={x∈Rn:Ax≤b}. A vertex of PPP is an extreme point; for a polyhedron this is a point of PPP at which nnn linearly independent inequalities are tight. Two vertices u≠vu \ne vu=v are adjacent if the segment [u,v][u,v][u,v] is an edge (a one-dimensional face) of PPP. This gives the polyhedral graph GP=(V,E)G_P = (V, E)GP​=(V,E), and the diameter of PPP is at most BBB if every two vertices are joined by a walk of at most BBB edges.

AAA has sub-determinants bounded by Δ\DeltaΔ if every k×kk\times kk×k submatrix, for every k≥1k \ge 1k≥1, has determinant in [−Δ,Δ][-\Delta, \Delta][−Δ,Δ]. In particular every entry is at most Δ\DeltaΔ in absolute value.

For a vertex vvv the normal cone CvC_vCv​ is the set of objectives ccc for which vvv maximizes cTxc^T xcTx over PPP. With BnB_nBn​ the closed unit ball, the volume of a set U⊆VU \subseteq VU⊆V of vertices is

vol(U)=vol(⋃v∈UCv∩Bn),\mathrm{vol}(U) = \mathrm{vol}\Big(\bigcup_{v\in U} C_v \cap B_n\Big),vol(U)=vol(v∈U⋃​Cv​∩Bn​),

and the neighbourhood N(I)\mathcal N(I)N(I) of I⊆VI \subseteq VI⊆V is the set of vertices outside III adjacent to a vertex of III. A spherical cone is S=C∩BnS = C \cap B_nS=C∩Bn​ with CCC closed under non-negative scaling; its dockable surface D(S)D(S)D(S) is the (n−1)(n-1)(n−1)-dimensional measure of the part of its boundary inside the open ball. A cone of revolution of angle 0<θ≤π/20<\theta\le\pi/20<θ≤π/2 is {x∈Bn:vTx≥cos⁡θ ∥v∥ ∥x∥}\{x \in B_n : v^T x \ge \cos\theta\,\|v\|\,\|x\|\}{x∈Bn​:vTx≥cosθ∥v∥∥x∥}. PPP is non-degenerate if every vertex has exactly nnn tight inequalities.

Formalization targets

Goal: Theorem 2 (p. 105)

If A∈Zm×nA \in \mathbb{Z}^{m\times n}A∈Zm×n has all sub-determinants bounded by Δ\DeltaΔ and PPP is bounded, then

diam⁡(P)≤2⌊2π Δ2n5/2ln⁡ ⁣(2n n! nn/2 Δn)⌋+2  =  O(Δ2n3.5log⁡nΔ).\operatorname{diam}(P) \le 2\Big\lfloor \sqrt{2\pi}\,\Delta^2 n^{5/2}\ln\!\big(2^n\, n!\, n^{n/2}\,\Delta^n\big)\Big\rfloor + 2 \;=\; O(\Delta^2 n^{3.5}\log n\Delta).diam(P)≤2⌊2π​Δ2n5/2ln(2nn!nn/2Δn)⌋+2=O(Δ2n3.5lognΔ).

No non-degeneracy, full-dimensionality or rank condition is assumed, and the bound is uniform in mmm and bbb.

Milestones

  1. Lemma 3 (p. 108): for a vertex vvv of a non-degenerate polytope, D(Sv)≤Δ2n3 vol(Sv)D(S_v) \le \Delta^2 n^3\,\mathrm{vol}(S_v)D(Sv​)≤Δ2n3vol(Sv​), where Sv=Cv∩BnS_v = C_v \cap B_nSv​=Cv​∩Bn​.
  2. Lemma 4 (p. 109): among spherical cones of a given volume, a cone of revolution has minimum dockable surface.
  3. Lemma 5 (p. 110): for a cone of revolution, D(S)≥2n/π vol(S)D(S) \ge \sqrt{2n/\pi}\,\mathrm{vol}(S)D(S)≥2n/π​vol(S).
  4. Lemma 6 (p. 111): for every measurable spherical cone with vol(S)≤12vol(Bn)\mathrm{vol}(S) \le \frac12 \mathrm{vol}(B_n)vol(S)≤21​vol(Bn​), D(S)≥2n/π vol(S)D(S) \ge \sqrt{2n/\pi}\,\mathrm{vol}(S)D(S)≥2n/π​vol(S).
  5. Lemma 1 (p. 105): for a non-degenerate polytope and I⊆VI \subseteq VI⊆V with vol(I)≤12vol(Bn)\mathrm{vol}(I) \le \frac12\mathrm{vol}(B_n)vol(I)≤21​vol(Bn​),
vol(N(I))≥2π 1Δ2n2.5 vol(I).\mathrm{vol}(\mathcal N(I)) \ge \sqrt{\tfrac{2}{\pi}}\,\frac{1}{\Delta^2 n^{2.5}}\,\mathrm{vol}(I).vol(N(I))≥π2​​Δ2n2.51​vol(I).
  1. Eq. (1) (p. 105): if IjI_jIj​ is the set of vertices at graph distance at most jjj from a vertex vvv and vol(Ij)≤12vol(Bn)\mathrm{vol}(I_j) \le \frac12\mathrm{vol}(B_n)vol(Ij​)≤21​vol(Bn​), then j≤2π Δ2n2.5ln⁡(2n/vol(I0))j \le \sqrt{2\pi}\,\Delta^2 n^{2.5}\ln(2^n/\mathrm{vol}(I_0))j≤2π​Δ2n2.5ln(2n/vol(I0​)).

Significance

The result. Theorem 2 bounds the diameter of every integral polytope by a polynomial in the dimension and the largest sub-determinant, independently of the number of facets. For totally unimodular matrices, which cover network-flow, bipartite matching and transportation polytopes, it gives O(n3.5log⁡n)O(n^{3.5}\log n)O(n3.5logn), improving the Dyer–Frieze bound by a large polynomial factor. It shows that the obstruction to a polynomial Hirsch bound, if any, must come from matrices with large sub-determinants. The volume-expansion method — measuring breadth-first search by the volume of the normal fan it has covered — was later refined, for instance in the shadow-vertex analysis of Dadush–Hähnle, which improves the dependence on nnn.

Formalizing it. The theorem is proved (2012/2014); no machine-checked proof is known. A formal development needs, on top of Mathlib, the normal fan of a polytope and its relation to the vertex-edge graph, a Hausdorff-measure calculus for cones (surface of a cone in terms of its base), Lévy's isoperimetric inequality on the sphere in a measure-theoretic form, and explicit Gamma-function estimates. Each of these is reusable well beyond this paper.

Difficulty

The combinatorial side is short; the geometry is not. Lemma 4 is the spherical isoperimetric inequality of Lévy, which Mathlib does not have in any form, and which the paper cites rather than proves; the relations between the volume of a spherical cone, the area of its base, its lateral surface and the length of the base's boundary (Eq. (3), "basic integration") are also absent. Lemma 3 depends on the structure of the normal cone of a vertex of a non-degenerate polytope (full-dimensional, simplicial, generated by rows of AAA), none of which is available for Mathlib's extreme points. Lemma 1 depends on the normal fan of a polytope: the normal cones have pairwise disjoint interiors, cover Rn\mathbb{R}^nRn, and share a facet exactly when their vertices are adjacent. The step from non-degenerate to arbitrary polytopes perturbs bbb and needs the diameter not to decrease, a statement about the vertex-edge graph under perturbation. A shortcut through a finite graph abstraction is not available: the constant depends on the geometry of the normal cones, not only on the graph.

Formalization scope

The polyhedron is Hirsch.Hpoly (rowVec A) b, with rowVec A i the iii-th row of A∈A \inA∈ Matrix (Fin m) (Fin n) ℤ as a vector of EuclideanSpace ℝ (Fin n). Vertices are Set.extremePoints ℝ P, adjacency is Hirsch.Adj, "diameter at most BBB" is Hirsch.DiamLE P B, all from the published Hirsch_model. The normal cone is the published FirstOrderOpt.ConvexTheory.normalCone. Volumes are Lebesgue measure with values in [0,∞][0,\infty][0,∞]; the dockable surface uses μHE[n-1], the Hausdorff measure normalized to agree with Lebesgue measure on hyperplanes, applied to frontier S ∩ Metric.ball 0 1. Δ\DeltaΔ is a natural number and the sub-determinant bound ranges over all sizes k≥1k \ge 1k≥1.

Explicit constants. The paper writes O(Δ2n3.5log⁡nΔ)O(\Delta^2 n^{3.5}\log n\Delta)O(Δ2n3.5lognΔ) in Theorem 2; the proof on pp. 105–106 yields 2⌊K⌋+22\lfloor K\rfloor + 22⌊K⌋+2 with K=2π Δ2n5/2ln⁡(2nn! nn/2Δn)K = \sqrt{2\pi}\,\Delta^2 n^{5/2}\ln(2^n n!\, n^{n/2}\Delta^n)K=2π​Δ2n5/2ln(2nn!nn/2Δn), from Eq. (1), the bound vol(I0)≥1/(n! nn/2Δn)\mathrm{vol}(I_0) \ge 1/(n!\,n^{n/2}\Delta^n)vol(I0​)≥1/(n!nn/2Δn) and the fact that the diameter is at most twice the number of breadth-first-search iterations needed to cover more than half of BnB_nBn​. This explicit bound is the goal. The ratios D/volD/\mathrm{vol}D/vol of Lemmas 3, 5, 6 are stated in multiplicative form.

Non-degeneracy is a hypothesis of Lemma 3, Lemma 1 and Eq. (1) only, as in the paper's §1.1, and never of Theorem 2. The neighbourhood N(I)\mathcal N(I)N(I) excludes III; including it would make Lemma 1 trivial, since its constant is below 111. Lemma 4 is stated against every competitor: for every measurable spherical cone SSS and every cone of revolution S∗S^*S∗ of the same volume, D(S∗)≤D(S)D(S^*) \le D(S)D(S∗)≤D(S); it does not assert existence of a cone of a prescribed volume. The goal is Theorem 2 about the polytope and its graph, not an abstract statement about set families with a volume-expansion property; integrality of AAA and the bound on minors of every size are both essential (scaling a real matrix down makes Δ\DeltaΔ arbitrarily small), and the raw Hausdorff measure μH[n-1] would put Lemmas 3 and 6 on incompatible scales.

Contributions are welcome at every level: the normal fan and its adjacency structure, cone surface formulas, the Gamma estimate Γ(x+12)/Γ(x)≥x−14\Gamma(x+\frac12)/\Gamma(x) \ge \sqrt{x-\frac14}Γ(x+21​)/Γ(x)≥x−41​​, and a formal Lévy inequality.

Selected references

  • N. Bonifas, M. Di Summa, F. Eisenbrand, N. Hähnle, M. Niemeier, On Sub-determinants and the Diameter of Polyhedra, Discrete Comput Geom 52 (2014) 102–115. https://doi.org/10.1007/s00454-014-9601-x
  • M. Dyer, A. Frieze, Random walks, totally unimodular matrices, and a randomised dual simplex algorithm, Math. Program. 64 (1994) 1–16. https://doi.org/10.1007/BF01582563
  • G. Kalai, D. J. Kleitman, A quasi-polynomial bound for the diameter of graphs of polyhedra, Bull. Amer. Math. Soc. 26 (1992) 315–316. https://doi.org/10.1090/S0273-0979-1992-00285-9
  • F. Santos, A counterexample to the Hirsch conjecture, Annals of Math. 176 (2012) 383–412. https://doi.org/10.4007/annals.2012.176.1.7
  • T. Figiel, J. Lindenstrauss, V. Milman, The dimension of almost spherical sections of convex bodies, Acta Math. 139 (1977) 53–94 (Lévy's isoperimetric inequality, Theorem 2.1). https://doi.org/10.1007/BF02392234
  • D. Dadush, N. Hähnle, On the shadow simplex method for curved polyhedra, Discrete Comput Geom 56 (2016). https://arxiv.org/abs/1412.6705
11 thms2 active usersReviewed
Combinatorics·Captain: mysticflounder

Superlinear or exact bounds for planar distinct distancesOpen Problem

Superlinear or exact bounds for planar distinct distances

Motivation

This mission asks how restrictions on collinear and cocircular points limit the reuse of distances in the plane. Its central question is Erdős Problem 98: must the minimum number of distances grow faster than the number of points?

Setting

For each positive integer n, let h(n) be the minimum number of distinct positive Euclidean distances determined by an n-point set in the plane with no three collinear points and no four cocircular points. Write D(P) for the number of distinct positive Euclidean distances determined by P.

Target

The mission is to establish a superlinear lower bound, or determine this extremal function exactly. The superlinear target is Erdős Problem 98:

lim⁡n→∞h(n)/n=∞.\lim_{n\to\infty} h(n)/n=\infty.n→∞lim​h(n)/n=∞.

Concretely, for every real A > 0, prove that there is an integer n_A such that every general-position configuration P with |P| = n >= n_A satisfies D(P) > A n. A fixed improvement of the coefficient 1/3, or an additive sublinear improvement above n/3, does not complete this objective.

The alternative completion target is an exact determination of h(n), proved by a universal lower bound and general-position constructions attaining that bound. State the range of n explicitly. An asymptotic estimate or a counterexample to superlinearity alone must be labeled with its actual scope; neither is an exact determination of h(n).

Significance and supporting results

A formally verified supporting result, an n^(1/4) additive bound for distinct distances, is

D(P)≥n/3+cn1/4D(P)\ge n/3+c n^{1/4}D(P)≥n/3+cn1/4

for some absolute c > 0 and all sufficiently large n. A complete Lean proof of this supporting theorem has been accepted on this platform. The n^(1/4), n^(1/5), and n^(1/6) theorem targets and their existing milestones are supporting results, not the mission's terminal goal. Resolving the superlinear target would establish a lower bound above every fixed linear coefficient. Determining h(n) exactly would settle the corresponding extremal problem with matching constructions.

Difficulty and research priorities

The n^(1/4) proof constructs an algebraic curve from a small distinct-distance count, proves pair separation, and constructs intrinsic curve cells to obtain the curve bound D(S) >= c d^(-4/3) |S|^(4/3) for finite samples S on noncircular, irreducible and absolutely irreducible, maximally circular curves of degree d, under |S|^2 >= T d^5 for absolute constants c,T > 0. The final calculation yields an additive n^(1/4) term. Stronger additive bounds count as intermediate progress; they must not be reported as a superlinear lower bound.

  • Develop an argument that excludes D(P) <= A n for every fixed A > 0. A fixed-coefficient distance-energy gap is one sufficient route.
  • Investigate additional structure of the Newton carriers and interactions between their factors, or another geometric or combinatorial route that can control the superlinear target.
  • Investigate constructions and universal lower bounds together when pursuing an exact extremal determination.
  • Preserve and formalize useful intermediate theorems while keeping their statements and remaining premises explicit.

Formalization scope

Configurations are finite subsets of the Euclidean plane. In the indexed formalization, a configuration is an injective map p : Fin n → ℝ²; write P = image(p) before using D(P). Both general-position hypotheses apply to the image. D(P) counts distinct positive distance values, not pairs or ordered multiplicities. The superlinear quantifier ranges over every real A > 0 and every sufficiently large general-position configuration.

The superlinear target and exact-determination target remain open here. Distinguish conjectures, conditional reductions, audited prose proofs, and kernel-checked Lean results. A completed supporting formalization does not by itself complete this mission.

Selected references

  • Erdős Problem 98 — extremal question and bibliography.

  • Proof of a fourth-root additive bound for planar distinct distances (dated 2026-09-13; revision 2026-09-14 UTC; Git revision 9e63c1648017becfc717d7e0e4f2877925ace28d; file SHA-256 8ceaca7bad8c07ce3d56f8bf726a0d273ddc41fd300b777fae1eb6525e087d29), manuscript, revised 2026-09-14.

  • Independent audit of the fourth-root additive bound (dated 2026-09-14 UTC; Git revision 9e63c1648017becfc717d7e0e4f2877925ace28d; file SHA-256 911a06495542d5b4b6cf9f6a41418e4baf4ecd8800b80001b3cc3748f82b0e7b; target SHA-256 8ceaca7bad8c07ce3d56f8bf726a0d273ddc41fd300b777fae1eb6525e087d29), 2026-09-14.

31 thms2 active users
Group Theory·Captain: xuanji

230 space groupsTextbook

Motivation: classify three-dimensional periodic symmetry

A space group describes the rigid motions compatible with a periodic spatial symmetry. The classification concerns possible symmetry types, rather than the size or shape of a particular drawing of a crystal. The classical three-dimensional numbers are 230, 219 when mirror-related types are identified, and 65 for the orientation-preserving subfamily. These are the three numbers recorded in Oliver Knill’s survey, §94, “Crystallography,” p. 41. Keeping their conventions separate matters: changing which coordinate transformations are allowed changes what counts as the same type.

The target is the known classification result selected by LeanEval v1, not an unsolved classification conjecture. Its authoritative specification is the declaration LeanEval.Geometry.SpaceGroupsProblem.space_groups. The accompanying benchmark manifest attributes the classification independently to Fedorov and Schoenflies in 1891. The work requested here is a machine-checked proof of that fixed statement.

Setting: groups, transformations, and orientation

For a natural number ddd, let E(d)=RdE(d)=\mathbb R^dE(d)=Rd with its Euclidean inner product. A Euclidean isometry is an invertible affine distance-preserving transformation of this space. The objects being counted are subgroups GGG of this full motion group. In the benchmark’s definitions, GGG is discrete when, for every point xxx and every real ε>0\varepsilon>0ε>0, the set of its elements satisfying dist⁡(gx,x)≤ε\operatorname{dist}(gx,x)\leq\varepsilondist(gx,x)≤ε is finite.

Such a group is crystallographic when it also contains translations by the members of some linearly independent family of ddd vectors. Translation by vvv means exactly that the transformation sends every xxx to x+vx+vx+v. These conditions specify the underlying groups directly; they do not start with a list of previously classified examples. Although the structure field containing the translation condition is called cocompact, its actual content is the existence of these independent translations, not a separately assumed compact quotient.

An affine equivalence between two groups is an invertible affine map whose conjugation carries the first group’s set of transformations onto the second’s. It need not be an isometry. An orientation-preserving affine equivalence additionally requires the determinant of that affine map’s linear part to be positive. Separately, an individual isometry preserves orientation when its own linear part has positive determinant. The Sohncke subfamily restricts the groups themselves: every element of a group must preserve orientation. This is distinct from restricting the map used to compare two groups, as explicitly distinguished by the source definitions and notes.

Target: one conjunction, with all three exact counts

Write COP(d)C_{\mathrm{OP}}(d)COP​(d) for crystallographicCountOP d, C(d)C(d)C(d) for crystallographicCount d, and COP,only(d)C_{\mathrm{OP,only}}(d)COP,only​(d) for crystallographicCountOPOnly d. They count, respectively, orientation-preserving affine classes of all crystallographic groups, arbitrary affine classes of all crystallographic groups, and orientation-preserving affine classes within the all-elements-orientation-preserving subfamily. The sole goal is

COP(3)=230∧C(3)=219∧COP,only(3)=65.C_{\mathrm{OP}}(3)=230\quad\land\quad C(3)=219\quad\land\quad C_{\mathrm{OP,only}}(3)=65.COP​(3)=230∧C(3)=219∧COP,only​(3)=65.

This is the exact benchmark conjunction, in its original order. None of its three components is optional, and they are not separate theorem targets. There are no milestones or auxiliary theorem items.

Significance: exact cardinalities for the underlying groups

The result gives finite and exact answers for the specified spaces of symmetry types, while retaining the distinction between orientation of a coordinate change and orientation of every symmetry in a group. The difference between 230 and 219 reflects the identification of mirror-related types described in Knill, §94; the 65 count answers a different question, concerning the restricted subfamily. Neither a single count nor a list that silently merges the equivalence conventions establishes the full assertion.

The formalization would add a proof connecting these numerical claims to the actual groups and class subsets specified in Lean. The benchmark source currently supplies the statement with an unproved placeholder. This draft likewise supplies a statement, not a proof or a claim that the benchmark is solved.

Difficulty: a catalog is not a completeness theorem

A finite catalog can have 230 entries without representing every crystallographic group, and different entries can still represent the same affine class. Thus checking the length of a catalog alone does not establish the source’s cardinality assertion. The difficulty is the mathematical connection between concrete descriptions and all groups admitted by the definitions, with precisely the required equivalence relations. The restricted 65-count must also respect the condition on every group element, rather than just a label attached to an example.

Formalization scope: preserve the benchmark model

The Lean representation is EuclideanSpace ℝ (Fin d), with affine isometries and affine equivalences from Mathlib. The three source counting functions use Set.encard in the extended natural numbers N∪{∞}\mathbb N\cup\{\infty\}N∪{∞}. More precisely, each counts the set of subsets obtained as the class of some admissible group; it does not count representatives with multiplicity. Consequently, the displayed finite equalities include finiteness, which is not assumed beforehand.

One reusable definition bundle contains exactly the source’s model spaces, translation and discreteness predicates, crystallographic-group subtype, orientation predicate, two conjugacy relations, and three counting functions. Their declarations are preserved, including definitions for every natural dimension; only the theorem fixes d=3d=3d=3. Definitions for group actions, affine conjugation, and these class subsets can be used independently of this particular count. Contributions must establish the fixed goal with these meanings. Replacing the groups by a hard-coded finite type, defining a count to be its desired answer, or importing an unproved classification into the definition bundle would not establish this target.

Selected references

  • LeanEval contributors; problem submitted by Kim Morrison. LeanEval v1: 230 space groups, statement revision 1, 2026, repository commit 296b7491ec989d21bcf8636a9a69231a1e5d1d25. Exact Lean source; manifest with historical bibliography.
  • Oliver Knill. Some Fundamental Theorems in Mathematics, author-hosted expository survey, July 22, 2018; updated June 25, 2023. §94, “Crystallography,” p. 41. Full text. This provides background for the three counts; the exact formal conventions are those of LeanEval above.
17 thms2 active usersReviewed
Convex OptimizationOperations ResearchProbability·Captain: mikedeng1

Geometrical Bounds for Variance and Recentered Moments: Among Measures with Support Diameter at Most d, Equal Mass on the Vertices of a Regular n-Simplex Maximizes Every Recentered MomentResearch Paper

Why support geometry matters

A probability distribution can spread its mass while remaining inside a set of limited diameter. A bound on the diameter alone does not specify its center, shape, or number of support points. Lim and McCann ask how large a recentered moment can become under this sparse geometric information. Their answer identifies a regular simplex as the extremal configuration: in dimension two it is an equilateral triangle, and in dimension three a regular tetrahedron. This gives an exact higher-dimensional analogue of the familiar one-dimensional observation that a distribution on an interval maximizes variance by placing equal mass at the endpoints. The paper frames the question in connection with systems of particles or organisms that prefer separation from one another, where a diameter constraint limits how far their spatial distribution can spread (Lim and McCann, 2020, pp. 5–7).

The mission follows the first arXiv version of the paper. The theorem is already proved in the source. The formalization task is to express its sharp bound and its equality case precisely in Lean, together with the geometric statements that distinguish the simplex from other diameter-constrained supports.

The objects being compared

Work in Euclidean space En=RnE_n=\mathbb R^nEn​=Rn with its ordinary distance. For a Borel probability measure μ\muμ, its support spt⁡μ\operatorname{spt}\musptμ is the smallest closed set carrying all its mass. The diameter of the support is the supremum of the distances between its points. A bounded diameter places no restriction on where the support sits in space: translating the distribution leaves the diameter unchanged.

For a nonnegative convex function V:En→[0,∞)V:E_n\to[0,\infty)V:En​→[0,∞), the generalized variance allows the center to move:

Var⁡V(μ)=inf⁡z∈En∫EnV(x−z) dμ(x).\operatorname{Var}_V(\mu) =\inf_{z\in E_n}\int_{E_n}V(x-z)\,d\mu(x).VarV​(μ)=z∈En​inf​∫En​​V(x−z)dμ(x).

The center minimizing this integral can depend on VVV and need not be the ordinary barycenter. For V(x)=∥x∥2V(x)=\|x\|^2V(x)=∥x∥2 it is the barycenter and the expression is classical variance. This mission focuses on radial functions V(x)=v(∥x∥)V(x)=v(\|x\|)V(x)=v(∥x∥), where vvv is convex and increasing on nonnegative radii. The paper's definition and the relation to classical variance appear in equation (1.8) and Remark 1.6 (pp. 4–5).

A regular nnn-simplex in EnE_nEn​ has n+1n+1n+1 vertices, and every pair of distinct vertices has the same positive distance ddd. Its uniform vertex measure assigns mass 1/(n+1)1/(n+1)1/(n+1) to each vertex. Write

rn=n2n+2.r_n=\sqrt{\frac{n}{2n+2}}.rn​=2n+2n​​.

This is the distance from the centroid to each vertex when the simplex has side length one. The paper defines regular and unit simplices in Definition 1.8 and identifies this radius in Proposition 3.1 (pp. 6, 13).

Formalization targets

Sharp diameter bound

For a Borel probability measure μ\muμ with diam⁡(spt⁡μ)≤d\operatorname{diam}(\operatorname{spt}\mu)\le ddiam(sptμ)≤d, the target is

Var⁡v∘∥⋅∥(μ)≤v(rnd).\operatorname{Var}_{v\circ\|\cdot\|}(\mu)\le v(r_nd).Varv∘∥⋅∥​(μ)≤v(rn​d).

Here d≥0d\ge0d≥0, so the argument rndr_ndrn​d belongs to the domain on which vvv is specified. The bound fixes the optimal constant in every dimension and for each eligible radial cost (Theorem 1.10, pp. 6–7).

Equality and its geometric milestones

When d>0d>0d>0 and vvv is differentiable at rndr_ndrn​d, equality is characterized by a regular nnn-simplex with side length ddd:

Var⁡v∘∥⋅∥(μ)=v(rnd)⟺μ=1n+1∑i=0nδxi,∥xi−xj∥=d(i≠j).\operatorname{Var}_{v\circ\|\cdot\|}(\mu)=v(r_nd) \quad\Longleftrightarrow\quad \mu=\frac{1}{n+1}\sum_{i=0}^{n}\delta_{x_i}, \qquad \|x_i-x_j\|=d\quad(i\ne j).Varv∘∥⋅∥​(μ)=v(rn​d)⟺μ=n+11​i=0∑n​δxi​​,∥xi​−xj​∥=d(i=j).

The source's proof uses a maximization principle for generalized variance (Theorem 1.7), a constraint placing a minimizing center inside the support's closed convex hull, the two parts of Proposition 3.1 on spherical sets, and the uniqueness of barycentric weights in Remark 3.2. Those statements are the five milestones. Their order follows the progression from the measure problem to the limiting geometry and finally to equal mass (pp. 5, 13, 15).

What the result gives

The bound controls every convex radial recentered moment allowed by the hypotheses using only the diameter and dimension. Its equality clause identifies the entire extremal distribution, rather than merely the value of the maximum. For squared distance, this recovers a sharp variance statement; for powers of the distance, it supplies corresponding recentered moment bounds (Example 1.11, p. 7). The geometric threshold rnr_nrn​ also links the paper's result to Jung's enclosing-ball theorem, which the authors discuss separately (Theorem 1.12 and §3.1, pp. 7, 15–16).

A complete Lean proof would connect measure-theoretic optimization to finite-dimensional convex geometry. It would provide reusable formal statements for generalized variance, support-constrained probability measures, simplex geometry, and the equality mechanism. The source paper proves these results mathematically. In this mission, the Lean statements are drafted and compile with proof placeholders; they do not yet have machine-checked proofs.

Where the difficulty lies

The diameter restriction is imposed on the support, which need not be convex. Replacing it by a restriction on a convex set loses the particular geometry that makes the simplex extremal. The minimizing center is also part of the optimization: evaluating the integral at the barycenter would change the object for a general radial cost. In the equality case, knowing only that the support lies on a sphere is insufficient. The radius must reach exactly rnr_nrn​, the spherical support must be precisely the n+1n+1n+1 simplex vertices, and the mass at those vertices must be equal. Proposition 3.1 and Remark 3.2 isolate these separate assertions (pp. 13–15).

Formalization scope

Lean represents EnE_nEn​ as EuclideanSpace ℝ (Fin n), with the Euclidean norm and Borel measurable space. The support diameter in the goal is the extended metric diameter, which correctly excludes unbounded supports from a finite bound. The paper's conv⁡(K)\operatorname{conv}(K)conv(K) is its smallest closed convex hull, so the statements use closedConvexHull. Probability measures supported on KKK carry the condition μ(Kc)=0\mu(K^c)=0μ(Kc)=0. The generalized variance is the infimum over all centers z∈Enz\in E_nz∈En​; restricting centers to the support hull or replacing the infimum by the integral at the mean would give a different and potentially easier claim.

The function vvv is modeled on all real numbers because Lean functions are total, but every value used by the mathematical statements is at a nonnegative radius. The hypotheses require nonnegativity, convexity, and strict increase on [0,∞)[0,\infty)[0,∞); the latter makes the paper's equality conclusion meaningful. The equality clause additionally assumes differentiability at rndr_ndrn​d. Without it, a kinked strictly increasing convex vvv, with a change of slope at that radius, can have unequal masses attaining the same optimum. The manuscript states the clause without this condition. Theorem 1.7 explicitly assumes KKK nonempty, and the goal explicitly records d≥0d\ge0d≥0 and d>0d>0d>0 where they are needed. Theorem 1.7's support level set follows the sign proved on p. 12, V(x−z∗)=λV(x-z_*)=\lambdaV(x−z∗​)=λ; the printed translated set has the opposite sign.

The definitions and theorem files contain no custom probability or simplex infrastructure beyond these source objects. Contributions proving the optimization principle, the center-in-hull claim, either part of the spherical geometry proposition, or the barycentric uniqueness statement are within scope. Theorem 1.1, the paper's separate convex-envelope result, and the deduction of Jung's theorem are outside this mission.

Selected references

  • Tongseok Lim and Robert J. McCann, Geometrical bounds for the variance and recentered moments, arXiv:2001.11851v1, 2020; published in Mathematics of Operations Research, 2022. arXiv version.
7 thms1 active userReviewed
Harmonic AnalysisMathematical Physics·Captain: vatsj

Thomson problem: N = 9, 10 and 32Open Problem

The problem

Place NNN distinct points x1,…,xNx_1,\dots,x_Nx1​,…,xN​ on the unit sphere in R3\mathbb{R}^3R3 to minimise the Coulomb energy

E(x)=∑i<j1∥xi−xj∥.E(x) = \sum_{i<j} \frac{1}{\lVert x_i - x_j\rVert}.E(x)=i<j∑​∥xi​−xj​∥1​.

Good configurations are easy to find numerically, but proving that one is the global minimum is hard. Rigorous proofs exist for only a few NNN:

  • N=2,3N = 2, 3N=2,3: classical.
  • N=4,6,12N = 4, 6, 12N=4,6,12: the Delsarte–Yudin linear-programming bound is exact (Yudin; Andreev; Cohn–Kumar).
  • N=5N = 5N=5: Schwartz, by a large computer-assisted search.
  • N=8N = 8N=8: Kryvonos–Liehr–Taylor; in Lean by Tooby-Smith–Zughaid.
  • N=7N = 7N=7: a Lean-verified proof by Hung Tran (package, write-up).

This mission

  • N=10N = 10N=10 (capped square antiprism) and N=9N = 9N=9 (tricapped trigonal prism) are stated as "some member of this one-parameter symmetric family is a global minimiser over all configurations", because the optimal parameter has no known closed form. No symmetry is assumed of competing configurations, and uniqueness is not required. N=10N = 10N=10 comes first: like N=7N = 7N=7 it has an antipodal pair of poles, which makes the case split easier.
  • N=32N = 32N=32 (an icosahedron nested with its dual dodecahedron) is the challenge. It is the numerically best-known configuration, but the standard bounds do not close the gap there, so a new idea is needed.

The shared lemmas (the distance formula, two-point positivity behind the Delsarte–Yudin bound, and three-point positivity behind the Bachoc–Vallentin and Cohn–Woo bounds) are stated once for reuse, along with well-posedness of each configuration.

6 thms1 active userReviewed
AlgebraCombinatoricsInformation Theory·Captain: mikedeng1

A Large Family of Maximum Scattered Linear Sets of PG(1, q^n) and Their Associated MRD Codes 2: For Odd q = p^r and t > 4, a Lower Bound on the Number of PΓL-Classes of the Linear Sets L_{h,t}Research Paper

Motivation

A linear set of the projective line PG(1,qn)\mathrm{PG}(1, q^n)PG(1,qn) is the set of points ⟨u⟩Fqn\langle u\rangle_{\mathbb F_{q^n}}⟨u⟩Fqn​​ spanned by the nonzero vectors uuu of an Fq\mathbb F_qFq​-subspace U⊆Fqn2U \subseteq \mathbb F_{q^n}^2U⊆Fqn2​. The extremal ones, the maximum scattered linear sets, are given by scattered polynomials fff over Fqn\mathbb F_{q^n}Fqn​, and each such fff also gives an Fqn\mathbb F_{q^n}Fqn​-linear maximum rank-distance (MRD) code Cf={ax+bf(x):a,b∈Fqn}\mathcal C_f = \{ax + bf(x) : a, b \in \mathbb F_{q^n}\}Cf​={ax+bf(x):a,b∈Fqn​} of minimum distance n−1n-1n−1 (Sheekey, Adv. Math. Commun. 10, 2016). Rank-distance codes are used in network coding, and linear sets appear in the study of semifields, blocking sets and two-intersection sets.

Linear sets and codes are classified up to equivalence: two linear sets are the same object if a collineation of the line, an element of PΓL(2,qn)\mathrm{P\Gamma L}(2, q^n)PΓL(2,qn), maps one onto the other. The question is therefore not only whether new scattered polynomials exist but how many inequivalent linear sets they produce.

  • 2000–2001. The monomials xqsx^{q^s}xqs, gcd⁡(s,n)=1\gcd(s,n) = 1gcd(s,n)=1 (Blokhuis–Lavrauw 2000), all give a single linear set of PG(1,qn)\mathrm{PG}(1,q^n)PG(1,qn) (pseudoregulus type); the binomials δxqs+xqn−s\delta x^{q^s} + x^{q^{n-s}}δxqs+xqn−s (Lunardon–Polverino 2001) give at most φ(n)\varphi(n)φ(n) times a function of qqq classes.
  • 2018. Csajbók, Marino and Polverino (J. Combin. Theory Ser. A 157, 2018) develop the equivalence theory of linear sets of PG(1,qn)\mathrm{PG}(1, q^n)PG(1,qn), including coefficient identities that two qqq-polynomials defining the same linear set must satisfy.
  • 2021. Longobardi and Zanella (J. Algebraic Combin. 2021) give a third family, again with at most φ(n)/2\varphi(n)/2φ(n)/2 classes.
  • 2021. Longobardi, Marino, Trombetti and Zhou (arXiv:2102.08287v3) construct the family ψh,t\psi_{h,t}ψh,t​ (Theorem 3.1, the companion mission) and show (Theorem 5.1) that it contains a number of PΓL\mathrm{P\Gamma L}PΓL-classes growing exponentially in nnn. This mission formalizes that count.

Setting

Let ppp be an odd prime, r≥1r \ge 1r≥1, q=prq = p^rq=pr, t≥1t \ge 1t≥1, n=2tn = 2tn=2t, and let F=FqnF = \mathbb F_{q^n}F=Fqn​ be the finite field of order q2tq^{2t}q2t. Its subfield Fq\mathbb F_qFq​ is {c:cq=c}\{c : c^q = c\}{c:cq=c}.

A qqq-polynomial is a map f(x)=∑i=0n−1αixqif(x) = \sum_{i=0}^{n-1} \alpha_i x^{q^i}f(x)=∑i=0n−1​αi​xqi with αi∈F\alpha_i \in Fαi​∈F (qpoly). It is scattered (IsScatteredPoly) if for all z,y∈F∗z, y \in F^*z,y∈F∗, f(z)/z=f(y)/yf(z)/z = f(y)/yf(z)/z=f(y)/y forces z=cyz = cyz=cy with c∈Fqc \in \mathbb F_qc∈Fq​.

The graph of fff is Uf={(x,f(x)):x∈F}⊆F2U_f = \{(x, f(x)) : x \in F\} \subseteq F^2Uf​={(x,f(x)):x∈F}⊆F2 (graph). A point ⟨v⟩Fqn\langle v \rangle_{\mathbb F_{q^n}}⟨v⟩Fqn​​ of PG(1,qn)\mathrm{PG}(1, q^n)PG(1,qn) is the set of all FFF-multiples of a nonzero v∈F2v \in F^2v∈F2 (point), and the linear set of fff is the point set

Lf={⟨(x,f(x))⟩Fqn:x∈F∗}(linSet).L_f = \{\langle (x, f(x))\rangle_{\mathbb F_{q^n}} : x \in F^*\} \quad (\texttt{linSet}).Lf​={⟨(x,f(x))⟩Fqn​​:x∈F∗}(linSet).

For a field automorphism σ\sigmaσ and an invertible 2×22 \times 22×2 matrix AAA over FFF, the map v↦Avσv \mapsto A v^\sigmav↦Avσ (σ\sigmaσ acting on coordinates, act) is an element of ΓL(2,qn)\Gamma\mathrm L(2, q^n)ΓL(2,qn); with σ=id\sigma = \mathrm{id}σ=id it is in GL(2,qn)\mathrm{GL}(2, q^n)GL(2,qn). Two subsets U,W⊆F2U, W \subseteq F^2U,W⊆F2 are GL- or ΓL-equivalent (GLEquiv, GammaLEquiv) if such a map sends UUU onto WWW. Two point sets are PGL- or PΓL-equivalent (PGLEquiv, PGammaLEquiv) if the induced map on points, P↦{Avσ:v∈P}P \mapsto \{A v^\sigma : v \in P\}P↦{Avσ:v∈P}, sends one onto the other.

The family. For h∈Fh \in Fh∈F,

ψh,t(x)=xq+xqt−1−h1−qt+1xqt+1+h1−q2t−1xq2t−1(psi),\psi_{h,t}(x) = x^q + x^{q^{t-1}} - h^{1-q^{t+1}} x^{q^{t+1}} + h^{1-q^{2t-1}} x^{q^{2t-1}} \quad (\texttt{psi}),ψh,t​(x)=xq+xqt−1−h1−qt+1xqt+1+h1−q2t−1xq2t−1(psi),

H={h∈F:hqt+1=−1}H = \{h \in F : h^{q^t+1} = -1\}H={h∈F:hqt+1=−1} (hSet), Uh=Uψh,tU_h = U_{\psi_{h,t}}Uh​=Uψh,t​​ and Lh,t=Lψh,tL_{h,t} = L_{\psi_{h,t}}Lh,t​=Lψh,t​​. For t≥3t \ge 3t≥3 and h∈H∖Fqth \in H \setminus \mathbb F_{q^t}h∈H∖Fqt​, ψh,t\psi_{h,t}ψh,t​ is scattered and Lh,tL_{h,t}Lh,t​ is a maximum scattered linear set (Theorem 3.1).

Formalization targets

Goal: Theorem 5.1

For t>4t > 4t>4, the number MMM of PΓL(2,qn)\mathrm{P\Gamma L}(2,q^n)PΓL(2,qn)-classes met by the sets Lh,tL_{h,t}Lh,t​, h∈Hh \in Hh∈H, satisfies

M≥{⌊qt+18rt⌋,t≢2(mod4),⌊qt+14rt(q2+1)⌋,t≡2(mod4).M \ge \begin{cases} \left\lfloor \dfrac{q^t+1}{8rt} \right\rfloor, & t \not\equiv 2 \pmod 4,\\[1.5ex] \left\lfloor \dfrac{q^t+1}{4rt(q^2+1)} \right\rfloor, & t \equiv 2 \pmod 4.\end{cases}M≥⎩⎨⎧​⌊8rtqt+1​⌋,⌊4rt(q2+1)qt+1​⌋,​t≡2(mod4),t≡2(mod4).​

It is stated as: some S⊆HS \subseteq HS⊆H of at least that size has pairwise PΓL\mathrm{P\Gamma L}PΓL-inequivalent Lh,tL_{h,t}Lh,t​. The goal fixes the paper's constants; the proof on p. 26 in fact counts classes exactly in terms of class sizes, so a sharper bound is a natural follow-up.

Milestones

  1. Proposition 3.2 (p. 7): h∈H⇒hq2+1≠1h \in H \Rightarrow h^{q^2+1} \ne 1h∈H⇒hq2+1=1 and hqt−2≠−hh^{q^{t-2}} \ne -hhqt−2=−h.
  2. ∣H∣=qt+1|H| = q^t + 1∣H∣=qt+1 (p. 26).
  3. Display (29) (p. 16): for h,k∈Hh, k \in Hh,k∈H, Uh∼GLUkU_h \sim_{\mathrm{GL}} U_kUh​∼GL​Uk​ iff h=±kh = \pm kh=±k (t≢2t \not\equiv 2t≡2) or h=ℓkh = \ell kh=ℓk with ℓq2+1=1\ell^{q^2+1} = 1ℓq2+1=1 (t≡2 mod 4t \equiv 2 \bmod 4t≡2mod4).
  4. Theorem 4.2 (a)–(b) (p. 16), in the UhU_hUh​ form: Uh∼ΓLUkU_h \sim_{\Gamma\mathrm L} U_kUh​∼ΓL​Uk​ iff h=±kρh = \pm k^\rhoh=±kρ (resp. ℓkρ\ell k^\rhoℓkρ) for some ρ∈Aut(Fqn)\rho \in \mathrm{Aut}(\mathbb F_{q^n})ρ∈Aut(Fqn​).
  5. Proof of Corollary 4.3 (p. 18): the GL-class of UhU_hUh​ in the family has ξh=2\xi_h = 2ξh​=2 (resp. q2+1q^2+1q2+1) members, and its ΓL-class at most nrξhnr\xi_hnrξh​.
  6. Lemma 2.2 (p. 5, from Csajbók–Marino–Polverino): Lf=LgL_f = L_gLf​=Lg​ forces α0=β0\alpha_0 = \beta_0α0​=β0​ and two families of coefficient identities.
  7. Display (37) (pp. 20–21): PGL-equivalence of Lf,LgL_f, L_gLf​,Lg​ as an equality of sets of quotients.
  8. The case b=0b = 0b=0 of (37) (p. 21): c=0c = 0c=0.
  9. Lemma 5.2 (p. 21): the case b=0b = 0b=0 for ψh,t,ψk,t\psi_{h,t}, \psi_{k,t}ψh,t​,ψk,t​ gives (h/k)q2+1=1(h/k)^{q^2+1} = 1(h/k)q2+1=1 and Uh∼ΓLUkU_h \sim_{\Gamma\mathrm L} U_kUh​∼ΓL​Uk​.
  10. Lemma 5.3 (p. 22): the case b=1b = 1b=1, with its three conclusions.

Significance

The result. Theorem 5.1 shows that maximum scattered linear sets of PG(1,qn)\mathrm{PG}(1, q^n)PG(1,qn), and hence Fqn\mathbb F_{q^n}Fqn​-linear MRD codes, exist in numbers exponential in nnn up to equivalence, whereas each previously known infinite family yields at most φ(n)\varphi(n)φ(n) times a function of qqq classes (Table 1, p. 3). Remark 5.4 draws the consequence that the family is new. Its intermediate results, the classification of the GL\mathrm{GL}GL- and ΓL\Gamma\mathrm LΓL-equivalences between the UhU_hUh​ (Theorem 4.2) and the reduction of PGL\mathrm{PGL}PGL-equivalence to sets of quotients, are the standard tools for any later classification of this family.

Formalizing it. The theorem has a complete published proof; none of it is machine-checked. The platform has no action of GL(2,qn)\mathrm{GL}(2, q^n)GL(2,qn) or ΓL(2,qn)\Gamma\mathrm L(2, q^n)ΓL(2,qn) on the projective line, no linear sets of PG(1,qn)\mathrm{PG}(1,q^n)PG(1,qn) as point sets, and no coefficient identities for qqq-polynomials defining the same linear set. The mission produces these, then the class counts. Lemma 2.2 is a result of an earlier paper (Csajbók–Marino–Polverino 2018) cited without proof, so formalizing it means formalizing part of that paper too.

Difficulty

Equality of two linear sets is equality of sets of points, which says nothing directly about the polynomials: different qqq-polynomials can define the same linear set (all xqsx^{q^s}xqs with gcd⁡(s,n)=1\gcd(s,n)=1gcd(s,n)=1 do). The obvious idea, that PΓL\mathrm{P\Gamma L}PΓL-equivalent linear sets come from ΓL\Gamma\mathrm LΓL-equivalent subspaces UhU_hUh​, fails in general (p. 3), so the count cannot be read off Theorem 4.2. The step that needs the restriction t>4t > 4t>4 is the comparison of coefficients in the qqq-polynomial identities, which for t=3,4t = 3, 4t=3,4 is open (Remark 5.5). A second difficulty is bookkeeping: the exponents qt±1q^{t\pm1}qt±1, q2t−1q^{2t-1}q2t−1, the negative powers of hhh and the case split on t mod 4t \bmod 4tmod4 must all be tracked through every identity.

Formalization scope

  • Field. FFF is a finite field with CharP F p, p prime, Odd p ("qqq odd"), q = p ^ r with 0 < r, and Fintype.card F = q ^ (2 * t). rrr is the exponent of q=prq = p^rq=pr, as in Theorem 5.1, so ∣Aut(Fqn)∣=rn|\mathrm{Aut}(\mathbb F_{q^n})| = rn∣Aut(Fqn​)∣=rn. Lemma 2.2, (37) and the b=0b = 0b=0 step are stated for general n≥1n \ge 1n≥1 and Fintype.card F = q ^ n, without "qqq odd", as printed.
  • Exponents. xqix^{q^i}xqi is x ^ q ^ i; negative powers of hhh are field quotients (h / h ^ q ^ (t+1) for h1−qt+1h^{1-q^{t+1}}h1−qt+1), exact since h≠0h \ne 0h=0 on HHH. Natural-number subtractions t−1t-1t−1, 2t−12t-12t−1, 3t−23t-23t−2, n−kn-kn−k, qi−1q^i - 1qi−1 are exact in the stated ranges.
  • Geometry. The plane is Fin 2 → F; points are sets of scalar multiples; linear sets are sets of points. Equivalences are stated as images of sets under v↦Avσv \mapsto A v^\sigmav↦Avσ with det⁡A\det AdetA a unit and σ:F≃F\sigma : F \simeq Fσ:F≃F a ring automorphism.
  • Pinned readings.
    1. "Scattered" is the p. 2 quotient definition with Fq={c:cq=c}\mathbb F_q = \{c : c^q = c\}Fq​={c:cq=c}.
    2. Proposition 3.2's hqt−2h^{q^{t-2}}hqt−2 is read as hq3t−2h^{q^{3t-2}}hq3t−2 (equal for t≥2t \ge 2t≥2), and "qqq odd" is added.
    3. Theorem 4.2 (a)–(b), stated for codes Ch,t\mathcal C_{h,t}Ch,t​, is formalized as the ΓL\Gamma\mathrm LΓL-equivalence of Uh,UkU_h, U_kUh​,Uk​, which the proof establishes and which is the code equivalence by Theorem 2.1 (from Sheekey 2016). The automorphism-group sentence is omitted.
    4. Lemma 5.3's hypothesis "(41) holds" is read in the setting of p. 21: b=1b = 1b=1, cˉ=c−da≠0\bar c = c - da \ne 0cˉ=c−da=0, gˉ=∑γiyqi\bar g = \sum \gamma_i y^{q^i}gˉ​=∑γi​yqi inverse to x↦ax+g(x)x \mapsto ax + g(x)x↦ax+g(x) (47), and d+cˉγ0=0d + \bar c\gamma_0 = 0d+cˉγ0​=0 (40).
    5. In (37), ax+bg(x)≠0ax + bg(x) \ne 0ax+bg(x)=0 for x≠0x \ne 0x=0 is part of the condition (stated on p. 21 for b=1b = 1b=1).
    6. "The total number MMM of inequivalent Lh,tL_{h,t}Lh,t​" is the number of classes met by {Lh,t:h∈H}\{L_{h,t} : h \in H\}{Lh,t​:h∈H}, formalized through a pairwise-inequivalent subset of HHH. HHH is not restricted to h∉Fqth \notin \mathbb F_{q^t}h∈/Fqt​, as in the count ∣H∣=qt+1|H| = q^t+1∣H∣=qt+1 of p. 26; "maximum scattered" is descriptive.
  • Trivializing formalizations ruled out. Equivalence is stated on point sets, never on quotient sets {f(x)/x}\{f(x)/x\}{f(x)/x}, where Lean's y/0=0y/0 = 0y/0=0 could merge inequivalent sets; the goal requires pairwise inequivalence, not distinctness of the parameters hhh (note ψ−h,t=ψh,t\psi_{-h,t} = \psi_{h,t}ψ−h,t​=ψh,t​); and the hypotheses are satisfiable, with bound ⌊244/40⌋=6\lfloor 244/40 \rfloor = 6⌊244/40⌋=6 at q=3q = 3q=3, t=5t = 5t=5.
  • Citation basis. All indices and pages refer to arXiv:2102.08287v3 (printed page = PDF page). The Combinatorica version is not used.
  • Welcome contributions. The ΓL(2,qn)\Gamma\mathrm L(2,q^n)ΓL(2,qn) action on F2F^2F2 and on point sets, the cardinality of {h:hqt+1=−1}\{h : h^{q^t+1} = -1\}{h:hqt+1=−1}, and Lemma 2.2 are reusable beyond this mission. Proofs of any milestone, and of Theorem 3.1 in the companion mission, are welcome.

Selected references

  • G. Longobardi, G. Marino, R. Trombetti, Y. Zhou, A large family of maximum scattered linear sets of PG(1, q^n) and their associated MRD codes, arXiv:2102.08287v3, 2021 (Combinatorica 43, 2023). https://arxiv.org/abs/2102.08287v3
  • B. Csajbók, G. Marino, O. Polverino, Classes and equivalence of linear sets in PG(1, q^n), J. Combin. Theory Ser. A 157, 402–426, 2018.
  • J. Sheekey, A new family of linear maximum rank distance codes, Adv. Math. Commun. 10(3), 475–488, 2016.
  • G. Longobardi, C. Zanella, Linear sets and MRD-codes arising from a class of scattered linearized polynomials, J. Algebraic Combin., 2021. https://doi.org/10.1007/s10801-020-01011-9
  • A. Blokhuis, M. Lavrauw, Scattered spaces with respect to a spread in PG(n, q), Geometriae Dedicata 81(1), 231–243, 2000.
  • G. Lunardon, O. Polverino, Blocking sets and derivable partial spreads, J. Algebraic Combin. 14(1), 49–56, 2001.
  • B. Csajbók, C. Zanella, On the equivalence of linear sets, Des. Codes Cryptogr. 81(2), 269–281, 2016.
13 thms1 active userReviewed
PreviousPage 1 of 4Next

Get started

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

About Prove2Me

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

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me