Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Theoretical Computer Science

464 missions · 238 completed

The mathematical foundations of computation: which problems can be solved, by what algorithms, and at what cost in time, space, or communication. Distinguished by its emphasis on rigor and unconditional lower bounds, it spans computational complexity, algorithm design, automata and computability, cryptography, and the analysis of Boolean functions.

Missions

Open226Completed238All464
Complexity TheoryOperations Research·Captain: wurtle

A Polynomial-Time Algorithm for Three-Machine Unit-Job SchedulingResearch Paper

Motivation

A set of unit-length jobs with precedence constraints is to be run on a fixed number of identical machines so that all jobs finish as early as possible. For two machines this was solved in 1969–1972; for a number of machines given as part of the input the problem is NP-complete (Ullman, 1975). Whether the problem is polynomial-time solvable for three machines, written P3∣prec,pj=1∣Cmax⁡P3\mid\mathrm{prec},p_j=1\mid C_{\max}P3∣prec,pj​=1∣Cmax​, is problem OPEN8 in the appendix of Garey and Johnson's 1979 book and is one of the best-known open questions on the complexity of scheduling.

Timeline

  • 1969–1971. Fujii, Kasami and Ninomiya solve the two-processor case by matching (doi:10.1137/0117070, erratum doi:10.1137/0120018).
  • 1972. Coffman and Graham give an optimal list-scheduling algorithm for two processors (doi:10.1007/BF00288685).
  • 1975. Ullman proves NP-completeness when the number of machines is part of the input (doi:10.1016/S0022-0000(75)80008-0).
  • 1979. Garey and Johnson list the fixed-three-machine problem as OPEN8; Graham, Lawler, Lenstra and Rinnooy Kan introduce the standard notation (doi:10.1016/S0167-5060(08)70356-X).
  • 1983. Garey, Johnson, Tarjan and Yannakakis give a linear-time algorithm on three processors for opposing forests (doi:10.1137/0604011).
  • 1984. Dolev and Warmuth give an O(nh(m−1)+1)O(n^{h(m-1)+1})O(nh(m−1)+1) algorithm for precedence graphs of height hhh (doi:10.1016/0196-6774(84)90039-7).
  • 2016–2022. Approximation schemes for fixed machine counts: Levey and Rothvoß via LP hierarchies (doi:10.1145/2897518.2897532), Garg's quasi-PTAS (doi:10.4230/LIPIcs.ICALP.2018.59), Li (doi:10.1137/1.9781611976465.178), Das and Wiese (doi:10.4230/LIPIcs.ESA.2022.40).
  • 2025. Nederlof, Swennenhuis and Węgrzycki give an exact algorithm in time 2O(nlog⁡n)2^{O(\sqrt n\log n)}2O(n​logn) for three machines (doi:10.1137/1.9781611978322.16).

The source of this mission, an OpenAI preprint dated September 24, 2026, claims a deterministic polynomial-time algorithm.

Setting

The input is a directed acyclic graph G=(V,E)G=(V,E)G=(V,E) on V={1,…,n}V=\{1,\dots,n\}V={1,…,n}, n≥1n\ge1n≥1, given as an explicit list of arcs, and optionally an integer deadline TTT with 1≤T≤n1\le T\le n1≤T≤n. A feasible schedule with makespan at most TTT is a map

τ:V→{1,…,T}\tau:V\to\{1,\dots,T\}τ:V→{1,…,T}

such that each value is taken at most three times (three machines) and τ(u)<τ(v)\tau(u)<\tau(v)τ(u)<τ(v) for every arc (u,v)∈E(u,v)\in E(u,v)∈E (precedence). Slot ttt means execution during [t−1,t)[t-1,t)[t−1,t). There are no release dates, communication delays or other resources. The input length LLL is the length of the binary encoding of the graph and deadline.

Formalization targets

Goal: Theorem 1.1

There is a single deterministic multitape Turing machine MMM and a constant CCC such that, on every such input of length LLL, MMM halts within

C (L+2)150020C\,(L+2)^{150020}C(L+2)150020

steps and

  • without a deadline, outputs a feasible schedule of minimum makespan;
  • with a deadline TTT, outputs "infeasible" exactly when no feasible schedule with makespan at most TTT exists, and otherwise outputs one.

The exponent is the one stated in the paper and is not optimized; the content of the theorem is that a fixed polynomial bound exists. The Lean statement OAI.ThreeMachine.main_theorem is open on the platform.

Significance

The theorem resolves the polynomial-time side of Garey and Johnson's OPEN8: unit-job makespan scheduling with arbitrary precedence constraints on three identical machines lies in P. Previous exact algorithms were either restricted to special precedence structures (forests, bounded height) or subexponential, and approximation schemes did not decide the optimum. The structural statement behind the algorithm, that feasible schedules admit recursive decompositions whose interval job sets have descriptions of bounded size, may be of independent use for other fixed-machine scheduling problems. No practical running time is claimed.

The result is claimed in an OpenAI preprint; it has not been peer reviewed and no machine-checked proof exists. A formal proof would also certify the running-time analysis in an explicit machine model, which is rarely done for algorithmic results of this size.

Difficulty

Removing selected slots from a feasible schedule leaves gaps that can be rescheduled independently, which suggests a recursive search over the job sets of gaps. The obstacle is their number: naming a gap by its two boundary triples does not determine its job set, and intersecting descriptions inherited from ancestors can make the descriptions grow without bound along a recursion of linear depth. The subexponential decomposition of Nederlof–Swennenhuis–Węgrzycki controls the number of subproblems only up to 2O(nlog⁡n)2^{O(\sqrt n\log n)}2O(n​logn); a polynomial bound needs every subproblem to have a global description of constant size.

Formalization scope

  • Instance n is a list of arcs on Fin n; acyclic forbids a Relation.TransGen cycle. Jobs are Fin n (shifted to 1,…,n1,\dots,n1,…,n in the encoding).
  • Feasible G T τ: slots in [1,T][1,T][1,T], at most three jobs per slot, strict precedence.
  • CorrectOutput distinguishes the optimization mode (none) and the decision mode (some T).
  • The machine model is a custom multitape Turing machine Machine k q g over Fin (g+3) with k+1k+1k+1 tapes; input on tape 0 encoded by encodeInput (self-delimiting binary numbers), output read from tape 0 after halting. The time bound is C (L+2)150020C\,(L+2)^{150020}C(L+2)150020 with LLL the encoded length.
  • The quantifiers are ordered so that the machine and constant are fixed before the instance (uniformity).

A complete development needs the combinatorics of the decomposition (separators and bounded global descriptions), a dynamic program over polynomially many descriptions, and its implementation and time analysis on the given Turing machine model. Contributions formalizing the purely combinatorial structure theorem, independently of the machine model, are welcome.

Selected references

  • OpenAI, A Polynomial-Time Algorithm for Three-Machine Unit-Job Scheduling, preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/A-polynomial-time-algorithm-for-three-machine-unit-job-scheduling-September-24-2026/paper.pdf
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979.
  • E. G. Coffman Jr., R. L. Graham, Optimal scheduling for two-processor systems, Acta Inform., 1972. https://doi.org/10.1007/BF00288685
  • J. D. Ullman, NP-complete scheduling problems, J. Comput. System Sci., 1975. https://doi.org/10.1016/S0022-0000(75)80008-0
  • D. Dolev, M. K. Warmuth, Scheduling precedence graphs of bounded height, J. Algorithms, 1984. https://doi.org/10.1016/0196-6774(84)90039-7
  • J. Nederlof, C. M. F. Swennenhuis, K. Węgrzycki, A subexponential time algorithm for makespan scheduling of unit jobs with precedence constraints, SODA 2025. https://doi.org/10.1137/1.9781611978322.16
  • E. Levey, T. Rothvoß, A (1+ϵ)(1+\epsilon)(1+ϵ)-approximation for makespan scheduling with precedence constraints using LP hierarchies, STOC 2016. https://doi.org/10.1145/2897518.2897532
2 thms1 active userReviewed
Information TheoryProbability·Captain: wurtle

Quantitative lower bounds for trace reconstructionResearch Paper

Motivation: how many deletion traces does it take to recover a word?

In the trace reconstruction problem an unknown binary word x∈{0,1}nx\in\{0,1\}^nx∈{0,1}n is passed repeatedly through a deletion channel: each bit is deleted independently with known probability q∈(0,1)q\in(0,1)q∈(0,1) and the surviving bits are concatenated in order. Each output, called a trace, reveals its length and its bits but not which positions survived. The question is how many independent traces are needed to recover every word xxx exactly with a fixed success probability. The problem arises in DNA storage and sequencing, in multiple sequence alignment, and as a basic test of how much information a channel with synchronization errors destroys. Its central open question has been whether a polynomial number of traces suffices at a fixed deletion probability.

Timeline

  • 2001 — Levenshtein studies efficient reconstruction of sequences from corrupted copies, in combinatorial and memoryless-channel versions (IEEE Trans. Inf. Theory 2001).
  • 2004 — Batu, Kannan, Khanna and McGregor introduce the independent-deletion trace problem in connection with multiple sequence alignment and outline a linear lower bound at fixed deletion probability (SODA 2004).
  • 2008 — Holenstein, Mitzenmacher, Panigrahy and Wieder give the first worst-case upper bound exp⁡(O~(n))\exp(\widetilde O(\sqrt n))exp(O(n​)) at constant deletion probability (SODA 2008).
  • 2014 — McGregor, Price and Vorotnikova analyze lower bounds via a single displaced one among zeros (ESA 2014).
  • 2017 — De, O'Donnell and Servedio (STOC 2017) and, independently, Nazarov and Peres (STOC 2017) improve the upper bound to exp⁡(O(n1/3))\exp(O(n^{1/3}))exp(O(n1/3)).
  • 2020 — Holden and Lyons prove the lower bound Ω(n5/4/log⁡n)\Omega(n^{5/4}/\sqrt{\log n})Ω(n5/4/logn​) (Ann. Appl. Probab. 2020; erratum 2022).
  • 2021 — Chase improves the lower bound to Ω(n3/2/log⁡7n)\Omega(n^{3/2}/\log^7 n)Ω(n3/2/log7n) (Ann. Inst. Henri Poincaré Probab. Stat. 2021) and the upper bound to exp⁡(O(n1/5log⁡5n))\exp(O(n^{1/5}\log^5 n))exp(O(n1/5log5n)) (STOC 2021).
  • 2026 — Burudgunte, Valiant and Wang give a quasipolynomial sample upper bound exp⁡(p−7/3(log⁡2n)C0)\exp(p^{-7/3}(\log_2 n)^{C_0})exp(p−7/3(log2​n)C0​) with p=1−qp=1-qp=1−q (arXiv:2607.04073).
  • 2026 — An OpenAI preprint, Quantitative lower bounds for trace reconstruction (OpenAI Math Release, September 24, 2026), claims that nΩ(log⁡log⁡n)n^{\Omega(\log\log n)}nΩ(loglogn) traces are necessary at every fixed qqq, answering the polynomial-sample question negatively. The preprint has not been peer reviewed and its theorems are not formally verified.

Setting

Fix nnn and q∈(0,1)q\in(0,1)q∈(0,1). A deletion mask a∈{0,1}na\in\{0,1\}^na∈{0,1}n keeps bit iii when ai=1a_i=1ai​=1; it has probability ∏i(1−q)aiq1−ai\prod_i (1-q)^{a_i} q^{1-a_i}∏i​(1−q)ai​q1−ai​. The trace of xxx under aaa is the list of the bits xix_ixi​ with ai=1a_i=1ai​=1, in order. Write Dq(x)\mathcal D_q(x)Dq​(x) for the law of the trace of xxx under a random mask.

An estimator with budget mmm receives mmm independent traces of xxx and returns a random word; it may be randomized and its computation is unrestricted. For s∈(0,1]s\in(0,1]s∈(0,1], the sample complexity Tq,s(n)T_{q,s}(n)Tq,s​(n) is the least mmm such that some estimator recovers every x∈{0,1}nx\in\{0,1\}^nx∈{0,1}n with probability at least sss; it is +∞+\infty+∞ if no finite mmm works. The one-trace distance is

TV⁡(Dq(x),Dq(y))=12∑w∣Dq(x)(w)−Dq(y)(w)∣,\operatorname{TV}\bigl(\mathcal D_q(x),\mathcal D_q(y)\bigr)=\tfrac12\sum_{w}\bigl|\mathcal D_q(x)(w)-\mathcal D_q(y)(w)\bigr|,TV(Dq​(x),Dq​(y))=21​w∑​​Dq​(x)(w)−Dq​(y)(w)​,

the sum running over all finite binary words.

In Lean (namespace OAI.TraceReconstruction), words are Fin n → Bool, traceMass q x w is the probability of output w, an estimator is a function A : (Fin m → Trace) → Word n → ℝ giving a probability vector for each sample, sampleComplexity n q s : ℝ≥0∞ is the infimum of sufficient budgets, and minTraceTV n q is the minimum of traceTV q x y over x≠yx\neq yx=y.

Formalization targets

Goal: quantitative sample lower bound (Theorem 1.1)

Fix s∈(0,1]s\in(0,1]s∈(0,1] and 0<c<1/(4log⁡2)0<c<1/(4\log 2)0<c<1/(4log2). For every sequence of instances (nk,qk)(n_k,q_k)(nk​,qk​) with nk→∞n_k\to\inftynk​→∞, qk∈(0,1)q_k\in(0,1)qk​∈(0,1) and qk3log⁡nk→∞q_k^3\log n_k\to\inftyqk3​lognk​→∞,

Tqk,s(nk)  ≥  nk clog⁡(qk3log⁡nk)for all sufficiently large k.T_{q_k,s}(n_k)\;\ge\; n_k^{\,c\log(q_k^3\log n_k)}\qquad\text{for all sufficiently large }k.Tqk​,s​(nk​)≥nkclog(qk3​lognk​)​for all sufficiently large k.

At fixed qqq this gives Tq,s(n)≥nΩ(log⁡log⁡n)T_{q,s}(n)\ge n^{\Omega(\log\log n)}Tq,s​(n)≥nΩ(loglogn), so no polynomial number of traces suffices. The goal is published on the platform with status Open.

Milestones: superpolynomial indistinguishability (Theorem 1.2)

For every fixed q∈(0,1)q\in(0,1)q∈(0,1) and every A>0A>0A>0,

lim⁡n→∞nAmin⁡x≠y∈{0,1}nTV⁡(Dq(x),Dq(y))=0,lim⁡n→∞Tq,s(n)nA=∞(s∈(0,1]).\lim_{n\to\infty} n^{A}\min_{x\neq y\in\{0,1\}^n}\operatorname{TV}\bigl(\mathcal D_q(x),\mathcal D_q(y)\bigr)=0, \qquad \lim_{n\to\infty}\frac{T_{q,s}(n)}{n^{A}}=\infty\quad (s\in(0,1]).n→∞lim​nAx=y∈{0,1}nmin​TV(Dq​(x),Dq​(y))=0,n→∞lim​nATq,s​(n)​=∞(s∈(0,1]).

Significance

The result itself. Theorem 1.1 is the first superpolynomial lower bound for exact worst-case trace reconstruction at a fixed deletion probability; previous lower bounds were of order n3/2n^{3/2}n3/2 up to logarithms. It holds against unrestricted computation and any fixed positive success probability, and it settles in the negative the question of whether polynomially many traces suffice. Combined with the 2026 quasipolynomial upper bound of Burudgunte, Valiant and Wang, the sample complexity at fixed qqq is now known to lie between exp⁡(Θ(log⁡nlog⁡log⁡n))\exp(\Theta(\log n\log\log n))exp(Θ(lognloglogn)) and exp⁡((log⁡n)O(1))\exp((\log n)^{O(1)})exp((logn)O(1)). Theorem 1.2 isolates the underlying phenomenon: some pairs of distinct words have single-trace laws closer than every inverse polynomial.

Formalizing it. The definitions fix exactly what an estimator may see (only the complete traces, no positions or boundaries) and what "sample complexity" means, including the value +∞+\infty+∞. A formal proof would certify a lower bound whose argument is a long multi-scale construction with many quantitative error budgets, the kind of argument where an overlooked loss would invalidate the conclusion.

Difficulty

Earlier lower bounds compare two words that differ by a single local defect and compute their trace distance directly; such pairs are distinguishable from polynomially many traces, so no single-defect construction can give a superpolynomial bound. A superpolynomial bound needs pairs of words whose trace laws agree to every polynomial order, while the argument must remain uniform in the deletion probability (which may tend to 000 slowly), keep all constants below scale margins that shrink as the number of construction stages grows with nnn, and transfer a bound for one trace to a bound for mmm independent traces and for arbitrary randomized estimators.

Formalization scope

  • Traces are List Bool; traceTV sums over the finite set of all lists of length at most nnn, which contains the support of every trace law, so it agrees with the source's total variation over all finite words.
  • sampleComplexity is an sInf in ℝ≥0∞; the infimum of the empty set is ⊤\top⊤, matching Tq,s(n)=+∞T_{q,s}(n)=+\inftyTq,s​(n)=+∞. The lower bound is compared after ENNReal.ofReal, and the right-hand side nclog⁡(q3log⁡n)n^{c\log(q^3\log n)}nclog(q3logn) uses the real power.
  • Theorem 1.1 is encoded with sequences n q : ℕ → …, Tendsto n atTop atTop, and Tendsto (q k ^ 3 * log (n k)) atTop atTop; the conclusion is ∀ᶠ k in atTop. Logarithms are natural.
  • Estimators are arbitrary real-valued probability kernels from samples to words; success is the exact probability of returning the true word, required to be at least sss for every input.
  • minTraceTV n q is a real sInf, attained for n≥1n\ge1n≥1; at n=0n=0n=0 the set is empty and the junk value 000 does not affect the limit.
  • Needed infrastructure: the deletion channel and its trace law, total variation and Hellinger distance for finite laws and their tensorization, and the reduction from sample complexity to pairwise distinguishability. These are reusable for other deletion-channel problems.

Selected references

  • V. I. Levenshtein, Efficient reconstruction of sequences, IEEE Trans. Inform. Theory, 2001. https://doi.org/10.1109/18.904499
  • T. Batu, S. Kannan, S. Khanna and A. McGregor, Reconstructing strings from random traces, SODA 2004. https://people.cs.umass.edu/~mcgregor/papers/04-soda.pdf
  • T. Holenstein, M. Mitzenmacher, R. Panigrahy and U. Wieder, Trace reconstruction with constant deletion probability and related results, SODA 2008. https://www.eecs.harvard.edu/~michaelm/postscripts/soda2008c.pdf
  • A. McGregor, E. Price and S. Vorotnikova, Trace reconstruction revisited, ESA 2014. https://doi.org/10.1007/978-3-662-44777-2_57
  • A. De, R. O'Donnell and R. A. Servedio, Optimal mean-based algorithms for trace reconstruction, STOC 2017. https://doi.org/10.1145/3055399.3055450
  • F. Nazarov and Y. Peres, Trace reconstruction with exp⁡(O(n1/3))\exp(O(n^{1/3}))exp(O(n1/3)) samples, STOC 2017. https://doi.org/10.1145/3055399.3055494
  • N. Holden and R. Lyons, Lower bounds for trace reconstruction, Ann. Appl. Probab., 2020. https://doi.org/10.1214/19-AAP1506
  • Z. Chase, New lower bounds for trace reconstruction, Ann. Inst. Henri Poincaré Probab. Stat. 57 (2021), 627–643. https://doi.org/10.1214/20-AIHP1089
  • Z. Chase, Separating words and trace reconstruction, STOC 2021. https://doi.org/10.1145/3406325.3451118
  • A. Burudgunte, P. Valiant and H. Wang, Quasipolynomial trace reconstruction, preprint, 2026. https://arxiv.org/abs/2607.04073
  • OpenAI, Quantitative lower bounds for trace reconstruction, OpenAI Math Release preprint, September 24, 2026 (Theorem 1.1, p. 1; Theorem 1.2, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/quantitative-lower-bounds-for-trace-reconstruction-September-24-2026/paper.pdf
4 thms1 active userReviewed
Discrete GeometryOptimization·Captain: wurtle

Near-square-root logarithmic integrality gaps for uniform sparsest cutResearch Paper

Motivation

Sparsest cut asks for a partition of a weighted graph whose crossing capacity is small compared with the demand it separates. It is a basic graph-partitioning primitive (used for divide-and-conquer algorithms, clustering and expansion testing), and its relaxations are the main bridge between approximation algorithms and the geometry of finite metric spaces. In the uniform version every pair of vertices has demand one. The best known polynomial-time approximation, due to Arora, Rao and Vazirani, uses the Goemans–Linial semidefinite relaxation and achieves ratio O(log⁡n)O(\sqrt{\log n})O(logn​). How good this relaxation can be, i.e. its integrality gap on uniform instances, has been a long-standing question.

Timeline

  • 1995. Linial, London and Rabinovich develop the metric approach to cuts: ℓ1\ell_1ℓ1​ metrics are combinations of cut metrics, so embeddings into ℓ1\ell_1ℓ1​ round to cuts (doi:10.1007/BF01200757).
  • 1997 and 2002. Goemans and Linial propose the semidefinite relaxation with triangle inequalities and the conjecture that negative-type metrics embed into ℓ1\ell_1ℓ1​ with constant distortion (doi:10.1007/BF02614315).
  • 2003 and 2008. Rabinovich relates uniform sparsest cut to average distortion of embeddings into the line and into L1L_1L1​ (doi:10.1145/780542.780609, doi:10.1007/s00454-007-9047-5).
  • 2004/2009. Arora, Rao and Vazirani prove the O(log⁡n)O(\sqrt{\log n})O(logn​) upper bound for uniform demands (doi:10.1145/1502793.1502794).
  • 2005/2015. Khot and Vishnoi disprove the Goemans–Linial conjecture for general demands (doi:10.1145/2629614).
  • 2006. Devanur, Khot, Saket and Vishnoi give an Ω(log⁡log⁡n)\Omega(\log\log n)Ω(loglogn) uniform integrality gap (doi:10.1145/1132516.1132594). Lee and Naor propose the Heisenberg-group approach for general demands (doi:10.1109/FOCS.2006.47).
  • 2009–2011. Cheeger, Kleiner and Naor prove a (log⁡n)Ω(1)(\log n)^{\Omega(1)}(logn)Ω(1) general-demand gap (arXiv:0910.2024) and note that their examples cannot give a diverging uniform gap.
  • 2013. Kane and Meka improve the uniform gap to exp⁡(Ω(log⁡log⁡n))\exp(\Omega(\sqrt{\log\log n}))exp(Ω(loglogn​)) (doi:10.1145/2488608.2488610).
  • 2018 and 2025. Naor and Young obtain the sharp general-demand lower bound Ω(log⁡n)\Omega(\sqrt{\log n})Ω(logn​) (doi:10.4007/annals.2018.188.1.4); Chang, Naor and Ren prove the matching upper bound (doi:10.1145/3717823.3718285).

For uniform demands the gap between the lower bound exp⁡(Ω(log⁡log⁡n))\exp(\Omega(\sqrt{\log\log n}))exp(Ω(loglogn​)) and the upper bound O(log⁡n)O(\sqrt{\log n})O(logn​) remained. The source of this mission, an OpenAI preprint dated September 24, 2026, claims a lower bound matching the upper bound up to a power of log⁡log⁡n\log\log nloglogn.

Setting

Let C=(cij)C=(c_{ij})C=(cij​) be symmetric nonnegative capacities on pairs of distinct vertices of [n][n][n]. The uniform sparsest-cut value is

OPT(C)=min⁡∅≠B⊊[n]∑i∈B, j∉Bcij∣B∣ (n−∣B∣).\mathrm{OPT}(C)=\min_{\varnothing\ne B\subsetneq[n]}\frac{\sum_{i\in B,\,j\notin B}c_{ij}}{|B|\,(n-|B|)}.OPT(C)=∅=B⊊[n]min​∣B∣(n−∣B∣)∑i∈B,j∈/B​cij​​.

A negative-type semimetric is d(i,j)=∥xi−xj∥2d(i,j)=\|x_i-x_j\|^2d(i,j)=∥xi​−xj​∥2 for vectors xix_ixi​ in a real Hilbert space, satisfying all triangle inequalities d(i,k)≤d(i,j)+d(j,k)d(i,k)\le d(i,j)+d(j,k)d(i,k)≤d(i,j)+d(j,k); distinct points may have distance zero. The Goemans–Linial value is

GL(C)=min⁡{∑i<jcij d(i,j) : ∑i<jd(i,j)=1, d negative type}.\mathrm{GL}(C)=\min\Bigl\{\sum_{i<j}c_{ij}\,d(i,j)\ :\ \sum_{i<j}d(i,j)=1,\ d\ \text{negative type}\Bigr\}.GL(C)=min{i<j∑​cij​d(i,j) : i<j∑​d(i,j)=1, d negative type}.

Normalized cut metrics are feasible, so GL(C)≤OPT(C)\mathrm{GL}(C)\le\mathrm{OPT}(C)GL(C)≤OPT(C). The integrality gap of CCC is OPT(C)/GL(C)\mathrm{OPT}(C)/\mathrm{GL}(C)OPT(C)/GL(C).

Formalization targets

Goal: Theorem 1.1

There are an absolute c>0c>0c>0 and instances C(j)C^{(j)}C(j) on nj→∞n_j\to\inftynj​→∞ vertices with GL(C(j))>0\mathrm{GL}(C^{(j)})>0GL(C(j))>0 such that, for all large jjj,

OPT(C(j))GL(C(j)) ≥ c log⁡nj(log⁡log⁡nj)3.\frac{\mathrm{OPT}(C^{(j)})}{\mathrm{GL}(C^{(j)})}\ \ge\ c\,\frac{\sqrt{\log n_j}}{(\log\log n_j)^3}.GL(C(j))OPT(C(j))​ ≥ c(loglognj​)3lognj​​​.

The Lean statement OAI.UniformSparsestCut.mainGap is open on the platform.

Significance

The bound reaches the exponent 1/21/21/2 in log⁡n\log nlogn, so the Arora–Rao–Vazirani analysis of the Goemans–Linial SDP is tight for uniform demands up to the factor (log⁡log⁡n)3(\log\log n)^3(loglogn)3. It strengthens the Devanur–Khot–Saket–Vishnoi and Kane–Meka refutations of the uniform constant-gap conjecture. Geometrically, it gives a finite negative-type metric for which every 111-Lipschitz map into ℓ1\ell_1ℓ1​ preserves only a (log⁡log⁡n)3/log⁡n(\log\log n)^3/\sqrt{\log n}(loglogn)3/logn​ fraction of the average distance. The theorem concerns the basic Goemans–Linial SDP only, not strengthened hierarchies.

The result is claimed in an OpenAI preprint; it has not been peer reviewed and no machine-checked proof exists.

Difficulty

The sharp general-demand lower bounds (Heisenberg group, Naor–Young) control worst-pair distortion, and Cheeger–Kleiner–Naor note that those examples have constant average distortion into a line, so they cannot give a diverging uniform gap. A uniform-demand example must force every ℓ1\ell_1ℓ1​ contraction to lose most of the average distance over all pairs, not merely on a few pairs. At the same time the constructed distance must remain of negative type and satisfy all triangle inequalities exactly, and the passage from a probability measure on points to an exactly uniform demand with positive SDP value must be carried out without losing the bound.

Formalization scope

  • Capacity n stores cap : Fin n → Fin n → ℝ, nonnegative, symmetric, zero on the diagonal.
  • cutRatio C B is the crossing capacity of BBB divided by ∣B∣(n−∣B∣)|B|(n-|B|)∣B∣(n−∣B∣); OPT C is the sInf over nonempty proper Finsets BBB.
  • NegativeType d requires points x : Fin n → EuclideanSpace ℝ (Fin n) with d=∥xi−xj∥2d=\|x_i-x_j\|^2d=∥xi​−xj​∥2 and all triangle inequalities; Feasible d adds ∑i<jd(i,j)=1\sum_{i<j}d(i,j)=1∑i<j​d(i,j)=1; glValue C is the sInf of ∑i<jcijd(i,j)\sum_{i<j}c_{ij}d(i,j)∑i<j​cij​d(i,j) over feasible ddd.
  • The goal asks for a sequence nj≥2n_j\ge2nj​≥2, nj→∞n_j\to\inftynj​→∞, with GL>0\mathrm{GL}>0GL>0 for every jjj (ruling out a gap that is trivial because of division by zero) and the bound eventually.

A complete development needs cut-cone duality, a Hilbert-space kernel construction with rounded charts, Gaussian direction families, and contraction estimates for maps into ℓ1\ell_1ℓ1​. Cut-cone duality and the reduction to exactly uniform demands are reusable. Contributions formalizing Lemma 6.1 (cut-cone duality) or the weak duality GL≤OPT\mathrm{GL}\le\mathrm{OPT}GL≤OPT are welcome first steps.

Selected references

  • OpenAI, Near-square-root logarithmic integrality gaps for uniform sparsest cut, preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Near-square-root-logarithmic-integrality-gaps-for-uniform-sparsest-cut-September-24-2026/Near-square-root-logarithmic-integrality-gaps-for-uniform-sparsest-cut-September-24-2026.pdf
  • S. Arora, S. Rao, U. Vazirani, Expander flows, geometric embeddings and graph partitioning, J. ACM, 2009. https://doi.org/10.1145/1502793.1502794
  • N. R. Devanur, S. A. Khot, R. Saket, N. K. Vishnoi, Integrality gaps for sparsest cut and minimum linear arrangement problems, STOC 2006. https://doi.org/10.1145/1132516.1132594
  • D. M. Kane, R. Meka, A PRG for Lipschitz functions of polynomials with applications to sparsest cut, STOC 2013. https://doi.org/10.1145/2488608.2488610
  • S. A. Khot, N. K. Vishnoi, The Unique Games Conjecture, integrality gap for cut problems and embeddability of negative type metrics into ℓ1\ell_1ℓ1​, J. ACM, 2015. https://doi.org/10.1145/2629614
  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica, 1995. https://doi.org/10.1007/BF01200757
  • Y. Rabinovich, On average distortion of embedding metrics into the line, Discrete Comput. Geom., 2008. https://doi.org/10.1007/s00454-007-9047-5
  • J. Cheeger, B. Kleiner, A. Naor, A (log⁡n)Ω(1)(\log n)^{\Omega(1)}(logn)Ω(1) integrality gap for the Sparsest Cut SDP, FOCS 2009. https://arxiv.org/abs/0910.2024
  • A. Naor, R. Young, Vertical perimeter versus horizontal perimeter, Ann. of Math., 2018. https://doi.org/10.4007/annals.2018.188.1.4
2 thms1 active userReviewed
CombinatoricsMarkov Chain·Captain: wurtle

An FPRAS for Cell-Bounded Contingency TablesResearch Paper

Motivation

A cell-bounded contingency table is a nonnegative integer matrix with prescribed row sums, column sums and an upper bound on each entry. Bounds of zero forbid entries (structural zeros), and bounds of one give bipartite graphs with prescribed degrees, so the model contains degree-constrained subgraph counting and perfect matchings (the permanent). Counting such tables is central to exact conditional inference in statistics and is equivalent to counting integer flows in bipartite networks. Exact counting is #P-complete even for two-row ordinary tables (Dyer, Kannan and Mount, 1997), so the target is a fully polynomial randomized approximation scheme (FPRAS): a randomized algorithm that, for any ε,δ∈(0,1)\varepsilon,\delta\in(0,1)ε,δ∈(0,1), returns a (1±ε)(1\pm\varepsilon)(1±ε)-approximation with probability at least 1−δ1-\delta1−δ in time polynomial in the input length, 1/ε1/\varepsilon1/ε and log⁡(1/δ)\log(1/\delta)log(1/δ).

Timeline

  • 1986. Jerrum, Valiant and Vazirani relate approximate counting and sampling through self-reducibility.
  • 1995. Diaconis and Gangolli survey fixed-margin arrays (doi:10.1007/978-1-4612-0801-3_3).
  • 1997. Dyer, Kannan and Mount prove #P-completeness for two-row tables and give geometric methods for large margins; Morris improves the margin requirements in 2002 (doi:10.1002/rsa.10049).
  • 2000–2003. Dyer and Greenhill (two rows), Cryan and Dyer (fixed number of rows), Dyer (dynamic programming) give approximate counting for a bounded number of rows.
  • 2004. Jerrum, Sinclair and Vigoda's permanent algorithm yields an FPRAS for all 0/1 caps (doi:10.1145/1008731.1008738).
  • 2009–2010. Barvinok bounds weighted table counts within NO(m+n)N^{O(m+n)}NO(m+n) (doi:10.1093/imrn/rnn133); Barvinok, Luria, Samorodnitsky and Yong give quasipolynomial randomized approximations for smooth margins (doi:10.1002/rsa.20301).
  • 2010. Cryan, Dyer and Randall give an FPRAS for integral flows with large tight capacities and for cell-bounded tables with a fixed number of rows, and list unrestricted approximation as open (doi:10.1137/060650544).
  • 2023. Brändén, Leake and Pak give capacity lower bounds via Lorentzian polynomials (doi:10.1007/s11856-022-2364-9); Guo and Jerrum count vertices of certain 0/1 polytopes (doi:10.1007/s00454-022-00406-8).

The source of this mission, an OpenAI preprint dated September 24, 2026, claims an FPRAS with no restriction on dimensions, margins or caps.

Setting

For positive integers m,nm,nm,n, margins r∈Z≥0mr\in\mathbb Z_{\ge0}^mr∈Z≥0m​, c∈Z≥0nc\in\mathbb Z_{\ge0}^nc∈Z≥0n​ with equal totals, and bounds b∈Z≥0m×nb\in\mathbb Z_{\ge0}^{m\times n}b∈Z≥0m×n​,

Ω(r,c,b)={X∈Z≥0m×n: ∑jXij=ri, ∑iXij=cj, Xij≤bij},Z(r,c,b)=∣Ω(r,c,b)∣.\Omega(r,c,b)=\Bigl\{X\in\mathbb Z_{\ge0}^{m\times n}:\ \textstyle\sum_jX_{ij}=r_i,\ \sum_iX_{ij}=c_j,\ X_{ij}\le b_{ij}\Bigr\},\qquad Z(r,c,b)=|\Omega(r,c,b)|.Ω(r,c,b)={X∈Z≥0m×n​: ∑j​Xij​=ri​, ∑i​Xij​=cj​, Xij​≤bij​},Z(r,c,b)=∣Ω(r,c,b)∣.

All numbers are written in binary, so a capacity can be exponentially larger than its description. LLL denotes the total binary input length including the rational accuracy parameters ε,δ\varepsilon,\deltaε,δ. The algorithm uses independent unbiased random bits and is charged for every bit operation.

Formalization targets

Goal: Theorem 1.1 (FPRAS)

One randomized algorithm, given (r,c,b)(r,c,b)(r,c,b) and rational ε,δ∈(0,1)\varepsilon,\delta\in(0,1)ε,δ∈(0,1), outputs a nonnegative rational Z^\widehat ZZ with

Pr⁡[(1−ε)Z(r,c,b)≤Z^≤(1+ε)Z(r,c,b)]≥1−δ,\Pr\bigl[(1-\varepsilon)Z(r,c,b)\le\widehat Z\le(1+\varepsilon)Z(r,c,b)\bigr]\ge1-\delta,Pr[(1−ε)Z(r,c,b)≤Z≤(1+ε)Z(r,c,b)]≥1−δ,

outputs exactly 000 on every execution when Ω(r,c,b)=∅\Omega(r,c,b)=\emptysetΩ(r,c,b)=∅, and performs at most poly(L,ε−1,log⁡δ−1)\mathrm{poly}(L,\varepsilon^{-1},\log\delta^{-1})poly(L,ε−1,logδ−1) bit operations on every execution. Lean: OAI.ContingencyTables.counting, open on the platform.

Significance

The theorem would answer the question left open by Cryan, Dyer and Randall in 2010: approximate counting of cell-bounded tables (equivalently, integral flows in bipartite networks with arbitrary capacities) without any fixed-dimension, large-capacity, density or balance assumption, with polynomial cost on every execution and no oracle. It contains the 0/1 case handled by Jerrum–Sinclair–Vigoda and the fixed-row results as special cases, and via self-reduction it gives approximate samplers for integer flows. A companion preprint of the same family treats exact uniform sampling of uncapped tables. The result is claimed in an OpenAI preprint that has not been peer reviewed; no machine-checked proof exists.

Difficulty

Two regimes clash. Large capacities cannot be expanded into that many binary choices without exponential blow-up, so geometric (polytope volume) methods are needed for them; but those methods require every cell to be large. Small capacities need combinatorial Markov chains, which in turn must preserve exactly one unit of mass per table despite the multiplicities created by encoding a bounded integer by binary choices. Arbitrary mixtures of tiny and huge caps, zero caps and variable dimensions defeat each approach used alone.

Formalization scope

  • Algorithms are MatchingFPRAS.RandomMachines: a single finite transition table on an 8-symbol alphabet, one random bit per tick; time is counted in ticks (polynomially equivalent to bit operations).
  • Inputs are delimited binary encodings of m,n,r,c,bm,n,r,c,bm,n,r,c,b and of ε,δ\varepsilon,\deltaε,δ as reduced fractions (encodeCountingInput); outputs are encoded rationals.
  • count r c b is Fintype.card of bounded tables.
  • CountingStatement: one machine and constants C>0C>0C>0, ddd; with t=C(∣input∣+⌈1/ε⌉+⌈log⁡2⌈1/δ⌉⌉+1)dt=C(|\text{input}|+\lceil1/\varepsilon\rceil+\lceil\log_2\lceil1/\delta\rceil\rceil+1)^dt=C(∣input∣+⌈1/ε⌉+⌈log2​⌈1/δ⌉⌉+1)d, every tape of length ttt halts with a rational q≥0q\ge0q≥0; if count = 0 every tape outputs 000; and the fraction of tapes with (1−ε)Z≤q≤(1+ε)Z(1-\varepsilon)Z\le q\le(1+\varepsilon)Z(1−ε)Z≤q≤(1+ε)Z is at least 1−δ1-\delta1−δ.
  • m,n≥1m,n\ge1m,n≥1 and equal margin totals are hypotheses, as in the paper.

Needed infrastructure: randomized Turing machines, Markov chain mixing and Poincaré inequalities, log-concavity (Prékopa–Leindler), and network-flow feasibility. The machine model and encodings are shared with the exact-sampling mission of this family. Contributions formalizing Proposition 2.2 (contamination bound), Lemma 3.3 (log-concave tree marginals), Proposition 6.1 (weight evaluation) or Proposition 7.2 are welcome.

Selected references

  • OpenAI, An FPRAS for Cell-Bounded Contingency Tables, preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/An-FPRAS-for-Cell-Bounded-Contingency-Tables-September-24-2026/main.pdf
  • M. Cryan, M. Dyer, D. Randall, Approximately counting integral flows and cell-bounded contingency tables, SIAM J. Comput., 2010. https://doi.org/10.1137/060650544
  • M. Jerrum, A. Sinclair, E. Vigoda, A polynomial-time approximation algorithm for the permanent of a matrix with nonnegative entries, J. ACM, 2004. https://doi.org/10.1145/1008731.1008738
  • B. Morris, Improved bounds for sampling contingency tables, Random Structures Algorithms, 2002. https://doi.org/10.1002/rsa.10049
  • A. Barvinok, Asymptotic estimates for the number of contingency tables, integer flows, and volumes of transportation polytopes, IMRN, 2009. https://doi.org/10.1093/imrn/rnn133
  • A. Barvinok, Z. Luria, A. Samorodnitsky, A. Yong, An approximation algorithm for counting contingency tables, Random Structures Algorithms, 2010. https://doi.org/10.1002/rsa.20301
  • P. Brändén, J. Leake, I. Pak, Lower bounds for contingency tables via Lorentzian polynomials, Israel J. Math., 2023. https://doi.org/10.1007/s11856-022-2364-9
  • H. Guo, M. Jerrum, Counting vertices of integral polytopes defined by facets, Discrete Comput. Geom., 2023. https://doi.org/10.1007/s00454-022-00406-8
2 thms1 active userReviewed
CombinatoricsMarkov Chain·Captain: wurtle

Exact Uniform Sampling of Contingency Tables with Arbitrary MarginsResearch Paper

Motivation

A contingency table with margins r=(r1,…,rm)r=(r_1,\dots,r_m)r=(r1​,…,rm​) and c=(c1,…,cn)c=(c_1,\dots,c_n)c=(c1​,…,cn​) is a nonnegative integer matrix whose row sums are rrr and column sums are ccc. Such tables are the basic objects of categorical data analysis: conditional tests of independence (Fisher's exact test and its generalizations) compare an observed table with the uniform distribution on all tables with the same margins, and in practice this requires sampling uniformly from that set. The set is finite but can be exponentially large in the input size when margins are written in binary. Whether one can sample it (approximately or exactly) in time polynomial in both dimensions and in the bit length of the margins has been open since the 1990s.

Timeline

  • 1995. Diaconis and Gangolli survey enumeration and random generation of rectangular arrays with fixed margins (doi:10.1007/978-1-4612-0801-3_3).
  • 1997–2002. Dyer, Kannan and Mount sample dense tables via the transportation polytope; Morris improves the dense bounds (doi:10.1002/rsa.10049).
  • 2000. Dyer and Greenhill give polynomial-time approximate counting and sampling for two-rowed tables.
  • 2003. Kijima and Matsui give exact sampling for two rows by coupling from the past; Cryan and Dyer approximately count for a fixed number of rows (doi:10.1016/S0022-0000(03)00014-X); Dyer gives an exact dynamic-programming sampler for fixed row count.
  • 2006. Cryan, Dyer, Goldberg, Jerrum and Martin prove rapid mixing of the heat-bath chain for every fixed number of rows (doi:10.1137/S0097539703434243).
  • 2016. DeSalvo and Zhao give an exact divide-and-conquer sampler whose polynomial runtime is conditional on a conjecture (arXiv:1507.00070).
  • 2021. Arman, Gao and Wormald give linear-time exact generation for sparse margins with 5Δ4<N5\Delta^4<N5Δ4<N, and report that no polynomial-time approximately uniform sampler for arbitrary margins was known (arXiv:2104.09413).
  • 2024. Göbel, Liu, Manurangsi and Pappik turn sufficiently accurate samplers into perfect samplers (arXiv:2410.00882).

The source of this mission, an OpenAI preprint dated September 24, 2026, claims both an almost-uniform and an exactly uniform sampler for arbitrary margins.

Setting

For nonnegative integer vectors r∈Z≥0mr\in\mathbb Z_{\ge0}^mr∈Z≥0m​, c∈Z≥0nc\in\mathbb Z_{\ge0}^nc∈Z≥0n​ with common total N=∑iri=∑jcjN=\sum_ir_i=\sum_jc_jN=∑i​ri​=∑j​cj​, let

Ω(r,c)={X∈Z≥0m×n: ∑jXij=ri, ∑iXij=cj}.\Omega(r,c)=\Bigl\{X\in\mathbb Z_{\ge0}^{m\times n}:\ \textstyle\sum_jX_{ij}=r_i,\ \sum_iX_{ij}=c_j\Bigr\}.Ω(r,c)={X∈Z≥0m×n​: ∑j​Xij​=ri​, ∑i​Xij​=cj​}.

The input is (m,n,r,c)(m,n,r,c)(m,n,r,c) with margins in binary, so its length is polynomial in mmm, nnn and log⁡(N+1)\log(N+1)log(N+1). An algorithm uses unbiased random bits. For laws on a finite set, TV(μ,ν)=12∑x∣μ(x)−ν(x)∣\mathrm{TV}(\mu,\nu)=\tfrac12\sum_x|\mu(x)-\nu(x)|TV(μ,ν)=21​∑x​∣μ(x)−ν(x)∣.

Formalization targets

Milestone: Theorem 1.1(i) (almost-uniform, bounded time)

For each k≥1k\ge1k≥1, an algorithm halts on every execution in time poly(m,n,log⁡(N+1),k)\mathrm{poly}(m,n,\log(N+1),k)poly(m,n,log(N+1),k) and outputs X∈Ω(r,c)X\in\Omega(r,c)X∈Ω(r,c) with

TV(law(X), Unif Ω(r,c))≤2−k.\mathrm{TV}\bigl(\mathrm{law}(X),\ \mathrm{Unif}\,\Omega(r,c)\bigr)\le2^{-k}.TV(law(X), UnifΩ(r,c))≤2−k.

Goal: Theorem 1.1(ii) (exact uniform sampling)

An algorithm outputs X∈Ω(r,c)X\in\Omega(r,c)X∈Ω(r,c) with

Pr⁡[X=Y]=1∣Ω(r,c)∣(Y∈Ω(r,c)),\Pr[X=Y]=\frac1{|\Omega(r,c)|}\quad(Y\in\Omega(r,c)),Pr[X=Y]=∣Ω(r,c)∣1​(Y∈Ω(r,c)),

terminates almost surely, and has expected running time poly(m,n,log⁡(N+1))\mathrm{poly}(m,n,\log(N+1))poly(m,n,log(N+1)), with polynomials uniform over all margins. Lean: OAI.ContingencyTables.exactSampling, open on the platform.

Significance

The theorem would give the first unconditional polynomial-time exact uniform sampler for contingency tables with arbitrary dimensions and binary margins, with no positivity, sparsity, balance or fixed-dimension restriction, and with all random bits and integer arithmetic counted. It removes the fixed-row restriction of the Dyer–Greenhill and Cryan–Dyer line and the sparsity condition of Arman–Gao–Wormald. The exponents are deliberately large; the contribution is the existence of a polynomial algorithm. Exact time is polynomial in expectation, not on every run. A companion preprint of the same family gives an FPRAS for counting tables with cell bounds. The result is claimed in an OpenAI preprint that has not been peer reviewed; no machine-checked proof exists.

Difficulty

Markov-chain approaches mix rapidly for a fixed number of rows, but the bounds degrade with the row count. Geometric (polytope) methods need all margins to be large, and fail when some margins are tiny while others are huge, which is exactly the mixed regime that arbitrary binary margins allow. Making the sampler exact rather than approximate adds a separate difficulty: the residual-mixture correction requires that the actual output law of the approximate sampler can be computed, with cost polynomial in the accuracy bits.

Formalization scope

  • Algorithms are MatchingFPRAS.RandomMachines: a single finite transition table over an 8-symbol tape alphabet; each tick reads one random bit and either moves, writes or halts. Time is counted in ticks, polynomially equivalent to bit operations.
  • Inputs are binary encodings with delimiters (encodeNat, encodeMargins); outputs are row-major matrix encodings left on the tape (OutputsTable).
  • tableMass A input t X and haltMass A input t are the fractions of the 2t2^t2t random tapes of length ttt that halt with output XXX, respectively halt at all, within ttt ticks.
  • ExactSamplingStatement: one machine and constants C>0C>0C>0, ddd such that for all m,n,r,cm,n,r,cm,n,r,c with equal totals: every halted output is in Ω(r,c)\Omega(r,c)Ω(r,c); tableMass → 1/|Ω| for each table; haltMass → 1; and ∑t(1−haltMasst)\sum_t(1-\text{haltMass}_t)∑t​(1−haltMasst​), the expected running time, is at most C(m+n+⌈log⁡2(N+1)⌉+1)dC(m+n+\lceil\log_2(N+1)\rceil+1)^dC(m+n+⌈log2​(N+1)⌉+1)d.
  • BoundedSamplingStatement: every tape of length t=C(m+n+⌈log⁡2(N+1)⌉+k+1)dt=C(m+n+\lceil\log_2(N+1)\rceil+k+1)^dt=C(m+n+⌈log2​(N+1)⌉+k+1)d halts with a feasible table, and the TV distance to uniform is at most 2−k2^{-k}2−k.
  • Dimensions m,nm,nm,n vary freely (zero allowed); the machine is uniform, not one per size.

Needed infrastructure: randomized Turing machines and their output laws, Markov chain mixing (conductance, Poincaré inequalities), and log-concave sampling. Contributions formalizing Theorem 4.3 (graph Poincaré inequality), Theorem 5.1 (dense completion draw), Proposition 6.3 (law tabulation) or the Göbel–Liu–Manurangsi–Pappik correction are welcome.

Selected references

  • OpenAI, Exact Uniform Sampling of Contingency Tables with Arbitrary Margins, preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Exact-Uniform-Sampling-of-Contingency-Tables-with-Arbitrary-Margins-September-24-2026/main.pdf
  • P. Diaconis, A. Gangolli, Rectangular arrays with fixed margins, Discrete Probability and Algorithms, 1995. https://doi.org/10.1007/978-1-4612-0801-3_3
  • B. Morris, Improved bounds for sampling contingency tables, Random Structures Algorithms, 2002. https://doi.org/10.1002/rsa.10049
  • M. Cryan, M. Dyer, A polynomial-time algorithm to approximately count contingency tables when the number of rows is constant, J. Comput. System Sci., 2003. https://doi.org/10.1016/S0022-0000(03)00014-X
  • M. Cryan, M. Dyer, L. A. Goldberg, M. Jerrum, R. Martin, Rapidly mixing Markov chains for sampling contingency tables with a constant number of rows, SIAM J. Comput., 2006. https://doi.org/10.1137/S0097539703434243
  • S. DeSalvo, J. Y. Zhao, Random sampling of contingency tables via probabilistic divide-and-conquer, preprint, 2016. https://arxiv.org/abs/1507.00070
  • A. Arman, P. Gao, N. Wormald, Linear-time uniform generation of random sparse contingency tables with specified marginals, preprint, 2021. https://arxiv.org/abs/2104.09413
  • A. Göbel, J. Liu, P. Manurangsi, M. Pappik, Perfect sampling from rapidly mixing Markov chains, preprint, 2024. https://arxiv.org/abs/2410.00882
3 thms1 active userReviewed
CombinatoricsMarkov Chain·Captain: wurtle

Approximate counting of common bases of two matroidsResearch Paper

Motivation: counting common bases of two matroids

A matroid abstracts linear independence: a finite ground set with a family of "independent" subsets that is closed under subsets and satisfies the augmentation axiom. Its maximal independent sets, the bases, all have the same size, the rank. Matroid intersection, finding a common basis of two matroids, is one of the central polynomially solvable problems of combinatorial optimization (Edmonds, 1979). Counting common bases is much harder. It contains counting bases of a single matroid (take both matroids equal) and counting perfect matchings of a bipartite graph (take two partition matroids), and exact counting is #P-hard. The natural goal is therefore a fully polynomial randomized approximation scheme (FPRAS): an algorithm whose output is within relative error ϵ\epsilonϵ with probability at least 1−δ1-\delta1−δ, in time polynomial in the input size, ϵ−1\epsilon^{-1}ϵ−1 and log⁡δ−1\log\delta^{-1}logδ−1. The matroids should be accessible only through independence oracles.

Background and timeline

  • 1979 — Edmonds gives a polynomial-time algorithm for matroid intersection (Ann. Discrete Math. 4).
  • 1992 — Feder and Mihail give sampling and approximate counting for balanced matroids (STOC 1992).
  • 2004 — Jerrum, Sinclair and Vigoda give an FPRAS for the permanent of nonnegative matrices, covering bipartite perfect matchings (J. ACM 51).
  • 2007 — Barvinok and Samorodnitsky obtain coarse logarithmic estimates via random weighting (Israel J. Math. 158).
  • 2016–2017 — Anari, Oveis Gharan and Rezaei prove rapid mixing for strongly Rayleigh distributions (COLT 2016); Anari and Oveis Gharan obtain exponential-factor approximations from real-stable polynomials (arXiv:1702.02937).
  • 2020 — Brändén and Huh develop Lorentzian polynomials (Ann. of Math. 192).
  • 2021 — Anari, Oveis Gharan and Vinzant give a deterministic 2O(r)2^{O(r)}2O(r)-factor approximation for common bases of rank-rrr matroids (Duke Math. J. 170); Cryan, Guo and Mousa ask for a rapidly mixing chain for common bases (Ann. Probab. 49).
  • 2023–2024 — Liu's dissertation records the FPRAS question for two arbitrary oracle matroids (Section 13.2, Problem 1, thesis); Anari, Liu, Oveis Gharan and Vinzant give an FPRAS for bases of a single matroid (Ann. of Math. 199).
  • 2026 — An OpenAI preprint, Approximate counting of common bases of two matroids (OpenAI Math Release, September 23, 2026), claims an independence-oracle FPRAS for common bases of two arbitrary matroids of equal rank. It has not been peer reviewed, and its proof is not formally verified.

Setting

Let M1,M2M_1,M_2M1​,M2​ be matroids of the same rank rrr on the ground set [n][n][n], supplied by exact independence oracles: a query is an nnn-bit indicator of a subset, and the answer says whether it is independent. The quantity to estimate is

Z(M1,M2)=∣B(M1)∩B(M2)∣,Z(M_1,M_2)=\bigl|\mathcal B(M_1)\cap\mathcal B(M_2)\bigr|,Z(M1​,M2​)=​B(M1​)∩B(M2​)​,

the number of common bases. The input also contains rationals ϵ,δ∈(0,1)\epsilon,\delta\in(0,1)ϵ,δ∈(0,1); ℓ\ellℓ denotes the binary encoding length of n,r,ϵ,δn,r,\epsilon,\deltan,r,ϵ,δ. Oracle calls are counted separately from the bit operations used to build queries and process answers.

Formalization targets

Goal: an FPRAS for common bases in the independence-oracle model (Theorem 1.1)

There is a single randomized oracle algorithm and a fixed polynomial ppp such that, for all rank-rrr matroids M1,M2M_1,M_2M1​,M2​ on [n][n][n] and rational ϵ,δ∈(0,1)\epsilon,\delta\in(0,1)ϵ,δ∈(0,1), it outputs a nonnegative rational Z^\widehat ZZ with

Pr⁡[(1−ϵ)Z≤Z^≤(1+ϵ)Z] ≥ 1−δ,\Pr\bigl[(1-\epsilon)Z\le\widehat Z\le(1+\epsilon)Z\bigr]\ \ge\ 1-\delta,Pr[(1−ϵ)Z≤Z≤(1+ϵ)Z] ≥ 1−δ,

the output is always 000 when Z=0Z=0Z=0, and on every execution the numbers of oracle calls and bit operations are at most p(n,ℓ,ϵ−1,log⁡δ−1)p(n,\ell,\epsilon^{-1},\log\delta^{-1})p(n,ℓ,ϵ−1,logδ−1). The algorithm uses only independent fair random bits. The goal statement is published on the platform with status Open.

Significance

The result itself. Theorem 1.1 answers the question recorded by Liu (2023) and extends both the single-matroid FPRAS of Anari–Liu–Oveis Gharan–Vinzant and, in the bipartite-matching case, the permanent FPRAS of Jerrum–Sinclair–Vigoda. It needs no representation of either matroid and holds for unrestricted rank, whereas the previous oracle-model scheme was fully polynomial only for r=O(log⁡n)r=O(\log n)r=O(logn). The source derives FPRASs and almost-uniform samplers for common independent sets of prescribed, unrestricted or maximum size, even for matroids of different ranks (Corollary 9.1).

Formalizing it. The Lean goal fixes a concrete machine model, so it states a complete complexity-theoretic claim with no informal "polynomial-time" left. A proof would combine Lorentzian-polynomial inequalities (Brändén–Huh), a transport inequality, a Markov-chain Poincaré bound, and simulated annealing with explicit bit accounting. Each layer is reusable: Mathlib has matroids, but no Lorentzian polynomials, no mixing-time theory of this kind and no randomized-algorithm framework.

Difficulty

For one matroid, the bases-exchange walk mixes rapidly because the basis generating polynomial is log-concave. The intersection of two matroids has no such structure: its common bases need not be connected by short exchanges, and the natural chain on transversals of the paired construction passes through states of exponentially small weight, so positivity and connectivity alone give no polynomial bound. Earlier approaches either lost exponential factors or needed real-stable polynomials that an independence oracle does not provide. The proof must balance defect classes, as in Jerrum–Sinclair–Vigoda, and establish a Poincaré inequality for observables from matroid-theoretic inequalities alone.

Formalization scope

  • Matroids are Mathlib Matroid (Fin n) with ground set univ and eRank = r for both. Oracles are functions Bool → (Fin n → Bool) → Bool, required to answer independence queries exactly.
  • The algorithm is one fixed Program k s for a register machine with k+8k+8k+8 binary stacks: push/pop, fair coin, oracle query on an nnn-bit register, and halt with output (register 6)/(register 7) as a rational. Inputs n,r,ϵ,δn,r,\epsilon,\deltan,r,ϵ,δ are loaded in binary.
  • The resource bound is B=C(1+n+ℓ+⌈ϵ−1⌉+⌈log⁡2⌈δ−1⌉⌉)dB=C(1+n+\ell+\lceil\epsilon^{-1}\rceil+\lceil\log_2\lceil\delta^{-1}\rceil\rceil)^dB=C(1+n+ℓ+⌈ϵ−1⌉+⌈log2​⌈δ−1⌉⌉)d with constants C>0,dC>0,dC>0,d fixed before the instance. For every infinite bit stream, the run halts within BBB steps with at most BBB oracle calls and BBB bit operations, outputs a nonnegative rational, and outputs 000 when Z=0Z=0Z=0.
  • The success probability is expressed by counting: at least (1−δ)2B(1-\delta)2^B(1−δ)2B of the 2B2^B2B bit strings of length BBB lead to an output within relative error ϵ\epsilonϵ of ZZZ.
  • A trivializing reading is excluded: the program and constants come first and must work for every n,rn,rn,r, matroid pair, oracle and ϵ,δ\epsilon,\deltaϵ,δ.
  • Needed infrastructure: Lorentzian polynomials and the homogeneous Tutte inequalities, Markov chains on finite state spaces with Poincaré/variance bounds, annealing estimators, and correctness proofs for the register machine. Contributions on any layer are welcome.

Selected references

  • J. Edmonds, Matroid intersection, Ann. Discrete Math. 4 (1979), 39–49. https://doi.org/10.1016/S0167-5060(08)70817-3
  • T. Feder and M. Mihail, Balanced matroids, STOC 1992. https://doi.org/10.1145/129712.129716
  • M. Jerrum, A. Sinclair and E. Vigoda, A polynomial-time approximation algorithm for the permanent of a matrix with nonnegative entries, J. ACM 51 (2004), 671–697. https://doi.org/10.1145/1008731.1008738
  • A. Barvinok and A. Samorodnitsky, Random weighting, asymptotic counting, and inverse isoperimetry, Israel J. Math. 158 (2007). https://doi.org/10.1007/s11856-007-0008-8
  • P. Brändén and J. Huh, Lorentzian polynomials, Ann. of Math. 192 (2020), 821–891. https://doi.org/10.4007/annals.2020.192.3.4
  • N. Anari, S. Oveis Gharan and C. Vinzant, Log-concave polynomials I, Duke Math. J. 170 (2021). https://doi.org/10.1215/00127094-2020-0091
  • M. Cryan, H. Guo and G. Mousa, Modified log-Sobolev inequalities for strongly log-concave distributions, Ann. Probab. 49 (2021). https://doi.org/10.1214/20-AOP1453
  • N. Anari, K. Liu, S. Oveis Gharan and C. Vinzant, Log-concave polynomials II: high-dimensional walks and an FPRAS for counting bases of a matroid, Ann. of Math. 199 (2024). https://doi.org/10.4007/annals.2024.199.1.4
  • K. Liu, Spectral Independence: A New Tool to Analyze Markov Chains, PhD thesis, University of Washington, 2023. https://kuikuiliu.github.io/files/dissertation.pdf
  • OpenAI, Approximate counting of common bases of two matroids, OpenAI Math Release preprint, September 23, 2026 (Theorem 1.1, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Approximate-counting-of-common-bases-of-two-matroids-September-23-2026/main.pdf
2 thms1 active userReviewed
CombinatoricsInformation Theory·Captain: wurtle

Entropy and Face Dimension of the Perfect-Matching PolytopeResearch Paper

Motivation: how much randomness fits under prescribed edge probabilities?

Prescribing the probability xex_exe​ that each edge belongs to a random perfect matching does not determine the distribution of the whole matching. The largest entropy compatible with those edge probabilities, H(x)H(x)H(x), measures how many matching choices can coexist with them, and it links three topics: lower bounds on the number of perfect matchings, convex-optimization approaches to approximate counting, and the facial geometry of the perfect-matching polytope. For bipartite graphs, Schrijver's permanent inequality (Schrijver 1998) and Gurvits' weighted Bethe inequality (Gurvits 2011) give the clean comparison H≥F−BH\ge F-BH≥F−B between HHH and sums of one-coordinate entropies. For general graphs, odd cuts make such comparisons harder, and Anari, Oveis Gharan and Vinzant (2018) asked (FOCS version, Conjecture 6) whether the maximum binary-coordinate entropy over the polytope exceeds log⁡N(G)\log N(G)logN(G) by only a linear function of the number of vertices.

Timeline

  • 1965 — Edmonds describes the perfect-matching polytope by degree equations, nonnegativity and odd-cut inequalities (Edmonds 1965).
  • 1982 — Naddef (Math. Program. 1982) and Edmonds, Pulleyblank and Lovász (Combinatorica 1982) determine the dimension of the perfect-matching polytope.
  • 1986 — Chung, Graham, Frankl and Shearer: Shearer's entropy inequality (JCTA 1986), which gives H≤FH\le FH≤F.
  • 1998–2011 — Schrijver's lower bound for regular bipartite graphs and Gurvits' Bethe-approximation form give H≥F−BH\ge F-BH≥F−B in the bipartite case.
  • 2011 — Esperet, Kardoš, King, Král' and Norine: exponentially many perfect matchings in bridgeless cubic graphs (Adv. Math. 2011).
  • 2018 — Anari, Oveis Gharan and Vinzant propose the entropy framework and the linear-error question (FOCS 2018).
  • 2022 — Ebrahimnejad, Nagda and Oveis Gharan count perfect matchings in regular expanding nonbipartite graphs (ITCS 2022).
  • 2026 — Abdi, Cornuéjols, Dadush and Dalirrooyfard bound the rank of tight odd cuts by 3m−13m-13m−1 and give exponential lower counts for regular graphs of degree at least four with an odd-cut condition (Proc. LMS 2026). An OpenAI preprint, Entropy and Face Dimension of the Perfect-Matching Polytope (OpenAI Math Release, September 23, 2026), claims the pointwise bound F−(2−2/m)B≤H≤FF-(2-2/m)B\le H\le FF−(2−2/m)B≤H≤F and a positive answer to the Anari–Oveis Gharan–Vinzant question. The preprint has not been peer reviewed, and its main theorem is not formally verified.

Setting

Let GGG be a finite loopless multigraph with vertex set VVV, ∣V∣=2m≥2|V|=2m\ge2∣V∣=2m≥2, and edge set EEE (parallel edges are distinct coordinates). A perfect matching is a set of edges covering each vertex exactly once; M(G)\mathcal M(G)M(G) is the set of them and N(G)=∣M(G)∣N(G)=|\mathcal M(G)|N(G)=∣M(G)∣. The perfect-matching polytope is

P(G)=conv⁡{1M:M∈M(G)}⊆RE.P(G)=\operatorname{conv}\{\mathbf 1_M: M\in\mathcal M(G)\}\subseteq\mathbb R^E .P(G)=conv{1M​:M∈M(G)}⊆RE.

For x∈P(G)x\in P(G)x∈P(G) put, with natural logarithms and 0log⁡0=00\log 0=00log0=0,

H(x)=max⁡{−∑MpMlog⁡pM: p a probability law on M(G), ∑MpM1M=x},H(x)=\max\Big\{-\sum_M p_M\log p_M:\ p\text{ a probability law on }\mathcal M(G),\ \sum_M p_M\mathbf 1_M=x\Big\},H(x)=max{−M∑​pM​logpM​: p a probability law on M(G), M∑​pM​1M​=x}, F(x)=−∑exelog⁡xe,B(x)=−∑e(1−xe)log⁡(1−xe).F(x)=-\sum_{e}x_e\log x_e,\qquad B(x)=-\sum_e(1-x_e)\log(1-x_e).F(x)=−e∑​xe​logxe​,B(x)=−e∑​(1−xe​)log(1−xe​).

Formalization targets

The milestones are listed in attack order: the face of a triangle expansion (Lemma 5.1), the coefficient-eight nonlinear bound (Theorem B.1), deterministic approximate counting (Theorem 1.2) and its singleton-loop variant (Proposition C.1).

Goal: Theorem 1.1 (pointwise entropy bounds)

If GGG has at least one perfect matching, then for every x∈P(G)x\in P(G)x∈P(G)

F(x)−(2−2m)B(x)  ≤  H(x)  ≤  F(x).F(x)-\Big(2-\frac 2m\Big)B(x)\;\le\;H(x)\;\le\;F(x).F(x)−(2−m2​)B(x)≤H(x)≤F(x).

The bound holds on the whole polytope, boundary included. The goal statement is published on the platform with status Open: no machine-checked proof exists yet.

Significance

The result itself. The upper bound is Shearer's inequality; the content is the lower comparison with coefficient 2−2/m2-2/m2−2/m, valid in nonbipartite graphs where the bipartite coefficient 111 fails (the source gives an eight-vertex counterexample, and Corollary 5.3 shows any universal coefficient is at least ln⁡3/(2ln⁡(3/2))\ln3/(2\ln(3/2))ln3/(2ln(3/2))). Optimizing over xxx gives log⁡N(G)≤max⁡P(G)(F+B)≤log⁡N(G)+3m−2\log N(G)\le\max_{P(G)}(F+B)\le\log N(G)+3m-2logN(G)≤maxP(G)​(F+B)≤logN(G)+3m−2, answering the Anari–Oveis Gharan–Vinzant question. Plugging in uniform marginals gives N(G)≥e−2(m−1)kmN(G)\ge e^{-2(m-1)}k^mN(G)≥e−2(m−1)km for graphs with kkk edge-disjoint perfect matchings, and for kkk-regular graphs whose odd cuts all have at least kkk edges. The same entropy framework yields a deterministic polynomial-time algorithm approximating N(G)N(G)N(G) within a factor 512n512^n512n for graphs with binary edge multiplicities (Theorem 1.2).

Formalizing it. The goal is a finite-dimensional inequality about polytopes and entropies, accessible to Mathlib's convexity and analysis libraries. The milestones extend the mission to a certified polynomial-time counting algorithm in Mathlib's Turing-machine model, a rare example of a formally specified approximation algorithm for a #P-hard quantity.

Difficulty

Shearer's inequality gives the upper bound in one line, and in bipartite graphs the lower bound follows from permanent inequalities that have no analogue for general graphs. A direct convexity argument fails because HHH is not differentiable across faces of P(G)P(G)P(G): at a boundary mean the maximum-entropy law lives on a lower-dimensional face, and its covariance is singular in the ambient edge space. The lower bound must therefore be proved face by face, and it needs the sharp geometric estimate ∣supp⁡x∣−dim⁡Fx≤3m−2|\operatorname{supp}x|-\dim F_x\le 3m-2∣suppx∣−dimFx​≤3m−2 for the minimal face FxF_xFx​ containing xxx; the earlier rank bound 3m−13m-13m−1 would not give the coefficient 2−2/m2-2/m2−2/m.

Formalization scope

  • Graphs are LooplessGraph V E: endpoint maps left right : E → V with left e ≠ right e, on finite types with decidable equality, so parallel edges are distinct elements of E. A perfect matching is a Finset E with exactly one edge at each vertex.
  • polytope G is the convex hull of matching indicator vectors in E → ℝ; maxMatchingEntropy G x is sSup of the entropies of probability vectors on matchings with mean x. On the polytope this set is nonempty, compact and bounded, so the sSup is the attained maximum; off the polytope it would be junk, which is why the hypothesis x ∈ G.polytope is required.
  • entropy uses Real.negMulLog, so 0log⁡0=00\log0=00log0=0; complementEntropy x is B(x)B(x)B(x).
  • Hypotheses: 0 < m, Fintype.card V = 2 * m, and a perfect matching exists. The coefficient is the real number 2 - 2 / m.
  • The counting milestones use Turing.TM2ComputableInPolyTime with finite alphabets, binary encodings of the vertex count and of positive multiplicity records, and require polynomially bounded output length.
  • Infrastructure needed: Edmonds' polytope description, faces and minimal faces of polytopes, exponential families and covariance, and polynomial-time computability of rational arithmetic in TM2. The polytope and entropy definitions are reusable for other counting missions.

Selected references

  • J. Edmonds, Maximum matching and a polyhedron with 0,1-vertices, J. Res. Nat. Bur. Standards B 69B (1965), 125–130. https://doi.org/10.6028/jres.069B.013
  • L. G. Valiant, The complexity of computing the permanent, Theoret. Comput. Sci. 8 (1979), 189–201. https://doi.org/10.1016/0304-3975(79)90044-6
  • D. Naddef, Rank of maximum matchings in a graph, Math. Program. 22 (1982), 52–70. https://doi.org/10.1007/BF01581025
  • J. Edmonds, W. R. Pulleyblank, L. Lovász, Brick decompositions and the matching rank of graphs, Combinatorica 2 (1982), 247–274. https://doi.org/10.1007/BF02579233
  • F. R. K. Chung, R. L. Graham, P. Frankl, J. B. Shearer, Some intersection theorems for ordered sets and graphs, J. Combin. Theory Ser. A 43 (1986), 23–37. https://doi.org/10.1016/0097-3165(86)90019-1
  • A. Schrijver, Counting 1-factors in regular bipartite graphs, J. Combin. Theory Ser. B 72 (1998), 122–135. https://doi.org/10.1006/jctb.1997.1798
  • N. Linial, A. Samorodnitsky, A. Wigderson, A deterministic strongly polynomial algorithm for matrix scaling and approximate permanents, Combinatorica 20 (2000), 545–568. https://doi.org/10.1007/s004930070007
  • M. Jerrum, A. Sinclair, E. Vigoda, A polynomial-time approximation algorithm for the permanent of a matrix with nonnegative entries, J. ACM 51 (2004), 671–697. https://doi.org/10.1145/1008731.1008738
  • L. Esperet, F. Kardoš, A. D. King, D. Král', S. Norine, Exponentially many perfect matchings in cubic graphs, Adv. Math. 227 (2011), 1646–1664. https://doi.org/10.1016/j.aim.2011.03.015
  • L. Gurvits, Unleashing the power of Schrijver's permanental inequality with the help of the Bethe approximation (2011). https://arxiv.org/abs/1106.2844v11
  • M. Cygan, M. Pilipczuk, R. Škrekovski, A bound on the number of perfect matchings in Klee-graphs, Discrete Math. Theor. Comput. Sci. 15 (2013), 37–52. https://doi.org/10.46298/dmtcs.633
  • A. Barvinok, Approximating permanents and hafnians, Discrete Analysis (2017). https://doi.org/10.19086/da.1244
  • N. Anari, S. Oveis Gharan, C. Vinzant, Log-concave polynomials, entropy, and a deterministic approximation algorithm for counting bases of matroids, FOCS 2018, 35–46. https://doi.org/10.1109/FOCS.2018.00013
  • F. Ebrahimnejad, A. Nagda, S. Oveis Gharan, Counting and sampling perfect matchings in regular expanding non-bipartite graphs, ITCS 2022, 61:1–61:12. https://doi.org/10.4230/LIPIcs.ITCS.2022.61
  • A. Abdi, G. Cornuéjols, D. Dadush, M. Dalirrooyfard, Lower bounds for cube-ideal set-systems, Proc. London Math. Soc. 133 (2026), e70199. https://doi.org/10.1112/plms.70199
  • OpenAI, Entropy and Face Dimension of the Perfect-Matching Polytope, OpenAI Math Release preprint, September 23, 2026 (source of the goal; Theorem 1.1, p. 1). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Entropy-and-Face-Dimension-of-the-Perfect-Matching-Polytope-September-23-2026/main.pdf
2 thms1 active userReviewed
CombinatoricsMarkov Chain·Captain: wurtle

A Fully Polynomial Randomized Approximation Scheme for Perfect Matchings in General GraphsResearch Paper

Motivation: approximately counting perfect matchings

A perfect matching of a graph is a set of edges covering every vertex exactly once. Finding one is easy — Edmonds (1965) gave a polynomial-time algorithm for general graphs — but counting them is hard: Valiant (1979) proved exact counting #P-complete already for bipartite graphs, where the count is the permanent of a 0–1 matrix. Exact counting is tractable only in special classes, such as planar graphs via Pfaffians (Kasteleyn 1963). The natural relaxation is a fully polynomial randomized approximation scheme (FPRAS): an algorithm that, with probability at least 1−δ1-\delta1−δ, returns the count within relative error ε\varepsilonε, in time polynomial in the input size, 1/ε1/\varepsilon1/ε and log⁡(1/δ)\log(1/\delta)log(1/δ). Whether perfect matchings in general (nonbipartite) graphs admit an FPRAS was explicitly raised by Jerrum and Sinclair in 1989 and has been a benchmark question in approximate counting since. Perfect-matching counts are partition functions of the monomer–dimer model, so the question also connects to statistical physics.

Timeline

  • 1963 — Kasteleyn: exact counting of perfect matchings in planar graphs (J. Math. Phys. 1963).
  • 1965 — Edmonds: polynomial-time maximum matching in general graphs (Canad. J. Math. 1965).
  • 1979 — Valiant: computing the permanent is #P-complete (TCS 1979).
  • 1986 — Jerrum, Valiant and Vazirani: equivalence of approximate counting and almost-uniform sampling for self-reducible problems (TCS 1986).
  • 1989 — Jerrum and Sinclair: FPRAS for the number of all matchings in any graph, and for perfect matchings when the near-perfect/perfect ratio is polynomially bounded; they raise the general question (SIAM J. Comput. 1989, Section 7(i)).
  • 2004 — Jerrum, Sinclair and Vigoda: FPRAS for the permanent of every nonnegative matrix, i.e. bipartite perfect matchings (J. ACM 2004); accelerated by Bezáková, Štefankovič, Vazirani and Vigoda (SIAM J. Comput. 2008).
  • 2018 — Štefankovič, Vigoda and Wilmes: graphs on which every chain of Jerrum–Sinclair–Vigoda type either gives perfect matchings exponentially small mass or mixes exponentially slowly (LATIN 2018).
  • 2020–2025 — Cai and Liu relate the problem to the eight-vertex model (ICALP 2020); Fei, Goldberg and Lu locate it among two-state spin systems (Inf. Comput. 2025).
  • 2026 — Further permanent speedups by Chen, Vigoda and Yang (arXiv:2608.26599) and Chen, Guo, Vigoda and Yang (arXiv:2609.20717); deterministic hafnian approximation for dense graphs by Yi (arXiv:2609.04079). An OpenAI preprint, A Fully Polynomial Randomized Approximation Scheme for Perfect Matchings in General Graphs (OpenAI Math Release, September 23, 2026), claims an FPRAS for every finite simple graph. The preprint has not been peer reviewed, and its main theorem is not formally verified.

Setting

A finite simple undirected graph GGG on vertices {0,…,n−1}\{0,\dots,n-1\}{0,…,n−1} is given by a set of edges, each written once as a pair (u,v)(u,v)(u,v) with u<vu<vu<v. A set M⊆EM\subseteq EM⊆E is a perfect matching if every vertex lies in exactly one edge of MMM, and

Z(G)=#{M⊆E:M is a perfect matching of G}.Z(G)=\#\{M\subseteq E: M\text{ is a perfect matching of }G\}.Z(G)=#{M⊆E:M is a perfect matching of G}.

(The graph with no vertices has Z=1Z=1Z=1; any graph with an odd number of vertices has Z=0Z=0Z=0.)

A randomized algorithm is modelled as a single fixed machine with a finite transition table over a finite tape alphabet, which at each step reads the current tape symbol and one fresh fair random bit, then writes, moves, or halts. The input is the graph in binary (number of vertices, number of edges, the sorted edge list) followed by the rationals ε\varepsilonε and δ\deltaδ in reduced form; the output is a rational number written on the tape when the machine halts.

Formalization targets

Goal: Theorem 1.1 (FPRAS for perfect matchings)

There are a machine AAA and constants C>0C>0C>0, ddd such that for every graph GGG and rationals 0<ε<10<\varepsilon<10<ε<1, 0<δ<1/20<\delta<1/20<δ<1/2, with t=C (∣input∣+⌈ε−1⌉+⌈log⁡2⌈δ−1⌉⌉+1)dt=C\,(|\mathrm{input}|+\lceil\varepsilon^{-1}\rceil+\lceil\log_2\lceil\delta^{-1}\rceil\rceil+1)^dt=C(∣input∣+⌈ε−1⌉+⌈log2​⌈δ−1⌉⌉+1)d:

  • on every random string of length ttt, AAA halts within ttt steps with a nonnegative rational output Z^\widehat ZZ;
  • if Z(G)=0Z(G)=0Z(G)=0, the output is 000 on every random string;
Pr⁡bits[(1−ε)Z(G)≤Z^≤(1+ε)Z(G)]  ≥  1−δ.\Pr_{\text{bits}}\big[(1-\varepsilon)Z(G)\le\widehat Z\le(1+\varepsilon)Z(G)\big]\;\ge\;1-\delta.bitsPr​[(1−ε)Z(G)≤Z≤(1+ε)Z(G)]≥1−δ.

The goal statement is published on the platform with status Open: no machine-checked proof exists yet.

Significance

The result itself. The theorem answers the Jerrum–Sinclair question affirmatively: perfect matchings in arbitrary graphs can be approximately counted as efficiently, up to polynomial factors, as in the bipartite case. The worst-case time bound holds for every execution, not just in expectation. Via standard reductions, the source also derives counting schemes for edge subsets with prescribed vertex degrees and for matchings of a fixed size (Corollary 10.1). The accompanying samplers for those families always return a feasible object and approximate the uniform law in total variation.

Formalizing it. The statement fixes a concrete machine model, input encoding and time bound, so a formal proof would certify both correctness of the estimator and the polynomial running time of an explicit algorithm — a level of rigor rarely reached for Markov-chain Monte Carlo results, whose analyses involve many interacting constants.

Difficulty

The natural approach extends the Jerrum–Sinclair–Vigoda Markov chain on perfect and near-perfect matchings, with weights learned while edge activities are gradually lowered. In general graphs this approach provably fails: Štefankovič, Vigoda and Wilmes constructed graphs on which every chain of that type, for any hole-dependent weights, has either exponentially small stationary mass on perfect matchings or exponentially slow mixing. Odd cycles (blossoms) are the reason the bipartite canonical-path arguments do not carry over. A solution needs a different state space and a different mixing analysis.

Formalization scope

  • GraphInput is a vertex count n and a Finset (Fin n × Fin n) of edges with strictly increasing endpoints, so loops and duplicate edges are excluded. Perfect G M requires M ⊆ edges and that each vertex is incident to exactly one edge of M; Z G is the number of such M.
  • The machine model is RandomMachine: a finite transition table on Fin 8 symbols, whose step consumes one random bit and performs a Turing.TM0 write or move, or halts. run folds tick over a list of t bits; halted configurations are absorbing.
  • encodeInput writes the vertex count, edge count and lexicographically sorted edges in binary, then encodeRat ε and encodeRat δ (reduced numerator with sign, denominator). Outputs requires that the machine has halted and the tape to the right of the head is exactly the encoding of the output rational.
  • The time bound is a fixed polynomial with natural ceilings; the success probability is the fraction of good tapes among all 2t2^t2t tapes.
  • The scheme must be uniform: one machine for all inputs, quantified before the graph. A trivializing reading is excluded because every tape must halt within the bound with a nonnegative output, and zero must be returned exactly when Z(G)=0Z(G)=0Z(G)=0.
  • Infrastructure needed: a compiler-style layer turning structured randomized algorithms into RandomMachine transitions, Markov-chain spectral-gap bounds, and sampling-to-counting reductions. The machine model and encodings are reusable for other approximate-counting missions.

Selected references

  • P. W. Kasteleyn, Dimer statistics and phase transitions, J. Math. Phys. 4 (1963), 287–293. https://doi.org/10.1063/1.1703953
  • J. Edmonds, Paths, trees, and flowers, Canad. J. Math. 17 (1965), 449–467. https://doi.org/10.4153/CJM-1965-045-4
  • L. G. Valiant, The complexity of computing the permanent, Theoret. Comput. Sci. 8 (1979), 189–201. https://doi.org/10.1016/0304-3975(79)90044-6
  • M. R. Jerrum, L. G. Valiant, V. V. Vazirani, Random generation of combinatorial structures from a uniform distribution, Theoret. Comput. Sci. 43 (1986), 169–188. https://doi.org/10.1016/0304-3975(86)90174-X
  • M. Jerrum, A. Sinclair, Approximating the permanent, SIAM J. Comput. 18 (1989), 1149–1178. https://doi.org/10.1137/0218077
  • M. Jerrum, A. Sinclair, E. Vigoda, A polynomial-time approximation algorithm for the permanent of a matrix with nonnegative entries, J. ACM 51 (2004), 671–697. https://doi.org/10.1145/1008731.1008738
  • I. Bezáková, D. Štefankovič, V. V. Vazirani, E. Vigoda, Accelerating simulated annealing for the permanent and combinatorial counting problems, SIAM J. Comput. 37 (2008), 1429–1454. https://doi.org/10.1137/050644033
  • D. Štefankovič, E. Vigoda, J. Wilmes, On counting perfect matchings in general graphs, LATIN 2018, LNCS 10807, 873–885. https://doi.org/10.1007/978-3-319-77404-6_63
  • J.-Y. Cai, T. Liu, Counting perfect matchings and the eight-vertex model, ICALP 2020, LIPIcs 168, 23:1–23:18. https://doi.org/10.4230/LIPIcs.ICALP.2020.23
  • Y. Fei, L. A. Goldberg, P. Lu, Two-state spin systems with negative interactions, Inf. Comput. 307 (2025), 105340. https://doi.org/10.1016/j.ic.2025.105340
  • OpenAI, A Fully Polynomial Randomized Approximation Scheme for Perfect Matchings in General Graphs, OpenAI Math Release preprint, September 23, 2026 (source of the goal; Theorem 1.1, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/A-Fully-Polynomial-Randomized-Approximation-Scheme-for-Perfect-Matchings-in-General-Graphs-September-23-2026/main.pdf
2 thms1 active userReviewed
Complexity Theory·Captain: wurtle

Beyond the Square-Root Exponent for Depth-Three Boolean CircuitsResearch Paper

Motivation

Proving that explicit Boolean functions need large circuits is the central problem of circuit complexity. Even for very shallow circuits the known bounds are far from what counting arguments suggest. A depth-three OR–AND–OR circuit is an OR of CNFs, with unbounded fan-in and arbitrary sharing; most functions need 2Ω(n)2^{\Omega(n)}2Ω(n) gates in this model, yet for every explicit function in P the best known lower bounds had the form 2cn2^{c\sqrt n}2cn​ for a constant ccc. Håstad's switching lemma gave 2Ω(n)2^{\Omega(\sqrt n)}2Ω(n​) for parity (1989), and Paturi, Pudlák and Zane showed this is tight for parity up to a constant factor (1999), so parity cannot go further. Moving the exponent beyond any constant multiple of n\sqrt nn​ for an explicit language was identified as an open problem in the depth-three literature (Gurumukhani, Paturi, Pudlák, Saks and Talebanfard, CCC 2024; Gurumukhani, Kleber, Paturi, Rosin and Talebanfard, CCC 2026). Depth-three lower bounds are also tied to satisfiability algorithms and to Valiant's approach to lower bounds for logarithmic-depth circuits.

This mission asks for a formal proof that one polynomial-time language requires 2ω(n)2^{\omega(\sqrt n)}2ω(n​) depth-three gates, as stated in an OpenAI preprint dated September 23, 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

  • 1989 — Håstad's switching lemma yields 2Ω(n)2^{\Omega(\sqrt n)}2Ω(n​) depth-three lower bounds for parity (Almost optimal lower bounds for small depth circuits, Adv. Comput. Res. 5, 1989).
  • 1995 — Håstad, Jukna and Pudlák's top-down method improves the constants for parity and majority (Comput. Complexity 1995) and asks for 2n1/2+ε2^{n^{1/2+\varepsilon}}2n1/2+ε.
  • 1999–2000 — Paturi, Pudlák and Zane prove the sharp parity bound via the satisfiability coding lemma (Chicago J. TCS 1999); Paturi, Saks and Zane prove 2n−o(n)2^{n-o(n)}2n−o(n) bounds for bottom fan-in two (Comput. Complexity 2000).
  • 2001 — Impagliazzo, Paturi and Zane's sunflower sparsification (JCSS 2001).
  • 2005 — Paturi, Pudlák, Saks and Zane's code-membership bound with exponent (π/6)n(\pi/\sqrt6)\sqrt n(π/6​)n​ (JACM 2005).
  • 2024–2026 — Local enumeration and bottom-fan-in-three bounds (CCC 2024); optimal depth-three circuits for inner product (CCC 2026).
  • September 2026 — The OpenAI preprint claims a polynomial-time language with S3(fn)=2ω(n)S_3(f_n)=2^{\omega(\sqrt n)}S3​(fn​)=2ω(n​) (Theorem 1.1, p. 1).

Setting

A literal is a variable xix_ixi​ or its negation; negations are free. A depth-three OR–AND–OR circuit on nnn inputs has bottom OR gates (clauses) reading literals and optional constants, middle AND gates reading bottom gates, and one top OR gate reading middle gates. Fan-in and fan-out are unbounded and gates may be shared. Its size counts every bottom gate, every middle gate and the output gate. For f:{0,1}n→{0,1}f:\{0,1\}^n\to\{0,1\}f:{0,1}n→{0,1}, S3(f)S_3(f)S3​(f) is the least size of such a circuit computing fff exactly. For a language L⊆{0,1}∗L\subseteq\{0,1\}^*L⊆{0,1}∗, fnf_nfn​ is its membership function on {0,1}n\{0,1\}^n{0,1}n.

Formalization targets

Goal: Theorem 1.1 (p. 1)

There are a language L⊆{0,1}∗L\subseteq\{0,1\}^*L⊆{0,1}∗, a deterministic multi-tape Turing machine deciding LLL, and positive integers C,aC,aC,a such that the machine halts within C(n+1)aC(n+1)^aC(n+1)a steps on every input of length nnn, and for every real A>0A>0A>0 there is NAN_ANA​ with

S3(fn) > 2Anfor all n≥NA.S_3(f_n)\ >\ 2^{A\sqrt n}\qquad\text{for all } n\ge N_A .S3​(fn​) > 2An​for all n≥NA​.

The language and machine are fixed before AAA; no uniformity or bottom fan-in restriction is imposed on the circuits.

Significance

The result itself. It is the first explicit depth-three lower bound whose exponent is not O(n)O(\sqrt n)O(n​), crossing the barrier at which parity, majority and the earlier code-membership functions stop. It is weaker than the 2n1/2+ε2^{n^{1/2+\varepsilon}}2n1/2+ε bound asked by Håstad, Jukna and Pudlák and gives no fixed improvement of the exponent's power. Because every function has a single-CNF representation when bottom fan-in is unbounded, the result depends on counting all gates, as stated.

Formalizing it. Circuit lower bounds are combinatorial statements well suited to formal proof, but the libraries are thin: a complete development needs a Turing-machine time bound for a concrete polynomial-evaluation language, restriction arguments for CNFs, the Erdős–Rado sunflower lemma, and limited-independence estimates from hashing and polynomial interpolation. All of these are reusable in complexity theory. No machine-checked depth-three lower bound beyond small cases is known.

Difficulty

The classical switching lemma changes the normal form and loses a constant in the exponent, which keeps bounds at the cnc\sqrt ncn​ scale, and parity is provably no harder than 2O(n)2^{O(\sqrt n)}2O(n​), so a different hard function and a different argument are needed. The preprint needs a function whose sign has small correlation with every moderately wide CNF, independently of the number of clauses. Its key step (Lemma 3.1, p. 4) expresses each randomly restricted width-kkk CNF as a signed combination of width-bbb CNFs with expected coefficient mass at most two, with no dependence on the clause count. A uniform polynomial-time language cannot search for a hard function at each length, so hardness must come from inputs that carry the parameters selecting a hard slice (Section 5).

Formalization scope

  • Circuit3 (Fin n) has bottomCount raw clauses (lists of literals or Boolean constants), middleCount AND gates given as Finsets of bottom indices, and a top Finset of middle indices; gateCount = bottomCount + middleCount + 1. Computes requires exact agreement with L (List.ofFn x) on every input.
  • The machine model FiniteMultiTapeMachine has finitely many tapes, a finite alphabet, finite states, an input tape initialized with the input, and an accept predicate on the halting state. MultiTapeHaltsIn M w (L w) (C*(|w|+1)^a) requires halting within the bound with output L w.
  • The order of quantifiers is ∃L,M,C,a\exists L, M, C, a∃L,M,C,a before ∀A>0 ∃N ∀n≥N ∀\forall A>0\ \exists N\ \forall n\ge N\ \forall∀A>0 ∃N ∀n≥N ∀ circuits, with a strict inequality 2An<gateCount2^{A\sqrt n}<\mathrm{gateCount}2An​<gateCount in ℝ.
  • The statement is not vacuous: circuits exist for every function, so the claim is a genuine lower bound on all of them.

Selected references

  • OpenAI, Beyond the Square-Root Exponent for Depth-Three Boolean Circuits, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Beyond-the-Square-Root-Exponent-for-Depth-Three-Boolean-Circuits-September-23-2026/main.pdf
  • J. Håstad, Almost optimal lower bounds for small depth circuits, Randomness and Computation, Adv. Comput. Res. 5 (1989), 143–170.
  • J. Håstad, S. Jukna, P. Pudlák, Top-down lower bounds for depth-three circuits, Comput. Complexity 5 (1995), 99–112. https://doi.org/10.1007/BF01268140
  • R. Paturi, P. Pudlák, F. Zane, Satisfiability coding lemma, Chicago J. Theor. Comput. Sci. (1999), Article 11.
  • R. Paturi, M. E. Saks, F. Zane, Exponential lower bounds for depth three Boolean circuits, Comput. Complexity 9 (2000), 1–15. https://doi.org/10.1007/PL00001598
  • R. Impagliazzo, R. Paturi, F. Zane, Which problems have strongly exponential complexity?, J. Comput. Syst. Sci. 63 (2001), 512–530. https://doi.org/10.1006/jcss.2001.1774
  • R. Paturi, P. Pudlák, M. E. Saks, F. Zane, An improved exponential-time algorithm for k-SAT, J. ACM 52 (2005), 337–364. https://doi.org/10.1145/1066100.1066101
  • M. Gurumukhani, R. Paturi, P. Pudlák, M. Saks, N. Talebanfard, Local enumeration and majority lower bounds, CCC 2024. https://doi.org/10.4230/LIPIcs.CCC.2024.17
  • M. Gurumukhani, D. Kleber, R. Paturi, C. Rosin, N. Talebanfard, Optimal depth-three circuits for inner product, CCC 2026. https://doi.org/10.4230/LIPIcs.CCC.2026.27
  • H. Dell, T. Husfeldt, D. Marx, N. Taslaman, M. Wahlén, Exponential time complexity of the permanent and the Tutte polynomial, ACM Trans. Algorithms 10 (2014). https://doi.org/10.1145/2635812
  • P. Erdős, R. Rado, Intersection theorems for systems of sets, J. London Math. Soc. (1960). https://doi.org/10.1112/jlms/s1-35.1.85
2 thms1 active userReviewed
Complexity TheoryOperations Research·Captain: wurtle

Uniform computation of the squared-logarithmic k-server boundResearch Paper

Motivation

In the kkk-server problem, kkk labelled servers occupy points of a metric space, requests arrive online, and each request is served by moving one server onto it; the cost is the total distance moved, compared with the offline optimum through the competitive ratio. A companion OpenAI preprint claims that on every metric there exists a randomized policy with ratio O(log⁡2(k+1))O(\log^2(k+1))O(log2(k+1)) against oblivious requests, matching the known worst-case lower bound. That result is an existence statement: it gives no procedure for computing the policy's probabilities.

This mission asks whether the bound can be realized by one uniform algorithm with explicit bit complexity on all finite rational metrics, 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

  • 1988–1990 — Manasse, McGeoch and Sleator introduce the kkk-server problem and prove the deterministic lower bound kkk (STOC 1988; J. Algorithms 1990).
  • 1991 — Paging (the uniform metric): Fiat, Karp, Luby, McGeoch, Sleator and Young give the randomized marking algorithm and the harmonic lower bound (J. Algorithms 1991); McGeoch and Sleator attain HkH_kHk​ exactly (Algorithmica 1991). This suggested the randomized kkk-server conjecture, an O(log⁡k)O(\log k)O(logk) ratio on every metric.
  • 1995 — Koutsoupias and Papadimitriou prove that the work function algorithm is (2k−1)(2k-1)(2k−1)-competitive (J. ACM 1995).
  • 1996–2004 — Probabilistic tree embeddings: Bartal (FOCS 1996) and Fakcharoenphol–Rao–Talwar with distortion O(log⁡n)O(\log n)O(logn) (JCSS 2004).
  • 2012–2015 — Bansal, Buchbinder and Naor give O(log⁡k)O(\log k)O(logk) for weighted paging (J. ACM 2012); Bansal, Buchbinder, Mądry and Naor give the first polylogarithmic bound O(log⁡2klog⁡3nlog⁡log⁡n)O(\log^2k\log^3n\log\log n)O(log2klog3nloglogn) on arbitrary finite metrics (J. ACM 2015).
  • 2018 — Bubeck, Cohen, Lee, Lee and Mądry achieve O(log⁡2k)O(\log^2 k)O(log2k) on hierarchically separated trees, hence O(log⁡2klog⁡n)O(\log^2k\log n)O(log2klogn) on nnn-point metrics (STOC 2018).
  • 2023 — Bubeck, Coester and Rabani disprove the randomized kkk-server conjecture: some (k+1)(k+1)(k+1)-point metrics need ratio Ω(log⁡2k)\Omega(\log^2k)Ω(log2k) (STOC 2023).
  • 2026 — Coester and Cosson make the tree and finite-metric bounds polynomial-time, with ratio O(log⁡nlog⁡2k)O(\log n\log^2k)O(lognlog2k) on general metrics (ICALP 2026).
  • September 2026 — Two OpenAI preprints claim the matching upper bound O(log⁡2(k+1))O(\log^2(k+1))O(log2(k+1)) on every metric space (source) and a uniform bit-level implementation on finite rational metrics (source).

Setting

Fix n≥3n\ge3n≥3 and 2≤k<n2\le k<n2≤k<n. A rational metric on X={1,…,n}X=\{1,\dots,n\}X={1,…,n} is a table d(x,y)∈Qd(x,y)\in\mathbb Qd(x,y)∈Q that is nonnegative, zero exactly on the diagonal, symmetric, and satisfies the triangle inequality. An initial labelled tuple s∈Xks\in X^ks∈Xk may have repetitions. For a request word σ\sigmaσ, costA,s(σ)\mathrm{cost}_{A,s}(\sigma)costA,s​(σ) is the algorithm's total movement and OPTs(σ)\mathrm{OPT}_s(\sigma)OPTs​(σ) the minimum over all choices of serving labels with knowledge of σ\sigmaσ. LLL denotes the total binary length of the encoding of (n,k,d,s)(n,k,d,s)(n,k,d,s).

The computational model is a bit machine: a finite control with an input tape, finitely many work tapes over {0,1,blank}\{0,1,\text{blank}\}{0,1,blank}, an output channel, and one fresh unbiased random bit per step. At each request the input tape is replaced by the encoded request, the internal state is retained, and the output after the step budget names the server to move.

Formalization targets

Goal: Theorem 1.1 (p. 1)

There are an absolute constant CCC and a polynomial ppp such that a single uniform randomized online algorithm, on every instance (n,k,d,s)(n,k,d,s)(n,k,d,s) as above,

  • preprocesses the instance in at most p(L)p(L)p(L) bit operations;
  • serves request rtr_trt​, for every reachable history and internal state, in at most p(L+⌈log⁡2(t+1)⌉)p\bigl(L+\lceil\log_2(t+1)\rceil\bigr)p(L+⌈log2​(t+1)⌉) further bit operations, outputting a valid server label and moving exactly that server;
  • for every fixed finite request sequence σ\sigmaσ independent of its random bits satisfies
E costA,s(σ) ≤ C (log⁡(k+1))2 OPTs(σ)+Bd,k,s,\mathbb E\,\mathrm{cost}_{A,s}(\sigma)\ \le\ C\,\bigl(\log(k+1)\bigr)^2\,\mathrm{OPT}_s(\sigma)+B_{d,k,s},EcostA,s​(σ) ≤ C(log(k+1))2OPTs​(σ)+Bd,k,s​,

with Bd,k,s<∞B_{d,k,s}<\inftyBd,k,s​<∞ independent of σ\sigmaσ and of its length.

Significance

The result itself. It converts the existence theorem into an algorithm whose ratio is independent of nnn, using only unbiased random bits and polynomial work per request in the input length and the bit length of the request counter. The price is an additive constant Bd,k,sB_{d,k,s}Bd,k,s​ with no size bound: the algorithm may serve a very long initial stretch with a fixed label while its constructor runs. Coester and Cosson obtain stronger time and randomness guarantees but with ratio O(log⁡nlog⁡2k)O(\log n\log^2k)O(lognlog2k) on general metrics; the two results are incomparable.

Formalizing it. The goal is stated in an explicit machine model, so a formal proof certifies both the competitive bound and the bit-complexity bound. It needs the companion existence theorem (Theorem 2.1 here, p. 4) as an input; a full development therefore also requires formalizing that theorem or taking it as a separate milestone. Exact rational linear algebra (Fourier–Motzkin elimination), dyadic rounding of probabilities, and step-counted machine simulation are reusable.

Difficulty

Knowing that a good policy exists does not give its probabilities, and a finite horizon alone does not control a long input on which the optimum barely moves. The preprint (pp. 3–4) solves a rational linear system for each horizon, then filters requests so that only those missed by some strongly lazy service with bounded moves are retained, which bounds the retained word length; a marking fallback pays for trajectories exceeding the move cap. The constructor's running time is uncontrolled, so the algorithm simulates it for only ⌊log⁡2(t+1)⌋\lfloor\log_2(t+1)\rfloor⌊log2​(t+1)⌋ steps at request ttt (Proposition 4.2, p. 14), and the delayed activation is absorbed into Bd,k,sB_{d,k,s}Bd,k,s​.

Formalization scope

  • RationalMetric n is a Fin n → Fin n → ℚ table with the metric axioms; configurations are Fin k → Fin n; offlineCost is the sInf over label histories matching the request word (a finite, nonempty set of costs).
  • BitMachine has finitely many controls and tapes and a transition consuming one coin per step. Preprocessing (boot) runs c(L+1)ec(L+1)^ec(L+1)e steps with all coins 000, so it is deterministic; this is satisfied by the source's deterministic constructor, and Theorem 1.1 does not require randomness there.
  • Each request runs exactly requestBudget c e L t =c(L+⌈log⁡2(t+1)⌉+1)e=c(L+\lceil\log_2(t+1)\rceil+1)^e=c(L+⌈log2​(t+1)⌉+1)e steps (Nat.clog 2 (t+1)) and must have yielded with output value <k<k<k for every reachable state and every coin string; requests are encoded by a fixed self-delimiting binary code.
  • machineExpectedCost averages over all coin strings of the step budget, i.e. uniformly random bits; the bound is required for every request list w, with B≥0B\ge0B≥0 chosen after (d,s)(d,s)(d,s) and before w.
  • The machine, the polynomial (c,e)(c,e)(c,e) and CCC are fixed before nnn, kkk, ddd, sss: the algorithm is uniform.

Selected references

  • OpenAI, Uniform computation of the squared-logarithmic k-server bound, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Uniform-computation-of-the-squared-logarithmic-k-server-bound-September-24-2026/Uniform-computation-of-the-squared-logarithmic-k-server-bound-September-24-2026.pdf
  • OpenAI, Squared-logarithmic randomized k-server on arbitrary metrics, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Squared-logarithmic-randomized-k-server-on-arbitrary-metrics-September-24-2026/Squared-logarithmic-randomized-k-server-on-arbitrary-metrics-September-24-2026.pdf
  • D. Komm, R. Královič, R. Královič, T. Mömke, Randomized online computation with high probability guarantees, Algorithmica (2022). https://doi.org/10.1007/s00453-022-00925-z
  • G. B. Dantzig, B. C. Eaves, Fourier–Motzkin elimination and its dual, J. Combin. Theory Ser. A 14 (1973). https://doi.org/10.1016/0097-3165(73)90004-6
  • M. S. Manasse, L. A. McGeoch, D. D. Sleator, Competitive algorithms for server problems, J. Algorithms 11 (1990). https://doi.org/10.1016/0196-6774(90)90003-W
  • E. Koutsoupias, C. H. Papadimitriou, On the k-server conjecture, J. ACM 42 (1995). https://doi.org/10.1145/210118.210128
  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive paging algorithms, J. Algorithms 12 (1991). https://doi.org/10.1016/0196-6774(91)90041-V
  • N. Bansal, N. Buchbinder, A. Mądry, J. Naor, A polylogarithmic-competitive algorithm for the k-server problem, J. ACM 62 (2015). https://doi.org/10.1145/2783434
  • S. Bubeck, M. B. Cohen, J. R. Lee, Y. T. Lee, A. Mądry, k-server via multiscale entropic regularization, STOC 2018. https://doi.org/10.1145/3188745.3188798
  • S. Bubeck, C. Coester, Y. Rabani, The randomized k-server conjecture is false!, STOC 2023. https://doi.org/10.1145/3564246.3585132
  • C. Coester, R. Cosson, Randomized k-server in polynomial time, ICALP 2026. https://doi.org/10.4230/LIPIcs.ICALP.2026.65
  • J. Fakcharoenphol, S. Rao, K. Talwar, A tight bound on approximating arbitrary metrics by tree metrics, JCSS 69 (2004). https://doi.org/10.1016/j.jcss.2004.04.011
2 thms1 active userReviewed
Operations Research·Captain: wurtle

Squared-logarithmic randomized k-server on arbitrary metricsResearch Paper

Motivation

In the kkk-server problem, kkk servers sit at points of a metric space; requests arrive one at a time, and each request must be served by moving some server onto it, at a cost equal to the distance travelled. An online policy sees only the requests revealed so far, and its quality is measured by its competitive ratio: how much more it moves, in expectation, than the offline optimum that knows the whole sequence. The problem was introduced as a common framework for paging, caching and other online service problems, and the question of the best ratio achievable by randomized policies has driven the development of metric embeddings, online primal-dual methods and entropic regularization.

This mission asks for a formal proof that randomized kkk-server has competitive ratio O(log⁡2(k+1))O(\log^2(k+1))O(log2(k+1)) on every metric space, 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

  • 1988–1990 — Manasse, McGeoch and Sleator introduce the kkk-server problem and prove the deterministic lower bound kkk (STOC 1988; J. Algorithms 1990).
  • 1991 — Paging (the uniform metric): Fiat, Karp, Luby, McGeoch, Sleator and Young give the randomized marking algorithm and the harmonic lower bound (J. Algorithms 1991); McGeoch and Sleator attain HkH_kHk​ exactly (Algorithmica 1991). This suggested the randomized kkk-server conjecture, an O(log⁡k)O(\log k)O(logk) ratio on every metric.
  • 1995 — Koutsoupias and Papadimitriou prove that the work function algorithm is (2k−1)(2k-1)(2k−1)-competitive (J. ACM 1995).
  • 1996–2004 — Probabilistic tree embeddings: Bartal (FOCS 1996) and Fakcharoenphol–Rao–Talwar with distortion O(log⁡n)O(\log n)O(logn) (JCSS 2004).
  • 2012–2015 — Bansal, Buchbinder and Naor give O(log⁡k)O(\log k)O(logk) for weighted paging (J. ACM 2012); Bansal, Buchbinder, Mądry and Naor give the first polylogarithmic bound O(log⁡2klog⁡3nlog⁡log⁡n)O(\log^2k\log^3n\log\log n)O(log2klog3nloglogn) on arbitrary finite metrics (J. ACM 2015).
  • 2018 — Bubeck, Cohen, Lee, Lee and Mądry achieve O(log⁡2k)O(\log^2 k)O(log2k) on hierarchically separated trees, hence O(log⁡2klog⁡n)O(\log^2k\log n)O(log2klogn) on nnn-point metrics (STOC 2018).
  • 2023 — Bubeck, Coester and Rabani disprove the randomized kkk-server conjecture: some (k+1)(k+1)(k+1)-point metrics need ratio Ω(log⁡2k)\Omega(\log^2k)Ω(log2k) (STOC 2023).
  • 2026 — Coester and Cosson make the tree and finite-metric bounds polynomial-time, with ratio O(log⁡nlog⁡2k)O(\log n\log^2k)O(lognlog2k) on general metrics (ICALP 2026).
  • September 2026 — Two OpenAI preprints claim the matching upper bound O(log⁡2(k+1))O(\log^2(k+1))O(log2(k+1)) on every metric space (source) and a uniform bit-level implementation on finite rational metrics (source).

Setting

Let (X,d)(X,d)(X,d) be a metric space and s=(s1,…,sk)∈Xks=(s_1,\dots,s_k)\in X^ks=(s1​,…,sk​)∈Xk the initial positions of kkk labelled servers; repeated positions are allowed. To serve a request rtr_trt​, a policy chooses one label jjj and moves server jjj to rtr_trt​, leaving the others fixed (a server already at rtr_trt​ may be chosen at zero cost). For a finite request sequence σ\sigmaσ, costA,s(σ)\mathrm{cost}_{A,s}(\sigma)costA,s​(σ) is the total distance moved and OPTs(σ)\mathrm{OPT}_s(\sigma)OPTs​(σ) is the minimum of this cost over all label sequences, chosen with knowledge of σ\sigmaσ. A randomized online policy chooses, at each step, a probability distribution over labels as a function of the past requests, its own past choices and the current request. Requests are fixed in advance, independently of the policy's randomness (the oblivious adversary).

Formalization targets

Goal: Theorem 1.1 (p. 2)

There is an absolute constant C<∞C<\inftyC<∞ such that for every k≥2k\ge2k≥2, every metric space (X,d)(X,d)(X,d) with at least k+1k+1k+1 points, and every s∈Xks\in X^ks∈Xk, there are a randomized online policy AAA and a finite B≥0B\ge0B≥0 with

E costA,s(σ) ≤ C (log⁡(k+1))2 OPTs(σ)+B\mathbb E\,\mathrm{cost}_{A,s}(\sigma)\ \le\ C\,\bigl(\log(k+1)\bigr)^2\,\mathrm{OPT}_s(\sigma)+BEcostA,s​(σ) ≤ C(log(k+1))2OPTs​(σ)+B

for every finite request sequence σ\sigmaσ. The same policy serves all sequences and horizons; BBB may depend on (X,d)(X,d)(X,d), kkk and sss but not on σ\sigmaσ; and B=0B=0B=0 when the entries of sss are distinct.

Significance

The result itself. Combined with the Bubeck–Coester–Rabani lower bound, it identifies Θ(log⁡2k)\Theta(\log^2 k)Θ(log2k) as the worst-case order of the randomized competitive ratio over all metrics. The bound has no dependence on the number of points, the aspect ratio, finiteness or boundedness of XXX, removing the log⁡n\log nlogn factor left by static tree embeddings. No computational efficiency is claimed; a companion preprint turns the existence statement into a uniform algorithm on finite rational metrics.

Formalizing it. The statement is elementary to state, but its proof combines a rank-based allocation on separated trees, slowly varying simplex trackers, a compact potential for partition edits, online partitions with an embedded comparator, balanced rounding, and a compactness argument (Tychonoff) producing one policy for all finite inputs. None of these online-algorithm tools exist in a formal library, and a machine-checked competitive analysis of a randomized kkk-server algorithm on general metrics is not known.

Difficulty

Embedding the metric into a random hierarchically separated tree and running an O(log⁡2k)O(\log^2k)O(log2k) tree algorithm loses the embedding distortion, which is Θ(log⁡n)\Theta(\log n)Θ(logn) in the worst case and is not bounded in terms of kkk. Removing it requires changing the tree as requests arrive, and the cost of re-partitioning must be paid for. A previous dynamic approach (Lee's fusible HSTs) had its general claim withdrawn after a gap in that accounting (p. 3). In the preprint, the partition schedule is driven by the posterior measure of a hidden comparator, and the edit costs and the travel of parked mass are bounded separately (Sections 7–8). Passing from finite metrics and finite horizons to one policy on an arbitrary, possibly infinite and unbounded, space with a fixed constant requires a further compactness step (Section 10).

Formalization scope

  • X : Type u carries an arbitrary MetricSpace instance; "at least k+1k+1k+1 points" is an injection Fin (k+1) → X. No finiteness, separability or boundedness is assumed.
  • A Policy k X maps the history (list of past requests with the labels used) and the current request to a LabelDistribution k, i.e. a probability vector on Fin k. Every step moves exactly one server (Function.update), as in the source.
  • expectedCost is the finite sum over all label sequences of the path probability times the service cost; optimalCost is the sInf over all label sequences of the service cost (nonempty finite set, so the infimum is a minimum). Restricting the offline optimum to one-server-per-request service does not change it (triangle inequality).
  • The constant C>0C>0C>0 is chosen before kkk, XXX and sss; the policy and B≥0B\ge0B≥0 are chosen after sss and before the request list; Function.Injective s → B = 0 encodes the distinct-start clause.
  • Logarithms are natural (Real.log).

Selected references

  • OpenAI, Squared-logarithmic randomized k-server on arbitrary metrics, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Squared-logarithmic-randomized-k-server-on-arbitrary-metrics-September-24-2026/Squared-logarithmic-randomized-k-server-on-arbitrary-metrics-September-24-2026.pdf
  • OpenAI, Uniform computation of the squared-logarithmic k-server bound, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Uniform-computation-of-the-squared-logarithmic-k-server-bound-September-24-2026/Uniform-computation-of-the-squared-logarithmic-k-server-bound-September-24-2026.pdf
  • J. R. Lee, Fusible HSTs and the randomized k-server conjecture, FOCS 2018. https://doi.org/10.1109/FOCS.2018.00049
  • A. Tychonoff, Über einen Funktionenraum, Math. Ann. 111 (1935). https://doi.org/10.1007/BF01472255
  • M. S. Manasse, L. A. McGeoch, D. D. Sleator, Competitive algorithms for server problems, J. Algorithms 11 (1990). https://doi.org/10.1016/0196-6774(90)90003-W
  • E. Koutsoupias, C. H. Papadimitriou, On the k-server conjecture, J. ACM 42 (1995). https://doi.org/10.1145/210118.210128
  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive paging algorithms, J. Algorithms 12 (1991). https://doi.org/10.1016/0196-6774(91)90041-V
  • N. Bansal, N. Buchbinder, A. Mądry, J. Naor, A polylogarithmic-competitive algorithm for the k-server problem, J. ACM 62 (2015). https://doi.org/10.1145/2783434
  • S. Bubeck, M. B. Cohen, J. R. Lee, Y. T. Lee, A. Mądry, k-server via multiscale entropic regularization, STOC 2018. https://doi.org/10.1145/3188745.3188798
  • S. Bubeck, C. Coester, Y. Rabani, The randomized k-server conjecture is false!, STOC 2023. https://doi.org/10.1145/3564246.3585132
  • C. Coester, R. Cosson, Randomized k-server in polynomial time, ICALP 2026. https://doi.org/10.4230/LIPIcs.ICALP.2026.65
  • J. Fakcharoenphol, S. Rao, K. Talwar, A tight bound on approximating arbitrary metrics by tree metrics, JCSS 69 (2004). https://doi.org/10.1016/j.jcss.2004.04.011
2 thms1 active userReviewed
Algebraic GeometryComplexity Theory·Captain: wurtle

A cubic lower bound for border determinantal complexity of the permanentResearch Paper

Motivation

The permanent of an m×mm\times mm×m matrix X=(xij)X=(x_{ij})X=(xij​),

per⁡m(X)=∑σ∈Sm∏i=1mxi,σ(i),\operatorname{per}_m(X)=\sum_{\sigma\in S_m}\prod_{i=1}^m x_{i,\sigma(i)},perm​(X)=σ∈Sm​∑​i=1∏m​xi,σ(i)​,

has the same monomials as the determinant but no signs. The determinant can be computed in polynomial time; the permanent is #P-hard (Valiant, 1979). Valiant's algebraic theory makes this contrast precise: the determinant is complete for small formulas and branching programs, the permanent is complete for the class VNP, and VP ≠ VNP, the algebraic analogue of P ≠ NP, would follow if the permanent cannot be written as a determinant of polynomially larger size with affine-linear entries. Mulmuley and Sohoni's geometric complexity theory reformulates this through orbit closures, which leads to the border version of the question studied here.

Timeline

  • 1979. Valiant shows that formulas reduce to determinants and that the permanent is complete for its algebraic class (doi:10.1145/800135.804419), and that computing the permanent is #P-complete (doi:10.1016/0304-3975(79)90044-6).
  • 1987–1990. Linear lower bounds of order 2 m\sqrt2\,m2​m: von zur Gathen (Babai–Seress bound) (doi:10.1016/0024-3795(87)90337-5), Meshulam (doi:10.1016/0024-3795(89)90465-5), Cai (doi:10.1016/0890-5401(90)90036-H).
  • 2001 and 2008. Mulmuley and Sohoni propose geometric complexity theory (doi:10.1137/S009753970038715X, doi:10.1137/080718115).
  • 2004. Mignon and Ressayre prove dc(per⁡m)≥m2/2\mathrm{dc}(\operatorname{per}_m)\ge m^2/2dc(perm​)≥m2/2 in characteristic zero via Hessian rank (doi:10.1155/S1073792804142566).
  • 2010. Cai, Chen and Li extend the quadratic bound to odd characteristic (doi:10.1007/s00037-009-0284-2).
  • 2012. Grenet represents per⁡m\operatorname{per}_mperm​ as a determinant of order 2m−12^m-12m−1 (upper bound).
  • 2013. Landsberg, Manivel and Ressayre prove the border bound bdc(per⁡m)≥m2/2\mathrm{bdc}(\operatorname{per}_m)\ge m^2/2bdc(perm​)≥m2/2 (doi:10.4171/CMH/292).
  • 2017. Alper, Bogart and Velasco bound dc\mathrm{dc}dc via the codimension of the singular locus, which is only linear for the permanent (doi:10.1007/s10208-015-9300-x); Landsberg and Ressayre prove exponential bounds for representations with symmetry (doi:10.1016/j.difgeo.2017.03.017).

Until this work the best lower bounds for both exact and border determinantal complexity of the permanent were quadratic. The source of this mission, an OpenAI preprint dated September 24, 2026, claims a cubic lower bound in the border model.

Setting

For a complex polynomial f∈C[x1,…,xN]f\in\mathbb C[x_1,\dots,x_N]f∈C[x1​,…,xN​], the determinantal complexity dc(f)\mathrm{dc}(f)dc(f) is the least positive integer nnn with

f=det⁡(A0+∑j=1NxjAj),A0,…,AN∈Matn(C).f=\det\Bigl(A_0+\sum_{j=1}^N x_jA_j\Bigr),\qquad A_0,\dots,A_N\in\mathrm{Mat}_n(\mathbb C).f=det(A0​+j=1∑N​xj​Aj​),A0​,…,AN​∈Matn​(C).

The border determinantal complexity bdc(f)\mathrm{bdc}(f)bdc(f) is the least positive nnn such that fff is a coefficientwise limit of such determinants of size nnn (the size stays fixed along the sequence; the matrices are arbitrary). Clearly bdc(f)≤dc(f)\mathrm{bdc}(f)\le\mathrm{dc}(f)bdc(f)≤dc(f).

Formalization targets

Corollary 9.3 (milestone): smooth initial forms

If g∈C[x1,…,xd]g\in\mathbb C[x_1,\dots,x_d]g∈C[x1​,…,xd​], d≥2d\ge2d≥2, a∈Cda\in\mathbb C^da∈Cd, and the lowest nonzero homogeneous term GrG_rGr​ of g(a+x)g(a+x)g(a+x) has degree r≥2r\ge2r≥2 and defines a smooth projective hypersurface, then

bdc(g) ≥ (r−1)(d−1)4e,dc(g) ≥ (r−1)(d−1)4e.\mathrm{bdc}(g)\ \ge\ \frac{(r-1)(d-1)}{4e},\qquad \mathrm{dc}(g)\ \ge\ \frac{(r-1)(d-1)}{4e}.bdc(g) ≥ 4e(r−1)(d−1)​,dc(g) ≥ 4e(r−1)(d−1)​.

Goal: Theorem 1.1 and Corollary 9.1

For every m≥1408m\ge1408m≥1408,

bdc(per⁡m) ≥ m35 529 600 e,dc(per⁡m) ≥ m35 529 600 e.\mathrm{bdc}(\operatorname{per}_m)\ \ge\ \frac{m^3}{5\,529\,600\,e},\qquad \mathrm{dc}(\operatorname{per}_m)\ \ge\ \frac{m^3}{5\,529\,600\,e}.bdc(perm​) ≥ 5529600em3​,dc(perm​) ≥ 5529600em3​.

The constants are those stated in the paper and are not optimized. The Lean statement OAI.PermanentBorder.permanent_cubic_lower_bounds is open on the platform.

Significance

The theorem raises the exponent of the best lower bound for the permanent versus determinant problem from 222 (Mignon–Ressayre, Landsberg–Manivel–Ressayre) to 333, in the border model and hence also in the exact model. Through determinant reductions it also gives Ω(m3)\Omega(m^3)Ω(m3) lower bounds on the size of affine-linear algebraic branching programs for the permanent. It remains a polynomial bound: superpolynomial determinantal complexity of the permanent, which would separate VP from VNP-type classes, is still open. The smooth-form bound of Corollary 9.3 applies to any polynomial with a smooth initial form, not only the permanent.

The result is claimed in an OpenAI preprint; it has not been peer reviewed and no machine-checked proof exists. Formalizing it would certify the polar-degree counting argument and the coefficient-projection construction that extracts a smooth form of linear degree and quadratic dimension from the permanent.

Difficulty

The Hessian method stops at the quadratic scale: the Hessian of an order-nnn determinant has rank at most 2n2n2n at a singular matrix, while the permanent's Hessian has rank at most m2m^2m2, so it cannot give more than m2/2m^2/2m2/2. Singular-locus methods give only linear bounds for the permanent because its singular locus has small codimension. A cubic bound needs an invariant that grows multiplicatively in degree and dimension, and in the border model it must survive limits, so any intersection count used must be shown to persist for a single sufficiently close determinant.

Formalization scope

  • Polynomials are MvPolynomial σ ℂ; the permanent is permanentPolynomial m on variables Fin m × Fin m.
  • affineMatrixPolynomial n A₀ A is det⁡(A0+∑vxvAv)\det(A_0+\sum_v x_vA_v)det(A0​+∑v​xv​Av​) with complex matrices.
  • HasExactDeterminant n p: n>0n>0n>0 and ppp equals such a determinant. HasBorderDeterminant n p: n>0n>0n>0 and sequences A0(j),A(j)A_0(j),A(j)A0​(j),A(j) whose determinant polynomials converge to ppp coefficient by coefficient. Instead of defining bdc\mathrm{bdc}bdc as a minimum, the goal states that every admissible size nnn obeys the bound, which is equivalent.
  • For the milestone, affineSubstitution a I g is g(a+x)g(a+x)g(a+x), homogeneousComponent extracts GrG_rGr​, and smoothness is stated on the affine cone: the only common zero of GrG_rGr​ and all its partial derivatives is 000.
  • Real constants use Real.exp 1 for eee.

A complete development needs projective smoothness and polar intersections (Bézout-type counts with simple points), adjugate and kernel-line computations for determinant pencils, and dimension estimates for matrix fibers. Contributions formalizing Lemma 2.1 (affine substitution preserves the closure), Proposition 2.2 (equivalence of border models) or Theorem 2.3 are welcome.

Selected references

  • OpenAI, A cubic lower bound for border determinantal complexity of the permanent, preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/A-cubic-lower-bound-for-border-determinantal-complexity-of-the-permanent-September-24-2026/A-cubic-lower-bound-for-border-determinantal-complexity-of-the-permanent-September-24-2026.pdf
  • L. G. Valiant, Completeness classes in algebra, STOC 1979. https://doi.org/10.1145/800135.804419
  • T. Mignon, N. Ressayre, A quadratic bound for the determinant and permanent problem, IMRN, 2004. https://doi.org/10.1155/S1073792804142566
  • J. M. Landsberg, L. Manivel, N. Ressayre, Hypersurfaces with degenerate duals and the Geometric Complexity Theory Program, Comment. Math. Helv., 2013. https://doi.org/10.4171/CMH/292
  • J.-Y. Cai, X. Chen, D. Li, Quadratic lower bound for permanent vs. determinant in any characteristic, Comput. Complexity, 2010. https://doi.org/10.1007/s00037-009-0284-2
  • K. D. Mulmuley, M. Sohoni, Geometric Complexity Theory I, SIAM J. Comput., 2001. https://doi.org/10.1137/S009753970038715X
  • J. Alper, T. Bogart, M. Velasco, A lower bound for the determinantal complexity of a hypersurface, Found. Comput. Math., 2017. https://doi.org/10.1007/s10208-015-9300-x
  • R. Piene, Polar classes of singular varieties, Ann. Sci. ÉNS, 1978. https://doi.org/10.24033/asens.1346
2 thms1 active userReviewed
Complexity TheoryLinear algebra·Captain: wurtle

Staggered extraction for exact matrix multiplication over every fieldResearch Paper

Motivation: the matrix-multiplication exponent

Multiplying two n×nn\times nn×n matrices by the schoolbook method costs about n3n^3n3 operations. The matrix-multiplication exponent ωF\omega_FωF​ of a field FFF is the infimum of exponents τ\tauτ such that n×nn\times nn×n matrices over FFF can be multiplied with O(nτ+o(1))O(n^{\tau+o(1)})O(nτ+o(1)) additions, subtractions and multiplications. It governs the asymptotic cost of a central bilinear operation and of the many algorithms built from it (determinants, inversion, linear systems). Bounding ω\omegaω has driven algebraic complexity theory since 1969. In recent years the bounds have moved in the third and fourth decimal place, with each improvement coming from a finer analysis of the same tensor-power framework.

Background

  • 1969 — Strassen's seven-product algorithm gives ω≤log⁡27\omega\le\log_27ω≤log2​7 (Strassen 1969).
  • 1980–1987 — Approximate bilinear algorithms (Bini 1980), the asymptotic sum inequality (Schönhage 1981) and the laser method (Strassen 1987).
  • 1990 — Coppersmith and Winograd introduce the tensor family CWq\mathrm{CW}_qCWq​ and progression-free extraction (J. Symbolic Comput. 1990).
  • 2010–2014 — Higher tensor powers: Stothers (fourth power; with Davie, Proc. Roy. Soc. Edinburgh 2013), Vassilevska Williams (eighth power, STOC 2012), Le Gall (thirty-second power, ISSAC 2014).
  • 2023–2024 — The refined laser method gives ω<2.3728596\omega<2.3728596ω<2.3728596 (Alman–Vassilevska Williams, TheoretiCS 2024). Asymmetric hashing reduces combination loss (Duan–Wu–Zhou). Finer variable sharing follows (Vassilevska Williams–Xu–Xu–Zhou).
  • 2024–2026 — Successive asymmetric compatibility tests (Alman–Duan–Vassilevska Williams–Xu–Xu–Zhou, arXiv:2404.16349). Dupont et al. optimize the eighth-power instance of this framework to ω<2.371177\omega<2.371177ω<2.371177 (arXiv:2608.16884).
  • 2026 — An OpenAI preprint, Staggered extraction for exact matrix multiplication over every field (OpenAI Math Release, September 24, 2026), claims ωF<2.371054886006746\omega_F<2.371054886006746ωF​<2.371054886006746 for every fixed field FFF, including every positive characteristic. It has not been peer reviewed, and the claim has not been formally verified.

Setting

Fix a field FFF. An arithmetic program is a straight-line sequence of gates, each a constant from FFF, an input entry, or the sum, difference or product of two earlier gates. Its cost counts the addition, subtraction and multiplication gates. A program multiplies n×nn\times nn×n matrices if designated registers equal the entries of ABABAB for every pair of input matrices A,BA,BA,B over FFF. A real number τ\tauτ is admissible if for every ε>0\varepsilon>0ε>0 there is C>0C>0C>0 such that for every n≥1n\ge1n≥1 some correct program has cost at most C nτ+εC\,n^{\tau+\varepsilon}Cnτ+ε, and

ωF=inf⁡{τ∈R: τ admissible over F}.\omega_F=\inf\{\tau\in\mathbb R:\ \tau\text{ admissible over }F\}.ωF​=inf{τ∈R: τ admissible over F}.

In Lean this is MatrixAllFields.MatrixMultiplication.Arithmetic.omega F, defined for every F : Type u with [Field F].

Formalization targets

Goal: the exponent bound over every field (Theorem 1.1)

∀ F field:ωF<2.371054886006746.\forall\,F\ \text{field}:\qquad \omega_F<2.371054886006746 .∀F field:ωF​<2.371054886006746.

The Lean statement also contains the numerical comparison 2.371054886006746<2.3710562.371054886006746<2.3710562.371054886006746<2.371056, which carries no mathematical content (the source writes the comparison as <2.371054887<2.371054887<2.371054887). This is Theorem 1.1 of the source. The goal is published on the platform with status Open: no machine-checked proof exists yet.

Significance

The result itself. If the source's argument and certificate are correct, the theorem lowers the best stated bound on ω\omegaω from 2.3711772.3711772.371177 (Dupont et al.) to 2.3710552.3710552.371055, and it does so uniformly over all fields, with exact algorithms in every positive characteristic. Any improvement in ω\omegaω lowers the exponent of every algorithm reducible to matrix multiplication. The method — pooling directional capacities across extraction depths and staggering lots so that one lot sits at each splitting stage — shows that combination losses in the asymmetric framework can be traded across depths.

Formalizing it. Bounds of this type rest on long chains of tensor-power constructions, recovery arguments and large rational parameter tables, checked in the source by an accompanying verifier with interval certificates for logarithms. A machine-checked proof would certify the whole chain from the arithmetic model to the final rational inequality. It would also give Mathlib the asymptotic sum inequality and Coppersmith–Winograd machinery, which every future bound on ω\omegaω needs.

Difficulty

Every recent improvement fights extraction losses. Restricting a tensor power to a direct sum of independent matrix products wastes variables whenever whole blocks must be disjoint, while the final products only need disjoint subsets. Prescribed marginals may also allow several joint distributions. Earlier analyses take the minimum of three directional capacities separately at each depth, so a shortage on one side at one depth is lost. Using capacity across depths requires that variables deleted in earlier extractions (to resolve conflicts or impose type conditions) can still be recovered exactly, over the original field and with only subexponentially many extra copies. Positive characteristic rules out arguments that divide by small integers or rely on characteristic-zero identities.

Formalization scope

  • The model is division-free straight-line programs over an arbitrary field F : Type u, with free constants from FFF and unit cost for add, sub and mul. This matches the source's "scalar additions and multiplications".
  • AdmissibleExponent F τ asks, for each ε>0\varepsilon>0ε>0, for one constant C>0C>0C>0 that works for every n≥1n\ge1n≥1. omega F is the sInf of the admissible set, which is nonempty (schoolbook algorithm) and bounded below by the output size, so a junk value cannot make the strict inequality true.
  • The goal quantifies over every field, including finite fields, so a proof cannot use characteristic-zero arguments or embeddings into C\mathbb CC.
  • The second conjunct is a trivially true rational comparison and can be ignored.
  • Infrastructure needed: tensor rank and border rank over arbitrary fields, the asymptotic sum inequality, the Coppersmith–Winograd tensor and its powers, method-of-types counting, and verified rational/logarithm interval arithmetic for the parameter certificate. All of this is reusable for other matrix-multiplication results.

Selected references

  • V. Strassen, Gaussian elimination is not optimal, Numer. Math. 13 (1969), 354–356. https://doi.org/10.1007/BF02165411
  • D. Bini, Relations between exact and approximate bilinear algorithms. Applications, Calcolo 17 (1980), 87–97. https://doi.org/10.1007/BF02575865
  • A. Schönhage, Partial and total matrix multiplication, SIAM J. Comput. 10 (1981), 434–455. https://doi.org/10.1137/0210032
  • V. Strassen, Relative bilinear complexity and matrix multiplication, J. reine angew. Math. 375/376 (1987), 406–443. https://doi.org/10.1515/crll.1987.375-376.406
  • D. Coppersmith and S. Winograd, Matrix multiplication via arithmetic progressions, J. Symbolic Comput. 9 (1990), 251–280. https://doi.org/10.1016/S0747-7171(08)80013-2
  • A. M. Davie and A. J. Stothers, Improved bound for complexity of matrix multiplication, Proc. Roy. Soc. Edinburgh Sect. A 143 (2013), 351–369. https://doi.org/10.1017/S0308210511001648
  • V. Vassilevska Williams, Multiplying matrices faster than Coppersmith–Winograd, STOC 2012, 887–898. https://doi.org/10.1145/2213977.2214056
  • F. Le Gall, Powers of tensors and fast matrix multiplication, ISSAC 2014, 296–303. https://arxiv.org/abs/1401.7714v1
  • J. Alman and V. Vassilevska Williams, A refined laser method and faster matrix multiplication, TheoretiCS 3 (2024), Art. 21. https://doi.org/10.46298/theoretics.24.21
  • R. Duan, H. Wu and R. Zhou, Faster matrix multiplication via asymmetric hashing, arXiv:2210.10173. https://arxiv.org/abs/2210.10173v5
  • V. Vassilevska Williams, Y. Xu, Z. Xu and R. Zhou, New bounds for matrix multiplication: from alpha to omega, arXiv:2307.07970. https://arxiv.org/abs/2307.07970v2
  • J. Alman, R. Duan, V. Vassilevska Williams, Y. Xu, Z. Xu and R. Zhou, More asymmetry yields faster matrix multiplication, arXiv:2404.16349v3. https://arxiv.org/abs/2404.16349v3
  • E. Dupont et al., Improving the matrix multiplication exponent with modern optimization, arXiv:2608.16884 (2026). https://doi.org/10.48550/arXiv.2608.16884
  • OpenAI, Staggered extraction for exact matrix multiplication over every field, OpenAI Math Release preprint, September 24, 2026 (Theorem 1.1, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Staggered-extraction-for-exact-matrix-multiplication-over-every-field-September-24-2026/Staggered-extraction-for-exact-matrix-multiplication-over-every-field-September-24-2026.pdf
2 thms1 active userReviewed
Complexity TheoryGraph Theory·Captain: wurtle

Hardness of finding large independent sets in three-colorable graphsResearch Paper

Motivation: how hard is it to use the promise of three-colorability?

A graph is three-colorable if its vertices can be colored with three colors so that adjacent vertices get different colors. Deciding three-colorability is NP-complete. The approximate coloring problem asks how much a promise helps: given a graph promised to be three-colorable, can one efficiently find a coloring with a few more colors, or at least a large independent set? Every three-colorable graph on nnn vertices has an independent set of size n/3n/3n/3 (its largest color class). The best polynomial-time algorithms use n0.19…n^{0.19\ldots}n0.19… colors, while hardness results long reached only small palettes and needed extra conjectures for stronger gaps. Problems of this kind are central test cases for the PCP theorem, Label Cover and the theory of promise constraint satisfaction.

Background and timeline

  • 2000 — Khanna, Linial and Safra prove NP-hardness of 444-coloring three-colorable graphs (Combinatorica); Guruswami and Khanna give a different proof with bounded-degree variants (SIAM J. Discrete Math. 2004).
  • 2009 — Dinur, Mossel and Regev obtain three-colorable completeness with arbitrarily small independent-set soundness, conditional on a fish-shaped Label Cover conjecture (SIAM J. Comput.).
  • 2010 — Dinur, Khot, Perkins and Safra prove hardness of finding independent sets of density 1/91/91/9 in almost-three-colorable graphs (FOCS 2010).
  • 2020 — Guruswami and Sandeep show that the ddd-to-111 Games Conjecture with perfect completeness implies hardness of coloring three-colorable graphs with O(1)O(1)O(1) colors (ICALP 2020).
  • 2021 — Barto, Bulín, Krokhin and Opršal prove hardness of (2k−1)(2k-1)(2k−1)-coloring kkk-colorable graphs (five colors for k=3k=3k=3) via the algebraic theory of promise CSPs (J. ACM); Braverman, Khot, Lifshitz and Minzer obtain the independent-set gap from the Rich 222-to-111 Games Conjecture (FOCS 2021).
  • 2023 — Krokhin, Opršal, Wrochna and Živný give the arc-graph reduction to palette three (SIAM J. Comput.); Hecht, Minzer and Safra treat almost coloring almost-three-colorable graphs (APPROX/RANDOM 2023).
  • 2024–2026 — Kawarabayashi, Thorup and Yoneda (STOC 2024) and Bansal, Huang and Lee (arXiv:2602.05904) improve coloring algorithms to O(n0.19539)O(n^{0.19539})O(n0.19539) colors; Fei, Minzer and Wang prove the 444-to-111 Games Conjecture with perfect completeness, giving constant-palette hardness and small independent-set density with an eight-colorable promise (ECCC TR26-179).
  • 2026 — An OpenAI preprint, Hardness of finding large independent sets in three-colorable graphs (OpenAI Math Release, September 24, 2026), claims unconditional NP-hardness of distinguishing three-colorable graphs from graphs with independence ratio below any fixed δ>0\delta>0δ>0. It has not been peer reviewed, and its proof is not formally verified.

Setting

A 3-CNF formula is a conjunction of clauses, each a disjunction of at most three literals (variables or their negations); it is satisfiable if some truth assignment makes every clause true. For a finite simple undirected graph GGG with n=∣V(G)∣≥1n=|V(G)|\ge1n=∣V(G)∣≥1 vertices, α(G)\alpha(G)α(G) is the size of a largest independent set, a set of pairwise nonadjacent vertices. A polynomial-time reduction is a deterministic algorithm that, on an explicit binary encoding of its input, runs in time polynomial in the encoding length and outputs an explicit encoding (here the full adjacency matrix) of the output graph.

Formalization targets

Goal: NP-hardness of the three-colorable vs. small-independence-ratio gap (Theorem 1.1)

For every fixed real 0<δ<1/30<\delta<1/30<δ<1/3 there is a deterministic polynomial-time algorithm RδR_\deltaRδ​ mapping each explicitly encoded 3-CNF formula ϕ\phiϕ to a nonempty finite simple graph GϕG_\phiGϕ​ with

ϕ satisfiable ⟹ Gϕ three-colorable,ϕ unsatisfiable ⟹ α(Gϕ)<δ ∣V(Gϕ)∣.\phi\ \text{satisfiable}\ \Longrightarrow\ G_\phi\ \text{three-colorable}, \qquad \phi\ \text{unsatisfiable}\ \Longrightarrow\ \alpha(G_\phi)<\delta\,|V(G_\phi)| .ϕ satisfiable ⟹ Gϕ​ three-colorable,ϕ unsatisfiable ⟹ α(Gϕ​)<δ∣V(Gϕ​)∣.

The running time and output size are polynomial in the binary length of ϕ\phiϕ, with constants depending only on δ\deltaδ. The goal statement is published on the platform with status Open.

Significance

The result itself. Theorem 1.1 shows that, even with the full promise of three-colorability, it is NP-hard to find an independent set of any fixed positive density. It implies hardness of distinguishing kkk-colorable graphs from non-ccc-colorable graphs for all fixed 3≤k≤c3\le k\le c3≤k≤c (Corollary 1.2), the long-standing conjecture recorded by Barto–Bulín–Krokhin–Opršal. Unlike the routes through ddd-to-111 games and arc graphs, it keeps the independent-set soundness together with three-colorable completeness, and unlike Dinur–Khot–Perkins–Safra it does not delete vertices. Earlier results with this combination were conditional (Dinur–Mossel–Regev; Braverman–Khot–Lifshitz–Minzer).

Formalizing it. The statement is fully explicit: Mathlib's TM2ComputableInPolyTime fixes the machine model, and the formula and graph encodings are spelled out bit by bit. A formal proof would require a machine-checked PCP theorem with perfect completeness (Label Cover), a dimension-independent junta approximation theorem, and the paper's alignment and list-decoding lemmas. Any of these is a substantial and reusable formalization project; complexity-theoretic reductions are barely represented in Mathlib today.

Difficulty

Gap reductions from Label Cover to coloring traditionally use long-code or dictatorship tests whose soundness analysis loses control either of three-colorability of every vertex (completeness) or of independent sets of arbitrarily small density (soundness). Known unconditional routes reach constant-palette coloring hardness only by composing reductions that lose the independent-set guarantee, and routes keeping it required unproven games conjectures with structural assumptions. The proof must decode a labeling from any independent set of density δ\deltaδ in a graph that is genuinely three-colorable on every vertex when ϕ\phiϕ is satisfiable, without bounds on projection-fiber sizes.

Formalization scope

  • Formulas are lists of clauses of width at most 333 over literals named by natural numbers; satisfiability uses an assignment N→Bool\mathbb N\to\mathrm{Bool}N→Bool. Formulas and graphs are encoded as explicit bit strings (self-delimiting binary numbers; the graph encoding lists the vertex count and the full row-major adjacency matrix).
  • Graphs are Graph records: n≥1n\ge1n≥1 vertices Fin n, a symmetric loopless Boolean adjacency. Independence number is the maximum card of an independent Finset; three-colorability is a proper map to Fin 3.
  • RealGraphReduction δ packages a function reduce : List Bool → Graph with a Turing.TM2ComputableInPolyTime witness (identity input encoding, graphBits output encoding, finite stack alphabets), an explicit polynomial output-size bound, completeness and strict soundness α<δn\alpha<\delta nα<δn on encoded formulas. The reduction must be total and polynomial-time on all bit strings, which is harmless since malformed inputs can be sent to a fixed graph.
  • The goal quantifies over every real δ∈(0,1/3)\delta\in(0,1/3)δ∈(0,1/3); the source notes that rational δ\deltaδ suffices by monotonicity.
  • Needed infrastructure: polynomial-time Turing-machine composition, Label Cover with perfect completeness, Fourier analysis on product spaces and junta approximation. Contributions on any of these are welcome.

Selected references

  • S. Khanna, N. Linial and S. Safra, On the hardness of approximating the chromatic number, Combinatorica (2000). https://doi.org/10.1007/s004930070013
  • V. Guruswami and S. Khanna, On the hardness of 4-coloring a 3-colorable graph, SIAM J. Discrete Math. (2004). https://doi.org/10.1137/S0895480100376794
  • I. Dinur, E. Mossel and O. Regev, Conditional hardness for approximate coloring, SIAM J. Comput. (2009). https://doi.org/10.1137/07068062X
  • I. Dinur, S. Khot, W. Perkins and M. Safra, Hardness of finding independent sets in almost 3-colorable graphs, FOCS 2010. https://doi.org/10.1109/FOCS.2010.84
  • V. Guruswami and S. Sandeep, d-to-1 hardness of coloring 3-colorable graphs with O(1) colors, ICALP 2020. https://doi.org/10.4230/LIPIcs.ICALP.2020.62
  • L. Barto, J. Bulín, A. Krokhin and J. Opršal, Algebraic approach to promise constraint satisfaction, J. ACM (2021). https://doi.org/10.1145/3457606
  • M. Braverman, S. Khot, N. Lifshitz and D. Minzer, An invariance principle for the multi-slice, with applications, FOCS 2021. https://doi.org/10.1109/FOCS52979.2021.00030
  • A. Krokhin, J. Opršal, M. Wrochna and S. Živný, Topology and adjunction in promise constraint satisfaction, SIAM J. Comput. (2023). https://doi.org/10.1137/20M1378223
  • N. Bansal, N. Huang and E. Lee, Improved SDP-based algorithm for coloring 3-colorable graphs, 2026. https://arxiv.org/abs/2602.05904v1
  • OpenAI, Hardness of finding large independent sets in three-colorable graphs, OpenAI Math Release preprint, September 24, 2026 (Theorem 1.1, p. 1). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Hardness-of-finding-large-independent-sets-in-three-colorable-graphs-September-24-2026/Hardness-of-finding-large-independent-sets-in-three-colorable-graphs-September-24-2026.pdf
2 thms1 active userReviewed
Complexity Theory·Captain: wurtle

Perfect completeness for 2-to-1 gamesResearch Paper

Motivation

A projection game (label cover) asks for labels on the two sides of a bipartite multigraph so that, on each edge, the right label is the image of the left label under a prescribed map. In a 2-to-1 game the left alphabet has 2q2q2q labels, the right alphabet qqq, and every edge map has exactly two preimages of each right label. Khot introduced the 2-to-1 and Unique Games conjectures in 2002 (STOC 2002) as hypotheses that yield tight hardness of approximation. The 2-to-1 Games Conjecture with perfect completeness asserts that, for every fixed δ>0\delta>0δ>0, it is NP-hard to tell satisfiable 2-to-1 games from games in which every labeling satisfies at most a δ\deltaδ fraction of edges. Perfect completeness matters for problems whose YES instances must be exactly satisfiable, such as graph coloring: hardness with completeness 1−ε1-\varepsilon1−ε does not imply it.

This mission asks for a formal proof of the conjecture, as stated in an OpenAI preprint dated September 23, 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

  • 1998 — The PCP theorem (Arora–Safra, JACM 1998; Arora–Lund–Motwani–Sudan–Szegedy, JACM 1998) and Raz's parallel repetition theorem (SICOMP 1998).
  • 2002 — Khot formulates the 2-to-1 and Unique Games conjectures.
  • 2011 — Rao's parallel repetition for projection games (SICOMP 2011).
  • 2014 — Austrin, O'Donnell, Tan and Wright prove perfect-completeness 2-to-1 label cover hardness with alphabets of sizes 666 and 333 and soundness 23/24+ε23/24+\varepsilon23/24+ε (ACM TOCT 2014).
  • 2017–2023 — The Grassmann/shortcode programme: Khot–Minzer–Safra (STOC 2017; ToC 2025), Dinur–Khot–Kindler–Minzer–Safra (STOC 2018), Barak–Kothari–Steurer (ITCS 2019), and the Grassmann expansion theorem of Khot–Minzer–Safra (Annals 2023), giving 2-to-1 hardness with completeness arbitrarily close to 111.
  • 2026 — Fei, Minzer and Wang prove perfect-completeness hardness for 4-to-1 games at arbitrarily small soundness (ECCC TR26-179).
  • September 2026 — The OpenAI preprint claims perfect completeness for 2-to-1 games (Theorem 1.1, p. 1).

Setting

For an integer q≥2q\ge2q≥2, a 2-to-1 instance has finite left and right vertex sets and a nonempty list of edge occurrences e=(ue,ve,πe)e=(u_e,v_e,\pi_e)e=(ue​,ve​,πe​), where πe:[2q]→[q]\pi_e:[2q]\to[q]πe​:[2q]→[q] is given by an explicit table with ∣πe−1(b)∣=2|\pi_e^{-1}(b)|=2∣πe−1​(b)∣=2 for every b∈[q]b\in[q]b∈[q]. Parallel occurrences, possibly with different maps, are allowed and counted with multiplicity; there are no weights. A labeling (a,b)(a,b)(a,b) assigns a(u)∈[2q]a(u)\in[2q]a(u)∈[2q] and b(v)∈[q]b(v)\in[q]b(v)∈[q] and satisfies eee when πe(a(ue))=b(ve)\pi_e(a(u_e))=b(v_e)πe​(a(ue​))=b(ve​). The value is the maximum, over labelings, of the fraction of satisfied occurrences. 3-SAT inputs are bit strings that decode to conjunctions of three-literal clauses (repetitions allowed); malformed inputs count as unsatisfiable.

Formalization targets

Goal: Theorem 1.1 (p. 1)

For every fixed rational δ∈(0,1)\delta\in(0,1)δ∈(0,1) there are an integer q=q(δ)≥2q=q(\delta)\ge2q=q(δ)≥2 and a deterministic polynomial-time reduction φ↦Gφ\varphi\mapsto G_\varphiφ↦Gφ​ from 3-SAT to 2-to-1 instances with alphabets [2q][2q][2q] and [q][q][q] such that

φ satisfiable ⇒ val⁡(Gφ)=1,φ unsatisfiable ⇒ val⁡(Gφ)≤δ.\varphi\ \text{satisfiable}\ \Rightarrow\ \operatorname{val}(G_\varphi)=1,\qquad \varphi\ \text{unsatisfiable}\ \Rightarrow\ \operatorname{val}(G_\varphi)\le\delta .φ satisfiable ⇒ val(Gφ​)=1,φ unsatisfiable ⇒ val(Gφ​)≤δ.

The running time is polynomial in the binary input length, with degree and constants depending only on δ\deltaδ.

Significance

The result itself. It removes the conditional hypothesis from perfect-completeness reductions. The preprint derives NP-hardness for maximum kkk-colorable subgraph at soundness 1−1k+Clog⁡kk21-\frac1k+C\frac{\log k}{k^2}1−k1​+Ck2logk​ through Guruswami and Sinop's reduction (Corollary 1.2, p. 3), and recovers satisfiable Not-Two and three-query PCP hardness (Corollary 1.3, p. 4), the latter already known unconditionally by Håstad (SICOMP 2014). A companion preprint uses it for hardness of finding large independent sets in three-colorable graphs.

Formalizing it. No PCP-style hardness result of this strength has been machine-checked. A complete development needs a perfect-completeness PCP gap (Theorem 2.1, p. 5), projection-game parallel repetition (Theorem 2.2, p. 6), the Grassmann expansion theorem (Theorem 2.3, p. 7) and inverse shortcode estimates (Theorem 2.6, p. 9), on top of a polynomial-time Turing-machine framework. These components are shared with the Unique Games formalization and would be reusable.

Difficulty

Near-perfect completeness does not give perfect completeness: the 1−ε1-\varepsilon1−ε reductions need not produce any finite satisfiable game. Amplifying a constant-soundness perfect-completeness game by parallel repetition destroys the 2-to-1 structure, since kkk-fold repetition gives 2k2^k2k preimages per right label. The preprint therefore builds a new tree-structured construction (Section 3) whose joint-partition questions keep two-element fibers and satisfy exact identities, and the soundness analysis must decode labelings through a rank test whose tested rows depend on the decoded form itself (Sections 5–6), before projection-game repetition gives the contradiction (Section 7).

Formalization scope

  • ProjectionTable q is a Vector (Fin q) (2*q) with exactly two preimages of every b : Fin q. Instance q lists edges as a List with a nonemptiness proof, so value = maxSatisfied / edges.length never divides by zero; maxSatisfied is a Finset.sup over all labelings.
  • BinaryGapReduction δ fixes alphabet ≥ 2, a construct map, a Turing.TM2ComputableInPolyTime witness from raw input bits (identity encoding) to the explicit bit encoding gameBits, finite tape alphabets, completeness value = 1 on encodings of satisfiable formulas, and soundness value ≤ δ on all other inputs.
  • The goal quantifies over rational δ\deltaδ with 0<δ<10<\delta<10<δ<1, as in the source; the alphabet and machine are fixed before the input.
  • No vacuous reading: the 3-SAT language is nontrivial and value = 1 requires a labeling satisfying every edge.

Selected references

  • OpenAI, Perfect completeness for 2-to-1 games, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Perfect-completeness-for-2-to-1-games-September-23-2026/paper.pdf
  • S. Khot, On the power of unique 2-prover 1-round games, STOC 2002. https://doi.org/10.1145/509907.510017
  • P. Austrin, R. O'Donnell, L.-Y. Tan, J. Wright, New NP-hardness results for 3-coloring and 2-to-1 label cover, ACM Trans. Comput. Theory 6 (2014).
  • S. Khot, D. Minzer, M. Safra, Pseudorandom sets in Grassmann graph have near-perfect expansion, Ann. of Math. 198 (2023). https://doi.org/10.4007/annals.2023.198.1.1
  • I. Dinur, S. Khot, G. Kindler, D. Minzer, M. Safra, Towards a proof of the 2-to-1 games conjecture?, STOC 2018. https://doi.org/10.1145/3188745.3188804
  • B. Barak, P. K. Kothari, D. Steurer, Small-set expansion in shortcode graph and the 2-to-2 conjecture, ITCS 2019. https://doi.org/10.4230/LIPIcs.ITCS.2019.9
  • Y. Fei, D. Minzer, S. Wang, On the hardness of 4-to-1 games with perfect completeness, ECCC TR26-179 (2026). https://eccc.weizmann.ac.il/report/2026/179/
  • A. Rao, Parallel repetition in projection games and a concentration bound, SIAM J. Comput. 40 (2011). https://doi.org/10.1137/080734042
  • R. Raz, A parallel repetition theorem, SIAM J. Comput. (1998). https://doi.org/10.1137/S0097539795280895
  • V. Guruswami, A. K. Sinop, Improved inapproximability results for maximum k-colorable subgraph, Theory Comput. 9 (2013). https://doi.org/10.4086/toc.2013.v009a011
  • M. Braverman, S. Khot, D. Minzer, On rich 2-to-1 games, ITCS 2021. https://doi.org/10.4230/LIPIcs.ITCS.2021.27
  • J. Håstad, On the NP-hardness of Max-Not-2, SIAM J. Comput. (2014). https://doi.org/10.1137/120882718
2 thms1 active userReviewed
Complexity TheoryGraph Theory·Captain: wurtle

Constant-factor hardness of directed feedback vertex setResearch Paper

Motivation

A directed feedback vertex set of a finite digraph GGG is a set FFF of vertices meeting every directed cycle, equivalently one whose removal leaves G−FG-FG−F acyclic; DFVS(G)\mathrm{DFVS}(G)DFVS(G) is the minimum size of such a set. The problem arises wherever cyclic dependencies must be broken — deadlock resolution, circuit testing, scheduling — and a topological order of G−FG-FG−F certifies a solution in polynomial time.

The best general approximation guarantee is O(log⁡nlog⁡log⁡n)O(\log n\log\log n)O(lognloglogn), originating in Seymour's work on fractional packing of directed circuits (Combinatorica 1995) and made constructive by Even, Naor, Schieber and Sudan (Algorithmica 1998). On the hardness side, NP-hardness was known only for particular constants (via Vertex Cover), while ruling out every constant factor required the Unique Games Conjecture. Whether a constant-factor approximation exists under P ≠ NP alone was open.

This mission asks for a machine-checked proof that approximating DFVS within any fixed factor A≥1A\ge1A≥1 is NP-hard on unweighted digraphs, as claimed in an OpenAI preprint dated September 23, 2026 (source, Theorem 1.1, p. 2). The preprint has not been peer reviewed and its claim has not been independently verified; on the platform the Lean goal is open.

Timeline

  • 1995–1998 — Seymour's fractional packing bound (Combinatorica 1995) and the algorithm of Even–Naor–Schieber–Sudan (Algorithmica 1998) give O(log⁡nlog⁡log⁡n)O(\log n\log\log n)O(lognloglogn) approximation.
  • 2005 — Dinur and Safra prove Vertex Cover NP-hard below 105−21≈1.3610\sqrt5-21\approx1.36105​−21≈1.36 (Annals 2005); doubling each undirected edge transfers this to DFVS.
  • 2008–2011 — Guruswami, Manokaran and Raghavendra show every constant factor is Unique-Games-hard via maximum acyclic subgraph and feedback arc set (FOCS 2008); expanded with Håstad and Charikar (SICOMP 2011).
  • 2013 — Svensson's direct Unique-Games-based construction for vertex deletion problems (ToC 2013).
  • 2016 — Guruswami and Lee give a simpler UG-based proof (ToC 2016).
  • 2018–2025 — The imperfect-completeness 2-to-1 Games Theorem is completed (Dinur–Khot–Kindler–Minzer–Safra, ToC 2025; Barak–Kothari–Steurer, ITCS 2019; Khot–Minzer–Safra, ECCC TR18-006), giving Vertex Cover, hence DFVS, hardness below 2\sqrt22​.
  • 2026 — Ghorbani and Mnich record the gap for general digraphs (ICALP 2026); the OpenAI preprint claims NP-hardness for every constant factor.

Setting

A digraph is finite and loopless, given by a vertex count nnn and a duplicate-free list of arcs (u,v)(u,v)(u,v) with u≠vu\ne vu=v; opposite arcs (u,v),(v,u)(u,v),(v,u)(u,v),(v,u) are allowed. A directed cycle is a cyclic sequence of r+1≥2r+1\ge2r+1≥2 distinct vertices with every successive arc present, including the closing one. FFF is a feedback set if it meets every directed cycle, and DFVS(G)\mathrm{DFVS}(G)DFVS(G) is the least size of one.

A language L⊆{0,1}∗L\subseteq\{0,1\}^*L⊆{0,1}∗ is in NP if there are a polynomial ppp and a deterministic polynomial-time verifier VVV (a multi-stack Turing machine with finite alphabets) with x∈Lx\in Lx∈L iff some certificate www with ∣w∣≤p(∣x∣)|w|\le p(|x|)∣w∣≤p(∣x∣) makes V(x,w)V(x,w)V(x,w) accept.

Formalization targets

Goal: constant-factor hardness (Theorem 1.1, p. 2)

For every real A≥1A\ge1A≥1 and every language L∈NPL\in\mathrm{NP}L∈NP there is a deterministic polynomial-time map x↦(Gx,kx)x\mapsto(G_x,k_x)x↦(Gx​,kx​) with kx≥1k_x\ge1kx​≥1 such that

x∈L ⇒ DFVS(Gx)≤kx,x∉L ⇒ DFVS(Gx)>A kx.x\in L\ \Rightarrow\ \mathrm{DFVS}(G_x)\le k_x,\qquad x\notin L\ \Rightarrow\ \mathrm{DFVS}(G_x)>A\,k_x .x∈L ⇒ DFVS(Gx​)≤kx​,x∈/L ⇒ DFVS(Gx​)>Akx​.

In Lean this is OAI.DirectedFeedback.main : DirectedFeedback.MainStatement. The preprint states Theorem 1.1 as a reduction from an NP-hard promise problem (the 2-to-1 game gap of Theorem 2.1, p. 5) to (G,k)(G,k)(G,k); composing with that promise problem's NP-hardness gives the Lean form, which quantifies over all NP languages directly. A deterministic polynomial-time AAA-approximation would then decide every NP language (proof of Theorem 1.1, p. 28).

The preprint also proves Corollary 8.3 (p. 28) for directed feedback arc set; it is not formalized on the platform.

Significance

The result itself. Theorem 1.1 shows that a constant-factor approximation for DFVS exists if and only if P = NP, removing the Unique Games hypothesis from the Guruswami–Manokaran–Raghavendra and Svensson results. It does not determine the true growth of the optimal factor between constants and O(log⁡nlog⁡log⁡n)O(\log n\log\log n)O(lognloglogn).

Formalizing it. No machine-checked inapproximability result of this kind is known to exist. A complete development needs the 2-to-1 Games Theorem with imperfect completeness (Theorem 2.1, p. 5, external input), Friedgut's junta theorem in a dyadic fixed-bias form (Theorem 2.2 and Corollary 2.3, p. 6), finite Ramsey theory and von Neumann's minimax theorem, and the preprint's rank-graph construction (Definition 3.1, p. 8; Proposition 3.7, p. 12), pivotality budget (Definition 6.1, p. 20) and weighted soundness (Proposition 7.4, p. 25). An NP and Karp-reduction library on TM2 machines (including Cook–Levin) would be reusable across the whole family of hardness missions.

Difficulty

The central difficulty (Introduction, pp. 3–4) is extracting a game labeling from an arbitrary acyclic remainder: a topological order exists, but it can interleave the many vertices representing one local test arbitrarily. Comparisons in that order must first be made consistent and then read as Boolean functions; Friedgut-type junta bounds must hold at many biases simultaneously for the same random function, so a separate influence bound at each scale does not suffice. The coordinate lists must also be fixed before the tested game edge is drawn, or they do not define a game labeling. Gadget reductions from Vertex Cover cannot go beyond the Vertex Cover constants.

Formalization scope

  • DirectedFeedback.Digraph: n : ℕ, a Nodup list of arcs in Fin n × Fin n, loopless; cycles are injective maps Fin (r+1) → Fin n with all successive arcs (including the wrap-around) present, so 2-cycles count and 1-cycles cannot occur. dfvs is the least cardinality of a feedback Finset.
  • NP is defined from scratch by NPVerifier: a polynomial certificate bound and a polynomial-time TM2 verifier with finite alphabets on the unary-framed pair encoding pairBits; no SAT-specific definition is used.
  • GapReduction A L is a polynomial-time TM2 machine (finite alphabets, raw-bit input) producing a GapInstance (digraph plus k > 0), encoded explicitly in unary-framed form by GapInstance.bits; completeness and soundness are as displayed, with soundness compared in ℝ.
  • The statement quantifies over all NP languages, so its proof must include (or reprove) NP-hardness of the source promise problem in this machine model — e.g. a Cook–Levin theorem for TM2 verifiers plus the 2-to-1 Games reduction.
  • k > 0 and the strict inequality rule out a trivial reading (e.g. empty graphs with k = 0).

Welcome contributions: Friedgut's junta theorem, finite Ramsey and minimax in usable forms, the 2-to-1 Games Theorem, and Cook–Levin for Turing.TM2ComputableInPolyTime.

Selected references

  • OpenAI, Constant-factor hardness of directed feedback vertex set, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Constant-factor-hardness-of-directed-feedback-vertex-set-September-23-2026/paper.pdf
  • P. D. Seymour, Packing directed circuits fractionally, Combinatorica 15 (1995). https://doi.org/10.1007/BF01200760
  • G. Even, J. Naor, B. Schieber, M. Sudan, Approximating minimum feedback sets and multicuts in directed graphs, Algorithmica (1998). https://doi.org/10.1007/PL00009191
  • I. Dinur, S. Safra, On the hardness of approximating minimum vertex cover, Ann. of Math. 162 (2005). https://doi.org/10.4007/annals.2005.162.439
  • V. Guruswami, R. Manokaran, P. Raghavendra, Beating the random ordering is hard: inapproximability of maximum acyclic subgraph, FOCS 2008. https://doi.org/10.1109/FOCS.2008.51
  • O. Svensson, Hardness of vertex deletion and project scheduling, Theory of Computing 9 (2013). https://doi.org/10.4086/toc.2013.v009a024
  • V. Guruswami, E. Lee, Simple proof of hardness of feedback vertex set, Theory of Computing 12 (2016). https://doi.org/10.4086/toc.2016.v012a006
  • I. Dinur, S. Khot, G. Kindler, D. Minzer, M. Safra, Towards a proof of the 2-to-1 games conjecture?, Theory of Computing 21 (2025). https://doi.org/10.4086/toc.2025.v021a011
  • S. Khot, D. Minzer, M. Safra, Pseudorandom sets in Grassmann graph have near-perfect expansion, ECCC TR18-006 (2018). https://eccc.weizmann.ac.il/report/2018/006/revision/2/download
  • E. Friedgut, Boolean functions with low average sensitivity depend on few coordinates, Combinatorica (1998). https://doi.org/10.1007/PL00009809
  • F. P. Ramsey, On a problem of formal logic, Proc. London Math. Soc. (1930). https://doi.org/10.1112/plms/s2-30.1.264
2 thms1 active userReviewed
Complexity TheoryGraph Theory·Captain: wurtle

Constant-factor hardness of Min-UnCutResearch Paper

Motivation

Min-UnCut asks, for a finite simple graph G=(V,E)G=(V,E)G=(V,E), for the minimum number of edges whose deletion makes GGG bipartite. Equivalently, for a bipartition V=S⊔(V∖S)V=S\sqcup(V\setminus S)V=S⊔(V∖S) let UncutG(S)\mathrm{Uncut}_G(S)UncutG​(S) count the edges with both endpoints on the same side, and set OPTuncut(G)=min⁡SUncutG(S)\mathrm{OPT}_{\mathrm{uncut}}(G)=\min_S\mathrm{Uncut}_G(S)OPTuncut​(G)=minS​UncutG​(S). The quantity measures distance from bipartiteness. Although it equals ∣E∣|E|∣E∣ minus the maximum cut, a good multiplicative approximation for Max-Cut gives nothing multiplicative for Min-UnCut on nearly bipartite graphs, which is exactly the regime of interest.

The best known algorithm (Agarwal, Charikar, Makarychev, Makarychev, STOC 2005) achieves an O(log⁡∣V∣)O(\sqrt{\log|V|})O(log∣V∣​) factor. Before the source preprint, NP-hardness was known only for specific constant factors below about 1.491.491.49, and hardness for every constant factor was known only under the Unique Games Conjecture.

This mission asks for a machine-checked proof that Min-UnCut has no polynomial-time approximation within any constant factor unless P = NP, as claimed in an OpenAI preprint dated September 23, 2026 (source, Theorem 1.1, p. 2). The preprint has not been peer reviewed and its claim has not been independently verified; on the platform the Lean goal is open.

Timeline

  • 2001 — Håstad's near-satisfiable gap for three-variable parity equations (JACM 2001); combined with parallel repetition (Raz, SICOMP 1998; Holenstein, ToC 2009) it is the preprint's starting point.
  • 2005 — Agarwal, Charikar, Makarychev and Makarychev give O(log⁡n)O(\sqrt{\log n})O(logn​) approximations for Min-UnCut and Min 2CNF Deletion (STOC 2005).
  • 2007 — Under the Unique Games Conjecture, the near-satisfiable Max-Cut gap of Khot, Kindler, Mossel and O'Donnell implies arbitrarily large constant gaps for weighted Min-UnCut (SICOMP 2007).
  • 2017 — Håstad, Huang, Manokaran, O'Donnell and Wright: NP-hard below 11/811/811/8 for the weighted deletion formulation (ToC 2017).
  • 2018 — Wiman: factors below 1.456851.456851.45685 (KTH master's thesis, 2018).
  • 2024 — Martinsson: factors below 73139148/49096883≈1.4896973139148/49096883\approx1.4896973139148/49096883≈1.48969 for violated two-variable parity equations, transferring to weighted Min-UnCut (APPROX/RANDOM 2024).
  • 2026 — The OpenAI preprint claims NP-hardness for every constant factor on simple unweighted graphs.

Setting

A graph is finite, simple, undirected and unweighted, given by a vertex count NNN and a symmetric loopless Boolean adjacency table. A cut is a map c:{0,…,N−1}→{0,1}c:\{0,\dots,N-1\}\to\{0,1\}c:{0,…,N−1}→{0,1}; the parts may be empty or unbalanced. Uncut(c)\mathrm{Uncut}(c)Uncut(c) counts unordered edges {u,v}\{u,v\}{u,v} with c(u)=c(v)c(u)=c(v)c(u)=c(v), and OPTuncut(G)\mathrm{OPT}_{\mathrm{uncut}}(G)OPTuncut​(G) is its minimum over all cuts.

3SAT formulas are conjunctions of three-slot clauses over natural-number variable names (repeated variables allowed), given to the reduction by a canonical prefix-free binary encoding.

Formalization targets

Goal: constant-factor hardness (Theorem 1.1, p. 2)

For every fixed integer K≥2K\ge2K≥2 there is a deterministic polynomial-time reduction φ↦(G,k)\varphi\mapsto(G,k)φ↦(G,k) to an explicit simple graph GGG and a binary integer k≥1k\ge1k≥1 with

φ satisfiable⇒OPTuncut(G)≤k,φ unsatisfiable⇒OPTuncut(G)>K k,\varphi\ \text{satisfiable}\Rightarrow \mathrm{OPT}_{\mathrm{uncut}}(G)\le k,\qquad \varphi\ \text{unsatisfiable}\Rightarrow \mathrm{OPT}_{\mathrm{uncut}}(G)>K\,k,φ satisfiable⇒OPTuncut​(G)≤k,φ unsatisfiable⇒OPTuncut​(G)>Kk,

with graph size and running time polynomial in the bit length of φ\varphiφ for each fixed KKK. Choosing K≥CK\ge CK≥C shows that a deterministic polynomial-time CCC-approximation for any real C>1C>1C>1 would imply P = NP. The positive threshold k≥1k\ge1k≥1 is what separates this from bipartiteness testing. In Lean this is OAI.MinUncut.main : Nonempty MinUncut.Reduction.

The preprint also derives Corollary 1.2 (p. 3): minimum 2CNF clause deletion is NP-hard to approximate within every fixed factor C>1C>1C>1. It is not formalized on the platform.

Significance

The result itself. Theorem 1.1 settles the constant-factor approximability of Min-UnCut in the hardness direction without the Unique Games Conjecture, and by an exact cost-preserving transformation the same holds for Min 2CNF Deletion. Together with the O(log⁡n)O(\sqrt{\log n})O(logn​) algorithm, it places these problems strictly between constant-factor approximable and inapproximable within logarithmic factors (the remaining gap concerns super-constant factors).

Formalizing it. No machine-checked proof of super-constant-type inapproximability for any natural graph problem is known to exist. A complete development needs Håstad's parity gap (Theorem 2.2, p. 6) and uniform parallel repetition (Theorem 2.3, p. 7), used as external inputs, plus the preprint's inner decoding theorem (Theorem 3.2, p. 9), the composition gap for bit comparisons (Proposition 6.1, p. 27) and the finite constructions (Lemmas 7.1 and 7.2, pp. 31 and 34). PCP and parallel-repetition infrastructure would be reusable across many hardness missions.

Difficulty

A natural route is to compose a long-code test with a projection game, as in Håstad's work. The obstacle is extracting a label from an arbitrary folded bit proof with a probability and coefficient size independent of the label alphabet: the alphabet grows as the outer soundness is driven down, and every inner constant must be fixed before it (Introduction, pp. 3–4). In addition, the decoder's choices depend on shared background affine functions, so the outer game's soundness must survive such shared hints. Gadget-based approaches yield only isolated constants (up to about 1.491.491.49); unbounded constants previously required the Unique Games hypothesis.

Formalization scope

  • MinUncut.Output packages the graph (vertex count, symmetric loopless Bool adjacency) with a threshold threshold ≥ 1; uncut counts pairs u < v that are adjacent and on the same side, and opt is the minimum over all Fin N → Bool cuts (nonempty, so always defined).
  • MinUncut.Reduction is a single TM2 machine taking the pair (K,φ)(K,\varphi)(K,φ) (encoded by inputBits) with finite tape alphabets, a polynomial time K bounding the running time in the formula's bit length, and a polynomial size K bounding vertices plus output length; yes/no hold for every K ≥ 2.
  • The reduction acts on decoded formulas, so malformed inputs do not arise; runtime is measured in formulaBits φ, the canonical encoding including variable names.
  • The hypotheses are satisfiable and the conclusion demands an actual strict multiplicative gap with a positive integer threshold, so outputting a trivially bipartite graph cannot satisfy it.

Welcome contributions: Fourier and Gaussian analysis on F2\mathbb F_2F2​-affine spaces, the long code with folding, parallel repetition for projection games, and Turing-machine complexity infrastructure.

Selected references

  • OpenAI, Constant-factor hardness of Min-UnCut, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Constant-factor-hardness-of-Min-UnCut-September-23-2026/paper.pdf
  • A. Agarwal, M. Charikar, K. Makarychev, Y. Makarychev, O(√log n) approximation algorithms for Min UnCut, Min 2CNF Deletion, and directed cut problems, STOC 2005. https://doi.org/10.1145/1060590.1060675
  • S. Khot, G. Kindler, E. Mossel, R. O'Donnell, Optimal inapproximability results for MAX-CUT and other 2-variable CSPs?, SIAM J. Comput. 37 (2007). https://doi.org/10.1137/S0097539705447372
  • J. Håstad, S. Huang, R. Manokaran, R. O'Donnell, J. Wright, Improved NP-inapproximability for 2-variable linear equations, Theory of Computing (2017). https://doi.org/10.4086/toc.2017.v013a019
  • M. Wiman, Improved inapproximability of Max-Cut through Min-Cut, master's thesis, KTH, 2018. https://www.diva-portal.org/smash/get/diva2:1229467/FULLTEXT01.pdf
  • B. Martinsson, On the NP-hardness approximation curve for Max-2Lin(2), APPROX/RANDOM 2024. https://doi.org/10.4230/LIPIcs.APPROX/RANDOM.2024.11
  • J. Håstad, Some optimal inapproximability results, J. ACM 48 (2001). https://doi.org/10.1145/502090.502098
  • R. Raz, A parallel repetition theorem, SIAM J. Comput. (1998). https://doi.org/10.1137/S0097539795280895
  • T. Holenstein, Parallel repetition: simplifications and the no-signaling case, Theory of Computing (2009). https://doi.org/10.4086/toc.2009.v005a008
  • M. Bellare, O. Goldreich, M. Sudan, Free bits, PCPs, and nonapproximability — towards tight results, SIAM J. Comput. (1998). https://doi.org/10.1137/S0097539796302531
2 thms1 active userReviewed
Complexity TheoryGraph Theory·Captain: wurtle

The Factor-Two Hardness Threshold for Vertex CoverResearch Paper

Motivation

A vertex cover of a graph is a set of vertices meeting every edge; τ(G)\tau(G)τ(G) denotes the minimum size of one. Computing τ(G)\tau(G)τ(G) exactly is NP-hard, but a two-line algorithm — take both endpoints of a maximal matching — always returns a cover of size at most 2τ(G)2\tau(G)2τ(G). Whether any polynomial-time algorithm achieves a fixed factor 2−ε2-\varepsilon2−ε is one of the best-known questions in approximation algorithms. Unconditional NP-hardness was known below 105−21≈1.3610\sqrt5-21\approx1.36105​−21≈1.36 and, after the 2-to-2 Games Theorem, below 2\sqrt22​; factor-two hardness was known only under Khot's Unique Games Conjecture.

This mission asks for a machine-checked proof that every fixed approximation factor below two is NP-hard, as claimed in an OpenAI preprint dated September 23, 2026 (source, Corollary 1.2, p. 1). The preprint's direct reduction starts from ordinary perfect-completeness Label Cover. It has not been peer reviewed and its claim has not been independently verified; on the platform the Lean goal is open.

Timeline

  • 1972 — Karp proves Node Cover NP-complete (Karp 1972); the maximal-matching algorithm gives factor 222.
  • 1992–1998 — Feige–Goldwasser–Lovász–Safra–Szegedy (JACM 1996) link proof checking to approximation; the PCP theorem (Arora–Safra, JACM 1998; Arora–Lund–Motwani–Sudan–Szegedy, JACM 1998) gives a constant gap; Raz's parallel repetition (SICOMP 1998) yields Label Cover with arbitrarily small soundness.
  • 2001 — Håstad: NP-hard below factor 7/67/67/6 (JACM 2001).
  • 2005 — Dinur–Safra: NP-hard below 105−21≈1.3606810\sqrt5-21\approx1.36068105​−21≈1.36068 (Annals 2005).
  • 2002–2008 — Khot's Unique Games Conjecture (STOC 2002); Khot–Regev prove factor-(2−ε)(2-\varepsilon)(2−ε) hardness assuming it (JCSS 2008).
  • 2009 — Karakostas's algorithm achieves 2−Θ(1/log⁡n)2-\Theta(1/\sqrt{\log n})2−Θ(1/logn​) (TALG 2009), a saving that vanishes with nnn.
  • 2017–2023 — The 2-to-2 Games programme (Khot–Minzer–Safra, ToC 2025; Dinur–Khot–Kindler–Minzer–Safra, ToC 2025; Barak–Kothari–Steurer, ITCS 2019; Khot–Minzer–Safra expansion, Annals 2023) gives NP-hardness below 2\sqrt22​.
  • 2026 — The OpenAI preprint claims NP-hardness for every fixed factor below 222.

Setting

A graph is finite, simple, undirected and unweighted, given explicitly by its vertex count nnn and a duplicate-free list of edges {u,v}\{u,v\}{u,v} with u<vu<vu<v. A set SSS of vertices is a cover if every edge has an endpoint in SSS, and τ(G)\tau(G)τ(G) is the least size of a cover. An α\alphaα-approximation algorithm is a deterministic polynomial-time machine that, on every graph, outputs a duplicate-free list of vertices that covers all edges and has length at most α τ(G)\alpha\,\tau(G)ατ(G).

3SAT inputs are bit strings decoding canonically to a conjunction of three-slot clauses; malformed strings are NO instances. A 3SAT decider is a deterministic polynomial-time machine answering membership correctly on every input.

Formalization targets

Milestone: the explicit cover gap (Theorem 1.1, p. 1)

For every integer m≥4m\ge4m≥4 there is a deterministic polynomial-time reduction φ↦Gφ\varphi\mapsto G_\varphiφ↦Gφ​ from 3SAT to explicit simple graphs with

φ satisfiable⇒τ(Gφ)<(12+1m)∣V(Gφ)∣,φ unsatisfiable⇒τ(Gφ)>(1−1m)∣V(Gφ)∣.\varphi\text{ satisfiable}\Rightarrow \tau(G_\varphi)<\Bigl(\tfrac12+\tfrac1m\Bigr)|V(G_\varphi)|,\qquad \varphi\text{ unsatisfiable}\Rightarrow \tau(G_\varphi)>\Bigl(1-\tfrac1m\Bigr)|V(G_\varphi)|.φ satisfiable⇒τ(Gφ​)<(21​+m1​)∣V(Gφ​)∣,φ unsatisfiable⇒τ(Gφ​)>(1−m1​)∣V(Gφ​)∣.

The ratio of the thresholds is 2−6/(m+2)→22-6/(m+2)\to22−6/(m+2)→2. In Lean: OAI.VertexCover.explicit_cover_gap.

Goal: factor-two threshold (Corollary 1.2, p. 1)

For every real α\alphaα with 1≤α<21\le\alpha<21≤α<2, an α\alphaα-approximation algorithm for Vertex Cover yields a polynomial-time 3SAT decider:

∀ α∈[1,2):Approximation(α)  ⟹  ThreeSATDecision.\forall\,\alpha\in[1,2):\quad \mathrm{Approximation}(\alpha)\;\Longrightarrow\;\mathrm{ThreeSATDecision}.∀α∈[1,2):Approximation(α)⟹ThreeSATDecision.

In Lean: OAI.VertexCover.every_fixed_factor_below_two. The preprint derives it from Theorem 1.1 in its final section (pp. 23–24) by choosing mmm from α\alphaα and comparing the returned cover with an exact rational threshold. Together with the matching algorithm, it identifies 222 as the approximation threshold unless P = NP.

Significance

The result itself. The statement removes the Unique Games hypothesis from the Khot–Regev theorem and closes the gap between the elementary factor-222 algorithm and hardness. The density gap of Theorem 1.1 is also a statement about independent sets: YES graphs have independent sets of density nearly 1/21/21/2, NO graphs have none of density above 1/m1/m1/m.

Formalizing it. No machine-checked proof of any super-constant Vertex Cover hardness bound is known to exist. A complete development requires the perfect-completeness Label Cover theorem (Theorem 2.1, p. 4; PCP theorem plus parallel repetition, used as an external input), the explicit graph construction and its polynomial-time implementation (Proposition 3.3, p. 7), and the preprint's analytic soundness argument. The Label Cover input and a library for polynomial-time gap reductions on TM2 machines are reusable for many other hardness missions.

Difficulty

Soundness requires extracting short lists of candidate labels from an arbitrary large independent set, where each list may depend only on its own position's query. If a list also sees the opposite endpoint of a constraint or its occurrence index, the pair-decoding argument against Label Cover soundness breaks down (Introduction, pp. 3–4). Classical approaches either need the stronger Unique Games source promise (Khot–Regev) or only reach densities near 1−1/21-1/\sqrt21−1/2​ (the 2-to-2 route), which gives factor 2\sqrt22​ rather than 222. Averaging to impose the information restriction can erase all variance; the analytic core is retaining it without any dependence on the label alphabet sizes.

Formalization scope

  • Graphs are VertexCover.ExplicitGraph: n : ℕ and a Nodup list of pairs in Fin n × Fin n with e.1 < e.2; coverNumber is the least cardinality of a covering Finset (always defined since univ covers).
  • Encodings: inputs are raw bits; graphs are encoded by ExplicitGraph.bits (vertex count, all vertex names, edge count, edge endpoints, each with prefix-free binary framing); an approximation's output list is encoded by natListBits.
  • Approximation α bundles a polynomial-time TM2 program on graph encodings with finite tape alphabets, whose output is duplicate-free, in range, covers every edge and has length at most α * coverNumber. ThreeSATDecision is a polynomial-time TM2 decider correct on every input.
  • The goal is an implication from an algorithm to a decider: no assumption P ≠ NP is built in, and the hypothesis 1≤α<21\le\alpha<21≤α<2 is satisfied by real algorithms only if 3SAT is easy, which is exactly the content. It cannot be discharged trivially because Approximation α must be correct on every graph.
  • The milestone's thresholds are compared in ℝ with strict inequalities as in the source.

Welcome contributions: Label Cover hardness (PCP theorem + parallel repetition), Efron–Stein decompositions and variance inequalities on product spaces, and general complexity-theoretic infrastructure for Turing.TM2ComputableInPolyTime.

Selected references

  • OpenAI, The Factor-Two Hardness Threshold for Vertex Cover, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-Factor-Two-Hardness-Threshold-for-Vertex-Cover-September-23-2026/paper.pdf
  • R. M. Karp, Reducibility among combinatorial problems, Complexity of Computer Computations, 1972. https://doi.org/10.1007/978-1-4684-2001-2_9
  • J. Håstad, Some optimal inapproximability results, J. ACM 48 (2001). https://doi.org/10.1145/502090.502098
  • I. Dinur, S. Safra, On the hardness of approximating minimum vertex cover, Ann. of Math. 162 (2005). https://doi.org/10.4007/annals.2005.162.439
  • S. Khot, O. Regev, Vertex cover might be hard to approximate to within 2−ε, J. Comput. Syst. Sci. 74 (2008). https://doi.org/10.1016/j.jcss.2007.06.019
  • G. Karakostas, A better approximation ratio for the Vertex Cover problem, ACM Trans. Algorithms (2009). https://doi.org/10.1145/1597036.1597045
  • R. Raz, A parallel repetition theorem, SIAM J. Comput. (1998). https://doi.org/10.1137/S0097539795280895
  • S. Khot, D. Minzer, M. Safra, Pseudorandom sets in Grassmann graph have near-perfect expansion, Ann. of Math. 198 (2023). https://doi.org/10.4007/annals.2023.198.1.1
  • I. Dinur, S. Khot, G. Kindler, D. Minzer, M. Safra, Towards a proof of the 2-to-1 games conjecture?, Theory of Computing 21 (2025). https://doi.org/10.4086/toc.2025.v021a011
  • S. Khot, On the power of unique 2-prover 1-round games, STOC 2002. https://doi.org/10.1145/509907.510017
3 thms1 active userReviewed
Complexity TheoryOperations Research·Captain: wurtle

A Direct Proof of Optimal Max-Cut HardnessResearch Paper

Motivation

Max-Cut asks for a partition of the vertices of a graph that maximizes the number of edges crossing it. Its decision version is one of Karp's original NP-complete problems, so the natural question is how well it can be approximated in polynomial time. Goemans and Williamson's semidefinite-programming algorithm with random-hyperplane rounding achieves ratio approaching

αGW=min⁡−1≤ρ<12arccos⁡ρπ(1−ρ)=0.878567…\alpha_{\mathrm{GW}}=\min_{-1\le\rho<1}\frac{2\arccos\rho}{\pi(1-\rho)}=0.878567\ldotsαGW​=−1≤ρ<1min​π(1−ρ)2arccosρ​=0.878567…

(GW 1995). Whether any polynomial-time algorithm can beat this constant is a central question of approximation algorithms: unconditional NP-hardness was known only above 16/1716/1716/17, and optimality of αGW\alpha_{\mathrm{GW}}αGW​ was known only under Khot's Unique Games Conjecture.

This mission asks for a machine-checked proof that approximating Max-Cut within any fixed ratio α∈(αGW,1]\alpha\in(\alpha_{\mathrm{GW}},1]α∈(αGW​,1] is NP-hard, even on simple unweighted graphs, as claimed in an OpenAI preprint dated September 23, 2026 (source, Theorem 1.1, p. 1). The preprint gives a direct reduction from the 2-to-1 Games Theorem, not via the Unique Games Conjecture. It has not been peer reviewed and its claim has not been independently verified; on the platform the Lean goal is open.

Timeline

  • 1972 — Karp lists Max-Cut among the original NP-complete problems (Karp 1972).
  • 1995 — Goemans and Williamson give the SDP algorithm with ratio αGW\alpha_{\mathrm{GW}}αGW​ (JACM 1995).
  • 2000–2001 — Håstad's optimal inapproximability results (JACM 2001) combined with the gadgets of Trevisan, Sorkin, Sudan and Williamson (SICOMP 2000) rule out ratios above 16/1716/1716/17.
  • 2002 — Feige and Schechtman show the SDP integrality gap approaches αGW\alpha_{\mathrm{GW}}αGW​ (RSA 2002); Khot formulates the Unique Games Conjecture (STOC 2002).
  • 2007–2010 — Khot, Kindler, Mossel and O'Donnell prove that αGW\alpha_{\mathrm{GW}}αGW​ is optimal assuming the UGC, conditional on Majority Is Stablest (SICOMP 2007), which Mossel, O'Donnell and Oleszkiewicz prove (Annals 2010) via Borell's Gaussian noise-stability bound (1985). O'Donnell–Wu (STOC 2008) and Raghavendra (STOC 2008) extend the UG-based picture.
  • 2017–2025 — The Grassmann-graph programme establishes the 2-to-1 Games Theorem with imperfect completeness: Khot–Minzer–Safra (ToC 2025), Dinur–Khot–Kindler–Minzer–Safra (ToC 2025), Barak–Kothari–Steurer (ITCS 2019), Khot–Minzer–Safra expansion (Annals 2023).
  • 2026 — The OpenAI preprint claims unconditional NP-hardness at every ratio above αGW\alpha_{\mathrm{GW}}αGW​.

Setting

A simple unweighted graph on NNN vertices is given by a symmetric, loopless Boolean adjacency table. A cut is a map ccc from the vertices to {0,1}\{0,1\}{0,1}; its size is the number of unordered edges {u,v}\{u,v\}{u,v} with c(u)≠c(v)c(u)\ne c(v)c(u)=c(v), and MaxCut(G)\mathrm{MaxCut}(G)MaxCut(G) is the largest cut size. An algorithm is an α\alphaα-approximation if it always returns a cut of size at least α⋅MaxCut(G)\alpha\cdot\mathrm{MaxCut}(G)α⋅MaxCut(G).

3SAT inputs are bit strings decoding canonically to a conjunction of three-slot clauses (repeated variables and literals allowed); malformed strings are NO instances. A gap reduction for α\alphaα consists of rationals Y>0Y>0Y>0 and N≥0N\ge0N≥0 with N<αYN<\alpha YN<αY, and a deterministic polynomial-time map φ↦(Gφ,Qφ)\varphi\mapsto(G_\varphi,Q_\varphi)φ↦(Gφ​,Qφ​) outputting a graph together with a positive integer scale QφQ_\varphiQφ​ such that

φ∈3SAT⇒MaxCut(Gφ)Qφ≥Y,φ∉3SAT⇒MaxCut(Gφ)Qφ≤N.\varphi\in\mathrm{3SAT}\Rightarrow \frac{\mathrm{MaxCut}(G_\varphi)}{Q_\varphi}\ge Y,\qquad \varphi\notin\mathrm{3SAT}\Rightarrow \frac{\mathrm{MaxCut}(G_\varphi)}{Q_\varphi}\le N.φ∈3SAT⇒Qφ​MaxCut(Gφ​)​≥Y,φ∈/3SAT⇒Qφ​MaxCut(Gφ​)​≤N.

Given such a reduction, an α\alphaα-approximation algorithm would decide 3SAT by comparing its output with NQφN Q_\varphiNQφ​.

Formalization targets

Goal: optimal Max-Cut hardness (Theorem 1.1, p. 1)

For every real α\alphaα with αGW<α≤1\alpha_{\mathrm{GW}}<\alpha\le1αGW​<α≤1, a gap reduction for α\alphaα from 3SAT to simple unweighted Max-Cut exists:

∀ α∈(αGW,1]:GapReduction(α)≠∅.\forall\,\alpha\in(\alpha_{\mathrm{GW}},1]:\quad \mathrm{GapReduction}(\alpha)\ne\varnothing.∀α∈(αGW​,1]:GapReduction(α)=∅.

In Lean this is OAI.OptimalMaxCut.main : OptimalMaxCut.MainStatement. The preprint proves it from a sharper weighted gap (Theorem 1.2, p. 2): for rational t∈(0,1)t\in(0,1)t∈(0,1), B(t)=2πarcsin⁡tB(t)=\tfrac{2}{\pi}\arcsin tB(t)=π2​arcsint and 0<ε<(t−B(t))/40<\varepsilon<(t-B(t))/40<ε<(t−B(t))/4, it is NP-hard to distinguish Val≥1+t2−ε\mathrm{Val}\ge\frac{1+t}{2}-\varepsilonVal≥21+t​−ε from Val≤1+B(t)2+ε\mathrm{Val}\le\frac{1+B(t)}{2}+\varepsilonVal≤21+B(t)​+ε on graphs with a rational edge distribution; Appendix B (p. 35) converts weighted graphs to simple unweighted ones.

Significance

The result itself. Theorem 1.1 shows that the Goemans–Williamson algorithm is optimal among polynomial-time algorithms unless P = NP, removing the Unique Games hypothesis from the result of Khot, Kindler, Mossel and O'Donnell. It is the most prominent instance of an SDP-tight threshold, and the gap form (Theorem 1.2) gives a whole curve of completeness/soundness pairs.

Formalizing it. No machine-checked proof of any optimal inapproximability result for Max-Cut exists. A full development requires the 2-to-1 Games Theorem in the affine form stated as Proposition 2.2 (p. 6) and derived in Appendix A, the Majority Is Stablest theorem (Lemma 2.1, p. 5), and the paper's new ingredients: tree-code Fourier analysis (Section 5), a vector-valued affine hashing lemma (Lemma 6.1, p. 14), dimension-independent extraction (Proposition 7.1, p. 18) and decoding (Proposition 8.4, p. 24). Several of these — Majority Is Stablest, Boolean Fourier analysis, polynomial-time gap reductions on Turing machines — are reusable libraries in their own right.

Difficulty

The classical long-code test for Max-Cut, analysed via Majority Is Stablest, identifies an influential coordinate of an arbitrary cut. Under the Unique Games Conjecture that coordinate is directly a label of the source game. Starting instead from 2-to-1 games, the influential coordinate is indexed by a gate input of the test, not by a source-game label, and the two endpoints of a source edge cannot in general reconstruct a common context without knowing their projection. The decoding loss must also be independent of the source alphabet, because the alphabet dimension is chosen after the target soundness; a dimension-dependent loss would make the parameter choice circular (Introduction, pp. 3–4).

Formalization scope

  • Graphs are OptimalMaxCut.Graph: a vertex count with a symmetric loopless Bool adjacency function; maxCut maximizes cutSize over all Fin n → Bool, counting each unordered crossing edge once.
  • The output is a ScaledGraph with an explicit positive integer scale; the YES/NO bounds compare maxCut / scale in ℝ, and 0 < scale rules out division by zero.
  • alphaGW is defined as the sInf of 2arccos⁡ρ/(π(1−ρ))2\arccos\rho/(\pi(1-\rho))2arccosρ/(π(1−ρ)) over ρ∈[−1,1)\rho\in[-1,1)ρ∈[−1,1); this set is nonempty and bounded below, so the infimum is the true constant.
  • The reduction is a Turing.TM2ComputableInPolyTime machine with finite tape alphabets, from raw input bits (identity encoding) to ScaledGraph.bits (delimited vertex count and scale, then the full adjacency table); the rational bounds and the machine are fixed before the input.
  • The hypothesis αGW<α≤1\alpha_{\mathrm{GW}}<\alpha\le1αGW​<α≤1 is satisfiable, and the conclusion requires an actual strict gap noBound < α * yesBound with 0 < yesBound, so there is no trivial reading. No assumption P ≠ NP is built in.

Welcome contributions: Boolean and Gaussian Fourier analysis (noise stability, influences, Borell's theorem), Majority Is Stablest, the 2-to-1 Games Theorem, and a reusable library for polynomial-time gap reductions on TM2 machines.

Selected references

  • OpenAI, A Direct Proof of Optimal Max-Cut Hardness, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/A-Direct-Proof-of-Optimal-Max-Cut-Hardness-September-23-2026/paper.pdf
  • M. X. Goemans, D. P. Williamson, Improved approximation algorithms for maximum cut and satisfiability problems using semidefinite programming, J. ACM 42 (1995). https://doi.org/10.1145/227683.227684
  • J. Håstad, Some optimal inapproximability results, J. ACM 48 (2001). https://doi.org/10.1145/502090.502098
  • L. Trevisan, G. Sorkin, M. Sudan, D. Williamson, Gadgets, approximation, and linear programming, SIAM J. Comput. 29 (2000). https://doi.org/10.1137/S0097539797328847
  • S. Khot, G. Kindler, E. Mossel, R. O'Donnell, Optimal inapproximability results for MAX-CUT and other 2-variable CSPs?, SIAM J. Comput. 37 (2007). https://doi.org/10.1137/S0097539705447372
  • E. Mossel, R. O'Donnell, K. Oleszkiewicz, Noise stability of functions with low influences: invariance and optimality, Ann. of Math. 171 (2010). https://doi.org/10.4007/annals.2010.171.295
  • U. Feige, G. Schechtman, On the optimality of the random hyperplane rounding technique for MAX CUT, Random Structures Algorithms 20 (2002). https://doi.org/10.1002/rsa.10036
  • I. Dinur, S. Khot, G. Kindler, D. Minzer, M. Safra, Towards a proof of the 2-to-1 games conjecture?, Theory of Computing 21 (2025). https://doi.org/10.4086/toc.2025.v021a011
  • S. Khot, D. Minzer, M. Safra, Pseudorandom sets in Grassmann graph have near-perfect expansion, Ann. of Math. 198 (2023). https://doi.org/10.4007/annals.2023.198.1.1
  • S. Khot, On the power of unique 2-prover 1-round games, STOC 2002. https://doi.org/10.1145/509907.510017
2 thms1 active userReviewed
Complexity Theory·Captain: wurtle

The Unique Games TheoremResearch Paper

Motivation

A Unique Game is a system of constraints between pairs of variables, each constraint saying that the label of one endpoint is a fixed permutation of the label of the other. Deciding whether all constraints can be satisfied is easy (propagate labels along a spanning tree), so the interesting question is approximate: given an instance in which almost all constraints can be satisfied, can a polynomial-time algorithm find a labeling that satisfies even a small fraction?

Khot's Unique Games Conjecture (UGC, 2002) asserts that it cannot, unless P = NP. The conjecture matters because many optimal approximation thresholds have been proved conditionally on it: the Goemans–Williamson ratio for Max-Cut (KKMO 2007), the factor 222 for Vertex Cover (Khot–Regev 2008), the basic-SDP threshold for every constraint satisfaction problem (Raghavendra 2008), and hardness for ordering, multicut/sparsest-cut and correlation-clustering problems. An unconditional proof would turn all of these into NP-hardness results.

This mission asks for a machine-checked proof of the Unique Games Theorem as stated in an OpenAI preprint dated September 23, 2026 (source), which claims a positive resolution of the conjecture. The preprint has not been peer reviewed, and its claim has not been independently verified; on the platform the Lean goal below is open.

Timeline

  • 2001 — Håstad proves optimal inapproximability for linear equations mod 2 with near-perfect completeness (JACM 2001); this gap is the hardness input used by the preprint.
  • 2002 — Khot introduces Unique Games and formulates the conjecture (STOC 2002).
  • 2005–2008 — Conditional optimal hardness for Max-Cut (Khot, Kindler, Mossel, O'Donnell, SICOMP 2007), Vertex Cover (Khot–Regev, JCSS 2008) and every CSP (Raghavendra, STOC 2008). Khot–Vishnoi build SDP integrality gaps (JACM 2015, FOCS 2005).
  • 2005–2010 — Algorithms that delimit the conjecture's quantifiers: Trevisan (ToC 2008), Charikar–Makarychev–Makarychev with value 1−O(ϵlog⁡q)1-O(\sqrt{\epsilon\log q})1−O(ϵlogq​) (STOC 2006), and the subexponential algorithm of Arora–Barak–Steurer (FOCS 2010).
  • 2017–2018 — The 2-to-2 Games programme: Khot–Minzer–Safra (STOC 2017; ToC 2025), Dinur–Khot–Kindler–Minzer–Safra (STOC 2018; ToC 2025), Barak–Kothari–Steurer's shortcode formulation (ITCS 2019), and the Grassmann expansion theorem of Khot–Minzer–Safra (Annals 2023). Together these give NP-hardness of Unique Games with completeness near 1/21/21/2 and arbitrarily small soundness.
  • 2026 — The OpenAI preprint claims completeness arbitrarily close to 111, i.e. the full conjecture (Theorem 1.1, p. 2).

Setting

Fix a finite alphabet KKK. An instance GGG has a finite vertex set and a nonempty list of oriented constraints e=(ue,ve,πe)e=(u_e,v_e,\pi_e)e=(ue​,ve​,πe​) with πe\pi_eπe​ a permutation of KKK; repeated list entries count with multiplicity, and no weights are allowed. A labeling aaa assigns a label to every vertex and satisfies eee when a(ve)=πe(a(ue))a(v_e)=\pi_e(a(u_e))a(ve​)=πe​(a(ue​)). The value of aaa is the fraction of satisfied constraints and val(G)\mathrm{val}(G)val(G) is the maximum over labelings.

An instance is simple bipartite if its vertices split into two sides, every constraint goes from the first side to the second, and no two list entries join the same ordered pair. It is a translation instance over K=F2sK=\mathbb F_2^sK=F2s​ if every constraint has the form a(ve)=a(ue)+cea(v_e)=a(u_e)+c_ea(ve​)=a(ue​)+ce​.

3SAT inputs are bit strings decoding to a conjunction of clauses with exactly three literal slots (repeated variables and literals allowed); the input is a YES instance when it decodes to a satisfiable formula.

Formalization targets

Goal: the Unique Games Theorem (Theorem 1.1, p. 2)

For every ε,δ∈(0,1/2)\varepsilon,\delta\in(0,1/2)ε,δ∈(0,1/2) there exist s≥1s\ge1s≥1 and a deterministic polynomial-time map φ↦Gφ\varphi\mapsto G_\varphiφ↦Gφ​ to simple bipartite translation instances over F2s\mathbb F_2^sF2s​ such that

φ∈3SAT  ⟹  val(Gφ)≥1−ε,φ∉3SAT  ⟹  val(Gφ)≤δ.\varphi\in\mathrm{3SAT}\;\Longrightarrow\;\mathrm{val}(G_\varphi)\ge 1-\varepsilon,\qquad \varphi\notin\mathrm{3SAT}\;\Longrightarrow\;\mathrm{val}(G_\varphi)\le \delta .φ∈3SAT⟹val(Gφ​)≥1−ε,φ∈/3SAT⟹val(Gφ​)≤δ.

The alphabet, the program and the running-time polynomial depend only on (ε,δ)(\varepsilon,\delta)(ε,δ); the two errors are chosen independently, before the alphabet. In Lean this is OAI.UniqueGamesTheorem.theorem11, which asserts that the structure BinaryGapReduction ε δ is inhabited.

Significance

The result itself. The conjecture is the hypothesis behind a large body of tight approximation thresholds; the preprint's Section 8 (pp. 39–43) records how Theorem 1.1 feeds the earlier reductions for Max-Cut, Vertex Cover, Raghavendra's CSP theorem, ordering CSPs, multicut/sparsest cut and correlation clustering. Those reductions and thresholds are due to their original authors. The statement also fixes the conjecture's quantifier order (errors first, then alphabet), which is exactly what the known algorithms above do not contradict.

Formalizing it. No machine-checked proof of the UGC, or of the 2-to-2 theorem it builds on, is known to exist. A complete development would have to formalize Håstad's parity gap, classical parallel repetition for projection games, the Khot–Minzer–Safra Grassmann/shortcode inverse theorem (used as an external input via Theorem 5.1, p. 22 of the preprint), and the preprint's new latent-alphabet gadget and matrix test. Each of these is reusable well beyond this mission.

Difficulty

The central obstacle, identified in the preprint's introduction (pp. 2–4), is completeness: in a unique constraint every label permits exactly one reply, so the natural codeword tests lose a constant fraction of constraints even on honest proofs. For example, on the rank-one matrix shortcode the evaluation M↦MzM\mapsto MzM↦Mz is preserved by a random rank-one update only with probability about 1/21/21/2. The 2-to-2 theorem gives soundness but only near-1/21/21/2 completeness, and pushing completeness to 1−ε1-\varepsilon1−ε while keeping soundness δ\deltaδ independent of the alphabet is a separate problem that the earlier programme left open.

Formalization scope

  • Instances are Foundations.Target.Instance q over alphabet Fin q, with constraints given by explicit forward/inverse permutation tables and a nonemptiness proof; the value is countSatisfied / constraints.length in ℝ, so division by zero cannot occur.
  • The reduction is a Turing.TM2ComputableInPolyTime program from raw input bits (identity encoding, so runtime is measured in the raw input length including sparse variable names) to the bit encoding gameBits; its tape alphabets must be finite.
  • coordinates : Fin alphabet ≃ (Fin s → ZMod 2) fixes the F2s\mathbb F_2^sF2s​ structure, and translations requires every constraint to be a translation in these coordinates; simpleBipartite gives the orientation and the no-parallel-edge condition on list positions.
  • The 3SAT language is BinaryLanguage.language: inputs that decode (canonically) to a satisfiable three-slot formula. Malformed inputs are NO instances and must be mapped to low-value games.
  • The hypotheses 0<ε,δ<1/20<\varepsilon,\delta<1/20<ε,δ<1/2 are satisfiable and the conclusion is a data-carrying structure, so there is no vacuous reading; the completeness clause requires an actual labeling, and soundness quantifies over all labelings.

Welcome contributions: a reusable Lean library for PCP-style gap reductions on TM2 machines, Fourier analysis over F2n\mathbb F_2^nF2n​, parallel repetition for projection games, and the Grassmann expansion theorem.

Selected references

  • OpenAI, The Unique Games Theorem, OpenAI Math Release preprint, September 23, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-Unique-Games-Theorem-September-23-2026/paper.pdf
  • S. Khot, On the power of unique 2-prover 1-round games, STOC 2002. https://doi.org/10.1145/509907.510017
  • S. Khot, On the Unique Games Conjecture, CCC 2010. https://doi.org/10.1109/CCC.2010.19
  • S. Khot, G. Kindler, E. Mossel, R. O'Donnell, Optimal inapproximability results for MAX-CUT and other 2-variable CSPs?, SIAM J. Comput. 37 (2007). https://doi.org/10.1137/S0097539705447372
  • S. Khot, O. Regev, Vertex cover might be hard to approximate to within 2−ε, J. Comput. Syst. Sci. 74 (2008). https://doi.org/10.1016/j.jcss.2007.06.019
  • P. Raghavendra, Optimal algorithms and inapproximability results for every CSP?, STOC 2008. https://doi.org/10.1145/1374376.1374414
  • J. Håstad, Some optimal inapproximability results, J. ACM 48 (2001). https://doi.org/10.1145/502090.502098
  • B. Barak, P. Kothari, D. Steurer, Small-set expansion in shortcode graph and the 2-to-2 conjecture, ITCS 2019. https://doi.org/10.4230/LIPIcs.ITCS.2019.9
  • S. Khot, D. Minzer, M. Safra, Pseudorandom sets in Grassmann graph have near-perfect expansion, Ann. of Math. 198 (2023). https://doi.org/10.4007/annals.2023.198.1.1
  • S. Arora, B. Barak, D. Steurer, Subexponential algorithms for Unique Games and related problems, FOCS 2010. https://doi.org/10.1109/FOCS.2010.59
2 thms1 active userReviewed
CombinatoricsFunctional Analysis·Captain: wurtle

Finite-Circle Obstructions, Binary Codes, and Histogram Embeddings for Edit DistanceResearch Paper

Motivation

Edit distance ED(x,y)\mathrm{ED}(x,y)ED(x,y), the least number of single-symbol insertions, deletions and substitutions turning xxx into yyy, is the standard distance on strings. Embedding strings into ℓ1\ell_1ℓ1​ with small distortion would make edit distance amenable to the fast nearest-neighbour and sketching tools available for ℓ1\ell_1ℓ1​, so the least achievable distortion is a central quantity in metric embedding theory. A single insertion shifts every later symbol, yet cyclically shifting a word costs only two edits; this tension between position and alignment is what makes edit distance hard to embed.

This preprint gives two independent constructions of binary word sets whose least ℓ1\ell_1ℓ1​ distortion is exp⁡(Ω(log⁡dlog⁡log⁡d))\exp(\Omega(\sqrt{\log d\log\log d}))exp(Ω(logdloglogd​)), two binary coding arguments that transfer large-alphabet constructions to binary, and a complete histogram embedding attaining the matching upper scale. Together they determine the exponential scale of the distortion.

Background

  • 1966, 1974. Levenshtein introduces insertion/deletion codes; Wagner and Fischer describe alignments by increasing traces.
  • 1995. Linial, London and Rabinovich develop finite metric embeddings into ℓ1\ell_1ℓ1​.
  • 1999. Schulman and Zuckerman construct asymptotically good codes for insertions and deletions.
  • 2003. Andoni, Deza, Gupta, Indyk and Raskhodnikova give binary sets with ℓ1\ell_1ℓ1​ distortion approaching 3/23/23/2.
  • 2006. Khot and Naor prove a (log⁡d)1/2−o(1)(\log d)^{1/2 - o(1)}(logd)1/2−o(1) lower bound by Fourier analysis of cuts.
  • 2007. Ostrovsky and Rabani prove the upper bound exp⁡(O(log⁡dlog⁡log⁡d))\exp(O(\sqrt{\log d\log\log d}))exp(O(logdloglogd​)) for fixed-length binary strings.
  • 2009. Krauthgamer and Rabani prove an Ω(log⁡d)\Omega(\log d)Ω(logd) lower bound.
  • 2026. The OpenAI preprint Finite-Circle Obstructions, Binary Codes, and Histogram Embeddings for Edit Distance (dated September 27, 2026), a companion to Edit Distance in ℓ1\ell_1ℓ1​: Matching Bounds up to Constants in the Exponent, proves the results below. It has not been peer reviewed; the Lean goal is open on this platform.

Setting

For a finite alphabet Σ\SigmaΣ with ∣Σ∣≥2|\Sigma| \ge 2∣Σ∣≥2 and an integer d≥1d \ge 1d≥1, let Σ≤d\Sigma^{\le d}Σ≤d be the words of length at most ddd, including the empty word. For an injective f:Σ≤d→ℓ1f : \Sigma^{\le d} \to \ell_1f:Σ≤d→ℓ1​,

dist(f)=(max⁡x≠y∥f(x)−f(y)∥1ED(x,y))(max⁡x≠yED(x,y)∥f(x)−f(y)∥1),EΣ(d)=inf⁡fdist(f).\mathrm{dist}(f) = \Big(\max_{x\ne y}\frac{\|f(x)-f(y)\|_1}{\mathrm{ED}(x,y)}\Big)\Big(\max_{x\ne y}\frac{\mathrm{ED}(x,y)}{\|f(x)-f(y)\|_1}\Big), \qquad E_\Sigma(d) = \inf_f \mathrm{dist}(f).dist(f)=(x=ymax​ED(x,y)∥f(x)−f(y)∥1​​)(x=ymax​∥f(x)−f(y)∥1​ED(x,y)​),EΣ​(d)=finf​dist(f).

The longest common subsequence LCS(x,y)\mathrm{LCS}(x,y)LCS(x,y) and the deficit rL(x,y)=L−LCS(x,y)r_L(x,y) = L - \mathrm{LCS}(x,y)rL​(x,y)=L−LCS(x,y) for words of length LLL are used in the coding results. A binary substitution c:Σ→{0,1}wc : \Sigma \to \{0,1\}^wc:Σ→{0,1}w replaces each letter by a codeword.

In Lean (OAI.FiniteCircle), editDistance is the least script length, lcs and deficit are as above, optimalDistortion is the infimum of distortion over injective maps into lp (fun _ : ℕ => ℝ) 1, E α d is the least distortion of words of length at most ddd, and maskWord, excludedWord are the two finite-circle constructions with prime-period parameters.

Formalization targets

Goal: finiteCircle_fullMain

The conjunction of the paper's main statements, with common absolute constants:

  • Theorem 1.1. For d≥d0d \ge d_0d≥d0​ and every finite Σ\SigmaΣ with ∣Σ∣≥2|\Sigma| \ge 2∣Σ∣≥2,
exp⁡(clog⁡dlog⁡log⁡d)≤EΣ(d)≤exp⁡(Clog⁡dlog⁡log⁡d),\exp\big(c\sqrt{\log d\log\log d}\big) \le E_\Sigma(d) \le \exp\big(C\sqrt{\log d\log\log d}\big),exp(clogdloglogd​)≤EΣ​(d)≤exp(Clogdloglogd​),

and the same bounds for sup⁡∣Σ∣≥2EΣ(d)\sup_{|\Sigma| \ge 2} E_\Sigma(d)sup∣Σ∣≥2​EΣ​(d).

  • Lower constructions (Sections 5–7): for each d≥d0d \ge d_0d≥d0​, both the mask construction and the excluded-prime construction, encoded in binary, give a set of words of one common length n≤dn \le dn≤d with least distortion at least exp⁡(clog⁡dlog⁡log⁡d)\exp(c\sqrt{\log d\log\log d})exp(clogdloglogd​).
  • Theorem 4.1 (two binary substitutions): for every alphabet of A≥2A \ge 2A≥2 letters there are codes c:Σ→{0,1}wc : \Sigma \to \{0,1\}^wc:Σ→{0,1}w with w=O(log⁡2A)w = O(\log 2A)w=O(log2A) and a w rL(x,y)≤rwL(c(x),c(y))≤w rL(x,y)a\,w\,r_L(x,y) \le r_{wL}(c(x),c(y)) \le w\,r_L(x,y)awrL​(x,y)≤rwL​(c(x),c(y))≤wrL​(x,y) for equal-length words, ED(c(x),c(y))≤w ED(x,y)\mathrm{ED}(c(x),c(y)) \le w\,\mathrm{ED}(x,y)ED(c(x),c(y))≤wED(x,y), plus the local separation property each construction needs.
  • Theorem 8.1 (uniform upper bound): a deterministic finite-dimensional map FFF with ED(x,y)≤∥F(x)−F(y)∥1≤exp⁡(Clog⁡dlog⁡log⁡d) ED(x,y)\mathrm{ED}(x,y) \le \|F(x)-F(y)\|_1 \le \exp(C\sqrt{\log d\log\log d})\,\mathrm{ED}(x,y)ED(x,y)≤∥F(x)−F(y)∥1​≤exp(Clogdloglogd​)ED(x,y).

Significance

The result. The lower bounds reach the Ostrovsky–Rabani scale, so the order of log⁡EΣ(d)\log E_\Sigma(d)logEΣ​(d) is log⁡dlog⁡log⁡d\sqrt{\log d\log\log d}logdloglogd​; this improves the previous Ω(log⁡d)\Omega(\log d)Ω(logd) lower bound to exp⁡(Ω(log⁡dlog⁡log⁡d))\exp(\Omega(\sqrt{\log d\log\log d}))exp(Ω(logdloglogd​)). The histogram embedding gives an explicit, alphabet-uniform upper bound including all shorter words.

Formalizing it. No machine-checked proof exists; the source is an unrefereed preprint. Unlike the companion paper, the upper bound here is proved in full (Section 8) rather than imported, so the goal is self-contained relative to the preprint. The edit-distance and distortion definitions are reusable across the edit-distance missions.

Difficulty

Lower bounds for edit distance must control every increasing matching of the two words, including matchings crossing the row boundaries used to construct them. A cheap translation of phases moves row lists at low edit cost, while a simultaneous half-turn at every leaf must be expensive; proving the latter for arbitrary alignments, and then showing that every ℓ1\ell_1ℓ1​ image must distort one of these comparisons, requires combining a combinatorial separation estimate with a Fourier argument on the finite phase space that forces large total frequency.

Formalization scope

  • Words are List α with Fintype α; the domain Σ≤d\Sigma^{\le d}Σ≤d is a subtype by length; the empty word is included.
  • distortion uses sSup over finite sets of ratios; optimalDistortion is sInf over injective maps into ℓ1(N)\ell_1(\mathbb N)ℓ1​(N); exponentScale d = sqrt(log d · log log d) with natural logarithms.
  • The lower witnesses fix the paper's constructions (prime assignments MaskPrimes, ExcludedPrimes, depth parameters) and require a common positive length n≤dn \le dn≤d.
  • UniformUpper asks for a map into Fin M → ℝ with the ℓ1\ell_1ℓ1​ norm written as a finite sum.

The goal is not vacuous: a set with one element would have least distortion 000, so witnesses must be nontrivial, and both bounds are required uniformly.

Needed infrastructure: LCS and alignment combinatorics, cut decompositions of finite ℓ1\ell_1ℓ1​ metrics, Fourier analysis on finite cyclic groups, prime counting estimates, and random partition arguments for the upper bound. Contributions toward any component are welcome.

Selected references

  • OpenAI, Finite-Circle Obstructions, Binary Codes, and Histogram Embeddings for Edit Distance, OpenAI Math Release preprint, September 27, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Finite-Circle-Obstructions-Binary-Codes-and-Histogram-Embeddings-for-Edit-Distance-September-27-2026/paper.pdf
  • OpenAI, Edit Distance in ℓ1\ell_1ℓ1​: Matching Bounds up to Constants in the Exponent, OpenAI Math Release preprint, September 27, 2026. https://github.com/openai/math/blob/main/preprints/Edit-Distance-in-l1-Matching-Bounds-up-to-Constants-in-the-Exponent-September-27-2026/paper.pdf
  • R. Ostrovsky, Y. Rabani, Low distortion embeddings for edit distance, J. ACM 54(5), 2007. https://doi.org/10.1145/1284320.1284322
  • R. Krauthgamer, Y. Rabani, Improved lower bounds for embeddings into L1L_1L1​, SIAM J. Comput. 38(6), 2009.
  • S. Khot, A. Naor, Nonembeddability theorems via Fourier analysis, Math. Ann. 334(4), 2006. https://doi.org/10.1007/s00208-005-0745-0
  • A. Andoni, M. Deza, A. Gupta, P. Indyk, S. Raskhodnikova, Lower bounds for embedding edit distance into normed spaces, SODA 2003.
  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15(2), 1995. https://doi.org/10.1007/BF01200757
  • L. J. Schulman, D. Zuckerman, Asymptotically good codes correcting insertions, deletions, and transpositions, IEEE Trans. Inform. Theory 45(7), 1999. https://doi.org/10.1109/18.796406
  • V. I. Levenshtein, Binary codes capable of correcting deletions, insertions, and reversals, Soviet Physics Doklady 10(8), 1966.
  • R. A. Wagner, M. J. Fischer, The string-to-string correction problem, J. ACM 21(1), 1974. https://doi.org/10.1145/321796.321811
2 thms1 active userReviewed
CombinatoricsFunctional Analysis·Captain: wurtle

Edit Distance in l1: Matching Bounds up to Constants in the ExponentResearch Paper

Motivation

Edit distance ED(x,y)\mathrm{ED}(x,y)ED(x,y) is the least number of single-symbol insertions, deletions and substitutions transforming a string xxx into yyy. It is the basic similarity measure for strings, but it is expensive to compute and awkward to index. A standard way around this is to embed strings into ℓ1\ell_1ℓ1​ so that ℓ1\ell_1ℓ1​ distances approximate edit distances; the quality of such a map is its distortion. Low-distortion embeddings give fast approximate nearest-neighbour search and sketching for edit distance, so the best possible distortion has been studied for two decades.

This preprint determines the exponential scale of the optimal distortion: for strings of length at most ddd, it is exp⁡(Θ(log⁡d log⁡log⁡d))\exp(\Theta(\sqrt{\log d\,\log\log d}))exp(Θ(logdloglogd​)), uniformly over finite alphabets.

Background

  • 1966. Levenshtein studies codes correcting insertions, deletions and substitutions.
  • 1974. Wagner and Fischer describe edit scripts through increasing traces (alignments).
  • 1995. Linial, London and Rabinovich develop the geometry of finite metric embeddings into ℓ1\ell_1ℓ1​.
  • 2003. Andoni, Deza, Gupta, Indyk and Raskhodnikova give binary subsets with ℓ1\ell_1ℓ1​ distortion approaching 3/23/23/2.
  • 2006. Khot and Naor prove an Ω(log⁡d/log⁡log⁡d)\Omega(\sqrt{\log d/\log\log d})Ω(logd/loglogd​) lower bound via cuts and Fourier analysis of noisy shifts.
  • 2007. Ostrovsky and Rabani prove the upper bound exp⁡(O(log⁡dlog⁡log⁡d))\exp(O(\sqrt{\log d\log\log d}))exp(O(logdloglogd​)) for binary strings of length ddd.
  • 2009. Krauthgamer and Rabani prove an Ω(log⁡d)\Omega(\log d)Ω(logd) lower bound for binary strings of length ddd.
  • 2026. The OpenAI preprint Edit Distance in ℓ1\ell_1ℓ1​: Matching Bounds up to Constants in the Exponent (dated September 27, 2026) proves a matching lower bound exp⁡(Ω(log⁡dlog⁡log⁡d))\exp(\Omega(\sqrt{\log d\log\log d}))exp(Ω(logdloglogd​)). It has not been peer reviewed; the Lean goal is open on this platform.

Setting

For a finite alphabet Σ\SigmaΣ and d≥1d \ge 1d≥1, let Σ≤d\Sigma^{\le d}Σ≤d be the set of strings of length at most ddd, including the empty string. Edit scripts may pass through strings of any length. Let ℓ1\ell_1ℓ1​ be the space of absolutely summable real sequences. Define

EΣ(d)=inf⁡f:Σ≤d↪ℓ1(max⁡x≠y∥f(x)−f(y)∥1ED(x,y))(max⁡x≠yED(x,y)∥f(x)−f(y)∥1),E_\Sigma(d) = \inf_{f:\Sigma^{\le d}\hookrightarrow \ell_1}\Big(\max_{x\ne y}\frac{\|f(x)-f(y)\|_1}{\mathrm{ED}(x,y)}\Big)\Big(\max_{x\ne y}\frac{\mathrm{ED}(x,y)}{\|f(x)-f(y)\|_1}\Big),EΣ​(d)=f:Σ≤d↪ℓ1​inf​(x=ymax​ED(x,y)∥f(x)−f(y)∥1​​)(x=ymax​∥f(x)−f(y)∥1​ED(x,y)​),

the infimum over injective maps with no restriction on dimension or computability.

In Lean (OAI.EditDistortion), edit x y is the least length of a Script true n x y of unit insertions, deletions and substitutions on List Alpha; L1 = lp (fun _ : ℕ => ℝ) 1; distortion ρ f is the product of the two maximal ratios for an embedding f : X ↪ L1; leastDistortion is the infimum over embeddings; E Alpha d applies this to {x : List Alpha // x.length ≤ d}; exponentScale c d = exp(c·sqrt(log d · log log d)).

Formalization targets

Goal: Theorem 1.1

There are absolute c,C>0c, C > 0c,C>0 and d0d_0d0​ such that for every d≥d0d \ge d_0d≥d0​ and every finite alphabet with ∣Σ∣≥2|\Sigma| \ge 2∣Σ∣≥2,

exp⁡(clog⁡d log⁡log⁡d)≤EΣ(d)≤exp⁡(Clog⁡d log⁡log⁡d).\exp\big(c\sqrt{\log d\,\log\log d}\big) \le E_\Sigma(d) \le \exp\big(C\sqrt{\log d\,\log\log d}\big).exp(clogdloglogd​)≤EΣ​(d)≤exp(Clogdloglogd​).

The lower bound is witnessed by a finite set of binary strings of one common length at most ddd, and the same two bounds hold for sup⁡2≤∣Σ∣<∞EΣ(d)\sup_{2 \le |\Sigma| < \infty} E_\Sigma(d)sup2≤∣Σ∣<∞​EΣ​(d).

Significance

The result. Theorem 1.1 shows that the Ostrovsky–Rabani embedding is optimal up to the constant in the exponent, closing the gap between Ω(log⁡d)\Omega(\log d)Ω(logd) (Krauthgamer–Rabani) and exp⁡(O(log⁡dlog⁡log⁡d))\exp(O(\sqrt{\log d\log\log d}))exp(O(logdloglogd​)). It rules out any ℓ1\ell_1ℓ1​ embedding of edit distance with distortion exp⁡(o(log⁡dlog⁡log⁡d))\exp(o(\sqrt{\log d\log\log d}))exp(o(logdloglogd​)), which limits ℓ1\ell_1ℓ1​-based approaches to approximate edit-distance search.

Formalizing it. No machine-checked proof exists; the source is an unrefereed preprint. The upper bound relies on Ostrovsky and Rabani's fixed-length theorem (Theorem 5.1 in the preprint, taken from the literature) plus reductions proved in the preprint, so a full formalization also requires formalizing the Ostrovsky–Rabani embedding. The lower bound (Sections 2–4) is the new contribution.

Difficulty

Any lower bound must control every alignment of two strings, including matchings that cross the boundaries of the blocks used to build them, because an insertion or deletion shifts an entire suffix. Earlier Fourier-analytic arguments (Khot–Naor) and recursive constructions (Krauthgamer–Rabani) lose too much to reach the log⁡dlog⁡log⁡d\sqrt{\log d\log\log d}logdloglogd​ exponent; reaching it requires a construction whose combinatorial separation and an ℓ1\ell_1ℓ1​ displacement inequality match at that scale.

Formalization scope

  • Strings are List Alpha with Fintype Alpha, 2 ≤ card Alpha; the domain includes the empty string and all lengths up to ddd.
  • Edit distance includes substitutions; edit is an sInf over script lengths, which is attained since scripts always exist.
  • distortion uses sSup over finitely many ratios (the domain is finite), so no default values arise; leastDistortion is sInf over all injective maps into ℓ1(N)\ell_1(\mathbb{N})ℓ1​(N).
  • BinaryWitness c d asks for a finite nonempty set of Boolean lists of common length n≤dn \le dn≤d whose least distortion is at least exponentScale c d; a singleton has least distortion 000, so the witness must be nontrivial.

The statement cannot be met trivially: both bounds are required with constants uniform in the alphabet.

Needed infrastructure: cut decompositions of finite ℓ1\ell_1ℓ1​ metrics, Fourier analysis on finite abelian groups, prime-number estimates (Rosser–Schoenfeld), and the Ostrovsky–Rabani construction. Contributions toward the displacement inequality (Proposition 2.2) or the recursive construction (Proposition 3.1) are welcome.

Selected references

  • OpenAI, Edit Distance in ℓ1\ell_1ℓ1​: Matching Bounds up to Constants in the Exponent, OpenAI Math Release preprint, September 27, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Edit-Distance-in-l1-Matching-Bounds-up-to-Constants-in-the-Exponent-September-27-2026/paper.pdf
  • R. Ostrovsky, Y. Rabani, Low distortion embeddings for edit distance, J. ACM 54(5), Article 23, 2007. https://doi.org/10.1145/1284320.1284322
  • R. Krauthgamer, Y. Rabani, Improved lower bounds for embeddings into L1L_1L1​, SIAM J. Comput. 38(6), 2009.
  • S. Khot, A. Naor, Nonembeddability theorems via Fourier analysis, Math. Ann. 334(4), 2006. https://doi.org/10.1007/s00208-005-0745-0
  • A. Andoni, M. Deza, A. Gupta, P. Indyk, S. Raskhodnikova, Lower bounds for embedding edit distance into normed spaces, SODA 2003.
  • N. Linial, E. London, Y. Rabinovich, The geometry of graphs and some of its algorithmic applications, Combinatorica 15(2), 1995. https://doi.org/10.1007/BF01200757
  • V. I. Levenshtein, Binary codes capable of correcting deletions, insertions, and reversals, Soviet Physics Doklady 10(8), 1966.
  • R. A. Wagner, M. J. Fischer, The string-to-string correction problem, J. ACM 21(1), 1974. https://doi.org/10.1145/321796.321811
2 thms1 active userReviewed
CombinatoricsDiscrete Geometry·Captain: wurtle

The Euclidean Steinitz–Bergström theoremResearch Paper

Motivation: keeping partial sums of vectors small

Given vectors v1,…,vNv_1,\dots,v_Nv1​,…,vN​ in the Euclidean unit ball of Rd\mathbb R^dRd that sum to zero, can they always be reordered so that every partial sum stays in a ball whose radius depends only on ddd? Steinitz's lemma, from his work on rearrangements of conditionally convergent vector series, says yes: in any norm, radius ddd suffices (Grinberg–Sevast'yanov 1980). For the Euclidean norm, the expected answer is O(d)O(\sqrt d)O(d​), the Euclidean Steinitz–Bergström conjecture. A closely related prefix discrepancy problem asks for signs εi∈{±1}\varepsilon_i\in\{\pm1\}εi​∈{±1}, chosen once and for all, such that every signed prefix ∑i≤kεivi\sum_{i\le k}\varepsilon_iv_i∑i≤k​εi​vi​ is O(d)O(\sqrt d)O(d​), independently of NNN. Both questions are basic in discrepancy theory and combinatorial vector balancing, with applications to scheduling, rounding and the analysis of online and streaming algorithms.

Timeline

  • 1913 — Steinitz proves his rearrangement lemma for conditionally convergent vector series (historical account in Ambrus–Heck 2026, Section 2).
  • 1954 — Behrend discusses the expected square-root growth in Euclidean space (Canad. J. Math., p. 108).
  • 1980 — Grinberg and Sevast'yanov prove that partial sums can be kept within ddd times the unit ball, for every norm (Funct. Anal. Appl.).
  • 1994–2023 — Chobanyan's transference principle relates rearrangements to signed sums (Probability in Banach Spaces 9, 1994), later in a finite positive-forward, negative-reverse form (Chobanyan et al., Bull. TICMI 2023).
  • 1998 / 2012 — Banaszczyk's Gaussian-measure balancing theorem (Random Structures Algorithms 1998) and his prefix bound O(d+log⁡N)O(\sqrt d+\sqrt{\log N})O(d​+logN​) for signed series and rearrangements (Random Structures Algorithms 2012).
  • 2021 — Bansal, Jiang, Meka, Singla and Sinha state the O(d)O(\sqrt d)O(d​) prefix-signing question explicitly (Conjecture 6.3, arXiv:2111.07049).
  • 2026 — Ambrus and Heck record the conjecture and its attribution to Bergström (Conjecture 5) and reduce it to a relaxed problem (Mathematika); Dutta, Jha and Jiang give efficient bounds O(d+d1/4log⁡7/4N)O(\sqrt d+d^{1/4}\log^{7/4}N)O(d​+d1/4log7/4N) (arXiv:2604.13355).
  • 2026 — An OpenAI preprint, The Euclidean Steinitz–Bergström theorem (OpenAI Math Release, September 24, 2026), claims the O(d)O(\sqrt d)O(d​) bound for both problems with an absolute constant. It has not been peer reviewed, and its proof is not formally verified.

Setting

Work in Euclidean space Rd\mathbb R^dRd with norm ∥⋅∥2\|\cdot\|_2∥⋅∥2​, and let v1,…,vNv_1,\dots,v_Nv1​,…,vN​ satisfy ∥vi∥2≤1\|v_i\|_2\le1∥vi​∥2​≤1; repetitions and zero vectors are allowed. A prefix sum is ∑i=1k\sum_{i=1}^k∑i=1k​ for 0≤k≤N0\le k\le N0≤k≤N (the empty sum is 000). For a zero-sum family, the Steinitz quantity

β(v1,…,vN)=min⁡π∈SN max⁡0≤k≤N∥∑i=1kvπ(i)∥2\beta(v_1,\dots,v_N)=\min_{\pi\in\mathfrak S_N}\ \max_{0\le k\le N}\Bigl\|\sum_{i=1}^kv_{\pi(i)}\Bigr\|_2β(v1​,…,vN​)=π∈SN​min​ 0≤k≤Nmax​​i=1∑k​vπ(i)​​2​

is the best achievable maximal prefix norm over all orderings. The Euclidean Steinitz constant S2(d)S_2(d)S2​(d) is its supremum over all such families.

Formalization targets

Goal: prescribed-order signing and Euclidean Steinitz with one absolute constant (Theorems 1.1 and 1.2)

There is an absolute constant CCC such that for all d,N≥1d,N\ge1d,N≥1 and all v1,…,vN∈Rdv_1,\dots,v_N\in\mathbb R^dv1​,…,vN​∈Rd with ∥vi∥2≤1\|v_i\|_2\le1∥vi​∥2​≤1:

  1. (Theorem 1.1) there are signs εi∈{−1,1}\varepsilon_i\in\{-1,1\}εi​∈{−1,1} with
max⁡0≤k≤N∥∑i=1kεivi∥2≤Cd;\max_{0\le k\le N}\Bigl\|\sum_{i=1}^k\varepsilon_iv_i\Bigr\|_2\le C\sqrt d ;0≤k≤Nmax​​i=1∑k​εi​vi​​2​≤Cd​;
  1. (Theorem 1.2) if moreover ∑ivi=0\sum_iv_i=0∑i​vi​=0, some permutation π\piπ satisfies
max⁡0≤k≤N∥∑i=1kvπ(i)∥2≤Cd.\max_{0\le k\le N}\Bigl\|\sum_{i=1}^kv_{\pi(i)}\Bigr\|_2\le C\sqrt d .0≤k≤Nmax​​i=1∑k​vπ(i)​​2​≤Cd​.

The constant is not fixed. The simplex example in the source shows S2(d)≥12dS_2(d)\ge\tfrac12\sqrt dS2​(d)≥21​d​, so the order is optimal. The goal statement is published on the platform with status Open.

Significance

The result itself. Theorem 1.2 settles the Euclidean Steinitz–Bergström conjecture: S2(d)=Θ(d)S_2(d)=\Theta(\sqrt d)S2​(d)=Θ(d​). Theorem 1.1 settles the ℓ2\ell_2ℓ2​ prefix-discrepancy question of Bansal et al., removing the log⁡N\sqrt{\log N}logN​ term in Banaszczyk's bound. A single choice of signs controls all prefixes at once, unlike terminal balancing results that bound one sum. The source also derives finite-dimensional ℓp\ell_pℓp​ Steinitz bounds and a colorful Steinitz theorem (Corollaries 8.1–8.2). The result is existential and gives no efficient algorithm.

Formalizing it. The statement is finite and elementary. The deduction of Theorem 1.2 from Theorem 1.1 is a short transference argument (Section 1.1 of the source), suitable as a first formalization step. The proof of Theorem 1.1 uses Gaussian measures of convex bodies, Dirichlet energies, Steiner symmetrization and Brownian survival estimates, which would make substantial reusable additions to Mathlib's probability and convex-geometry libraries.

Difficulty

Signing each prefix separately is easy, by terminal balancing results, but one common signing must work for all NNN prefixes. Union bounds over prefixes and Banaszczyk's Gaussian-measure method lose a log⁡N\sqrt{\log N}logN​ factor, because the probability that a random walk stays in a ball of radius CdC\sqrt dCd​ for NNN steps decays with NNN. The proof must therefore control a single body of coefficients enforcing all prefix constraints, with bounds independent of the number of steps.

Formalization scope

  • Vectors live in EuclideanSpace ℝ (Fin d), indexed by Fin N, with d,N≥1d,N\ge1d,N≥1 and ∥vi∥≤1\|v_i\|\le1∥vi​∥≤1.
  • Prefixes are sums over {i:i<k}\{i : i < k\}{i:i<k} for every natural k≤Nk\le Nk≤N, so the empty prefix and the full sum are included.
  • SignedPrefixBound C asks for signs εi∈{−1,1}\varepsilon_i\in\{-1,1\}εi​∈{−1,1} (as reals) with every signed prefix at most CdC\sqrt dCd​. OrderingPrefixBound C asks, for zero-sum families, for a permutation π\piπ with every reordered prefix at most CdC\sqrt dCd​. The goal is ∃C, SignedPrefixBound C∧OrderingPrefixBound C\exists C,\ \mathrm{SignedPrefixBound}\,C\wedge\mathrm{OrderingPrefixBound}\,C∃C, SignedPrefixBoundC∧OrderingPrefixBoundC; the same constant serves both, as stated in the source.
  • The quantifier order is the source's: CCC is chosen before ddd, NNN and the vectors.
  • Needed infrastructure: Gaussian measures of symmetric convex bodies, Steiner symmetrization, Dirichlet energies of densities, Ornstein–Uhlenbeck/Brownian survival estimates, and the Chobanyan transference. Contributions proving the transference step (Theorem 1.1 ⇒ Theorem 1.2) or the lower-bound example are welcome.

Selected references

  • V. S. Grinberg and S. V. Sevast'yanov, Value of the Steinitz constant, Funct. Anal. Appl. (1980).
  • F. A. Behrend, The Steinitz–Gross theorem on sums of vectors, Canad. J. Math. (1954).
  • S. Chobanyan, Convergence a.s. of rearranged random series in Banach space and associated inequalities, Probability in Banach Spaces 9 (1994).
  • W. Banaszczyk, Balancing vectors and Gaussian measures of n-dimensional convex bodies, Random Structures Algorithms (1998).
  • W. Banaszczyk, On series of signed vectors and their rearrangements, Random Structures Algorithms (2012).
  • N. Bansal, H. Jiang, R. Meka, S. Singla and M. Sinha, Prefix discrepancy, smoothed analysis, and combinatorial vector balancing, 2021. https://arxiv.org/abs/2111.07049
  • G. Ambrus and R. Heck, A note on the Steinitz Lemma, Mathematika (2026). https://doi.org/10.1112/mtk.70085
  • K. Dutta, A. V. Jha and H. Jiang, Near-optimal constructive bounds for ℓ2\ell_2ℓ2​ prefix discrepancy and Steinitz problems via affine spectral independence, 2026. https://arxiv.org/abs/2604.13355
  • OpenAI, The Euclidean Steinitz–Bergström theorem, OpenAI Math Release preprint, September 24, 2026 (Theorems 1.1 and 1.2, p. 2). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-Euclidean-Steinitz-Bergstrom-theorem-September-24-2026/The-Euclidean-Steinitz-Bergstrom-theorem-September-24-2026.pdf
2 thms1 active userReviewed
Discrete GeometryProbability·Captain: wurtle

The Gaussian propeller bound in every dimensionResearch Paper

Motivation

Split Euclidean space into finitely many measurable cells and record, for each cell, its Gaussian first moment, the integral of xxx against the standard Gaussian measure over that cell. How large can the sum of the squared lengths of these vectors be? The propeller conjecture predicts that, with at least three cells in dimension at least two, the answer is 9/(8π)9/(8\pi)9/(8π), attained by three planar sectors of angle 2π/32\pi/32π/3 (a "propeller") times the orthogonal complement. The question arose in Khot and Naor's work on approximate kernel clustering (Mathematika 2009), where this Gaussian partition parameter determines both the approximation ratio of Gaussian rounding and the matching Unique-Games hardness threshold. The surprise is that more cells or more dimensions do not help: the optimum uses only three cells in a plane.

This mission asks for a formal proof of the propeller bound in every dimension, 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

  • 1983–2003 — The Gaussian Minkowski (Ehrhard) inequality: Ehrhard for convex sets (Math. Scand. 1983), Latała when one set is convex (Studia Math. 1996), Borell for Borel sets (C. R. Acad. Sci. 2003); later van Handel's game proof (PTRF 2018).
  • 2009 — Khot and Naor introduce the problem, reduce maximizers to conical partitions in the span of their centroids, and compute the three-cell value (Mathematika 2009); see also their sharp kernel-clustering paper (RSA 2013).
  • 2013 — Heilman, Jagannath and Naor prove the conjecture in R3\mathbb R^3R3 with a computer-assisted argument (DCG 2013).
  • 2022 — Heilman gives a conditional route for three or four cells in any dimension, under a stability hypothesis on noise-stability maximizers (arXiv:2209.11216).
  • September 2026 — The OpenAI preprint claims the bound for every dimension and every number of cells (Theorem 1.1, p. 1).

Setting

Let γd\gamma_dγd​ be the standard Gaussian probability measure on Rd\mathbb R^dRd, with density (2π)−d/2e−∥x∥2/2(2\pi)^{-d/2}e^{-\|x\|^2/2}(2π)−d/2e−∥x∥2/2. For a measurable set A⊆RdA\subseteq\mathbb R^dA⊆Rd, its centroid is the unnormalized vector

z(A)=∫Ax dγd(x)∈Rd.z(A)=\int_A x\,d\gamma_d(x)\in\mathbb R^d .z(A)=∫A​xdγd​(x)∈Rd.

A measurable partition (A1,…,Ak)(A_1,\dots,A_k)(A1​,…,Ak​) of Rd\mathbb R^dRd consists of measurable sets such that almost every point (for γd\gamma_dγd​) lies in exactly one of them; empty cells and arbitrary cell probabilities are allowed. Its value is F(A)=∑i=1k∥z(Ai)∥2F(A)=\sum_{i=1}^k\|z(A_i)\|^2F(A)=∑i=1k​∥z(Ai​)∥2. The propeller is the partition of Rd\mathbb R^dRd, d≥2d\ge2d≥2, into the three sets Sj×Rd−2S_j\times\mathbb R^{d-2}Sj​×Rd−2, where S1,S2,S3S_1,S_2,S_3S1​,S2​,S3​ are the planar sectors of opening 2π/32\pi/32π/3 bounded by the rays at angles ±π/3\pm\pi/3±π/3 and π\piπ, with any further cells empty.

Formalization targets

Goal: Theorem 1.1 (p. 1)

For all positive integers d,kd,kd,k and every measurable partition (A1,…,Ak)(A_1,\dots,A_k)(A1​,…,Ak​) of Rd\mathbb R^dRd,

∑i=1k∥∫Aix dγd(x)∥2 ≤ 98π,\sum_{i=1}^k\Bigl\|\int_{A_i}x\,d\gamma_d(x)\Bigr\|^2\ \le\ \frac{9}{8\pi},i=1∑k​​∫Ai​​xdγd​(x)​2 ≤ 8π9​,

and for d≥2d\ge2d≥2, k≥3k\ge3k≥3 the propeller is a measurable partition attaining equality.

Significance

The result itself. It settles the propeller conjecture in all dimensions and for all numbers of cells, without constraints on cell masses. Through Khot and Naor's analysis, the constant 9/(8π)9/(8\pi)9/(8π) fixes the loss factor αk=8π9(1−1k)\alpha_k=\frac{8\pi}{9}(1-\frac1k)αk​=98π​(1−k1​) of Gaussian rounding for identity-target kernel clustering; combined with the companion Unique Games preprint, the preprint deduces NP-hardness of any better fixed loss factor (Theorem 6.1, p. 18). The bound also gives a sharp inequality for expected Gaussian maxima (Corollary 5.1, p. 17).

Formalizing it. A formal proof needs Gaussian measures on Euclidean space (available in Mathlib as stdGaussian), a compactness argument for extremal partitions, the convex case of the Ehrhard inequality, and the three-dimensional theorem of Heilman, Jagannath and Naor, which the preprint uses as an established input and which itself relies on a finite computer verification. Each of these is a reusable component; the three-dimensional case alone would be a substantial formalization.

Difficulty

The optimization runs over all measurable partitions, with both shapes and probabilities free. Khot and Naor's reduction gives conical extremizers whose mmm nonzero centroids span an (m−1)(m-1)(m−1)-dimensional space, so the three-dimensional theorem disposes of m≤4m\le4m≤4 active cells. The remaining cases m≥5m\ge5m≥5 cannot be handled by any fixed finite computation. The preprint constrains them in two ways: deleting one score costs a definite amount of expected maximum, via concavity from the Ehrhard inequality (Section 3), and competition against two residual scores constrains a covariance determinant (Section 4, Lemma 4.1 and Corollary 4.2, pp. 11–12). Combining these must rule out every configuration with five or more cells.

Formalization scope

  • Space d is EuclideanSpace ℝ (Fin d) and the measure is stdGaussian. IsPartition A asks each cell to be measurable and almost every point to lie in exactly one cell, matching the paper's convention that partitions are taken up to null sets.
  • centroid A is the Bochner integral of the identity over A; value sums squared Euclidean norms. Integrability holds because Gaussian measures have finite first moments.
  • The goal is a conjunction: the bound for all d,k≥1d,k\ge1d,k≥1, and, for d≥2d\ge2d≥2, k≥3k\ge3k≥3, that the explicit three-sector propeller d k (cells ≥3\ge3≥3 empty) is a partition with value exactly 9/(8π)9/(8\pi)9/(8π). The sectors are written with closed inequalities, so they overlap only on null rays.
  • No hypothesis is vacuous: the case k=1k=1k=1 is included (value 000), and the equality clause pins the constant.

Selected references

  • OpenAI, The Gaussian propeller bound in every dimension, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/The-Gaussian-Propeller-Bound-in-Every-Dimension-September-24-2026/main.pdf
  • S. Khot, A. Naor, Approximate kernel clustering, Mathematika 55 (2009), 129–165. https://doi.org/10.1112/S002557930000098X
  • S. Khot, A. Naor, Sharp kernel clustering algorithms and their associated Grothendieck inequalities, Random Struct. Algorithms 42 (2013), 269–300. https://doi.org/10.1002/rsa.20398
  • S. Heilman, A. Jagannath, A. Naor, Solution of the Propeller Conjecture in R3\mathbb R^3R3, Discrete Comput. Geom. 50 (2013), 263–305. https://doi.org/10.1007/s00454-013-9530-0
  • S. Heilman, Hyperstable sets with voting and algorithmic hardness applications, arXiv:2209.11216 (2022). https://arxiv.org/abs/2209.11216v1
  • A. Ehrhard, Symétrisation dans l'espace de Gauss, Math. Scand. 53 (1983), 281–301. https://doi.org/10.7146/math.scand.a-12035
  • R. Latała, A note on the Ehrhard inequality, Studia Math. 118 (1996), 169–174. https://doi.org/10.4064/sm-118-2-169-174
  • C. Borell, The Ehrhard inequality, C. R. Acad. Sci. Paris 337 (2003), 663–666. https://doi.org/10.1016/j.crma.2003.09.031
  • R. van Handel, The Borell–Ehrhard game, Probab. Theory Relat. Fields 170 (2018), 555–585. https://doi.org/10.1007/s00440-017-0762-4
  • E. Mossel, R. O'Donnell, K. Oleszkiewicz, Noise stability of functions with low influences: invariance and optimality, Ann. of Math. 171 (2010), 295–341. https://doi.org/10.4007/annals.2010.171.295
2 thms1 active userReviewed
Discrete GeometryGraph Theory·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 thms1 active userReviewed
PreviousPage 8 of 10Next

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