Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Graph Theory

343 missions · 141 completed

Missions

Open202Completed141All343
🏆Completed
Convex OptimizationOptimization·Captain: mikedeng1

Shortest Paths in Graphs of Convex Sets I: The Perspective Mixed-Integer Convex Program (5.5) Has the Same Optimal Value as the Shortest-Path Problem in a Graph of Convex SetsResearch Paper

Motivation

Many planning problems combine a discrete choice with a continuous one. A robot moving through a cluttered room chooses which free regions to traverse and where exactly to cross between them; a vehicle chooses a sequence of operating modes and a trajectory within each mode. Marcucci, Umenberger, Parrilo and Tedrake (arXiv:2101.11565, SIAM J. Optim. 2024, doi:10.1137/22M1523790) abstract this as the shortest-path problem in a graph of convex sets (SPP in GCS): a shortest-path problem in which each vertex is not a point but a convex set, and the length of an edge depends on where in its two sets the path enters and leaves. The paper shows that the problem is NP-hard even in simple cases, and then gives a mixed-integer convex formulation that is exact and whose convex relaxation is strong in practice. That formulation underlies the GCS motion planners developed later by the same group.

Setting

A graph of convex sets consists of a finite directed graph G=(V,E)G = (\mathcal V, \mathcal E)G=(V,E), a source vertex sss, a target vertex t≠st \neq st=s, a nonempty compact convex set Xv⊆Rn\mathcal X_v \subseteq \mathbb R^nXv​⊆Rn for each vertex vvv, and for each edge e=(u,v)e = (u,v)e=(u,v) an edge length ℓe:Rn×Rn→R≥0∪{∞}\ell_e : \mathbb R^n \times \mathbb R^n \to \mathbb R_{\ge 0}\cup\{\infty\}ℓe​:Rn×Rn→R≥0​∪{∞}, which is proper, closed and convex. No symmetry or triangle inequality is assumed. An sss-ttt path is a sequence p=(v0,…,vK)p = (v_0,\dots,v_K)p=(v0​,…,vK​) of distinct vertices with v0=sv_0 = sv0​=s, vK=tv_K = tvK​=t and (vk,vk+1)∈E(v_k, v_{k+1}) \in \mathcal E(vk​,vk+1​)∈E; Ep\mathcal E_pEp​ is the set of its edges. The SPP in GCS (2.1) is

SPP⁡(G):=inf⁡{∑e=(u,v)∈Epℓe(xu,xv):p an s-t path, xv∈Xv for v∈p}.\operatorname{SPP}(G) := \inf\Big\{ \sum_{e=(u,v)\in\mathcal E_p} \ell_e(x_u,x_v) : p \text{ an } s\text{-}t\text{ path},\ x_v \in \mathcal X_v \text{ for } v \in p \Big\}.SPP(G):=inf{e=(u,v)∈Ep​∑​ℓe​(xu​,xv​):p an s-t path, xv​∈Xv​ for v∈p}.

The perspective of a set X\mathcal XX is the closed cone X~:=cl⁡{(x,λ):λ≥0, x∈λX}\tilde{\mathcal X} := \operatorname{cl}\{(x,\lambda) : \lambda \ge 0,\ x \in \lambda\mathcal X\}X~:=cl{(x,λ):λ≥0, x∈λX} (Definition 4.1). The perspective of a function fff is f~(x,λ):=inf⁡{σ:(x,σ,λ)∈epi⁡f~}\tilde f(x,\lambda) := \inf\{\sigma : (x,\sigma,\lambda) \in \widetilde{\operatorname{epi} f}\}f~​(x,λ):=inf{σ:(x,σ,λ)∈epif​}, the function whose epigraph is the perspective of the epigraph of fff (Definition 4.4). For λ>0\lambda > 0λ>0 it equals λf(x/λ)\lambda f(x/\lambda)λf(x/λ).

With Evin\mathcal E^{\mathrm{in}}_vEvin​ and Evout\mathcal E^{\mathrm{out}}_vEvout​ the edges entering and leaving vvv, the MICP (5.5) has a flow ye∈{0,1}y_e \in \{0,1\}ye​∈{0,1} and two vectors ze,ze′∈Rnz_e, z'_e \in \mathbb R^nze​,ze′​∈Rn per edge:

minimize ∑e∈Eℓ~e(ze,ze′,ye)s.t. ∑Esoutye=1, ∑Etinye=1,  ∑Evoutye≤1,  ∑Evin(ze′,ye)=∑Evout(ze,ye)  (v≠s,t),(ze,ye)∈X~u, (ze′,ye)∈X~v, ye∈{0,1}  (e=(u,v)∈E).\begin{aligned} \text{minimize}\ & \textstyle\sum_{e\in\mathcal E}\tilde\ell_e(z_e,z'_e,y_e)\\ \text{s.t.}\ & \textstyle\sum_{\mathcal E^{\mathrm{out}}_s} y_e = 1,\ \sum_{\mathcal E^{\mathrm{in}}_t} y_e = 1,\ \ \sum_{\mathcal E^{\mathrm{out}}_v} y_e \le 1,\ \ \sum_{\mathcal E^{\mathrm{in}}_v}(z'_e,y_e) = \sum_{\mathcal E^{\mathrm{out}}_v}(z_e,y_e)\ \ (v\ne s,t),\\ & (z_e,y_e)\in\tilde{\mathcal X}_u,\ (z'_e,y_e)\in\tilde{\mathcal X}_v,\ y_e\in\{0,1\}\ \ (e=(u,v)\in\mathcal E). \end{aligned}minimize s.t. ​∑e∈E​ℓ~e​(ze​,ze′​,ye​)∑Esout​​ye​=1, ∑Etin​​ye​=1,  ∑Evout​​ye​≤1,  ∑Evin​​(ze′​,ye​)=∑Evout​​(ze​,ye​)  (v=s,t),(ze​,ye​)∈X~u​, (ze′​,ye​)∈X~v​, ye​∈{0,1}  (e=(u,v)∈E).​

Its optimal value is MICP⁡(G)\operatorname{MICP}(G)MICP(G). The paper derives (5.5) from a biconvex program (5.2), the shortest-path LP (5.1) together with vertex positions xv∈Xvx_v \in \mathcal X_vxv​∈Xv​ and the bilinear constraints ze=yexuz_e = y_e x_uze​=ye​xu​, ze′=yexvz'_e = y_e x_vze′​=ye​xv​. Lemma 5.4 turns every valid linear inequality on the flows at a vertex into a valid perspective constraint.

Formalization targets

Goal: Theorem 5.7 (p. 10)

Under the standing assumptions of §2 and the assumption of p. 7 that no edge enters sss and no edge leaves ttt,

MICP⁡(G)=SPP⁡(G).\operatorname{MICP}(G) = \operatorname{SPP}(G).MICP(G)=SPP(G).

Both sides are infima in the extended reals and equal +∞+\infty+∞ when there is no sss-ttt path.

Milestones

  1. Remark 4.2 (p. 5): for bounded closed X\mathcal XX, X~={(x,λ):λ≥0, x∈λX}\tilde{\mathcal X} = \{(x,\lambda): \lambda\ge0,\ x\in\lambda\mathcal X\}X~={(x,λ):λ≥0, x∈λX}.
  2. Remark 4.5 (p. 6): f~(x,λ)=λf(x/λ)\tilde f(x,\lambda) = \lambda f(x/\lambda)f~​(x,λ)=λf(x/λ) for λ>0\lambda>0λ>0, f~=∞\tilde f = \inftyf~​=∞ for λ<0\lambda<0λ<0, f~(0,0)=0\tilde f(0,0)=0f~​(0,0)=0 for proper closed convex fff.
  3. Lemma 5.4 (p. 9): a valid inequality ∑e∈Evceye+d≥0\sum_{e\in\mathcal E_v} c_e y_e + d \ge 0∑e∈Ev​​ce​ye​+d≥0 for (5.2) lifts to the valid constraint (∑Evinceze′+∑Evoutceze+dxv, ∑Evceye+d)∈X~v\big(\sum_{\mathcal E^{\mathrm{in}}_v} c_e z'_e + \sum_{\mathcal E^{\mathrm{out}}_v} c_e z_e + d x_v,\ \sum_{\mathcal E_v} c_e y_e + d\big)\in\tilde{\mathcal X}_v(∑Evin​​ce​ze′​+∑Evout​​ce​ze​+dxv​, ∑Ev​​ce​ye​+d)∈X~v​.
  4. Remark 5.5 (p. 9): the equality version of Lemma 5.4.
  5. Path and cycles (proof of Theorem 5.7): binary flows satisfying (5.5b)–(5.5d) form the vertex-disjoint union of an sss-ttt path and cycles.
  6. Simplification of (5.5e) (proof of Theorem 5.7): ye=0y_e = 0ye​=0 forces ze=ze′=0z_e = z'_e = 0ze​=ze′​=0 at cost 000; ye=1y_e = 1ye​=1 forces ze∈Xuz_e\in\mathcal X_uze​∈Xu​, ze′∈Xvz'_e\in\mathcal X_vze′​∈Xv​ at cost ℓe(ze,ze′)\ell_e(z_e,z'_e)ℓe​(ze​,ze′​).
  7. Flow along the path (proof of Theorem 5.7): for consecutive path edges e,fe, fe,f, ze′=zfz'_e = z_fze′​=zf​.

A companion item states the recovery part of Theorem 5.7 (second and third sentences): an optimal MICP solution yields an optimal path and optimal positions xu=zex_u = z_exu​=ze​, xv=ze′x_v = z'_exv​=ze′​.

Significance

Theorem 5.7 is what makes (5.5) usable: it says that the mixed-integer program, whose only nonconvexity is ye∈{0,1}y_e\in\{0,1\}ye​∈{0,1}, solves the SPP in GCS exactly. The SPP in GCS is NP-hard by Theorem 3.1, so exactness is the most one can ask of a polynomial-size formulation, and dropping integrality gives a convex program whose value is a lower bound for branch and bound. The perspective construction of Lemma 5.4 is a general recipe, closely related to the reformulation-linearization technique, for convexifying products of flows and positions.

On the formal side, the paper's proof is half a page and leaves several points to the reader: why the integer flows decompose into a path and cycles, why the closure in the perspective cone disappears for compact sets, and why the perspective function vanishes at the origin. A machine-checked proof makes these steps explicit. The result has not been formalized before, and the mission also produces reusable statements about perspectives of sets and functions in extended-real arithmetic and about the decomposition of unit sss-ttt flows with degree bounds.

Difficulty

The inequality MICP⁡(G)≤SPP⁡(G)\operatorname{MICP}(G) \le \operatorname{SPP}(G)MICP(G)≤SPP(G) is the routine half: a path and its positions give a feasible point of (5.5) of the same cost. The other half requires reading a path out of an arbitrary feasible point. Integer flows satisfying (5.5b)–(5.5d) need not be a path: they may also contain disjoint cycles, and the degree constraint (5.5c) is what prevents a cycle from sharing a vertex with the path (Example 5.10 shows the value drops without it). The cost of a cycle must be shown to be nonnegative, and the perspective constraints must be shown to collapse to ze=ze′=0z_e = z'_e = 0ze​=ze′​=0 on zero-flow edges. That step uses compactness of Xv\mathcal X_vXv​. For unbounded sets the zero slice of X~v\tilde{\mathcal X}_vX~v​ is the recession cone, and the argument fails. The cost at zero flow is ℓ~e(0,0,0)\tilde\ell_e(0,0,0)ℓ~e​(0,0,0), defined by a closure, and its value 000 uses properness and convexity of ℓe\ell_eℓe​, not just the formula λℓe(⋅/λ)\lambda\ell_e(\cdot/\lambda)λℓe​(⋅/λ).

Formalization scope

  • Rn\mathbb R^nRn is Fin n → ℝ. Edge lengths and all optimal values are EReal; optimal values are infima (⨅/sInf), equal to +∞+\infty+∞ on an empty feasible set. Attainment is not assumed.
  • The graph is a structure with a finite edge set E : Finset (V × V) over a finite vertex type, a source, a target, sets X v and lengths ℓ e. The standing assumptions of §2 (s≠ts\ne ts=t; Xv\mathcal X_vXv​ nonempty, compact, convex; ℓe≥0\ell_e \ge 0ℓe​≥0 proper closed convex) are the predicate IsGCS. The assumption ∣Esin∣=∣Etout∣=0|\mathcal E^{\mathrm{in}}_s| = |\mathcal E^{\mathrm{out}}_t| = 0∣Esin​∣=∣Etout​∣=0 of p. 7 is the separate predicate NoInSNoOutT, a hypothesis of the goal, of the path-and-cycles milestone and of the recovery item.
  • Paths are duplicate-free vertex lists from sss to ttt along edges; Ep\mathcal E_pEp​ is the list of consecutive pairs.
  • The perspective of a function is Definition 4.4's infimum over the closed perspective cone of the epigraph. The formula λf(x/λ)\lambda f(x/\lambda)λf(x/λ) is not used as a definition: in Lean it would evaluate to 000 at λ=0\lambda = 0λ=0 for every argument and erase the content of (5.5e). The integrality constraint ye∈{0,1}y_e \in\{0,1\}ye​∈{0,1} is part of the feasible set; without it the goal is false (§9.4 reports relaxation gaps).
  • Disclosed additions: Lemma 5.4 and Remark 5.5 assume no self-loop at vvv, which the page's partition of Ev\mathcal E_vEv​ into incoming and outgoing edges presupposes. The recovery item takes the sss-ttt path contained in {e:ye=1}\{e : y_e = 1\}{e:ye​=1}, since the page's claim that this set is itself a path fails when an optimal solution carries a zero-length cycle.
  • Needed infrastructure: closures of cones over compact sets, extended-real infima, and the combinatorics of unit flows on finite digraphs. Proofs of any milestone, and reusable lemmas on perspectives, are welcome.

Selected references

  • T. Marcucci, J. Umenberger, P. A. Parrilo, R. Tedrake, Shortest Paths in Graphs of Convex Sets, SIAM Journal on Optimization 34(1), 2024; arXiv:2101.11565v5. https://arxiv.org/abs/2101.11565
  • R. T. Rockafellar, Convex Analysis, Princeton University Press, 1970 (perspectives and recession cones, §8). https://doi.org/10.1515/9781400873173
  • R. K. Ahuja, T. L. Magnanti, J. B. Orlin, Network Flows: Theory, Algorithms, and Applications, Prentice Hall, 1993 (flow decomposition, shortest-path LP).
  • T. Marcucci, M. Petersen, D. von Wrangel, R. Tedrake, Motion Planning around Obstacles with Convex Optimization, Science Robotics 8(84), 2023. https://doi.org/10.1126/scirobotics.adf7843
11 thms2 active usersReviewed
🏆Completed
CombinatoricsLinear algebra·Captain: mikedeng1

Three Conjectures in Extremal Spectral Graph Theory 2: For Large n, K₂ + Pₙ₋₂ Is the Unique Planar Graph on n Vertices of Maximum Spectral Radius (Boots–Royle / Cao–Vince Conjecture)Research Paper

Motivation

The spectral radius of a graph, the largest eigenvalue of its adjacency matrix, measures how densely and how unevenly the edges of the graph are distributed. Asking which graph in a class maximizes it is a spectral analogue of the extremal question "which graph in the class has the most edges". For planar graphs the question has a long history in spectral graph theory, going back at least to Schwenk and Wilson. It also arose in geography, where the spectral radius serves as a measure of the connectivity of planar networks. To compare real networks with a theoretical upper bound, geographers wanted the planar graph of maximum spectral radius (arXiv:1606.01916v2, p. 3).

Timeline, as recorded by Tait and Tobin (pp. 1–3):

  • 1978, Schwenk and Wilson: early study of the eigenvalues of planar graphs.
  • 1991, Boots and Royle, and independently 1993, Cao and Vince: conjecture that the planar graph on n≥9n \ge 9n≥9 vertices of maximum spectral radius is P2+Pn−2P_2 + P_{n-2}P2​+Pn−2​ (Conjecture 1 of the paper). The conjecture is false for n∈{7,8}n \in \{7, 8\}n∈{7,8}, so some lower bound on nnn is necessary (p. 2).
  • 1988–2000: successive improvements of the upper bound on the spectral radius of planar graphs, by Hong (1988, 1995, 1998), Cao and Vince (1993), Guiduli (1996) and Ellingham and Zha (2000). These bounds do not identify the extremal graph.
  • 2010, Dvořák and Mohar: an upper bound on the spectral radius of planar graphs of given maximum degree.
  • Guiduli and Hayes were reported by Ellingham and Zha to have proved the conjecture for large nnn. The preprint never appeared (p. 3).
  • 2017, Tait and Tobin: proof of the conjecture for all sufficiently large nnn (Theorem 15).

Setting

All graphs are finite and simple. A graph on nnn vertices has vertex set {0,1,…,n−1}\{0, 1, \dots, n-1\}{0,1,…,n−1}. Its adjacency matrix A(G)A(G)A(G) is the n×nn \times nn×n real matrix with entry 111 in position (u,w)(u, w)(u,w) when u∼wu \sim wu∼w and 000 otherwise. Its eigenvalues are real, and the largest is the spectral radius λ1(G)\lambda_1(G)λ1​(G).

A graph is planar if it can be drawn in the plane with vertices as distinct points and edges as simple arcs joining their ends, such that an arc meets no vertex other than its ends and two arcs meet only at common ends.

For graphs GGG and HHH on disjoint vertex sets, the join G+HG + HG+H is the graph on the union of the two vertex sets that keeps every edge of GGG and of HHH and adds every edge between a vertex of GGG and a vertex of HHH. Write K2K_2K2​ for a single edge, PmP_mPm​ for the path on mmm vertices, and K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​ for their join: two adjacent vertices, both joined to every vertex of a path on the remaining n−2n - 2n−2 vertices. The paper also writes it P2+Pn−2P_2 + P_{n-2}P2​+Pn−2​.

A planar graph GGG on nnn vertices is a planar spectral maximizer if λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every planar graph G′G'G′ on the same nnn vertices.

The milestones refer to a leading eigenvector v\mathbf vv of a maximizer, scaled so that every entry is positive and the largest entry is 111, and to a vertex xxx with vx=1\mathbf v_x = 1vx​=1. For a parameter ε>0\varepsilon > 0ε>0 they use the sets L={z:vz>ε}L = \{z : \mathbf v_z > \varepsilon\}L={z:vz​>ε} and S=V(G)∖LS = V(G) \setminus LS=V(G)∖L, the edge counts e(X)e(X)e(X) (edges inside XXX) and e(X,Y)e(X, Y)e(X,Y) (edges between XXX and YYY), and a second vertex www with vw>1−24ε\mathbf v_w > 1 - 24\varepsilonvw​>1−24ε that is adjacent to all but 94εn94\varepsilon n94εn vertices of SSS.

Formalization targets

Goal: Theorem 15

∃N0 ∀n≥N0 ∀G on n vertices:G is a planar spectral maximizer  ⟺  G≅K2+Pn−2.\exists N_0\ \forall n \ge N_0\ \forall G \text{ on } n \text{ vertices}: \quad G \text{ is a planar spectral maximizer} \iff G \cong K_2 + P_{n-2}.∃N0​ ∀n≥N0​ ∀G on n vertices:G is a planar spectral maximizer⟺G≅K2​+Pn−2​.

The threshold N0N_0N0​ is not fixed. The paper proves the statement only for nnn beyond an unspecified bound, and the conjecture's bound n≥9n \ge 9n≥9 remains open. The equivalence encodes "the unique maximizer is K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​": every maximizer is isomorphic to it, and it is itself a maximizer.

Milestones, in the order of the proof

  1. If HHH is a proper spanning subgraph of a connected graph GGG, then λ1(H)<λ1(G)\lambda_1(H) < \lambda_1(G)λ1​(H)<λ1​(G) (§1.2, already proved on the platform).
  2. A planar graph on n≥3n \ge 3n≥3 vertices has no K3,3K_{3,3}K3,3​ subgraph and at most 3n−63n - 63n−6 edges, and each of its bipartite subgraphs has at most 2n−42n - 42n−4 edges (p. 8).
  3. Lemma 8: 2n−4<λ1<6n\sqrt{2n-4} < \lambda_1 < \sqrt{6n}2n−4​<λ1​<6n​.
  4. The counts ∣L∣≤32n−4/ε|L| \le 3\sqrt{2n-4}/\varepsilon∣L∣≤32n−4​/ε, e(S,L)≤2n−4e(S, L) \le 2n - 4e(S,L)≤2n−4, e(S)≤3n−6e(S) \le 3n - 6e(S)≤3n−6, e(L)≤92n−4/εe(L) \le 9\sqrt{2n-4}/\varepsilone(L)≤92n−4​/ε (p. 8).
  5. Lemma 9: ∑z∈Lvz≤ε2n−4+18/ε\sum_{z \in L} \mathbf v_z \le \varepsilon\sqrt{2n-4} + 18/\varepsilon∑z∈L​vz​≤ε2n−4​+18/ε and ∑z∈Svz≤(1+3ε)2n−4\sum_{z \in S} \mathbf v_z \le (1 + 3\varepsilon)\sqrt{2n-4}∑z∈S​vz​≤(1+3ε)2n−4​.
  6. Eq. (7): ∑y∼u, y∈Svy≥(vu−ε)2n−4−18/ε\sum_{y \sim u,\, y \in S} \mathbf v_y \ge (\mathbf v_u - \varepsilon)\sqrt{2n-4} - 18/\varepsilon∑y∼u,y∈S​vy​≥(vu​−ε)2n−4​−18/ε for u∈Lu \in Lu∈L.
  7. Lemma 10: vz>1/6n\mathbf v_z > 1/\sqrt{6n}vz​>1/6n​ for every vertex zzz.
  8. ∣{y∈S:y≁x}∣≤14εn|\{y \in S : y \not\sim x\}| \le 14\varepsilon n∣{y∈S:y∼x}∣≤14εn (p. 10).
  9. Lemma 11: a vertex w∈Lw \in Lw∈L, w≠xw \ne xw=x, with vw>1−24ε\mathbf v_w > 1 - 24\varepsilonvw​>1−24ε and ∣{y∈S:y≁w}∣≤94εn|\{y \in S : y \not\sim w\}| \le 94\varepsilon n∣{y∈S:y∼w}∣≤94εn.
  10. Lemma 12: vu<1/10\mathbf v_u < 1/10vu​<1/10 for every u∉{x,w}u \notin \{x, w\}u∈/{x,w}.
  11. Lemma 13: x∼wx \sim wx∼w.
  12. Lemma 14: every vertex other than xxx and www is adjacent to both.
  13. A planar graph with two distinct vertices adjacent to all other vertices is a spanning subgraph of K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​ (proof of Theorem 15, p. 13).

Significance

The result. Theorem 15 settles the Boots–Royle / Cao–Vince conjecture for all large nnn. It identifies the planar graph of maximum spectral radius exactly, so the theoretical upper bound that the geographic application asks for is λ1(K2+Pn−2)\lambda_1(K_2 + P_{n-2})λ1​(K2​+Pn−2​). Along the way, Lemmas 8–14 describe the eigenvector of the extremal graph: two vertices carry almost all of the eigenvector weight, and every other entry is below 1/101/101/10. The same method proves the outerplanar analogue (Theorem 7, the companion mission of this series).

Formalizing it. The theorem is proved on paper and has no machine-checked proof. Formalizing it needs several pieces that are absent from Mathlib or the platform and are reusable: the edge bound 3n−63n - 63n−6 and the bipartite bound 2n−42n - 42n−4 for topologically planar graphs, the non-planarity of K3,3K_{3,3}K3,3​ and K5K_5K5​ and of their subdivisions, Rayleigh-quotient comparison of spectral radii, and the face structure of maximal planar graphs. Each milestone is a self-contained target.

Difficulty

Counting edges does not separate the candidates. A maximizer can be taken to be a maximal planar graph, so it has exactly 3n−63n - 63n−6 edges. Every maximal planar graph on nnn vertices has the same number of edges, and the standard lower bound λ1≥2∣E∣/n\lambda_1 \ge 2|E|/nλ1​≥2∣E∣/n is the same for all of them. The extremal graph is pinned down only through its eigenvector. One has to show that the weight concentrates on two vertices, which are then forced to be adjacent to each other (Lemma 13) and to every other vertex (Lemma 14). Lemma 13 is the step specific to planarity. It uses the triangular faces of a maximal planar graph to find an edge switch that keeps the graph planar. In the formal setting, planarity is topological, so every planarity fact of the argument (the edge counts, the K3,3K_{3,3}K3,3​ and K5K_5K5​ obstructions, the faces of a triangulation) has to be derived from the drawing.

Formalization scope

  • Graphs are SimpleGraph (Fin n). "GGG is K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​" is a graph isomorphism G ≃g bookPath n, where bookPath n lives on Fin 2 ⊕ Fin (n - 2). It is not an equality on a fixed labelling.
  • λ1\lambda_1λ1​ is the published WangKangXue.SpectralTuran.specRad, the largest eigenvalue of the real adjacency matrix.
  • Planarity is the published topological RobertsonSeymour1986.GM5.IsPlanar: a drawing by simple arcs. It is not replaced by K5K_5K5​/K3,3K_{3,3}K3,3​-minor-freeness or by an edge count, either of which would change the class the maximum ranges over.
  • Extremality is IsSpecMax IsPlanar G: GGG is planar and λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every planar G′G'G′ on Fin n. Every milestone from Lemma 8 to Lemma 14 assumes it. These lemmas are false for an arbitrary planar graph, so a statement that drops the hypothesis, or replaces it by a property implying the conclusion, is not acceptable.
  • The eigenvector is universally quantified: every v\mathbf vv with Av=λ1vA\mathbf v = \lambda_1 \mathbf vAv=λ1​v, all entries positive, all entries at most 111, and every xxx with vx=1\mathbf v_x = 1vx​=1. Such a vector exists because the maximizer is connected (Perron–Frobenius). A statement that only asserts the existence of some eigenvector is weaker and is ruled out.
  • Quantifiers. "For nnn sufficiently large" is ∃N ∀n≥N\exists N\ \forall n \ge N∃N ∀n≥N, and NNN may depend on ε\varepsilonε. Milestones 4–8 hold for every ε>0\varepsilon > 0ε>0. Lemmas 11–14 hold for 0<ε<ε00 < \varepsilon < \varepsilon_00<ε<ε0​ with an unspecified ε0\varepsilon_0ε0​, because the paper chooses ε\varepsilonε small. The quantifier on ε\varepsilonε is always outermost, so ε\varepsilonε never depends on nnn or GGG. Lemmas 12–14 hold for every vertex www satisfying the conclusion of Lemma 11 (the predicate IsSecondHub), with the constants 242424 and 949494.
  • Added hypotheses. n≥3n \ge 3n≥3 in the planar facts and in Lemma 8, where the paper's standing "nnn large enough" is needed only to that extent. The page names the vertex of Lemma 12 vvv, the same letter as the eigenvector; it is renamed uuu. Lemma 14's set AAA is stated in unfolded form.
  • Infrastructure that a complete development needs: Euler's formula for topologically planar graphs, the non-planarity of K3,3K_{3,3}K3,3​ and K5K_5K5​, Perron–Frobenius for connected graphs, and the Rayleigh quotient. Proofs of the planar facts (milestone 2) are reusable well beyond this mission. Contributions toward any of them are welcome.

Selected references

  • M. Tait and J. Tobin, Three conjectures in extremal spectral graph theory, J. Combin. Theory Ser. B (2017). Preprint arXiv:1606.01916v2. https://arxiv.org/abs/1606.01916, DOI https://doi.org/10.1016/j.jctb.2017.04.006
  • B. N. Boots and G. F. Royle, A conjecture on the maximum value of the principal eigenvalue of a planar graph, Geographical Analysis 23(3) (1991), 276–282.
  • D. Cao and A. Vince, The spectral radius of a planar graph, Linear Algebra Appl. 187 (1993), 251–257.
  • A. J. Schwenk and R. J. Wilson, On the eigenvalues of a graph, in: Selected Topics in Graph Theory, Academic Press, 1978, 307–336.
  • Hong Yuan, A bound on the spectral radius of graphs, Linear Algebra Appl. 108 (1988), 135–139.
  • M. N. Ellingham and X. Zha, The spectral radius of graphs on surfaces, J. Combin. Theory Ser. B 78(1) (2000), 45–56.
  • Z. Dvořák and B. Mohar, Spectral radius of finite and infinite planar graphs and of graphs of bounded genus, J. Combin. Theory Ser. B 100(6) (2010), 729–739.

Bibliographic data as given in the bibliography of Tait and Tobin (pp. 21–23).

21 thms2 active usersReviewed
🏆Completed
CombinatoricsLinear algebra·Captain: mikedeng1

Three Conjectures in Extremal Spectral Graph Theory 1: For Large n, K₁ + Pₙ₋₁ Is the Unique Outerplanar Graph on n Vertices of Maximum Spectral Radius (Cvetković–Rowlinson Conjecture)Research Paper

Motivation

The spectral radius λ1(G)\lambda_1(G)λ1​(G) of a graph GGG, the largest eigenvalue of its adjacency matrix, is a standard measure of how densely and how centrally a graph is connected. A classical theme of spectral graph theory asks, for a natural family of graphs, which member of the family on nnn vertices maximizes λ1\lambda_1λ1​. For planar graphs this question was raised by geographers comparing road and transport networks, and for outerplanar graphs it was posed by Cvetković and Rowlinson in their 1990 survey of the largest eigenvalue of a graph, where they relate it to the study of subfamilies of Hamiltonian graphs.

Conjecture (Cvetković–Rowlinson, 1990). The outerplanar graph on nnn vertices of maximum spectral radius is K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​.

Timeline:

  • 1990. Cvetković and Rowlinson state the conjecture (Linear and Multilinear Algebra 28).
  • 1990. Rowlinson makes partial progress, on the index of certain outerplanar graphs (Ars Combinatoria 29).
  • 1993. Cao and Vince work on the spectral radius of planar graphs and the related conjecture (Linear Algebra Appl. 187).
  • 2001. Zhou, Lin and Hu study the spectral radius of Hamiltonian planar and outerplanar graphs (Tsinghua Science and Technology 6).
  • 2017. Tait and Tobin prove the conjecture for all sufficiently large nnn (arXiv:1606.01916, J. Combin. Theory Ser. B, doi:10.1016/j.jctb.2017.04.006), Theorem 7.

This mission formalizes Section 2 of Tait and Tobin.

Setting

All graphs are finite and simple. An nnn-vertex graph GGG has vertex set V(G)={0,…,n−1}V(G) = \{0, \dots, n-1\}V(G)={0,…,n−1}, adjacency matrix A(G)A(G)A(G), and spectral radius λ1(G)\lambda_1(G)λ1​(G), the largest eigenvalue of the real symmetric matrix A(G)A(G)A(G). Write dud_udu​ for the degree of a vertex uuu and N(u)N(u)N(u) for its neighbourhood.

A plane drawing of GGG places the vertices at distinct points of R2\mathbb R^2R2 and each edge as a simple arc between its ends; an arc passes through no other vertex, and two arcs meet only at a common end. The union of the vertex points and arcs is a compact set, so its complement has exactly one unbounded connected component, the outer face. A graph is outerplanar if it has a plane drawing in which every vertex lies on the boundary of the outer face.

The join G+HG + HG+H of two graphs is their disjoint union together with all edges between them. The fan K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​ is a single vertex joined to every vertex of a path on n−1n-1n−1 vertices; it is outerplanar and has 2n−32n-32n−3 edges.

An outerplanar graph GGG on nnn vertices is a spectral maximizer if λ1(G′)≤λ1(G)\lambda_1(G') \le \lambda_1(G)λ1​(G′)≤λ1​(G) for every outerplanar graph G′G'G′ on the same nnn vertices.

Several milestones refer to a Perron vector: an eigenvector v\mathbf vv of A(G)A(G)A(G) for λ1\lambda_1λ1​ with all entries positive, normalized so that its maximum entry is 111, together with a vertex xxx with vx=1\mathbf v_x = 1vx​=1. For a connected graph such a vector exists by the Perron–Frobenius theorem.

Formalization targets

Goal: Theorem 7

There is NNN such that for every n≥Nn \ge Nn≥N and every graph GGG on nnn vertices,

G is an outerplanar spectral maximizer  ⟺  G≅K1+Pn−1.G \text{ is an outerplanar spectral maximizer} \iff G \cong K_1 + P_{n-1}.G is an outerplanar spectral maximizer⟺G≅K1​+Pn−1​.

The forward direction is the paper's Theorem 7; the converse says that the fan is the maximizer, as the conjecture reads. The threshold NNN is left unspecified, as in the paper.

Milestones

The milestones follow the proof in Section 2, in order:

  1. a strict subgraph of a connected graph has strictly smaller λ1\lambda_1λ1​ (already proved on the platform);
  2. an outerplanar graph on n≥2n \ge 2n≥2 vertices has at most 2n−32n-32n−3 edges and no K2,3K_{2,3}K2,3​ subgraph;
  3. Lemma 3: for n≥3n \ge 3n≥3 a spectral maximizer has λ1>n−1\lambda_1 > \sqrt{n-1}λ1​>n−1​;
  4. Lemma 4: du>vun−11nd_u > \mathbf v_u n - 11\sqrt ndu​>vu​n−11n​ for every vertex uuu;
  5. Lemma 5: for nnn large, dx>n−11nd_x > n - 11\sqrt ndx​>n−11n​ and vu<C1/n\mathbf v_u < C_1/\sqrt nvu​<C1​/n​ for every u≠xu \ne xu=x, with an absolute constant C1C_1C1​;
  6. Lemma 6: for nnn large, ∑z∈Bvz<C2/n\sum_{z\in B}\mathbf v_z < C_2/\sqrt n∑z∈B​vz​<C2​/n​, where B=V(G)∖(N(x)∪{x})B = V(G)\setminus(N(x)\cup\{x\})B=V(G)∖(N(x)∪{x});
  7. for nnn large, B=∅B = \emptysetB=∅: xxx is adjacent to every other vertex;
  8. an outerplanar graph with a vertex adjacent to all others is a spanning subgraph of K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​.

Significance

The theorem settles, for all large nnn, a conjecture that had been open since 1990 and had resisted several partial attacks. The same paper settles, in the same framework, the planar analogue (the Boots–Royle and Cao–Vince conjecture, K2+Pn−2K_2 + P_{n-2}K2​+Pn−2​) and a conjecture on spectral irregularity. The outerplanar case is the shortest of the three.

On the formal side, the result is proved on paper but has no machine-checked proof. Completing the mission produces a formal theory of outerplanar plane drawings (edge bound, K2,3K_{2,3}K2,3​-freeness, outerplanarity of the fan), eigenvector estimates for extremal graphs, and the full extremal theorem. The edge bound and the K2,3K_{2,3}K2,3​ exclusion are reusable for any later work on outerplanar or series–parallel graphs. Determining the exact threshold NNN, or proving the statement for every n≥2n \ge 2n≥2, would be a strengthening beyond the paper.

Difficulty

There are two separate difficulties. The first is topological: the facts about outerplanar graphs that the spectral argument uses (at most 2n−32n-32n−3 edges, no K2,3K_{2,3}K2,3​ subgraph, no cycle in G−xG - xG−x when xxx dominates, outerplanarity of the fan) are folklore on paper but require a Jordan-curve-type argument about plane drawings, of which Mathlib currently contains very little. The second is spectral: the edge-switching step that proves B=∅B = \emptysetB=∅ compares λ1\lambda_1λ1​ of two graphs through the Rayleigh quotient, and needs a positive Perron vector and the variational characterization of the largest eigenvalue for symmetric matrices. A purely combinatorial edge-maximization argument does not work: the fan is not the only outerplanar graph with 2n−32n-32n−3 edges, so maximizing the number of edges does not identify the spectral maximizer.

Formalization scope

  • Graphs are SimpleGraph (Fin n); λ1\lambda_1λ1​ is the published WangKangXue.SpectralTuran.specRad, the largest eigenvalue of G.adjMatrix ℝ.
  • Outerplanarity is topological: a published RobertsonSeymour1986.GM5.PlaneDrawing with every vertex point in the closure of the unbounded component of the complement of the drawing. It is not replaced by a forbidden-minor characterization; such characterizations are theorems about outerplanarity, and would change which graphs compete in the maximization.
  • Extremality is IsSpecMax IsOuterplanar G: GGG is outerplanar and no outerplanar graph on the same vertex set has larger λ1\lambda_1λ1​. It is never weakened to a comparison with the fan.
  • "GGG is K1+Pn−1K_1 + P_{n-1}K1​+Pn−1​" is a graph isomorphism G ≃g fan n, with the fan on Fin 1 ⊕ Fin (n - 1).
  • Perron-vector statements quantify over every positive eigenvector for λ1\lambda_1λ1​ with maximum entry 111 and every vertex xxx with vx=1\mathbf v_x = 1vx​=1; existential quantification over v\mathbf vv would be weaker.
  • "For some absolute constant" is an existential quantifier placed before the threshold, the graph and the vector; "for nnn sufficiently large" is ∃ N, ∀ n ≥ N.
  • Added hypotheses, each forced by a degenerate case: n≥2n \ge 2n≥2 for the edge bound (the bound 2n−32n-32n−3 is negative below), n≥3n \ge 3n≥3 for Lemma 3 (it fails for n=2n = 2n=2, where the maximizer K2K_2K2​ has λ1=1=n−1\lambda_1 = 1 = \sqrt{n-1}λ1​=1=n−1​).

A complete development needs the Perron–Frobenius theorem for symmetric nonnegative matrices and the Rayleigh quotient characterization of λ1\lambda_1λ1​, both of general use; the theory of outerplanar drawings; and the counting arguments of the lemmas. Contributions to any of these layers are welcome, as are independent proofs of the outerplanar facts from the drawing definition.

Selected references

  • M. Tait and J. Tobin, Three conjectures in extremal spectral graph theory, J. Combin. Theory Ser. B (2017), https://doi.org/10.1016/j.jctb.2017.04.006; preprint arXiv:1606.01916v2, https://arxiv.org/abs/1606.01916
  • D. Cvetković and P. Rowlinson, The largest eigenvalue of a graph: a survey, Linear and Multilinear Algebra 28 (1990), 3–33.
  • P. Rowlinson, On the index of certain outerplanar graphs, Ars Combinatoria 29 (1990), 221–225.
  • D. Cao and A. Vince, The spectral radius of a planar graph, Linear Algebra Appl. 187 (1993), 251–257.
  • J. Zhou, C. Lin and G. Hu, Spectral radius of Hamiltonian planar graphs and outerplanar graphs, Tsinghua Science and Technology 6 (2001), 350–354.
12 thms2 active usersReviewed
🏆Completed
CombinatoricsProbabilityStochastic Systems·Captain: mikedeng1

Separation of Time-Scales and Model Reduction for Stochastic Reaction Networks 1: The Balance Condition for Linear Combinations of Species Decomposes over Strongly Connected ComponentsResearch Paper

Why reaction-network balance matters

A stochastic reaction network tracks the counts of chemical species as reactions occur. Different species and reaction rates may grow at different powers of a large system-size parameter. Kang and Kurtz's model-reduction paper uses these powers to identify candidate time scales and combinations of species that can remain bounded. Before a reduced model can be justified, the dominant reactions that increase a combination must be compared with those that decrease it. Condition 3.2 makes that comparison.

Checking the condition for every nonnegative linear combination is a large task even for a finite network: the weights range over a continuum, and changing a weight can change which reactions increase or decrease the combination. The paper's Lemma 3.4 relates this task to the directed graph of species. It says that combinations supported on different strongly connected components can be checked component by component. This mission formalizes that result and the claims in its proof that carry the reduction.

The network and its graph

There are finitely many species SiS_iSi​ and finitely many indexed reactions kkk. Reaction kkk consumes νik\nu_{ik}νik​ molecules of species iii and produces νik′\nu'_{ik}νik′​ molecules of that species. The stoichiometric change is ζik=νik′−νik\zeta_{ik}=\nu'_{ik}-\nu_{ik}ζik​=νik′​−νik​. Species iii has a nonnegative scaling exponent αi\alpha_iαi​, and reaction kkk has a real exponent βk\beta_kβk​. The effective reaction exponent is ρk=βk+∑iνikαi\rho_k=\beta_k+\sum_i\nu_{ik}\alpha_iρk​=βk​+∑i​νik​αi​.

For a vector of nonnegative weights θ\thetaθ, the scalar θ⋅ζk\theta\cdot\zeta_kθ⋅ζk​ is the change in its linear combination when reaction kkk occurs. Write Γθ+\Gamma^+_\thetaΓθ+​ for reactions with positive change and Γθ−\Gamma^-_\thetaΓθ−​ for reactions with negative change. The balance equation equates the largest ρk\rho_kρk​ on these two sets. If that equation fails, the time-scale inequality permits a chosen exponent γ\gammaγ only when it is no larger than γθ\gamma_\thetaγθ​, the largest αi\alpha_iαi​ on the support of θ\thetaθ minus the largest ρk\rho_kρk​ on Γθ+∪Γθ−\Gamma^+_\theta\cup\Gamma^-_\thetaΓθ+​∪Γθ−​. Condition 3.2 is the disjunction of these two alternatives.

The species graph has an edge Si→SjS_i\to S_jSi​→Sj​ when one reaction consumes SiS_iSi​ and produces SjS_jSj​. Its strongly connected components are the maximal sets of species that can reach one another by directed paths. A single species is a component when it shares no larger mutual-reachability class. A nonnegative vector supported in one such component is a component-level test vector.

Formalization targets

The goal is Lemma 3.4, p. 15. Suppose θ=∑j=1mθj\theta=\sum_{j=1}^m\theta^jθ=∑j=1m​θj, each θj\theta^jθj is nonnegative and supported within one maximal strongly connected component, and distinct pieces belong to distinct components. For fixed γ\gammaγ, the lemma asserts all three implications:

(∀j,  Cond3.2(γ,θj))⟹Cond3.2(γ,θ),(∀j,  Balance(θj))⟹Balance(θ),(∀j,  γ≤γθj)⟹γ≤γθ.\begin{aligned} \bigl(\forall j,\;\mathrm{Cond}_{3.2}(\gamma,\theta^j)\bigr)&\Longrightarrow\mathrm{Cond}_{3.2}(\gamma,\theta),\\ \bigl(\forall j,\;\mathrm{Balance}(\theta^j)\bigr)&\Longrightarrow\mathrm{Balance}(\theta),\\ \bigl(\forall j,\;\gamma\le\gamma_{\theta^j}\bigr)&\Longrightarrow\gamma\le\gamma_\theta. \end{aligned}(∀j,Cond3.2​(γ,θj))(∀j,Balance(θj))(∀j,γ≤γθj​)​⟹Cond3.2​(γ,θ),⟹Balance(θ),⟹γ≤γθ​.​

The companion consequence on the same page quantifies over all nonnegative vectors: if any one of these properties holds for every vector supported in a strongly connected subgraph, then it holds for every vector. The milestones mirror the paper's proof: a chain of component pieces cannot revisit a piece; under Condition 3.2 the maximal positive exponent either satisfies the time-scale bound or is bounded by the maximal negative exponent; and the two uniform cases yield their respective alternatives.

What the result establishes

Lemma 3.4 replaces a global family of tests with tests inside individual strongly connected components. In the paper this narrows the balance checks needed before the authors study separated time scales and reduced stochastic models. The uniform balance and uniform time-scale clauses matter separately: the first preserves equality of dominant positive and negative exponents, while the second preserves the same admissible γ\gammaγ across a sum of pieces. The general Condition 3.2 clause allows different pieces to satisfy different alternatives.

The mathematical result was proved in the 2010 arXiv preprint and published in 2013; it is not an open conjecture. This proposal supplies Lean statements and a reusable finite reaction-network vocabulary, with proofs left for solvers. The local sanity file checks the empty-maximum and zero-vector behavior without using sorry. The draft theorem files compile as statements, but their six proofs are not yet machine checked.

Where the reasoning is difficult

The obvious inequality for a sum says that a reaction increasing θ\thetaθ increases at least one piece. It does not say that the piece's dominant increasing reaction also increases the sum: contributions from other pieces can cancel it. Consequently, taking a maximum separately on each piece does not directly give the balance equation for the sum. The graph restriction is decisive because a sequence of compensating reactions would otherwise be able to cycle among pieces. The milestones isolate that no-repetition property and the two inequalities the paper obtains from it.

Formalization scope

Species, reactions, and pieces are represented by Fin s, Fin r, and Fin m; these are zero-based versions of the paper's one-based indices. The stoichiometric coefficients are natural numbers, the changes ζk\zeta_kζk​ and all scaling exponents are real, and αi≥0\alpha_i\ge0αi​≥0. Following the standing restriction on p. 5, each reaction consumes at most two molecules. The paper's ρk\rho_kρk​ is computed from βk\beta_kβk​, νk\nu_kνk​, and α\alphaα, not supplied independently.

For a reaction set AAA, max⁡k∈Aρk\max_{k\in A}\rho_kmaxk∈A​ρk​ uses WithBot ℝ, so the maximum of an empty set is −∞-\infty−∞, as the paper explicitly uses on p. 17. The support maximum in γθ\gamma_\thetaγθ​ is an ordinary real maximum when relevant. The value γθ\gamma_\thetaγθ​ is represented in EReal: when Γθ+∪Γθ−\Gamma^+_\theta\cup\Gamma^-_\thetaΓθ+​∪Γθ−​ is empty, it is +∞+\infty+∞ and the time-scale inequality holds for every real γ\gammaγ. Subtraction is performed in the reals before conversion to EReal. Thus zero vectors do not acquire a spurious bound from a default-valued maximum or extended-real subtraction.

An edge uses positive consumption and positive production in the same indexed reaction, matching the wording on pp. 4 and 15. Mutual reachability defines a component. Each piece is required to lie in one component, and different pieces in different components. These conditions prevent splitting a single component into separately balanced pieces, which would make the goal false. An empty family of pieces is permitted; its sum is zero. The definition file is useful for the other balance results in the paper, while the graph and component lemmas can support later work on reaction-network decomposition.

Selected references

  • Hye-Won Kang and Thomas G. Kurtz, Separation of time-scales and model reduction for stochastic reaction networks, arXiv:1011.1672v1, 2010, especially §2 and §§3.2–3.4, pp. 4–17. Preprint.
  • Hye-Won Kang and Thomas G. Kurtz, Separation of time-scales and model reduction for stochastic reaction networks, Annals of Applied Probability 23(2), 2013. DOI. The numbered statements in this mission follow the pinned 2010 preprint.
6 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: mikedeng1

Cycle Lengths in Expanding Graphs 2: Every α-Expander on n Vertices Has Ω(α³/log(1/α))·n Distinct Cycle LengthsResearch Paper

Motivation

An expander graph has many edges leaving each small vertex set. This condition is useful when a graph must remain connected after local losses, and it also constrains the cycles that can occur. Knowing that one long cycle exists leaves open how many different lengths its cycles can have. For example, a bipartite graph has no odd cycle, so expansion alone does not make every length available. The question here is quantitative: how large must the set of available lengths be in terms of the graph's order and expansion parameter? Friedman and Krivelevich answered it with a lower bound whose dependence on the expansion parameter is polynomial up to a logarithm. Their Theorem 2 is the target of this mission. Friedman and Krivelevich, Cycle lengths in expanding graphs, pp. 2–3

The same paper separates two questions. Theorem 1 says that cycle lengths are well spread over a range, but its resulting lower bound on their number has a much smaller dependence on the expansion parameter. Theorem 2 obtains a stronger count without claiming that those lengths form an interval. This distinction matters: counting many available lengths and prescribing an arbitrary length are different demands on a graph. The paper's third theorem studies the latter demand under a different graph condition. Friedman and Krivelevich, pp. 2–4

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite simple graph with n=∣V∣n=|V|n=∣V∣ vertices. A simple graph has no loops or parallel edges. For a vertex set U⊆VU\subseteq VU⊆V, its external neighborhood is

NG(U)={v∈V∖U:v is adjacent to some vertex in U}.N_G(U)=\{v\in V\setminus U:\text{$v$ is adjacent to some vertex in $U$}\}.NG​(U)={v∈V∖U:v is adjacent to some vertex in U}.

The graph is an α\alphaα-expander if α>0\alpha>0α>0 and, for every U⊆VU\subseteq VU⊆V with ∣U∣≤⌈n/2⌉|U|\le\lceil n/2\rceil∣U∣≤⌈n/2⌉,

∣NG(U)∣≥α∣U∣.|N_G(U)|\ge\alpha|U|.∣NG​(U)∣≥α∣U∣.

The ceiling is part of the paper's formal Definition on p. 2. It matters at odd nnn, where the largest tested set has (n+1)/2(n+1)/2(n+1)/2 vertices. A related (k,α)(k,\alpha)(k,α)-expander uses the same inequality for every UUU with ∣U∣≤k|U|\le k∣U∣≤k; here kkk may be real and the size of UUU is still an integer. This broader notion appears in the supporting long-path lemma. Friedman and Krivelevich, pp. 1–2, 4

A cycle is a closed path with no repeated vertex apart from its starting and ending vertex. Its length is its number of edges, hence at least three. Write

L(G)={ℓ∈N:G contains a cycle of length ℓ}.L(G)=\{\ell\in\mathbb N:\text{$G$ contains a cycle of length $\ell$}\}.L(G)={ℓ∈N:G contains a cycle of length ℓ}.

The quantity ∣L(G)∣|L(G)|∣L(G)∣ counts lengths, rather than the number of cycles. Two cycles of the same length contribute just one element. Since a simple cycle visits distinct vertices, L(G)L(G)L(G) is finite and its cardinality is an ordinary natural number. Friedman and Krivelevich, p. 2

Formalization targets

The goal is the polynomial lower bound from Theorem 2. There is one absolute constant c>0c>0c>0 such that, for each 0<α≤10<\alpha\le10<α≤1 and every sufficiently large nnn, every α\alphaα-expander GGG on nnn vertices satisfies

∣L(G)∣≥c α3log⁡2(2/α) n.|L(G)|\ge c\,\frac{\alpha^3}{\log_2(2/\alpha)}\,n.∣L(G)∣≥clog2​(2/α)α3​n.

The paper writes Ω(α3/log⁡(1/α))n\Omega(\alpha^3/\log(1/\alpha))nΩ(α3/log(1/α))n. The denominator above records the endpoint convention used in this mission; it agrees up to an absolute factor for α≤1/2\alpha\le1/2α≤1/2. The constant ccc is uniform across all α\alphaα, graph orders, and graphs. The threshold for nnn may depend on α\alphaα. Friedman and Krivelevich, p. 3, Theorem 2

Four intermediate targets reproduce claims used by the paper. Lemma 2.5 gives a path of a prescribed length from a prescribed start vertex when every kkk-vertex set has enough external neighbors and the start vertex's component is large enough. In the proof of Theorem 2, deleting ⌊αn/4⌋\lfloor\alpha n/4\rfloor⌊αn/4⌋ vertices leaves a path of length ⌈αn/4⌉−1\lceil\alpha n/4\rceil-1⌈αn/4⌉−1. The chosen integer

k(α)=2⌈log⁡2(3/α)log⁡2(1+α/2)⌉+1k(\alpha)=2\left\lceil\frac{\log_2(3/\alpha)}{\log_2(1+\alpha/2)}\right\rceil+1k(α)=2⌈log2​(1+α/2)log2​(3/α)​⌉+1

has order at most log⁡2(2/α)/α\log_2(2/\alpha)/\alphalog2​(2/α)/α with an absolute constant. Finally, the paper's set X0X_0X0​ of path vertices connected to the deleted set by short paths satisfies ∣X0∣≥α2n/12|X_0|\ge\alpha^2n/12∣X0​∣≥α2n/12 once the rounding-sensitive small cases are excluded. These are the milestone statements, in the order in which they enter the target. Friedman and Krivelevich, pp. 6, 13–14

Significance

The theorem guarantees linearly many distinct cycle lengths for every fixed positive expansion parameter, and gives an explicit polynomial-scale dependence on that parameter. A long cycle alone proves only that one length occurs. Even a statement placing cycles near many target lengths can lose substantial strength in its dependence on α\alphaα. The count in Theorem 2 supplies a separate guarantee for graph families where the expansion parameter varies. It does not imply that every length in an interval occurs. Friedman and Krivelevich, p. 3

The result is proved in the cited paper. The formalization target is a machine-checked proof of that known statement and its selected intermediate claims. The reusable output includes an exact interface for external vertex expansion, simple cycle lengths, and the prescribed-start path lemma. The local Lean statements currently compile as open targets; they do not supply machine-checked proofs of the claims.

Difficulty

Expansion gives large neighborhoods and forces connectedness, but these facts do not by themselves distinguish cycle lengths. A tempting first argument finds a long cycle and then treats its subpaths as shorter cycles; the subpaths do not close into cycles without additional edges. A second temptation is to count cycles rather than lengths: a graph may have many cycles while all of them share only a few lengths. The central challenge is to obtain a family of cycles whose lengths can be certified as different while keeping the graph's expansion parameter visible in the count. The proof also has integer rounding at the scale αn\alpha nαn, where claims written for large graphs can fail on small ones. Friedman and Krivelevich, pp. 13–14

Formalization scope

Graphs are SimpleGraph (Fin n) for the goal and proof-specific milestones; the general path lemma allows any finite vertex type. Vertex sets are Set V, with Set.ncard for finite cardinality. An α\alphaα-expander tests every set up to ⌈n/2⌉=(n+1)/2\lceil n/2\rceil=(n+1)/2⌈n/2⌉=(n+1)/2 in natural-number arithmetic. Cycles use SimpleGraph.Walk.IsCycle, and paths use Walk.IsPath; their length is the number of edges. This rules out a two-edge closed walk being counted as a cycle. The formal development also uses the explicit k(α)k(\alpha)k(α) and the paper's X0X_0X0​, including its condition on the internal vertices of a short connecting path.

All logarithms are binary. The shifted log⁡2(2/α)\log_2(2/\alpha)log2​(2/α) replaces log⁡2(1/α)\log_2(1/\alpha)log2​(1/α) only inside asymptotic bounds, because the latter vanishes at α=1\alpha=1α=1. The constant in the lower bound precedes α\alphaα in the quantifier order. The goal includes n≥n0(α)n\ge n_0(\alpha)n≥n0​(α), a necessary disclosed addition: the two-vertex complete graph is an α\alphaα-expander for α≤1\alpha\le1α≤1 but has no cycle. The X0X_0X0​ milestone adds αn≥12\alpha n\ge12αn≥12 to make the paper's suppressed rounding valid; otherwise the initial vertex set can be empty. The deletion milestone states the path claim for any deleted set of the indicated size, since the displayed estimate uses only that size. A formal proof needs finite graph paths, component sizes, vertex neighborhoods, and cardinality estimates. Contributions to these general lemmas are useful beyond this mission. No claim is discharged by letting an asymptotic constant depend on α\alphaα, by counting closed walks, or by exploiting division by zero.

Selected references

  • L. Friedman and M. Krivelevich, Cycle lengths in expanding graphs, arXiv:1912.11011v2, 2020; published in Combinatorica 41 (2021). arXiv; DOI.
7 thms2 active usersReviewed
🏆Completed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

A c^k n 5-Approximation Algorithm for Treewidth 2: X Is a Balanced S-Separator iff V(G)∖X Splits into Three Mutually Non-Adjacent Parts Each Holding at Most |S|/2 Vertices of SResearch Paper

Motivation

The paper of Bodlaender, Drange, Dregi, Fomin, Lokshtanov and Pilipczuk (2016) has an algorithmic headline. This mission concerns one of its finite graph lemmas, Lemma 2.9. The lemma translates a condition phrased separately for every connected component into a condition on three vertex sets. A reader of graph decompositions encounters both forms: components express what remains connected after deleting vertices, while a small fixed number of parts gives a compact way to state constraints on a partition. Lemma 2.9 identifies exactly when the two forms agree for the half-balance threshold. The same result appears again in the paper as Lemma 6.2, where the authors use it in a later section.

The result does not require a tree decomposition or a treewidth hypothesis. It holds for every finite graph. That generality makes it a useful standalone target: the definitions of deletion, connected components, and the weighted three-partition can be reused wherever graph separators are studied. The mission follows the published SIAM article, with printed page 328 for the main lemma and pages 363–364 for its restatement and the supporting combinatorial observation.

Setting

Let GGG be a finite simple graph with vertex set V(G)V(G)V(G). Let SSS and XXX be subsets of its vertices. The deleted graph G∖XG\setminus XG∖X retains precisely the vertices outside XXX and the edges joining them. Two remaining vertices are in the same connected component when a walk joins them without passing through any vertex of XXX. Write CX(u)C_X(u)CX​(u) for the vertex set of the component containing u∉Xu\notin Xu∈/X. This interpretation follows the paper's notation on page 320.

A balanced SSS-separator is a set XXX for which each component of G∖XG\setminus XG∖X contains at most half the vertices of SSS:

∣CX(u)∩S∣≤∣S∣2for every u∈V(G)∖X.|C_X(u)\cap S|\le\frac{|S|}{2}\qquad\text{for every }u\in V(G)\setminus X.∣CX​(u)∩S∣≤2∣S∣​for every u∈V(G)∖X.

The paper introduces this phrase on page 321. Here balance is a property of XXX relative to SSS; it imposes no bound on the size of XXX. The denominator counts all vertices of SSS, including any that lie in XXX.

A three-part partition of V(G)∖XV(G)\setminus XV(G)∖X is a triple (M1,M2,M3)(M_1,M_2,M_3)(M1​,M2​,M3​) of pairwise disjoint sets whose union is V(G)∖XV(G)\setminus XV(G)∖X. A part may be empty. Two parts have no edge between them if no edge of GGG joins a vertex of one to a vertex of the other. Since the graph is undirected, checking the three unordered pairs of distinct parts covers the requirement for every i≠ji\ne ji=j.

Formalization targets

Lemma 2.9: balanced separators and three parts

The goal is the equivalence stated on page 328, also restated as Lemma 6.2 on page 363:

X is a balanced S-separator⟺V(G)∖X=M1∪˙M2∪˙M3,no edge joins distinct Mi,∣Mi∩S∣≤∣S∣/2(i=1,2,3).X\text{ is a balanced }S\text{-separator} \quad\Longleftrightarrow\quad \begin{gathered} V(G)\setminus X=M_1\mathbin{\dot\cup}M_2\mathbin{\dot\cup}M_3,\\ \text{no edge joins distinct }M_i,\\ |M_i\cap S|\le |S|/2\quad(i=1,2,3). \end{gathered}X is a balanced S-separator⟺V(G)∖X=M1​∪˙M2​∪˙M3​,no edge joins distinct Mi​,∣Mi​∩S∣≤∣S∣/2(i=1,2,3).​

The equivalence is the target, including both directions. It applies even when SSS or V(G)∖XV(G)\setminus XV(G)∖X is empty. It does not include a bound on ∣X∣|X|∣X∣.

Lemma 6.3: three bounded sums

The main supporting result on page 363 says that nonnegative integers a1,…,apa_1,\ldots,a_pa1​,…,ap​ with total qqq and ai≤q/2a_i\le q/2ai​≤q/2 can be divided into three classes, each with sum at most q/2q/2q/2. The milestone list also records the paper's merge inequality and its assertion that a component cannot cross between parts having no edge between them. These are the specific claims used in the published discussion of Lemma 2.9.

Significance

Lemma 2.9 lets a componentwise balance condition be represented by three sets whose coverage, disjointness, edge separation, and cardinalities can be inspected directly. It gives an exact equivalence, so the three-part formulation loses no balanced separators and introduces none that fail the original definition. That is the mathematical content needed when a later argument represents a separator by a bounded number of regions rather than by a list of all its components. The result is already proved in the 2016 paper; the open work in this mission is a machine-checked formalization of that known result and its cited supporting claims.

Formalizing the equivalence also establishes reusable interfaces for vertex deletion and component membership in a finite simple graph. A future development can state a separator condition using IsBalancedSSep, or use the three-part conclusion without redefining what a component of G∖XG\setminus XG∖X means. The independent integer-partition lemma has uses whenever indivisible nonnegative weights, each no larger than half the total, must be grouped into three bounded classes. Nothing in that statement depends on graphs.

Difficulty

There may be arbitrarily many components of G∖XG\setminus XG∖X, although the conclusion permits only three parts. Simply assigning all components to one part can violate the half bound, while splitting a component between different parts may create forbidden edges. A component's weight is the number of its vertices in SSS, and those weights can be uneven or zero. The theorem has to accommodate all such patterns without changing the threshold ∣S∣/2|S|/2∣S∣/2. In the converse direction, the absence of edges between the three parts must control entire walks, not only their first edge; ordinary connectivity in GGG would allow a walk through XXX and would describe the wrong components.

Formalization scope

The Lean development uses SimpleGraph V with V : Type, [Fintype V], and [DecidableEq V]. Vertex sets are Finset V. The shared walk-based AvoidReach G X u v requires every vertex on a walk to lie outside XXX, including both endpoints; avoidComp G X u collects the vertices reachable in this way. It is used only for u∉Xu\notin Xu∈/X when quantifying over components. This shared definition matches the paper's G∖XG\setminus XG∖X. The graph itself need not have decidable adjacency, since the statements concern propositions rather than computations.

Half bounds are written as 2c≤∣S∣2c\le |S|2c≤∣S∣ or 2c≤q2c\le q2c≤q in natural-number arithmetic. This is equivalent to the paper's rational half bound and avoids any interpretation of natural-number division. The three parts may be empty, so X=V(G)X=V(G)X=V(G) is a genuine boundary case. The set SSS may meet XXX; the right-hand denominator remains ∣S∣|S|∣S∣. The graph may have no vertices. There is no treewidth parameter or size condition on XXX: adding one would change Lemma 2.9 and can make its reverse implication false. Nor may components be computed using connectivity in the original graph, where paths through deleted vertices could merge them.

The required infrastructure is finite graph walks, finite sets, cardinalities, sums over finite index types, and the three original claims in the milestone list. The graph-component definitions and the integer grouping theorem are reusable beyond this mission. Contributions that prove any milestone or the full equivalence under exactly these conventions are in scope.

Selected references

  • H. L. Bodlaender, P. G. Drange, M. S. Dregi, F. V. Fomin, D. Lokshtanov and M. Pilipczuk, A c^k n 5-approximation algorithm for treewidth, SIAM Journal on Computing 45(2):317–378, 2016. DOI: 10.1137/130947374.
6 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: mikedeng1

Generation and Properties of Snarks 4: Every 3-Edge-Colourable Cubic Graph with m Edges Has a Cycle Cover of Total Length at Most 4m/3Research Paper

Motivation

A cycle cover of a graph is a collection of cycles that together use every edge at least once. Covering the edges of a network by closed tours appears in routing and inspection problems (the Chinese postman problem is the closed-walk version), and the minimum total length of a cycle cover is a classical extremal parameter of bridgeless graphs. Alon and Tarsi (Covering multigraphs by simple circuits, SIAM J. Algebraic Discrete Methods 6, 1985, pp. 345–350; reference [2] of the paper) conjectured that every bridgeless graph with mmm edges has a cycle cover of length at most 7m/57m/57m/5 (Conjecture 7.1 of the paper).

The question is closely tied to cycle double covers (CDCs): families of cycles covering every edge exactly twice. The Cycle Double Cover Conjecture of Szekeres and Seymour asserts that every bridgeless graph has one, and it suffices to prove it for cubic graphs. Brinkmann, Goedgebeur, Hägglund and Markström (Generation and properties of snarks, arXiv:1206.6690v3, 2013; J. Combin. Theory Ser. B 103, 2013) generated all snarks on up to 36 vertices and tested these conjectures on them. Section 7 of their paper derives cycle-cover bounds from double covers. For colourable cubic graphs (those with a proper 3-edge-colouring) the argument gives the bound 4m/34m/34m/3, and the authors conjecture (Conjecture 7.4) that 4m/3+o(m)4m/3 + o(m)4m/3+o(m) holds for all snarks. This mission formalizes the colourable bound and the chain of statements behind it.

Setting

All graphs are finite and simple. A graph is cubic if every vertex has degree 3. It is colourable if its edges can be coloured with three colours so that edges sharing a vertex receive different colours.

A cycle is a 2-regular connected graph. Inside a graph GGG a cycle is recorded by its edge set C={v1v2,…,vℓ−1vℓ,vℓv1}C = \{v_1v_2, \dots, v_{\ell-1}v_\ell, v_\ell v_1\}C={v1​v2​,…,vℓ−1​vℓ​,vℓ​v1​}, with v1,…,vℓv_1, \dots, v_\ellv1​,…,vℓ​ distinct and ℓ≥3\ell \ge 3ℓ≥3; its length is ∣C∣=ℓ|C| = \ell∣C∣=ℓ. A 2-regular subgraph of GGG is an edge set DDD such that every vertex lies on none or exactly two edges of DDD: a union of vertex-disjoint cycles, possibly empty and not necessarily spanning. A 2-factor is a spanning 2-regular subgraph. An even subgraph is an edge set SSS such that every vertex lies on an even number of edges of SSS.

A cycle double cover of GGG is a multiset of cycles such that every edge lies in exactly two of them. It is a kkk-CDC if the cycles can be coloured with kkk colours so that cycles sharing an edge get different colours. Each colour class is then an even subgraph, and following the paper the kkk-CDC is written as a family D0,…,Dk−1D_0, \dots, D_{k-1}D0​,…,Dk−1​ of even subgraphs with every edge in exactly two of them, counted with multiplicity.

A cycle cover of GGG is a set F\mathcal FF of cycles of GGG such that every edge lies in at least one of them, and its length is ℓ(F)=∑C∈F∣C∣\ell(\mathcal F) = \sum_{C \in \mathcal F} |C|ℓ(F)=∑C∈F​∣C∣. Throughout, mmm denotes the number of edges of GGG.

Formalization targets

Goal: the 4m/34m/34m/3 bound (§7, p. 24)

For every colourable cubic graph GGG with mmm edges there is a cycle cover F\mathcal FF with

ℓ(F)≤4m3.\ell(\mathcal F) \le \frac{4m}{3}.ℓ(F)≤34m​.

Equivalently, the shortest cycle cover of GGG has length at most 4m/34m/34m/3.

Milestones

  1. Cycle decompositions (§5, p. 15). An edge set SSS is an even subgraph of GGG if and only if it is the union of pairwise edge-disjoint cycles of GGG.
  2. Proposition 5.3 (p. 15). If GGG is a colourable cubic graph and DDD is any 2-regular subgraph of GGG, then GGG has a 4-CDC (D0,D1,D2,D3)(D_0, D_1, D_2, D_3)(D0​,D1​,D2​,D3​) with D0=DD_0 = DD0​=D.
  3. Lemma 7.2 (p. 24). If a cubic graph GGG has a 2-regular subgraph CCC and a CDC D\mathcal DD with C∈DC \in \mathcal DC∈D, then GGG has a cycle cover F\mathcal FF with
ℓ(F)≤2m−∣C∣.\ell(\mathcal F) \le 2m - |C|.ℓ(F)≤2m−∣C∣.

Significance

The bound 4m/34m/34m/3 is far below the Alon–Tarsi bound 7m/57m/57m/5 and is attained up to lower-order terms on the snarks the paper enumerates (Corollary 7.3, which rests on a computer search and is not part of this mission). It isolates the colourable case of Conjecture 7.4, so any proof of that conjecture must at least reproduce it. Proposition 5.3 is a strong form of the Strong Cycle Double Cover Conjecture for colourable cubic graphs: every 2-regular subgraph, not only every cycle, is a colour class of a 4-CDC. Lemma 7.2 is the general device that turns a double cover containing a long 2-regular subgraph into a short cycle cover; the paper applies it again with long cycles of snarks.

On the formal side, none of these objects exists in Mathlib or on the platform: there is no cycle double cover, no cycle cover, and no decomposition of even graphs into cycles (Veblen's theorem). Mathlib has walks, cycles as walks (Walk.IsCycle), Eulerian trails and the predicate IsCycles, but the step from an even edge set to a family of edge-disjoint cycles is missing. The results of this mission are all proved in the paper or classical; none has a machine-checked proof.

Difficulty

The arithmetic of the goal is immediate once the pieces exist: a 3-edge-colouring of a cubic graph on nnn vertices consists of three perfect matchings, two of which form a 2-factor with n=2m/3n = 2m/3n=2m/3 edges, and 2m−2m/3=4m/32m - 2m/3 = 4m/32m−2m/3=4m/3. The work lies elsewhere.

First, Veblen's theorem must be formalized: an even edge set is a union of edge-disjoint cycles. This requires extracting a cycle (a walk with no repeated vertex) from a nonempty even edge set and inducting on the number of edges; the bookkeeping between walks, their edge lists and finite edge sets is where formal proofs of this kind spend their effort. Lemma 7.2 needs it, because the colour classes of a CDC are even subgraphs, while a cycle cover consists of cycles.

Second, Proposition 5.3 needs the passage from a proper 3-edge-colouring (a colouring of the line graph) to three perfect matchings and to the three 2-factors formed by pairs of colours, and then a parity check that the symmetric differences with DDD are even and cover every edge exactly twice. A tempting shortcut fails: the three pairwise unions of colour classes already form a 3-CDC, but a prescribed 2-regular subgraph DDD need not be one of its classes.

Formalization scope

Graphs are SimpleGraph V on a finite type V with decidable equality and decidable adjacency. Cubic is G.IsRegularOfDegree 3. Colourable is G.lineGraph.Colorable 3, a proper vertex colouring of the line graph, never a vertex colouring of G. Edges are Sym2 V, and every subgraph is recorded by its finite edge set. A cycle is the edge set of a closed walk satisfying Walk.IsCycle, of length its number of edges. A kkk-CDC is a function Fin k → Finset (Sym2 V) of even edge sets with every edge of G in exactly two indices, so repeated and empty classes count correctly. A cycle cover is a Finset of cycles covering every edge, of length ∑ C ∈ F, C.card.

"The shortest cycle cover has length at most XXX" is stated as "some cycle cover has length at most XXX"; the two agree because a finite graph has finitely many cycle covers. Bounds are kept in N\mathbb NN without subtraction or division: the goal is 3 ℓ(F)≤4m3\,\ell(\mathcal F) \le 4m3ℓ(F)≤4m, and Lemma 7.2 is ℓ(F)+∣C∣≤2m\ell(\mathcal F) + |C| \le 2mℓ(F)+∣C∣≤2m. Lemma 7.2 is stated for a CDC given as a family of even subgraphs, the form in which the paper applies it (with CCC a 2-factor that is a colour class of a 4-CDC); a CDC of single cycles is a special case. The cubic hypothesis of Lemma 7.2 is kept as printed although the argument does not use it.

The bound becomes trivial if a "cycle" is weakened to an even edge set or to the edge set of a closed walk, since then the single set E(G)E(G)E(G) is a cover of length m≤4m/3m \le 4m/3m≤4m/3. The cycles here are genuine cycles of G.

Out of scope: Corollary 7.3 (the cycle-cover lengths of all snarks on at most 36 vertices) and every other statement resting on the exhaustive enumeration, Conjectures 7.1 and 7.4, and the Strong Cycle Double Cover material of §5.1 and §6.

Reusable infrastructure: cycles as edge sets, even subgraphs, kkk-CDCs, cycle covers and Veblen's theorem are general graph-theoretic objects, useful well beyond this paper (for flows, the Cycle Double Cover Conjecture, and Eulerian subgraph arguments). Proofs of Veblen's theorem and of the passage from a 3-edge-colouring to perfect matchings are particularly welcome.

Selected references

  • G. Brinkmann, J. Goedgebeur, J. Hägglund, K. Markström, Generation and properties of snarks, arXiv:1206.6690v3, 2013; J. Combin. Theory Ser. B 103 (2013). https://arxiv.org/abs/1206.6690
  • N. Alon, M. Tarsi, Covering multigraphs by simple circuits, SIAM J. Algebraic Discrete Methods 6 (1985), no. 3, 345–350; MR 791163. Cited as [2] in https://arxiv.org/abs/1206.6690
10 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: marwahaha

A logarithmic independence bound for clique-free graphsOpen Problem

Motivation

Excluding a fixed clique should force independent sets larger than those guaranteed by average degree alone. The source asks for the full logarithmic improvement uniformly across graph sizes. The source is OpenAI's September 2026 manuscript.

Setting

An independent set contains no adjacent pair. Write nnn for the number of vertices, d=2∣E∣/nd=2|E|/nd=2∣E∣/n for average degree, and α(G)α(G)α(G) for the largest independent-set size. A graph is clique-free of order rrr when it contains no rrr pairwise adjacent vertices.

Formalization targets

∀r≥4 ∃cr>0:α(G)≥crnlog⁡dd(d≥2, G is Kr-free).\forall r\ge4\ \exists c_r>0:\quad \alpha(G)\ge c_r\frac{n\log d}{d}\qquad(d\ge2,\ G\text{ is }K_r\text{-free}).∀r≥4 ∃cr​>0:α(G)≥cr​dnlogd​(d≥2, G is Kr​-free).

For every natural number r ≥ 4, there is a real constant c > 0 such that the following holds for every finite simple graph G on a vertex type V (a Fintype). Here the average degree d(G) is defined as 2|E(G)|/|V|, twice the number of edges divided by the number of vertices (as a real number). If G contains no clique of size r (it is r-clique-free) and d(G) ≥ 2, then the independence number α(G), the maximum size of an independent set, satisfies α(G) ≥ c · |V| · log(d(G)) / d(G), where log is the natural logarithm of real numbers. The constant c may depend on r but not on G.

The goal is OAI.CliqueFreeLog.logarithmic_independence_bound.

Significance

The bound improves the general inverse-degree scale by a factor of natural logarithmic size. The constant may depend on the forbidden clique but must work for every input graph. The manuscript presents a mathematical argument; the attached goal currently has status Open, so the requested contribution is a Lean proof of the displayed formal target.

Difficulty

The condition forbids ordinary cliques rather than imposing bounded complexity on every neighborhood. Estimates using maximum degree alone do not establish the stated average-degree bound.

Formalization scope

The representation is Mathlib SimpleGraph on a finite type. Real.log is the natural logarithm and the hypothesis d≥2d\ge2d≥2 keeps its argument and the denominator positive. The empty graph is excluded by that hypothesis.

The shared definitions are supplied by CliqueFreeLog. The target uses Lean 4.33.1 with the pinned Mathlib environment. Contributions should preserve the existing statement and develop reusable supporting lemmas in its definitions.

Selected references

  • OpenAI, A logarithmic independence bound for clique-free graphs, preprint, September 2026. Manuscript.
  • OpenAI, accompanying Lean statement. Pinned source.
25 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: marwahaha

Paired states and Hamiltonian cycles in cubic bipartite planar graphsOpen Problem

Motivation

A spanning subgraph in which every vertex has degree two may split into several cycles. Barnette’s question asks when planarity, bipartiteness, and connectivity force those vertices onto one cycle. The source is OpenAI's September 2026 manuscript.

Setting

A cubic graph has degree three at every vertex. A bipartite graph has two vertex classes with every edge crossing between them. A Hamiltonian cycle visits all vertices exactly once before returning to its starting point. Planarity is represented by a drawing using injective continuous edge arcs in the real plane.

Formalization targets

G finite, simple, cubic, bipartite, planar, and 3-vertex-connected ⟹ G has a Hamiltonian cycle.G\text{ finite, simple, cubic, bipartite, planar, and 3-vertex-connected}\ \Longrightarrow\ G\text{ has a Hamiltonian cycle}.G finite, simple, cubic, bipartite, planar, and 3-vertex-connected ⟹ G has a Hamiltonian cycle.

Every finite simple graph with decidable vertex equality and adjacency has a Hamiltonian cycle if each vertex has degree three, the vertices can be partitioned into two classes with every edge joining different classes, the graph is planar, and it is three-vertex-connected. Here planarity means that vertices can be placed at distinct points of the real plane and edges drawn as injective continuous arcs whose interiors contain no vertices and whose interiors are disjoint for distinct undirected edges; reversing an edge reverses its arc. Three-vertex-connectivity means that the graph has at least four vertices and remains connected after deletion of any set of at most two vertices. The conclusion is the existence of one closed cycle visiting every vertex exactly once before returning to its starting vertex.

The goal is OAI.Barnette.main.

Significance

A proof would turn the source’s Hamiltonicity claim into a result available for arbitrary finite graphs satisfying these exact hypotheses. It would also provide a reusable bridge between geometric embeddings and combinatorial walks. The manuscript presents a mathematical argument; the attached goal currently has status Open, so the requested contribution is a Lean proof of the displayed formal target.

Difficulty

Degree constraints do not guarantee that a spanning cycle cover is connected. The target requires a single cycle and must handle unrestricted face sizes, so verifying local degree conditions is insufficient.

Formalization scope

The published goal unfolds MainStatement. Three-vertex-connectivity requires at least four vertices and connectedness after deletion of any set of at most two vertices. The Hamiltonian object is a closed graph walk satisfying the cycle and vertex-coverage conditions. Stronger edge-avoidance conclusions from the paper are not attached as targets.

The shared definitions are supplied by BarnetteHamiltonian. The target uses Lean 4.33.1 with the pinned Mathlib environment. Contributions should preserve the existing statement and develop reusable supporting lemmas in its definitions.

Selected references

  • OpenAI, Paired states and Hamiltonian cycles in cubic bipartite planar graphs, preprint, September 2026. Manuscript.
  • OpenAI, accompanying Lean statement. Pinned source.
53 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: marwahaha

A counterexample to Sidorenko's conjectureOpen Problem

Motivation

Graph homomorphism inequalities compare the frequency of a fixed pattern with what edge density predicts. The mission asks for one finite host graph violating the proposed bound for a specified bipartite pattern. The source manuscript presents the surrounding research claim.

Setting

The pattern H is the incidence graph between thirteen points and twenty-two explicitly listed triples, giving thirty-five vertices and sixty-six edges. Homomorphism density counts all edge-preserving maps, including noninjective maps, divided by the number of all vertex maps.

Formalization target

Construct a positive-size finite simple host G with at least one edge such that

t(H,G)<p(G)66,p(G)=2∣E(G)∣∣V(G)∣2.t(H,G)<p(G)^{66},\qquad p(G)=\frac{2|E(G)|}{|V(G)|^2}.t(H,G)<p(G)66,p(G)=∣V(G)∣22∣E(G)∣​.

The selected formal target is OAI.SidorenkoCounterexample.main.

Significance and status

The counterexample uses the exact supplied pattern rather than an unspecified bipartite graph. A finite host makes the density comparison an ordinary finite counting statement. The published target is currently Open; the manuscript's argument and a checked Lean proof are separate deliverables.

Difficulty

The required strict inequality concerns all homomorphisms, not only embeddings. The host must have an edge so the assertion cannot be replaced by a zero-density degeneracy.

Formalization scope

H is built from the published table of twenty-two triples. Hosts use Fin size, and both graph counts and real density normalizations are fixed. The manuscript’s forcing-conjecture consequence is not an attached formal target.

Selected references

  • OpenAI, A counterexample to Sidorenko's conjecture, preprint, 2026. Pinned manuscript.
  • OpenAI, accompanying formal statements, commit adc7f1241b42. Goal source.
5 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: marwahaha

A linear list-coloring bound in terms of the Hadwiger numberOpen Problem

Motivation

List coloring permits each vertex to have different allowable colors. The mission asks whether the size of the largest clique minor still controls the required list size up to a universal linear factor. The source manuscript presents the surrounding research claim.

Setting

A clique minor of order t is represented by disjoint connected vertex sets, with an edge between every pair. The Hadwiger number h(G) is the largest possible t. The list chromatic number is the least k for which every assignment of at least k colors per vertex admits a proper selection.

Formalization target

Prove one integer C≥1 works for every finite nonempty simple graph:

χlist(G)≤C h(G).\chi_{list}(G)\le C\,h(G).χlist​(G)≤Ch(G).

The constant is chosen before the graph, the vertex type and all color lists.

The selected formal target is OAI.LinearListHadwiger.main_theorem.

Significance and status

The assertion controls list coloring rather than only a common-color-set chromatic number. It is a linear-factor statement, without coefficient-one optimality. The published target is currently Open; the manuscript's argument and a checked Lean proof are separate deliverables.

Difficulty

A coloring from one chosen palette need not respect independently assigned lists. The proof must cover arbitrary color types and list assignments.

Formalization scope

The definitions use finite graph vertex types, arbitrary finite color sets, induced-subgraph connectivity and natural-number infima/suprema. Nonemptiness of the vertex type is explicit. Only the main inequality and its definition block are attached.

Selected references

  • OpenAI, A linear list-coloring bound in terms of the Hadwiger number, preprint, 2026. Pinned manuscript.
  • OpenAI, accompanying formal statements, commit adc7f1241b42. Goal source.
9 thms2 active usersReviewed
🏆Completed
Number Theory·Captain: marwahaha

Bounded-Step Walks on Gaussian PrimesOpen Problem

Motivation

A planar prime walk can go around a gap that would block a walk on a line. The Gaussian moat problem therefore asks for an obstruction to every possible bounded-step route through distinct Gaussian primes. The source manuscript presents the surrounding research claim.

Setting

A Gaussian prime is an irreducible element of ℤ[i]. For D ∈ ℝ, two distinct primes are adjacent when their Euclidean distance is at most D. A component contains every prime reachable by a finite path.

Formalization target

Prove the conjunction

no injective infinite D-step prime walk exists,∀D  ∃B  ∀p, ∣component⁡D(p)∣≤B.\text{no injective infinite D-step prime walk exists},\qquad\forall D\;\exists B\;\forall p,\ |\operatorname{component}_D(p)|\le B.no injective infinite D-step prime walk exists,∀D∃B∀p, ∣componentD​(p)∣≤B.

The second clause also bounds the length of every finite injective D-step prime chain by B and asserts component finiteness.

The selected formal target is OAI.GaussianMoat.fullMain.

Significance and status

The same B must work for every starting prime at a fixed step bound. The theorem makes no numerical claim about this bound. The published target is currently Open; the manuscript's argument and a checked Lean proof are separate deliverables.

Difficulty

Showing that each component is finite does not give a uniform bound across starting points. Every route in two dimensions and primes on the axes remain within the quantification.

Formalization scope

The formal graph uses Mathlib Gaussian integers and complex Euclidean distance. Step bounds range over all real numbers. MainEndpoint and UniformEndpoint are both mandatory clauses of fullMain. Only that published target and its definitions are supplied.

Selected references

  • OpenAI, Bounded-Step Walks on Gaussian Primes, preprint, 2026. Pinned manuscript.
  • OpenAI, accompanying formal statements, commit adc7f1241b42. Goal source.
2 thms2 active usersReviewed
🏆Completed
Combinatorics·Captain: wurtle

A linear cycle-and-edge decomposition of every graphResearch Paper

Motivation

How efficiently can the edges of a graph be split into cycles? Some single edges must be allowed, since a forest has no cycles at all, and a tree on nnn vertices already needs n−1n-1n−1 single-edge parts. The Erdős–Gallai cycle decomposition conjecture asserts that this is essentially the only obstruction: every graph on nnn vertices can have its edge set partitioned into O(n)O(n)O(n) simple cycles and single edges. Equivalently, every graph with all degrees even partitions into O(n)O(n)O(n) cycles. The problem is a basic question about edge decompositions, and the partial bounds O(nlog⁡n)O(n\log n)O(nlogn), O(nlog⁡log⁡n)O(n\log\log n)O(nloglogn) and O(nlog⁡∗n)O(n\log^*n)O(nlog∗n) trace several decades of progress on it.

This mission asks for a formal proof of the linear bound, as stated in an OpenAI preprint dated September 24, 2026 (source). The preprint has not been peer reviewed and its claims have not been independently verified; on the platform the Lean goal is open.

Background

  • 1966 — Erdős, Goodman and Pósa record the problem and an O(nlog⁡n)O(n\log n)O(nlogn) bound, distinguishing covers from edge-disjoint partitions (Canad. J. Math. 1966).
  • 1968 — Lovász: every graph on nnn vertices partitions into at most ⌊n/2⌋\lfloor n/2\rfloor⌊n/2⌋ paths and cycles (On covering of graphs, 1968).
  • 1985 — Pyber shows that n−1n-1n−1 cycles and edges suffice to cover (not partition) every graph (Combinatorica 1985).
  • 2014 — Conlon, Fox and Sudakov prove O(nlog⁡log⁡n)O(n\log\log n)O(nloglogn) in general and linear bounds for random graphs and graphs of linear minimum degree (RSA 2014).
  • 2015–2021 — Asymptotically sharp counts for random graphs (Korándi–Krivelevich–Sudakov, CPC 2015); (3/2+δ)n(3/2+\delta)n(3/2+δ)n for dense graphs (Girão–Granet–Kühn–Osthus, J. LMS 2021).
  • 2024 — Bucić and Montgomery prove O(nlog⁡∗n)O(n\log^*n)O(nlog∗n) (Adv. Math. 2024).
  • 2025 — Akbari, Aloni, Beikmohammadi and Clow prove n−1n-1n−1 for maximum degree at most four (arXiv:2509.01901).
  • September 2026 — The OpenAI preprint claims the linear bound (Theorem 1.1, p. 1).

Setting

A finite simple graph GGG has vertex set of size nnn and edge set E(G)E(G)E(G) of unordered pairs of distinct vertices. A cycle is a closed walk of length at least 333 with no repeated vertices except the endpoints; its edge set is the set of edges it traverses. A cycle-and-edge decomposition of GGG into kkk parts is a family E1,…,EkE_1,\dots,E_kE1​,…,Ek​ of pairwise disjoint subsets of E(G)E(G)E(G) with union E(G)E(G)E(G), where each EiE_iEi​ is either the edge set of a cycle of GGG or a single edge of GGG. Parts may share vertices; isolated vertices need not be covered; for an edgeless graph k=0k=0k=0 is allowed.

Formalization targets

Goal: Theorem 1.1 (p. 1)

There is an absolute constant C>0C>0C>0 such that for every nnn and every simple graph GGG on nnn vertices,

G has a cycle-and-edge decomposition into at most Cn parts.G\ \text{has a cycle-and-edge decomposition into at most } Cn \text{ parts}.G has a cycle-and-edge decomposition into at most Cn parts.

The value of CCC is left unspecified, which keeps the goal stable under future improvements of the constant (the best possible constant is at least 3/23/23/2, by complete bipartite examples).

Significance

The result itself. The order nnn is optimal (trees). The theorem implies that every Eulerian graph (all degrees even) partitions into O(n)O(n)O(n) cycles (Corollary 1.2, p. 2), and via known equivalences it gives the sharp Δ(G)/2+δn\Delta(G)/2+\delta nΔ(G)/2+δn cycle bound for Eulerian graphs satisfying a large-cut condition. It removes the last iterated-logarithm loss from the Bucić–Montgomery bound.

Formalizing it. The statement uses only Mathlib's SimpleGraph.Walk.IsCycle and edge sets, so it is a clean combinatorial target. A complete development needs Lovász's path-cycle decomposition theorem, the Aharoni–Haxell hypergraph matching theorem, expansion and routing lemmas, and binomial tail bounds; all are reusable in extremal graph theory. No machine-checked proof of any o(nlog⁡n)o(n\log n)o(nlogn) bound is known.

Difficulty

Covering is easy (Pyber's n−1n-1n−1 bound), but a cover does not yield a partition: deleting repeated edges can break cycles into many paths. For partitions, the standard strategy decomposes the graph into expanding pieces at a sequence of scales, and paying a cost proportional to the full order at every scale accumulates a factor from the number of scales — this is where the log⁡log⁡n\log\log nloglogn and log⁡∗n\log^*nlog∗n losses come from. The preprint (Introduction, pp. 2–3) charges the work on a prefix of scales to a single expanding layer and saves a fixed fraction of vertices by identifying pairs, then lifts quotient cycles back with reserved edge-disjoint routes (Lemma 3.1, p. 5); keeping a separate assignment for each original edge throughout is what makes the induction close.

Formalization scope

  • Graphs are SimpleGraph (Fin n) for arbitrary n (including 0 and edgeless graphs).
  • CycleOrSingleEdge G s: either sss is p.edgeSet for some p : G.Walk v v with p.IsCycle (Mathlib cycles have length at least 333), or s={e}s=\{e\}s={e} for an edge e∈Ge\in Ge∈G.
  • EdgeDecomposition G k: parts indexed by Fin k, each a cycle or single edge, pairwise disjoint, with union G.edgeSet.
  • MainStatement: ∃C>0, ∀n ∀G, ∃k,\exists C>0,\ \forall n\ \forall G,\ \exists k,∃C>0, ∀n ∀G, ∃k, EdgeDecomposition G k ∧ k≤Cn\wedge\ k\le Cn∧ k≤Cn.
  • Trivializations are ruled out: parts must be nonempty edge sets of the required shape and must exactly partition the edge set, so the count kkk is a genuine count.

Selected references

  • OpenAI, A linear cycle-and-edge decomposition of every graph, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/A-linear-cycle-and-edge-decomposition-of-every-graph-September-24-2026/main.pdf
  • P. Erdős, A. W. Goodman, L. Pósa, The representation of a graph by set intersections, Canad. J. Math. 18 (1966). https://doi.org/10.4153/CJM-1966-014-3
  • L. Lovász, On covering of graphs, Theory of Graphs (Proc. Colloq. Tihany, 1966), 1968.
  • L. Pyber, An Erdős–Gallai conjecture, Combinatorica 5 (1985). https://doi.org/10.1007/BF02579444
  • D. Conlon, J. Fox, B. Sudakov, Cycle packing, Random Structures Algorithms (2014). https://doi.org/10.1002/rsa.20574
  • D. Korándi, M. Krivelevich, B. Sudakov, Decomposing random graphs into few cycles and edges, Combin. Probab. Comput. (2015). https://doi.org/10.1017/S0963548314000844
  • A. Girão, B. Granet, D. Kühn, D. Osthus, Path and cycle decompositions of dense graphs, J. London Math. Soc. (2021). https://doi.org/10.1112/jlms.12455
  • M. Bucić, R. Montgomery, Towards the Erdős–Gallai cycle decomposition conjecture, Adv. Math. (2024). https://doi.org/10.1016/j.aim.2023.109434
  • R. Aharoni, P. Haxell, Hall's theorem for hypergraphs, J. Graph Theory 35 (2000). https://doi.org/10.1002/1097-0118(200010)35:2<83::AID-JGT2>3.0.CO;2-V
  • S. Akbari, J. Aloni, A. Beikmohammadi, A. Clow, Tight bounds for cycle-edge decompositions and covers, arXiv:2509.01901 (2025). https://arxiv.org/abs/2509.01901v2
8 thms2 active usersReviewed
🏆Completed
CombinatoricsProbability·Captain: wurtle

The second Kahn–Kalai conjectureResearch Paper

Motivation

When does the Erdős–Rényi random graph G(n,p)G(n,p)G(n,p) contain a copy of a prescribed graph HHH? A necessary condition is visible from first moments: if some subgraph F⊆HF\subseteq HF⊆H has expected copy count below a constant, then HHH is unlikely to appear. Kahn and Kalai conjectured in 2007 that this elementary obstruction determines the containment threshold up to a logarithmic factor, uniformly over all target graphs, including spanning ones such as perfect matchings, Hamilton cycles or cubes (CPC 2007). Their abstract conjecture for general increasing families (the "first" conjecture) was proved by Park and Pham; the graph-specific "second" conjecture compares the threshold with the smaller, purely subgraph-count-based graph expectation threshold, and is not implied by the abstract theorem.

This mission asks for a formal proof of the second Kahn–Kalai conjecture with explicit constants, as stated in an OpenAI preprint dated September 24, 2026 (source). The preprint has not been peer reviewed and its claims have not been independently verified; on the platform the Lean goal is open.

Background

  • 1960–1981 — Erdős and Rényi's evolution of random graphs (1960) and perfect matchings (1966); Bollobás determines thresholds for fixed subgraphs (Math. Proc. Camb. Phil. Soc. 1981).
  • 2007 — Kahn and Kalai state both conjectures (CPC 2007).
  • 2010 — Talagrand introduces fractional expectation thresholds (STOC 2010).
  • 2021 — Frankston, Kahn, Narayanan and Park prove pc=O(qflog⁡ℓ)p_c=O(q_f\log\ell)pc​=O(qf​logℓ) (Annals 2021), building on the sunflower methods of Alweiss–Lovett–Wu–Zhang (Annals 2021).
  • 2022–2025 — Mossel, Niles-Weed, Sun and Zadik obtain a single-log bound for a modified graph threshold (arXiv:2209.03326) and a Bayesian proof of the spread lemma (RSA 2025).
  • 2024 — Park and Pham prove the abstract Kahn–Kalai conjecture (JAMS 2024).
  • 2025 — Dubroff, Kahn and Park prove pc=O(pElog⁡3n)p_c=O(p_E\log^3n)pc​=O(pE​log3n) and special cases (arXiv:2508.14269, arXiv:2510.02005).
  • 2026 — Tran reduces the loss to O(log⁡2(2e(H)))O(\log^2(2e(H)))O(log2(2e(H))) and proves the conjecture for trees and other classes (arXiv:2609.20546); the OpenAI preprint claims the general case (Theorem 1.1, p. 1).

Setting

For finite simple graphs G,FG,FG,F let N(G,F)N(G,F)N(G,F) be the number of subgraphs of GGG isomorphic to FFF (ordinary, non-induced copies). G(n,p)G(n,p)G(n,p) is the random graph on [n][n][n] with each of the (n2)\binom n2(2n​) edges present independently with probability ppp. For a graph HHH with at least one edge and at most nnn vertices,

pc(n,H)=inf⁡{p∈[0,1]: Pr⁡[N(G(n,p),H)>0]≥12},p_c(n,H)=\inf\{p\in[0,1]:\ \Pr[N(G(n,p),H)>0]\ge\tfrac12\},pc​(n,H)=inf{p∈[0,1]: Pr[N(G(n,p),H)>0]≥21​}, pE(n,H)=inf⁡{p∈[0,1]: E N(G(n,p),F)≥12 for every subgraph F⊆H}.p_E(n,H)=\inf\{p\in[0,1]:\ \mathbb E\,N(G(n,p),F)\ge\tfrac12\ \text{for every subgraph }F\subseteq H\}.pE​(n,H)=inf{p∈[0,1]: EN(G(n,p),F)≥21​ for every subgraph F⊆H}.

The first is the containment threshold, the second the graph expectation threshold; always pE≤pcp_E\le p_cpE​≤pc​. Write h=e(H)h=e(H)h=e(H) for the number of edges.

Formalization targets

Goal: Theorem 1.1 (p. 1)

For every n≥2n\ge2n≥2 and every finite simple graph HHH with h≥1h\ge1h≥1 edges and at most nnn vertices,

pc(n,H) ≤ min⁡{1, 2048 e50 pE(n,H) (1+log⁡2h)},pc(n,H) ≤ 6144 e50 pE(n,H)log⁡2n.p_c(n,H)\ \le\ \min\bigl\{1,\ 2048\,e^{50}\,p_E(n,H)\,(1+\log_2 h)\bigr\},\qquad p_c(n,H)\ \le\ 6144\,e^{50}\,p_E(n,H)\log_2 n .pc​(n,H) ≤ min{1, 2048e50pE​(n,H)(1+log2​h)},pc​(n,H) ≤ 6144e50pE​(n,H)log2​n.

The constants are explicit and universal; HHH may vary with nnn.

Significance

The result itself. The bound says that, up to a universal constant times log⁡h\log hlogh, the only obstruction to containing HHH is a subgraph with too few expected copies. The logarithm cannot be removed: for a perfect matching pE=Θ(1/n)p_E=\Theta(1/n)pE​=Θ(1/n) while pc=Θ(log⁡n/n)p_c=\Theta(\log n/n)pc​=Θ(logn/n), because isolated vertices must disappear. This settles the problem left by the abstract Park–Pham theorem, whose bound is in terms of the possibly larger integral expectation threshold of the family of graphs containing HHH.

Formalizing it. Both thresholds are finite, explicit optimization problems over Bernoulli product measures on the edges of KnK_nKn​, so the statement is fully elementary. Mathlib has SimpleGraph.copyCount but no random-graph or threshold theory; the spread-family resampling lemma and the tree covering theorem would be reusable for other threshold results (Park–Pham, Frankston–Kahn–Narayanan–Park). No machine-checked proof of any Kahn–Kalai-type theorem is known.

Difficulty

The abstract theorem compares pcp_cpc​ with the integral expectation threshold of the family FH\mathcal F_HFH​ of edge sets containing a copy of HHH, and covers of FH\mathcal F_HFH​ may be much cheaper than the copies of a single subgraph, so it does not give a bound in terms of pEp_EpE​. Previous graph-specific approaches lost extra logarithms because each conditional extension step paid the sampling cost again. The preprint (Section 1.2, pp. 3) builds a tree of nested subgraph extensions with geometrically growing, pairwise disjoint labels along each path, and must reuse one random set at every level of the tree while keeping the spread property under conditioning (Theorem 3.1, p. 6).

Formalization scope

  • G(n,p)G(n,p)G(n,p) is encoded as a Bernoulli product weight on Finset (Edge n), where Edge n is the type of non-diagonal Sym2 (Fin n); expectations and probabilities are finite sums, defined for every real ppp but used only on [0,1][0,1][0,1].
  • expectedCopies and containmentProbability use Mathlib's copyCount of the graph built from the edge set; both are ordinary (non-induced) copies.
  • criticalThreshold and expectationThreshold are sInf over p∈[0,1]p\in[0,1]p∈[0,1]; the sets are nonempty (they contain p=1p=1p=1 when HHH has at most nnn vertices), so the infima are the true thresholds. The expectation constraint ranges over all H.Subgraphs, including subgraphs with isolated vertices, whose constraints are automatically satisfied.
  • HHH is a SimpleGraph on an arbitrary finite type V with Fintype.card V ≤ n, 2 ≤ n and at least one edge; logTwo is log⁡x/log⁡2\log x/\log 2logx/log2.
  • The constants 2048e502048e^{50}2048e50 and 6144e506144e^{50}6144e50 are hard-coded, as in the source.

Selected references

  • OpenAI, The second Kahn–Kalai conjecture, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-second-Kahn-Kalai-conjecture-September-24-2026/paper.pdf
  • J. Kahn, G. Kalai, Thresholds and expectation thresholds, Combin. Probab. Comput. 16 (2007), 495–502. https://doi.org/10.1017/S0963548307008474
  • J. Park, H. T. Pham, A proof of the Kahn–Kalai conjecture, J. Amer. Math. Soc. 37 (2024), 235–243. https://doi.org/10.1090/jams/1028
  • K. Frankston, J. Kahn, B. Narayanan, J. Park, Thresholds versus fractional expectation-thresholds, Ann. of Math. 194 (2021), 475–495. https://doi.org/10.4007/annals.2021.194.2.2
  • R. Alweiss, S. Lovett, K. Wu, J. Zhang, Improved bounds for the sunflower lemma, Ann. of Math. 194 (2021), 795–815. https://doi.org/10.4007/annals.2021.194.3.5
  • M. Talagrand, Are many small sets explicitly small?, STOC 2010. https://doi.org/10.1145/1806689.1806693
  • E. Mossel, J. Niles-Weed, N. Sun, I. Zadik, On the second Kahn–Kalai conjecture, arXiv:2209.03326 (2022). https://arxiv.org/abs/2209.03326v1
  • E. Mossel, J. Niles-Weed, N. Sun, I. Zadik, A Bayesian proof of the spread lemma, Random Structures Algorithms 66 (2025), e70008. https://doi.org/10.1002/rsa.70008
  • Q. Dubroff, J. Kahn, J. Park, On the "second" Kahn–Kalai conjecture, arXiv:2508.14269 (2025). https://arxiv.org/abs/2508.14269v1
  • T. Tran, Fractional expectation thresholds and the "second" Kahn–Kalai conjecture, arXiv:2609.20546 (2026). https://arxiv.org/abs/2609.20546v2
  • B. Bollobás, Threshold functions for small subgraphs, Math. Proc. Cambridge Philos. Soc. 90 (1981), 197–206. https://doi.org/10.1017/S0305004100058655
47 thms2 active usersReviewed
🏆Completed
CombinatoricsTheoretical Computer Science·Captain: wurtle

The strong thin tree conjectureResearch Paper

Motivation

A spanning tree keeps a graph connected with as few edges as possible. A spanning tree is thin if, in addition, it uses only a small fraction of the edges of every cut. Thin trees were introduced as a tool for the asymmetric traveling salesman problem (ATSP): Asadpour, Goemans, Mądry, Oveis Gharan and Saberi showed that thin trees, together with cost control, can be augmented cheaply into tours, which gave an O(log⁡n/log⁡log⁡n)O(\log n/\log\log n)O(logn/loglogn)-approximation (doi:10.1287/opre.2017.1603). Goddyn's thin tree conjecture asks whether high edge connectivity alone forces an ε\varepsilonε-thin spanning tree, independently of the number of vertices; the strong form asks for thinness C/kC/kC/k in a kkk-edge-connected graph.

Background. Goddyn recorded the conjecture in a 2004 problem list. The Nash–Williams–Tutte theorem (1961) gives ⌊k/2⌋\lfloor k/2\rfloor⌊k/2⌋ edge-disjoint spanning trees in a kkk-edge-connected graph (doi:10.1112/jlms/s1-36.1.445, doi:10.1112/jlms/s1-36.1.221), so each cut is used about 2/k2/k2/k times on average, but no single tree need be good for every cut. Oveis Gharan and Saberi proved the bound for planar and bounded-genus graphs (2011, doi:10.1137/1.9781611973082.75). Anari and Oveis Gharan obtained thinness poly(log⁡log⁡n)/k\mathrm{poly}(\log\log n)/kpoly(loglogn)/k in general (2015, arXiv:1411.4613), using the Marcus–Spielman–Srivastava interlacing-families method (doi:10.4007/annals.2015.182.1.8). Klein and Olver handled any prescribed laminar family of cuts (2023, doi:10.1109/FOCS57990.2023.00011), and Klein, Olver and Yeoh controlled all near-minimum cuts (2026, doi:10.4230/LIPIcs.ICALP.2026.129). Constant-factor ATSP approximations were obtained by other means (Svensson–Tarnawski–Végh 2020, doi:10.1145/3424306; Traub–Vygen 2022, doi:10.1137/20M1339313), but the thin tree question stayed unresolved.

The source of this mission, an OpenAI preprint dated September 23, 2026, claims the strong thin tree conjecture.

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite undirected multigraph without loops; parallel edges are counted separately. For a nonempty proper subset ∅≠S⊊V\varnothing\ne S\subsetneq V∅=S⊊V, let δG(S)\delta_G(S)δG​(S) be the set of edges with exactly one endpoint in SSS. GGG is kkk-edge-connected if ∣δG(S)∣≥k|\delta_G(S)|\ge k∣δG​(S)∣≥k for every such SSS. For an edge set T⊆ET\subseteq ET⊆E, δT(S)=δG(S)∩T\delta_T(S)=\delta_G(S)\cap TδT​(S)=δG​(S)∩T. A spanning tree TTT is α\alphaα-thin if ∣δT(S)∣≤α∣δG(S)∣|\delta_T(S)|\le\alpha|\delta_G(S)|∣δT​(S)∣≤α∣δG​(S)∣ for every nonempty proper SSS.

Formalization targets

Goal: Theorem 1.1 (strong thin trees)

There is a universal constant C>0C>0C>0 such that for every k≥1k\ge1k≥1, every finite loopless kkk-edge-connected multigraph with at least two vertices has a spanning tree TTT with

∣δT(S)∣≤Ck ∣δG(S)∣(∅≠S⊊V).|\delta_T(S)|\le\frac{C}{k}\,|\delta_G(S)|\qquad(\varnothing\ne S\subsetneq V).∣δT​(S)∣≤kC​∣δG​(S)∣(∅=S⊊V).

Lean: OAI.StrongThinTree.strongThinTree, open on the platform.

Significance

The order 1/k1/k1/k is optimal (two vertices joined by kkk parallel edges). The theorem would settle both Goddyn's thin tree conjecture and its strong form, removing the poly(log⁡log⁡n)\mathrm{poly}(\log\log n)poly(loglogn) loss of Anari–Oveis Gharan. It shows that connectivity alone controls all cuts at once, and gives a structural route to ATSP integrality-gap bounds through thin trees. An algorithmic version (deterministic polynomial-time construction) is the subject of the companion mission in this family. The result is claimed in an OpenAI preprint that has not been peer reviewed; no machine-checked proof exists.

Difficulty

Averaging over a Nash–Williams–Tutte packing controls each cut on average but not all cuts simultaneously: there are exponentially many cuts and a union bound over them loses a factor depending on nnn. Spectral approaches control all cuts via a matrix inequality, but spectrally thin trees need not exist when some edges have large effective resistance, which is why previous spectral arguments lost log⁡log⁡n\log\log nloglogn factors. One must instead keep a large tree packing on low-resistance edges through a sequence of sparsifications whose relative losses stay bounded.

Formalization scope

  • MultiGraph n m: vertices Fin n, edges Fin m with endpoint maps and a looplessness proof; parallel edges are distinct edges.
  • cut G T S: the edges of T with exactly one endpoint in S. Connected G T: every nonempty proper S is crossed by T. SpanningTree G T: T connected and every proper deletion disconnects.
  • EdgeConnected G k: every nonempty proper cut has at least k edges.
  • MainStatement: ∃ C > 0, ∀ k ≥ 1, ∀ n ≥ 2, ∀ G k-edge-connected, ∃ T spanning tree, ∀ S nonempty proper, |cut T S| ≤ (C/k)|cut G S| (real arithmetic).

Needed infrastructure: graph Laplacians and effective resistance, the Nash–Williams–Tutte theorem, the Marcus–Spielman–Srivastava mixed-characteristic-polynomial bound, and a Brouwer-type fixed point theorem (the paper proves the finite-dimensional case via Sperner's lemma). Contributions formalizing Theorem 2.3 (Nash–Williams–Tutte), Theorem 3.2, Theorem 4.1, Proposition 6.2 or Theorem A.1 are welcome.

Selected references

  • OpenAI, The strong thin tree conjecture, preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-strong-thin-tree-conjecture-September-23-2026/paper.pdf
  • L. A. Goddyn, Some open problems I like, problem list, 2004.
  • A. Asadpour, M. X. Goemans, A. Mądry, S. Oveis Gharan, A. Saberi, An O(log n/log log n)-approximation algorithm for the asymmetric traveling salesman problem, Oper. Res., 2017. https://doi.org/10.1287/opre.2017.1603
  • N. Anari, S. Oveis Gharan, Effective-resistance-reducing flows, spectrally thin trees, and asymmetric TSP, FOCS, 2015. https://arxiv.org/abs/1411.4613
  • A. W. Marcus, D. A. Spielman, N. Srivastava, Interlacing families II: mixed characteristic polynomials and the Kadison–Singer problem, Ann. of Math., 2015. https://doi.org/10.4007/annals.2015.182.1.8
  • S. Oveis Gharan, A. Saberi, The asymmetric traveling salesman problem on graphs with bounded genus, SODA, 2011. https://doi.org/10.1137/1.9781611973082.75
  • N. Klein, N. Olver, Thin trees for laminar families, FOCS, 2023. https://doi.org/10.1109/FOCS57990.2023.00011
  • N. Klein, N. Olver, Z. S. Yeoh, Thin trees for near minimum cuts, ICALP, 2026. https://doi.org/10.4230/LIPIcs.ICALP.2026.129
  • C. St. J. A. Nash-Williams, Edge-disjoint spanning trees of finite graphs, J. London Math. Soc., 1961. https://doi.org/10.1112/jlms/s1-36.1.445
5 thms2 active usersReviewed
🏆Completed
Complexity TheoryTheoretical Computer Science·Captain: wurtle

Parity lifts and bounded-treewidth witnesses for Weisfeiler–Leman equivalenceResearch Paper

Motivation

The Weisfeiler–Leman (WL) method is the standard combinatorial heuristic for graph isomorphism: it repeatedly refines colors of kkk-tuples of vertices, and two graphs it cannot tell apart are called kkk-WL equivalent. It is equivalent to counting-logic indistinguishability and, by results of Dvořák (2010) and Dell, Grohe and Rattan (2018), to equality of homomorphism counts from graphs of treewidth below kkk (bags of size at most k+1k+1k+1). For fixed kkk, direct refinement decides kkk-WL equivalence in nO(k)n^{O(k)}nO(k) time; whether the growing exponent is necessary is a natural complexity question.

This preprint gives an exact reduction from a simple constraint problem (choosing compatible elements from finite domains) to kkk-WL equivalence of two explicit uncolored graphs, and uses it to derive nΩ(k)n^{\Omega(k)}nΩ(k) conditional time lower bounds.

Background

  • 1968. Weisfeiler and Leman introduce the refinement method for graph canonization.
  • 1992. Cai, Fürer and Immerman construct parity (CFI) graph pairs of bounded degree that need WL dimension linear in their order.
  • 1999. Grohe shows equivalence in finite-variable logics is complete for polynomial time.
  • 2001. Impagliazzo and Paturi formulate the Exponential Time Hypothesis; Impagliazzo, Paturi and Zane give sparsification.
  • 2010, 2018. Dvořák, and Dell–Grohe–Rattan, characterize kkk-WL equivalence by homomorphism counts from bounded-treewidth graphs.
  • 2019. Atserias, Mančinska, Roberson, Šámal, Severini and Varvitsiotis encode solutions of binary linear systems as graph vertices.
  • 2025. Grohe, Lichter, Neuen and Schweitzer prove Ω(nk/2)\Omega(n^{k/2})Ω(nk/2) round lower bounds for joint refinement; Lichter, Raßmann and Schweitzer conjecture that WL equivalence and identification do not admit no(k)n^{o(k)}no(k) time.
  • 2026. The OpenAI preprint Parity lifts and bounded-treewidth witnesses for Weisfeiler–Leman equivalence (dated September 25, 2026) proves the parity reduction below. It has not been peer reviewed; the Lean goal is open on this platform.

Setting

Fix k≥4k \ge 4k≥4. The atomic type of a tuple x∈V(X)k\mathbf{x} \in V(X)^kx∈V(X)k records all equalities and adjacencies among its entries. In joint refinement a round records the old color of x\mathbf xx and the multiset over y∈V(X)y \in V(X)y∈V(X) of the vectors (Cr(x[1←y]),…,Cr(x[k←y]))(C_r(\mathbf x[1\leftarrow y]), \dots, C_r(\mathbf x[k\leftarrow y]))(Cr​(x[1←y]),…,Cr​(x[k←y])). In separate refinement it records the old color and the kkk multisets { ⁣{Cr(x[i←y]):y} ⁣}\{\!\{C_r(\mathbf x[i\leftarrow y]) : y\}\!\}{{Cr​(x[i←y]):y}}. Two graphs are kkk-WL equivalent (≡k\equiv_k≡k​) when their tuple-color histograms agree at every round. The bag parameter is t=k+1t = k+1t=k+1 (joint) or t=kt = kt=k (separate).

A choice system on ttt domains consists of finite sets D1,…,DtD_1, \dots, D_tD1​,…,Dt​ and, for each pair i<ji < ji<j, a finite label set LijL_{ij}Lij​ with maps λiji:Di→Lij\lambda^i_{ij} : D_i \to L_{ij}λiji​:Di​→Lij​, λijj:Dj→Lij\lambda^j_{ij} : D_j \to L_{ij}λijj​:Dj​→Lij​. A successful choice is (d1,…,dt)∈∏iDi(d_1, \dots, d_t) \in \prod_i D_i(d1​,…,dt​)∈∏i​Di​ with λiji(di)=λijj(dj)\lambda^i_{ij}(d_i) = \lambda^j_{ij}(d_j)λiji​(di​)=λijj​(dj​) for all i<ji < ji<j. Empty domains are allowed.

In Lean (OAI.ParityWL), ChoiceSystem t packages the domains and labels, template t is the ttt-clique with helpers, baseGraph C replaces template types by domains and labels, and zeroLift C, starLift C are the two parity lifts. Colors are defined by recursion on rounds (jointColor, separateColor) with histograms compared as multisets.

Formalization targets

Goal: Theorem 1.2 (Parity reduction)

For k≥4k \ge 4k≥4, either convention and the corresponding ttt, for every choice system on ttt domains the two explicit uncolored simple graphs X0,X⋆X_0, X_\starX0​,X⋆​ have equal order and

X0≡kX⋆  ⟺  the choice system has no successful choice,X_0 \equiv_k X_\star \iff \text{the choice system has no successful choice},X0​≡k​X⋆​⟺the choice system has no successful choice,

both with common color names on separate graphs and with pure-tuple comparison in the disjoint union. If every domain and label set has size at most A≥1A \ge 1A≥1, the common order is at most AMtA M_tAMt​ with

Mt=t 2 t−2+3(t−12)+12(t3).M_t = t\, 2^{\,t-2+3\binom{t-1}{2}} + 12\binom{t}{3}.Mt​=t2t−2+3(2t−1​)+12(3t​).

When a successful choice exists, the histograms already differ after two joint rounds or one separate round.

Significance

The result. The reduction turns compatible-choice problems, which encode sparse satisfiability, into WL equivalence of uncolored graphs whose size grows only linearly in the domain size. Combined with sparsification it gives Theorem 1.3: under the positive-rate ETH there are c>0c > 0c>0, K≥4K \ge 4K≥4 such that no deterministic algorithm decides kkk-WL equivalence of nnn-vertex graphs in time O(nck)O(n^{ck})O(nck) for fixed k≥Kk \ge Kk≥K, even when only early-round histograms are compared. This addresses the expectation stated by Lichter, Raßmann and Schweitzer.

Formalizing it. No machine-checked proof exists; the source is an unrefereed preprint. A formalization would also give Lean definitions of both WL conventions, reusable by the companion missions in this family on identification and unconditional lower bounds. The ETH consequence (Theorem 1.3) is not part of the Lean goal.

Difficulty

The forward direction (a successful choice yields a distinguishing homomorphism test detected in few rounds) is constructive. The converse is the hard part: one must show that every bounded-treewidth homomorphism test that distinguishes the lifts yields a successful choice, even when bags do not induce cliques, projected images repeat, and some weights vanish. A naive argument that counts only type-preserving homomorphisms fails because non-type-preserving maps also contribute to hom⁡(F,X0)−hom⁡(F,X⋆)\hom(F, X_0) - \hom(F, X_\star)hom(F,X0​)−hom(F,X⋆​).

Formalization scope

  • Graphs are SimpleGraph on finite types; the lifts' vertex types are sigma types of base vertices with parity tags.
  • WL colors are nested Multiset-valued types defined by recursion on the round number; equivalence quantifies over all rounds rrr.
  • The order bound uses orderFactor t = t * 2^(t-2+3*(t-1).choose 2) + 12 * t.choose 3 with natural-number subtraction (harmless since t≥4t \ge 4t≥4).
  • Successful is an existential over dependent tuples; empty domains are allowed, in which case no successful choice exists and the lifts must be equivalent.

The goal is not trivial: the biconditional pins down equivalence exactly, and the detection clause fixes the round.

Needed infrastructure: tree decompositions and homomorphism counts, the Dvořák/Dell–Grohe–Rattan interface for both conventions, and F2\mathbb{F}_2F2​ linear algebra for the parity lifts. Contributions toward either direction of the equivalence are welcome.

Selected references

  • OpenAI, Parity lifts and bounded-treewidth witnesses for Weisfeiler–Leman equivalence, OpenAI Math Release preprint, September 25, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Parity-lifts-and-bounded-treewidth-witnesses-for-Weisfeiler-Leman-equivalence-September-25-2026/paper.pdf
  • B. Weisfeiler, A. Leman, The Reduction of a Graph to Canonical Form and the Algebra Which Appears Therein, 1968. https://www.iti.zcu.cz/wl2018/pdf/wl_paper_translation.pdf
  • J.-Y. Cai, M. Fürer, N. Immerman, An Optimal Lower Bound on the Number of Variables for Graph Identification, Combinatorica 12 (1992). https://doi.org/10.1007/BF01305232
  • Z. Dvořák, On Recognizing Graphs by Numbers of Homomorphisms, J. Graph Theory (2010). https://doi.org/10.1002/jgt.20461
  • H. Dell, M. Grohe, G. Rattan, Lovász Meets Weisfeiler and Leman, ICALP 2018. https://doi.org/10.4230/LIPIcs.ICALP.2018.40
  • A. Atserias et al., Quantum and Non-Signalling Graph Isomorphisms, J. Combin. Theory Ser. B (2019). https://doi.org/10.1016/j.jctb.2018.11.002
  • R. Impagliazzo, R. Paturi, On the Complexity of k-SAT, J. Comput. Syst. Sci. (2001). https://doi.org/10.1006/jcss.2000.1727
  • M. Grohe, M. Lichter, D. Neuen, P. Schweitzer, Compressing CFI Graphs and Lower Bounds for the Weisfeiler–Leman Refinements, J. ACM (2025). https://doi.org/10.1145/3727978
  • M. Lichter, S. Raßmann, P. Schweitzer, Computational Complexity of the Weisfeiler–Leman Dimension, CSL 2025. https://doi.org/10.4230/LIPIcs.CSL.2025.13
31 thms2 active usersReviewed
🏆Completed
Discrete GeometryTheoretical Computer Science·Captain: wurtle

L1 Embeddings of Graphs of Bounded TreewidthResearch Paper

Motivation: L1L_1L1​ embeddings, sparsest cut and the GNRS conjecture

A finite metric embeds into L1L_1L1​ exactly when it is a nonnegative combination of cut metrics δS(u,v)=∣1S(u)−1S(v)∣\delta_S(u,v)=|\mathbf 1_S(u)-\mathbf 1_S(v)|δS​(u,v)=∣1S​(u)−1S​(v)∣. This makes L1L_1L1​ distortion the right measure for the flow–cut gap: given edge capacities and pairwise demands, the ratio between the sparsest cut value ϕ∗\phi_*ϕ∗​ and the maximum concurrent flow λ∗\lambda_*λ∗​. Linial, London and Rabinovich (1995) and Aumann and Rabani (1998) developed this connection, and Gupta, Newman, Rabinovich and Sinclair (2004) proved that the worst flow–cut gap on a fixed graph equals the worst L1L_1L1​ distortion of its weighted shortest-path metrics. They conjectured that every proper minor-closed family embeds into L1L_1L1​ with uniformly bounded distortion (the GNRS conjecture). Graphs of bounded treewidth form one of its basic cases.

Timeline

  • 1986 — Robertson and Seymour: excluding a fixed planar graph as a minor bounds the treewidth (JCTB 1986).
  • 2004 — Gupta, Newman, Rabinovich and Sinclair prove the treewidth-two (series-parallel) case, show that some series-parallel metrics need expected tree distortion Ω(log⁡n)\Omega(\log n)Ω(logn), and state the GNRS conjecture (Combinatorica 2004).
  • 2006 — Chekuri, Gupta, Newman, Rabinovich and Sinclair: bounded distortion for kkk-outerplanar graphs (SIDMA 2006).
  • 2008–2010 — Chakrabarti, Jaffe, Lee and Vincent obtain the optimal bound 222 for treewidth two (FOCS 2008), matched by Lee and Raghavendra (DCG 2010).
  • 2010 — Chlamtáč, Krauthgamer and Raghavendra: bounded distortion for bounded treewidth under a local-consistency assumption on small vertex sets (APPROX/RANDOM 2010).
  • 2013 — Lee and Sidiropoulos settle bounded pathwidth via stochastic tree embeddings (Combinatorica 2013); Lee and Poore treat 2-sums of a fixed finite family (SoCG 2013).
  • 2022 — Abraham, Filtser, Gupta and Neiman: O(p)O(\sqrt p)O(p​) distortion for pathwidth ppp (SICOMP 2022).
  • 2025 — Filtser et al.: optimal padded decompositions for bounded treewidth and an O(log⁡(2+t)log⁡n)O(\sqrt{\log(2+t)\log n})O(log(2+t)logn​) bound (TheoretiCS 2025).
  • 2026 — An OpenAI preprint, L1L_1L1​ Embeddings of Graphs of Bounded Treewidth (OpenAI Math Release, September 23, 2026), claims a uniform bound for every fixed treewidth. The preprint has not been peer reviewed, and its main theorem is not formally verified.

Setting

Let G=(V,E)G=(V,E)G=(V,E) be a finite connected graph with edge lengths ℓ:E→(0,∞)\ell:E\to(0,\infty)ℓ:E→(0,∞). The shortest-path metric dG,ℓ(u,v)d_{G,\ell}(u,v)dG,ℓ​(u,v) is the least total length of a path from uuu to vvv.

A tree decomposition of GGG is a finite tree TTT with bags Bt⊆VB_t\subseteq VBt​⊆V (t∈V(T)t\in V(T)t∈V(T)) such that the bags cover VVV, each edge has both endpoints in some bag, and for each vertex vvv the nodes whose bag contains vvv induce a connected subtree. In this mission, as in the source, kkk bounds the bag cardinality ∣Bt∣≤k|B_t|\le k∣Bt​∣≤k, so the treewidth bound is k−1k-1k−1.

A map F:V→ℓ1mF:V\to\ell_1^mF:V→ℓ1m​ has distortion at most CCC if dG,ℓ(u,v)≤∥F(u)−F(v)∥1≤C dG,ℓ(u,v)d_{G,\ell}(u,v)\le\|F(u)-F(v)\|_1\le C\,d_{G,\ell}(u,v)dG,ℓ​(u,v)≤∥F(u)−F(v)∥1​≤CdG,ℓ​(u,v).

Formalization targets

Goal: Theorem 1.1 (bounded-treewidth graph metrics in L1L_1L1​)

For every integer k≥2k\ge2k≥2 there is a constant Ck≥1C_k\ge1Ck​≥1 such that for every nonempty finite connected graph GGG with a tree decomposition of bag size at most kkk and every ℓ:E→(0,∞)\ell:E\to(0,\infty)ℓ:E→(0,∞), there are mmm and F:V→ℓ1mF:V\to\ell_1^mF:V→ℓ1m​ with

dG,ℓ(u,v)  ≤  ∥F(u)−F(v)∥1  ≤  Ck dG,ℓ(u,v)(u,v∈V).d_{G,\ell}(u,v)\;\le\;\|F(u)-F(v)\|_1\;\le\;C_k\,d_{G,\ell}(u,v)\qquad(u,v\in V).dG,ℓ​(u,v)≤∥F(u)−F(v)∥1​≤Ck​dG,ℓ​(u,v)(u,v∈V).

CkC_kCk​ is independent of the number of vertices, the decomposition and all ratios of edge lengths; the goal asserts only its existence. The goal statement is published on the platform with status Open: no machine-checked proof exists yet.

Significance

The result itself. The theorem settles the bounded-treewidth case of the GNRS conjecture. Through the GNRS equivalence, it gives λ∗≤ϕ∗≤Ckλ∗\lambda_*\le\phi_*\le C_k\lambda_*λ∗​≤ϕ∗​≤Ck​λ∗​ for every multicommodity instance on graphs of treewidth less than kkk. By the Robertson–Seymour excluded-planar-minor theorem it also covers every family excluding a fixed planar minor. Combined with known approximation algorithms for the minimum L1L_1L1​ distortion (Gupta–Talwar–Witmer 2013; Cohen-Addad–Mömke–Verdugo 2024), it yields fixed-parameter-time embeddings of distortion (2+ε)Ck(2+\varepsilon)C_k(2+ε)Ck​. With the planar companion preprint it covers the almost-embeddable families of Lee and Sidiropoulos.

Formalizing it. The proof constructs cut measures directly on a tree decomposition and combines them with padded random partitions from the literature. A formal proof would certify that the constants depend only on kkk, the point at which earlier approaches lost uniformity.

Difficulty

The successful pathwidth approach embeds the graph stochastically into dominating trees, but that route cannot work here: some series-parallel metrics require expected tree distortion Ω(log⁡n)\Omega(\log n)Ω(logn) (GNRS 2004, Theorem 5.6). Local-consistency approaches assume that small vertex sets admit compatible isometric cut representations, which a general bounded-treewidth metric need not provide. Multiscale constructions from padded partitions give good separation at each scale, but summing their contributions loses a factor growing with the number of scales. The construction must handle arbitrary branching of the decomposition tree and arbitrary edge-length ratios without such accumulation.

Formalization scope

  • The vertex type V is any nonempty Fintype with decidable equality; G : SimpleGraph V is connected.
  • HasTreeDecomposition G k: a tree T on Fin n, bags Fin n → Finset V covering all vertices and all edges, with the bags containing each vertex inducing a connected subgraph of T, and all bags of cardinality at most k.
  • Edge lengths are ℓ : G.edgeSet → ℝ, strictly positive. shortestPathDistance is the sInf of lengths of paths (IsPath) from uuu to vvv; this set is finite and nonempty for a connected finite graph, so the infimum is the true distance.
  • The target is the finite-dimensional space PiLp 1 (Fin m → ℝ) with mmm chosen per graph; this is the source's remark (p. 2) that the target can be taken to be finite-dimensional ℓ1\ell_1ℓ1​, and it is at least as strong as an L1(μ)L_1(\mu)L1​(μ) target.
  • The constant C depends only on k and is quantified before the graph.
  • Infrastructure needed: tree decompositions, cut-measure representations of ℓ1\ell_1ℓ1​ metrics, Markov flows on layered graphs, padded decompositions. Tree-decomposition and cut-metric infrastructure is reusable.

Selected references

  • N. Robertson, P. D. Seymour, Graph minors. V. Excluding a planar graph, J. Combin. Theory Ser. B 41 (1986), 92–114. https://doi.org/10.1016/0095-8956(86)90030-4
  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15 (1995), 215–245. https://doi.org/10.1007/BF01200757
  • Y. Aumann, Y. Rabani, An O(log⁡k)O(\log k)O(logk) approximate min-cut max-flow theorem and approximation algorithm, SIAM J. Comput. 27 (1998), 291–301. https://doi.org/10.1137/S0097539794285983
  • A. Gupta, I. Newman, Y. Rabinovich, A. Sinclair, Cuts, trees and ℓ1\ell_1ℓ1​-embeddings of graphs, Combinatorica 24 (2004), 233–269. https://doi.org/10.1007/s00493-004-0015-x
  • A. Chakrabarti, A. Jaffe, J. R. Lee, J. Vincent, Embeddings of topological graphs: lossy invariants, linearization, and 2-sums, FOCS 2008, 761–770. https://doi.org/10.1109/FOCS.2008.79
  • J. R. Lee, P. Raghavendra, Coarse differentiation and multi-flows in planar graphs, Discrete Comput. Geom. 43 (2010), 346–362. https://doi.org/10.1007/s00454-009-9172-4
  • E. Chlamtáč, R. Krauthgamer, P. Raghavendra, Approximating sparsest cut in graphs of bounded treewidth, APPROX/RANDOM 2010, LNCS 6302, 124–137. https://doi.org/10.1007/978-3-642-15369-3_10
  • J. R. Lee, A. Sidiropoulos, Pathwidth, trees, and random embeddings, Combinatorica 33 (2013), 349–374. https://doi.org/10.1007/s00493-013-2685-8
  • I. Abraham, A. Filtser, A. Gupta, O. Neiman, Metric embedding via shortest path decompositions, SIAM J. Comput. 51 (2022), 290–314. https://doi.org/10.1137/19M1296021
  • A. Filtser, T. Friedrich, D. Issac, N. Kumar, H. Le, N. Mallek, Z. Zeif, Optimal padded decomposition for bounded treewidth graphs, TheoretiCS 4 (2025), Article 22. https://doi.org/10.46298/theoretics.25.22
  • OpenAI, L1L_1L1​ Embeddings of Graphs of Bounded Treewidth, OpenAI Math Release preprint, September 23, 2026 (source of the goal; Theorem 1.1, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/L1-Embeddings-of-Graphs-of-Bounded-Treewidth-September-23-2026/paper.pdf
2 thms2 active usersReviewed
🏆Completed
OptimizationStatistics·Captain: mikedeng1

Random Graphs with a Given Degree Sequence 1: The Fixed-Point Iteration x ↦ φ(x) Converges Geometrically to the Unique β-Model MLE, and Diverges When No MLE ExistsResearch Paper

Motivation

The β\betaβ-model is the simplest exponential random graph model in which every vertex has its own parameter. It was studied by Holland and Leinhardt (1981) in the directed case and by Park and Newman (2004) and Blitzstein and Diaconis (2011) in the undirected case, and it is a close relative of the Bradley–Terry model for paired comparisons. Its sufficient statistic is the degree sequence of the observed graph, so fitting the model to data means solving a system of nnn nonlinear equations in nnn unknowns, one per vertex, where nnn can be in the thousands.

Chatterjee, Diaconis and Sly (arXiv:1005.1136, Ann. Appl. Probab. 21 (2011)) give a simple iteration for this system and prove that it converges geometrically fast whenever a solution exists, at a rate that does not deteriorate with the number of vertices, and that it detects when no solution exists. The same estimate is then used in the paper's consistency theorem for the maximum likelihood estimator (Theorem 1.3) and in its graph-limit theorem (Theorem 1.1). This mission formalizes that algorithmic result, Theorem 1.5, together with the chain of numbered displays in Section 2 that proves it.

Setting

Let n≥3n\ge 3n≥3 and label the vertices 1,…,n1,\dots,n1,…,n. For β=(β1,…,βn)∈Rn\beta=(\beta_1,\dots,\beta_n)\in\mathbb R^nβ=(β1​,…,βn​)∈Rn, the β\betaβ-model puts an edge between distinct vertices iii and jjj independently with probability

pij(β)=eβi+βj1+eβi+βj.p_{ij}(\beta)=\frac{e^{\beta_i+\beta_j}}{1+e^{\beta_i+\beta_j}}.pij​(β)=1+eβi​+βj​eβi​+βj​​.

Given an observed graph with degrees d1,…,dnd_1,\dots,d_nd1​,…,dn​, the maximum likelihood estimate β^\hat\betaβ^​ must satisfy the ML equations

di=∑j≠ieβ^i+β^j1+eβ^i+β^j,i=1,…,n.(3)d_i=\sum_{j\ne i}\frac{e^{\hat\beta_i+\hat\beta_j}}{1+e^{\hat\beta_i+\hat\beta_j}},\qquad i=1,\dots,n.\tag{3}di​=j=i∑​1+eβ^​i​+β^​j​eβ^​i​+β^​j​​,i=1,…,n.(3)

The sup norm of x∈Rnx\in\mathbb R^nx∈Rn is ∣x∣∞=max⁡i∣xi∣|x|_\infty=\max_i|x_i|∣x∣∞​=maxi​∣xi​∣. For i≠ji\ne ji=j put rij(x)=1/(e−xj+exi)r_{ij}(x)=1/(e^{-x_j}+e^{x_i})rij​(x)=1/(e−xj​+exi​), and define the map φ:Rn→Rn\varphi:\mathbb R^n\to\mathbb R^nφ:Rn→Rn by

φi(x)=log⁡di−log⁡∑j≠irij(x).\varphi_i(x)=\log d_i-\log\sum_{j\ne i}r_{ij}(x).φi​(x)=logdi​−logj=i∑​rij​(x).

Starting from any x0∈Rnx_0\in\mathbb R^nx0​∈Rn, the fixed-point iteration is xk+1=φ(xk)x_{k+1}=\varphi(x_k)xk+1​=φ(xk​). The fixed points of φ\varphiφ are exactly the solutions of (3).

For an n×nn\times nn×n matrix A=(aij)A=(a_{ij})A=(aij​), the L∞L^\inftyL∞ operator norm is ∣A∣∞=max⁡i∑j∣aij∣|A|_\infty=\max_i\sum_j|a_{ij}|∣A∣∞​=maxi​∑j​∣aij​∣. For δ>0\delta>0δ>0, the class Ln(δ)\mathcal L_n(\delta)Ln​(δ) consists of the matrices with ∣A∣∞≤1|A|_\infty\le 1∣A∣∞​≤1, aii≥δa_{ii}\ge\deltaaii​≥δ and aij≤−δ/(n−1)a_{ij}\le-\delta/(n-1)aij​≤−δ/(n−1) for i≠ji\ne ji=j. For x,y∈Rnx,y\in\mathbb R^nx,y∈Rn, J(x,y)J(x,y)J(x,y) is the matrix with entries Jij(x,y)=∫01∂φi∂xj(tx+(1−t)y) dtJ_{ij}(x,y)=\int_0^1\frac{\partial\varphi_i}{\partial x_j}(tx+(1-t)y)\,dtJij​(x,y)=∫01​∂xj​∂φi​​(tx+(1−t)y)dt.

Formalization targets

Goal: Theorem 1.5

There are functions Θ(a,b)∈[0,1)\Theta(a,b)\in[0,1)Θ(a,b)∈[0,1) and C(a,b)C(a,b)C(a,b), CCC continuous, independent of nnn and ddd, such that for every n≥3n\ge3n≥3 and every d∈(0,∞)nd\in(0,\infty)^nd∈(0,∞)n: if (3) has a solution β^\hat\betaβ^​, then β^=φ(β^)\hat\beta=\varphi(\hat\beta)β^​=φ(β^​), it is the unique solution of (3), and for all x0x_0x0​ and kkk

∣xk−β^∣∞≤Θ(∣β^∣∞,∣x0∣∞)⌊k/2⌋∣x0−β^∣∞,∣x0−β^∣∞≤C(∣β^∣∞,∣x0∣∞) ∣x0−x1∣∞;|x_k-\hat\beta|_\infty\le\Theta(|\hat\beta|_\infty,|x_0|_\infty)^{\lfloor k/2\rfloor}|x_0-\hat\beta|_\infty,\qquad |x_0-\hat\beta|_\infty\le C(|\hat\beta|_\infty,|x_0|_\infty)\,|x_0-x_1|_\infty;∣xk​−β^​∣∞​≤Θ(∣β^​∣∞​,∣x0​∣∞​)⌊k/2⌋∣x0​−β^​∣∞​,∣x0​−β^​∣∞​≤C(∣β^​∣∞​,∣x0​∣∞​)∣x0​−x1​∣∞​;

and if (3) has no solution, then every orbit {xk}\{x_k\}{xk​} is unbounded.

Milestones, in the order the proof uses them

  1. p. 9: φ(x)=x\varphi(x)=xφ(x)=x if and only if xxx solves (3).
  2. (6)–(7), p. 12: on ∣x∣∞≤K|x|_\infty\le K∣x∣∞​≤K, −e2Kn−1≤∂φi∂xj≤−e−4K2(n−1)-\frac{e^{2K}}{n-1}\le\frac{\partial\varphi_i}{\partial x_j}\le-\frac{e^{-4K}}{2(n-1)}−n−1e2K​≤∂xj​∂φi​​≤−2(n−1)e−4K​ for i≠ji\ne ji=j, and 12e−4K≤∂φi∂xi≤e2K\frac12e^{-4K}\le\frac{\partial\varphi_i}{\partial x_i}\le e^{2K}21​e−4K≤∂xi​∂φi​​≤e2K.
  3. (8), pp. 12–13: φ(x)−φ(y)=J(x,y)(x−y)\varphi(x)-\varphi(y)=J(x,y)(x-y)φ(x)−φ(y)=J(x,y)(x−y); off-diagonal partials are negative, diagonal ones positive, and every row of the Jacobian has absolute sum 111.
  4. Lemma 2.1, p. 11: for A,B∈Ln(δ)A,B\in\mathcal L_n(\delta)A,B∈Ln​(δ),
∣AB∣∞≤1−2(n−2)δ2n−1.|AB|_\infty\le 1-\frac{2(n-2)\delta^2}{n-1}.∣AB∣∞​≤1−n−12(n−2)δ2​.
  1. (9), p. 13: with KKK the largest of ∣x∣∞,∣y∣∞,∣φ(x)∣∞,∣φ(y)∣∞|x|_\infty,|y|_\infty,|\varphi(x)|_\infty,|\varphi(y)|_\infty∣x∣∞​,∣y∣∞​,∣φ(x)∣∞​,∣φ(y)∣∞​ and δ=12e−4K\delta=\frac12e^{-4K}δ=21​e−4K,
∣φ(φ(x))−φ(φ(y))∣∞≤(1−2(n−2)δ2n−1)∣x−y∣∞,|\varphi(\varphi(x))-\varphi(\varphi(y))|_\infty\le\Big(1-\frac{2(n-2)\delta^2}{n-1}\Big)|x-y|_\infty,∣φ(φ(x))−φ(φ(y))∣∞​≤(1−n−12(n−2)δ2​)∣x−y∣∞​,

and ∣φ(x)−φ(y)∣∞≤∣x−y∣∞|\varphi(x)-\varphi(y)|_\infty\le|x-y|_\infty∣φ(x)−φ(y)∣∞​≤∣x−y∣∞​. 6. (10), p. 13: a single θ=Θ(∣β^∣∞,∣x0∣∞)∈[0,1)\theta=\Theta(|\hat\beta|_\infty,|x_0|_\infty)\in[0,1)θ=Θ(∣β^​∣∞​,∣x0​∣∞​)∈[0,1), continuous in its arguments, with ∣xk+3−xk+2∣∞≤θ∣xk+1−xk∣∞|x_{k+3}-x_{k+2}|_\infty\le\theta|x_{k+1}-x_k|_\infty∣xk+3​−xk+2​∣∞​≤θ∣xk+1​−xk​∣∞​ and ∣xk+2−β^∣∞≤θ∣xk−β^∣∞|x_{k+2}-\hat\beta|_\infty\le\theta|x_k-\hat\beta|_\infty∣xk+2​−β^​∣∞​≤θ∣xk​−β^​∣∞​.

Significance

The theorem turns an implicit statistical object, the MLE of an nnn-parameter model, into the limit of an explicit iteration with three guarantees: a convergence rate that depends on the size of the parameters but not on nnn; an a posteriori error bound in terms of the computable step ∣x0−x1∣∞|x_0-x_1|_\infty∣x0​−x1​∣∞​; and a certificate of non-existence. The uniqueness of the MLE comes out of the same estimate. In the paper, the a priori bound ∣x0−β^∣∞≤C∣x0−x1∣∞|x_0-\hat\beta|_\infty\le C|x_0-x_1|_\infty∣x0​−β^​∣∞​≤C∣x0​−x1​∣∞​ is the device that proves consistency of the MLE (Theorem 1.3, applied with x0=βx_0=\betax0​=β), and the two-step contraction is reused in the proof of the graph-limit theorem.

The result is proved in the paper; no machine-checked proof of it is known. The work here is to formalize the known proof. The matrix inequality of Lemma 2.1 and the derivative bounds (6)–(7) are self-contained and can be attacked independently of the iteration.

Difficulty

The obvious approach is to show that φ\varphiφ is a contraction and apply Banach's fixed-point theorem. This fails: every row of the Jacobian of φ\varphiφ has absolute sum exactly 111, so φ\varphiφ is only non-expansive in the sup norm, never strictly contracting. Contraction appears only for the two-step map φ∘φ\varphi\circ\varphiφ∘φ, through the sign pattern of the Jacobian (positive diagonal, uniformly negative off-diagonal), and only on bounded sets, with a factor that tends to 111 as the norms grow. A second difficulty is uniformity: the factor must be bounded away from 111 independently of nnn, which needs (n−2)/(n−1)≥12(n-2)/(n-1)\ge\frac12(n−2)/(n−1)≥21​, hence n≥3n\ge3n≥3. For the converse, there is no fixed point to anchor the orbit, so boundedness of the orbit has to be converted into convergence before a fixed point exists.

Formalization scope

  • Vertices are Fin n ={0,…,n−1}=\{0,\dots,n-1\}={0,…,n−1}; vectors are Fin n → ℝ, whose Mathlib norm is the sup norm ∣x∣∞|x|_\infty∣x∣∞​. Matrix norms are the explicit maximal absolute row sum.
  • The degrees did_idi​ are arbitrary positive reals. This is a generalization of the page (degrees of a graph), and it is all the proof uses; positivity is required because log⁡di\log d_ilogdi​ is junk in Lean at 000.
  • n≥3n\ge3n≥3 is a hypothesis of the goal and of (10). It is not on the page but is necessary: for n=2n=2n=2, (3) reads d1=d2=p12d_1=d_2=p_{12}d1​=d2​=p12​, whose solutions form a line, so uniqueness fails, and the factor in Lemma 2.1 equals 111.
  • Partial derivatives are Fréchet derivatives applied to unit vectors; J(x,y)J(x,y)J(x,y) uses the interval integral over [0,1][0,1][0,1].
  • "Geometrically fast, with rate depending only on (∣β^∣∞,∣x0∣∞)(|\hat\beta|_\infty,|x_0|_\infty)(∣β^​∣∞​,∣x0​∣∞​)" is stated with the function Θ\ThetaΘ quantified before nnn, ddd, β^\hat\betaβ^​ and x0x_0x0​, and "CCC is a continuous function of the pair" likewise. A version that chooses the rate after fixing nnn, ddd and x0x_0x0​ is strictly weaker and nearly says only that a convergent sequence converges; it is not the target.
  • "A divergent subsequence" is stated as unboundedness of {∣xk∣∞}\{|x_k|_\infty\}{∣xk​∣∞​}, which is equivalent for sequences in Rn\mathbb R^nRn.
  • The paper's remark that θ\thetaθ is "uniformly bounded away from 1 on subsets of Rn×Rn\mathbb R^n\times\mathbb R^nRn×Rn" holds only on bounded subsets; it is not stated separately, and its quantitative content is (10).

Needed infrastructure: calculus of φ\varphiφ (derivatives of log⁡\loglog of a sum of exponentials), the integral form of the mean value theorem for maps Rn→Rn\mathbb R^n\to\mathbb R^nRn→Rn, and elementary matrix-norm estimates. The class Ln(δ)\mathcal L_n(\delta)Ln​(δ) and Lemma 2.1 are reusable for other diagonally dominant iterations. Proofs of any milestone are welcome, in any order.

Selected references

  • S. Chatterjee, P. Diaconis and A. Sly, Random graphs with a given degree sequence, Ann. Appl. Probab. 21(4), 1400–1435, 2011. arXiv:1005.1136v5, DOI 10.1214/10-AAP728
  • P. W. Holland and S. Leinhardt, An exponential family of probability distributions for directed graphs, J. Amer. Statist. Assoc. 76, 33–65, 1981. DOI 10.1080/01621459.1981.10477598
  • J. Park and M. E. J. Newman, Statistical mechanics of networks, Phys. Rev. E 70, 066117, 2004. arXiv:cond-mat/0405566
  • J. Blitzstein and P. Diaconis, A sequential importance sampling algorithm for generating random graphs with prescribed degrees, Internet Mathematics 6(4), 489–522, 2011. DOI 10.1080/15427951.2010.557277
9 thms2 active usersReviewed
🏆Completed
CombinatoricsDiscrete GeometryLinear algebra·Captain: mikedeng1

The Geometry of Graphs and Some of Its Algorithmic Applications 4: The m-Cube Has Isometric Dimension mResearch Paper

Why graph dimension matters

A finite graph has distances defined by the lengths of its shortest paths. Representing those distances inside a normed vector space lets geometric language describe a combinatorial object, but the number of coordinates required can depend sharply on the graph. The isometric dimension asks for the smallest dimension of any real normed space that preserves every graph distance exactly. Linial, London, and Rabinovich study this invariant alongside metric embeddings and algorithmic uses of graph geometry in their 1995 paper. The cube is a useful test case: its vertices already have an obvious coordinate description, while the theorem must show that no different norm compresses those exact distances into fewer coordinates.

The paper's Section 5.3 relates isometric dimension to the rank of a matrix recording vertex distances in a maximum norm. It then computes the dimension of an even cycle and uses an isometric cycle inside the cube. These are known results of the paper, not an open conjecture. The mission isolates the rank characterization, the cycle calculation, and the cube calculation as reusable formal targets.

Graphs, norms, and the cube

For a finite connected simple graph G=(V,E)G=(V,E)G=(V,E), let dG(x,y)d_G(x,y)dG​(x,y) be the minimum number of edges in a path from xxx to yyy. A norm NNN on Rd\mathbb R^dRd is nonnegative, vanishes only at zero, satisfies N(av)=∣a∣N(v)N(av)=|a|N(v)N(av)=∣a∣N(v), and obeys N(u+v)≤N(u)+N(v)N(u+v)\leq N(u)+N(v)N(u+v)≤N(u)+N(v). An isometric embedding is a map φ:V→Rd\varphi:V\to\mathbb R^dφ:V→Rd such that N(φ(x)−φ(y))=dG(x,y)N(\varphi(x)-\varphi(y))=d_G(x,y)N(φ(x)−φ(y))=dG​(x,y) for every pair of vertices. The number dim⁡(G)\dim(G)dim(G) is the least ddd for which some norm and some such map exist. The paper states the more general metric-space definition in Definition 2.1 and adopts connected graphs as the standing convention in Section 5.1 (pp. 218–219, 229).

The mmm-cube has binary strings x=(x1,…,xm)x=(x_1,\ldots,x_m)x=(x1​,…,xm​) as vertices. Two strings are adjacent when exactly one coordinate differs. Its graph distance is the Hamming distance, the number of differing coordinates. For m≥2m\geq2m≥2, the even cycle C2mC_{2m}C2m​ has 2m2m2m vertices in cyclic order, each adjacent to its two neighbors. A graph JJJ is an isometric subgraph of GGG when its distances between its own vertices agree with their distances in GGG.

An ∣V∣×r|V|\times r∣V∣×r real matrix MMM implements GGG when the maximum absolute coordinate difference between rows iii and jjj is dG(i,j)d_G(i,j)dG​(i,j), for every pair. The number of columns rrr is allowed to vary; the matrix rank is its ordinary linear-algebraic rank. This is the matrix object in Theorem 5.7, p. 232.

Formalization targets

The first main intermediate target is the rank characterization

dim⁡(G)=min⁡M implements Grank⁡(M).\dim(G)=\min_{M\text{ implements }G}\operatorname{rank}(M).dim(G)=M implements Gmin​rank(M).

The next is the even-cycle result, stated for m≥2m\geq2m≥2:

dim⁡(C2m)=m.\dim(C_{2m})=m.dim(C2m​)=m.

The goal is Corollary 5.12, p. 233, for every positive integer mmm:

dim⁡(m-Cube)=m.\boxed{\dim(m\text{-Cube})=m.}dim(m-Cube)=m.​

The milestone list also records the isometric-subgraph lower bound, the nonsingularity of the distance matrix of a path, an isometric embedding of the cube into ℓ1m\ell_1^mℓ1m​, and an isometric C2mC_{2m}C2m​ inside the cube. The paper's proof of the corollary names the last two facts in consecutive sentences (p. 233).

What the result establishes

The upper bound says the familiar mmm binary coordinates suffice: the cube embeds in ℓ1m\ell_1^mℓ1m​ with all distances preserved. The lower bound is stronger than a claim about that particular norm or embedding. It rules out every real normed space of dimension below mmm, regardless of its norm. The even cycle provides a separate exact dimension calculation, while Theorem 5.7 gives a matrix-rank way to study other finite graph metrics. These statements separate the graph invariant from whichever embedding happens to be easiest to draw (Section 5.3, pp. 232–233).

The paper proves these claims mathematically. The remaining formalization work is to build machine-checked proofs of their precise statements, including the finite-dimensional norm interface, graph-distance calculations, implementing matrices, and the rank lower bound. The definitions of normed realization and implementing matrix can be reused for other graphs; the path-distance determinant can be used independently of graph embedding questions. The mission's theorem files are open proof targets, not claims that these proofs have already been checked in Lean.

Where the difficulty lies

Placing binary cube vertices in ℓ1m\ell_1^mℓ1m​ gives an immediate candidate for an upper bound, but it cannot by itself exclude embeddings into lower-dimensional spaces with a different norm. The lower bound has to control all norms. Theorem 5.7 converts that unrestricted geometric question into a rank statement for every implementing matrix, and Proposition 5.10 supplies the exact cycle dimension needed by the cube argument (pp. 232–233).

The paper's printed proof of Proposition 5.10 first obtains only rank⁡(A)≥m−1\operatorname{rank}(A)\geq m-1rank(A)≥m−1, then says a more careful analysis yields mmm without spelling it out. That final rank step is a genuine proof obligation for solvers; the stated proposition and its use in Corollary 5.12 are retained exactly (p. 233). The path-distance matrix cited in that proof is nonsingular when its size is at least two; size one would be a zero matrix.

Formalization scope

Lean represents an arbitrary ddd-dimensional real normed space by a norm function on Fin d → ℝ. This includes every finite-dimensional real normed space up to a choice of coordinates. The isoDim definition takes the least natural-number dimension admitting an isometric embedding, with equality of distances for every vertex pair. The graph metric comes from SimpleGraph.dist; connectedness is explicit where the paper assumes it. The cube is the graph on Fin m → Bool with Hamming-distance-one adjacency, and C2mC_{2m}C2m​ is Mathlib's cycleGraph (2 * m). The goal assumes m≥1m\geq1m≥1; its m=1m=1m=1 case is direct, while the cycle milestones use m≥2m\geq2m≥2. A zero-dimensional cube consists of one vertex and is covered by the auxiliary upper-bound statement, though the goal follows the paper's positive-dimensional convention.

An implementing matrix has arbitrary finite column count. Instead of taking a maximum over an empty index set, the formal definition requires all column differences to be bounded by the graph distance and a column to attain that distance for distinct rows. This agrees with the paper's maximum formula whenever columns are present and handles a one-vertex graph with zero columns. Matrix rank is taken over R\mathbb RR; the path-distance matrix has entries ∣r−s∣|r-s|∣r−s∣ for indices 0,…,m−10,\ldots,m-10,…,m−1. The infimum over achievable dimensions or ranks is a minimum, since these are nonempty subsets of N\mathbb NN for finite connected graphs. No asymptotic O(⋅)O(\cdot)O(⋅) or Ω(⋅)\Omega(\cdot)Ω(⋅) bound occurs in these targets, so there is no hidden constant to instantiate.

The statements require genuine norm axioms and exact graph distances. Allowing a seminorm, assuming the desired dimension bound, or replacing the cube by a path would empty the main claim. Useful contributions include the cube-distance calculation, the path determinant, an exact even-cycle lower bound, and the general rank characterization. The paper makes no running-time claim for these particular results; no algorithmic complexity assertion is included.

Selected references

  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15(2), 215–245, 1995. DOI: 10.1007/BF01200757.
  • J. H. van Lint, R. M. Wilson, A Course in Combinatorics, Cambridge University Press, 1992, pp. 64–65, cited as [44] by Linial et al. for the path distance-matrix fact. Book record.
11 thms2 active usersReviewed
🏆Completed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Multicommodity Max-Flow Min-Cut Theorems and Their Use in Designing Approximation Algorithms V: On a 3-Regular c-Expander the Uniform Max-Flow Is at Most 6𝒮/(c(⌊log n⌋ − 2))Research Paper

Motivation

A multicommodity flow problem asks to route many source–sink demands through one capacitated network at the same time. Its min-cut (sparsest cut) is the smallest ratio, over all cuts, of the capacity crossing the cut to the demand separated by it. Routing is always limited by every cut, so the max-flow never exceeds the min-cut. For a single commodity the two are equal (Ford–Fulkerson); for many commodities they are not, and the size of the gap governs a large family of approximation algorithms for graph partitioning, layout and routing problems.

Leighton and Rao (J. ACM 46 (1999), extended abstract FOCS 1988) proved that for uniform problems, where every pair of nodes demands one unit, the gap is at most O(log⁡n)O(\log n)O(logn) (their Theorem 2). This mission formalizes the other half: Section 2.1 of the paper, A bad example, which shows that the log⁡n\log nlogn is necessary. On a 3-regular expander graph the max-flow is a Θ(log⁡n)\Theta(\log n)Θ(logn) factor smaller than the min-cut, so Theorem 2 cannot be improved beyond the constant.

Timeline. Ford and Fulkerson (1956) and Elias–Feinstein–Shannon (1956): single-commodity max-flow equals min-cut. Hu (1963): equality for two commodities. Leighton and Rao (1988, journal 1999): uniform and product problems have an O(log⁡n)O(\log n)O(logn) gap, and expanders show it is tight. Klein, Agrawal, Ravi and Rao (1990) and later Linial–London–Rabinovich (Combinatorica 1995) and Aumann–Rabani (1998): an O(log⁡k)O(\log k)O(logk) gap for general demands with kkk commodities, tight by the same expander example.

Setting

Let VVV be a finite set of nnn nodes. A network assigns a capacity C(u,v)≥0C(u,v)\ge0C(u,v)≥0 to every pair, symmetric and zero on the diagonal. For U⊆VU\subseteq VU⊆V with complement Uˉ\bar UUˉ, the cut capacity C(U,Uˉ)C(U,\bar U)C(U,Uˉ) is the total capacity of the pairs with one end in UUU and the other in Uˉ\bar UUˉ.

In the uniform multicommodity flow problem (UMFP) every pair of distinct nodes is a commodity with unit demand. Its max-flow fff is the largest fraction such that fff units of every commodity can be routed simultaneously without exceeding any capacity (concurrent flow). Its min-cut is

S=min⁡∅≠U⊊VC(U,Uˉ)∣U∣ ∣Uˉ∣.\mathcal S=\min_{\emptyset\ne U\subsetneq V}\frac{C(U,\bar U)}{|U|\,|\bar U|}.S=∅=U⊊Vmin​∣U∣∣Uˉ∣C(U,Uˉ)​.

Section 2.1 takes a 3-regular nnn-node graph GGG with unit edge capacities, so that C(U,Uˉ)=∣⟨U,Uˉ⟩∣C(U,\bar U)=|\langle U,\bar U\rangle|C(U,Uˉ)=∣⟨U,Uˉ⟩∣ is the number of edges between UUU and Uˉ\bar UUˉ, and assumes that GGG is a ccc-edge-expander:

∣⟨U,Uˉ⟩∣ ≥ c min⁡{∣U∣,∣Uˉ∣}for all U⊆V,|\langle U,\bar U\rangle|\ \ge\ c\,\min\{|U|,|\bar U|\}\qquad\text{for all }U\subseteq V,∣⟨U,Uˉ⟩∣ ≥ cmin{∣U∣,∣Uˉ∣}for all U⊆V,

for a constant c>0c>0c>0. Families of such graphs exist for a sufficiently small ccc (Margulis 1973; a random 3-regular graph has the property with high probability). Throughout, log⁡\loglog is log⁡2\log_2log2​.

Formalization targets

Goal: the gap on an expander

For a 3-regular ccc-edge-expander with unit capacities and n≥8n\ge 8n≥8,

f ≤ 6 Sc (⌊log⁡n⌋−2).f\ \le\ \frac{6\,\mathcal S}{c\,\bigl(\lfloor\log n\rfloor-2\bigr)} .f ≤ c(⌊logn⌋−2)6S​.

This is the final display of §2.1 (pp. 794–795), with the floor that the argument actually supports (see Formalization scope). It leaves the expander family unspecified: it holds for every graph with the two properties.

Milestones, in the order the argument uses them

  1. S≥c/(n−1)\mathcal S\ge c/(n-1)S≥c/(n−1) for the expander (p. 794).
  2. In a 3-regular graph, at most n/2n/2n/2 nodes lie within distance ⌊log⁡n⌋−3\lfloor\log n\rfloor-3⌊logn⌋−3 of any node.
  3. Hence at least half of the source–sink pairs are at distance at least ⌊log⁡n⌋−2\lfloor\log n\rfloor-2⌊logn⌋−2.
  4. A uniform flow of fraction fff uses capacity at least 12(n2)f(⌊log⁡n⌋−2)\tfrac12\binom n2 f(\lfloor\log n\rfloor-2)21​(2n​)f(⌊logn⌋−2), and the total capacity is 3n/23n/23n/2.
  5. Hence f≤3n/((n2)(⌊log⁡n⌋−2))=6/((n−1)(⌊log⁡n⌋−2))f\le 3n/\bigl(\binom n2(\lfloor\log n\rfloor-2)\bigr)=6/\bigl((n-1)(\lfloor\log n\rfloor-2)\bigr)f≤3n/((2n​)(⌊logn⌋−2))=6/((n−1)(⌊logn⌋−2)).

Companions: Theorem 1

Theorem 1 (p. 795) states that for any nnn there is an nnn-node UMFP with f≤O(S/log⁡n)f\le O(\mathcal S/\log n)f≤O(S/logn). It is stated twice: conditionally on a family of 3-regular ccc-expanders (the form §2.1 proves), and unconditionally, with one constant KKK for all n≥2n\ge2n≥2.

Significance

The result shows that the O(log⁡n)O(\log n)O(logn) max-flow min-cut gap of Leighton and Rao's Theorem 2 is tight up to the constant, and, through the same example, that the O(log⁡k)O(\log k)O(logk) bound for general demands is tight. Every approximation algorithm that rounds a multicommodity flow, or its dual, to a cut loses a Θ(log⁡n)\Theta(\log n)Θ(logn) factor in the worst case for this reason; improving on it, as Arora, Rao and Vazirani did (J. ACM 2009) with an O(log⁡n)O(\sqrt{\log n})O(logn​) approximation for sparsest cut, requires a stronger relaxation than the flow LP.

The result is classical and proved. To our knowledge it has no machine-checked proof. Formalizing it requires the concurrent-flow LP, the uniform sparsest cut, graph distance balls in bounded-degree graphs and a counting argument that converts path lengths into capacity usage; each of these is reusable. The unconditional Theorem 1 additionally needs the existence of bounded-degree expanders on every number of nodes, which is a substantial formalization project in its own right.

Difficulty

The min-cut bound and the ball bound are elementary. The central step is milestone 4: turning "the shortest path of a commodity has at least ℓ\ellℓ edges" into "the commodity uses at least fℓf\ellfℓ units of capacity". In an arc-flow formulation a commodity is a conservation-respecting function on arcs, not a set of paths, so the step needs either a path decomposition of each commodity's flow or a weak-duality argument with unit edge lengths. The obvious attempt, summing the conservation equations directly, gives only the net flow across one cut and loses the factor ℓ\ellℓ.

The unconditional Theorem 1 is hard for a different reason: the paper cites the existence of 3-regular expanders and does not prove it, and the theorem asks for every nnn, including odd nnn, where no 3-regular graph exists.

Formalization scope

Nodes are a Fintype VVV with decidable equality and n=∣V∣n=|V|n=∣V∣. A network is a symmetric nonnegative real capacity function with zero diagonal; the unit-capacity network of a simple graph GGG has C(u,v)=1C(u,v)=1C(u,v)=1 on edges and 000 elsewhere. Commodities are ordered pairs (s,t)(s,t)(s,t), s≠ts\ne ts=t, each with demand 12\tfrac1221​ (the paper's footnote 2), so each unordered pair carries one unit; a flow is a nonnegative arc flow per commodity, and both directions of an edge share its capacity. The max-flow is the supremum of feasible fractions. The min-cut ranges over nonempty proper UUU only: at U=∅U=\emptysetU=∅ or VVV the ratio is 0/00/00/0, which Lean reads as 000 and which would make S=0\mathcal S=0S=0.

Conventions and added hypotheses:

  • The floor. The page writes log⁡n−2\log n-2logn−2. Its argument shows that far pairs have distance at least ⌊log⁡n⌋−2\lfloor\log n\rfloor-2⌊logn⌋−2, which differs from log⁡n−2\log n-2logn−2 for non-integer log⁡n\log nlogn. The goal and milestones state ⌊log⁡n⌋−2\lfloor\log n\rfloor-2⌊logn⌋−2 (Nat.log 2 n), a different constant in the same Θ(log⁡n)\Theta(\log n)Θ(logn).
  • n≥8n\ge 8n≥8 makes ⌊log⁡n⌋−2≥1\lfloor\log n\rfloor-2\ge1⌊logn⌋−2≥1, so the bound divides by a positive number.
  • Graph distance is the extended distance, ∞\infty∞ between components, never Lean's junk value 000.
  • Regularity is SimpleGraph.IsRegularOfDegree 3.
  • Theorem 1 requires the network to be connected (the paper's standing assumption, p. 789); a disconnected network has f=S=0f=\mathcal S=0f=S=0 and would make the statement empty. Its constant is quantified before nnn.
  • Polynomial-time claims are not formalized.

The hypotheses of the goal are satisfiable: the 3-cube Q3Q_3Q3​ (n=8n=8n=8) is 3-regular and a 111-edge-expander. Contributions welcome: the arc-flow path-length lemma, the ball bound for bounded-degree graphs, and an existence theorem for 3-regular expanders.

Selected references

  • T. Leighton and S. Rao, Multicommodity max-flow min-cut theorems and their use in designing approximation algorithms, J. ACM 46(6):787–832, 1999. https://doi.org/10.1145/331524.331526
  • G. A. Margulis, Explicit constructions of concentrators, Problems of Information Transmission 9(4):325–332, 1973.
  • N. Linial, E. London and Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15(2):215–245, 1995. https://doi.org/10.1007/BF01200757
  • S. Arora, S. Rao and U. Vazirani, Expander flows, geometric embeddings and graph partitioning, J. ACM 56(2):5, 2009. https://doi.org/10.1145/1502793.1502794
7 thms2 active usersReviewed
🏆Completed
CombinatoricsOperations Research·Captain: mikedeng1

An Introduction to Timetabling 3: Course Scheduling with Unavailabilities and Preassignments in p Periods Is Equivalent to p-Coloring the Nodes of a GraphResearch Paper

Motivation

Universities schedule lectures and examinations every term. When students choose their courses freely, with no fixed curriculum, the basic requirement is that no student is asked to attend two lectures at the same time. On top of this come local constraints: a course cannot take place at some periods (its teacher or a suitable room is unavailable), and some lectures are fixed in advance at a given period. D. de Werra's invited review An introduction to timetabling (Eur. J. Oper. Res. 19 (1985) 151–162) surveys the graph models of these problems. Its §3.1 shows that course scheduling, with unavailabilities and preassignments, is exactly a node colouring problem. This reduction is the standard starting point for graph colouring heuristics in timetabling, including the examination timetabling literature.

Setting

A course scheduling instance has qqq courses KbK_bKb​, course KbK_bKb​ consisting of kbk_bkb​ lectures lal_ala​, and rrr students; student iii takes a set TiT_iTi​ of courses. There are ppp periods.

The lecture graph GGG has one lecture-node mabm_{ab}mab​ for each lecture lal_ala​ of each course KbK_bKb​. Two lecture-nodes of the same course are adjacent. Two lecture-nodes of distinct courses KbK_bKb​, KbˉK_{\bar b}Kbˉ​ are adjacent whenever some student takes both courses.

A feasible course schedule in ppp periods assigns a period s(l)∈{1,…,p}s(l) \in \{1,\dots,p\}s(l)∈{1,…,p} to each lecture so that two lectures of one course are at different periods and no student has two lectures at the same period.

A course scheduling problem with unavailabilities and preassignments (CSUP) adds, for each course KbK_bKb​, a set UbU_bUb​ of forbidden periods, and for some lectures lll a prescribed period kˉ(l)\bar k(l)kˉ(l). A solution is a feasible schedule with s(l)∉Ubs(l) \notin U_bs(l)∈/Ub​ for every lecture lll of KbK_bKb​ and s(l)=kˉ(l)s(l) = \bar k(l)s(l)=kˉ(l) for every preassigned lecture.

The graph G^\hat GG^ of a CSUP adds one period-node kkk per period to GGG. All pairs of period-nodes are adjacent. A lecture-node of KbK_bKb​ is adjacent to period-node kkk if k∈Ubk \in U_bk∈Ub​. A lecture-node of a lecture preassigned to kˉ\bar kkˉ is adjacent to every period-node k≠kˉk \neq \bar kk=kˉ.

A node colouring with ppp colours assigns one of ppp colours to each node so that adjacent nodes get different colours.

Formalization targets

Goal: Proposition 3.1

the CSUP P has a solution in p periods  ⟺  G^(P) has a node colouring with p colours.\text{the CSUP } P \text{ has a solution in } p \text{ periods} \iff \hat G(P) \text{ has a node colouring with } p \text{ colours.}the CSUP P has a solution in p periods⟺G^(P) has a node colouring with p colours.

Milestones

  1. §3.1 (p. 157). Without unavailabilities and preassignments, the feasible schedules in ppp periods are exactly the node colourings of GGG with ppp colours.
  2. Figure 4 (p. 157). For courses K1,K2,K3K_1,K_2,K_3K1​,K2​,K3​ with 1,2,21,2,21,2,2 lectures, student 1 taking K1,K2K_1,K_2K1​,K2​ and student 2 taking K2,K3K_2,K_3K2​,K3​, the printed schedule in 444 periods is feasible.
  3. Proof of Proposition 3.1 (p. 158). Every colouring of G^\hat GG^ with ppp colours gives a solution in which lecture lal_ala​ of KbK_bKb​ is at period kkk exactly when mabm_{ab}mab​ has the colour of period-node kkk, and every solution arises in this way from some colouring.
  4. Proposition 3.2 (p. 158). Let HHH be a graph on lecture-nodes and a clique of ppp period-nodes, with some lecture-nodes precoloured consistently with the colours of period-nodes. Removing the precoloured nodes and joining their neighbours to the matching period-nodes gives a graph G^′\hat G'G^′ with
H has a p-colouring respecting the precolouring  ⟺  G^′ is p-colourable.H \text{ has a } p\text{-colouring respecting the precolouring} \iff \hat G' \text{ is } p\text{-colourable.}H has a p-colouring respecting the precolouring⟺G^′ is p-colourable.

Significance

The result. Proposition 3.1 turns a family of constrained timetabling problems into one well-studied problem. Any exact or heuristic method for node colouring (sequential colouring, DSATUR, implicit enumeration; §3.2 of the paper) applies directly to course scheduling with forbidden and fixed periods, and lower bounds on the chromatic number become lower bounds on the number of periods. Proposition 3.2 adds that preassignments cost nothing for node colouring: a precoloured instance is again a plain colouring instance. For edge colouring the paper contrasts this with the NP-completeness of precoloured edge colouring.

The formalization. The results are proved on the page in a few lines and are not machine-checked anywhere. This mission makes the reduction precise. It writes out the graph G^\hat GG^ as a definition, isolates the consistency condition that Proposition 3.2 needs, and states the correspondence between colourings and schedules in a form other formal timetabling and scheduling developments can use. The definitions of course instances, CSUP and G^\hat GG^ can serve formal work on examination timetabling and on the complexity of timetabling problems.

Difficulty

The forward direction of Proposition 3.1 is direct: colour each lecture-node with its period and each period-node with itself. The converse is where care is needed. A colouring of G^\hat GG^ need not give period-node kkk the colour kkk; the colours of the period-nodes can be any permutation of the palette. The schedule must therefore be read through the colours of the period-nodes, and the fact that the ppp period-nodes form a clique is what makes this well defined.

For Proposition 3.2, the obvious reading of the page is false. If two adjacent lecture-nodes are precoloured with the same colour, no colouring respects the precolouring, yet removing both nodes deletes the edge between them and the reduced graph may still be colourable. The same happens when a node is precoloured with the colour of a period-node it is adjacent to. The proposition therefore needs a consistency condition the page leaves implicit. The page's "1–1 correspondence" between colourings fails when precolours are fixed colours of the palette, since G^′\hat G'G^′'s colourings can permute the colours of the period clique; the target is the equivalence of the proposition.

Formalization scope

Courses, students, lectures of a course and periods are indexed from 000 (Fin q, Fin r, Fin (lectures b), Fin p). A lecture is a dependent pair Σ b, Fin (lectures b), and the nodes of G^\hat GG^ are Lecture ⊕ Fin p. Graphs and colourings are Mathlib's SimpleGraph, SimpleGraph.Coloring and SimpleGraph.Colorable. Unavailability is per course and preassignment per lecture, as on the page. The set-valued constraint of the Figure 5 example ("1 lecture of K3K_3K3​ at period 1 or 2") is not modelled. Feasibility includes the clique on the lectures of each course even when no student takes the course, because the page's graph contains it. p=0p = 0p=0 is allowed.

The paper states Proposition 3.1 as "one can construct a graph G^\hat GG^ such that …". The existential form is trivially true (pick a complete graph on p+1p+1p+1 nodes or an edgeless graph according to the answer), so the goal is stated about the specific graph csupGraph P built in the proof.

Proposition 3.2 is stated for any graph whose period-nodes form a clique, with two disclosed consistency hypotheses: no precoloured node is adjacent to its own period-node, and no two adjacent nodes are precoloured with the same period-node's colour. A sanity file shows that the equivalence fails without them. "Precoloured" means "receives the colour of period-node kkk", the reading of the proof of Proposition 3.1.

Not targeted: the NP-completeness remarks (cited, not proved), the heuristics of §3.2–§3.4, and the hypergraph classroom model of §3.5, whose "1–1 correspondence" has the same palette-permutation issue as Proposition 3.2.

Contributions welcome: proofs of the milestones and the goal, and general lemmas about colourings of graphs containing a clique of size equal to the number of colours.

Selected references

  • D. de Werra, An introduction to timetabling, European Journal of Operational Research 19 (1985) 151–162. https://doi.org/10.1016/0377-2217(85)90167-5
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979 (NP-completeness of graph ppp-colourability, reference [13] of the paper).
  • D. Brélaz, New methods to color the vertices of a graph, Communications of the ACM 22 (1979) 251–256. https://doi.org/10.1145/359094.359101
7 thms2 active usersReviewed
🏆Completed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

A Theorem on Boolean Matrices: Warshall's Construction M* Equals the Boolean Sum M ∨ M² ∨ ⋯ ∨ Mᵈ of the Boolean Powers of MResearch Paper

Motivation

A finite directed graph on ddd vertices is recorded by its adjacency matrix, a d×dd \times dd×d array of zeros and ones. Many questions about the graph are questions about reachability: can vertex jjj be reached from vertex iii by following one or more edges? In the late 1950s this question arose in the analysis of programs. Prosser (Applications of Boolean matrices to the analysis of flow diagrams, 1959) represented the flow diagram of a program by a boolean matrix and read off properties such as which blocks can follow which; Marimont (A new method of checking the consistency of precedence matrices, J. ACM 1959) used the same matrices to check precedence relations for consistency.

In both works the reachability matrix is described as a boolean sum of boolean powers of the adjacency matrix. Computed that way, it costs ddd boolean matrix products, that is on the order of d3d^3d3 operations per product. Stephen Warshall's two-page note (A theorem on Boolean matrices, J. ACM 9 (1962) 11–12) proves that a much shorter procedure, ddd passes of row updates performed in place, gives the same matrix. The same procedure was published in the same year by Floyd as Algorithm 96: Ancestor (Comm. ACM 5 (1962) 344), and its weighted version is Floyd's Algorithm 97: Shortest Path. The procedure is now taught as the Floyd–Warshall or Roy–Warshall algorithm.

Timeline. 1959: Prosser and Marimont describe reachability as a sum of boolean powers. 1959: Roy publishes the same closure construction in French (C. R. Acad. Sci. Paris 249). 1962: Warshall's note proves the in-place construction correct; Floyd's Algorithm 96 publishes it as code, and Algorithm 97 extends it to shortest paths.

Setting

Fix a natural number ddd and let M=(mij)M = (m_{ij})M=(mij​) be a d×dd \times dd×d matrix each of whose entries is 000 or 111.

  • The boolean product of two such matrices AAA and BBB is the matrix A∧BA \wedge BA∧B with (A∧B)ij=⋁k(aik∧bkj)(A \wedge B)_{ij} = \bigvee_k (a_{ik} \wedge b_{kj})(A∧B)ij​=⋁k​(aik​∧bkj​): the entry is 111 when some index kkk has aik=bkj=1a_{ik} = b_{kj} = 1aik​=bkj​=1.
  • The boolean sum A∨BA \vee BA∨B is the entrywise disjunction, (A∨B)ij=aij∨bij(A \vee B)_{ij} = a_{ij} \vee b_{ij}(A∨B)ij​=aij​∨bij​.
  • The boolean powers are M1=MM^1 = MM1=M and Mp+1=Mp∧MM^{p+1} = M^p \wedge MMp+1=Mp∧M.
  • The introduction's matrix is the power sum
M′  =  ⋁p=1dMp.M' \;=\; \bigvee_{p=1}^{d} M^p .M′=p=1⋁d​Mp.
  • The THEOREM defines M′M'M′ instead by chains: mij′=1m'_{ij} = 1mij′​=1 if and only if either mij=1m_{ij} = 1mij​=1 or there are indices k1,…,knk_1, \dots, k_nk1​,…,kn​ with mik1=mk1k2=⋯=mknj=1m_{ik_1} = m_{k_1k_2} = \cdots = m_{k_nj} = 1mik1​​=mk1​k2​​=⋯=mkn​j​=1. Footnote 2 states that the two definitions agree.
  • Warshall's construction M∗M^*M∗: set M∗=MM^* = MM∗=M; for i=1,…,di = 1, \dots, di=1,…,d in turn, for every row jjj with mji∗=1m^*_{ji} = 1mji∗​=1 and every column kkk, set mjk∗=mjk∗∨mik∗m^*_{jk} = m^*_{jk} \vee m^*_{ik}mjk∗​=mjk∗​∨mik∗​. Entries are overwritten during a pass and later steps of the same pass read the new values.
  • Footnote 3 gives a recursive version: (mij)0=mij(m_{ij})_0 = m_{ij}(mij​)0​=mij​ and (mij)n+1=(mij)n∨((mi,n+1)n∧(mn+1,j)n)(m_{ij})_{n+1} = (m_{ij})_n \vee \bigl((m_{i,n+1})_n \wedge (m_{n+1,j})_n\bigr)(mij​)n+1​=(mij​)n​∨((mi,n+1​)n​∧(mn+1,j​)n​), and asserts mij∗=(mij)dm^*_{ij} = (m_{ij})_dmij∗​=(mij​)d​.

In Lean the matrix is M : Fin d → Fin d → Bool. The power sum is powerSum M, the chain relation ChainRel M i j, the recursion warshallRec M n, all in the namespace WarshallBool.Closure. The construction M∗M^*M∗ is the published FloydAlgorithms.ShortestPath.algorithm96 M.

Formalization targets

Goal: the construction computes the power sum

M∗  =  ⋁p=1dMpfor every d and every d×d 0-1 matrix M.M^* \;=\; \bigvee_{p=1}^{d} M^p \qquad \text{for every } d \text{ and every } d \times d \text{ 0-1 matrix } M .M∗=p=1⋁d​Mpfor every d and every d×d 0-1 matrix M.

This is the transformation the introduction asks for, combined with the THEOREM's assertion M∗=M′M^* = M'M∗=M′ and footnote 2. The goal has no hypothesis beyond the shape of MMM.

Milestone 1 (Proved, referenced): the THEOREM

mij∗=1  ⟺  mij=1 or mik1=⋯=mknj=1 for some k1,…,kn.m^*_{ij} = 1 \iff m_{ij} = 1 \text{ or } m_{ik_1} = \cdots = m_{k_nj} = 1 \text{ for some } k_1, \dots, k_n .mij∗​=1⟺mij​=1 or mik1​​=⋯=mkn​j​=1 for some k1​,…,kn​.

It is already on the platform as FloydAlgorithms.ShortestPath.algorithm96_eq_transGen, stated with Mathlib's transitive closure Relation.TransGen, and it is not posed again.

Milestone 2: footnote 2

mij=1 or a chain i→k1→⋯→kn→j of 1 entries  ⟺  (⋁p=1dMp)ij=1.m_{ij} = 1 \text{ or a chain } i \to k_1 \to \cdots \to k_n \to j \text{ of 1 entries} \iff \Bigl(\bigvee_{p=1}^{d} M^p\Bigr)_{ij} = 1 .mij​=1 or a chain i→k1​→⋯→kn​→j of 1 entries⟺(p=1⋁d​Mp)ij​=1.

Companion: footnote 3

M∗  =  ((mij)d)i,j.M^* \;=\; \bigl((m_{ij})_d\bigr)_{i,j} .M∗=((mij​)d​)i,j​.

This is not a milestone of the goal; it identifies the in-place construction with its stage-by-stage recursion.

Significance

The result. The theorem says that reachability in a directed graph on ddd vertices, defined through matrix powers, can be computed by ddd passes of in-place row operations, on the order of d2d^2d2 row updates instead of ddd full boolean matrix products. The same scheme, with (min⁡,+)(\min, +)(min,+) in place of (∨,∧)(\vee, \wedge)(∨,∧), is the Floyd–Warshall shortest-path algorithm; with other semirings it computes regular expressions of finite automata (Kleene, McNaughton–Yamada) and path problems in general. Footnote 3's recursion is the form in which the algorithm is usually proved correct today, and its equivalence with the in-place version is what justifies running it without a second copy of the matrix.

Formalizing it. The chain form of the THEOREM is already machine-checked on Prove2Me (algorithm96_eq_transGen). What this mission adds is the statement in the paper's own terms: the introduction's power sum ⋁p=1dMp\bigvee_{p=1}^d M^p⋁p=1d​Mp, footnote 2's equivalence between powers and chains, which requires that chains can be shortened to length at most ddd, and footnote 3's in-place versus recursive equivalence. None of these three is on the platform. The paper itself proves only the THEOREM; footnotes 2 and 3 are stated without proof.

Difficulty

Footnote 2 calls the equivalence "trivial", but one direction needs a pigeonhole argument: a chain i→k1→⋯→kn→ji \to k_1 \to \cdots \to k_n \to ji→k1​→⋯→kn​→j of arbitrary length witnesses some power Mn+1M^{n+1}Mn+1, and the power sum stops at MdM^dMd. A chain with repeated indices has to be shortened, and the shortening must keep at least one step, since there is no reflexive part (mii′=1m'_{ii} = 1mii′​=1 only through a cycle). For a ddd-cycle the diagonal entries appear only in MdM^dMd, so the range 1..d1..d1..d is exact.

Footnote 3 compares two computations that read different data. The construction overwrites entries during a pass, and later steps of the same pass read the overwritten values; the recursion reads only the previous stage. The two agree only because of what a pass does and does not overwrite, and the published algorithm96 is three nested left folds over List.finRange, so the comparison has to be made against that concrete loop structure rather than an idealized simultaneous update.

Formalization scope

  • Matrices are Fin d → Fin d → Bool, with 111 as true. Indices are 0,…,d−10, \dots, d-10,…,d−1 instead of the page's 1,…,d1, \dots, d1,…,d; the pivot of pass iii (1-based) is the 0-based index i−1i-1i−1.
  • d=0d = 0d=0 is allowed; all statements are then trivially true, as on the page.
  • The boolean product boolProd A B i j is decide (∃ k, A i k = true ∧ B k j = true). Mathlib's matrix product on Matrix (Fin d) (Fin d) Bool is not the boolean product: Bool is a boolean ring there, its addition is exclusive or, and that product counts walks modulo 2. A formalization through * or ^ on Matrix _ _ Bool states a different, false theorem and is ruled out.
  • boolPow M p is the page's MpM^pMp for p≥1p \ge 1p≥1; the base boolPow M 0 is the identity matrix, which is not on the page and does not enter the sum. powerSum M i j holds iff boolPow M p i j for some p∈{1,…,d}p \in \{1, \dots, d\}p∈{1,…,d}. No reflexive closure is included anywhere.
  • ChainRel M i j is the existence of a list [k1,…,kn][k_1, \dots, k_n][k1​,…,kn​] (possibly empty) making i,k1,…,kn,ji, k_1, \dots, k_n, ji,k1​,…,kn​,j a List.IsChain of true entries.
  • warshallRec M n uses the 0-based pivot nnn at stage n+1n+1n+1 and is constant after stage ddd, which footnote 3 never uses.
  • M∗M^*M∗ is the published algorithm96, which fixes the inner loops over rows and columns in ascending order; Warshall leaves that order open, and it does not affect the result.
  • The goal must not be stated through TransGen, ChainRel or warshallRec, and powerSum is defined from boolPow, not from the closure or the algorithm.

Reusable beyond this mission: the boolean matrix product and powers, the power-sum characterization of the transitive closure of a relation on Fin d (footnote 2), and the in-place/recursive equivalence (footnote 3), which carries over to Floyd's Algorithm 97. Proofs of milestone 2 and of the companion are welcome independently of the goal.

Selected references

  • S. Warshall, A theorem on Boolean matrices, J. ACM 9(1) (1962) 11–12. https://doi.org/10.1145/321105.321107
  • R. W. Floyd, Algorithm 96: Ancestor, Comm. ACM 5(6) (1962) 344–345. https://doi.org/10.1145/367766.368166
  • R. W. Floyd, Algorithm 97: Shortest path, Comm. ACM 5(6) (1962) 345. https://doi.org/10.1145/367766.368168
  • R. T. Prosser, Applications of Boolean matrices to the analysis of flow diagrams, Proc. Eastern Joint Computer Conference (1959) 133–138. https://doi.org/10.1145/1460299.1460314
  • R. B. Marimont, A new method of checking the consistency of precedence matrices, J. ACM 6(2) (1959) 164–171. https://doi.org/10.1145/320964.320975
  • B. Roy, Transitivité et connexité, C. R. Acad. Sci. Paris 249 (1959) 216–218.
5 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations Research·Captain: mikedeng1

Graphs and Cooperation in Games: The Unique Fair Allocation Rule, the Shapley Value of the Graph-Restricted Game, Is Totally Stable for Superadditive GamesResearch Paper

Motivation

Classical cooperative game theory assumes that any coalition of players can form and coordinate. In many settings cooperation is mediated by pairwise relationships: communication lines, contracts, alliances, or trade links. Roger Myerson's discussion paper Graphs and Cooperation in Games (Northwestern University, 1976; published in Mathematics of Operations Research 2(3), 1977, doi:10.1287/moor.2.3.225) models such a cooperation structure as a graph on the set of players and asks how the players should share the worth of the game when only linked players can coordinate directly.

The answer, now called the Myerson value, is the starting point of the literature on communication situations and network games. It underlies later work on network formation (Jackson and Wolinsky, A strategic model of social and economic networks, JET 1996) and on allocation rules for networks, and it is a standard textbook example of an axiomatic solution concept.

Timeline:

  • 1953. Shapley defines the Shapley value φ\varphiφ for transferable-utility games by efficiency, symmetry, a carrier axiom and additivity (Shapley 1953).
  • 1976–1977. Myerson introduces cooperation graphs, the fairness (equal gains from each link) condition, proves existence and uniqueness of the fair allocation rule, identifies it as the Shapley value of the graph-restricted game, proves its total stability for superadditive games, and extends existence and uniqueness to games without transferable utility in graph function form.

Setting

Let NNN be a nonempty finite set of players and CLCLCL the set of nonempty coalitions S⊆NS\subseteq NS⊆N. A game in characteristic function form is a vector v∈RCLv\in\mathbb R^{CL}v∈RCL; vSv_SvS​ is the transferable wealth coalition SSS can divide.

A link n:mn{:}mn:m is an unordered pair of distinct players and a graph ggg is a set of links; GRGRGR is the set of all graphs, gˉN\bar g^Ngˉ​N the complete graph, g∖n:mg\setminus n{:}mg∖n:m the graph with one link removed, and ∣g∣|g|∣g∣ the number of links. Players a,ba,ba,b are connected in SSS by ggg if a=b∈Sa=b\in Sa=b∈S or a path of links of ggg joins them through players of SSS only; the classes form the partition S/gS/gS/g of SSS, and N/gN/gN/g is the set of connected components of ggg. The graph-restricted game is

(v/g)S=∑T∈S/gvT.(v/g)_S=\sum_{T\in S/g}v_T .(v/g)S​=T∈S/g∑​vT​.

An allocation rule is a function Y:GR→RNY:GR\to\mathbb R^NY:GR→RN, with Yn(g)Y_n(g)Yn​(g) the payoff of player nnn under cooperation structure ggg. It is a fair allocation rule for vvv if

  1. (efficiency, (7)) ∑n∈SYn(g)=vS\sum_{n\in S}Y_n(g)=v_S∑n∈S​Yn​(g)=vS​ for every ggg and every component S∈N/gS\in N/gS∈N/g, and
  2. (equity, (10)) Yn(g)−Yn(g∖n:m)=Ym(g)−Ym(g∖n:m)Y_n(g)-Y_n(g\setminus n{:}m)=Y_m(g)-Y_m(g\setminus n{:}m)Yn​(g)−Yn​(g∖n:m)=Ym​(g)−Ym​(g∖n:m) for every ggg and every link n:m∈gn{:}m\in gn:m∈g.

YYY is totally stable if Yn(g)≥Yn(g∖n:m)Y_n(g)\ge Y_n(g\setminus n{:}m)Yn​(g)≥Yn​(g∖n:m) for every ggg and every link n:m∈gn{:}m\in gn:m∈g, and vvv is superadditive if vS∪T≥vS+vTv_{S\cup T}\ge v_S+v_TvS∪T​≥vS​+vT​ for disjoint S,T∈CLS,T\in CLS,T∈CL.

A game in graph function form assigns to every pair (S,g)(S,g)(S,g) with S∈N/gS\in N/gS∈N/g a closed, comprehensive (closed under decreasing coordinates), proper (∅≠W≠RS\emptyset\ne W\ne\mathbb R^S∅=W=RS) set w(S,g)⊆RSw(S,g)\subseteq\mathbb R^Sw(S,g)⊆RS of feasible payoffs; its fair allocation rule replaces (7) by (14), (Yn(g))n∈S∈∂w(S,g)(Y_n(g))_{n\in S}\in\partial w(S,g)(Yn​(g))n∈S​∈∂w(S,g).

Formalization targets

Goal: Theorem 3 (p. 8)

If vvv is superadditive, then the fair allocation rule YYY for vvv is totally stable:

Yn(g)≥Yn(g∖n:m)for all g∈GR, n:m∈g.Y_n(g)\ge Y_n(g\setminus n{:}m)\qquad\text{for all } g\in GR,\ n{:}m\in g .Yn​(g)≥Yn​(g∖n:m)for all g∈GR, n:m∈g.

Theorem 1 (p. 7) and Theorem 4 (p. 10)

Every v∈RCLv\in\mathbb R^{CL}v∈RCL has a unique fair allocation rule; every game www in graph function form has a unique YYY satisfying (14) and (10).

Theorem 2 (p. 8)

The fair allocation rule is

Y(g)=φ(v/g)for all g∈GR,in particular Y(gˉN)=φ(v).Y(g)=\varphi(v/g)\quad\text{for all }g\in GR,\qquad\text{in particular } Y(\bar g^N)=\varphi(v).Y(g)=φ(v/g)for all g∈GR,in particular Y(gˉ​N)=φ(v).

The milestones follow the paper's Section 7: the steps of the proof of Theorem 4 (displays (15)–(17)), Theorem 4, Theorem 1, the decomposition of v/gv/gv/g and the efficiency and equity of φ(v/g)\varphi(v/g)φ(v/g) (proof of Theorem 2), Theorem 2, and the monotonicity of v/gv/gv/g and φ(v/g)\varphi(v/g)φ(v/g) in a link (proof of Theorem 3). Theorem 3 is the goal because its proof uses the formula of Theorem 2, whose proof uses the uniqueness of Theorem 1, which is the transferable-utility case of Theorem 4: all four numbered results lie on one path.

Significance

The theorems give an axiomatic foundation for allocation in networks: two natural requirements, component-wise efficiency and equal gains from each link, single out one rule, and that rule is an explicit formula built from the Shapley value. At the complete graph it recovers the Shapley value itself, so the result is also a new characterization of the Shapley value. Total stability says that, for superadditive games, no player has an incentive to break a link, which connects the allocation rule to the question of which networks form.

The results are proved in the paper. A search of the platform found no formalization of cooperation graphs, the graph-restricted game, the Myerson value or total stability; the Shapley value is available as a platform definition and is reused. This mission produces machine-checked statements of all four theorems and of the intermediate steps of their proofs, together with a reusable definition layer for communication situations.

Difficulty

The equity condition (10) relates the payoff at a graph to the payoff at a graph with one fewer link, and it constrains only pairs of linked players; it is not evident that these conditions, together with efficiency on components, determine the rule at every graph, nor that they are consistent. Uniqueness and existence must hold simultaneously for all 2∣gˉN∣2^{|\bar g^N|}2∣gˉ​N∣ graphs, and the constraints couple each graph to all of its subgraphs and each player to every other player of its component. In the graph function form, efficiency becomes a boundary condition on a closed comprehensive set, and the existence and uniqueness of the relevant extremal point depend on all three properties of w(S,g)w(S,g)w(S,g).

Identifying the rule with φ(v/g)\varphi(v/g)φ(v/g) requires the efficiency of the Shapley value on each component of ggg, which is not the efficiency of φ\varphiφ on NNN, and equity requires tracking how the partitions S/gS/gS/g change when a link is removed. Total stability is false for general games: without superadditivity a player can gain by cutting a link, so the hypothesis cannot be dropped.

Formalization scope

Players are Fin n, with 0<n0<n0<n in every theorem as the paper assumes NNN nonempty; the paper's player kkk is index k−1k-1k−1. The general game carrier has its own definition module, and the graph-specific definitions build on it. Graphs are SimpleGraph (Fin n): gˉN\bar g^Ngˉ​N is ⊤, the empty graph ⊥, h⊆gh\subseteq gh⊆g is h ≤ g and h⊂gh\subset gh⊂g is h < g. A game is a function Finset (Fin n) → ℝ; its value at ∅\emptyset∅ is a junk coordinate that no definition reads, and no theorem assumes v∅=0v_\emptyset=0v∅​=0. Superadditivity quantifies over nonempty coalitions only. The Shapley value is the platform definition Supermodularity.Cooperative.ShapleyValue, applied after resetting the value at ∅\emptyset∅ to 000. RS\mathbb R^SRS is ↥S → ℝ with the product topology, and ∂\partial∂ is the topological frontier.

Explicit readings of the typescript:

  • The printed definition of "connected in SSS by ggg" (p. 3) asks ni∈Sn^i\in Sni∈S only for i≥1i\ge1i≥1; the formalization requires every vertex of the path, including the start, to lie in SSS, which makes S/gS/gS/g the partition of SSS the paper describes.
  • Comprehensiveness (p. 10) prints "∀n∈N\forall n\in N∀n∈N" for coordinates of vectors in RS\mathbb R^SRS; read as n∈Sn\in Sn∈S.
  • Display (17) writes Yn(h)Y_n(h)Yn​(h) under the index m∈Sm\in Sm∈S; read as Ym(h)Y_m(h)Ym​(h), as in (18b).
  • The proof headed "PROOF OF THEOREM 4." on p. 13 is the proof of Theorem 3.
  • "The unique fair allocation rule" in Theorems 2 and 3 is stated for every fair allocation rule; existence and uniqueness are Theorem 1.
  • The paper's NNN is nonempty; every theorem carries the corresponding hypothesis 0<n0<n0<n.

The equity and stability conditions quantify over links of ggg only. Quantifying over all pairs of players would make Theorem 1 false and the goal vacuous, so that reading is ruled out. Likewise (14) is boundary membership, not membership in w(S,g)w(S,g)w(S,g).

Contributions welcome: general lemmas on the partitions S/gS/gS/g (refinement under link deletion, restriction to components), Möbius inversion over the subgraph lattice of a finite simple graph, the carrier and null-player properties of the Shapley value, and the extremal-point lemma for closed comprehensive sets. These are reusable beyond this mission.

Selected references

  • R. B. Myerson, Graphs and Cooperation in Games, Discussion Paper No. 246, CMS-EMS, Northwestern University, September 1976; published in Mathematics of Operations Research 2(3):225–229, 1977. https://doi.org/10.1287/moor.2.3.225
  • L. S. Shapley, A Value for n-Person Games, in Contributions to the Theory of Games II, Annals of Mathematics Studies 28, Princeton University Press, 1953, pp. 307–317. https://doi.org/10.1515/9781400881970-018
  • M. O. Jackson and A. Wolinsky, A Strategic Model of Social and Economic Networks, Journal of Economic Theory 71(1):44–74, 1996. https://doi.org/10.1006/jeth.1996.0108
18 thms2 active usersReviewed
🏆Completed
CombinatoricsOperations Research·Captain: mikedeng1

Zur allgemeinen Kurventheorie I: Menger's Theorem for Finite Graphs, If No Fewer Than n Points Separate Two Finite Sets P and Q, Then n Disjoint Paths Join P to QResearch Paper

Motivation

Menger's theorem is the oldest of the min–max theorems of combinatorial optimization. It says that the largest number of disjoint paths joining two sets of vertices of a finite graph equals the smallest number of vertices whose removal disconnects them. The max-flow/min-cut theorem of Ford and Fulkerson (1956), the Kőnig–Egerváry theorem on bipartite matchings (1931), Whitney's characterization of kkk-connected graphs (1932) and the integrality of network-flow polytopes all belong to its family. In operations research it underlies network reliability (how many vertex failures a network survives) and disjoint routing.

Karl Menger proved it in 1927, in a paper on curves. In Zur allgemeinen Kurventheorie (Fund. Math. 10, pp. 96–115) the graph statement is a lemma, Satz δ (p. 101), on the way to characterizing the order of a point of a regular curve by the number of arcs ending there (the subject of mission II of this series).

Timeline.

  • 1927. Menger states Satz β for compact regular one-dimensional spaces and proves its finite case, Satz δ, for finite unions of arcs, by induction on the number of arcs (Menger 1927).
  • 1931. Kőnig observes that one step of the induction, called "offenbar" on p. 102, is in the bipartite case exactly his matching theorem, and supplies it (Kőnig 1931, see the references). Kőnig's 1936 book gives the first complete graph-theoretic proof.
  • 1932. Whitney derives the two-vertex form for kkk-connected graphs.
  • 1956. Ford and Fulkerson's max-flow/min-cut theorem gives the flow-theoretic proof; Dantzig and Fulkerson connect it to LP duality.

Setting

Let GGG be a finite simple graph on a vertex set VVV, and let P,Q⊆VP, Q \subseteq VP,Q⊆V be disjoint sets of vertices.

A set S⊆VS \subseteq VS⊆V separates PPP and QQQ if every walk of GGG from a vertex of PPP to a vertex of QQQ passes through a vertex of SSS. The set SSS may contain vertices of PPP and of QQQ; in particular S=PS = PS=P always separates.

GGG is nnn-point connected between PPP and QQQ ("n-punktig zusammenhängend", p. 100) if no set of fewer than nnn vertices separates PPP and QQQ:

∣S∣≥nfor every S separating P and Q.|S| \ge n \quad\text{for every } S \text{ separating } P \text{ and } Q .∣S∣≥nfor every S separating P and Q.

A family of nnn disjoint PPP–QQQ paths consists of paths W1,…,WnW_1, \dots, W_nW1​,…,Wn​ of GGG, each from a vertex of PPP to a vertex of QQQ, whose vertex sets are pairwise disjoint, end vertices included.

The degree (Grad) of GGG is its number of edges. A part of GGG is a spanning subgraph H≤GH \le GH≤G; it is proper if it lacks an edge of GGG. GGG is irreducibly nnn-point connected between PPP and QQQ if it is nnn-point connected and no proper part is. For a separating set SSS, the side K1K_1K1​ of SSS towards PPP is formed by the components of G−SG - SG−S that meet P−SP - SP−S, together with their edges to S−PS - PS−P; the side K2K_2K2​ towards QQQ is defined symmetrically.

In Lean these are Menger27.Graphs.Separates, NPointConnected, HasDisjointPaths, IrreduciblyNPointConnected and sidePart.

Formalization targets

Goal: Satz δ (p. 101)

G is n-point connected between disjoint P and Q⟹G contains n disjoint P–Q paths.G \text{ is } n\text{-point connected between disjoint } P \text{ and } Q \quad\Longrightarrow\quad G \text{ contains } n \text{ disjoint } P\text{–}Q \text{ paths.}G is n-point connected between disjoint P and Q⟹G contains n disjoint P–Q paths.

This is the direction Menger states; the converse (disjoint paths force every separator to be large) is immediate and is not part of the mission.

Milestones: the proof of Satz δ, in the paper's order (pp. 101–102)

  1. A graph nnn-point connected between PPP and QQQ has at least nnn edges.
  2. If it has exactly nnn edges, it consists of nnn disjoint PPP–QQQ paths.
  3. It contains an irreducibly nnn-point connected part K′K'K′.
  4. If K′K'K′ is irreducible and has more than nnn edges, it has a vertex s∉P∪Qs \notin P \cup Qs∈/P∪Q on one of its edges.
  5. Through such an sss passes a set SSS of exactly nnn vertices separating PPP and QQQ in K′K'K′.
  6. For an nnn-vertex separating set SSS with ∣S∩P∣=p|S \cap P| = p∣S∩P∣=p, the side K1K_1K1​ is (n−p)(n-p)(n−p)-point connected between P−SP - SP−S and S−PS - PS−P.
  7. n−pn - pn−p disjoint paths in K1K_1K1​ and n−qn - qn−q disjoint paths in K2K_2K2​ (with q=∣S∩Q∣q = |S \cap Q|q=∣S∩Q∣) combine into nnn disjoint PPP–QQQ paths in GGG.

Significance

The result. Satz δ is the vertex form of Menger's theorem for two sets, from which the other standard forms follow: the two-vertex form, Whitney's theorem, the edge form (through line graphs), Kőnig's theorem (for bipartite GGG with PPP, QQQ the two sides) and, through a standard reduction, integral max-flow/min-cut for unit capacities. Menger used it as the finite step of his theorem on nnn-legs at points of regular curves.

Formalizing it. The theorem is proved, with many published proofs. Mathlib at the pinned revision has edge connectivity of simple graphs but no vertex separators between sets and no form of Menger's theorem; nothing on Prove2Me states it in set form. The mission produces a statement of Satz δ about Mathlib's SimpleGraph, together with the intermediate claims of Menger's own induction, so that the original argument, including the step Kőnig completed, can be checked. Related posed items on the platform are the two-vertex Whitney theorem Balinski61.Whitney.whitney_theorem (internally disjoint paths between two vertices) and Kőnig's theorem MatousekLP.Integrality.konig; neither is the set form stated here, and neither is proved yet.

Difficulty

The obvious argument, choosing disjoint paths one at a time, fails: a first path chosen badly can block every completion, so the paths must be chosen together. Menger's induction avoids this by passing to an irreducible part and cutting it at an nnn-vertex separator, but it needs a vertex outside P∪QP \cup QP∪Q to cut at (milestone 4). When every vertex of the irreducible part lies in P∪QP \cup QP∪Q the part is bipartite between PPP and QQQ, and the claim is then Kőnig's matching theorem rather than an observation; this is the gap of 1927. The gluing step (milestone 7) is also more delicate than on the page: the paths of the two sides must meet the separator only at their ends, which depends on how the sides are cut out.

Formalization scope

Representation. Menger's gewöhnlich eindimensionaler Raum, a finite union of arcs meeting at most in end points, is read as a finite simple graph G : SimpleGraph V with [Fintype V]; PPP and QQQ are Finset V. Vertices are the points of P∪QP \cup QP∪Q and the end and branch points (Menger's punktförmige Stücke); edges are the arcs between them, subdivided so that there are no loops or parallel edges. Subdivision changes neither separation nor the number of disjoint paths, and a separating point inside an arc can be moved to an end of the arc. Menger's degree (the number of arcs between consecutive punktförmige Stücke) becomes the number of edges, which a subdivision only increases; every milestone stays true.

Conventions and explicit readings of the paper's phrases.

  • Separation is phrased through walks: every walk from PPP to QQQ meets SSS. Separators may meet PPP and QQQ.
  • Disjointness of PPP and QQQ ("fremden", p. 100) is a hypothesis of every theorem, not part of the definition. Paths are disjoint including their end points; with P∩Q=∅P \cap Q = \emptysetP∩Q=∅ every path has an edge.
  • "Teil" is a spanning subgraph H ≤ G; "echter Teil" is H < G.
  • "Offenbar enthält K′ ein punktförmiges Stück s …" (p. 102) is read as: a vertex s∉P∪Qs \notin P \cup Qs∈/P∪Q lying on an edge of K′K'K′.
  • "In K′ − s sind n − 1 Punkte … so dass … S = {s, s₂, … sₙ}" is read as: a separating set of exactly nnn vertices containing sss.
  • "Besteht aus n Bögen" (milestone 2) is read as: nnn disjoint PPP–QQQ paths covering every edge.
  • The side K1K_1K1​ is constructed from GGG, PPP and SSS (sidePart G P S), not quantified over; it has no edges at the vertices of S∩PS \cap PS∩P, matching Menger's K1⊆K′−SK_1 \subseteq K' - SK1​⊆K′−S considered between P−S⋅PP - S\cdot PP−S⋅P and S−S⋅PS - S\cdot PS−S⋅P.
  • The subtractions n−pn - pn−p and n−qn - qn−q are exact, since ∣S∣=n|S| = n∣S∣=n.
  • On p. 101 the induction hypothesis is printed as "Satz γ"; it is Satz δ.

Trivializations ruled out. Forcing separators to avoid P∪QP \cup QP∪Q, counting internally disjoint or edge-disjoint paths, or allowing two paths to share an end vertex would each change the theorem; none is done. No hypothesis restricts nnn, and n=0n = 0n=0 is the trivial case, not an excluded one.

Needed infrastructure. Walk surgery (prefixes up to the first visit of a set, concatenation, conversion of walks to paths), transfer of walks along SimpleGraph inclusions, and well-founded induction on the number of edges. Separators and disjoint path families between sets are reusable beyond this mission, for Whitney's theorem, Kőnig's theorem and flow integrality. Proofs by any method are welcome, including a proof of the goal through max-flow/min-cut or an independent argument (Göring's) bypassing the milestones.

Selected references

  • K. Menger, Zur allgemeinen Kurventheorie, Fundamenta Mathematicae 10 (1927), 96–115. https://doi.org/10.4064/fm-10-1-96-115
  • D. Kőnig, Gráfok és mátrixok, Matematikai és Fizikai Lapok 38 (1931), 116–119 (no DOI; English discussion in Diestel, Graph Theory, §2.1 and §3.3).
  • H. Whitney, Congruent graphs and the connectivity of graphs, American Journal of Mathematics 54 (1932), 150–168. https://doi.org/10.2307/2371086
  • L. R. Ford, D. R. Fulkerson, Maximal flow through a network, Canadian Journal of Mathematics 8 (1956), 399–404. https://doi.org/10.4153/CJM-1956-045-5
  • F. Göring, Short proof of Menger's theorem, Discrete Mathematics 219 (2000), 295–296. https://doi.org/10.1016/S0012-365X(00)00088-1
  • R. Diestel, Graph Theory, 5th ed., Springer, 2017, §3.3. https://doi.org/10.1007/978-3-662-53622-3
10 thms2 active usersReviewed
🏆Completed
CombinatoricsLinear OptimizationOperations Research·Captain: mikedeng1

Integer Programming Formulation of Traveling Salesman Problems: The Miller–Tucker–Zemlin Constraints Exactly Describe the Itineraries with at Most p Cities per TourResearch Paper

Why this formulation matters

The traveling salesman problem asks for a shortest closed route through a set of cities. Writing it as an integer program requires linear constraints whose integer solutions are exactly the routes. The assignment constraints (each city entered once and left once) are not enough: their integer solutions include collections of disjoint cycles, called subtours, that never pass through the starting city. Some extra family of constraints must rule these out.

In 1960 C. E. Miller, A. W. Tucker and R. A. Zemlin gave such a family using one auxiliary real variable per city (J. ACM 7(4), 326–329). Their MTZ constraints have polynomially many rows, n2−nn^2-nn2−n for nnn cities, whereas the subtour elimination constraints of Dantzig, Fulkerson and Johnson (1954) (Oper. Res. 2(4)) have one row per subset of cities. They also handle a generalization: a salesman who returns to base several times, visiting at most ppp cities per trip. This is the basic form of a capacitated vehicle routing problem. MTZ-type constraints remain standard in vehicle routing models and in textbook formulations of the TSP. Their linear relaxations are known to be weaker than those of the subtour elimination formulation (Padberg and Sung, 1991).

The paper attributes the model to Tucker and reports five machine experiments with Gomory's all-integer algorithm. The computational part is outside this mission.

Setting

There are nnn cities 1,…,n1,\dots,n1,…,n and a base city 000. A tour is a succession of visits to cities without stopping at city 000. In problem (1) the salesman leaves 000, visits each of the nnn cities exactly once and returns to 000, possibly several times, with at most ppp cities per tour. Formally, an itinerary is a list I=(T1,…,Tt)I=(T_1,\dots,T_t)I=(T1​,…,Tt​) of tours. Each tour is a nonempty list of cities other than 000 of length at most ppp, and each city 1,…,n1,\dots,n1,…,n occurs exactly once in the concatenation of the tours. The number ttt of tours is the number of returns to 000. Tour (c1,…,ck)(c_1,\dots,c_k)(c1​,…,ck​) travels the arcs (0,c1),(c1,c2),…,(ck,0)(0,c_1),(c_1,c_2),\dots,(c_k,0)(0,c1​),(c1​,c2​),…,(ck​,0). aij(I)a_{ij}(I)aij​(I) denotes the number of times III travels the arc (i,j)(i,j)(i,j). For real distances dijd_{ij}dij​, the length of III is the sum of dijd_{ij}dij​ over the arcs it travels.

Problem (2) has non-negative integer variables xijx_{ij}xij​ (0≤i≠j≤n0\le i\ne j\le n0≤i=j≤n) and real variables u1,…,unu_1,\dots,u_nu1​,…,un​:

min⁡∑∑0≤i≠j≤ndijxijs.t.∑i=0i≠jnxij=1 (j≥1),∑j=0j≠inxij=1 (i≥1),ui−uj+p xij≤p−1 (1≤i≠j≤n).\min \sum\sum_{0\le i\ne j\le n} d_{ij}x_{ij}\quad\text{s.t.}\quad \sum_{\substack{i=0\\ i\ne j}}^{n} x_{ij}=1\ (j\ge1),\quad \sum_{\substack{j=0\\ j\ne i}}^{n} x_{ij}=1\ (i\ge1),\quad u_i-u_j+p\,x_{ij}\le p-1\ (1\le i\ne j\le n).min∑0≤i=j≤n∑​dij​xij​s.t.i=0i=j​∑n​xij​=1 (j≥1),j=0j=i​∑n​xij​=1 (i≥1),ui​−uj​+pxij​≤p−1 (1≤i=j≤n).

When the number of tours ttt is fixed, the relation ∑i=1nxi0=t\sum_{i=1}^n x_{i0}=t∑i=1n​xi0​=t is added. With t=1t=1t=1 and p≥np\ge np≥n, (1) is the standard traveling salesman problem.

Formalization targets

Goal: the equivalence of (1) and (2)

For every n≥0n\ge0n≥0 and every p≥1p\ge1p≥1:

{ x:∃u, (x,u) feasible for (2) }  =  { a(I):I a legitimate itinerary of (1) },\{\,x : \exists u,\ (x,u)\ \text{feasible for (2)}\,\} \;=\; \{\,a(I) : I \text{ a legitimate itinerary of (1)}\,\},{x:∃u, (x,u) feasible for (2)}={a(I):I a legitimate itinerary of (1)},

with ∑i=1nxi0\sum_{i=1}^n x_{i0}∑i=1n​xi0​ equal to the number of tours t(I)t(I)t(I), and with the uiu_iui​ in the converse taken to be non-negative integers. In addition, ∑∑i≠jdij aij(I)\sum\sum_{i\ne j} d_{ij}\,a_{ij}(I)∑∑i=j​dij​aij​(I) equals the length of III for all real ddd. The paper says that (2) "will be shown to be equivalent to (1)". Its "burden of proof" is that the two feasible sets correspond and the two objectives agree, and the goal states exactly that.

Milestones, in the order of the paper's proof

  1. The constraints of (2) force xij∈{0,1}x_{ij}\in\{0,1\}xij​∈{0,1} (p. 327).
  2. Along a path r0,…,rmr_0,\dots,r_mr0​,…,rm​ of non-base cities with all xrkrk+1=1x_{r_kr_{k+1}}=1xrk​rk+1​​=1, ur0−urm≤−mu_{r_0}-u_{r_m}\le-mur0​​−urm​​≤−m (p. 328).
  3. A feasible xxx has no cycle that avoids city 000, so all tours include city 0 (pp. 327–328).
  4. For p≥1p\ge1p≥1, no tour of a feasible xxx visits more than ppp cities (p. 328).
  5. The converse: ui=u_i=ui​= the position of city iii in its tour makes (a(I),u)(a(I),u)(a(I),u) feasible, with 1≤ui≤p1\le u_i\le p1≤ui​≤p (p. 328).
  6. Under x=a(I)x=a(I)x=a(I) the objective of (2) is the distance travelled in (1) (p. 327).
  7. An itinerary with ttt tours of at most ppp cities requires tp≥ntp\ge ntp≥n (p. 326).

Significance

The result. The theorem justifies using (2) in place of (1). Because the feasible sets match and the objectives agree, an optimal solution of (2) is an optimal itinerary of (1), and the same holds with the number of tours fixed. It also shows that a compact, polynomial-size family of linear inequalities in auxiliary potentials excludes every subtour. The same idea, with the uiu_iui​ read as loads or arrival times, underlies the standard capacity and time-window constraints of vehicle routing.

Formalizing it. The result is classical and its proof is short, but the paper argues informally. It follows successor chains, treats "a tour" without defining an itinerary, and leaves several steps implicit, including the hypothesis p≥1p\ge1p≥1 and the distinctness of rp+1r_{p+1}rp+1​ and r1r_1r1​. No machine-checked statement of the MTZ equivalence is recorded on the platform. The remaining work is the whole formal proof: the combinatorics that turns a 0/1 matrix with unit row and column sums and no subtours into a list of tours, the converse labelling, and the bookkeeping identity between the two objectives.

Difficulty

The inequalities themselves are elementary. Each milestone about (2) is a short computation with the potentials uiu_iui​. The work is in the passage from a feasible matrix xxx to an itinerary. The successor walk from city 000 must be shown to terminate at 000 and to cover every arc of xxx. The resulting tours must be listed so that their arc counts reproduce xxx exactly, including the number of arcs leaving 000. That number is not constrained directly in (2). It equals ∑ixi0\sum_{i}x_{i0}∑i​xi0​ only through a counting argument on the row and column sums. A proof that only shows some itinerary exists does not establish the correspondence the paper claims.

Formalization scope

Cities are Fin (n + 1) with 0 the base city, which is the paper's indexing. The commitments, each also stated in the items' Formalization Notes, are:

  • An itinerary is a List of tours, each tour a List of cities. The paper's "legitimate itinerary" is read as: every tour nonempty, free of 0, of length at most ppp; the concatenation duplicate-free; every city 1,…,n1,\dots,n1,…,n present.
  • "A feasible solution to (2) has xijx_{ij}xij​ which do define a legitimate itinerary" is read as x=a(I)x=a(I)x=a(I), the arc counts of III, together with ∑i=1nxi0=t(I)\sum_{i=1}^n x_{i0}=t(I)∑i=1n​xi0​=t(I) ("The number of returns to city 0 is given by ∑xi0\sum x_{i0}∑xi0​"). Arc counts are counts, so xij∈{0,1}x_{ij}\in\{0,1\}xij​∈{0,1} is a theorem, not an encoding choice.
  • "Together with appropriate uiu_iui​" is made explicit as the existence of natural-number uiu_iui​, following the paper's remark that the uiu_iui​ may be taken non-negative integers. The labelling ui=ju_i=jui​=j is a definition.
  • Problem (2) has no variables xiix_{ii}xii​. The matrix xxx is indexed by all pairs and the clause xii=0x_{ii}=0xii​=0 encodes their absence. u0u_0u0​ exists in the encoding but appears in no constraint.
  • xij∈Nx_{ij}\in\mathbb Nxij​∈N ("non-negative integers"), ui∈Ru_i\in\mathbb Rui​∈R ("arbitrary real numbers"), and p−1p-1p−1 and p xijp\,x_{ij}pxij​ are computed in R\mathbb RR.
  • Added hypothesis p≥1p\ge1p≥1 in the goal and in milestone 4. The paper does not state it, and the claim is false without it: for n=1n=1n=1, p=0p=0p=0, x01=x10=1x_{01}=x_{10}=1x01​=x10​=1 is feasible for (2) but no itinerary exists.
  • The distances dijd_{ij}dij​ are arbitrary reals, with no sign, symmetry or triangle condition, as on the page.
  • Milestone 2 states the exact summed bound −m-m−m. The page prints the weaker j+1−kj+1-kj+1−k.
  • Milestone 3 excludes every cycle of xxx avoiding 000. The page argues for the cycle met while following successors from 000, by an argument that applies to any such cycle.

Not formalized: the constraint and variable counts and the elimination of xi0,x0jx_{i0},x_{0j}xi0​,x0j​ "to produce an equivalent problem with n2+nn^2+nn2+n inequalities and n2n^2n2 variables" (p. 328), which the paper does not make precise, and the machine experiments (pp. 328–329).

Trivializing readings are ruled out. The goal identifies xxx with the arc counts of an itinerary, not merely the existence of some itinerary. The hypotheses are satisfiable: for n=3n=3n=3, p=2p=2p=2, the itinerary [[1,2],[3]][[1,2],[3]][[1,2],[3]] with u=(1,2,1)u=(1,2,1)u=(1,2,1) is feasible for (2), while the subtour x01=x10=x23=x32=1x_{01}=x_{10}=x_{23}=x_{32}=1x01​=x10​=x23​=x32​=1 is not, for any uuu. Contributions welcome: proofs of the milestones, a reusable library of "functional digraph with unit in- and out-degree" decompositions, and the optimization corollary (equal optimal values of (1) and (2)).

Selected references

  • C. E. Miller, A. W. Tucker, R. A. Zemlin, Integer Programming Formulation of Traveling Salesman Problems, Journal of the ACM 7(4), 326–329, 1960. https://doi.org/10.1145/321043.321046
  • G. Dantzig, R. Fulkerson, S. Johnson, Solution of a Large-Scale Traveling-Salesman Problem, Operations Research 2(4), 393–410, 1954. https://doi.org/10.1287/opre.2.4.393
  • M. Padberg, T.-Y. Sung, An analytical comparison of different formulations of the travelling salesman problem, Mathematical Programming 52, 315–357, 1991. https://doi.org/10.1007/BF01582894
  • R. E. Gomory, Outline of an algorithm for integer solutions to linear programs, Bulletin of the AMS 64(5), 275–278, 1958. https://doi.org/10.1090/S0002-9904-1958-10224-4
9 thms2 active usersReviewed
PreviousPage 4 of 6Next

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