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 · 237 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

Open227Completed237All464
AlgebraAnalysisCombinatorics+3·Captain: Lucas

Formal Conjectures Portfolio: Bateman-Horn and CompanionsOpen Problem

1. Motivation

Wikipedia's pages on open problems are, for many mathematicians, the first contact with a conjecture: a one-paragraph statement, a short history, a list of partial results. The Formal Conjectures library (Google DeepMind, Apache-2.0) turned a large part of that material into Lean 4 statements, so that the conjectures can be attacked — and, just as importantly, stated unambiguously — by machine.

This mission ports a coherent slice of that material to Prove2Me. It is deliberately a portfolio mission: the goal theorem is the Bateman–Horn conjecture, the strongest single statement in the collection, and the milestone list gathers the other conjectures and the landmark theorems that surround them. Some milestones are genuine steps toward the goal (the Bunyakovsky conjecture is literally the one-polynomial case); most are independent open problems from other fields, grouped here because they share a source, a level of difficulty, and a need for faithful formal statements. A reader should not assume that proving a milestone advances the goal theorem. The mission's value is that every statement in it has been written against the same Mathlib revision, checked to compile, and documented well enough to be attacked.

A rough timeline of the collection's landmarks:

  • 1947 — Mills: a real A>1A>1A>1 with ⌊A3n⌋\lfloor A^{3^n}\rfloor⌊A3n⌋ always prime.
  • 1962 — Radó: the busy beaver function outgrows every computable function.
  • 1971 — Davies: planar Kakeya sets have Hausdorff dimension 222.
  • 1978 — Apéry: ζ(3)\zeta(3)ζ(3) is irrational.
  • 1985 — Read (after Enflo, 1981): an operator on ℓ1\ell^1ℓ1 with no nontrivial closed invariant subspace.
  • 2001 — Zudilin: one of ζ(5),ζ(7),ζ(9),ζ(11)\zeta(5),\zeta(7),\zeta(9),\zeta(11)ζ(5),ζ(7),ζ(9),ζ(11) is irrational.
  • 2002 — Mihăilescu: 888 and 999 are the only consecutive perfect powers (Catalan's conjecture).
  • 2009 / 2021 — Dvir; Bukh–Chao: the finite-field Kakeya bound and its sharp density constant.
  • 2021 — Gardam: Kaplansky's unit conjecture is false (its zero-divisor and idempotent companions remain open).
  • 2024 — Saito: Mills' constant is irrational; bbchallenge: BB(5)=47 176 870\mathrm{BB}(5)=47\,176\,870BB(5)=47176870.
  • 2025 — Wang–Zahl: the Kakeya set conjecture in R3\mathbb{R}^3R3.

2. Setting

The goal theorem concerns prime values of polynomials. Fix a finite set S={f1,…,fk}⊆Z[X]S=\{f_1,\dots,f_k\}\subseteq\mathbb{Z}[X]S={f1​,…,fk​}⊆Z[X] of distinct polynomials. Say that fff satisfies the Bunyakovsky condition if its leading coefficient is positive, deg⁡f≥1\deg f\ge 1degf≥1, and fff is irreducible over Z\mathbb{Z}Z; say that SSS satisfies the Schinzel condition if for every prime ppp there is an integer nnn with p∤f1(n)⋯fk(n)p\nmid f_1(n)\cdots f_k(n)p∤f1​(n)⋯fk​(n) — i.e. no fixed prime divides the product at every argument.

For a prime ppp let ωp(S)\omega_p(S)ωp​(S) be the number of residue classes n mod pn \bmod pnmodp at which some fif_ifi​ vanishes, let D=∏ideg⁡fiD=\prod_i \deg f_iD=∏i​degfi​, and let

πS(x)=#{ n≤x:∣fi(n)∣ is prime for every i }.\pi_S(x)=\#\{\,n\le x : |f_i(n)| \text{ is prime for every } i\,\}.πS​(x)=#{n≤x:∣fi​(n)∣ is prime for every i}.

The Bateman–Horn constant is the (conditionally convergent) Euler product

C=lim⁡N→∞ ∏p<N(1−1p)−k(1−ωp(S)p).C=\lim_{N\to\infty}\ \prod_{p<N}\Big(1-\tfrac1p\Big)^{-k}\Big(1-\tfrac{\omega_p(S)}{p}\Big).C=N→∞lim​ p<N∏​(1−p1​)−k(1−pωp​(S)​).

The other groups use their own vocabulary, each fixed in a definition item of this mission: Kakeya sets in Rn\mathbb{R}^nRn and over Fq\mathbb{F}_qFq​; Mills' property ⌊A3n⌋∈P\lfloor A^{3^n}\rfloor \in \mathbb{P}⌊A3n⌋∈P; Wagstaff primes and Catalan–Mersenne numbers; polynomial self-maps and their Jacobian matrix; nontrivial closed invariant subspaces; linear extensions of a finite poset; Catalan's constant; and an explicit two-symbol Turing machine model with its maximum-shifts function BB\mathrm{BB}BB.

3. Target

The goal theorem is the Bateman–Horn asymptotic: under the Bunyakovsky and Schinzel hypotheses, CCC exists and is positive and

πS(x) ∼ CD x(log⁡x)k(x→∞).\pi_S(x)\ \sim\ \frac{C}{D}\,\frac{x}{(\log x)^{k}}\qquad (x\to\infty).πS​(x) ∼ DC​(logx)kx​(x→∞).

Weaker statements in the same direction appear as milestones, first of all Bunyakovsky's conjecture: under the same hypotheses with k=1k=1k=1, fff takes prime values infinitely often. The remaining milestones are listed in the milestone panel and are grouped by subject: Diophantine equations (Brocard, Pillai, Lebesgue–Nagell, Catalan/Mihăilescu), Mersenne-type primality (New Mersenne, infinitude of Mersenne primes, Catalan–Mersenne), prime-representing constants (Mills), geometric measure theory (Kakeya in Rn\mathbb{R}^nRn, Kakeya over Fq\mathbb{F}_qFq​, Falconer), operator theory (invariant subspace problem and Read's ℓ1\ell^1ℓ1 counterexample), group algebras (Kaplansky's zero-divisor and idempotent conjectures), affine algebraic geometry (the two-variable Jacobian conjecture), irrationality and transcendence (ζ(5)\zeta(5)ζ(5), all odd zeta values, Zudilin's theorem, e+πe+\pie+π, eπe\pieπ, γ\gammaγ, Catalan's constant), order theory (the 1/31/31/3–2/32/32/3 conjecture), and computability (Radó's theorem).

4. Significance

The results themselves. Bateman–Horn is the quantitative form of Schinzel's hypothesis H: it contains the twin prime conjecture, the infinitude of primes of the form n2+1n^2+1n2+1, and Bunyakovsky as special cases, and it is the standard heuristic behind prime-counting predictions. The other targets are each the headline question of their area: whether every bounded Hilbert-space operator has an invariant subspace; whether group algebras of torsion-free groups are domains; whether Kakeya sets must have full dimension. The solved milestones (Mihăilescu, Davies, Dvir, Zudilin, Read, Saito, Radó) are landmarks whose formal proofs would be significant library contributions in their own right.

Formalizing them. None of the open statements is expected to fall here; the concrete deliverable is a set of faithful, compiling, reusable statements plus formal proofs of the solved milestones, most of which are not in Mathlib today. Several are realistically in reach: the finite-field Kakeya bound (Dvir's polynomial method is short), the elementary fact that π+e\pi+eπ+e and πe\pi eπe cannot both be algebraic, and Radó's diagonal argument.

5. Difficulty

For Bateman–Horn, the obstruction is visible already for k=1k=1k=1, deg⁡f=2\deg f = 2degf=2: sieve methods bound πS(x)\pi_S(x)πS​(x) from above by a constant times the conjectured main term and produce almost-primes, but the parity problem blocks every known sieve from producing a single prime value of an irreducible quadratic. The conditional convergence of the Euler product is a second, smaller trap: the product over p<Np<Np<N must be taken in order, so any reformulation as an unordered infinite product changes the statement.

Each other group has its own obstruction, and they do not transfer: the parity problem says nothing about Kakeya, where the difficulty is that dimension is not stable under the natural compactness arguments, nor about the invariant subspace problem, where the known counterexamples on ℓ1\ell^1ℓ1 show that no soft argument can work.

6. Formalization scope

Conventions this mission commits to, all fixed in the definition items:

  • Polynomials are elements of ℤ[X]; primality of a polynomial value is primality of its absolute value, and the counting function ranges over natural numbers n≤⌊x⌋n \le \lfloor x\rfloorn≤⌊x⌋.
  • The Bateman–Horn constant is the limit of the ordered partial products over p<Np<Np<N, not an unordered infinite product.
  • Kakeya sets carry no compactness or measurability hypothesis, matching the source; the conjecture is stated as an equality of Hausdorff dimensions in [0,∞][0,\infty][0,∞].
  • Falconer's hypothesis is written d<2dim⁡HEd < 2\dim_H Ed<2dimH​E to avoid division in [0,∞][0,\infty][0,∞].
  • Torsion-freeness of a group is spelled out as "every element of finite order is the identity", which is the hypothesis the source intends (it is weaker than Mathlib's IsMulTorsionFree).
  • Linear extensions are order-preserving bijections onto {0,…,∣P∣−1}\{0,\dots,|P|-1\}{0,…,∣P∣−1}, and probabilities are quotients of set cardinalities in Q\mathbb{Q}Q.
  • The busy beaver model is an explicit nnn-state, 222-symbol machine with a bi-infinite Boolean tape; BB\mathrm{BB}BB counts transitions performed (maximum shifts), the halting transition included, and BB(0)=0\mathrm{BB}(0)=0BB(0)=0.
  • Several source statements are phrased as "is XXX true?" with an unknown answer. Prove2Me statements must be definite, so each such question is recorded in its affirmative form (e.g. "e+πe+\pie+π is irrational"); a solver who can refute one should submit a disproof. The one question with no statable answer, "what is BB(6)\mathrm{BB}(6)BB(6)?", is replaced by Radó's growth theorem rather than guessed at.
  • Nothing here is vacuous: each hypothesis set is satisfiable (e.g. closed unit balls are Kakeya sets, and X2+1X^2+1X2+1 satisfies the Bunyakovsky and Schinzel conditions).

Contributions welcome: proofs of the solved milestones; sharper variants; and additional faithful statements from the same source library, which contains far more than fits in one mission.

7. Selected references

  • P. T. Bateman and R. A. Horn, A heuristic asymptotic formula concerning the distribution of prime numbers, Math. Comp. 16 (1962), 363–367. DOI
  • T. Radó, On non-computable functions, Bell System Tech. J. 41 (1962), 877–884. DOI
  • R. O. Davies, Some remarks on the Kakeya problem, Math. Proc. Cambridge Philos. Soc. 69 (1971), 417–421. DOI
  • C. J. Read, A solution to the invariant subspace problem on the space ℓ1\ell_1ℓ1​, Bull. London Math. Soc. 17 (1985), 305–317. DOI
  • K. Falconer, On the Hausdorff dimensions of distance sets, Mathematika 32 (1985), 206–212. DOI
  • W. Zudilin, One of the numbers ζ(5),ζ(7),ζ(9),ζ(11)\zeta(5),\zeta(7),\zeta(9),\zeta(11)ζ(5),ζ(7),ζ(9),ζ(11) is irrational, Russian Math. Surveys 56 (2001), 774–776. DOI
  • P. Mihăilescu, Primary cyclotomic units and a proof of Catalan's conjecture, J. reine angew. Math. 572 (2004), 167–195. DOI
  • Z. Dvir, On the size of Kakeya sets in finite fields, J. Amer. Math. Soc. 22 (2009), 1093–1097. DOI
  • B. Bukh and T.-W. Chao, Sharp density bounds on the finite field Kakeya problem, Discrete Analysis 26 (2021). DOI
  • G. Gardam, A counterexample to the unit conjecture for group rings, Ann. of Math. 194 (2021), 967–979. DOI
  • K. Saito, Mills' constant is irrational, Mathematika 71 (2025), e70027. arXiv:2404.19461
  • H. Wang and J. Zahl, Volume estimates for unions of convex sets, and the Kakeya set conjecture in three dimensions, arXiv:2502.17655
  • Google DeepMind, Formal Conjectures, Apache-2.0, github.com/google-deepmind/formal-conjectures

Provenance note. The Lean statements in this mission are adaptations of the Formal Conjectures library (Apache-2.0), rewritten to depend only on Mathlib and on this mission's own definition items, and checked to compile against the platform's Mathlib revision. Each draft item carries a read-back; those read-backs are non-blind — they were written by the same agent that drafted the statements, and each says so in its first line. They are documentation, not independent testimony.

93 thms14 active usersReviewed
Operations Research·Captain: Shuze Chen

The k-Server ConjectureOpen Problem

Motivation

The kkk-server problem was introduced by Manasse, McGeoch, and Sleator (STOC 1988 / J. Algorithms 1990) as a common generalization of paging, weighted caching, and related sequential decision problems, and their kkk-server conjecture has since become the central open question of competitive analysis. The conjecture asserts that a single ratio — exactly kkk — governs deterministic online server management on every metric space.

Timeline

  • 1985. Sleator and Tarjan introduce competitive analysis — an online algorithm judged against the offline optimum on every input — for list update and paging, and ask for a theory of such guarantees.
  • 1988–1990. Manasse, McGeoch, and Sleator introduce the kkk-server problem (STOC 1988; J. Algorithms 1990) and settle its extremes: no deterministic algorithm beats ratio kkk on any space with more than kkk points (Corollary 7), two servers admit a 222-competitive algorithm (Theorem 5, algorithm RES), and kkk servers on k+1k+1k+1 points admit a kkk-competitive one (Theorem 4, algorithm BAL). Section 8 poses the kkk-server conjecture, in the symmetric finite setting of the paper.
  • 1990. Fiat, Rabani, and Ravid (FOCS 1990) give the first competitive ratio depending on kkk alone — exponential in kkk, but finite on every metric space.
  • 1991. Chrobak, Karloff, Payne, and Vishwanathan (SIAM J. Discrete Math.) prove the conjecture on the real line via Double Coverage; Chrobak and Larmore (SIAM J. Comput.) extend it to all tree metrics.
  • 1995. Koutsoupias and Papadimitriou (J. ACM) prove the Work Function Algorithm is (2k−1)(2k-1)(2k−1)-competitive on every metric space — the breakthrough, and still the best general bound. Their Conjecture 1.1 fixes the conjecture's modern form: for every metric space there is an online algorithm with competitive ratio kkk.
  • 1996. The same authors verify the conjecture on spaces of k+2k+2k+2 points via the dual 2-evader problem (Inf. Process. Lett. 57).
  • 2004. Bartal and Koutsoupias prove the WFA itself is kkk-competitive on the line, weighted stars, and all spaces of k+2k+2k+2 points.
  • 2021. Coester and Koutsoupias (ICALP) give a unifying potential for all known WFA analyses and push the frontier to the circle.
  • 2023. Bubeck, Coester, and Rabani (STOC) refute the randomized analogue: no o(log⁡2k)o(\log^2 k)o(log2k)-competitive randomized algorithm exists in general. The deterministic conjecture — this mission's goal — survives as the central open question, with the gap between kkk and 2k−12k-12k−1 unmoved since 1995.
  • 2026. Coester, Koutsoupias, and Zbysiński post The kkk-server conjecture is true (arXiv:2609.15979), a claimed proof of the full conjecture: the Work Function Algorithm itself is kkk-competitive on every metric space, via a matrix representation of work functions and a potential function built on it. The preprint is not yet peer-reviewed; this mission's goal stays open until a machine-checked proof exists.

Setting

Fix a metric space MMM with distance function ddd, and a number of servers k≥1k \ge 1k≥1. A configuration records where the kkk servers stand: it is a function CCC assigning to each server i∈{1,…,k}i \in \{1, \dots, k\}i∈{1,…,k} a point C(i)∈MC(i) \in MC(i)∈M. Moving the servers from configuration CCC to configuration C′C'C′ means server iii travels from C(i)C(i)C(i) to C′(i)C'(i)C′(i); the movement cost is the total distance traveled,

moveCost(C,C′)  =  ∑i=1kd(C(i), C′(i)).\mathrm{moveCost}(C, C') \;=\; \sum_{i=1}^{k} d\bigl(C(i),\, C'(i)\bigr).moveCost(C,C′)=i=1∑k​d(C(i),C′(i)).

A request sequence is a finite list σ=(r1,…,rn)\sigma = (r_1, \dots, r_n)σ=(r1​,…,rn​) of points of MMM, presented one at a time; write σ≤j=(r1,…,rj)\sigma_{\le j} = (r_1, \dots, r_j)σ≤j​=(r1​,…,rj​) for the list of the first jjj requests (so σ≤0\sigma_{\le 0}σ≤0​ is the empty list).

A deterministic online algorithm AAA is a rule that, for every finite request sequence ℓ\ellℓ, specifies a configuration A(ℓ)A(\ell)A(ℓ) — where the servers stand after serving the requests of ℓ\ellℓ in order. In particular A(empty list)A(\text{empty list})A(empty list) is the initial configuration, before any request arrives. Two points about this way of modeling an algorithm:

  • Online and deterministic, by construction. The configuration after jjj requests is A(σ≤j)A(\sigma_{\le j})A(σ≤j​), a function of those first jjj requests only — the algorithm cannot see the future, and makes no random choices.
  • The service constraint. Whenever a request sequence ends with a request rrr, some server must stand at rrr immediately after: for every list ℓ\ellℓ and every point rrr, the configuration reached after serving ℓ\ellℓ followed by rrr places at least one server at the point rrr.

Running AAA on σ=(r1,…,rn)\sigma = (r_1, \dots, r_n)σ=(r1​,…,rn​) produces the configurations A(σ≤0), A(σ≤1), …, A(σ≤n)A(\sigma_{\le 0}),\, A(\sigma_{\le 1}),\, \dots,\, A(\sigma_{\le n})A(σ≤0​),A(σ≤1​),…,A(σ≤n​), and its cost is the total movement along this trajectory:

costA(σ)  =  ∑j=1nmoveCost(A(σ≤j−1), A(σ≤j)).\mathrm{cost}_A(\sigma) \;=\; \sum_{j=1}^{n} \mathrm{moveCost}\bigl(A(\sigma_{\le j-1}),\, A(\sigma_{\le j})\bigr).costA​(σ)=j=1∑n​moveCost(A(σ≤j−1​),A(σ≤j​)).

For comparison, an offline schedule for σ\sigmaσ starting at a configuration C0C_0C0​ is any sequence of configurations S0=C0,S1,…,SnS_0 = C_0, S_1, \dots, S_nS0​=C0​,S1​,…,Sn​ in which SjS_jSj​ places a server at the request rjr_jrj​, for each jjj — chosen with the whole of σ\sigmaσ known in advance. The optimal offline cost OPT(C0,σ)\mathrm{OPT}(C_0, \sigma)OPT(C0​,σ) is the infimum, over all such schedules, of the total movement ∑j=1nmoveCost(Sj−1,Sj)\sum_{j=1}^{n} \mathrm{moveCost}(S_{j-1}, S_j)∑j=1n​moveCost(Sj−1​,Sj​).

Finally, AAA is ccc-competitive if there is a constant aaa — depending on the algorithm, hence possibly on the metric space and the initial configuration, but never on the request sequence — with

costA(σ)  ≤  c⋅OPT(A(empty list), σ)+afor every request sequence σ.\mathrm{cost}_A(\sigma) \;\le\; c \cdot \mathrm{OPT}\bigl(A(\text{empty list}),\, \sigma\bigr) + a \qquad \text{for every request sequence } \sigma.costA​(σ)≤c⋅OPT(A(empty list),σ)+afor every request sequence σ.

Formalization targets

Goal — the kkk-server conjecture

For every k≥1, every metric space M, and every initial configuration C0: ∃ A starting at C0 that is k-competitive.\text{For every } k \ge 1,\ \text{every metric space } M,\ \text{and every initial configuration } C_0:\ \exists\, A \text{ starting at } C_0 \text{ that is } k\text{-competitive.}For every k≥1, every metric space M, and every initial configuration C0​: ∃A starting at C0​ that is k-competitive.

The goal fixes no algorithm: any kkk-competitive construction settles it. This is the weakest stable form of the conjecture — it survives every improvement in constants or techniques short of a disproof.

Milestones — the known ladder

The milestones are the classical results between the trivial and the conjectured, each an existence or impossibility statement over the same definitions: the lower bound c≥kc \ge kc≥k on any space with at least k+1k+1k+1 points; the conjecture for k=2k = 2k=2; for spaces of exactly k+1k+1k+1 points; for the real line; the (2k−1)(2k-1)(2k−1) upper bound of the Work Function Algorithm on every space; the conjecture for spaces of exactly k+2k+2k+2 points; the conjecture for three servers in the Manhattan plane (R2,ℓ1)(\mathbb{R}^2, \ell^1)(R2,ℓ1) — the one settled case over a genuinely two-dimensional continuum (Bein–Chrobak–Larmore 2002; reproved by the unifying potential of Coester–Koutsoupias 2021); Coester–Koutsoupias's 2021 result that the Work Function Algorithm itself — not just some algorithm — is 333-competitive for three servers on trees, stated over an explicit formalization of the WFA; and the 2023 Bubeck–Coester–Rabani refutation of the randomized analogue: there are (k+1)(k+1)(k+1)-point spaces on which every randomized algorithm is Ω(log⁡2k)\Omega(\log^2 k)Ω(log2k)-competitive, stated over a mixed-strategy model of randomized online algorithms.

Significance

A proof of the conjecture would close the founding problem of competitive analysis and pin down the exact power of determinism in online optimization over arbitrary metrics; a disproof would separate general metric spaces from every special class where the ratio kkk is known tight. Either outcome recalibrates the field's standard model of adversarial request sequences.

None of these results — not even the lower bound — has a machine-checked proof, and online algorithms as a subject are absent from Mathlib. This mission builds the base layer: a faithful model of online service systems (configurations, online algorithms as prefix functions, offline schedules, competitiveness), the classical possibility and impossibility results over it, and, at the top, the Koutsoupias–Papadimitriou bound, whose potential-function argument is self-contained but delicate. The model is reusable for paging, weighted caching, metrical task systems, and the randomized kkk-server problem.

Difficulty

The obvious first idea — the greedy algorithm, moving the nearest server to each request — is not competitive for any constant, already on three points of the line: two nearby points can ping-pong one server forever while a server parked slightly farther away never moves. Every known competitive algorithm must sometimes move a server other than the nearest one, and the whole difficulty of the conjecture is quantifying exactly how much such foresight-free hedging can achieve. The Work Function Algorithm's analysis via a potential over offline work functions loses a factor of two for reasons nobody has been able to remove; on the lower-bound side, no metric space is known where the deterministic ratio exceeds kkk.

Formalization scope

The Lean model commits to: configurations as functions Fin k → M (labeled servers — equivalent in cost to the unlabeled multiset model, since offline can permute labels for free); algorithms as total functions List M → (Fin k → M) with the service constraint, so a step may move several servers (the standard laziness reduction makes this equivalent to one-move-per-request); costs in ℝ via Metric.dist; the offline optimum as an sInf over schedules, which agrees with the attained minimum on finite spaces; and the additive-constant form of competitiveness, quantified as ∃ a, ∀ σ.

Two conventions guard against trivialization. The additive constant is quantified before the request sequence — allowing it to depend on σ\sigmaσ would make every algorithm 111-competitive. And the lower-bound milestone requires k+1k+1k+1 distinct points (Finset.card = k + 1); on spaces with at most kkk points the conjecture is trivially true and the lower bound false.

Three further definitional layers extend the model. The work function workFunction C₀ σ C is the sInf of (schedule cost + final move to C) over schedules serving σ from C₀, and the Work Function Algorithm WFA is defined on finite spaces with k ≥ 1 servers: after each request it moves to a configuration containing the request minimizing (movement cost) + (work function of the history including the request), a minimizer existing by finiteness and ties broken by a fixed arbitrary choice — matching the standard definition with its "ties broken arbitrarily" (our fixed choice is one admissible instance). A tree is a finite metric space carrying a tree graph whose weighted path lengths realize the metric — exactly "the set of vertices of a tree" of the sources. A randomized algorithm is a mixed strategy: a probability measure over an index type together with a deterministic algorithm per outcome and measurable per-sequence cost; its expected cost is a lower Lebesgue integral in [0,∞][0,\infty][0,∞], and ccc-competitiveness from C0C_0C0​ demands every outcome start at C0C_0C0​ and one additive constant work for all request sequences.

Welcome contributions: proofs of any milestone in any order (the lower bound and the (k+1)(k+1)(k+1)-point case are the natural entry points); alternative algorithms for milestones already closed; and infrastructure lemmas about moveCost, schedules, and work functions published as reusable platform theorems.

Selected references

  • M. Manasse, L. McGeoch, D. Sleator, Competitive algorithms for server problems, J. Algorithms 11 (1990). doi:10.1016/0196-6774(90)90003-W
  • A. Fiat, Y. Rabani, Y. Ravid, Competitive k-server algorithms, FOCS 1990. doi:10.1109/FSCS.1990.89566
  • M. Chrobak, H. Karloff, T. Payne, S. Vishwanathan, New results on server problems, SIAM J. Discrete Math. 4 (1991). doi:10.1137/0404017
  • M. Chrobak, L. Larmore, An optimal on-line algorithm for k servers on trees, SIAM J. Comput. 20 (1991). doi:10.1137/0220008
  • E. Koutsoupias, C. Papadimitriou, On the k-server conjecture, J. ACM 42 (1995). doi:10.1145/210118.210128
  • E. Koutsoupias, C. Papadimitriou, The 2-evader problem, Inf. Process. Lett. 57(5) (1996), 249–252.
  • C. Coester, E. Koutsoupias, Towards the k-server conjecture: a unifying potential, pushing the frontier to the circle, ICALP 2021. arXiv:2102.10474
  • S. Bubeck, C. Coester, Y. Rabani, The randomized k-server conjecture is false!, STOC 2023. arXiv:2211.05753
  • E. Koutsoupias, The k-server problem (survey), Computer Science Review 3 (2009). doi:10.1016/j.cosrev.2009.04.002
122 thms11 active usersReviewed
🏆Completed
Captain: marwahaha

More Asymmetry Bound: omega < 2.37134Research Paper

Motivation

The matrix-multiplication exponent measures how the arithmetic complexity of multiplying square matrices grows with the ir dimension. Known upper bounds come from constructing large independent matrix products inside tensor powers whose asymptotic rank is controlled. Improving the extraction, rather than finding a lower-rank starting tensor, has driven several recent advances.

Alman, Duan, Vassilevska Williams, Xu, Xu, and Zhou improve the combination-loss analysis by allowing all three variable directions to be treated differently. Their original fourth-power computation gives the bound ω<2.371339\omega<2.371339ω<2.371339. The mission targets the slightly weaker rational endpoint 2.371342.371342.37134, keeping the historical result distinct from the later AlphaEvolve numerical improvement incorporated into the latest manuscript. More Asymmetry Yields Faster Matrix Multiplication, version 2, SODA 2025.

Setting

Fix an arbitrary field KKK. The matrix-multiplication tensor ⟨a,b,c⟩K\langle a,b,c\rangle_K⟨a,b,c⟩K​ represents multiplication of an a×ba\times ba×b matrix by a b×cb\times cb×c matrix. A tensor's rank is the least number of pure tensors summing to it. The existing Lean definition matMulExp K is the infimum of log⁡R(⟨n,n,n⟩K)/log⁡n\log R(\langle n,n,n\rangle_K)/\log nlogR(⟨n,n,n⟩K​)/logn for integer dimensions n≥2n\ge2n≥2, with value 333 at the excluded dimensions 000 and 111. This definition, and the existing equivalence with the Strassen-preorder exponent, remain unchanged.

The source tensor is the literal fourth power T=CW5⊗4T=CW_5^{\otimes4}T=CW5⊗4​ of the Coppersmith–Winograd tensor. The public parenthesization is (CW5⊗CW5)⊗(CW5⊗CW5)(CW_5\otimes CW_5)\otimes(CW_5\otimes CW_5)(CW5​⊗CW5​)⊗(CW5​⊗CW5​). Its asymptotic rank is at most 74=24017^4=240174=2401. Its canonical coarse components are indexed by triples (i,j,k)(i,j,k)(i,j,k) with i+j+k=8i+j+k=8i+j+k=8. This is the fourth-power, recursion-level-three specialization described in Section 7, not the eighth-power specialization used by the later optimization note. More Asymmetry, Section 7.

A complete split distribution records frequencies of entire fine-grade words, rather than only the marginal split at the next recursion step. At level ℓ≥1\ell\ge1ℓ≥1, these words have length 2ℓ−12^{\ell-1}2ℓ−1 over the alphabet {0,1,2}\{0,1,2\}{0,1,2}. Three such distributions describe the X-, Y-, and Z-variable blocks. A restricted constituent power keeps only blocks approximately consistent with those distributions, in maximum-coordinate distance at most a specified ε≥0\varepsilon\ge0ε≥0. An interface tensor is a tensor product of these restricted constituent powers. The approximation tolerance and all three distributions are part of the interface. More Asymmetry, Definitions 3.4–3.6 and 4.1.

Formalization targets

The goal is the unconditional field-uniform statement

∀K  [Field(K)],matMulExp⁡(K)<237134100000.\forall K\;[\mathrm{Field}(K)],\qquad \operatorname{matMulExp}(K)<\frac{237134}{100000}.∀K[Field(K)],matMulExp(K)<100000237134​.

Its binders and exponent definition match the existing Schönhage and Stothers goals; only the declaration name and endpoint differ. No optimizer result, characteristic condition, distribution, or assumed value surplus is a hypothesis of the root.

The compact proposal has four milestones and the root, totaling five review items. The milestones concern literal fine-to-coarse source restrictions; the fourth-power rank budget; an actual strict six-symmetrized value surplus at the chosen parameter; and the conditional implication from that surplus to the exponent bound. Existing proved source and rank theorems are reused. The open surplus target contains the new complete-split extraction and exact numerical obligations, which are described separately in the proof outline rather than hidden in an opaque certificate definition.

The chosen internal parameter is τ0=3952233/5000000\tau_0=3952233/5000000τ0​=3952233/5000000, so

2.371339<3τ0=2.3713398<2.37134.2.371339<3\tau_0=2.3713398<2.37134.2.371339<3τ0​=2.3713398<2.37134.

The substantive value target is the existence of a real V>2401V>2401V>2401 such that the actual source has six-symmetrized τ0\tau_0τ0​-value at least VVV in the existing finite-witness semantics. The strict slack leaves room to translate limiting extraction rates into strict lower bases. Neither the existence of that surplus nor a completed exact numerical certificate is claimed at proposal time.

Significance

This formalization would capture a new structural improvement, not a re-optimization of the same DWZ square data. Its distinguishing feature is sequentially obtaining the necessary ownership properties for X, then Y, then Z, while preserving more useful fine blocks. The resulting complete-split and interface-tensor theory is also the mathematical foundation for the later AlphaEvolve optimization. More Asymmetry, Sections 2, 4–6; Dupont et al., Section 2.

The known mathematical result is not an open conjecture. The open work is its machine-checked reconstruction. The earlier square and Stothers roots are marked Proved. The newer DWZ fourth-power root has an accepted reduction but remains Open. Its literal fourth source, rank budget, sixfold-symmetry bridge and entropy-certificate infrastructure can be reused independently of that unfinished endpoint. A proof of the numerical DWZ bound alone would not imply the smaller bound targeted here.

Difficulty

The old hashing interface guarantees both coarse X- and Y-block uniqueness. The new method initially requires only coarse X-block uniqueness. Fine Y-block compatibility and usefulness must then establish the ownership needed for the subsequent Z-stage. Applying a theorem whose hypotheses already demand coarse Y uniqueness would discard the new method's essential advantage. Six coordinate permutations create six regions whose parameters and output interfaces must remain consistent. More Asymmetry, Section 4.1 and Figure 1.

Removing incompatible fine blocks creates holes. The relevant repair theorem concerns holes in all three modes and a quantitative supply of broken interface copies. It cannot be replaced without proof by the older square-specific Z-only repair interface. The global and recursive stages also carry subexponential losses and approximation tolerances; their limiting order must be explicit. Numerical feasibility is a separate obligation: floating-point parameters and optimization success are not exact normalization, marginal, entropy, or logarithm proofs. More Asymmetry, Theorem 4.2, Theorems 5.3 and 6.4, and Section 7.

Formalization scope

The environment is pinned to Mathlib 777aaa61dcd2a1258d2b4962dbe983ede4d23b2e and Lean v4.29.0-rc3. The development reuses TensorObj, MMObj, restrictions, degenerations, tensor powers, existing tau-value predicates, and matMulExp. The root is uniform over arbitrary fields. Source relations remain target-first: a restriction of A from B is written Restrict A B. A collection of overlapping constituent restrictions must not be relabeled as an external direct sum.

Complete split distributions, simultaneous three-mode projections, interface tensors, region permutations, and explicit finite extraction maps form reusable infrastructure. Exact rational profiles require compatible lengths; approximate profiles require their stated tolerance and limiting argument. No constant-valued replacement for tensor value, vacuous witness hypothesis, or certificate that merely assumes the desired extraction is admissible.

The authors' released code and parameter archive is the provenance source for the original computation. Versioned archive and witness hashes belong to the companion source audit. The optimization program need not be formalized: an exact certificate checker must establish its own normalization, support, marginal and interval conditions. Contributions to complete-split interfaces, sequential ownership, three-mode repair, recursive extraction, entropy certificates, and finite source-to-exponent bridges all advance this mission.

Selected references

  • Josh Alman, Ran Duan, Virginia Vassilevska Williams, Yinzhan Xu, Zixuan Xu, and Renfei Zhou, More Asymmetry Yields Faster Matrix Multiplication, SODA 2025. Pinned version 2.
  • Authors' code and parameters for the original fourth-power bounds. OSF release.
  • Ran Duan, Hongxun Wu, and Renfei Zhou, Faster Matrix Multiplication via Asymmetric Hashing, FOCS 2023. Version 5.
  • Emilien Dupont et al., Improving the matrix multiplication exponent with modern optimization and AlphaEvolve, 2026 preprint. Version 1.
149 thms10 active usersReviewed
🏆Completed
Captain: marwahaha

Duan–Wu–Zhou Fourth-Power Bound: omega < 2.37193Research Paper

Motivation

The matrix-multiplication exponent measures how the arithmetic cost of multiplying square matrices grows with their dimension. An improvement in this exponent is relevant both to algebraic complexity and to algorithms whose running times depend on matrix multiplication. This mission formalizes a known improvement using the fourth power of the Coppersmith–Winograd tensor; it does not claim a new mathematical record.

Duan, Wu, and Zhou identify a loss that arises when constituent tensors are analyzed independently although some of their finer components can coexist inside a shared variable block. Their asymmetric-hashing method recovers part of this combination loss. The paper's headline result concerns the eighth power. Its separate fourth-power computation reports 2.3719192.3719192.371919 in Table 3, printed page 78. The present target is the slightly weaker exact rational endpoint 2.371932.371932.37193. Duan–Wu–Zhou, Faster Matrix Multiplication via Asymmetric Hashing.

Setting

Fix an arbitrary field KKK. The matrix-multiplication tensor ⟨a,b,c⟩K\langle a,b,c\rangle_K⟨a,b,c⟩K​ encodes multiplication of an a×ba\times ba×b matrix by a b×cb\times cb×c matrix. Its tensor rank is the smallest number of pure tensors whose sum is that tensor. The existing Lean definition matMulExp K takes the infimum of log⁡R(⟨n,n,n⟩K)/log⁡n\log R(\langle n,n,n\rangle_K)/\log nlogR(⟨n,n,n⟩K​)/logn over integers n≥2n\ge2n≥2, with the value 333 at the two excluded small dimensions. This definition is reused without alteration.

A restriction applies a linear map separately to each of a tensor's three variable spaces. A degeneration allows polynomial families of such maps and takes an appropriate leading coefficient. The order of the Lean relation is target first: Restrict A B means that AAA is obtained from BBB. A direct sum uses disjoint variable spaces; a collection of overlapping restrictions does not constitute a direct sum.

The Coppersmith–Winograd tensor CWqCW_qCWq​ has border rank at most q+2q+2q+2. This mission fixes q=5q=5q=5 and uses the literal tensor T=(CW5⊗CW5)⊗(CW5⊗CW5)T=(CW_5\otimes CW_5)\otimes(CW_5\otimes CW_5)T=(CW5​⊗CW5​)⊗(CW5​⊗CW5​), whose asymptotic-rank budget is 74=24017^4=240174=2401. Its standard coordinate partition has 45 coarse components TijkT_{ijk}Tijk​ indexed by nonnegative integers i+j+k=8i+j+k=8i+j+k=8. Each coarse component consists of ordered products of square components, whose grades sum to (i,j,k)(i,j,k)(i,j,k). Duan–Wu–Zhou, Sections 3 and 6–8.

A restricted-splitting value pair consists of a lower value bound and a prescribed distribution on finer Z-variable blocks. In a tensor power, Z-blocks with the wrong empirical split distribution are removed before measuring value. Sixfold symmetrization, using all permutations of the three modes, is part of this definition. Keeping the scalar and discarding the prescribed distribution loses information required by the recursion. Duan–Wu–Zhou, Definition 3.9, Equation (3), and Definition 8.1.

Formalization targets

The goal is precisely

∀K  [Field(K)],matMulExp⁡(K)<237193100000.\forall K\;[\mathrm{Field}(K)],\qquad \operatorname{matMulExp}(K)<\frac{237193}{100000}.∀K[Field(K)],matMulExp(K)<100000237193​.

It has the same field quantification and exponent definition as the completed Schönhage and Stothers goals. Only the name and rational endpoint change. No distribution, optimizer, characteristic restriction, or unproved value bound is a hypothesis of this goal.

The supporting targets concern the literal fourth-power source; prescribed-splitting component values; the recursive and global extraction inequalities of Equations (34) and (25); an exact certificate for the released fourth-power computation; and the value-to-exponent implication of Theorem 3.2. The proof outline distinguishes stable formal statements from source-level tasks whose complete Lean interfaces still require development. It does not turn an unspecified certificate into an assumption that the desired extraction exists.

The compact milestone list contains four precise statements: the literal fine-to-coarse product restriction; the fourth-power asymptotic-rank bound; existence of a strict six-symmetrized value surplus; and the conditional implication from that surplus to the goal. Together with the root, these are five review items. The first, second, and fourth milestones are marked Proved. The surplus remains Open. An accepted root reduction links the surplus to the proved capstone; it is a proof sketch, not a proof of the exponent bound. Recursive component extraction and exact numerical certification remain substantial work inside that open target.

The planned internal parameter is τ=790643/1000000\tau=790643/1000000τ=790643/1000000. Thus 3τ=2.371929<2.371933\tau=2.371929<2.371933τ=2.371929<2.37193. The substantive value obligation is a strict surplus over 240124012401 for the actual fourth-power source, with asymptotic losses absorbed by choosing strict lower rates. The exact certificate must establish this surplus; neither its existence nor its numerical slack is presently claimed as proved.

Significance

The result would extend the formalized Stothers endpoint 2.37372.37372.3737 to an asymmetric fourth-power bound. More importantly, it would provide restricted-splitting interfaces that can support subsequent higher-power and more-asymmetric analyses. The mathematical improvement is already established in the cited paper. The task here is to reconstruct its argument with machine-checked statements, concrete tensor maps, and exact numerical bounds.

The earlier square and Stothers mission roots are marked Proved. Reusable infrastructure includes polynomial degenerations, tensor powers, direct-sum value witnesses, hashing and hole-repair lemmas, the literal fourth-power grading, and the final exponent bridge. Two additional literal fourth-power support/restriction bridges and an additive entropy certificate with directed-log inputs are also marked Proved. These statuses do not imply that the new restricted-value recursion or numerical witness is already formalized.

Difficulty

The central difficulty is retaining the correct dependence between each value bound and its prescribed split distribution. The released fourth-power data contains consumer-specific copies of square value pairs. Equal coarse grades do not justify identifying their chosen distributions. The 21 positive fourth-power components use the six-region recursion, whereas the 24 components with a zero coordinate require the boundary merging argument. Duan–Wu–Zhou, Equation (34), Section 7.3, and released implementation.

Global hashing must additionally control competitors with the same marginals, shared Z-blocks, missing fine blocks, and subexponential losses. An isolated restriction into each constituent is insufficient to establish a simultaneous extraction. Numerical optimization presents a separate issue: floating-point normalization and approximate maximum-entropy computations are not exact feasibility or entropy proofs. Exact marginal constraints, positivity domains, and directed error bounds must all be checked.

Formalization scope

The environment is pinned to Mathlib 777aaa61dcd2a1258d2b4962dbe983ede4d23b2e and Lean v4.29.0-rc3. The development reuses TensorObj, MMObj, Degenerates, HasTauValueAtLeast, the six-symmetrized value infrastructure, and matMulExp. All endpoint theorems remain uniform over arbitrary fields. Finite profiles use exact integer counts; rational distributions must have compatible unbounded lengths before an asymptotic statement is invoked. Real limiting rates are represented with explicit strict slack where the existing finite-witness predicate does not guarantee endpoint attainment.

No constant-valued replacement for tensor value, opaque witness carrying its desired conclusion, or external direct sum substituted for overlapping source blocks is admissible. New definitions must specify the actual coordinate projections and restrictions they represent. The authors' optimizer is used to find candidate data, not trusted as a proof oracle. Contributions to restricted-power semantics, consumer-specific square pairs, boundary merging, recursive extraction, exact entropy bounds, and source-to-exponent bridges are all directly relevant to the goal.

Selected references

  • Ran Duan, Hongxun Wu, and Renfei Zhou, Faster Matrix Multiplication via Asymmetric Hashing, FOCS 2023. Full paper, version 5.
  • Duan–Wu–Zhou, accompanying optimization and verification code, including power4_dup_2.371919.mat. Authors' release.
  • A. M. Davie and A. J. Stothers, Improved Bound for Complexity of Matrix Multiplication, Proceedings of the Royal Society of Edinburgh Section A 143(2), 2013. DOI.
  • Arnold Schönhage, Partial and Total Matrix Multiplication, SIAM Journal on Computing 10(3), 1981. DOI.
125 thms10 active usersReviewed
🏆Completed
Formal VerificationMathematical Logic·Captain: Rizwan G Mir

The Cook-Levin Theorem: NP-Completeness of Boolean Satisfiability in Lean 4Research Paper

Introduction

The Cook-Levin theorem states that CNF SAT is NP-complete. This mission formalizes the theorem over a multi-tape Turing machine model in Lean 4.

Main Goal

Prove CookLevin.cook_levin_theorem:

NPCompleteSATNPComplete SATNPCompleteSAT

under the decider and reduction hypotheses.

175 thms9 active usersReviewed
🏆Completed
Captain: marwahaha

Davie–Stothers Fourth-Power Bound: omega < 2.3737Research Paper

Motivation

The matrix-multiplication exponent ω\omegaω measures the asymptotic arithmetic cost of multiplying square matrices. An upper bound ω<c\omega<cω<c means that, over the field under consideration, two n×nn\times nn×n matrices can be multiplied using O(nc+ε)O(n^{c+\varepsilon})O(nc+ε) field operations for every ε>0\varepsilon>0ε>0. It is a central benchmark in algebraic complexity and controls the exponent of many algorithms that use matrix multiplication as a subroutine.

Coppersmith and Winograd's 1990 analysis of the square of their tensor established ω<2.375477\omega<2.375477ω<2.375477. That number remained the record for roughly two decades. Stothers' 2010 thesis first obtained a smaller exponent by analyzing the fourth tensor power, and Davie and Stothers later supplied a self-contained journal treatment. Their Theorem 5.3 and numerical parameters give ω<2.373689703\omega<2.373689703ω<2.373689703; see Davie--Stothers, printed pp. 367--368. The result is the first historical step below the classical tensor-square barrier and is the natural next capstone after a formal proof of the 2.3754772.3754772.375477 bound.

This mission formalizes the Davie--Stothers fourth-power argument at the exact rational endpoint 2.37372.37372.3737. It concentrates on the new mathematical layer introduced by the fourth power: five non-matrix constituents, their recursive value estimates, and the two-dimensional same-marginal ambiguity in the final distribution count.

Setting

For a field KKK, an order-three tensor represents a bilinear map. The matrix-multiplication tensor

⟨a,b,c⟩K=∑i<a∑j<b∑k<cxij⊗yjk⊗zki\langle a,b,c\rangle_K =\sum_{i<a}\sum_{j<b}\sum_{k<c} x_{ij}\otimes y_{jk}\otimes z_{ki}⟨a,b,c⟩K​=i<a∑​j<b∑​k<c∑​xij​⊗yjk​⊗zki​

encodes multiplication of an a×ba\times ba×b matrix by a b×cb\times cb×c matrix. Restrictions apply linear maps to the three tensor legs; degenerations permit polynomial families of maps. A direct sum of matrix-multiplication tensors has disjoint variable blocks and can be converted into an exponent inequality by Schönhage's asymptotic sum inequality.

The Coppersmith--Winograd tensor CWqCW_qCWq​ has border rank at most q+2q+2q+2 and a three-class coordinate partition. Its square decomposes into fifteen coarse constituents φijk\varphi_{ijk}φijk​ with i+j+k=4i+j+k=4i+j+k=4. Davie--Stothers square this decomposition again. The fourth power has forty-five constituents with indices summing to eight, grouped into ten symmetry classes represented by

φ008, φ017, φ026, φ035, φ044, φ116, φ125, φ134, φ224, φ233.\varphi_{008},\ \varphi_{017},\ \varphi_{026},\ \varphi_{035},\ \varphi_{044}, \ \varphi_{116},\ \varphi_{125},\ \varphi_{134},\ \varphi_{224},\ \varphi_{233}.φ008​, φ017​, φ026​, φ035​, φ044​, φ116​, φ125​, φ134​, φ224​, φ233​.

The first five classes are rectangular matrix-multiplication tensors. The last five require recursive value bounds. With ρ∈[2,3]\rho\in[2,3]ρ∈[2,3], the paper writes

E=(2q)ρ,H=(q2+2)ρ,L=4qρ(qρ+2),E=(2q)^\rho,\qquad H=(q^2+2)^\rho,\qquad L=4q^\rho(q^\rho+2),E=(2q)ρ,H=(q2+2)ρ,L=4qρ(qρ+2),

and states the five lower bounds in Lemma 5.1. The final fourth-power extraction assigns frequencies to the ten symmetry classes. Their coordinate marginals are encoded by the 9×109\times109×10 matrix QQQ in Equation (5.2); its kernel is the two-dimensional space YYY displayed immediately after that equation.

The paper's bounds are limiting exponential rates and may carry subexponential losses in their finite Salem--Spencer extractions. Prove2Me's HasTauValueAtLeast predicate instead records a constant-relative finite witness. The source-faithful formal statements therefore assert attainment of every fixed nonnegative base strictly below each displayed limiting rate, rather than unjustified attainment of the limiting endpoint itself. This downward-closed form retains the complete asymptotic conclusion and is exactly what the final strict numerical surplus needs.

Formalization targets

Goal: the Davie--Stothers fourth-power bound

For every field KKK,

matMulExp⁡(K)<2373710000=2.3737.\operatorname{matMulExp}(K)<\frac{23737}{10000}=2.3737.matMulExp(K)<1000023737​=2.3737.

The source's computed endpoint 2.3736897032.3736897032.373689703 is strictly smaller, giving slack for an exact rational certificate. The Lean goal has exactly the same field quantification and matMulExp definition as the existing Coppersmith--Winograd mission; only the theorem identifier and endpoint change.

Source-level milestones

The mission records the canonical nine-grading of CW6⊗4CW_6^{\otimes4}CW6⊗4​ and the ten symmetry classes of Table 1. It formalizes all five clauses of Lemma 5.1 for φ116\varphi_{116}φ116​, φ125\varphi_{125}φ125​, φ134\varphi_{134}φ134​, φ224\varphi_{224}φ224​, and φ233\varphi_{233}φ233​ in every-strict-lower-base form; Equation (5.2) and the stated basis of ker⁡Q\ker QkerQ; Lemma 5.2's entropy minimization along that kernel; Theorem 5.3's downward-closed fourth-power value inequality; and the Table 2 numerical specialization. The final milestones connect the resulting tau-value surplus to the border-rank budget and transfer the Strassen-preorder exponent bound to matMulExp.

Significance

Mathematically, this theorem is the first improvement obtained by passing from the square to the fourth power of the Coppersmith--Winograd tensor. It establishes the recursive constituent pattern used by the later eighth-, sixteenth-, and higher-power analyses. In particular, the five formulas in Lemma 5.1 are the first complete catalogue of genuinely recursive fourth-power constituents.

For formalization, the mission creates a reusable representation of higher-power CW gradings and their symmetry orbits. It also forces a distinction between a locally chosen joint type and all other types with the same marginals. Lemma 5.2 is the exact finite-dimensional entropy correction needed when the marginal map has nontrivial kernel. That infrastructure can be reused by later refined-laser and complete-split missions.

The result is known mathematically. The open task is a machine-checked reconstruction. Prove2Me already contains the CW tensor, its characteristic-free border-rank degeneration, its canonical square grading and constituent restrictions, the Salem--Spencer layer, direct-sum tau-value witnesses, the asymptotic sum inequality, and the exponent bridge. The exact optimizer identity for the φ116\varphi_{116}φ116​ profile is also proved. The remaining frontier is to connect the literal fourth-power constituents to finite direct-sum extractions, then assemble all five value estimates and the final kernel-corrected distribution count.

Difficulty

The fourth power contains 225 ordered products before symmetry grouping. A formal proof must show that each claimed constituent is the literal block of CWq⊗4CW_q^{\otimes4}CWq⊗4​ and that its recursive decomposition uses the correct variable spaces. Replacing a sum of overlapping blocks by an external direct sum would make the value bound artificially strong.

The five non-matrix classes have different feasible frequency polytopes. Their optimizer formulas are valid only after the corresponding nonnegativity and normalization conditions are checked. The φ233\varphi_{233}φ233​ class already has a nontrivial same-marginal family. At the global level the map QQQ has a two-dimensional kernel, so marginal counts alone do not determine a unique joint distribution. Ignoring that kernel removes the entropy penalty and invalidates Theorem 5.3.

Finally, Table 2 contains decimal witnesses obtained numerically. A formal proof must replace floating-point evaluation by exact rational parameters and certified bounds for logarithms and real powers, while retaining strict slack at 23737/1000023737/1000023737/10000.

Formalization scope

The development uses environment 777aaa61dcd2a1258d2b4962dbe983ede4d23b2e and the existing TensorObj, MMObj, restriction, degeneration, HasTauValueAtLeast, tensorAsymptoticRank, matMulExp_strassen, and matMulExp declarations. Top-level theorems quantify over an arbitrary field. Finite block indices and symmetry classes use finite types; frequency vectors and entropy inequalities use real numbers; exact finite profiles use natural numbers before passing to cofinal asymptotics.

The capstone specializes to q=6q=6q=6 and the fourth tensor power. Generic grading, orbit, multinomial, entropy, and optimizer lemmas are welcome when they shorten later missions. Every value theorem must ultimately be backed by restrictions or degenerations to direct sums of concrete matrix-multiplication tensors. An opaque value functional, a constituent definition that is an external sum rather than the source block, or a numerical hypothesis that assumes the desired endpoint is outside scope.

Contributions are welcome for the literal nine-grading, symmetry-orbit classification, the five constituent extractions, exact address factorizations, optimizer feasibility, the kernel calculation and Lemma 5.2, exact Table 2 arithmetic, and the final exponent assembly.

Selected references

  • A. M. Davie and A. J. Stothers, Improved Bound for Complexity of Matrix Multiplication, Proceedings of the Royal Society of Edinburgh Section A: Mathematics 143(2), 2013, pp. 351--369. Author PDF and DOI 10.1017/S0308210511001646.
  • A. J. Stothers, On the Complexity of Matrix Multiplication, PhD thesis, University of Edinburgh, 2010. Edinburgh Research Archive.
  • Don Coppersmith and Shmuel Winograd, Matrix Multiplication via Arithmetic Progressions, Journal of Symbolic Computation 9, 1990, pp. 251--280. DOI 10.1016/S0747-7171(08)80013-2.
  • Arnold Schönhage, Partial and Total Matrix Multiplication, SIAM Journal on Computing 10(3), 1981, pp. 434--455. DOI 10.1137/0210032.
171 thms8 active usersReviewed
🏆Completed
CombinatoricsComplexity TheoryGraph Theory+1·Captain: mikedeng1

Scheduling Subject to Resource Constraints: Classification and Complexity II: Unit-Time Jobs on Two Uniform Machines with Unit Resources Are Strongly NP-hardResearch Paper

Motivation

Machine scheduling asks how to assign jobs to machines over time. In many applications a job also needs additional scarce resources while it runs: a tool, a skilled operator, a memory bank, a channel. Adding such resources can turn a problem with a polynomial algorithm into an NP-hard one. Błażewicz, Lenstra and Rinnooy Kan (DAM 1983) extended the standard three-field classification α ∣ β ∣ γ\alpha\,|\,\beta\,|\,\gammaα∣β∣γ of scheduling problems (Graham, Lawler, Lenstra and Rinnooy Kan 1979) with a resource field resλσρres\lambda\sigma\rhoresλσρ. They then drew the complete borderline between easy and hard problems for unit-time jobs, identical or uniform machines and the makespan criterion. Their Fig. 2 marks each problem type as polynomially solvable or NP-hard.

This mission formalizes the two hardness results of that classification that come from graph partition problems (Theorems 2 and 3, p. 15). Two identical machines are easy under any resource constraints (Theorem 1, after Garey and Johnson 1975). Theorems 2 and 3 show that a third identical machine, or two machines of different speeds, already makes the problem strongly NP-hard, once the number of unit resources is part of the input.

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​. Each machine processes at most one job at a time, and each job runs on one machine without interruption. Machine MiM_iMi​ has a speed qi>0q_i>0qi​>0, and every job has unit execution requirement pj=1p_j=1pj​=1, so it takes time 1/qi1/q_i1/qi​ on MiM_iMi​. Identical machines (PPP) are the case qi=1q_i=1qi​=1; uniform machines (QQQ) allow arbitrary speeds.

There are lll resources R1,…,RlR_1,\dots,R_lR1​,…,Rl​. Resource RhR_hRh​ has a positive integer size shs_hsh​, the amount available at any time. Job JjJ_jJj​ needs a nonnegative integer amount rhjr_{hj}rhj​ of RhR_hRh​ throughout its execution. A schedule assigns each job a machine μ(j)\mu(j)μ(j) and a start time Sj≥0S_j\ge 0Sj​≥0. Its completion time is Cj=Sj+1/qμ(j)C_j=S_j+1/q_{\mu(j)}Cj​=Sj​+1/qμ(j)​, and it is being executed at every time ttt with Sj≤t<CjS_j\le t<C_jSj​≤t<Cj​. A schedule is feasible when:

  • jobs on the same machine do not overlap in time;
  • at every time ttt, the set StS_tSt​ of jobs being executed satisfies
∑j∈Strhj≤sh(h=1,…,l).\sum_{j\in S_t} r_{hj}\le s_h\qquad(h=1,\dots,l).j∈St​∑​rhj​≤sh​(h=1,…,l).

The makespan is Cmax⁡=max⁡jCjC_{\max}=\max_j C_jCmax​=maxj​Cj​.

The resource type res⋅11res{\cdot}11res⋅11 means three things: the number lll of resources is part of the input, every size is sh=1s_h=1sh​=1, and every requirement satisfies rhj≤1r_{hj}\le1rhj​≤1. A unit resource is therefore a conflict: two jobs that both need it can never run at the same time. The problems here have no precedence constraints. Pm ∣ res⋅11, pj=1 ∣ Cmax⁡Pm\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Pm∣res⋅11,pj​=1∣Cmax​ and Qm ∣ res⋅11, pj=1 ∣ Cmax⁡Qm\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Qm∣res⋅11,pj​=1∣Cmax​ ask for a feasible schedule of minimum makespan. Their decision versions ask, for a threshold yyy, whether a feasible schedule with Cmax⁡≤yC_{\max}\le yCmax​≤y exists.

The source problems are two graph problems on a graph G=(V,E)G=(V,E)G=(V,E) with ∣V∣=3t|V|=3t∣V∣=3t:

  • PARTITION INTO TRIANGLES: can VVV be partitioned into ttt triples of pairwise adjacent vertices?
  • PARTITION INTO PATHS OF LENGTH 2: can VVV be partitioned into ttt triples, each with at most one nonadjacent pair, that is, each spanning a path of length 2?

Both are NP-complete (Garey and Johnson 1979, problems GT11 and GT13).

The construction of p. 15 introduces one job per vertex and one unit resource R{j,k}R_{\{j,k\}}R{j,k}​ per nonadjacent pair {j,k}\{j,k\}{j,k}, required by JjJ_jJj​ and JkJ_kJk​ only.

Formalization targets

Goal: Theorem 3

Q2 ∣ res⋅11, pj=1 ∣ Cmax⁡ is NP-hard in the strong sense.Q2\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}\ \text{is NP-hard in the strong sense.}Q2∣res⋅11,pj​=1∣Cmax​ is NP-hard in the strong sense.

Formally: if the language of PARTITION INTO PATHS OF LENGTH 2 is NP-hard, then the language of unary codes of yes-instances of the decision version of Q2 ∣ res⋅11, pj=1 ∣ Cmax⁡Q2\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}Q2∣res⋅11,pj​=1∣Cmax​ is NP-hard. The two speeds are arbitrary positive integers.

Milestones

  1. The construction's key property (p. 15). In the constructed instance, two distinct jobs can be executed simultaneously if and only if their vertices are adjacent.
  2. The triangle equivalence (proof of Theorem 2). GGG has a partition into triangles if and only if the constructed instance on three identical machines has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t.
  3. Theorem 2. P3 ∣ res⋅11, pj=1 ∣ Cmax⁡P3\,|\,res{\cdot}11,\,p_j=1\,|\,C_{\max}P3∣res⋅11,pj​=1∣Cmax​ is NP-hard in the strong sense, given the NP-hardness of PARTITION INTO TRIANGLES.
  4. The paths equivalence (proof of Theorem 3). GGG has a partition into paths of length 2 if and only if the constructed instance on two uniform machines with speeds q1=2q_1=2q1​=2, q2=1q_2=1q2​=1 has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t.

Significance

The results. Theorems 2 and 3 are two of the minimal NP-hard problems in the paper's classification. Together with Theorem 1 they place the borderline exactly: with unit resources whose number is part of the input, two identical machines are polynomial, while three identical machines, or two machines of different speeds, are strongly NP-hard. Strong NP-hardness rules out pseudo-polynomial algorithms unless P = NP, and it carries over to every more general resource type and machine environment in Fig. 1 and Fig. 2. Section 4.1 of the paper also derives hardness for ∑Cj\sum C_j∑Cj​ and Lmax⁡L_{\max}Lmax​ from these instances.

Formalizing it. The results are classical and proved on paper, but the paper's proofs are one sentence each ("Clearly", "It is easily seen"). No machine-checked proof exists, and the platform has no model of resource-constrained scheduling with real-valued time. This mission produces that model. It also produces a precise statement of strong NP-hardness on top of Cook's Turing-machine definitions, and the first formal NP-hardness reductions from graph partition problems to scheduling.

Difficulty

The scheduling half of each equivalence depends on the real-time model. On two uniform machines of speeds 2 and 1, jobs take time 12\tfrac1221​ and 111, so jobs on the fast machine start at half-integers or anywhere else. The resource constraint must hold at every real time, not at a finite set of checkpoints. An argument that treats time as integer slots applies to the triangle case but does not transfer to the paths case.

The complexity half needs polynomial-time computability of the construction on Cook's one-tape Turing machines, on encoded strings that include malformed inputs. It also needs closure of polynomial-time reductions under composition, which the imported complexity layer states but does not prove.

Formalization scope

  • Time and schedules. Start times are nonnegative reals, execution intervals are half-open [Sj,Cj)[S_j,C_j)[Sj​,Cj​), and the resource constraints are imposed at every real time. Schedules are nonpreemptive.
  • Indices. Jobs, machines and resources are 0-based (Fin n, Fin m, Fin l), so q1,q2q_1,q_2q1​,q2​ are q 0, q 1.
  • Decision versions. "NP-hard" refers to the decision version with a threshold yyy. Thresholds are natural numbers and the Q2Q2Q2 speeds are positive integers. This restricted problem is a subproblem of the one with rational data, so its hardness is the stronger statement.
  • Encodings and strong NP-hardness. Instances are strings over a two-letter alphabet with every number in unary. Graphs are ttt in unary followed by the 3t×3t3t\times 3t3t×3t adjacency matrix, so ∣V∣=3t|V|=3t∣V∣=3t is part of the instance. Languages contain only codes of yes-instances. Strong NP-hardness is NP-hardness of the unary code language. With unary numbers, Max(I)≤Length(I)\mathrm{Max}(I)\le\mathrm{Length}(I)Max(I)≤Length(I), so this is equivalent to Garey and Johnson's definition. The complexity layer is the published module CookPvsNP_defs.
  • Cited hypothesis. Each hardness theorem takes as its only hypothesis the NP-hardness of its source problem, which the paper cites from Garey and Johnson rather than proves. The hypothesis is a true statement about a nonempty, non-universal language. The statements are not weakened to a reduction between languages, and they assume nothing about P versus NP.
  • Source problems. The paper's phrase "three vertices, at most two of which are nonadjacent" is read as "at most one nonadjacent pair", which is Garey and Johnson's GT13. Reading it as "at most two nonadjacent pairs" would admit triples with a single edge and change the problem. PARTITION INTO PATHS OF LENGTH 2 reuses the published definition CubicP3Partition.P3Factor, a spanning non-induced P3P_3P3​-factor.
  • Construction. Resources are indexed by the nonadjacent pairs j<kj<kj<k in lexicographic order, one per unordered pair and none for a pair {j,j}\{j,j\}{j,j}. A diagonal resource would make every job infeasible.
  • Not trivial. A model that checks resources only at integer times, or only at start times, would make the paths equivalence false. A hypothesis on the target problem would make the goal circular. The definitions rule out both.

Welcome contributions: proofs of the two equivalences, polynomial-time computability of the construction on Cook's machines, and a general composition lemma for polynomial-time reductions. The composition lemma is reusable for every hardness mission built on CookPvsNP_defs.

Selected references

  • J. Błażewicz, J. K. Lenstra, A. H. G. Rinnooy Kan, Scheduling subject to resource constraints: classification and complexity, Discrete Applied Mathematics 5 (1983) 11–24. https://doi.org/10.1016/0166-218X(83)90012-4
  • M. R. Garey, D. S. Johnson, Complexity results for multiprocessor scheduling under resource constraints, SIAM Journal on Computing 4 (1975) 397–411. https://doi.org/10.1137/0204035
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, San Francisco, 1979, ISBN 0-7167-1045-5.
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Annals of Discrete Mathematics 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • S. Cook, The P versus NP problem, Clay Mathematics Institute Millennium Problems. https://www.claymath.org/wp-content/uploads/2022/06/pvsnp.pdf
41 thms7 active usersReviewed
🏆Completed
CombinatoricsGraph Theory·Captain: hao jia

Immune High-Girth Bipartite Graphs (Feghali-Lucke-Paulusma-Ries 2025)Research Paper

Motivation

A matching cut is a vertex bipartition whose crossing edges form a matching. The property was introduced under the name decomposability and has links to graph algorithms, stable cutsets in line graphs, and several graph-labeling problems. An Open Problem Garden question asked whether sufficiently large girth forces a matching cut once average degree is bounded.

Feghali, Lucke, Paulusma, and Ries answered that question negatively. Their conference paper appeared at ISAAC 2023, and the version of record was published in Algorithmica in 2025. The paper proves NP-completeness for bipartite graphs of arbitrarily prescribed girth and bounded maximum degree. A central input, Lemma 5, is a stronger structural existence statement: for every girth threshold there is an immune 141414-regular bipartite graph of at least that girth, and it has a perfect matching.

This is therefore a ResearchPaper mission, not a new open-problem mission. Its goal is to formalize the published theorem and its graph-theoretic consequence. Repository candidate constructions and finite arithmetic audits remain separate and are not credited as solving the problem.

Setting

For a finite simple graph GGG and a vertex set A⊆V(G)A\subseteq V(G)A⊆V(G), the associated cut consists of all edges with one endpoint in AAA and one in V(G)∖AV(G)\setminus AV(G)∖A. The cut is nontrivial when both shores are nonempty. It is a matching cut when each vertex is incident with at most one crossing edge. A graph is called immune in the cited paper when it has no matching cut.

The girth is the length of a shortest simple cycle; forests have infinite girth. A graph is 141414-regular when every vertex has exactly fourteen neighbors. Bipartiteness is witnessed by a partition into two independent sides. A perfect matching pairs every vertex with one adjacent partner.

The original OPG wording has a literal one-vertex boundary ambiguity: with nonempty shores required, K1K_1K1​ has no matching cut, average degree zero, and infinite girth. The research-paper target avoids that vacuity by constructing connected graphs with at least two vertices, exact degree fourteen, and arbitrarily large finite girth.

Formalization targets

Lemma 5 — immune high-girth graphs

The main theorem follows the paper's structural lemma:

∀g≥3 ∃G,G is finite, connected, bipartite, and 14-regular,\forall g\ge3\ \exists G, \quad G\text{ is finite, connected, bipartite, and $14$-regular}, ∀g≥3 ∃G,G is finite, connected, bipartite, and 14-regular, girth⁡(G)≥g,G has no matching cut,G has a perfect matching. \operatorname{girth}(G)\ge g, \qquad G\text{ has no matching cut}, \qquad G\text{ has a perfect matching}.girth(G)≥g,G has no matching cut,G has a perfect matching.

The graph may depend on ggg. The existence quantifier does not request an efficient algorithm or a numerical order bound.

Negative OPG consequence

A supporting theorem removes the perfect-matching and bipartite fields and records the direct substantive counterexample family:

∀g≥3 ∃G,d‾(G)=14<15,girth⁡(G)≥g,G has no matching cut.\forall g\ge3\ \exists G, \qquad \overline d(G)=14<15, \quad \operatorname{girth}(G)\ge g, \quad G\text{ has no matching cut}.∀g≥3 ∃G,d(G)=14<15,girth(G)≥g,G has no matching cut.

Thus choosing d=15d=15d=15 refutes the intended universal assertion that some girth threshold works for every graph of average degree below ddd.

Significance

The theorem shows that large girth and bounded degree do not force matching cuts. The examples are highly nontrivial: they are connected, regular, bipartite, and can have arbitrarily large girth. This separates local tree-like structure from the global expansion that prevents a matching cut.

Within the paper, the immune graphs serve as gadgets for hardness reductions. The journal theorem states that, for every g≥3g\ge3g≥3, Matching Cut is NP-complete even for bipartite graphs of girth at least ggg and maximum degree at most 606060. Formalizing Lemma 5 supplies the graph-theoretic core needed to reconstruct that result without forcing this mission to formalize an entire complexity-theory reduction in its first stage.

Difficulty

Large girth alone makes bounded neighborhoods look like trees, and trees have many matching cuts. Immunity must therefore come from global expansion rather than short local cycles. The paper obtains the required family from Lubotzky–Phillips–Sarnak Ramanujan graphs and combines spectral and isoperimetric bounds to show that every nontrivial cut has too many crossing incidences to be a matching.

A formal proof must bridge several exact interfaces: existence of suitable primes, the finite Cayley-graph construction, bipartiteness and regularity, the girth lower bound, the spectral-to-isoperimetric inequality, and Hall's theorem for the perfect matching. None of these can be replaced by a finite sample or an asymptotic slogan.

Formalization scope

Graphs are finite and simple. Connectedness is nonempty mutual graph reachability. A simple cycle is a cyclic list of at least three distinct vertices; girth at least ggg means every such cycle has length at least ggg, so forests satisfy every threshold. A matching cut requires both shores nonempty and is encoded by the condition that every vertex has at most one crossing neighbor. A perfect matching is represented by an adjacent involution.

The main theorem explicitly requires at least two vertices, although exact 14-regularity already forces nontrivial order; the redundant bound documents exclusion of the K1K_1K1​ ambiguity. The mission does not claim that the frozen OPG contract was well-posed at order one. It formalizes the paper's substantive counterexample family and the consequence for the intended question.

Candidate files in the associated repository explore alternative bounded-degree constructions and integer counts. They are candidate_only and are not proof dependencies. Contributions should follow the published Lemma 5 and its cited inputs, or provide a separately sourced proof of the same declaration. The later maximum-degree-60 NP-completeness theorem is welcome as a future extension after the finite complexity framework is fixed.

Selected references

  • C. Feghali, F. Lucke, D. Paulusma, and B. Ries, Matching Cuts in Graphs of High Girth and H-Free Graphs, Algorithmica 87 (2025), 1199–1221. https://doi.org/10.1007/s00453-025-01318-8
  • C. Feghali, F. Lucke, D. Paulusma, and B. Ries, ISAAC 2023 version. https://drops.dagstuhl.de/entities/document/10.4230/LIPIcs.ISAAC.2023.31
  • A. Lubotzky, R. Phillips, and P. Sarnak, Ramanujan graphs, Combinatorica 8 (1988), 261–277. https://doi.org/10.1007/BF02126799
  • Open Problem Garden, Matching cut and girth. https://www.openproblemgarden.org/op/matching_cut_and_girth
20 thms6 active usersReviewed
Linear OptimizationOperations ResearchOptimization·Captain: ORdos

Smale's Ninth Problem: Strongly Polynomial Linear ProgrammingOpen Problem

The problem of solving linear inequalities

The linear feasibility problem takes a matrix A∈Rm×nA \in \mathbb{R}^{m\times n}A∈Rm×n and a vector b∈Rmb \in \mathbb{R}^mb∈Rm and asks whether the system of mmm linear inequalities in nnn real unknowns

{ x∈Rn∣Ax≥b }  ≠  ∅\{\,x \in \mathbb{R}^n \mid Ax \ge b\,\} \;\ne\; \emptyset{x∈Rn∣Ax≥b}=∅

has a solution. By linear programming duality, optimizing a linear objective over such a set reduces to feasibility, so this decision problem carries the whole complexity of linear programming.

What "polynomial time" means here depends on the machine. In the bit model the input is a list of rational numbers, its size LLL counts the bits of all numerators and denominators, and an algorithm is polynomial if it runs in time poly(m,n,L)\mathrm{poly}(m, n, L)poly(m,n,L). In the real-number model the input is a list of mn+mmn + mmn+m exact real numbers, each arithmetic operation (+,−,×,÷+, -, \times, \div+,−,×,÷), comparison, or memory move costs one unit, and a running time may only depend on mmm and nnn. An algorithm polynomial in this second sense is what Smale asks for; the closely related bit-model notion — poly(m,n)\mathrm{poly}(m,n)poly(m,n) arithmetic operations and polynomially bounded intermediate bit sizes — is called strongly polynomial. This mission fixes the real-number model precisely as a Blum–Shub–Smale (BSS) machine (Blum–Shub–Smale 1989): a finite program of instructions acting on a bi-infinite tape Z→R\mathbb{Z} \to \mathbb{R}Z→R of real registers — loads of arbitrary real machine constants, exact field arithmetic at fixed addresses, two-sided tape shifts, a sign-test branch, and accept/reject — with cost equal to the number of executed instructions. The convention that costs something: the program must be uniform, one finite instruction list serving every mmm, nnn, and every real instance. Uniformity is exactly what separates the question from point-location tricks available to non-uniform families of decision trees.

Why it matters

For optimization, the question is the last gap in the complexity of its central problem. Linear programs with combinatorial structure already admit strongly polynomial algorithms — Tardos (1986) solved every LP whose running time may depend on the entries of AAA but not on bbb or ccc, covering network flows and all {0,±1}\{0,\pm1\}{0,±1}-constraint problems — and a positive answer for general LP would extend that unification to the whole class, while explaining why simplex-type methods behave so well in practice (Spielman–Teng 2004).

For the theory of computation over the reals, the problem is a benchmark for what unit-cost exact arithmetic can do: it is Problem 9 on Smale's list of mathematical problems for the twenty-first century (Smale 1998), posed in the BSS model as the real-number analogue of the P-versus-NP style questions of that program, and it interacts with polyhedral combinatorics through the polynomial Hirsch conjecture: a polynomial bound on polytope diameters is a necessary condition for any polynomial pivot rule. A problem that calibrates both the practice of optimization and the foundations of real computation is a subject, not a special case.

The question and what is known

Question (Smale’s 9th).Is there a uniform BSS program deciding {x∣Ax≥b}≠∅ in poly(m,n) steps?\textbf{Question (Smale's 9th).}\quad \text{Is there a uniform BSS program deciding } \{x \mid Ax \ge b\} \ne \emptyset \text{ in } \mathrm{poly}(m,n) \text{ steps?}Question (Smale’s 9th).Is there a uniform BSS program deciding {x∣Ax≥b}=∅ in poly(m,n) steps?

The timeline splits into a negative branch (lower bounds against algorithm classes) and a positive branch (polynomial algorithms in weaker senses).

Lower bounds. Klee–Minty (1972) constructed a deformed cube on which Dantzig's largest-coefficient simplex rule visits all 2n2^n2n vertices; analogous exponential examples were later found for essentially every deterministic pivot rule, and randomized rules were driven to subexponential lower bounds by Friedmann–Hansen–Zwick (2011) — against upper bounds of exp⁡(O(nlog⁡n))\exp(O(\sqrt{n \log n}))exp(O(nlogn​)) from Kalai (1992) and Matoušek–Sharir–Welzl (1996). On the interior-point side, Allamigeon–Benchimol–Gaubert–Joswig (2018) showed by tropical methods that log-barrier path following is not strongly polynomial, and Allamigeon–Gaubert–Vandame (2022) extended this to every self-concordant barrier: no interior-point method of that class can settle the question positively.

Polynomial algorithms in weaker senses. Khachiyan (1979/80) proved LP feasibility is polynomial in the bit model via the ellipsoid method; Karmarkar (1984) and then Renegar (1988) brought interior-point methods to O(n L)O(\sqrt{n}\,L)O(n​L) iterations. Megiddo (1984) solved LP in linear time for every fixed dimension; Tardos (1986) gave the combinatorial strongly polynomial class; Vavasis–Ye (1996) and Dadush–Huiberts–Natura–Végh (2020) replaced the bit size by condition measures of AAA alone; Ye (2011) proved policy iteration strongly polynomial for fixed-discount Markov decision processes.

The central difficulty is visible in every positive result: each known iteration count is controlled by a scale-dependent quantity — bit length, condition number, barrier curvature — that is unbounded over the real instances with m,nm, nm,n fixed. The naive plan, "run the ellipsoid method and round", fails at its first step in the real model: the number of iterations needed to separate a feasible system from an infeasible one grows with the thinness of the feasible set, which is not a function of (m,n)(m, n)(m,n); no data-independent perturbation ε\varepsilonε exists when the data are arbitrary reals. All results above are proved on paper only; none has a machine-checked proof in the literature. What is already formalized, on this platform, is the substrate this mission builds on: the simplex iteration (mission Introduction to Linear Optimization IV), the ellipsoid method with its volume-halving correctness theorem (XI), interior-point path following (XII), and self-concordance with the barrier method (Convex Optimization VI).

A hierarchy of formalization targets

The mission's milestone list realizes this hierarchy in order; each level states what it deliberately leaves open.

Level 0 — the model works. A uniform BSS program decides one-variable feasibility in linear time:

∃ P, C  ∀m, ∀(a,b)∈Rm×Rm: P decides {x∈R∣aix≥bi ∀i}≠∅ within C(m+1) steps.\exists\,P,\,C\ \ \forall m,\ \forall (a,b) \in \mathbb{R}^m \times \mathbb{R}^m:\ P \text{ decides } \{x \in \mathbb{R} \mid a_i x \ge b_i\ \forall i\} \ne \emptyset \text{ within } C(m{+}1) \text{ steps}.∃P,C  ∀m, ∀(a,b)∈Rm×Rm: P decides {x∈R∣ai​x≥bi​ ∀i}=∅ within C(m+1) steps.

It fixes nothing about n≥2n \ge 2n≥2; its role is to certify that the machine model and cost semantics of the goal are non-vacuous.

Level 1 — the classical method is exponential. On the Klee–Minty cube, Dantzig's rule admits a run of

2n−1 pivots2^n - 1 \text{ pivots}2n−1 pivots

from the all-slack basis to the optimum. It leaves open all other pivot rules — extensions to further rules are welcome as strengthenings.

Level 2 — the bit model succeeds. Through the Cramer–Hadamard solution bound ∣xj∣≤n! Un|x_j| \le n!\,U^n∣xj​∣≤n!Un and the perturbation estimates, Khachiyan's theorem: for integer data bounded by UUU, every admissible ellipsoid run decides feasibility within

t∗≤106 (n+2)4(log⁡2U+n+2) iterations.t^* \le 10^6\,(n{+}2)^4(\log_2 U + n + 2) \text{ iterations}.t∗≤106(n+2)4(log2​U+n+2) iterations.

The generous constants are deliberate — only the polynomial order is load-bearing. This level leaves open exactly the dependence on log⁡U\log UlogU.

Level 3 — the goal (open). A uniform program with data-independent polynomial cost:

∃ P, C, d  ∀m,n,A,b: P decides {x∣Ax≥b}≠∅ within C (mn+m+2)d steps.\exists\,P,\,C,\,d\ \ \forall m, n, A, b:\ P \text{ decides } \{x \mid Ax \ge b\} \ne \emptyset \text{ within } C\,(mn + m + 2)^d \text{ steps}.∃P,C,d  ∀m,n,A,b: P decides {x∣Ax≥b}=∅ within C(mn+m+2)d steps.

The statement asserts only the shape of the truth — no hard-coded degree or constant — so it is stable under every future quantitative improvement. These levels do not exhaust the project: Tardos' combinatorial LP theorem, Ye's fixed-discount MDP result, and impossibility statements for restricted program classes in the style of Allamigeon–Gaubert–Vandame are natural later milestones.

Formalization scope

Polyhedra, simplex states, pivots, and ellipsoid runs are the platform's existing LinearOptimization development over Matrix (Fin m) (Fin n) ℝ, with {x∣Ax≥b}\{x \mid Ax \ge b\}{x∣Ax≥b} as polyhedron A b; algorithms with data-dependent iteration counts are formalized as run predicates, as in the parent missions. The new SmaleNinth definitions supply what the goal genuinely needs and the run-predicate style cannot express: a concrete inductive type of BSS programs with operational semantics and unit-cost accounting, the Klee–Minty data with Dantzig's rule, and the explicit Khachiyan constants. One convention closes the degenerate escape hatch: the goal quantifies over finite BSSProgram terms under the fixed encodeLP input convention — formalizing "algorithm" as an arbitrary function Rmn+m→Bool\mathbb{R}^{mn+m} \to \mathrm{Bool}Rmn+m→Bool would make the statement trivially true and is not the theorem. Division is totalized as x/0=0x/0 = 0x/0=0 and the branch test is xi≤0x_i \le 0xi​≤0; both are benign for the class of programs quantified over.

The machine module is infrastructure beyond this mission — any real-number complexity statement (other Smale problems, sums-of-square-roots, BSS-completeness) can reuse it, as can any pivot-rule lower bound reuse the Klee–Minty module. Formalization forces distinctions the literature leaves informal: which machine variant carries the unit-cost claim, how ties in Dantzig's rule are resolved, and which of the interchangeable Khachiyan constants each estimate actually needs. Welcome contributions include proofs of any milestone, alternative exponential instances for other pivot rules, sharper constants in the Khachiyan module, and ports of the known strongly polynomial special cases.

Selected references

  • L. Blum, M. Shub, S. Smale, On a theory of computation and complexity over the real numbers, Bull. AMS 21(1):1–46, 1989. DOI
  • S. Smale, Mathematical problems for the next century, Math. Intelligencer 20(2):7–15, 1998. DOI
  • V. Klee, G. J. Minty, How good is the simplex algorithm?, in Inequalities III, Academic Press, 1972, pp. 159–175.
  • L. G. Khachiyan, Polynomial algorithms in linear programming, USSR Comput. Math. Math. Phys. 20:53–72, 1980. DOI
  • N. Karmarkar, A new polynomial-time algorithm for linear programming, Combinatorica 4:373–395, 1984. DOI
  • J. Renegar, A polynomial-time algorithm, based on Newton's method, for linear programming, Math. Programming 40:59–93, 1988. DOI
  • É. Tardos, A strongly polynomial algorithm to solve combinatorial linear programs, Oper. Res. 34(2):250–256, 1986. DOI
  • N. Megiddo, Linear programming in linear time when the dimension is fixed, J. ACM 31(1):114–127, 1984. DOI
  • G. Kalai, A subexponential randomized simplex algorithm, STOC 1992. DOI
  • O. Friedmann, T. D. Hansen, U. Zwick, Subexponential lower bounds for randomized pivoting rules for the simplex algorithm, STOC 2011. DOI
  • D. A. Spielman, S.-H. Teng, Smoothed analysis of algorithms: why the simplex algorithm usually takes polynomial time, J. ACM 51(3):385–463, 2004. DOI
  • S. A. Vavasis, Y. Ye, A primal-dual interior point method whose running time depends only on the constraint matrix, Math. Programming 74:79–120, 1996. DOI
  • Y. Ye, The simplex and policy-iteration methods are strongly polynomial for the Markov decision problem with a fixed discount rate, Math. Oper. Res. 36(4):593–603, 2011. DOI
  • X. Allamigeon, P. Benchimol, S. Gaubert, M. Joswig, Log-barrier interior point methods are not strongly polynomial, SIAM J. Appl. Algebra Geom. 2(1):140–178, 2018. DOI
  • X. Allamigeon, S. Gaubert, N. Vandame, No self-concordant barrier interior point method is strongly polynomial, STOC 2022. arXiv
  • D. Dadush, S. Huiberts, B. Natura, L. A. Végh, A scaling-invariant algorithm for linear programming whose running time depends only on the constraint matrix, STOC 2020. arXiv
  • D. Bertsimas, J. N. Tsitsiklis, Introduction to Linear Optimization, Athena Scientific, 1997 (Chapters 3, 8, 9 — formalized in the Introduction to Linear Optimization mission series).
  • B. Korte, J. Vygen, Combinatorial Optimization: Theory and Algorithms, 6th ed., Springer, 2018, §4.1–4.5.
29 thms6 active usersReviewed
🏆Completed
Captain: joe

The Sipser–Gács–Lautemann TheoremResearch Paper

Randomness appears to enlarge efficient computation, but the Sipser–Gács–Lautemann theorem places every bounded-error probabilistic polynomial-time language at the second level of the polynomial hierarchy, giving one of complexity theory’s foundational limits on the power of randomization.

66 thms6 active usersReviewed
🏆Completed
Machine LearningOperations Research·Captain: mikedeng1

The Weighted Majority Algorithm 8: With η Anomalies Every Deterministic Learner Can Be Forced to opt(F, 0) + 2η MistakesResearch Paper

Learning with anomalies

In the on-line prediction model studied by Littlestone and Warmuth (Inform. and Comput. 108, 1994), a learner receives instances one at a time, predicts a binary label for each, and is then told the correct label. No probabilistic assumption is made: performance is measured by the number of mistakes on a worst-case sequence. When every label is given by some function of a known target class FFF, the best attainable worst-case number of mistakes is a combinatorial invariant of FFF, characterized by Littlestone (Mach. Learn. 2, 1988) through mistake trees, now called the Littlestone dimension.

Real data are rarely labelled exactly by a member of the class. Littlestone and Warmuth measure the departure by the number of anomalies: the number of trials inconsistent with the best function in the class. Their Weighted Majority algorithm makes a number of mistakes that grows linearly in the number of anomalies, and the rate of growth can be pushed towards 222 mistakes per anomaly by the choice of the update factor β\betaβ. Theorem 8.1 of the paper shows that no deterministic learner does better: rate 222 is forced for every class with at least two members. This mission formalizes that lower bound.

Setting

Let XXX be an arbitrary set (the domain) and FFF a set of functions X→{0,1}X \to \{0,1\}X→{0,1}, possibly infinite. A trial is a pair (x,y)∈X×{0,1}(x,y) \in X \times \{0,1\}(x,y)∈X×{0,1}, and a sequence of trials S=((x1,y1),…,(xT,yT))\mathcal S = ((x_1,y_1),\dots,(x_T,y_T))S=((x1​,y1​),…,(xT​,yT​)) is finite; an instance may occur several times with different labels.

  • The number of anomalies of S\mathcal SS with respect to f:X→{0,1}f : X \to \{0,1\}f:X→{0,1} is anomf(S)=#{t≤T:f(xt)≠yt}\mathrm{anom}_f(\mathcal S) = \#\{t \le T : f(x_t) \ne y_t\}anomf​(S)=#{t≤T:f(xt​)=yt​}.
  • SηS_\etaSη​ is the collection of sequences with at most η\etaη anomalies with respect to FFF: those S\mathcal SS with min⁡f∈Fanomf(S)≤η\min_{f\in F} \mathrm{anom}_f(\mathcal S) \le \etaminf∈F​anomf​(S)≤η.
  • A deterministic on-line algorithm AAA maps the list of earlier trials and the current instance to a predicted label. MA(S)M_A(\mathcal S)MA​(S) is the number of trials ttt at which AAA's prediction, computed from (x1,y1),…,(xt−1,yt−1)(x_1,y_1),\dots,(x_{t-1},y_{t-1})(x1​,y1​),…,(xt−1​,yt−1​) and xtx_txt​, differs from yty_tyt​.
  • opt(F,η)=min⁡Asup⁡S∈SηMA(S)∈{0,1,… }∪{∞}\mathrm{opt}(F,\eta) = \min_A \sup_{\mathcal S \in S_\eta} M_A(\mathcal S) \in \{0,1,\dots\} \cup \{\infty\}opt(F,η)=minA​supS∈Sη​​MA​(S)∈{0,1,…}∪{∞}.

A complete kkk-mistake tree for FFF is a complete binary tree of height kkk whose internal nodes carry instances of XXX such that every labelling of a root-to-leaf path is realized by some f∈Ff \in Ff∈F; in Lean this is the predicate ShattersTree F k v.

Formalization targets

Goal: Theorem 8.1 (p. 250)

For every class FFF with ∣F∣>1|F| > 1∣F∣>1 and every η≥0\eta \ge 0η≥0,

opt(F,η) ≥ opt(F,0)+2η.\mathrm{opt}(F,\eta) \ \ge\ \mathrm{opt}(F,0) + 2\eta .opt(F,η) ≥ opt(F,0)+2η.

The constant 222 is the paper's, and the statement holds for infinite classes, including those with opt(F,0)=∞\mathrm{opt}(F,0) = \inftyopt(F,0)=∞.

Milestones

  1. The [Lit88, Lit89b] mistake tree (§8, p. 251): for any pool FFF there is a complete opt(F,0)\mathrm{opt}(F,0)opt(F,0)-mistake tree. This is published and proved on the platform as UnderstandingML.soa_mistake_bound (MSOA(F)≤Ldim(F)M_{\mathrm{SOA}}(F) \le \mathrm{Ldim}(F)MSOA​(F)≤Ldim(F)), referenced together with UnderstandingML.ldim_le_mistake_bound; together they identify opt(F,0)\mathrm{opt}(F,0)opt(F,0) with the Littlestone dimension.
  2. The first stage of the adversary (§8, p. 251): if FFF has a complete kkk-mistake tree with k≥1k \ge 1k≥1, then against any deterministic algorithm there are k−1k-1k−1 trials on which it errs every time, an instance xxx, and two functions f1,f2∈Ff_1, f_2 \in Ff1​,f2​∈F consistent with all these trials with f1(x)≠f2(x)f_1(x) \ne f_2(x)f1​(x)=f2​(x).

Significance

The theorem settles the optimal rate at which mistakes must grow with the number of anomalies for deterministic learners: 222 per anomaly, uniformly over all non-trivial classes. Together with the upper bounds of §2 of the paper, it shows that the deterministic Weighted Majority algorithm is within an arbitrarily small factor of optimal in its dependence on the anomalies, and it explains the factor of 222 that separates deterministic from randomized prediction (Theorem 8.2, a separate mission of this series, gives rate 111 for randomized learners).

The result is proved in the paper; to the best of a search of the platform it has not been formalized. The mission produces a machine-checked lower bound in a general setting: arbitrary domain, arbitrary (possibly infinite) class, extended-natural-valued optimal bounds. The noise-free side, the Littlestone dimension and the Standard Optimal Algorithm, is already published on the platform from Shalev-Shwartz and Ben-David's textbook, and this mission connects it to the noisy setting.

Difficulty

The obvious argument forces one mistake per trial and stops: an adversary that always contradicts the learner leaves no function of FFF with few anomalies once the trials outnumber the mistake tree. The difficulty is to spend the anomaly budget so that each anomaly buys two mistakes while the sequence stays in SηS_\etaSη​, and to do so against an arbitrary algorithm, which may behave adversarially itself. The definition of SηS_\etaSη​ through a minimum over the class (not a fixed target) and the possibility that opt(F,0)\mathrm{opt}(F,0)opt(F,0) is infinite both have to be handled.

Formalization scope

  • The domain is any type X; the class is F : Set (X → Bool); the hypothesis ∣F∣>1|F| > 1∣F∣>1 is F.Nontrivial, which is meaningful for infinite classes (a cardinality ncard would be 000 there).
  • Algorithms, histories and mistake counts are those of the published module UnderstandingML_Online (OnlineAlg X Bool, history, mistakes, ShattersTree); an algorithm is an arbitrary function of the history and the current instance, so every deterministic learner is covered.
  • Sequences are Fin T → X × Bool for every finite length T, with repeated instances and conflicting labels allowed.
  • opt F η : ℕ∞ is an infimum over algorithms of a supremum over sequences in SηS_\etaSη​; η\etaη is a natural number. Membership in SηS_\etaSη​ is "some f∈Ff \in Ff∈F has at most η\etaη anomalies", which is the paper's minimum over FFF.
  • Swapping the quantifiers (one sequence hard for all algorithms) or replacing opt\mathrm{opt}opt by a finite toNat truncation would give a different statement; neither is used.

Needed infrastructure: the identification opt(F,0)=min⁡AMA(F)\mathrm{opt}(F,0) = \min_A M_A(F)opt(F,0)=minA​MA​(F) (sequences with no anomaly are exactly the realizable ones), the existence of a shattered tree of depth opt(F,0)\mathrm{opt}(F,0)opt(F,0) from the published Littlestone-dimension results, and the two-stage adversary. Contributions of general lemmas about ShattersTree (subtrees, truncation to smaller depth) are reusable for the randomized mission and beyond.

Selected references

  • N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 212–261, 1994. https://doi.org/10.1006/inco.1994.1009
  • N. Littlestone, Learning quickly when irrelevant attributes abound: a new linear-threshold algorithm, Machine Learning 2, 285–318, 1988. https://doi.org/10.1007/BF00116827
  • S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 21. https://doi.org/10.1017/CBO9781107298019
6 thms5 active usersReviewed
🏆Completed
CombinatoricsGraph Theory·Captain: mikedeng1

Eigenvalues and Expanders: A Regular Bipartite Graph Is a Strong Expander If and Only If λ(G) Is Bounded Away from 0Research Paper

Motivation

Expander graphs are sparse graphs in which every set of vertices has many neighbours. Families of them with bounded degree and expansion bounded away from zero are a basic tool of theoretical computer science: they are the main component of the sorting network of Ajtai, Komlós and Szemerédi (AKS 1983), the building block of superconcentrators and other graphs with strong connectivity properties, and an ingredient of many later constructions in coding theory, derandomization and complexity.

Expansion is hard to certify. Checking that every set of vertices has many neighbours means looking at exponentially many sets, and computing the exact expansion of a graph is coNP-complete. A spectral quantity, by contrast, is computable in polynomial time. N. Alon's paper Eigenvalues and expanders (Combinatorica 6 (1986) 83–96) proves that for regular bipartite graphs the two notions are equivalent: a graph is a strong expander if and only if the second-smallest eigenvalue of its Laplacian is bounded away from 0, with explicit constants in both directions. This is a discrete counterpart of Cheeger's inequality for Riemannian manifolds.

Timeline.

  • 1983: Ajtai, Komlós and Szemerédi use bounded-degree bipartite expanders to build sorting networks of depth O(log⁡n)O(\log n)O(logn).
  • 1984: Tanner (SIAM J. Alg. Disc. Meth. 5) bounds the neighbourhood size of a set in a regular bipartite graph by its second eigenvalue, the direction "eigenvalue gap implies expansion".
  • 1985: Alon and Milman (J. Combin. Theory Ser. B 38) prove isoperimetric inequalities for graphs in terms of λ(G)\lambda(G)λ(G) and introduce enlargers.
  • 1986: Alon proves the converse direction, "expansion implies an eigenvalue gap" (Lemma 2.4 and Theorem 3.4 of the paper).

Setting

All graphs are finite and simple. For a graph G=(V,E)G = (V, E)G=(V,E) and a set X⊆VX \subseteq VX⊆V, N(X)={v∈V:vx∈E for some x∈X}N(X) = \{v \in V : vx \in E \text{ for some } x \in X\}N(X)={v∈V:vx∈E for some x∈X} is the set of neighbours of XXX; it may meet XXX.

The Laplacian of GGG is QG=diag(d(v))v∈V−AGQ_G = \mathrm{diag}(d(v))_{v \in V} - A_GQG​=diag(d(v))v∈V​−AG​, where AGA_GAG​ is the 0–1 adjacency matrix and d(v)d(v)d(v) the degree of vvv. It is symmetric with eigenvalues 0=λ0≤λ1≤⋯≤λn−10 = \lambda_0 \le \lambda_1 \le \dots \le \lambda_{n-1}0=λ0​≤λ1​≤⋯≤λn−1​, counted with multiplicity, and λ(G)=λ1\lambda(G) = \lambda_1λ(G)=λ1​ is its second-smallest eigenvalue. It is positive exactly when GGG is connected.

  • An (n,d,c)(n, d, c)(n,d,c)-magnifier is a graph on nnn vertices with maximal degree ddd in which every X⊆VX \subseteq VX⊆V with ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 satisfies ∣N(X)−X∣≥c∣X∣|N(X) - X| \ge c|X|∣N(X)−X∣≥c∣X∣.
  • An (n,d,ε)(n, d, \varepsilon)(n,d,ε)-enlarger is a graph on nnn vertices with maximal degree ddd and λ(G)≥ε\lambda(G) \ge \varepsilonλ(G)≥ε.
  • A bipartite graph G=(I,O;E)G = (I, O; E)G=(I,O;E) has inputs III, outputs OOO and edges only between III and OOO. It is a strong (n,d,c)(n, d, c)(n,d,c)-expander if ∣I∣=∣O∣=n|I| = |O| = n∣I∣=∣O∣=n, the maximal degree is ddd, and for every X⊆IX \subseteq IX⊆I
∣N(X)∣≥(1+c(1−∣X∣n))∣X∣.|N(X)| \ge \Bigl(1 + c\Bigl(1 - \frac{|X|}{n}\Bigr)\Bigr)|X|.∣N(X)∣≥(1+c(1−n∣X∣​))∣X∣.

Formalization targets

Goal: Theorem 3.4

Let G=(I,O;E)G = (I, O; E)G=(I,O;E) be a ddd-regular bipartite graph with ∣I∣=∣O∣=n|I| = |O| = n∣I∣=∣O∣=n and λ=λ(G)\lambda = \lambda(G)λ=λ(G).

  1. If GGG is a strong (n,d,c)(n, d, c)(n,d,c)-expander then
λ≥c21024+2c2.\lambda \ge \frac{c^2}{1024 + 2c^2}.λ≥1024+2c2c2​.
  1. If λ≥ε\lambda \ge \varepsilonλ≥ε then GGG is a strong (n,d,c)(n, d, c)(n,d,c)-expander with
c=2dε−ε2d2.c = \frac{2d\varepsilon - \varepsilon^2}{d^2}.c=d22dε−ε2​.

Milestones, in the order of the paper

  • Lemma 2.2 (Alon–Milman, already on the platform): for disjoint sets A,BA, BA,B at distance ϱ>1\varrho > 1ϱ>1, b≤(1−a)/(1+(λ/d)aϱ2)b \le (1-a)/(1 + (\lambda/d)a\varrho^2)b≤(1−a)/(1+(λ/d)aϱ2).
  • Corollary 2.3: every (n,d,ε)(n, d, \varepsilon)(n,d,ε)-enlarger is an (n,d,2ε/(d+2ε))(n, d, 2\varepsilon/(d+2\varepsilon))(n,d,2ε/(d+2ε))-magnifier.
  • Eq. (2.1): if fff is an eigenvector of QGQ_GQG​ for λ(G)\lambda(G)λ(G) and ggg its positive part, then ∑uv∈E(g(u)−g(v))2≤λ∑vg2(v)\sum_{uv \in E}(g(u)-g(v))^2 \le \lambda \sum_v g^2(v)∑uv∈E​(g(u)−g(v))2≤λ∑v​g2(v).
  • Lemma 2.4: every (n,d,c)(n, d, c)(n,d,c)-magnifier has λ(G)≥c2/(4+2c2)\lambda(G) \ge c^2/(4 + 2c^2)λ(G)≥c2/(4+2c2).
  • Lemma 3.1: a strong (n,d,c)(n, d, c)(n,d,c)-expander is a (2n,d,c/16)(2n, d, c/16)(2n,d,c/16)-magnifier.
  • Proof of Lemma 3.3, spectrum: the two largest eigenvalues of CTCC^TCCTC, with CCC the I×OI \times OI×O biadjacency matrix, are d2d^2d2 and (d−λ)2(d - \lambda)^2(d−λ)2.
  • Proof of Lemma 3.3, Tanner's bound: ∣N(X)∣≥d2∣X∣/(α(d2−(d−λ)2)+(d−λ)2)|N(X)| \ge d^2|X| / \bigl(\alpha(d^2 - (d-\lambda)^2) + (d-\lambda)^2\bigr)∣N(X)∣≥d2∣X∣/(α(d2−(d−λ)2)+(d−λ)2) with α=∣X∣/n\alpha = |X|/nα=∣X∣/n.
  • Lemma 3.3: a ddd-regular bipartite graph is a strong (n,d,(2dλ−λ2)/d2)(n, d, (2d\lambda - \lambda^2)/d^2)(n,d,(2dλ−λ2)/d2)-expander.

Part (1) of the goal combines Lemmas 3.1 and 2.4; part (2) follows from Lemma 3.3.

Significance

The result. Theorem 3.4 makes expansion of regular bipartite graphs checkable in polynomial time up to a constant-factor loss. A random regular bipartite graph can be generated and its expansion certified by computing one eigenvalue. Lemma 2.4 is one of the first discrete Cheeger inequalities. Together with Corollary 2.3 it shows that magnifiers and enlargers are the same graphs up to the constants, and it underlies the later theory of spectral expanders, including the Alon–Boppana bound and Ramanujan graphs.

The formalization. The results are proved in the paper. The work is to formalize the known proofs. That includes Tanner's eigenvalue bound, which the paper only cites, and a max-flow min-cut argument, for which Mathlib has no general theorem. No machine-checked version of Lemma 2.4, Lemma 3.1, Lemma 3.3 or Theorem 3.4 is known to exist. Lemma 2.2 is already stated and proved on the platform as part of the Alon–Milman mission.

Difficulty

The direction "eigenvalue gap implies expansion" is a variational argument on the spectrum of CTCC^TCCTC. The converse is the hard one. A first attempt bounds λ(G)\lambda(G)λ(G) from below by testing the Rayleigh quotient on indicator vectors of sets. That only gives upper bounds on λ\lambdaλ: any one test vector does. A lower bound has to control every vector orthogonal to the constants at once. The paper first reduces to the positive part of an eigenvector (Eq. (2.1)). It then turns the combinatorial expansion of the graph into an analytic inequality for that function, using a network flow whose existence comes from the max-flow min-cut theorem. Step (ii) of the flow conditions printed on p. 87 is false as stated: the arcs (u,u)(u, u)(u,u) of the network absorb part of the flow. The flow argument has to be repaired before it can be formalized.

Lemma 3.1 has its own obstacle: one-sided expansion of inputs must be converted into expansion of arbitrary vertex sets that mix inputs and outputs. This needs the strong form of expansion; for ordinary expanders the lemma is false.

Formalization scope

Graphs are SimpleGraph V on a Fintype. A bipartite graph lives on the sum type I ⊕ O, and IsIOBipartite forbids edges inside I and inside O. Cardinalities ∣N(X)∣|N(X)|∣N(X)∣ are Set.ncard. "Maximal degree ddd" is read as G.maxDegree ≤ d; every statement is monotone in ddd or fixes ddd by regularity (G.IsRegularOfDegree d). The condition ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 is written 2∣X∣≤n2|X| \le n2∣X∣≤n in N\mathbb{N}N. All constants are real, and every subtraction and division is taken in R\mathbb{R}R.

Reused published items:

  • λ(G)\lambda(G)λ(G) is AlonMilman.Diameter.lambda1, the second-smallest eigenvalue of G.lapMatrix ℝ, which is 000 by convention on fewer than two vertices.
  • N(X)N(X)N(X) is AKSSorting.Core.neighbours.
  • Lemma 2.2 is AlonMilman.Diameter.theorem_2_5, which carries Alon–Milman's standing hypotheses that GGG is connected and n≥2n \ge 2n≥2; outside them the inequality is trivial.

The page omits a few degenerate cases, and the following hypotheses are added for them. Each is necessary, with a counterexample recorded in the item's statement:

  • n≥1n \ge 1n≥1 and c≥0c \ge 0c≥0 in Theorem 3.4 (1);
  • ε>0\varepsilon > 0ε>0 in Theorem 3.4 (2);
  • n≥2n \ge 2n≥2 and c≥0c \ge 0c≥0 in Lemma 2.4;
  • n≥2n \ge 2n≥2 in Lemma 3.1 and in the CTCC^TCCTC statement;
  • d≥1d \ge 1d≥1 in Lemma 3.3;
  • ε≥0\varepsilon \ge 0ε≥0 in Corollary 2.3.

Eq. (2.1) is stated in multiplied form, so no quotient by ∑g2\sum g^2∑g2 appears.

The closing sentences of Theorem 2.5 and Theorem 3.4 ("Thus … one can prove efficiently …") are not formalized. Read as implications between expanders they reduce to monotonicity in ccc, because c′≤cc' \le cc′≤c; their content is algorithmic.

Several encodings would trivialize the mission and are excluded:

  • a λ\lambdaλ other than the published second-smallest Laplacian eigenvalue, in particular one defined as the best constant of a quotient;
  • expansion or magnifier conditions with a negative constant in a hypothesis;
  • ∣X∣≤n/2|X| \le n/2∣X∣≤n/2 with truncating natural-number division.

Contributions are welcome:

  • Tanner's bound in Lean, which is reusable for any regular bipartite graph;
  • a max-flow min-cut theorem for finite networks;
  • the corrected flow lemma behind Eqs. (2.2)–(2.3);
  • the spectral facts about λ(G)\lambda(G)λ(G) for bipartite graphs (λ≤d\lambda \le dλ≤d for n≥2n \ge 2n≥2, and λ=d−σ2(C)\lambda = d - \sigma_2(C)λ=d−σ2​(C)).

Selected references

  • N. Alon, Eigenvalues and expanders, Combinatorica 6 (1986) 83–96. https://doi.org/10.1007/BF02579166
  • N. Alon and V. D. Milman, λ₁, isoperimetric inequalities for graphs, and superconcentrators, J. Combin. Theory Ser. B 38 (1985) 73–88. https://doi.org/10.1016/0095-8956(85)90092-9
  • R. M. Tanner, Explicit concentrators from generalized N-gons, SIAM J. Algebraic Discrete Methods 5 (1984) 287–293. https://doi.org/10.1137/0605030
  • M. Ajtai, J. Komlós and E. Szemerédi, Sorting in c log n parallel steps, Combinatorica 3 (1983) 1–19. https://doi.org/10.1007/BF02579338
16 thms5 active usersReviewed
🏆Completed
Linear OptimizationOperations ResearchProbability·Captain: mikedeng1

Competitive Randomized Algorithms for Nonuniform Problems III: The Optimal Randomized Two-Server Ratio on the 1-d-d Isosceles TriangleResearch Paper

Motivation

The k-server problem of Manasse, McGeoch and Sleator (J. Algorithms 11 (1990)) asks how kkk mobile servers in a metric space should respond, on-line, to a sequence of requests at points of the space, each of which must be covered by a server. It is the central model of on-line computation: paging is the special case of a uniform metric, and many caching and scheduling problems reduce to it. For two servers the deterministic picture is complete: the optimal competitive ratio is 222 on every metric space with at least three points.

Randomization changes the picture, and the smallest nontrivial case already shows how. On the equilateral triangle the optimal randomized ratio against an oblivious adversary is 3/23/23/2. Karlin, Manasse, McGeoch and Owicki (Algorithmica 11 (1994) 542–571) computed the exact optimal randomized ratio for several nonuniform triangles, where the distances differ, and showed that it depends on the geometry. Their Theorem 12 settles the whole family of isosceles triangles with edge lengths 111, ddd, ddd. These exact values are among the few known optimal randomized ratios for server problems.

Timeline:

  • 1990: Manasse, McGeoch and Sleator introduce the kkk-server problem and prove the deterministic two-server ratio is 222.
  • 1990–1994: Karlin, Manasse, McGeoch and Owicki submit this paper (received August 1990, revised September 1991) and publish it in Algorithmica in 1994, with the isosceles-triangle ratios of Theorem 12 and the 3-4-5 triangle ratio 1652/10691652/10691652/1069 of Theorem 13.
  • Later: Karloff, Rabani and Ravid extend the technique to Ω(log⁡log⁡k)\Omega(\log\log k)Ω(loglogk) and Ω(log⁡k)\Omega(\log k)Ω(logk) randomized lower bounds (cited on p. 564); Bubeck, Coester and Rabani (STOC 2023) refute the randomized kkk-server conjecture.

Setting

Fix an integer d≥1d\ge1d≥1. The isosceles triangle MMM has three points aaa, bbb, ccc with

dist⁡(a,b)=1,dist⁡(a,c)=dist⁡(b,c)=d.\operatorname{dist}(a,b)=1,\qquad \operatorname{dist}(a,c)=\operatorname{dist}(b,c)=d.dist(a,b)=1,dist(a,c)=dist(b,c)=d.

A configuration C:{0,1}→MC:\{0,1\}\to MC:{0,1}→M places two labelled servers on points of MMM. A deterministic on-line algorithm assigns to every finite request sequence σ=(r1,…,rn)\sigma=(r_1,\dots,r_n)σ=(r1​,…,rn​) a configuration, computed from σ\sigmaσ alone and covering the last request; its value on the empty sequence is its initial configuration. Its cost CA(σ)C_A(\sigma)CA​(σ) is the total distance its servers move while serving σ\sigmaσ request by request. The off-line optimum Copt(σ)C_{opt}(\sigma)Copt​(σ) from an initial configuration C0C_0C0​ is the least total movement of any schedule that starts at C0C_0C0​ and covers each request in turn, knowing σ\sigmaσ in advance.

A randomized algorithm is a probability distribution on deterministic on-line algorithms; its expected cost is ECA(σ)\mathbf{E}C_A(\sigma)ECA​(σ). It is ρ\rhoρ-competitive against an oblivious adversary from C0C_0C0​ if every algorithm in its support starts at C0C_0C0​ and there is a constant aaa such that

ECA(σ)≤ρ⋅Copt(σ)+afor every request sequence σ.\mathbf{E}C_A(\sigma)\le\rho\cdot C_{opt}(\sigma)+a\qquad\text{for every request sequence }\sigma.ECA​(σ)≤ρ⋅Copt​(σ)+afor every request sequence σ.

The request sequence is fixed in advance and does not react to the algorithm's coin flips.

Write ep=(1+1/p)pe_p=(1+1/p)^pep​=(1+1/p)p and

αd=e2d−1+1/4d(e2d−1−1)+1/2d,e2d−1=(2d2d−1)2d−1.\alpha_d=\frac{e_{2d-1}+1/4d}{(e_{2d-1}-1)+1/2d},\qquad e_{2d-1}=\left(\frac{2d}{2d-1}\right)^{2d-1}.αd​=(e2d−1​−1)+1/2de2d−1​+1/4d​,e2d−1​=(2d−12d​)2d−1.

In Lean this is NonuniformCompetitive.Isosceles.isoscelesRatio d.

Formalization targets

Goal: Theorem 12

For every d≥1d\ge1d≥1 and every initial configuration C0C_0C0​:

∀A, ∀ρ,A is ρ-competitive from C0 ⟹ ρ≥αd,\forall A,\ \forall\rho,\quad A\text{ is }\rho\text{-competitive from }C_0\ \Longrightarrow\ \rho\ge\alpha_d,∀A, ∀ρ,A is ρ-competitive from C0​ ⟹ ρ≥αd​, ∃A: A is αd-competitive from C0.\exists A:\ A\text{ is }\alpha_d\text{-competitive from }C_0.∃A: A is αd​-competitive from C0​.

The two claims are also milestones of their own (no_better_ratio, ratio_attained).

The phase LP (§5, pp. 565–566)

For free real π1,…,π2d−1\pi_1,\dots,\pi_{2d-1}π1​,…,π2d−1​ and real α\alphaα with

(πk)2d+∑i=1k(1−πi)≤αk  (1≤k<2d),2d+∑i=12d−1(1−πi)+12≤α⋅2d,(\pi_k)2d+\sum_{i=1}^k(1-\pi_i)\le\alpha k\ \ (1\le k<2d),\qquad 2d+\sum_{i=1}^{2d-1}(1-\pi_i)+\tfrac12\le\alpha\cdot2d,(πk​)2d+i=1∑k​(1−πi​)≤αk  (1≤k<2d),2d+i=1∑2d−1​(1−πi​)+21​≤α⋅2d,

one has α≥αd\alpha\ge\alpha_dα≥αd​ (lp_lower_bound); and πk=(αd−1)((2d/(2d−1))k−1)\pi_k=(\alpha_d-1)\big((2d/(2d-1))^k-1\big)πk​=(αd​−1)((2d/(2d−1))k−1), π2d=1\pi_{2d}=1π2d​=1 is nondecreasing from π1≥0\pi_1\ge0π1​≥0 to 111 and makes every constraint an equality (lp_attained).

The limit remark (§5, p. 566)

α1<α2<α3<⋯ ,lim⁡d→∞αd=ee−1\alpha_1<\alpha_2<\alpha_3<\cdots,\qquad \lim_{d\to\infty}\alpha_d=\frac{e}{e-1}α1​<α2​<α3​<⋯,d→∞lim​αd​=e−1e​

(ratio_increases_to_e_ratio).

Significance

The theorem gives an exact optimal randomized ratio for an infinite family of metric spaces. It shows that the optimal randomized two-server ratio is not a constant: it runs from 3/23/23/2 on the equilateral triangle to e/(e−1)≈1.582e/(e-1)\approx1.582e/(e−1)≈1.582 as the triangle becomes long and thin, where the problem resembles ski rental. With the deterministic ratio 222, it quantifies exactly how much randomization gains on these spaces.

The results are proved in the paper; none is formalized on Prove2Me, and no machine-checked proof of them is known. A formal proof would require the paper's phase framework (Theorems 1–3 and the appendix's Theorem 15) for server problems, which this mission does not state separately, and a concrete randomized algorithm as a measurable mixed strategy. Both would be reusable for Theorem 13 (the 3-4-5 triangle) and for other exact ratios on small metric spaces.

Difficulty

The phase LP milestones are finite real arithmetic. The difficulty is the passage between them and the goal. The lower bound must hold for every randomized algorithm, not only phase-based lazy ones: an arbitrary algorithm may condition on the whole history, move non-lazily, and randomize in ways that do not reduce to the probabilities πk\pi_kπk​. The paper handles this with Theorem 3, which says that the LP bound of phase-based algorithms bounds the competitive factor of all algorithms; its proof uses an averaging argument over histories that must be made rigorous. The upper bound needs a mixed strategy over infinitely many phases, with measurable costs, an explicit additive constant covering the first partial phase from an arbitrary initial configuration, and an accounting of CoptC_{opt}Copt​ across phase boundaries.

Formalization scope

The model is the platform's published KServer_model and KServer_randomized (reference items): labelled servers Fin 2 → M; a deterministic on-line algorithm as a map from request prefixes to configurations; a randomized algorithm as a probability measure over deterministic algorithms, with the cost of each fixed sequence measurable in the random outcome; expected cost as a lower Lebesgue integral in [0,∞][0,\infty][0,∞]; the off-line optimum as a real infimum over schedules from C0C_0C0​ (nonempty and bounded below by 000); and IsCompetitiveFrom A C₀ c with a real additive constant.

Committed conventions:

  • The triangle is any metric space whose points are exactly a,b,ca,b,ca,b,c at distances 1,d,d1,d,d1,d,d, with ddd a natural number and d≥1d\ge1d≥1. Every such space is isometric to the paper's triangle; at d=0d=0d=0 it would not be a triangle.
  • Both claims are stated for every initial configuration, including both servers on one point. The paper treats the initial state {a,b}\{a,b\}{a,b} separately and absorbs the first partial phase into the additive constant.
  • The lower bound quantifies over all randomized algorithms (deterministic ones are point masses), never over phase-based ones only.
  • In the LP milestones the πk\pi_kπk​ are free reals, as printed; no box 0≤πk≤10\le\pi_k\le10≤πk​≤1 is imposed.
  • "Grows" in the limit remark is read as strictly increasing.
  • The paper prints the recurrence on p. 565 as πk=α−1+(πk−1)2d−12d\pi_k=\frac{\alpha-1+(\pi_{k-1})2d-1}{2d}πk​=2dα−1+(πk−1​)2d−1​; the equations (∗)(*)(∗) give πk=α−1+2d πk−12d−1\pi_k=\frac{\alpha-1+2d\,\pi_{k-1}}{2d-1}πk​=2d−1α−1+2dπk−1​​. The recurrence is not used; the closed form printed on p. 566 is correct and is the one stated.

Without the measurability field of a randomized algorithm the lower integral would under-report expected cost and the attainment claim would become easier than the paper's; the published definition includes it. The lower bound is not vacuous: the triangle hypotheses are satisfiable for every d≥1d\ge1d≥1.

Welcome contributions: a formal version of the phase framework (Theorems 1–3, 15) for finite metric spaces, reusable across missions III and IV; a measurable construction of phase-based randomized algorithms; and proofs of the LP milestones.

Selected references

  • A. R. Karlin, M. S. Manasse, L. A. McGeoch, S. Owicki, Competitive Randomized Algorithms for Nonuniform Problems, Algorithmica 11 (1994) 542–571. https://doi.org/10.1007/BF01189993
  • M. S. Manasse, L. A. McGeoch, D. D. Sleator, Competitive Algorithms for Server Problems, J. Algorithms 11 (1990) 208–230. https://doi.org/10.1016/0196-6774(90)90003-W
  • H. Karloff, Y. Rabani, Y. Ravid, Lower Bounds for Randomized k-Server and Motion-Planning Algorithms, SIAM J. Comput. 23 (1994) 293–312. https://doi.org/10.1137/S0097539792224838
  • S. Bubeck, C. Coester, Y. Rabani, The Randomized k-Server Conjecture Is False!, STOC 2023. https://arxiv.org/abs/2211.05753
9 thms5 active usersReviewed
🏆Completed
Captain: xuanji

AlphaEvolve Eighth-Power Bound: omega < 2.371177Research Paper

Formalize ω<2.371177\omega < 2.371177ω<2.371177, the current record from Dupont, Eisenberger, Kozlovskii, Mehrabian, Ruiz, See, Zhou, Alman, Vassilevska Williams and Balog (arXiv:2608.16884, August 2026).

The paper applies the combination-loss laser method of Alman et al. (arXiv:2404.16349), formalized here as the ω<2.37134\omega < 2.37134ω<2.37134 entry, to CW5⊗8CW_5^{\otimes 8}CW5⊗8​. That is level ℓ∗=4\ell^* = 4ℓ∗=4 instead of 3, with an exact rational certificate of about 7⋅1067\cdot 10^67⋅106 parameters found by gradient-based optimization and AlphaEvolve. The authors state that the solution and verification code are being prepared for release.

9 thms5 active usersReviewed
🏆Completed
CombinatoricsOperations Research·Captain: mikedeng1

Competitive Paging Algorithms IV: An Algorithm Competitive against Several Others Exists iff the Reciprocal Ratios Sum to at Most 1Research Paper

Motivation

Paging is the problem of managing a fast memory that holds kkk pages out of nnn: when a requested page is not in fast memory (a page fault), some resident page must be evicted, and the cost of an algorithm is its number of faults. Practitioners have many eviction rules. Least-recently-used (LRU) performs well on real workloads but can be kkk times worse than the optimal off-line schedule; the randomized marking algorithm of the same paper is 2Hk2H_k2Hk​-competitive and so has better worst-case guarantees. Fiat, Karp, Luby, McGeoch, Sleator and Young asked in 1991 whether one on-line algorithm can combine the advantages of several given ones, and answered the question exactly: the attainable combinations of ratios are characterized by one inequality (arXiv:cs/0205038, §6).

The question of combining on-line algorithms has since become a theme of its own: combining heuristics with worst-case-safe algorithms, and, more recently, combining machine-learned predictions with robust fallbacks, both ask for the same kind of guarantee against several reference algorithms at once.

Setting

A type (k,n)(k,n)(k,n) consists of kkk servers and a finite set MMM of nnn vertices with the uniform metric: two distinct vertices are at distance 111. This is paging: vertices are pages, the vertices covered by servers are the pages in fast memory, and a server move is a page fault.

A deterministic on-line algorithm AAA of type (k,n)(k,n)(k,n) has an initial configuration of its kkk servers and, after each request r∈Mr\in Mr∈M, moves servers so that some server covers rrr; its configuration after a request sequence depends only on that sequence. Its cost CA(σ)C_A(\sigma)CA​(σ) on a request sequence σ\sigmaσ is the total distance its servers travel, i.e. the number of server moves.

For algorithms AAA and BBB of the same type and a constant ccc, AAA is ccc-competitive against BBB if there is a constant aaa such that for every request sequence σ\sigmaσ

CA(σ)≤c⋅CB(σ)+a.C_A(\sigma)\le c\cdot C_B(\sigma)+a .CA​(σ)≤c⋅CB​(σ)+a.

A sequence c∗=(c(1),…,c(m))c^*=(c(1),\dots,c(m))c∗=(c(1),…,c(m)) of positive reals is realizable if for every type (k,n)(k,n)(k,n) and every mmm deterministic on-line algorithms B(1),…,B(m)B(1),\dots,B(m)B(1),…,B(m) of that type there is a deterministic on-line algorithm AAA of the same type that is c(i)c(i)c(i)-competitive against B(i)B(i)B(i) for every iii.

Formalization targets

Goal: Theorem 6

For m≥1m\ge1m≥1 and positive reals c(1),…,c(m)c(1),\dots,c(m)c(1),…,c(m),

c∗ is realizable  ⟺  ∑1≤i≤m1c(i)≤1.c^*\ \text{is realizable}\iff \sum_{1\le i\le m}\frac1{c(i)}\le 1 .c∗ is realizable⟺1≤i≤m∑​c(i)1​≤1.

Milestones

In the order of the paper's proof:

  1. Punishments are paid for. If AAA punishes BBB at a time step (an AAA-interval on a vertex vvv ends at that step and contains the end of a BBB-interval on vvv that began no later), then BBB has moved a server; the number of such steps is at most CB(σ)C_B(\sigma)CB​(σ).
  2. A fault leaves room to punish. If ∣SA∣=k|S_A|=k∣SA​∣=k, ∣SB∣≤k|S_B|\le k∣SB​∣≤k, x∈SBx\in S_Bx∈SB​ and x∉SAx\notin S_Ax∈/SA​, then some u∈SAu\in S_Au∈SA​ is not in SBS_BSB​.
  3. The greedy quota claim. If ∑i1/c(i)≤1\sum_i 1/c(i)\le 1∑i​1/c(i)≤1 and each unit of cost punishes the B(i)B(i)B(i) minimizing c(i)(PUN(i)+1)c(i)(\mathrm{PUN}(i)+1)c(i)(PUN(i)+1) (other algorithms may be punished incidentally), then after cost rrr every B(i)B(i)B(i) has been punished at least ⌊r/c(i)⌋\lfloor r/c(i)\rfloor⌊r/c(i)⌋ times.
  4. Shuttle algorithms. With 2m−12m-12m−1 servers on 2m2m2m vertices there are mmm algorithms, each keeping all vertices outside its own pair covered, no two of which move at the same step; in particular their total cost on any σ\sigmaσ is at most ∣σ∣|\sigma|∣σ∣.
  5. A forcing adversary. With 2m−12m-12m−1 servers on 2m2m2m vertices every algorithm can be forced to move at each of NNN steps, so CA(τ(N))≥NC_A(\tau(N))\ge NCA​(τ(N))≥N.

Significance

The result. Theorem 6 is an exact characterization, not a bound: the region of simultaneously attainable ratios against arbitrary deterministic paging algorithms is {c:∑1/c(i)≤1}\{c:\sum 1/c(i)\le 1\}{c:∑1/c(i)≤1}. For example, any two paging algorithms can be combined into one that is 222-competitive against each, and no better symmetric pair is possible in general. Combined with Theorem 7 of the same paper (not part of this mission), the same region is attainable against randomized algorithms, which is how LRU's practical behaviour and the marking algorithm's 2Hk2H_k2Hk​ worst-case guarantee can be obtained within constant factors by one algorithm.

Formalizing it. The theorem has been proved since 1991; no machine-checked proof is known. A formal proof produces a reusable notion of competitiveness of one on-line algorithm against another, built on the published KServer_model definitions, and a formal account of the scheduling fact at the core of the sufficiency proof.

Difficulty

Sufficiency looks like an averaging argument, but the combined algorithm cannot simulate the B(i)B(i)B(i) and follow one of them: switching between their configurations costs up to kkk per switch, which no additive constant absorbs. The accounting has to charge each of AAA's faults to a specific move of a specific B(i)B(i)B(i), and the charge must be injective; the paper's claim that CB(σ)C_B(\sigma)CB​(σ) is at least the number of punishments is where this happens, and it depends on how server intervals are matched. The allocation of faults to algorithms is then a deadline-scheduling problem whose feasibility is exactly ∑1/c(i)≤1\sum 1/c(i)\le 1∑1/c(i)≤1, and the floor functions make the counting delicate at the boundary. The paper's own definition of punishment only counts intervals that start with a move, so the first kkk faults of AAA (servers on their initial vertices) need separate treatment; they are absorbed by the additive constant.

Necessity needs the right family of hard instances: the mmm algorithms must never move at the same step, which pins the type to (2m−1,2m)(2m-1,2m)(2m−1,2m).

Formalization scope

The Lean development works in the namespace CompetitivePaging.Combining and imports the published KServer_model definitions: KServer.OnlineAlgorithm k M (a configuration map from request prefixes to Fin k → M with a serving condition) and OnlineAlgorithm.cost. Committed conventions:

  • a type (k,n)(k,n)(k,n) is any k : ℕ and any finite M : Type with a metric in which distinct points are at distance 111; realizability quantifies over all of them, never over one fixed type;
  • servers are labelled; each algorithm has its own initial configuration, and the additive constant aaa is chosen before the request sequence;
  • c(i)>0c(i)>0c(i)>0 and m≥1m\ge1m≥1 are hypotheses of the goal, as in the paper; without positivity, 1/0=01/0=01/0=0 in Lean would make a zero ratio free;
  • time ttt is the step processing the ttt-th request; the paper's PUN\mathrm{PUN}PUN counts time steps.

Trivializing encodings are ruled out: realizability is not stated for a single fixed type, the metric is not the metric of Fin n, and the competitive constant is not allowed to depend on the request sequence.

A complete proof needs the construction of the punishing algorithm as a KServer.OnlineAlgorithm (a lazy, injective algorithm whose moves depend on the prefix and on the B(i)B(i)B(i)'s configurations), the injective charging argument, the scheduling lemma, and the explicit shuttle algorithms. The scheduling lemma and the charging lemma are independent of paging and reusable. Proofs of any milestone are welcome, as are alternative statements of the sufficiency construction.

Selected references

  • A. Fiat, R. M. Karp, M. Luby, L. A. McGeoch, D. D. Sleator, N. E. Young, Competitive Paging Algorithms, J. Algorithms 12(4):685–699, 1991. doi:10.1016/0196-6774(91)90041-V; preprint arXiv:cs/0205038.
  • D. D. Sleator, R. E. Tarjan, Amortized efficiency of list update and paging rules, Comm. ACM 28(2):202–208, 1985. doi:10.1145/2786.2793
  • M. S. Manasse, L. A. McGeoch, D. D. Sleator, Competitive algorithms for server problems, J. Algorithms 11(2):208–230, 1990. doi:10.1016/0196-6774(90)90003-W
9 thms5 active usersReviewed
🏆Completed
CombinatoricsComplexity TheoryOperations Research·Captain: mikedeng1

Scheduling Subject to Resource Constraints: Classification and Complexity I: Unit-Time Chains on Two Identical Machines with One Unit Resource Are Strongly NP-hardResearch Paper

Resource constraints and the easy/hard borderline in scheduling

Machine scheduling asks how to assign jobs to machines over time so that a criterion such as the makespan Cmax⁡C_{\max}Cmax​, the time at which the last job completes, is as small as possible. In practice jobs also compete for scarce resources beyond the machines themselves: tools, operators, memory, power. Błażewicz, Lenstra and Rinnooy Kan (DAM 1983) extended the standard three-field classification α∣β∣γ\alpha\mid\beta\mid\gammaα∣β∣γ of Graham, Lawler, Lenstra and Rinnooy Kan (1979) by a resource field resλσρres\lambda\sigma\rhoresλσρ. They then settled the complexity of every problem with parallel identical or uniform machines, unit-time jobs, precedence constraints and the Cmax⁡C_{\max}Cmax​ criterion. Their Fig. 2 separates the maximal polynomially solvable problems from the minimal NP-hard ones, and it has been the reference map for resource-constrained scheduling since.

Brief timeline of the problems involved:

  • 1975. Garey and Johnson (SIAM J. Comput. 4) show that P2∣res⋯ ,pj=1∣Cmax⁡P2\mid res\cdots, p_j=1\mid C_{\max}P2∣res⋯,pj​=1∣Cmax​ is solvable in polynomial time via matchings, and that P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ and P2∣res1⋅⋅,tree,pj=1∣Cmax⁡P2\mid res1\cdot\cdot, tree, p_j=1\mid C_{\max}P2∣res1⋅⋅,tree,pj​=1∣Cmax​ are NP-hard in the strong sense, by reduction from 3-PARTITION.
  • 1976. Ullman (Complexity of sequencing problems, in Coffman, ed., Computer & Job/Shop Scheduling Theory, Wiley) gives strong NP-hardness of P2∣res111,prec,pj=1∣Cmax⁡P2\mid res111, prec, p_j=1\mid C_{\max}P2∣res111,prec,pj​=1∣Cmax​ under arbitrary precedence constraints.
  • 1983. Błażewicz, Lenstra and Rinnooy Kan prove Theorem 7: chains suffice. Two identical machines, one resource of size one, requirements in {0,1}\{0,1\}{0,1} and chain-like precedence already give a strongly NP-hard problem. The result dominates both earlier two-machine results.

Setting

There are nnn jobs J1,…,JnJ_1,\dots,J_nJ1​,…,Jn​ and mmm machines M1,…,MmM_1,\dots,M_mM1​,…,Mm​. Every job has processing time 111 on every machine, each machine handles at most one job at a time, and jobs are not preempted. There are lll resources; resource RhR_hRh​ has a positive integer size shs_hsh​, the amount available at any time, and job JjJ_jJj​ has a nonnegative integer requirement rhjr_{hj}rhj​, the amount it holds throughout its execution. A directed acyclic graph HHH on the jobs gives the precedence constraints: if HHH has a path from jjj to kkk (Jj→JkJ_j\to J_kJj​→Jk​), then JjJ_jJj​ must complete before JkJ_kJk​ starts. The precedence is chain-like when every vertex of HHH has indegree and outdegree at most one.

A schedule gives each job a machine and a real start time SjS_jSj​; the job occupies [Sj,Sj+1)[S_j, S_j+1)[Sj​,Sj​+1) and completes at Cj=Sj+1C_j = S_j+1Cj​=Sj​+1. It is feasible if jobs on one machine do not overlap, precedence is respected, and at every time ttt the jobs running at ttt require at most shs_hsh​ of each resource RhR_hRh​. The makespan is Cmax⁡=max⁡jCjC_{\max} = \max_j C_jCmax​=maxj​Cj​.

The problem P2∣res111,chain,pj=1∣Cmax⁡P2\mid res111, chain, p_j=1\mid C_{\max}P2∣res111,chain,pj​=1∣Cmax​ restricts this to m=2m=2m=2, one resource (λ=1\lambda=1λ=1) of size 111 (σ=1\sigma=1σ=1), every requirement at most 111 (ρ=1\rho=1ρ=1), and chain-like precedence. The problem P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ has m=3m=3m=3, one resource of arbitrary size and requirements, and no precedence.

3-PARTITION: given ttt, a positive integer bbb and positive integers a1,…,a3ta_1,\dots,a_{3t}a1​,…,a3t​ with ∑jaj=tb\sum_j a_j = tb∑j​aj​=tb and 14b<aj<12b\tfrac14 b<a_j<\tfrac12 b41​b<aj​<21​b, can {1,…,3t}\{1,\dots,3t\}{1,…,3t} be split into ttt disjoint 3-element sets SiS_iSi​ with ∑j∈Siaj=b\sum_{j\in S_i}a_j=b∑j∈Si​​aj​=b?

A problem is NP-hard in the strong sense if it remains NP-hard when every number of the instance is written in unary.

Formalization targets

Goal: Theorem 7

3-PARTITION is NP-hard in the strong sense  ⟹  P2∣res111, chain, pj=1∣Cmax⁡ is NP-hard in the strong sense.\text{3-PARTITION is NP-hard in the strong sense} \;\Longrightarrow\; P2\mid res111,\ chain,\ p_j=1\mid C_{\max}\ \text{is NP-hard in the strong sense.}3-PARTITION is NP-hard in the strong sense⟹P2∣res111, chain, pj​=1∣Cmax​ is NP-hard in the strong sense.

The hypothesis is Garey and Johnson's theorem on 3-PARTITION, which the paper cites and does not prove. The conclusion concerns the decision version: given an instance and y∈Ny\in\mathbb Ny∈N, is there a feasible schedule with Cmax⁡≤yC_{\max}\le yCmax​≤y?

Milestones

  1. Proof of Theorem 4, the saturation equivalence. For positive bbb, aja_jaj​ with ∑jaj=tb\sum_j a_j=tb∑j​aj​=tb, the P3∣res1⋅⋅P3\mid res1\cdot\cdotP3∣res1⋅⋅ instance with 3t3t3t unit jobs, resource size bbb and requirements aja_jaj​ has a feasible schedule with Cmax⁡≤tC_{\max}\le tCmax​≤t iff the 3-PARTITION instance has a solution.
  2. Theorem 4 (Garey and Johnson). Under the same hypothesis as the goal, P3∣res1⋅⋅,pj=1∣Cmax⁡P3\mid res1\cdot\cdot, p_j=1\mid C_{\max}P3∣res1⋅⋅,pj​=1∣Cmax​ is NP-hard in the strong sense.
  3. Proof of Theorem 7, "if". A 3-PARTITION solution yields a feasible schedule of the constructed two-machine instance with Cmax⁡=2tbC_{\max}=2tbCmax​=2tb.
  4. Proof of Theorem 7, "only if". A feasible schedule of the constructed instance with Cmax⁡≤2tbC_{\max}\le 2tbCmax​≤2tb yields a 3-PARTITION solution.

Significance

Theorem 7 is the sharpest hardness result of the paper's classification. Without resources, two-machine unit-time scheduling with arbitrary precedence is polynomial (Coffman and Graham, Acta Inform. 1972); without precedence, it is polynomial under arbitrary resources (Theorem 1 of the paper). The theorem shows that combining the weakest nontrivial versions of both constraints, chains and one unit resource, already crosses the borderline. The paper's §4.1 extends the same reduction to the ∑Cj\sum C_j∑Cj​ and Lmax⁡L_{\max}Lmax​ criteria.

On the formal side, the mission provides a machine-checked model of resource-constrained scheduling with real start times, a definition of NP-hardness in the strong sense on top of the platform's Turing-machine formalization of P\mathrm PP and NP\mathrm{NP}NP, and 3-PARTITION as a reusable source problem. As far as the platform's corpus shows, none of Theorems 4 and 7, 3-PARTITION, or strong NP-hardness has been formalized before. Both theorems are proved in the literature; what remains is to formalize the reductions and their polynomial running time.

Difficulty

The combinatorial heart is the "only if" direction: a schedule of length 2tb2tb2tb must be shown to be rigid. Start times are arbitrary reals, so the first obstacle is to show that both machines are busy throughout [0,2tb)[0,2tb)[0,2tb), that the chain LLL forces unit spacing, and that the primed jobs of the chains Kj′K'_jKj′​ can only run in the intervals the chain LLL leaves free of the resource. Only after this is established can the index sets SiS_iSi​ be read off. Arguing on integer time slots from the start is not enough: the model allows fractional start times, and ruling them out is part of the proof.

The second obstacle is the complexity layer. NP-hardness is stated with respect to polynomial-time many-one reductions computed by one-tape Turing machines. The reduction from 3-PARTITION therefore has to be implemented and its running time bounded on unary codes. The constructed instance has 4tb4tb4tb jobs, which is polynomial in the unary length of the 3-PARTITION instance; this is exactly why the reduction proves hardness in the strong sense.

Formalization scope

  • Model. Jobs are Fin n and machines Fin m, 0-based. Only identical machines with unit processing times are modelled. Start times are real, execution intervals are half-open, and the resource constraint is imposed at every real time. Precedence is the transitive closure of the arc list of HHH. Cmax⁡=0C_{\max}=0Cmax​=0 for an empty instance.
  • Decision version. Thresholds yyy are natural numbers; this narrower class makes the hardness statement stronger.
  • Encoding. An instance is described by its list of numbers (n,m,ln,m,ln,m,l, the sizes, the requirements row by row, the number of arcs and the arcs, then yyy). The unary language is the set of unary codes of yes-instances over a two-letter alphabet. No pairing function is used. The class conditions (two machines, one unit resource, requirements at most one, chain-like acyclic HHH) are part of the yes-predicate.
  • Strong sense. Strong NP-hardness is NP-hardness of the unary language. This is equivalent to Garey and Johnson's definition, which bounds the largest number by a polynomial in the instance length.
  • Cited hypothesis. The goal and Theorem 4 assume strong NP-hardness of 3-PARTITION (with 14b<aj<12b\tfrac14 b<a_j<\tfrac12 b41​b<aj​<21​b) and nothing else. Stating the goal as a bare reduction between the two languages, or adding P≠NP\mathrm P\ne\mathrm{NP}P=NP, would not be Theorem 7.
  • Constructions. The two scheduling instances built from a 3-PARTITION instance are explicit definitions following the page, not arbitrary instances with a property.
  • Reuse. The scheduling model and the strong-NP-hardness layer are shared with the other missions of this series; 3-PARTITION serves any strong NP-hardness proof by number partitioning.

Welcome contributions: proofs of the four milestones; a formalized polynomial-time implementation of the reduction on unary codes; general lemmas about composing polynomial-time reductions on the one-tape machine model.

Selected references

  • J. Błażewicz, J. K. Lenstra, A. H. G. Rinnooy Kan, Scheduling subject to resource constraints: classification and complexity, Discrete Applied Mathematics 5 (1983) 11–24. https://doi.org/10.1016/0166-218X(83)90012-4
  • M. R. Garey, D. S. Johnson, Complexity results for multiprocessor scheduling under resource constraints, SIAM J. Comput. 4 (1975) 397–411. https://doi.org/10.1137/0204035
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, Freeman, 1979.
  • R. L. Graham, E. L. Lawler, J. K. Lenstra, A. H. G. Rinnooy Kan, Optimization and approximation in deterministic sequencing and scheduling: a survey, Ann. Discrete Math. 5 (1979) 287–326. https://doi.org/10.1016/S0167-5060(08)70356-X
  • J. D. Ullman, Complexity of sequencing problems, in: E. G. Coffman, Jr., ed., Computer & Job/Shop Scheduling Theory, Wiley, 1976, 139–164.
  • E. G. Coffman, Jr., R. L. Graham, Optimal scheduling for two-processor systems, Acta Informatica 1 (1972) 200–213. https://doi.org/10.1007/BF00288685
  • S. Cook, The P versus NP problem, Clay Mathematics Institute official problem description.
11 thms5 active usersReviewed
Computational GeometryLinear OptimizationOperations Research·Captain: mikedeng1

Linear Programming in Linear Time When the Dimension Is Fixed: Fixed-Dimension LP Feasibility Decided in Linear Time on the Real RAMResearch Paper

Motivation

A linear program asks for a point x∈Rdx\in\mathbb{R}^dx∈Rd minimizing cTxc^TxcTx subject to nnn linear inequalities ∑j=1daijxj≥bi\sum_{j=1}^d a_{ij}x_j\ge b_i∑j=1d​aij​xj​≥bi​. Many problems in computational geometry and statistics are linear programs with few variables and very many constraints: separating two point sets by a line or plane, fitting a line in the Chebyshev (L∞L_\inftyL∞​) norm, finding the smallest disk or ball containing a point set (a related convex problem). For these problems the number of variables ddd is a small constant, and what matters is how the running time grows with nnn.

Nimrod Megiddo showed that for every fixed ddd the problem can be solved in time C(d)⋅nC(d)\cdot nC(d)⋅n (J. ACM 31(1), 1984).

Timeline.

  • 1983. Megiddo (SIAM J. Comput. 12) and, independently, Dyer (SIAM J. Comput. 13 (1984)) give linear-time algorithms for d=2d=2d=2 and d=3d=3d=3.
  • 1984. Megiddo extends the method to every fixed ddd, with C(d)<22d+2C(d)<2^{2^{d+2}}C(d)<22d+2 (the paper formalized here).
  • 1988–1991. Clarkson (J. ACM 42 (1995), conference version 1988) gives a randomized algorithm with expected time O(d2n)+dO(d)log⁡nO(d^2n)+d^{O(\sqrt d)}\log nO(d2n)+dO(d​)logn. Seidel (Discrete Comput. Geom. 6 (1991)) gives a simple randomized O(d! n)O(d!\,n)O(d!n) algorithm.
  • 1992–1996. Matoušek, Sharir and Welzl and, independently, Kalai give subexponential randomized bounds. Chazelle and Matoušek derandomize the linear dependence with C(d)=dO(d)C(d)=d^{O(d)}C(d)=dO(d) (J. Algorithms 21 (1996)).

Setting

Fix ddd. An instance is a matrix A∈Rn×dA\in\mathbb{R}^{n\times d}A∈Rn×d and a vector b∈Rnb\in\mathbb{R}^nb∈Rn, and its feasible region is the polyhedron P(A,b)={x∈Rd:Ax≥b}P(A,b)=\{x\in\mathbb{R}^d: Ax\ge b\}P(A,b)={x∈Rd:Ax≥b}. Here nnn is the number of constraints and ddd the number of variables.

The model of computation is the real RAM. A program is a finite list of instructions acting on real registers, integer pointer registers and a memory Z→R\mathbb{Z}\to\mathbb{R}Z→R. It performs exact +,−,×,/+,-,\times,/+,−,×,/ on reals at unit cost, tests the sign of a real, sets, copies, increments, decrements and compares pointers, and loads and stores through pointers. The input is the standard encoding of (A,b)(A,b)(A,b) in memory: the numbers nnn and ddd, then AAA row by row, then bbb. A program decides an instance within TTT steps with output β∈{accept,reject}\beta\in\{\text{accept},\text{reject}\}β∈{accept,reject} if it halts on that output after at most TTT steps.

Megiddo's method rests on multidimensional search. There is an unknown point x∗∈Rdx^*\in\mathbb{R}^dx∗∈Rd and an oracle that, for any hyperplane {x:aTx=b}\{x: a^Tx=b\}{x:aTx=b}, answers whether aTx∗<ba^Tx^*<baTx∗<b, =b=b=b or >b>b>b. Given hyperplanes Hi={aiTx=bi}H_i=\{a_i^Tx=b_i\}Hi​={aiT​x=bi​} with ai≠0a_i\ne0ai​=0, the question is how many oracle calls determine the position of x∗x^*x∗ relative to all of them. A search strategy is a ternary decision tree: inner nodes are hyperplane queries, leaves carry outputs, and the tree is built from the data alone. For linear programming, x∗x^*x∗ is an optimal solution, or a minimizer of the infeasibility function f(x)=max⁡i(bi−aiTx)f(x)=\max_i(b_i-a_i^Tx)f(x)=maxi​(bi​−aiT​x) when the system is infeasible. The oracle is implemented by solving problems in d−1d-1d−1 variables.

Formalization targets

Goal: linear-time feasibility on the real RAM

∀d ∃R ∃C ∀n ∀A∈Rn×d, b∈Rn:R decides within C (n+1) steps whether {x:Ax≥b}≠∅.\forall d\ \exists R\ \exists C\ \forall n\ \forall A\in\mathbb{R}^{n\times d},\,b\in\mathbb{R}^n:\quad R\text{ decides within }C\,(n+1)\text{ steps whether } \{x: Ax\ge b\}\neq\emptyset.∀d ∃R ∃C ∀n ∀A∈Rn×d,b∈Rn:R decides within C(n+1) steps whether {x:Ax≥b}=∅.

The program and the constant depend on ddd only. No explicit form of C(d)C(d)C(d) is fixed.

Milestones

  1. One query settles half of nnn hyperplanes on the line (A(1)=1A(1)=1A(1)=1, B(1)=12B(1)=\tfrac12B(1)=21​).
  2. v(ϵ)=(1,ϵ,…,ϵd−1)v(\epsilon)=(1,\epsilon,\dots,\epsilon^{d-1})v(ϵ)=(1,ϵ,…,ϵd−1) is orthogonal to some aia_iai​ for at most n(d−1)n(d-1)n(d−1) values of ϵ\epsilonϵ, so there is a basis in which all aij≠0a_{ij}\ne0aij​=0.
  3. For hyperplanes of opposite slopes in the (x1,x2)(x_1,x_2)(x1​,x2​) plane, the answers for Hik(1)H^{(1)}_{ik}Hik(1)​ and Hik(2)H^{(2)}_{ik}Hik(2)​ settle one of HiH_iHi​, HkH_kHk​.
  4. A linearly dependent pair of opposite slopes has ai1=ak1=0a_{i1}=a_{k1}=0ai1​=ak1​=0, and the middle hyperplane settles one of them.
  5. Approach I: 2d−12^{d-1}2d−1 queries settle at least ⌊21−2dn⌋\lfloor 2^{1-2^d}n\rfloor⌊21−2dn⌋ hyperplanes.
  6. C(d)log⁡nC(d)\log nC(d)logn queries settle all nnn hyperplanes.
  7. If a hyperplane contains no optimal point, all optimal points lie on one side of it.
  8. The oracle, Case I: at an optimum relative to {xd=0}\{x_d=0\}{xd​=0}, two auxiliary systems decide the side or certify global optimality.
  9. The oracle, Case II: at a minimizer of fff on {xd=0}\{x_d=0\}{xd​=0}, systems (1) and (2) decide the side or certify infeasibility.

Significance

The result. For every fixed dimension, linear programming is solvable in time linear in the number of constraints. The algorithm is also strongly polynomial in fixed dimension: its operation count does not depend on the bit size of the data. Deciding whether the optimum is at most ttt is feasibility of Ax≥bAx\ge bAx≥b together with −cTx≥−t-c^Tx\ge-t−cTx≥−t, so the goal also covers the decision form of optimization. The prune-and-search technique of the paper, which discards a constant fraction of the constraints per round, became a standard tool of computational geometry.

Formalizing it. The result is proved and classical. The platform already has the cases d=1d=1d=1 (linear time) and d=2d=2d=2 (quadratic time, by Fourier–Motzkin elimination) on the same machine and input encoding (SmaleNinth.real_ram_decides_one_variable_lp_linear, SmaleNinth.real_ram_decides_two_variable_lp_quadratic). No machine-checked proof of the general statement is known. The work consists of the query-complexity layer (milestones 1–6), the convex-analytic correctness of the oracle (milestones 7–9), and a real-RAM implementation with a step count linear in nnn, including linear-time median selection. Alternative proofs, for example through Clarkson's or Seidel's algorithms made deterministic, are welcome for the goal.

Difficulty

The obvious approach is to find the optimum by testing constraints one by one or by eliminating variables. Fourier–Motzkin elimination produces Θ(n2)\Theta(n^2)Θ(n2) constraints after one step. Pivoting methods have no known bound linear in nnn. The key difficulty is to discard a constant fraction of the constraints using only a constant number of recursive calls in dimension d−1d-1d−1, when no single hyperplane test gives information about more than one constraint. The multidimensional search layer gives this, and it is where the pairing of hyperplanes by slope and the degenerate cases (dependent pairs, zero coefficients) have to be handled exactly. At the machine level, the step count must stay linear in nnn for a fixed program, so every median selection and every recursive call must be implemented within the budget, with the recursion depth depending on ddd only.

Formalization scope

  • Machine and input. The machine is the platform's real RAM SmaleNinth.RAMProgram with RAMDecidesInTime, and the input convention is SmaleNinth.encodeLP (published definitions, reused unchanged). No instruction is added: there is no LP, median, floor or sort primitive. Time is the number of machine steps.
  • Quantifier order. ∀d ∃R ∃C ∀n,A,b\forall d\ \exists R\ \exists C\ \forall n, A, b∀d ∃R ∃C ∀n,A,b. The bound is C(n+1)C(n+1)C(n+1) in the number nnn of constraints, so that the machine can halt at n=0n=0n=0. The paper's C(d)<22d+2C(d)<2^{2^{d+2}}C(d)<22d+2 counts unspecified units of "effort" with an unquantified θ(nd)\theta(nd)θ(nd) term, and it is not transferred to machine steps. Where a milestone's proof fixes a constant exactly, the constant is stated: 2d−12^{d-1}2d−1 queries and ⌊n/22d−1⌋\lfloor n/2^{2^d-1}\rfloor⌊n/22d−1⌋ settled hyperplanes in milestone 5.
  • Feasibility only. The machine outputs accept or reject. Returning an optimizer, "unbounded", or a minimizer of fff is not part of the goal. The case d=0d=0d=0 is included.
  • Query trees. Nodes are queries compare (a ⬝ᵥ x) b and nothing else, leaves hold fixed values, and correctness is required for every xxx. A tree over arbitrary tests of xxx would make milestones 5 and 6 empty, and it is excluded by the definition.
  • Indices. The paper's x1,x2x_1,x_2x1​,x2​ are indices 0, 1 of Fin (d + 2), and its xdx_dxd​ is Fin.last d of Fin (d + 1).
  • Corrections. Two passages of §4 are stated in corrected form. The Case I auxiliary objective includes the ±cd\pm c_d±cd​ term of the direction. In Case II, feasibility of (1) puts improvement in {xd>0}\{x_d>0\}{xd​>0}, where the page's last sentence says {xd<0}\{x_d<0\}{xd​<0}. The pairing claim carries ak1ai2−ak2ai1≠0a_{k1}a_{i2}-a_{k2}a_{i1}\ne0ak1​ai2​−ak2​ai1​=0, the hypothesis its argument uses, since linear independence alone does not give it.
  • Not included. Approach II and its bound O(n(log⁡n)d2)O(n(\log n)^{d^2})O(n(logn)d2), the remarks on slowly growing ddd, the randomized variants, and the applications of §1.
  • Reusable parts. The query-tree definition and milestones 1–6 apply to any prune-and-search problem with a hyperplane oracle. The oracle lemmas (7–9) are statements about convex piecewise-linear functions and polyhedra.

Selected references

  • N. Megiddo, Linear programming in linear time when the dimension is fixed, J. ACM 31(1):114–127, 1984. https://doi.org/10.1145/2422.322418
  • N. Megiddo, Linear-time algorithms for linear programming in R3R^3R3 and related problems, SIAM J. Comput. 12(4):759–776, 1983. https://doi.org/10.1137/0212052
  • M. E. Dyer, Linear time algorithms for two- and three-variable linear programs, SIAM J. Comput. 13(1):31–45, 1984. https://doi.org/10.1137/0213003
  • K. L. Clarkson, Las Vegas algorithms for linear and integer programming when the dimension is small, J. ACM 42(2):488–499, 1995. https://doi.org/10.1145/201019.201036
  • R. Seidel, Small-dimensional linear programming and convex hulls made easy, Discrete Comput. Geom. 6:423–434, 1991. https://doi.org/10.1007/BF02574699
  • B. Chazelle, J. Matoušek, On linear-time deterministic algorithms for optimization problems in fixed dimension, J. Algorithms 21(3):579–597, 1996. https://doi.org/10.1006/jagm.1996.0046
16 thms5 active usersReviewed
Machine LearningProbabilityStatistics·Captain: mikedeng1

Foundations of Machine Learning X: Regression and Rademacher Complexity BoundsTextbook

Motivation

Every generalization bound presented so far in this series is for classification, where the error of a prediction is binary (correct or not). Regression asks a different question: predictions are real-valued, and error is measured by the magnitude of the deviation from the true label, via a loss function L. Chapter 11 develops generalization theory for bounded regression, showing that the same two complexity measures used for classification — Rademacher complexity and a VC-dimension analogue — extend naturally, once the loss function itself is folded into the machinery via a Lipschitz-contraction argument (Rademacher route) or a reduction to classification via level-set thresholding (pseudo-dimension route).

Setting

A regression hypothesis h:X→ℝ is scored by a loss L:ℝ×ℝ→ℝ against a joint distribution D on X×ℝ (the stochastic scenario, since regression labels are rarely exactly reproducible); R(h) = E_{(x,y)~D}[L(h(x),y)] (Eq. 11.1) and R̂_S(h) = (1/m)∑L(h(x_i),y_i) (Eq. 11.2). For a finite hypothesis set, Theorem 11.1 gives a Hoeffding/union-bound guarantee directly, the regression analogue of chunk 02-pac's finite-hypothesis bound. For infinite H, §11.2.2 develops a Rademacher-complexity route: Proposition 11.2 shows that if L is µ-Lipschitz in its first (predicted-value) argument, the Rademacher complexity of the loss-composed family G = {(x,y)↦L(h(x),y) : h∈H} is controlled by µ times H's own Rademacher complexity, via Talagrand's contraction lemma (chunk 05-svm's Lemma 5.7); Theorem 11.3 combines this with chunk 03's Theorem 3.3 to give the chapter's headline bound. §11.2.3 develops an independent, purely combinatorial route: pseudo-dimension (Definition 11.5), a real-valued analogue of VC-dimension defined via threshold-witnessed shattering (Definition 11.4, restated via its own Eq. 11.3 as the VC-dimension of a thresholded indicator family); Theorem 11.8 gives a pseudo-dimension generalization bound by reducing regression to a family of classification problems (one per threshold t), using the tail-integral identity Eq. 11.5.

Formalization targets

Theorem 11.1 (milestone). For L bounded by M and H finite: for any δ>0, with probability at least 1-δ, for all h∈H: R(h) ≤ R̂_S(h) + M√((log|H|+log(1/δ))/(2m)).

Proposition 11.2 (milestone). For L non-negative, bounded by M, µ-Lipschitz in its first argument: for any sample S, R̂_S(G) ≤ µR̂_S(H).

Theorem 11.3 — the mission's goal. Under Proposition 11.2's hypotheses on L: for any δ>0, with probability at least 1-δ, for all h∈H: E[L(h(x),y)] ≤ (1/m)∑L(h(x_i),y_i) + 2µR_m(H) + M√(log(1/δ)/(2m)), and also with 2µR̂_S(H) + 3M√(log(2/δ)/(2m)).

Theorem 11.8 (milestone). For Pdim(G)=d, L non-negative bounded by M: for any δ>0, with probability at least 1-δ over a sample of size m, for all h∈H: R(h) ≤ R̂_S(h) + M√(2d log(em/d)/m) + M√(log(1/δ)/(2m)).

Significance

Theorem 11.3 is the chapter's own choice of headline result (§11.2's stated goal is to show "how the Rademacher complexity bounds of theorem 3.3 can be used to derive generalization bounds for regression"), and its proof genuinely reuses two pieces of prior machinery from this series — chunk 03's Theorem 3.3 and chunk 05's Talagrand's-lemma-style contraction — combined via a new observation (Proposition 11.2) specific to loss-composed families, not a restatement of either. Theorem 11.8 is the chapter's second, structurally independent technique: its em/d bound parallels chunk 03's Corollary 3.19 (both ultimately reduce to a VC-dimension-style growth-function argument), but the reduction itself — regression to a continuum of threshold classification problems, via the Lebesgue-integral tail identity Eq. (11.5) applied to |R(h)-R̂_S(h)| — is genuinely new content for this book, and pseudo-dimension has no prior art on the platform or in Mathlib. No prior art exists for this chapter's overall content either: GET /theorems?q=generalization%20bound%20regression and GET /theorems?q=pseudo-dimension both return zero hits.

Difficulty

Proposition 11.2's proof needs Talagrand's contraction lemma applied with the Lipschitz constant taken in the first argument of L only — the predicted value h(x_i), holding the true label y_i fixed — exactly the pitfall BRIEF.md names: a loss Lipschitz in the wrong argument, or in both arguments jointly, would not license this step. Theorem 11.8's proof is the chapter's most involved: it defines, for every h∈H and threshold t≥0, a classifier c(h,t):(x,y)↦1_{L(h(x),y)>t}, bounds |R(h)-R̂_S(h)| by M·sup_{t∈[0,M]}|R(c(h,t))- R̂_S(c(h,t))| via the tail-integral identity, and then applies a VC-dimension-style classification bound (Corollary 3.19) to the family of thresholded classifiers — whose VC-dimension is, by Eq. (11.3), exactly Pdim(G) by construction. A formalization that conflated pseudo-dimension with ordinary VC-dimension, or reused chunk 03's HasVCDim definition by relabeling, would misrepresent this chapter's genuinely different (real-valued, threshold-witnessed) combinatorial notion — precisely the pitfall BRIEF.md flags.

Formalization scope

Y := ℝ throughout (the book's own "Y a measurable subset of ℝ"), a harmless simplification consistent with every hypothesis, loss and Lipschitz condition in this chapter being stated for real-valued scores and labels. EmpiricalRademacherComplexity/ RademacherComplexity restate chunk 03-rademacher-vc's Definitions 3.1/3.2 locally, since a draft item cannot import another chunk's draft module. Shatters/PseudoDim are formalized via the book's own equivalent reformulation (Eq. 11.3, the thresholded-indicator form), rather than the sign-function form of Definition 11.4 directly, since the two coincide except at a measure-zero boundary the book itself does not address; PseudoDim mirrors chunk 03's HasVCDim Prop-valued pattern (does not cover Pdim(G)=+∞; every consuming theorem takes it as an explicit hypothesis) but is a structurally distinct definition built on Shatters, never a relabeling of HasVCDim, per BRIEF.md's pitfall note. Proposition 11.2's and Theorem 11.3's Lipschitz hypothesis (hLlip) is stated with the true label y' universally quantified outside the two-point comparison y1, y2 (the predicted values), matching "for any fixed y' ∈ Y, y ↦ L(y,y') is µ-Lipschitz" exactly — Lipschitzness in the first argument only, per BRIEF.md's pitfall note. RademacherComplexity (Measure.map Prod.fst D) H m gives the book's R_m(H) (H's Rademacher complexity under the marginal sampling distribution of the inputs x, i.e. D's first marginal). No numerical constant is altered from the book in any of the four theorems.

Not formalized: the L_p-loss worked example following Theorem 11.3's proof (an instantiation of the general theorem for a specific loss family, not a separate numbered theorem); Theorem 11.6 and Theorem 11.7 (worked pseudo-dimension examples for hyperplanes and vector spaces, background/illustration rather than the chapter's general machinery — drafting only these examples instead of the general Theorem 11.8 would be this chapter's trivializing formalization); the two-sided variant of Theorem 11.1 mentioned immediately after its proof (an unnumbered remark, not a separately displayed/numbered theorem); and all of §11.3 (linear regression, kernel ridge regression, SVR, Lasso and their online variants), which is applications-heavy per BRIEF.md's chapter restriction to §11.1-11.2.

Selected references

  • M. Mohri, A. Rostamizadeh, A. Talwalkar, Foundations of Machine Learning, 2nd ed., MIT Press, 2018, Chapter 11 (§11.1-11.2).
  • D. Haussler, "Decision theoretic generalizations of the PAC model for neural net and other learning applications," Information and Computation 100(1), 1992 (pseudo-dimension's origin).
  • D. Pollard, Convergence of Stochastic Processes, Springer, 1984 (the tail-integral identity Eq. 11.5's classical antecedent).
15 thms5 active users
🏆Completed
Machine LearningProbabilityStatistics·Captain: mikedeng1

Foundations of Machine Learning V: Kernel Methods and the Representer TheoremTextbook

Motivation

Linear methods like SVMs work only when the classes are linearly separable, but most real data is not. Chapter 6 shows how to get non-linear decision boundaries for free: replace the input space's inner product with a kernel KKK that implicitly computes an inner product in a (possibly very high- or infinite-dimensional) feature space, without ever explicitly computing the feature mapping. This works for any positive definite symmetric (PDS) kernel — and the chapter's central theorem shows that such a kernel always induces a genuine Hilbert space (the reproducing kernel Hilbert space, RKHS) in which the kernel is literally an inner product. The chapter's capstone, the representer theorem, then shows that a broad class of optimization problems over this (possibly infinite-dimensional) Hilbert space always has a solution expressible as a finite linear combination of kernel evaluations at the training points — turning an infinite-dimensional problem into a finite, mmm-dimensional one.

Setting

A kernel K:X×X→RK:X\times X\to\mathbb RK:X×X→R is PDS (Definition 6.3) if for every finite sample {x1,…,xm}⊆X\{x_1,\dots,x_m\}\subseteq X{x1​,…,xm​}⊆X, the Gram matrix [K(xi,xj)][K(x_i,x_j)][K(xi​,xj​)] is symmetric positive semidefinite. Theorem 6.8 shows every PDS kernel is an inner product K(x,x′)=⟨Φ(x),Φ(x′)⟩K(x,x')=\langle \Phi(x),\Phi(x')\rangleK(x,x′)=⟨Φ(x),Φ(x′)⟩ in some Hilbert space HHH (the RKHS), which further has the reproducing property h(x)=⟨h,K(x,⋅)⟩h(x)=\langle h,K(x,\cdot)\rangleh(x)=⟨h,K(x,⋅)⟩ for every h∈Hh\in Hh∈H — evaluating hhh at a point is itself an inner product with the kernel section at that point. Theorem 6.10 shows PDS kernels are closed under sum, product, tensor product, pointwise limit, and power-series composition, letting complex kernels (Gaussian, and many others) be built from simple ones (polynomial kernels) without re-verifying positive-semidefiniteness from scratch. Section 6.3's representer theorem (Theorem 6.11) then considers minimizing, over h∈Hh\in Hh∈H, an objective F(h)=G(∥h∥H)+L(h(x1),…,h(xm))F(h)=G(\|h\|_H)+L(h(x_1),\dots,h(x_m))F(h)=G(∥h∥H​)+L(h(x1​),…,h(xm​)) that depends on hhh only through its norm and its values at mmm fixed points.

Formalization targets

Theorem 6.8 (RKHS existence, milestone). For a PDS kernel KKK, there exist a Hilbert space HHH and Φ:X→H\Phi:X\to HΦ:X→H with K(x,x′)=⟨Φ(x),Φ(x′)⟩K(x,x')=\langle\Phi(x),\Phi(x')\rangleK(x,x′)=⟨Φ(x),Φ(x′)⟩, and HHH has the reproducing property h(x)=⟨h,K(x,⋅)⟩h(x)=\langle h,K(x,\cdot)\rangleh(x)=⟨h,K(x,⋅)⟩ for all h∈Hh\in Hh∈H, x∈Xx\in Xx∈X.

Theorem 6.10 (closure properties, milestone). PDS kernels are closed under sum, product, tensor product, pointwise limit, and power-series composition with non-negative coefficients.

Theorem 6.11 — the mission's goal. For any non-decreasing G:R→RG:\mathbb R\to\mathbb RG:R→R and any loss L:Rm→R∪{+∞}L:\mathbb R^m\to\mathbb R\cup\{+\infty\}L:Rm→R∪{+∞}, argminh∈HG(∥h∥H)+L(h(x1),…,h(xm))\mathrm{argmin}_{h\in H} G(\|h\|_H)+ L(h(x_1),\dots,h(x_m))argminh∈H​G(∥h∥H​)+L(h(x1​),…,h(xm​)) admits a solution h⋆=∑i=1mαiK(xi,⋅)h^\star=\sum_{i=1}^m\alpha_i K(x_i,\cdot)h⋆=∑i=1m​αi​K(xi​,⋅); if GGG is increasing, every solution has this form.

Significance

Theorem 6.11 is the chapter's payoff and one of the most widely used structural results in kernel methods: it explains, in one general statement covering SVMs, kernel ridge regression, Gaussian process MAP estimation and many other algorithms simultaneously, why the dual (finite, mmm-coefficient) formulation always suffices — the RKHS's infinite dimensionality never has to be confronted directly. Theorem 6.8 is the structural fact the whole chapter (and every later kernelized algorithm in the book, chapters 9-11, 15) depends on: without it, "PDS kernel" would be a purely combinatorial condition on Gram matrices with no guarantee it corresponds to any actual inner product. No prior art on the Prove2Me platform is faithful to any of this chapter's content: GET /theorems?q=Representer theorem and q=reproducing kernel return no faithful match (one unrelated hit concerns a Gaussian-measure reproducing kernel in a different, probabilistic context, not this chapter's PDS-kernel/RKHS construction). All six items are drafted fresh.

Difficulty

Theorem 6.8's proof is a genuine construction: define H0H_0H0​ as finite linear combinations of kernel sections Φ(x)=K(x,⋅)\Phi(x)=K(x,\cdot)Φ(x)=K(x,⋅), define an inner product on H0H_0H0​ using KKK itself, verify it is well-defined (independent of the representation), positive semidefinite (via the PDS hypothesis), and — via the Cauchy-Schwarz-for-PDS-kernels lemma (Lemma 6.7) — actually positive definite, then complete H0H_0H0​ to a genuine Hilbert space HHH in which it is dense, and finally extend the reproducing property from the dense subspace H0H_0H0​ to all of HHH by a continuity argument. This is substantial analysis, not a restatement. Theorem 6.11's proof uses the orthogonal decomposition H=H1⊕H1⊥H=H_1\oplus H_1^\perpH=H1​⊕H1⊥​ (where H1=span{K(xi,⋅)}H_1=\mathrm{span}\{K(x_i, \cdot)\}H1​=span{K(xi​,⋅)}) and the reproducing property to show the orthogonal component h⊥h^\perph⊥ never helps and, when GGG is strictly increasing, strictly hurts — a short argument, but one that depends essentially on Theorem 6.8's reproducing property holding for the specific HHH constructed, not just any Hilbert space with the kernel as its inner product.

Formalization scope

IsPDS uses the book's own second SPSD characterization (c^T K c ≥ 0 for every finite sample and coefficient vector c) rather than the non-negative-eigenvalues characterization, avoiding spectral theory for a Prop-valued definition; the book states the two are equivalent. IsRKHSOf and IsMinimizer are formalization scaffolding, not book-numbered definitions: IsRKHSOf packages Theorem 6.8's own two displayed equations (6.8, 6.9) as a reusable predicate shared between Theorem 6.8 (its conclusion) and Theorem 6.11 (its "H its corresponding RKHS" hypothesis), using an explicit evaluation map ev : H → X → ℝ to stand in for "elements of H are functions on X," since Mathlib's abstract Hilbert spaces are not themselves spaces of functions; IsMinimizer packages argmin. Both X in Theorem 6.8's existential and Theorem 6.11's ambient type are plain Type rather than Type*, avoiding universe-polymorphic quantification over the constructed Hilbert space's own type — a harmless simplification, since every application in this book instantiates X at a concrete, small type (typically RN\mathbb R^NRN or a finite set). Theorem 6.11's loss codomain ℝ ∪ {+∞} is WithTop ℝ, not EReal (which would also admit -∞, an unstated generalization the book's own display does not license, since EReal's ⊤+⊥=⊥ collapse is a genuine faithfulness risk the book's own L:\mathbb R^m\to\mathbb R\cup\{+\infty\} avoids by construction). WithTop ℝ on its own does not avoid every collapse, though: an unconstrained L may be the constant function ⊤ (a legal instance of ℝ∪\{+\infty\}), forcing the objective identically ⊤ and every point to vacuously minimize it, which would make the theorem's second conjunct false. An added hypothesis, ∃ h₀, F h₀ ≠ ⊤, makes explicit the book's own implicit assumption that the objective is finite somewhere — see the Formalization scope note below. Theorem 6.11's two clauses are otherwise kept exactly as distinct as the book states them: existence needs only Monotone G (non-decreasing); "any solution has this form" needs StrictMono G (increasing) as an added hypothesis on the second conjunct only — per this chunk's own BRIEF.md, the crux of the theorem, and the trap this mission is most careful to avoid collapsing. Theorem 6.10's five closure clauses are stated as one conjunction (matching the book's single theorem, not five separate items); the power-series clause keeps the book's own radius-of-convergence domain restriction and adds an explicit summability hypothesis guarding the ∑' term.

Not formalized: Theorem 6.2 (Mercer's condition) — not needed by the goal's own proof chain (it is an equivalent characterization of PDS mentioned before the RKHS construction, not a premise Theorem 6.8's or 6.11's proof invokes) and its own hypotheses (compact X⊂RNX\subset \mathbb R^NX⊂RN, continuous KKK, an eigenfunction expansion of a compact self-adjoint integral operator) are real analytic content this mission's budget does not include; Lemma 6.7 (Cauchy-Schwarz for PDS kernels) and Lemma 6.9 (normalized PDS kernels) — supporting lemmas for Theorem 6.8's proof, not independently numbered results the goal cites; Theorem 6.12/Corollary 6.13 (Rademacher complexity/margin bounds for kernel-based hypotheses) — the chapter's optional further milestone, connecting to chunks 03/05's machinery, cut for budget; §6.5-6.8 (sequence kernels, weighted transducers, rational kernels, Bochner's theorem, approximate feature maps) — explicitly out of scope per this chunk's own BRIEF.md, a distinct, applications-heavy topic.

Selected references

  • M. Mohri, A. Rostamizadeh, A. Talwalkar, Foundations of Machine Learning, 2nd ed., MIT Press, 2018, Chapter 6, §6.1-6.4.
  • B. Schölkopf, R. Herbrich, A. J. Smola, "A generalized representer theorem," COLT 2001, Lecture Notes in Computer Science 2111, 2001, 416-426.
  • N. Aronszajn, "Theory of reproducing kernels," Transactions of the American Mathematical Society 68(3), 1950, 337-404.
7 thms5 active users
🏆Completed
Captain: marwahaha

Coppersmith–Winograd Bound: omega < 2.376Research Paper

AI generated but i think correct. I think the milestones make it really annoying but the central theorem looks correct.

Motivation

The matrix-multiplication exponent measures the asymptotic number of field operations needed to multiply two square matrices. A bound ω<c\omega<cω<c means that, for every ε>0\varepsilon>0ε>0, two n×nn\times nn×n matrices can be multiplied using O(nc+ε)O(n^{c+\varepsilon})O(nc+ε) arithmetic operations. Matrix multiplication is a central benchmark in algebraic complexity and a primitive for many algorithms in linear algebra, graph theory, and symbolic computation.

After Strassen showed that ω<3\omega<3ω<3, a sequence of tensor constructions reduced the exponent further. Schönhage's asymptotic sum inequality made it possible to exploit simultaneous matrix products rather than a single square product. In 1990, Don Coppersmith and Shmuel Winograd combined an explicit low-border-rank tensor with a block extraction argument based on Salem--Spencer sets. Their basic analysis gave ω<2.38719\omega<2.38719ω<2.38719; coupling the random weights in the tensor square sharpened this to ω<2.375477\omega<2.375477ω<2.375477, hence the exact rational consequence ω<2.376\omega<2.376ω<2.376.

This mission formalizes that historical Coppersmith--Winograd result. It follows the source tensor and its actual block restrictions, while excluding placeholder “laser values” that are not backed by extracted direct sums of matrix-multiplication tensors.

Setting

For a field KKK, an order-three tensor is represented by three finite-dimensional KKK-vector spaces and an element of their tensor product. The matrix-multiplication tensor

⟨a,b,c⟩K=∑i<a∑j<b∑k<cxij⊗yjk⊗zki\langle a,b,c\rangle_K =\sum_{i<a}\sum_{j<b}\sum_{k<c} x_{ij}\otimes y_{jk}\otimes z_{ki}⟨a,b,c⟩K​=i<a∑​j<b∑​k<c∑​xij​⊗yjk​⊗zki​

encodes multiplication of an a×ba\times ba×b matrix by a b×cb\times cb×c matrix. A restriction applies one linear map to each tensor leg. A degeneration permits those maps to depend polynomially on a formal parameter and selects their first nonzero coefficient. Thus a degeneration from the diagonal tensor IrI_rIr​ is a border-rank certificate R‾(T)≤r\underline R(T)\le rR​(T)≤r.

The Coppersmith--Winograd tensor with parameter qqq is

Tq=∑i=1q(x0yizi+xiy0zi+xiyiz0)+x0y0zq+1+x0yq+1z0+xq+1y0z0.T_q= \sum_{i=1}^{q} (x_0y_i z_i+x_i y_0z_i+x_i y_i z_0) +x_0y_0z_{q+1}+x_0y_{q+1}z_0+x_{q+1}y_0z_0.Tq​=i=1∑q​(x0​yi​zi​+xi​y0​zi​+xi​yi​z0​)+x0​y0​zq+1​+x0​yq+1​z0​+xq+1​y0​z0​.

It has border rank at most q+2q+2q+2. Its coordinates carry three classes, indexed by 0,1,20,1,20,1,2, and its six nonzero block types are

(0,1,1), (1,0,1), (1,1,0), (0,0,2), (0,2,0), (2,0,0).(0,1,1),\ (1,0,1),\ (1,1,0),\ (0,0,2),\ (0,2,0),\ (2,0,0).(0,1,1), (1,0,1), (1,1,0), (0,0,2), (0,2,0), (2,0,0).

The first three blocks are matrix-multiplication tensors with dimensions (1,1,q)(1,1,q)(1,1,q), (q,1,1)(q,1,1)(q,1,1), and (1,q,1)(1,q,1)(1,q,1); the other three are scalar products. Tensor powers therefore contain many typed rectangular matrix products. The laser method selects a large family with disjoint coordinate blocks and applies Schönhage's asymptotic sum inequality to all surviving products simultaneously.

Formalization targets

Goal: the 1990 Coppersmith--Winograd bound

For every field KKK,

matMulExp⁡(K)<297125=2.376.\operatorname{matMulExp}(K)<\frac{297}{125}=2.376.matMulExp(K)<125297​=2.376.

The Lean goal has the same quantified proposition and the same matMulExp definition as the existing Schönhage-bound mission; only the theorem identifier and rational endpoint change.

Tensor and block foundations

The development records the characteristic-free order-three degeneration

Tq⊴Iq+2T_q\unlhd I_{q+2}Tq​⊴Iq+2​

and the exact matrix-product dimensions associated with every supported type sequence in Tq⊗NT_q^{\otimes N}Tq⊗N​. These statements identify the algebraic input before any asymptotic counting is used.

Coupled-weight extraction

For q=6q=6q=6, the tensor-square grading and the coupled-weight pruning must produce the direct sums and asymptotic inequality stated in Section 8 and in the coupled-constituent lemma on journal pp. 270--272. The final numerical milestone certifies the rational endpoint 297/125297/125297/125 from exact inequalities, rather than treating the decimal 2.3754772.3754772.375477 as a proof object.

Significance

The result was the strongest matrix-multiplication bound for roughly two decades and introduced the tensor family that underlies the classical laser-method line of work. A formal proof supplies a checked bridge from an explicit border-rank identity to an exponent bound whose combinatorial extraction is substantially more delicate than the earlier Schönhage examples.

The formalization also produces reusable infrastructure. The order-three CW degeneration is an explicit polynomial-family test case over arbitrary fields. The six block identifications and type-count formulas can be reused in analyses of tensor powers. A faithful extraction predicate, stated through actual restrictions to direct sums of MMObj tensors, separates sound laser arguments from formulas that count incompatible or coordinate-sharing blocks as independent.

The mathematical bound is known. The open work is its machine-checked reconstruction in Lean. The border-rank theorem, per-type matrix-product restriction layer, tensor-square support invariant, balanced block calculation, Salem--Spencer set theorem, and exact q=6q=6q=6 numerical endpoint are already proved. The unrestricted value/rank bridge, the coupled-constituent extraction, and the full Section 8 auxiliary inequality remain the substantive frontier.

Difficulty

The main difficulty is not expanding TqT_qTq​ or evaluating a decimal logarithm. A tensor power contains exponentially many typed terms, but most share variables. They cannot all be placed in a direct sum, and counting all joint type sequences overestimates the usable matrix products. The source hashes coordinate blocks into a large progression-free set and prunes collisions so that the surviving blocks are genuinely independent.

The 2.3762.3762.376 improvement adds a second layer. It begins with Tq⊗2T_q^{\otimes2}Tq⊗2​, regroups variables into five classes, couples weights that were independent in the simpler analysis, and estimates a nontrivial central block by a further extraction. A formal proof must track the direction of every restriction, the exact multiplicities of all block types, and the loss introduced by pruning. Replacing exponential surviving-block counts by a polynomial number of blocks, or using joint entropy without the marginal compatibility constraints, changes the mathematical claim and is outside the mission.

Formalization scope

The mission uses the existing TensorObj, MMObj, TensorObj.Restrict, Degenerates, tensorAsymptoticRank, matMulExp, and matMulExp_strassen declarations in the Mathlib environment pinned by the earlier matrix-multiplication mission. Tensor dimensions and type counts are natural numbers; exponent and optimization inequalities are real-valued. All top-level bounds quantify over an arbitrary field, matching the integral polynomial identities used by the construction.

Laser statements must exhibit, directly or through a faithful reusable predicate, restrictions from a tensor power to a finite direct sum of concrete matrix-multiplication tensors. The number and dimensions of the summands remain part of the witness. A constant-valued “laser functional,” a vacuous witness hypothesis, or a capacity definition that discards the exponential number of surviving blocks does not satisfy the mission.

Welcome contributions include restriction composition lemmas, tensor-power block equivalences, multinomial and entropy estimates with all marginal constraints, formal Salem--Spencer pruning, exact real-inequality certificates, and the coupled central-block value lemma. Every milestone should cite the corresponding equation, table, or lemma in the primary paper.

Selected references

  • Don Coppersmith and Shmuel Winograd, Matrix Multiplication via Arithmetic Progressions, Journal of Symbolic Computation 9, 1990, pp. 251--280. ScienceDirect.
  • Arnold Schönhage, Partial and Total Matrix Multiplication, SIAM Journal on Computing 10(3), 1981, pp. 434--455. DOI 10.1137/0210032.
  • Avi Wigderson and Jeroen Zuiddam, Asymptotic Spectra: Theory, Applications and Extensions, 2023, for the tensor restriction and asymptotic-rank framework used by the Lean development. Author manuscript.
72 thms5 active usersReviewed
Complexity Theory·Captain: avi

Integer multiplication: κ = 0.00012612978530233 (Jain round eight)Open Problem

Motivation

Multiplying two binary integers is a basic task in computational complexity. The relevant resource here is the number of individual steps made by a machine whose tapes can only be accessed through moving heads. Harvey and van der Hoeven established an O(nlog⁡n)O(n\log n)O(nlogn) bound in this model. The OpenAI manuscript proposes a positive saving in the exponent of the logarithm. Subsequent constructions explore larger savings. The formalization target is the complete computation and its time bound, including the analytic and implementation arguments needed to connect a finite construction to every input length.

Setting

A deterministic multitape Turing machine has a fixed finite alphabet, a fixed finite state set, and a fixed finite number of tapes. Each tape extends in one direction and has one head. A step reads the symbols under the heads, changes state, writes at those positions, and moves each head by at most one cell. The alphabet contains distinct blank, start-marker, zero, one, and separator symbols. Tape zero is the read-only input tape; tape one is the output tape; further tapes are available for work. The halting state freezes the configuration.

An input consists of two bit strings x,yx,yx,y of the same positive length nnn, separated by #\##. The strings are read most significant bit first, with leading zeros permitted. If val⁡(x)\operatorname{val}(x)val(x) denotes the represented natural number, the required output is exactly the 2n2n2n-bit binary representation of val⁡(x)val⁡(y)\operatorname{val}(x)\operatorname{val}(y)val(x)val(y), padded with leading zeros and followed by blanks. All tapes begin with a preserved start marker. Define L(n)=max⁡(⌈log⁡2n⌉,1)L(n)=\max(\lceil\log_2 n\rceil,1)L(n)=max(⌈log2​n⌉,1).

Target

The full target is

IntMul.KappaBound(0.00012612978530233).\mathrm{IntMul.KappaBound}(0.00012612978530233).IntMul.KappaBound(0.00012612978530233).

This asserts the existence of one machine that computes the exact product correctly for every positive input length, together with constants c>0c>0c>0 and n0n_0n0​ such that every pair of nnn-bit inputs, for n≥n0n\ge n_0n≥n0​ and n≥1n\ge1n≥1, halts within

c n L(n)1−0.00012612978530233c\,n\,L(n)^{1-0.00012612978530233}cnL(n)1−0.00012612978530233

steps. Constants, alphabet, states, and tape count are independent of the input length. The number is an exact real value. Larger κ\kappaκ gives a stronger asymptotic saving because L(n)≥1L(n)\ge1L(n)≥1.

Each entry supplies its own exact numerical value and identifies the construction that motivates it. Historical and new entries use the same machine definition. A proposed value may be recorded while its complete proof is still open; it becomes a proved record only when the goal has a verified proof. Numerical certificates are supporting results with their own scope, and do not establish a complete machine time bound.

This entry: Jain, round eight

The value

κ=126129785302331017=0.00012612978530233  (>2−13)\kappa=\frac{12612978530233}{10^{17}}=0.00012612978530233\;(>2^{-13})κ=101712612978530233​=0.00012612978530233(>2−13)

is the headline conditional witness of S. Jain's integer-mult-kappa, round eight, at revision 5868dc7. It is bound by the complex side (ac=3945999/31250000000a_c=3945999/31250000000ac​=3945999/31250000000).

  • Bit side: the frozen round-seven witness (deferred readouts, VVV leaves, lifted frames and late copies at h=23h=23h=23; Avi Eisenberg's producer from PR #62, compiled by eumemic's joint frame compiler), recompiled with opposite bank orders and a stopped interchange after IceKylin (PR #104).
  • Complex side: a signed reclaim with completed-core sharing and one fix-up child per group.

The repository states that the claim is conditional: it assumes the full upstream multiplication theorem (the OpenAI manuscript and Colkitt's framework) and that stage two is the complement time-reversal of stage one, and it rests on three further interfaces that are stated but not machine-checked there — IceKylin's normal form for alternating residuals, the cost of binary address adapters for the exterior fix-up children, and the row stock of the stopped bit interchange. Its Lean files check only finite rational statements (moment certificates, the stopped mix and the parameter assembly), not the analytic inequalities or the multiplication theorem. It has not had independent review.

Earlier claims appear as milestones: Harvey–van der Hoeven (κ=0\kappa=0κ=0), OpenAI Theorem 1 (2−1822^{-182}2−182), Colkitt's 2−302^{-30}2−30 checkpoint, Jain rounds six and seven, and community PR #62. Each is implied by this goal because L(n)1−κL(n)^{1-\kappa}L(n)1−κ decreases in κ\kappaκ.

Significance

A proof for any positive value would give a strict asymptotic improvement over nlog⁡nn\log nnlogn. Such an improvement would contradict the proposed Ω(nlog⁡n)\Omega(n\log n)Ω(nlogn) lower bound, while remaining consistent with the known upper bound. The tracked value concerns asymptotic growth; it does not measure the runtime of a practical implementation. No claim about a currently fastest implementation follows from a larger saving.

Formalizing the full target connects exact finite calculations to a uniform algorithm. The community repository explicitly treats the original framework and remaining all-size interfaces as assumptions of its current evidence. In this mission those steps must be proved before the goal can be marked Proved. They are not added as hypotheses that would weaken the full target. Earlier published theorem statements and finite certificates may be reused with their original attribution and proof status.

Difficulty

A finite network certificate can establish a concrete role count or a rational margin, but does not itself supply a machine for every input length. A complete proof must establish correctness of the recursive construction, admissibility of its parameters, precision and error guarantees, and the cost of moving data on tapes. Counts of algebraic operations alone do not account for this last resource. The difficult connection is between these separate contracts and the single uniform machine quantified by the goal.

Formalization scope

The goal imports the existing public Definitions.Def_IntMul_MultitapeModel bundle and uses IntMul.KappaBound directly. No alternative machine or cost model is introduced. Natural numbers describe input lengths, head positions, and step counts; real numbers describe the asymptotic constant and exponent. The real power has base at least one. Correctness is required even below the asymptotic threshold. Inputs of length zero and malformed inputs are outside the specified task.

The input and output conventions exclude a machine that merely decides a predicate about the product. The entire output tape must contain the required product and then blanks. Neither failure to halt nor correctness on a finite collection of test inputs meets the target. The work tapes and final head positions are unconstrained. Reusable contributions include machine composition with explicit costs, exact encodings, analytic estimates, and proofs that the finite certificates imply the contracts required by the recursive algorithm. Every such contribution should state its own assumptions and conclusion explicitly.

Selected references

  • D. Harvey and J. van der Hoeven, Integer multiplication in time O(n log n), Annals of Mathematics 193(2), 563–617 (2021). Paper.
  • OpenAI, Integer multiplication below n log n, preprint (23 September 2026), pinned manuscript revision adc7f1241b42e322a6451854ab7e4b4c146bf78a. Source.
  • D. Colkitt and community contributors, A sharper exponent for integer multiplication, research repository (2026), revision 56b66d58297deca1d7dd130247d720e960f77a37. Checkpoint and attribution.
  • S. Jain, Integer multiplication: a conditional witness above 2−132^{-13}2−13 (integer-mult-kappa), research draft (2026), revision 5868dc701a9bec375c9dc5391d992ab3b8c35f3f. Round eight.
16 thms4 active usersReviewed
Complexity Theory·Captain: avi

Integer multiplication: κ = 0.00006397891967576 (Jain round seven)Open Problem

Motivation

Multiplying two binary integers is a basic task in computational complexity. The relevant resource here is the number of individual steps made by a machine whose tapes can only be accessed through moving heads. Harvey and van der Hoeven established an O(nlog⁡n)O(n\log n)O(nlogn) bound in this model. The OpenAI manuscript proposes a positive saving in the exponent of the logarithm. Subsequent constructions explore larger savings. The formalization target is the complete computation and its time bound, including the analytic and implementation arguments needed to connect a finite construction to every input length.

Setting

A deterministic multitape Turing machine has a fixed finite alphabet, a fixed finite state set, and a fixed finite number of tapes. Each tape extends in one direction and has one head. A step reads the symbols under the heads, changes state, writes at those positions, and moves each head by at most one cell. The alphabet contains distinct blank, start-marker, zero, one, and separator symbols. Tape zero is the read-only input tape; tape one is the output tape; further tapes are available for work. The halting state freezes the configuration.

An input consists of two bit strings x,yx,yx,y of the same positive length nnn, separated by #\##. The strings are read most significant bit first, with leading zeros permitted. If val⁡(x)\operatorname{val}(x)val(x) denotes the represented natural number, the required output is exactly the 2n2n2n-bit binary representation of val⁡(x)val⁡(y)\operatorname{val}(x)\operatorname{val}(y)val(x)val(y), padded with leading zeros and followed by blanks. All tapes begin with a preserved start marker. Define L(n)=max⁡(⌈log⁡2n⌉,1)L(n)=\max(\lceil\log_2 n\rceil,1)L(n)=max(⌈log2​n⌉,1).

Target

The full target is

IntMul.KappaBound(0.00006397891967576).\mathrm{IntMul.KappaBound}(0.00006397891967576).IntMul.KappaBound(0.00006397891967576).

This asserts the existence of one machine that computes the exact product correctly for every positive input length, together with constants c>0c>0c>0 and n0n_0n0​ such that every pair of nnn-bit inputs, for n≥n0n\ge n_0n≥n0​ and n≥1n\ge1n≥1, halts within

c n L(n)1−0.00006397891967576c\,n\,L(n)^{1-0.00006397891967576}cnL(n)1−0.00006397891967576

steps. Constants, alphabet, states, and tape count are independent of the input length. The number is an exact real value. Larger κ\kappaκ gives a stronger asymptotic saving because L(n)≥1L(n)\ge1L(n)≥1.

Each entry supplies its own exact numerical value and identifies the construction that motivates it. Historical and new entries use the same machine definition. A proposed value may be recorded while its complete proof is still open; it becomes a proved record only when the goal has a verified proof. Numerical certificates are supporting results with their own scope, and do not establish a complete machine time bound.

This entry: Jain, round seven

The value

κ=79973649594712500000000000000=0.00006397891967576  (>2−14)\kappa=\frac{799736495947}{12500000000000000}=0.00006397891967576\;(>2^{-14})κ=12500000000000000799736495947​=0.00006397891967576(>2−14)

is the conditional witness of S. Jain's integer-mult-kappa, round seven, at revision c05adfe, with the sharper moment certificate merged there (research/round7-moment-refinement/). The README's headline for the same construction is the slightly smaller 1599247689723/(2.5⋅1016)≈6.39699×10−51599247689723/(2.5\cdot10^{16})\approx6.39699\times10^{-5}1599247689723/(2.5⋅1016)≈6.39699×10−5; the refinement keeps the frozen round-seven bit network, histogram, complex saving and assembly unchanged, and replaces the reciprocal bound eu≤1/(1−u)e^u\le1/(1-u)eu≤1/(1−u) in the moment certificate by a degree-eight Taylor enclosure, which certifies the bit saving ab=63983013240044/1018a_b=63983013240044/10^{18}ab​=63983013240044/1018.

The construction builds on the OpenAI manuscript and on Colkitt's integer-mult-bounds (it uses the PR #62 producer with interval strips and core-aware pair assembly). Its new ingredients are lifted frames with late copies, deferred readouts, and per-use VVV leaves on the bit side; the complex side is round six's. The repository states that the claim is conditional on the OpenAI manuscript and Colkitt's framework, that stage two (the complement time-reversal of stage one) is assumed rather than machine-checked, and that its Lean files check only finite rational statements (moment certificates and parameter assembly), not the analytic tail inequalities or the multiplication theorem. It has not had independent review.

Earlier claims appear as milestones: Harvey–van der Hoeven (κ=0\kappa=0κ=0), OpenAI Theorem 1 (2−1822^{-182}2−182), Colkitt's 2−302^{-30}2−30 checkpoint, Jain round six (3666565558019/10173666565558019/10^{17}3666565558019/1017) and community PR #62 (0.0000510169201700780.0000510169201700780.000051016920170078). Each is implied by this goal because L(n)1−κL(n)^{1-\kappa}L(n)1−κ decreases in κ\kappaκ.

Significance

A proof for any positive value would give a strict asymptotic improvement over nlog⁡nn\log nnlogn. Such an improvement would contradict the proposed Ω(nlog⁡n)\Omega(n\log n)Ω(nlogn) lower bound, while remaining consistent with the known upper bound. The tracked value concerns asymptotic growth; it does not measure the runtime of a practical implementation. No claim about a currently fastest implementation follows from a larger saving.

Formalizing the full target connects exact finite calculations to a uniform algorithm. The community repository explicitly treats the original framework and remaining all-size interfaces as assumptions of its current evidence. In this mission those steps must be proved before the goal can be marked Proved. They are not added as hypotheses that would weaken the full target. Earlier published theorem statements and finite certificates may be reused with their original attribution and proof status.

Difficulty

A finite network certificate can establish a concrete role count or a rational margin, but does not itself supply a machine for every input length. A complete proof must establish correctness of the recursive construction, admissibility of its parameters, precision and error guarantees, and the cost of moving data on tapes. Counts of algebraic operations alone do not account for this last resource. The difficult connection is between these separate contracts and the single uniform machine quantified by the goal.

Formalization scope

The goal imports the existing public Definitions.Def_IntMul_MultitapeModel bundle and uses IntMul.KappaBound directly. No alternative machine or cost model is introduced. Natural numbers describe input lengths, head positions, and step counts; real numbers describe the asymptotic constant and exponent. The real power has base at least one. Correctness is required even below the asymptotic threshold. Inputs of length zero and malformed inputs are outside the specified task.

The input and output conventions exclude a machine that merely decides a predicate about the product. The entire output tape must contain the required product and then blanks. Neither failure to halt nor correctness on a finite collection of test inputs meets the target. The work tapes and final head positions are unconstrained. Reusable contributions include machine composition with explicit costs, exact encodings, analytic estimates, and proofs that the finite certificates imply the contracts required by the recursive algorithm. Every such contribution should state its own assumptions and conclusion explicitly.

Selected references

  • D. Harvey and J. van der Hoeven, Integer multiplication in time O(n log n), Annals of Mathematics 193(2), 563–617 (2021). Paper.
  • OpenAI, Integer multiplication below n log n, preprint (23 September 2026), pinned manuscript revision adc7f1241b42e322a6451854ab7e4b4c146bf78a. Source.
  • D. Colkitt and community contributors, A sharper exponent for integer multiplication, research repository (2026), revision 56b66d58297deca1d7dd130247d720e960f77a37. Checkpoint and attribution.
  • S. Jain, Integer multiplication: a conditional witness above 2−142^{-14}2−14 (integer-mult-kappa), research draft (2026), revision c05adfe1dfe85326dd6dcf12bfd53f435bd49bf3. Round seven and the sharper moment certificate.
16 thms4 active usersReviewed
🏆Completed
Machine LearningOperations Research·Captain: mikedeng1

The Weighted Majority Algorithm 9: With η Anomalies Every Randomized Learner Can Be Forced to ½ opt(F, 0) + η Expected MistakesResearch Paper

Motivation

Prediction with expert advice asks a learner to predict a binary label on each of a sequence of trials, after seeing the predictions of a pool of candidate functions, and to make few more mistakes than the best member of the pool. Littlestone and Warmuth's The Weighted Majority Algorithm (Inform. and Comput. 108, 1994) gave the deterministic algorithm WM and its randomized variant WMR, and proved worst-case bounds that hold for every sequence of trials, with no probabilistic assumption on the data. Their §8 asks how good such bounds can possibly be when the data contain anomalies: disagreements with a candidate function, counted relative to the best-fitting member of the pool. Such disagreements can arise from noise in the instances or labels, or from instances that do not determine the label.

For the randomized algorithm WMR the expected number of mistakes grows, by Corollary 6.1 of the paper, at rate ln⁡(1/β)/(1−β)\ln(1/\beta)/(1-\beta)ln(1/β)/(1−β) per anomaly, which tends to 111 as β→1\beta \to 1β→1. Theorem 8.2 shows that no randomized learner, however it is built, can do better than rate 111: an adversary who fixes the sequence in advance can force half the noise-free optimum plus one expected mistake per anomaly. Together with the deterministic counterpart, Theorem 8.1 (rate 222), it locates the cost of label noise in online learning exactly. Mistake-bound lower bounds of this kind are a standard reference point for the analysis of online classification and of the expert-advice algorithms descended from WM (Cesa-Bianchi and Lugosi, Prediction, Learning, and Games, 2006).

Timeline: Littlestone (1988) characterized the optimal noise-free mistake bound opt(F,0)\mathrm{opt}(F,0)opt(F,0) by complete mistake trees (the Littlestone dimension); Littlestone and Warmuth (FOCS 1989; journal version 1994) introduced WM and WMR and proved Theorems 8.1 and 8.2.

Setting

Let XXX be an arbitrary domain and FFF a pool of functions X→{0,1}X \to \{0,1\}X→{0,1}, finite or infinite. A trial is a pair (x,y)(x,y)(x,y) with x∈Xx\in Xx∈X and y∈{0,1}y\in\{0,1\}y∈{0,1}; a sequence of trials is a finite list S=((x1,y1),…,(xT,yT))\mathcal S = ((x_1,y_1),\dots,(x_T,y_T))S=((x1​,y1​),…,(xT​,yT​)), in which an instance may repeat with different labels. The number of anomalies of S\mathcal SS with respect to fff is #{t:f(xt)≠yt}\#\{t : f(x_t)\ne y_t\}#{t:f(xt​)=yt​}, and SηS_\etaSη​ is the collection of sequences for which some f∈Ff\in Ff∈F has at most η\etaη anomalies.

A deterministic algorithm maps the past trials and the current instance to a predicted label; MA(S)M_A(\mathcal S)MA​(S) counts its wrong predictions, and

opt(F,η)=inf⁡Asup⁡S∈SηMA(S).\mathrm{opt}(F,\eta) = \inf_A \sup_{\mathcal S\in S_\eta} M_A(\mathcal S).opt(F,η)=Ainf​S∈Sη​sup​MA​(S).

A randomized algorithm uses random choices independent of the sequence. On a fixed sequence its behaviour at trial ttt is summarized by pt∈[0,1]p_t\in[0,1]pt​∈[0,1], the probability that it predicts 111 given the preceding instances and labels and the current instance; its expected number of mistakes is ∑t∣pt−yt∣\sum_t |p_t - y_t|∑t​∣pt​−yt​∣, and

optRAND(F,η)=inf⁡Asup⁡S∈Sη∑t=1T∣pt−yt∣,\mathrm{opt}_{\mathrm{RAND}}(F,\eta) = \inf_A \sup_{\mathcal S\in S_\eta} \sum_{t=1}^T |p_t-y_t|,optRAND​(F,η)=Ainf​S∈Sη​sup​t=1∑T​∣pt​−yt​∣,

the infimum ranging over all maps AAA from (past trials, current instance) to pt∈[0,1]p_t \in [0,1]pt​∈[0,1].

Formalization targets

Goal: Theorem 8.2 (p. 252)

For every pool FFF with ∣F∣>1|F|>1∣F∣>1 and every integer η≥0\eta\ge 0η≥0,

optRAND(F,η) ≥ 12 opt(F,0)+η.\mathrm{opt}_{\mathrm{RAND}}(F,\eta)\ \ge\ \tfrac12\,\mathrm{opt}(F,0)+\eta .optRAND​(F,η) ≥ 21​opt(F,0)+η.

Both sides live in [0,∞][0,\infty][0,∞]; when opt(F,0)=∞\mathrm{opt}(F,0)=\inftyopt(F,0)=∞ the statement says optRAND(F,η)=∞\mathrm{opt}_{\mathrm{RAND}}(F,\eta)=\inftyoptRAND​(F,η)=∞.

Milestone: complete mistake trees (§8, p. 251)

"For any pool of functions FFF there exists a complete opt(F,0)\mathrm{opt}(F,0)opt(F,0)-mistake tree" ([Lit88]). On the platform this is the proved Lemma 21.7 of Shalev-Shwartz and Ben-David, MSOA(F)≤Ldim(F)M_{\mathrm{SOA}}(F)\le \mathrm{Ldim}(F)MSOA​(F)≤Ldim(F), together with Lemma 21.6, Ldim(F)≤MA(F)\mathrm{Ldim}(F) \le M_A(F)Ldim(F)≤MA​(F) for every AAA; a tree shattered by FFF is a complete mistake tree.

Significance

The theorem shows that the expected-mistake bound of WMR is optimal in its dependence on the number of anomalies: its coefficient of η\etaη approaches 111, and no randomized learner achieves a smaller one on every pool with two or more functions. It also quantifies the advantage of randomization: the deterministic rate is 222 (Theorem 8.1), the randomized rate is 111, and the noise-free term halves. The lower bound is obliviously adversarial, since the sequence is chosen before any coin is tossed, so it applies equally to the stronger model in which the adversary adapts to the learner's predictions.

The result has been proved since 1994. To our knowledge it has no machine-checked proof; the platform holds the Littlestone-dimension characterization of the noise-free optimum (Shalev-Shwartz and Ben-David, Chapter 21) but no statement about anomalies or randomized lower bounds. This mission formalizes the definitions of anomalies, opt(F,η)\mathrm{opt}(F,\eta)opt(F,η) and optRAND(F,η)\mathrm{opt}_{\mathrm{RAND}}(F,\eta)optRAND​(F,η), and the theorem.

Difficulty

An adversary against a randomized learner cannot react to its predictions: the whole sequence, instances and labels, must be fixed in advance, so the adaptive adversary that proves the deterministic bound is not directly available. The lower bound must instead come from the learner's prediction probabilities, which depend only on the sequence prefix; the content of the theorem is that this information is enough to construct, for every learner, one fixed sequence that stays inside SηS_\etaSη​ while forcing the stated expected count. Taking the infimum over learners inside an extended-real supremum over sequences of every length adds bookkeeping when opt(F,0)\mathrm{opt}(F,0)opt(F,0) is infinite.

Formalization scope

Sequences are Fin T → X × Bool for every T : ℕ. Deterministic algorithms, the mistake count, shattered trees and the Littlestone dimension are reused from the published definition UnderstandingML_Online (OnlineAlg, mistakes, ShattersTree, ldim, soa); randomized algorithms are its OnlineAlgR X := List (X × Bool) → X → ℝ restricted by IsRandAlg to values in [0,1][0,1][0,1], and the expected mistake count is its cumLoss. η\etaη is a natural number. opt(F,η)\mathrm{opt}(F,\eta)opt(F,η) takes values in ℕ∞ and optRAND(F,η)\mathrm{opt}_{\mathrm{RAND}}(F,\eta)optRAND​(F,η) in ℝ≥0∞; the paper's "minimum" is read as an infimum. ∣F∣>1|F|>1∣F∣>1 is F.Nontrivial, which admits infinite pools.

The reduction of randomized algorithms to probability maps is the main modelling choice: on a fixed sequence the expected number of mistakes of any algorithm whose randomization is independent of the sequence equals ∑t∣pt−yt∣\sum_t|p_t-y_t|∑t​∣pt​−yt​∣, and every map with values in [0,1][0,1][0,1] is realized by such an algorithm, so the two infima coincide. Dropping the constraint pt∈[0,1]p_t\in[0,1]pt​∈[0,1], or quantifying only over deterministic algorithms, gives a different quantity and is ruled out.

Needed: the definitions above, and the existence of a complete opt(F,0)\mathrm{opt}(F,0)opt(F,0)-mistake tree, available through the proved platform theorems on ldim and soa. Lemmas relating mistakes along a path of a shattered tree to ldim, and a reusable statement that an oblivious adversary can force expected loss 12\tfrac1221​ per trial, would serve other lower bounds in the Littlestone–Warmuth series and in online learning generally. Contributions of such lemmas are welcome.

Selected references

  • N. Littlestone, M. K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2):212–261, 1994. https://doi.org/10.1006/inco.1994.1009
  • N. Littlestone, Learning Quickly When Irrelevant Attributes Abound: A New Linear-Threshold Algorithm, Machine Learning 2:285–318, 1988. https://doi.org/10.1007/BF00116827
  • S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 21. https://doi.org/10.1017/CBO9781107298019
  • N. Cesa-Bianchi, G. Lugosi, Prediction, Learning, and Games, Cambridge University Press, 2006. https://doi.org/10.1017/CBO9780511546921
7 thms4 active usersReviewed
🏆Completed
CombinatoricsMachine LearningOperations Research·Captain: mikedeng1

The Weighted Majority Algorithm 7: On a Class Shattered by Its Domain, Every Deterministic Learner's Mistake Bounds Satisfy Σ 2^(−M_i) < 2Research Paper

Motivation

An online learner predicts a binary label for each instance before seeing that instance's correct label. Its guarantee is usually measured by its worst possible number of mistakes against a target function. Such guarantees can differ greatly across target functions: an algorithm might learn one target quickly while taking much longer on another. Littlestone and Warmuth's analysis of the Weighted Majority Algorithm asks which collections of individual mistake guarantees can coexist for one deterministic learner. In Section 7 of their 1994 paper, they isolate a condition under which these guarantees obey a strict numerical constraint. The condition concerns how the target functions vary together over the instance domain, so the result applies without a probability distribution on instances or labels.

The theorem also supports the paper's comparison between arbitrary deterministic learners and Weighted Majority. Given individual mistake bounds, its inequality is the numerical condition that makes a suitable initial weighting of the target functions possible in Theorem 7.1 and Corollary 7.1 of the same paper. The present mission centers on the inequality itself, which is independent of a particular implementation of Weighted Majority.

Setting

Let XXX be an instance domain, possibly finite or infinite, and let φ1,…,φn:X→{0,1}\varphi_1,\ldots,\varphi_n:X\to\{0,1\}φ1​,…,φn​:X→{0,1} be a finite family of target functions. A trial presents an instance x∈Xx\in Xx∈X to a learner, asks it to predict either 000 or 111, and then reveals a label. A finite sequence is consistent with φi\varphi_iφi​ when each revealed label equals φi(x)\varphi_i(x)φi​(x) at that trial's instance. The learner may use the current instance and the complete history of earlier labeled trials when predicting; it does not see the current label first. No sampling or independence assumption is made.

The family is shattered by its domain when every pattern b=(b1,…,bn)∈{0,1}nb=(b_1,\ldots,b_n)\in\{0,1\}^nb=(b1​,…,bn​)∈{0,1}n appears as the vector of function values at some instance:

∀b∈{0,1}n∃x∈X∀i∈{1,…,n}, φi(x)=bi.\forall b\in\{0,1\}^n\quad\exists x\in X\quad\forall i\in\{1,\ldots,n\},\ \varphi_i(x)=b_i.∀b∈{0,1}n∃x∈X∀i∈{1,…,n}, φi​(x)=bi​.

This is a dual notion of shattering. The familiar VC condition asks whether a class of functions realizes every labeling of chosen instances. Here the domain realizes every joint value pattern of a chosen family of functions. The two conditions should not be interchanged.

Let AAA be any deterministic online learner. Write MiM_iMi​ for the supremum of the number of mistakes made by AAA over every finite sequence consistent with φi\varphi_iφi​. Since mistake counts are natural numbers, a finite supremum is an attained natural-number maximum. The theorem assumes each MiM_iMi​ is finite. In the Lean representation, this is a separate finiteness condition on each target's extended-natural mistake bound.

Formalization targets

Strict constraint on individual mistake bounds

The goal is Theorem 7.2 of Littlestone and Warmuth, pp. 244–245. For every domain, every finite family shattered by that domain, and every deterministic online learner with finite individual mistake bounds,

∑i=1n2−Mi<2.\sum_{i=1}^{n}2^{-M_i}<2.i=1∑n​2−Mi​<2.

The inequality is strict. Replacing it by a weak inequality would lose the stated result and the threshold used in the paper's later comparison with Weighted Majority. The theorem quantifies over the learner itself; it does not select a specially behaved learning rule.

Arithmetic milestones

Lemma 7.1 concerns positive real numbers with integral adjacent quotients. When a positive target sss is no larger than their sum and s/r1s/r_1s/r1​ is integral, some prefix sums exactly to sss. Lemma 7.2 concerns positive powers of two with integer exponents. If a power 2ℓ2^\ell2ℓ lies at least as high as each individual weight and at most as high as the total weight, some subset sums exactly to 2ℓ2^\ell2ℓ. Both statements and their exact source wording are recorded as the mission's milestones; they are the numbered arithmetic results used in Theorem 7.2.

Significance

The theorem limits the set of claims that one deterministic learner can make simultaneously. Small MiM_iMi​ contribute large terms to the sum, so a learner cannot have arbitrarily small worst-case mistake counts for all members of a domain-shattered family. This is a worst-case statement about every finite consistent trial sequence, rather than an average-case statement about a chosen data distribution. The paper's Corollary 7.1 uses the resulting inequality to select initial weights for a Weighted Majority comparison, but that later algorithmic statement is outside this mission's goal.

Formalizing this result adds a reusable dual-shattering predicate and an exact constraint on per-target online mistake bounds to the existing online-learning definitions. The paper proves the theorem, while the mission's Lean goal and its two milestone lemmas are open proof obligations. The separate published online-learning module provides the definitions of deterministic predictors, trial histories, mistake counts, and extended-natural worst-case bounds; it does not already prove this theorem.

Difficulty

For a single target, an arbitrary learner may have a small finite mistake bound. The difficulty is relating those bounds across different targets without placing extra restrictions on the learner. Domain shattering gives the adversary control over the simultaneous predictions of all functions at an instance, but the labels of one trial sequence must still remain consistent with a target. A simple count of targets or of possible patterns does not directly produce the strict weighted sum, particularly when the different MiM_iMi​ are unequal. The two arithmetic lemmas address exact sums of structured weights; the main theorem must also connect such weights to the learner's worst-case mistakes.

Formalization scope

The function family is indexed by Fin n, with zero-based Lean indices corresponding to the paper's 1,…,n1,\ldots,n1,…,n. The domain XXX is an arbitrary type and may be infinite. Labels and predictions are Bool. ShatteredByDomain states that every Boolean pattern is realized by some domain point. The published UnderstandingML_Online definition gives OnlineAlg X Bool, finite trial histories, mistakes, and mistakeBound. The latter is valued in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞} and ranges over all finite sequences labeled by the target. Each MiM_iMi​ is converted to a natural number only after assuming its bound is finite. The powers 2−Mi2^{-M_i}2−Mi​ and their finite sum are real numbers. The n=0n=0n=0 conclusion is the empty sum 0<20<20<2; the shattering predicate itself then requires an instance to realize the unique empty pattern.

The two milestones use real-valued weights. Lemma 7.1 assumes n≥1n\ge1n≥1 because the source refers to r1r_1r1​. Lemma 7.2 represents the source's integral base-two logarithms by equations ri=2kir_i=2^{k_i}ri​=2ki​ with integer kik_iki​, which also enforce positivity. Its maximum hypothesis is written pointwise, so an empty family is harmless: the other inequality 2ℓ≤02^\ell\le02ℓ≤0 cannot hold. A coordinate family on X={0,1}nX=\{0,1\}^nX={0,1}n is verified locally to satisfy the shattering predicate, ruling out a vacuous encoding of that central assumption. Contributions proving either arithmetic milestone or the main theorem within these exact definitions are in scope.

Selected references

  • Nick Littlestone and Manfred K. Warmuth, The Weighted Majority Algorithm, Information and Computation 108(2), 212–261, 1994. DOI: 10.1006/inco.1994.1009.
  • Shai Shalev-Shwartz and Shai Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Chapter 21. DOI: 10.1017/CBO9781107298019.
5 thms4 active usersReviewed
🏆Completed
CombinatoricsOperations ResearchProbability·Captain: mikedeng1

An Optimal Algorithm for On-line Bipartite Matching 1: RANKING Finds a Matching of Expected Size at Least n(1 − 1/e) − o(n) on Every Graph with a Perfect MatchingResearch Paper

Motivation

Online matching describes allocation when requests must be answered as they arrive. A matching decision uses only the edges revealed so far and cannot be revised when later requests appear. In the bipartite setting studied by Karp, U. Vazirani, and V. Vazirani, the arriving vertices are girls and the possible partners are boys. Even when the full graph has a perfect matching, a fixed greedy priority can leave many girls unmatched. The paper introduced RANKING, which randomly chooses the boys' priority order once and then uses that order for every arrival.

The question is quantitative: how many pairs does RANKING guarantee in expectation against a graph and an arrival order chosen before its random ranking? The paper's target is a fraction approaching 1−1/e1-1/e1−1/e of the nnn pairs in a perfect matching. Its printed Theorem 1 concerns an auxiliary algorithm called EARLY; the RANKING statement follows the intended chain through Lemmas 3 and 5. The original EARLY analysis has a gap for general upper-triangular matrices, so this mission states the RANKING target and the earlier, unaffected lemmas separately. This distinction matters because an assertion about EARLY would be a different formalization target.

Setting

A bipartite graph has a boy side UUU and a girl side VVV, each with nnn vertices. An edge (u,v)(u,v)(u,v) means that boy uuu may be paired with girl vvv. A matching is a collection of edges in which no boy or girl occurs twice. The standing hypothesis for the performance guarantee is that the graph has a perfect matching: some bijection from boys to girls selects an edge for every boy. The graph is otherwise arbitrary.

Girls arrive in a predetermined order. When girl vvv arrives, only her incident edges are revealed. RANKING first chooses a uniformly random permutation π\piπ of the boys. For each arriving girl, it selects the highest-ranked adjacent boy who is still unmatched, if one exists. Write MR(G,π)M_{\mathrm R}(G,\pi)MR​(G,π) for the final matching. The expected size is the finite average over all n!n!n! rankings. The guarantee must hold for every graph and every predetermined arrival order; relabeling girls lets the formal statement fix their order to n,n−1,…,1n,n-1,\ldots,1n,n−1,…,1.

The paper also uses a dual rows-arrive view. Boys arrive in an order and choose among eligible girls, whose priority order is fixed. This is the same greedy rule with the sides exchanged. For its triangular reduction, the columns are numbered 1,…,n1,\ldots,n1,…,n and column nnn has highest priority. An upper-triangular matrix with unit diagonal is a graph with every edge (i,i)(i,i)(i,i) and with an edge (i,j)(i,j)(i,j) only when i≤ji\le ji≤j. On such a graph, the auxiliary algorithm EARLY declines to match row iii if column iii is already covered by EARLY's own matching. For any matching MMM, D(M)D(M)D(M) denotes the indices for which both row iii and column iii are covered.

Formalization targets

RANKING guarantee

For every ε>0\varepsilon>0ε>0, one threshold NNN must work for every size n≥Nn\ge Nn≥N and every graph GGG with a perfect matching:

Eπ∼Unif(Sn)∣MR(G,π)∣≥(1−e−1−ε)n.\mathbb E_{\pi\sim\mathrm{Unif}(S_n)}|M_{\mathrm R}(G,\pi)| \ge (1-e^{-1}-\varepsilon)n.Eπ∼Unif(Sn​)​∣MR​(G,π)∣≥(1−e−1−ε)n.

This is the uniform lower-bound reading of the paper's n(1−1/e)−o(n)n(1-1/e)-o(n)n(1−1/e)−o(n) target. It does not prescribe a finite-nnn additive constant. The mission goal is the RANKING assertion drawn from the paper's Theorem 1 and Lemmas 3 and 5. The six milestones formalize the paper's Lemmas 1–5 and the corollary to Lemma 4, in their source order. They cover the side-exchange identity, arbitrary refusal algorithms, the triangular reduction, the matching-count identity, its expected form, and the pointwise comparison with EARLY.

Significance

The guarantee gives a concrete worst-case floor for a simple randomized allocation rule: as the graph size grows, RANKING matches at least an asymptotic 1−1/e1-1/e1−1/e fraction of the pairs available in a perfect matching, in expectation. The order is selected before the random permutation, matching the paper's performance measure. The lower bound remains meaningful on sparse graphs and does not rely on a density assumption. The original paper also studies the limit on what any randomized online algorithm can guarantee; that upper-bound result is treated in a separate mission.

A formal development here would provide reusable finite definitions for online greedy matching, arbitrary state-dependent refusal, fixed-order duality, and uniform expectation over permutations. The goal is a theorem statement awaiting a machine-checked proof; compiling the draft declarations verifies their Lean syntax and types, not their truth. The local mission separates the valid early structural statements from later statements whose published EARLY argument does not justify them on general upper-triangular graphs.

Difficulty

The random permutation does not make the fate of different vertices independent. Matching one girl removes a boy who might be essential to a later girl, so a per-arrival probability estimate cannot simply be added across all arrivals. The dual and triangular views capture useful structure, but turning that structure into a uniform bound for every graph is the main obstacle. In particular, reasoning about EARLY as if its matched-column set had the same monotonicity as unrestricted RANKING fails on some upper-triangular matrices. A proof of the goal must establish the RANKING guarantee without treating those later EARLY statements as available facts.

Formalization scope

Boys and girls are both Fin n; an adjacency matrix is a relation between them. A published predicate represents a perfect matching as an edge-preserving bijection. The generic greedy run processes arrival times in increasing order and interprets a smaller priority index as higher rank. It makes a new matching decision only from the current matching and the arriving vertex. RANKING's returned edges are consistently ordered as (boy, girl), even though girls arrive. The dual run has rows arriving in an arbitrary permutation and takes column n−1n-1n−1 as highest priority in Lean's zero-based numbering. EARLY's refusal test consults the matching constructed by EARLY itself.

All expectations are finite averages over Equiv.Perm (Fin n), scaled by 1/n!1/n!1/n!. The formal o(n)o(n)o(n) claim is ∀ε>0, ∃N, ∀n≥N, ∀G\forall\varepsilon>0,\ \exists N,\ \forall n\ge N,\ \forall G∀ε>0, ∃N, ∀n≥N, ∀G with a perfect matching, the displayed lower bound. Placing NNN before GGG preserves the worst-case meaning. The perfect-matching condition is essential: the empty graph cannot satisfy a positive linear guarantee. The triangular matrix in Lemma 3 retains every diagonal edge, so a zero matrix cannot witness the reduction. These conventions rule out vacuous versions of the goal and reduction.

The definitions of partial matching, covered vertices, greedy run, RANKING, EARLY, and uniform average are part of the mission. A complete proof may build further finite counting and permutation machinery; such lemmas can be shared beyond this paper. Contributions toward the six milestone statements, the goal, and faithful supporting results are in scope. Lemmas 6–12 and the printed EARLY Theorem 1 are outside this mission because their analysis depends on claims that fail for some matrices allowed by their surrounding assumptions.

Selected references

  • Richard M. Karp, Umesh V. Vazirani, and Vijay V. Vazirani, An Optimal Algorithm for On-line Bipartite Matching, Proceedings of the 22nd Annual ACM Symposium on Theory of Computing, 1990, pp. 352–358. DOI: 10.1145/100216.100262.
10 thms4 active usersReviewed
PreviousPage 1 of 19Next

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