Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Number Theory

103 missions · 51 completed

The study of the integers and the structures built from them — prime numbers, and the rational, algebraic, and ppp-adic numbers that extend them. It reaches from analytic number theory, which uses the tools of analysis to understand the distribution of primes, to algebraic number theory, Diophantine equations, and the arithmetic of elliptic curves, modular forms, and LLL-functions.

Missions

Open52Completed51All103
Captain: xuanji

The irrationality measure of π is at most 7.606309 (Salikhov 2008)Research Paper

Motivation

The irrationality measure μ(π)\mu(\pi)μ(π) is the supremum of the μ\muμ for which ∣π−p/q∣<q−μ|\pi - p/q| < q^{-\mu}∣π−p/q∣<q−μ has infinitely many rational solutions p/qp/qp/q. Every irrational number has μ≥2\mu \ge 2μ≥2 (Dirichlet), almost every real number has μ=2\mu = 2μ=2, and it is conjectured that μ(π)=2\mu(\pi) = 2μ(π)=2. Known upper bounds:

  • Mahler (1953): 424242, the first proof that π\piπ is not a Liouville number.
  • Mignotte (1974): 20.620.620.6.
  • Chudnovsky (1982): 19.8899944…19.8899944\ldots19.8899944…
  • Rhin–Viola (1993): 14.79707414.79707414.797074.
  • Hata (1993): 8.016045…8.016045\ldots8.016045…
  • Salikhov (2008): 7.606308…7.606308\ldots7.606308…
  • Zeilberger–Zudilin (2020): 7.103205334137…7.103205334137\ldots7.103205334137…, the current record.

The campaign's first proved value is Mahler's 424242. This entry records Salikhov's bound.

Formalization target

The campaign template with the value 7.6063097.6063097.606309 filled in: PiIrrationality.UpperBound (7.606309 : ℝ), i.e. μ(π)≤7.606309\mu(\pi) \le 7.606309μ(π)≤7.606309.

Value. The bound is quoted as 7.606308…7.606308\ldots7.606308… (e.g. by Zeilberger–Zudilin), a truncation. This entry rounds the last digit up to 7.6063097.6063097.606309.

How the bound arises

Salikhov uses integrals of rational functions that are symmetric under a group of transformations (in the spirit of Rhin–Viola), which gives larger arithmetic savings in the common denominators than Hata's construction.

Significance

Each step down the list replaces Mahler's approximations with a sharper family. Formalizing 7.6063097.6063097.606309 would build reusable explicit machinery: integral constructions of rational approximations to π\piπ, bounds on their common denominators via prime-number estimates, and the standard lemma turning a sequence of good approximations into an irrationality-measure bound.

Selected references

  • V. Kh. Salikhov, On the irrationality measure of π\piπ, Russian Math. Surveys 63 (2008), no. 3, 570–572.
  • K. Mahler, On the approximation of π\piπ, Indag. Math. 15 (1953), 30–42.
  • F. Beukers, A rational approach to π\piπ, Nieuw Arch. Wiskd. (5) 1 (2000), 372–379.
  • Source table: https://teorth.github.io/optimizationproblems/constants/7a.html
11 thms3 active usersReviewed
Captain: xuanji

The irrationality measure of π is at most 8.016046 (Hata 1993)Research Paper

Motivation

The irrationality measure μ(π)\mu(\pi)μ(π) is the supremum of the μ\muμ for which ∣π−p/q∣<q−μ|\pi - p/q| < q^{-\mu}∣π−p/q∣<q−μ has infinitely many rational solutions p/qp/qp/q. Every irrational number has μ≥2\mu \ge 2μ≥2 (Dirichlet), almost every real number has μ=2\mu = 2μ=2, and it is conjectured that μ(π)=2\mu(\pi) = 2μ(π)=2. Known upper bounds:

  • Mahler (1953): 424242, the first proof that π\piπ is not a Liouville number.
  • Mignotte (1974): 20.620.620.6.
  • Chudnovsky (1982): 19.8899944…19.8899944\ldots19.8899944…
  • Rhin–Viola (1993): 14.79707414.79707414.797074.
  • Hata (1993): 8.016045…8.016045\ldots8.016045…
  • Salikhov (2008): 7.606308…7.606308\ldots7.606308…
  • Zeilberger–Zudilin (2020): 7.103205334137…7.103205334137\ldots7.103205334137…, the current record.

The campaign's first proved value is Mahler's 424242. This entry records Hata's bound.

Formalization target

The campaign template with the value 8.0160468.0160468.016046 filled in: PiIrrationality.UpperBound (8.016046 : ℝ), i.e. μ(π)≤8.016046\mu(\pi) \le 8.016046μ(π)≤8.016046.

Value. The paper computes the exponent as 7.016045…+17.016045\ldots + 17.016045…+1 and states the rounded bound 8.01618.01618.0161. This entry rounds the computed value's last digit up to 8.0160468.0160468.016046, which is still below the paper's 8.01618.01618.0161 and above 8.016045…8.016045\ldots8.016045….

How the bound arises

Hata obtains simultaneous approximations to 1,π,log⁡21, \pi, \log 21,π,log2 from complex contour integrals of Legendre type, with extra arithmetic savings from primes that divide the coefficients to a predictable extent. His linear-form measure is 7.016045…7.016045\ldots7.016045…, giving μ(π)≤8.016045…\mu(\pi) \le 8.016045\ldotsμ(π)≤8.016045… (stated in the paper as 8.01618.01618.0161). It stood as the record for about fifteen years.

Significance

Each step down the list replaces Mahler's approximations with a sharper family. Formalizing 8.0160468.0160468.016046 would build reusable explicit machinery: integral constructions of rational approximations to π\piπ, bounds on their common denominators via prime-number estimates, and the standard lemma turning a sequence of good approximations into an irrationality-measure bound.

Selected references

  • M. Hata, Rational approximations to π\piπ and some other numbers, Acta Arith. 63 (1993), no. 4, 335–349. http://matwbn.icm.edu.pl/ksiazki/aa/aa63/aa6344.pdf
  • K. Mahler, On the approximation of π\piπ, Indag. Math. 15 (1953), 30–42.
  • F. Beukers, A rational approach to π\piπ, Nieuw Arch. Wiskd. (5) 1 (2000), 372–379.
  • Source table: https://teorth.github.io/optimizationproblems/constants/7a.html
11 thms3 active usersReviewed
🏆Completed
Captain: marwahaha

Mahler's irrationality bound for π: 42Research Paper

Motivation

The irrationality of π\piπ rules out an exact representation as a rational number. A quantitative question asks how closely rational numbers can approximate it as their denominators grow. The irrationality measure records the threshold exponent for exceptionally accurate rational approximations. This mission formalizes the historical upper-bound milestone μ(π)≤42\mu(\pi)\le42μ(π)≤42 listed as C7aC_{7a}C7a​ in the optimization constants project.

Mahler's 1953 paper establishes a stronger, explicit inequality in Theorem 1. The mission extracts its consequence for the irrationality measure and states that consequence through a shared Lean predicate. The number 42 is the chosen historical milestone; it is neither a claim about the exact value of the measure nor a claim to the strongest bound mentioned anywhere in Mahler's paper. Mahler, original p. 33.

Setting

Write p∈Zp\in\mathbb Zp∈Z for a numerator and q∈Nq\in\mathbb Nq∈N for a positive denominator. The approximation error is the real number ∣π−p/q∣|\pi-p/q|∣π−p/q∣. The exponent BBB describes an upper bound on the irrationality measure through an eventual lower bound on this error.

The shared predicate PiIrrationality.UpperBound B means that for every real ε>0\varepsilon>0ε>0, some natural-number threshold QQQ satisfies

1qB+ε<∣π−pq∣\frac{1}{q^{B+\varepsilon}}<\left|\pi-\frac pq\right|qB+ε1​<​π−qp​​

for every integer ppp and every natural number q>0q>0q>0 with Q≤qQ\le qQ≤q. The threshold can depend on ε\varepsilonε and on the chosen bound BBB; it cannot depend on the later choices of ppp or qqq. Numerators may be negative, zero, or positive. Fractions need not be in lowest terms. This is the epsilon characterization used in the definition of C7aC_{7a}C7a​.

Formalization targets

The goal is

μ(π)≤42,\mu(\pi)\le42,μ(π)≤42,

represented by PiIrrationality.UpperBound (42 : ℝ). Expanded, the target is

∀ε>0  ∃Q∈N  ∀p∈Z  ∀q∈N,q>0 ∧ Q≤q ⟹ 1q42+ε<∣π−pq∣.\forall\varepsilon>0\;\exists Q\in\mathbb N\;\forall p\in\mathbb Z\;\forall q\in\mathbb N,\quad q>0\ \land\ Q\le q\ \Longrightarrow\ \frac1{q^{42+\varepsilon}}<\left|\pi-\frac pq\right|.∀ε>0∃Q∈N∀p∈Z∀q∈N,q>0 ∧ Q≤q ⟹ q42+ε1​<​π−qp​​.

The mission contains one shared definition and one goal theorem. The definition introduces the proposition without asserting any bound. The theorem has no additional hypotheses, and its proof is intentionally left open. Future historical-bound missions can import the same definition and state a different numeric bound without changing the quantity being tracked.

Significance

A finite upper bound restricts the quality of rational approximations to π\piπ and excludes approximation at arbitrarily large exponents. The formal result would supply a reusable quantitative fact beyond the assertion that π\piπ is irrational. The published mathematical result is known; the work requested here is a machine-checked proof of the stated consequence.

The shared definition also fixes the meaning of all entries in the accompanying campaign. A smaller bound makes a stronger claim. A proof of a stronger entry may establish this historical goal as a consequence, provided it uses the same definition and no extra hypotheses. The mission remains mathematically valid after further improvements to the numerical bound.

Difficulty

A proof that π\piπ is irrational only establishes nonzero approximation errors. The target requires a uniform lower estimate over every numerator once the denominator passes a threshold. Checking finitely many rational approximations cannot establish the quantified conclusion. A formal development must control the dependence of its estimates and thresholds, and justify every passage between an analytic estimate and the final rational-approximation inequality.

The exact theorem from Mahler should be distinguished from this goal: his explicit uniform inequality is stronger than the eventual epsilon statement recorded here. A proof may pass through that uniform result, but the goal does not require a particular proof method, a particular threshold, or a separate treatment of every auxiliary theorem in the original paper.

Formalization scope

The circle constant is Mathlib's Real.pi. Absolute value, division, and exponentiation in the displayed inequality are operations on the real numbers; in particular, q42+εq^{42+\varepsilon}q42+ε is a real power. The denominator is explicitly positive, so division by zero cannot satisfy the premises. Allowing Q=0Q=0Q=0 does not remove the positivity requirement on qqq.

The definition is stored in Definitions.Def_PiIrrationality_UpperBound. The goal imports this definition instead of introducing another version of it. The conclusion is not placed among the theorem's assumptions. The definition contains no proof placeholder; the sole sorry is the open proof of the goal theorem. Contributions may establish supporting estimates or a complete proof while preserving these conventions.

Selected references

  • K. Mahler, On the approximation of π, Nederl. Akad. Wetensch. Proc. Ser. A 56 = Indag. Math. 15 (1953), 30–42, Theorem 1, p. 33. EMS reprint.
  • Optimization problems project, The irrationality measure of π, constant C7aC_{7a}C7a​: definition and historical bounds. Source page.
3 thms3 active usersReviewed
🏆Completed
ProbabilityQuantum InformationTheoretical Computer Science·Captain: mikedeng1

Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer 3: The Success Probability of Quantum Order FindingResearch Paper

Motivation

The security of the RSA cryptosystem rests on the assumed difficulty of factoring large integers, and the best known classical algorithms for factoring run in super-polynomial time. In 1994 Peter Shor showed that a quantum computer can factor an nnn-digit integer in time polynomial in nnn (Shor, SIAM J. Comput. 1997; conference version FOCS 1994). The algorithm has two parts. A classical reduction, due to Miller (1976), turns factoring into order finding: given xxx coprime to nnn, find the least r≥1r \ge 1r≥1 with xr≡1(modn)x^r \equiv 1 \pmod nxr≡1(modn). The quantum part solves order finding.

This mission formalizes the quantum part as Shor analyzes it in §5 of the journal paper: the construction of the quantum state, the probability of each measurement outcome, and the classical post-processing that reads rrr off the measured value. The paper's claim is that one run of this procedure returns rrr with probability at least φ(r)/3r\varphi(r)/3rφ(r)/3r.

Timeline:

  • 1976: Miller reduces factoring to order finding (with randomization).
  • 1985–1994: Deutsch, Bernstein–Vazirani and Simon give the quantum Fourier sampling ideas the algorithm builds on.
  • 1994: Shor's FOCS paper introduces the factoring and discrete logarithm algorithms.
  • 1997: the SIAM J. Comput. version gives the analysis formalized here, with qqq the power of 222 in [n2,2n2)[n^2, 2n^2)[n2,2n2).

Setting

Fix an integer n≥2n \ge 2n≥2 and an integer xxx coprime to nnn. Its order rrr is the least r≥1r \ge 1r≥1 with xr≡1(modn)x^r \equiv 1 \pmod nxr≡1(modn); since xxx is a unit, r≤φ(n)<nr \le \varphi(n) < nr≤φ(n)<n. Let q=2lq = 2^lq=2l be the power of 222 with n2≤q<2n2n^2 \le q < 2n^2n2≤q<2n2.

A quantum state on two registers, the first holding 0≤a<q0 \le a < q0≤a<q and the second a residue y∈Z/ny \in \mathbb{Z}/ny∈Z/n, is a complex vector ψ(a,y)\psi(a, y)ψ(a,y) indexed by the basis states ∣a,y⟩|a, y\rangle∣a,y⟩. Measuring it returns ∣a,y⟩|a, y\rangle∣a,y⟩ with probability ∣ψ(a,y)∣2|\psi(a, y)|^2∣ψ(a,y)∣2.

The Fourier matrix AqA_qAq​ is the q×qq \times qq×q matrix with entries (Aq)a,c=q−1/2exp⁡(2πiac/q)(A_q)_{a,c} = q^{-1/2}\exp(2\pi i a c/q)(Aq​)a,c​=q−1/2exp(2πiac/q), with rows indexing inputs and columns outputs. The algorithm

  1. prepares 1q1/2∑a=0q−1∣a⟩∣xa mod n⟩\frac{1}{q^{1/2}}\sum_{a=0}^{q-1}|a\rangle|x^a \bmod n\rangleq1/21​∑a=0q−1​∣a⟩∣xamodn⟩ (eq. (5.2)),
  2. applies AqA_qAq​ to the first register, obtaining 1q∑a,cexp⁡(2πiac/q)∣c⟩∣xa mod n⟩\frac1q\sum_{a,c}\exp(2\pi iac/q)|c\rangle|x^a \bmod n\rangleq1​∑a,c​exp(2πiac/q)∣c⟩∣xamodn⟩ (eq. (5.4)),
  3. measures, obtaining some ∣c,y⟩|c, y\rangle∣c,y⟩,
  4. rounds c/qc/qc/q to the nearest fraction with denominator smaller than nnn.

The observed ccc gives us rrr if some fraction with lowest-terms denominator below nnn is within 1/2q1/2q1/2q of c/qc/qc/q, and every such fraction has lowest-terms denominator exactly rrr. In the Lean development these objects are preFourierState, finalState, outcomeProb and yieldsOrder, in the namespace ShorAlgorithms.OrderFinding, and the shared definition ShorAlgorithms.Shared.fourierMatrix.

Formalization targets

Goal: success probability at least φ(r)/3r\varphi(r)/3rφ(r)/3r

For all sufficiently large nnn, with xxx, rrr and qqq as above,

Pr⁡[the observed c gives us r]  =  ∑c gives r ∑y∈Z/n∣Ψ(c,y)∣2  ≥  φ(r)3r,\Pr\bigl[\text{the observed } c \text{ gives us } r\bigr] \;=\; \sum_{c\ \text{gives}\ r}\ \sum_{y \in \mathbb{Z}/n} |\Psi(c, y)|^2 \;\ge\; \frac{\varphi(r)}{3r},Pr[the observed c gives us r]=c gives r∑​ y∈Z/n∑​∣Ψ(c,y)∣2≥3rφ(r)​,

where Ψ\PsiΨ is the state (5.4). The threshold on nnn is uniform in xxx and qqq; it is the paper's "for sufficiently large nnn" from the per-state bound.

Milestones

  1. Eqs. (5.5)–(5.6). For 0≤k<r0 \le k < r0≤k<r, the probability of ∣c,xk⟩|c, x^k\rangle∣c,xk⟩ equals ∣1q∑b=0⌊(q−k−1)/r⌋exp⁡(2πi(br+k)c/q)∣2\left|\frac1q\sum_{b=0}^{\lfloor (q-k-1)/r\rfloor}\exp(2\pi i(br+k)c/q)\right|^2​q1​∑b=0⌊(q−k−1)/r⌋​exp(2πi(br+k)c/q)​2.
  2. Eq. (5.11). For nnn past a threshold, every ∣c,xk⟩|c, x^k\rangle∣c,xk⟩ with −r/2≤rc−dq≤r/2-r/2 \le rc - dq \le r/2−r/2≤rc−dq≤r/2 for some integer ddd has probability at least 1/3r21/3r^21/3r2.
  3. Eq. (5.13). If n2≤qn^2 \le qn2≤q, at most one fraction with denominator below nnn lies within 1/2q1/2q1/2q of c/qc/qc/q.
  4. p. 1500. Such a fraction is a convergent of the continued fraction of c/qc/qc/q.
  5. p. 1501. At least φ(r)\varphi(r)φ(r) values of ccc are within 1/2q1/2q1/2q of some d/rd/rd/r with gcd⁡(d,r)=1\gcd(d, r) = 1gcd(d,r)=1; with the rrr distinct values of xkx^kxk this gives at least rφ(r)r\varphi(r)rφ(r) states ∣c,xk⟩|c, x^k\rangle∣c,xk⟩, and each such ccc gives us rrr.

Significance

The goal is the quantitative statement behind "order finding is in bounded-error quantum polynomial time": since φ(r)/r≥δ/log⁡log⁡r\varphi(r)/r \ge \delta/\log\log rφ(r)/r≥δ/loglogr for a constant δ\deltaδ (Hardy and Wright, Thm. 328), O(log⁡log⁡r)O(\log\log r)O(loglogr) repetitions find rrr with high probability, and Miller's reduction then factors nnn. Without the bound, the algorithm is a procedure with no guarantee.

The result is proved, in the paper and in textbooks (Nielsen and Chuang, 2000, §5.3), usually with a phase-estimation analysis rather than Shor's direct count. What this mission adds is a machine-checked proof of Shor's own argument, with his choice of qqq and his constants, starting from the state built by applying AqA_qAq​ to (5.2). Formal proofs of idealized versions exist elsewhere, for instance in the exact-period model where rrr divides qqq and the output is uniform on rrr peaks, but that model removes the approximation that the 1/3r21/3r^21/3r2 bound is about. Legendre's theorem on continued fractions is already on the platform (FamousTheorems.legendre_continued_fraction_theorem) and is included as a reference item.

Difficulty

The obvious route is to compute the output distribution in closed form. That works only when rrr divides qqq; here qqq is a power of 222 and rrr is arbitrary, so the amplitudes are geometric sums of ⌊(q−k−1)/r⌋+1\lfloor (q-k-1)/r\rfloor + 1⌊(q−k−1)/r⌋+1 terms whose phases do not cancel exactly. The per-state bound 1/3r21/3r^21/3r2 requires a lower bound on such a sum that is uniform in rrr, ccc and kkk, with error terms of order 1/q1/q1/q controlled against a main term of order 1/r21/r^21/r2. The constant 1/31/31/3 leaves only a small margin below the limiting value 4/π2≈0.4054/\pi^2 \approx 0.4054/π2≈0.405, so the errors must be bounded explicitly, not merely shown to vanish.

The second difficulty is the counting: distinct coprime numerators ddd must give distinct outcomes ccc in [0,q)[0, q)[0,q), and each good ccc must determine rrr uniquely, which uses r<nr < nr<n and n2≤qn^2 \le qn2≤q.

Formalization scope

Conventions the statements commit to:

  • States are functions Fin q × ZMod n → ℂ; the matrix convention is row = input, so applying AqA_qAq​ to the first register gives the amplitude ∑aψ(a,y)(Aq)a,c\sum_a \psi(a, y)(A_q)_{a,c}∑a​ψ(a,y)(Aq​)a,c​ at (c,y)(c, y)(c,y).
  • The final state is built by applying AqA_qAq​ to the state (5.2); the closed forms (5.5) and (5.6) are theorems, not definitions. No normalization hypothesis is assumed.
  • Probabilities are squared moduli; the probability of the event "ccc gives us rrr" sums over all y∈Z/ny \in \mathbb{Z}/ny∈Z/n, which is exact because yyy that are not powers of xxx have probability zero.
  • xxx is a natural number with gcd⁡(x,n)=1\gcd(x, n) = 1gcd(x,n)=1; rrr is orderOf (x : ZMod n). qqq enters through the three hypotheses q=2lq = 2^lq=2l, n2≤qn^2 \le qn2≤q, q<2n2q < 2n^2q<2n2, not through a function of nnn.
  • Fractions are rationals, and "in lowest terms" is Rat.den.
  • Thresholds "for sufficiently large nnn" are ∃N, ∀n≥N\exists N,\ \forall n \ge N∃N, ∀n≥N, with NNN quantified before xxx, qqq, ccc and kkk.
  • Condition (5.11) is stated in its equivalent form (5.12), with an integer ddd.
  • Printed slip. Eq. (5.13)'s justification says "Because q>n2q > n^2q>n2", but qqq was chosen with n2≤qn^2 \le qn2≤q, and q=n2q = n^2q=n2 when nnn is a power of 222. The uniqueness claim holds under n2≤qn^2 \le qn2≤q, and that is what is stated.

Typing the closed form (5.4)–(5.6) in as the definition of the final state would make milestone 1 trivial and hide whether the probability model is the paper's; the definitions exclude this by construction.

Not stated: the polynomial running time of any step, the O(log⁡log⁡r)O(\log\log r)O(loglogr) repetition count (no explicit constant), the reversible modular exponentiation of §3, and the post-processing heuristics on p. 1501. Needed infrastructure: bounds on geometric exponential sums, Euler's totient, Diophantine approximation by fractions with bounded denominator, and Mathlib's continued fractions. Lemmas on geometric sums of roots of unity and on the order of units mod nnn are reusable in the companion discrete logarithm mission.

Selected references

  • P. W. Shor, Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer, SIAM J. Comput. 26(5):1484–1509, 1997. https://doi.org/10.1137/S0097539795293172
  • P. W. Shor, Algorithms for quantum computation: discrete logarithms and factoring, Proc. 35th FOCS, 1994. https://doi.org/10.1109/SFCS.1994.365700
  • G. L. Miller, Riemann's hypothesis and tests for primality, J. Comput. System Sci. 13(3):300–317, 1976. https://doi.org/10.1016/S0022-0000(76)80043-8
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford, 1979 (Ch. X, continued fractions; Thm. 328).
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge, 2000. https://doi.org/10.1017/CBO9780511976667
12 thms3 active usersReviewed
🏆Completed
CombinatoricsGraph TheoryTheoretical Computer Science·Captain: mikedeng1

Fast Algorithms for Finding Nearest Common Ancestors II: Nearest Common Ancestors in a Complete Binary Tree by Symmetric-Order ArithmeticResearch Paper

Motivation

The nearest common ancestor (nca) problem asks, for a fixed rooted tree and a sequence of vertex pairs (v,w)(v, w)(v,w), for the deepest vertex that is an ancestor of both. It is a basic step in suffix-tree string algorithms and is equivalent to range-minimum queries (Bender, Farach-Colton, 2000). Harel and Tarjan, Fast Algorithms for Finding Nearest Common Ancestors, SIAM J. Comput. 13 (1984) 338–355, gave the first algorithm answering each query on a static tree in constant time on a random-access machine after linear preprocessing.

Their construction reduces the general problem to the case of a complete binary tree, where §3 of the paper shows that nca queries can be answered "by direct calculation" on vertex numbers: multiplication, division, powers of two, the base-two logarithm and bitwise exclusive or. The later simplification of Schieber and Vishkin (1988) is built on the same in-order numbering of a complete binary tree. This mission formalizes that arithmetic core.

Timeline, as reviewed in the paper's §1 (pp. 338–340):

  • 1976: Aho, Hopcroft and Ullman (SIAM J. Comput. 5) give an O(n+mα(m+n,n))O(n + m\alpha(m+n, n))O(n+mα(m+n,n))-time off-line algorithm on a pointer machine, and for static trees a random-access algorithm with O(nlog⁡log⁡n)O(n \log\log n)O(nloglogn) preprocessing and O(log⁡log⁡n)O(\log\log n)O(loglogn) time per query.
  • 1976: van Leeuwen (unpublished report) gives an O(n+mlog⁡log⁡n)O(n + m \log\log n)O(n+mloglogn)-time algorithm for linking roots and static trees that runs on a pointer machine in O(n)O(n)O(n) space.
  • 1980: Harel (Proc. 21st FOCS) gives a preliminary version of the paper's results.
  • 1984: Harel and Tarjan prove that pointer machines need Ω(log⁡log⁡n)\Omega(\log\log n)Ω(loglogn) time per query on static trees (Theorem 1), and give the O(n)O(n)O(n)-preprocessing, O(1)O(1)O(1)-query random-access algorithm whose base case is the subject of this mission.

Setting

Fix d≥0d \ge 0d≥0 and let TTT be the complete binary tree of depth ddd. A vertex is identified with the path from the root to it, a word of at most ddd left or right turns; the root is the empty word and TTT has n=2d+1−1n = 2^{d+1} - 1n=2d+1−1 vertices. Following the paper's Appendix (pp. 354–355):

  • www is an ancestor of vvv (vvv a descendant of www) if the word www is a prefix of the word vvv; every vertex is its own ancestor. vvv and www are unrelated if neither is an ancestor of the other.
  • The depth of vvv is its distance to the root; its height h(v)h(v)h(v) is the length of the longest path from a leaf to vvv, which in TTT is d−depth⁡(v)d - \operatorname{depth}(v)d−depth(v).
  • nca⁡(v,w)\operatorname{nca}(v, w)nca(v,w) is the vertex of greatest depth that is an ancestor of both: the longest common prefix.

The vertices of TTT are numbered from 111 to nnn in symmetric order (in-order): at every vertex, first the left subtree, then the vertex, then the right subtree. sym(v)\mathrm{sym}(v)sym(v) is the number of vvv and sym−1(i)\mathrm{sym}^{-1}(i)sym−1(i) the vertex numbered iii. For d=4d = 4d=4 (Fig. 1 of the paper) the root is 161616, its children 888 and 242424, and the leaves 1,3,5,…,311, 3, 5, \dots, 311,3,5,…,31. i⊕ji \oplus ji⊕j denotes bitwise exclusive or and lg⁡\lglg the base-two logarithm.

Two procedures of §3 use only numbers, heights and ddd:

  • the nca depth algorithm: return d−h(v)d - h(v)d−h(v) if sym(w)∈[sym(v)−2h(v)+1,sym(v)+2h(v)−1]\mathrm{sym}(w) \in [\mathrm{sym}(v) - 2^{h(v)} + 1, \mathrm{sym}(v) + 2^{h(v)} - 1]sym(w)∈[sym(v)−2h(v)+1,sym(v)+2h(v)−1]; else d−h(w)d - h(w)d−h(w) if the same holds with v,wv, wv,w exchanged; else d−⌊lg⁡(sym(v)⊕sym(w))⌋d - \lfloor \lg(\mathrm{sym}(v) \oplus \mathrm{sym}(w)) \rfloord−⌊lg(sym(v)⊕sym(w))⌋;
  • the depth algorithm: given vvv and a depth d2≤depth⁡(v)d_2 \le \operatorname{depth}(v)d2​≤depth(v), with h=d−d2h = d - d_2h=d−d2​, return sym−1(2h+1⌊sym(v)/2h+1⌋+2h)\mathrm{sym}^{-1}\bigl(2^{h+1}\lfloor \mathrm{sym}(v)/2^{h+1}\rfloor + 2^h\bigr)sym−1(2h+1⌊sym(v)/2h+1⌋+2h).

Formalization targets

Goal: the nca algorithm is correct

The algorithm to compute nca⁡(v,w)\operatorname{nca}(v,w)nca(v,w) (p. 342) runs the nca depth algorithm to obtain d0d_0d0​ and then the depth algorithm on (v,d0)(v, d_0)(v,d0​). The goal states that it returns the nearest common ancestor: for all vertices v,wv, wv,w of TTT, with d0d_0d0​ the output of the nca depth algorithm and h=d−d0h = d - d_0h=d−d0​,

sym(nca⁡(v,w))=2h+1⌊sym(v)2h+1⌋+2h.\mathrm{sym}(\operatorname{nca}(v,w)) = 2^{h+1}\left\lfloor \frac{\mathrm{sym}(v)}{2^{h+1}} \right\rfloor + 2^h .sym(nca(v,w))=2h+1⌊2h+1sym(v)​⌋+2h.

Milestones

In the order the paper uses them:

  1. Numbers at height hhh (p. 341): the vertices of height hhh are numbered 2h,3⋅2h,5⋅2h,…2^h, 3\cdot 2^h, 5\cdot 2^h, \dots2h,3⋅2h,5⋅2h,… from left to right.
  2. Lemma 1: h(v)h(v)h(v) is the largest hhh with 2h∣sym(v)2^h \mid \mathrm{sym}(v)2h∣sym(v).
  3. Lemma 2: the descendants of vvv are the vertices numbered in [sym(v)−2h(v)+1,sym(v)+2h(v)−1][\mathrm{sym}(v) - 2^{h(v)} + 1, \mathrm{sym}(v) + 2^{h(v)} - 1][sym(v)−2h(v)+1,sym(v)+2h(v)−1].
  4. Lemma 3: for a height h≥h(v)h \ge h(v)h≥h(v), the height-hhh ancestor of vvv has number 2h+1⌊sym(v)/2h+1⌋+2h2^{h+1}\lfloor \mathrm{sym}(v)/2^{h+1}\rfloor + 2^h2h+1⌊sym(v)/2h+1⌋+2h.
  5. Lemma 4: for unrelated v,wv, wv,w,
h(nca⁡(v,w))=⌊lg⁡(sym(v)⊕sym(w))⌋.h(\operatorname{nca}(v,w)) = \lfloor \lg(\mathrm{sym}(v) \oplus \mathrm{sym}(w)) \rfloor .h(nca(v,w))=⌊lg(sym(v)⊕sym(w))⌋.
  1. The nca depth algorithm returns depth⁡(nca⁡(v,w))\operatorname{depth}(\operatorname{nca}(v,w))depth(nca(v,w)).
  2. The depth algorithm returns the number of the depth-d2d_2d2​ ancestor of vvv.

Two supporting statements pin the definitions to the paper: sym\mathrm{sym}sym is a bijection onto {1,…,2d+1−1}\{1, \dots, 2^{d+1} - 1\}{1,…,2d+1−1}, and the longest common prefix is the deepest common ancestor.

Significance

The constant-time nca computation on complete binary trees is the base case of the whole paper: §§4–5 embed an arbitrary tree into a moderately sized complete binary tree through a compressed tree and a balanced binary tree, and every query ends with the arithmetic of §3. The same idea, that in-order numbers encode ancestry in their low-order bits, underlies the Schieber–Vishkin algorithm. Lemma 1 identifies the height with the 2-adic valuation of the number, and Lemma 4 identifies the nca height with the position of the highest differing bit.

The results are proved in the paper, with the proofs left as "easy to verify". No machine-checked version of this numbering or of these four lemmas is known to exist in Mathlib or on this platform. A formal development supplies proofs of the four lemmas and the two algorithms, and a reusable library connecting in-order ranks of a complete binary tree to binary arithmetic (Nat.log, bitwise xor, 2-adic valuation).

Difficulty

The numbering is defined by a traversal order, while the lemmas speak about divisibility, floor division and exclusive or. The work lies in connecting the rank of a vertex in symmetric order to its closed form (2j+1)⋅2h(v)(2j+1)\cdot 2^{h(v)}(2j+1)⋅2h(v), where jjj is its left-to-right position. That counting argument sums the sizes of the subtrees that precede vvv and is where most of the effort goes. Lemma 4 then needs the observation that two unrelated numbers agree in all bits above the height of their nca and differ in the bit at that height. This is a statement about Nat.testBit of the exclusive or, and it fails for related vertices. The algorithm statements add a case analysis whose first two cases overlap when v=wv = wv=w.

Formalization scope

  • A vertex of the tree of depth ddd is a List Bool of length at most ddd (false = left). Ancestry is the prefix relation, nca⁡\operatorname{nca}nca the longest common prefix, depth the length, and height d−lengthd - \text{length}d−length. None of these structural notions uses the numbering.
  • sym(v)\mathrm{sym}(v)sym(v) is the number of vertices whose in-order sort key is lexicographically at most that of vvv. The key is the path with left ↦0\mapsto 0↦0, right ↦2\mapsto 2↦2, followed by 111. The numbering is not defined by the closed form or by a recursion on numbers: a definition of that kind would make the height-hhh numbering and Lemma 1 immediate and move the content of the mission into an uncheckable definition.
  • ⌊lg⁡x⌋\lfloor \lg x \rfloor⌊lgx⌋ is Nat.log 2 x, which agrees for x≥1x \ge 1x≥1. ⊕\oplus⊕ is ^^^ on N\mathbb NN, and floor division is / on N\mathbb NN.
  • Interval tests a∈[b−c+1,b+c−1]a \in [b - c + 1, b + c - 1]a∈[b−c+1,b+c−1] are written additively as b+1≤a+cb + 1 \le a + cb+1≤a+c and a+1≤b+ca + 1 \le b + ca+1≤b+c. The subtractions d−h(v)d - h(v)d−h(v) and d−d2d - d_2d−d2​ never truncate for heights and depths of vertices.
  • Lemma 3 states explicitly that h≤dh \le dh≤d ("hhh is a height") and that the ancestor exists. The depth algorithm assumes d2≤depth⁡(v)d_2 \le \operatorname{depth}(v)d2​≤depth(v), as printed.
  • sym−1\mathrm{sym}^{-1}sym−1 is not defined as a function. The goal and the depth algorithm state that a vertex has the computed number if and only if it is the nearest common ancestor (respectively the ancestor at depth d2d_2d2​), which says that sym−1\mathrm{sym}^{-1}sym−1 of that number is that vertex.
  • The O(1)O(1)O(1) time bounds are not formalized, since the random-access machine model is out of scope.

Proofs of any milestone are welcome.

Selected references

  • D. Harel and R. E. Tarjan, Fast Algorithms for Finding Nearest Common Ancestors, SIAM J. Comput. 13(2) (1984), 338–355. https://doi.org/10.1137/0213024
  • A. V. Aho, J. E. Hopcroft and J. D. Ullman, On Finding Lowest Common Ancestors in Trees, SIAM J. Comput. 5(1) (1976), 115–132. https://doi.org/10.1137/0205011
  • B. Schieber and U. Vishkin, On Finding Lowest Common Ancestors: Simplification and Parallelization, SIAM J. Comput. 17(6) (1988), 1253–1262. https://doi.org/10.1137/0217079
  • M. A. Bender and M. Farach-Colton, The LCA Problem Revisited, LATIN 2000, LNCS 1776, 88–94. https://doi.org/10.1007/10719839_9
11 thms3 active usersReviewed
Discrete GeometryOperations Research·Captain: mikedeng1

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

Motivation

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

Timeline.

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

Setting

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

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

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

Formalization targets

Goal: Theorem (5.5), corrected reading

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

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

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

Milestones

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

Significance

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

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

Difficulty

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

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

Formalization scope

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

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

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

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

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

Selected references

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

Erdős Problem 30: Sidon sets in {1,…,N} have size √N + O(N^ε)Open Problem

Motivation

A set of integers is a Sidon set if all of its pairwise sums a+ba+ba+b (a≤ba\le ba≤b) are different. Sidon, in connection with Fourier analysis, asked how dense such sets can be, and the question became one of the standard problems of additive combinatorics. Let

h(N)=max⁡{∣A∣:A⊆{1,…,N}, A Sidon}.h(N)=\max\{|A| : A\subseteq\{1,\dots,N\},\ A \text{ Sidon}\}.h(N)=max{∣A∣:A⊆{1,…,N}, A Sidon}.

A counting argument shows h(N)≤(1+o(1))2Nh(N)\le (1+o(1))\sqrt{2N}h(N)≤(1+o(1))2N​, and the true order was settled early: h(N)∼Nh(N)\sim\sqrt Nh(N)∼N​. What remains open is the size of the error term h(N)−Nh(N)-\sqrt Nh(N)−N​. Erdős and Turán asked whether it is smaller than every power of NNN; Erdős offered $1000 for this problem (Erdős Problem #30), and it is also Problem 31 on Green's list of open problems and problem C9 in Guy's Unsolved Problems in Number Theory.

Timeline.

  • 1938 — Singer constructs, for every prime power qqq, a set of q+1q+1q+1 residues modulo q2+q+1q^2+q+1q2+q+1 with all differences distinct. Combined with the density of primes this gives h(N)≥(1−o(1))Nh(N)\ge(1-o(1))\sqrt Nh(N)≥(1−o(1))N​ (Singer 1938).
  • 1941 — Erdős and Turán prove h(N)≤N1/2+O(N1/4)h(N)\le N^{1/2}+O(N^{1/4})h(N)≤N1/2+O(N1/4) (Erdős–Turán 1941).
  • 1969 — Lindström gives an alternative proof with the explicit bound h(N)≤N1/2+N1/4+1h(N)\le N^{1/2}+N^{1/4}+1h(N)≤N1/2+N1/4+1 (Lindström 1969).
  • 2021 — Balogh, Füredi and Roy lower the constant: h(N)≤N1/2+0.998N1/4h(N)\le N^{1/2}+0.998N^{1/4}h(N)≤N1/2+0.998N1/4 for large NNN (arXiv:2103.15850).
  • 2022 — O'Bryant: h(N)≤N1/2+0.99703N1/4h(N)\le N^{1/2}+0.99703N^{1/4}h(N)≤N1/2+0.99703N1/4 for large NNN (arXiv:2207.07800).
  • 2023 — Carter, Hunter and O'Bryant: h(N)≤N1/2+0.98183N1/4+O(1)h(N)\le N^{1/2}+0.98183N^{1/4}+O(1)h(N)≤N1/2+0.98183N1/4+O(1), with substantial computer assistance (arXiv:2310.20032).

No upper bound with an error exponent below 1/41/41/4 is known, and no lower bound of the form h(N)≥N−O(Nε)h(N)\ge\sqrt N-O(N^{\varepsilon})h(N)≥N​−O(Nε) for every ε>0\varepsilon>0ε>0 is known either.

Setting

A set AAA in an additive commutative monoid is Sidon if for all i1,j1,i2,j2∈Ai_1,j_1,i_2,j_2\in Ai1​,j1​,i2​,j2​∈A,

i1+i2=j1+j2 ⟹ (i1=j1∧i2=j2) ∨ (i1=j2∧i2=j1).i_1+i_2=j_1+j_2\ \Longrightarrow\ (i_1=j_1\wedge i_2=j_2)\ \vee\ (i_1=j_2\wedge i_2=j_1).i1​+i2​=j1​+j2​ ⟹ (i1​=j1​∧i2​=j2​) ∨ (i1​=j2​∧i2​=j1​).

For a finite set XXX, maxSidon⁡(X)\operatorname{maxSidon}(X)maxSidon(X) is the largest size of a Sidon subset of XXX (the empty set is Sidon, so this is well defined), and

h(N)=maxSidon⁡({1,2,…,N}),h(0)=0.h(N)=\operatorname{maxSidon}(\{1,2,\dots,N\}),\qquad h(0)=0.h(N)=maxSidon({1,2,…,N}),h(0)=0.

The first values are h(1),…,h(15)=1,2,2,3,3,3,4,4,4,4,4,5,5,5,5h(1),\dots,h(15)=1,2,2,3,3,3,4,4,4,4,4,5,5,5,5h(1),…,h(15)=1,2,2,3,3,3,4,4,4,4,4,5,5,5,5 (OEIS A143824).

Formalization targets

Goal (Erdős Problem #30)

∀ε>0:h(N)−N=O ⁣(Nε)(N→∞).\forall\varepsilon>0:\qquad h(N)-\sqrt N = O\!\left(N^{\varepsilon}\right)\quad(N\to\infty).∀ε>0:h(N)−N​=O(Nε)(N→∞).

This is a two-sided statement: it asks both for an upper bound h(N)≤N+CεNεh(N)\le\sqrt N+C_\varepsilon N^\varepsilonh(N)≤N​+Cε​Nε and for a matching lower bound h(N)≥N−CεNεh(N)\ge\sqrt N-C_\varepsilon N^\varepsilonh(N)≥N​−Cε​Nε for large NNN. Erdős asked it as a yes/no question; the goal fixes the conjectured answer yes, so a disproof on the platform settles the question negatively.

Milestones (known results, weakest to strongest)

  1. Singer's construction: h(q2+q+1)≥q+1h(q^2+q+1)\ge q+1h(q2+q+1)≥q+1 for every prime power qqq.
  2. Singer's lower bound: h(N)≥(1−ε)Nh(N)\ge(1-\varepsilon)\sqrt Nh(N)≥(1−ε)N​ for every ε>0\varepsilon>0ε>0 and all large NNN.
  3. Erdős–Turán / Lindström: h(N)≤N+N1/4+1h(N)\le\sqrt N+N^{1/4}+1h(N)≤N​+N1/4+1 for all NNN.
  4. Balogh–Füredi–Roy: h(N)≤N+0.998N1/4h(N)\le\sqrt N+0.998N^{1/4}h(N)≤N​+0.998N1/4 for all large NNN.
  5. O'Bryant: h(N)≤N+0.99703N1/4h(N)\le\sqrt N+0.99703N^{1/4}h(N)≤N​+0.99703N1/4 for all large NNN.
  6. Carter–Hunter–O'Bryant: h(N)≤N+0.98183N1/4+Ch(N)\le\sqrt N+0.98183N^{1/4}+Ch(N)≤N​+0.98183N1/4+C for an absolute constant CCC.

Significance

The result itself. An affirmative answer would pin h(N)h(N)h(N) down to N\sqrt NN​ up to a sub-polynomial error, in both directions; Erdős even speculated that h(N)=N+O(1)h(N)=\sqrt N+O(1)h(N)=N​+O(1) might hold, while remarking that this is perhaps too optimistic. A negative answer would show that the Singer-type constructions or the counting upper bounds are off by a power of NNN. Either answer would be the first change in the exponent of the error term since 1941.

Formalizing it. The goal is open. All milestones are published theorems. The platform already contains weaker related results in other formalizations (for example the order-of-magnitude bounds cN≤max⁡∣A∣≤2N+1c\sqrt N\le \max|A|\le\sqrt{2N}+1cN​≤max∣A∣≤2N​+1 for Sidon subsets of an initial segment, and the Erdős–Turán construction); the sharp bounds listed as milestones are not stated there for this hhh. The upper bounds of Balogh–Füredi–Roy and O'Bryant are elementary but delicate optimizations, and the Carter–Hunter–O'Bryant bound relies on a large computation, so formalizing it is a substantial verification task in its own right.

Difficulty

For the upper bound, every known argument counts differences a−a′a-a'a−a′ in short windows and loses at the scale N1/4N^{1/4}N1/4; improvements since 1941 only change the constant in front of N1/4N^{1/4}N1/4. For the lower bound, the constructions (Singer, Bose, Ruzsa) produce Sidon sets of size about p\sqrt pp​ in a modulus ppp of size about NNN, and the loss comes from the gap between NNN and the nearest admissible modulus; bringing it below NεN^\varepsilonNε requires either new constructions or information about primes in very short intervals that is far beyond current knowledge.

Formalization scope

Sidon sets are formalized for sets in an arbitrary additive commutative monoid, with the definition, the decidability instance and maxSidon⁡\operatorname{maxSidon}maxSidon transcribed from the formal-conjectures library (definitions IsSidon, Finset.maxSidonSubsetCard, and Erdos30.h in FormalConjectures/ErdosProblems/30.lean), placed in the namespace Erdos30. The value h(N)h(N)h(N) is a natural number cast to R\mathbb RR; ⋅\sqrt{\cdot}⋅​ is the real square root and NεN^{\varepsilon}Nε, N1/4N^{1/4}N1/4 are real powers of N≥0N\ge 0N≥0. The goal's O(⋅)O(\cdot)O(⋅) is Mathlib's Asymptotics.IsBigO along atTop on N\mathbb NN. The goal is not trivialized by any junk value: hhh is a genuine finite maximum, and the O(⋅)O(\cdot)O(⋅) statement concerns all large NNN.

Useful infrastructure: basic lemmas on Sidon sets (hereditary under subsets, translation invariance, distinct differences), finite projective geometry or Bose's construction for the lower bounds, and prime gaps (Bertrand's postulate suffices for h(N)≥cNh(N)\ge c\sqrt Nh(N)≥cN​ with c<1c<1c<1; a prime number theorem in short intervals is needed for 1−o(1)1-o(1)1−o(1)). Contributions of reusable Sidon-set lemmas are welcome.

Selected references

  • J. Singer, A theorem in finite projective geometry and some applications to number theory, Trans. Amer. Math. Soc. 43 (1938), 377–385. https://doi.org/10.1090/S0002-9947-1938-1501951-4
  • P. Erdős and P. Turán, On a problem of Sidon in additive number theory, and on some related problems, J. London Math. Soc. 16 (1941), 212–215. https://doi.org/10.1112/jlms/s1-16.4.212
  • B. Lindström, An inequality for B2B_2B2​-sequences, J. Combin. Theory 6 (1969), 211–212. https://doi.org/10.1016/S0021-9800(69)80124-9
  • J. Balogh, Z. Füredi and S. Roy, An upper bound on the size of Sidon sets, Amer. Math. Monthly (2023). https://arxiv.org/abs/2103.15850
  • K. O'Bryant, On the size of finite Sidon sets (2022). https://arxiv.org/abs/2207.07800
  • D. Carter, Z. Hunter and K. O'Bryant, On the diameter of finite Sidon sets (2023). https://arxiv.org/abs/2310.20032
  • K. O'Bryant, A complete annotated bibliography of work related to Sidon sequences, Electron. J. Combin. DS11 (2004). https://arxiv.org/abs/math/0407117
  • T. F. Bloom, Erdős Problem #30, https://www.erdosproblems.com/30
  • Google DeepMind, formal-conjectures, FormalConjectures/ErdosProblems/30.lean. https://github.com/google-deepmind/formal-conjectures
8 thms3 active usersReviewed
Combinatorics·Captain: Lucas

Erdős Problem 3: arithmetic progressions in sets with divergent reciprocal sumOpen Problem

Motivation

Which sets of positive integers are forced to contain long arithmetic progressions? Van der Waerden (1927) showed that in any finite colouring of N\mathbb NN some colour class does; Erdős and Turán (1936) asked for a density version, which became Szemerédi's theorem. Erdős then proposed the strongest natural size condition: divergence of the reciprocal sum. Erdős Problem #3 (erdosproblems.com/3) asks whether every A⊆NA\subseteq\mathbb NA⊆N with ∑n∈A1/n=∞\sum_{n\in A}1/n=\infty∑n∈A​1/n=∞ contains arbitrarily long arithmetic progressions. Erdős attached one of his largest prizes to it. The primes are the motivating example: ∑p1/p=∞\sum_p 1/p=\infty∑p​1/p=∞, so a positive answer would contain the Green–Tao theorem.

Timeline.

  • 1936 — Erdős and Turán conjecture that sets of positive density contain arbitrarily long progressions.
  • 1953 — Roth proves the case k=3k=3k=3 of the density conjecture by Fourier analysis.
  • 1975 — Szemerédi proves the density conjecture for all kkk.
  • 2001 — Gowers gives the first quantitative bounds for all kkk: rk(N)≪N/(log⁡log⁡N)ckr_k(N)\ll N/(\log\log N)^{c_k}rk​(N)≪N/(loglogN)ck​.
  • 2008 — Green and Tao prove that the primes contain arbitrarily long progressions.
  • 2020 — Bloom and Sisask prove r3(N)≪N/(log⁡N)1+cr_3(N)\ll N/(\log N)^{1+c}r3​(N)≪N/(logN)1+c, which settles the case k=3k=3k=3 of Erdős Problem #3.
  • 2023 — Kelley and Meka prove r3(N)≤Nexp⁡(−c(log⁡N)1/12)r_3(N)\le N\exp(-c(\log N)^{1/12})r3​(N)≤Nexp(−c(logN)1/12).
  • 2024 — Leng, Sah and Sawhney prove rk(N)≤Nexp⁡(−(log⁡log⁡N)ck)r_k(N)\le N\exp(-(\log\log N)^{c_k})rk​(N)≤Nexp(−(loglogN)ck​) for every k≥5k\ge5k≥5.

The problem is open for every k≥4k\ge4k≥4.

Setting

A set S⊆NS\subseteq\mathbb NS⊆N is an arithmetic progression of length kkk if ∣S∣=k|S|=k∣S∣=k and S={a,a+d,…,a+(k−1)d}S=\{a,a+d,\dots,a+(k-1)d\}S={a,a+d,…,a+(k−1)d} for some a,d∈Na,d\in\mathbb Na,d∈N (for k≥2k\ge2k≥2 the size condition forces d>0d>0d>0). For k,N∈Nk,N\in\mathbb Nk,N∈N, rk(N)r_k(N)rk​(N) denotes the largest size of a subset of {1,…,N}\{1,\dots,N\}{1,…,N} containing no arithmetic progression of length kkk. A set AAA has divergent reciprocal sum if ∑n∈A1/n=∞\sum_{n\in A}1/n=\infty∑n∈A​1/n=∞.

Formalization targets

Goal (Erdős Problem #3)

For every A⊆NA\subseteq\mathbb NA⊆N,

∑n∈A1n=∞ ⟹ A contains arithmetic progressions of arbitrarily large length.\sum_{n\in A}\frac1n=\infty\ \Longrightarrow\ A\ \text{contains arithmetic progressions of arbitrarily large length}.n∈A∑​n1​=∞ ⟹ A contains arithmetic progressions of arbitrarily large length.

This is the formal-conjectures statement erdos_3 with its answer(sorry) instantiated to the conjectured answer yes. A disproof of the goal on the platform settles the problem negatively.

Milestones

  • Szemerédi's theorem for sets of positive upper density (the density case), and the existing platform statement rk(N)=o(N)r_k(N)=o(N)rk​(N)=o(N).
  • The Green–Tao theorem (the case A=A=A= primes).
  • The Bloom–Sisask bound on r3(N)r_3(N)r3​(N) and its corollary, the case k=3k=3k=3 of the goal; the Kelley–Meka bound (existing platform statement).
  • The Leng–Sah–Sawhney bound for k≥5k\ge5k≥5.
  • The partial-summation reduction: bounds rk(N)≤N/(log⁡N)1+ckr_k(N)\le N/(\log N)^{1+c_k}rk​(N)≤N/(logN)1+ck​ for all k≥3k\ge3k≥3 imply the goal.

Significance

A positive answer would be a common strengthening of Szemerédi's theorem and the Green–Tao theorem, obtained from a single size condition with no arithmetic structure. Through the reduction milestone, it is closely tied to the quantitative theory of rk(N)r_k(N)rk​(N): bounds of the shape N/(log⁡N)1+cN/(\log N)^{1+c}N/(logN)1+c for every kkk would suffice. Formalizing the milestones would also give reusable Lean statements of Szemerédi-type theorems in a common language.

Difficulty

Divergence of ∑1/n\sum 1/n∑1/n is a very weak condition: such sets can have density zero, and the natural approach through rk(N)r_k(N)rk​(N) requires bounds just past N/log⁡NN/\log NN/logN. For k=3k=3k=3 this barrier was only broken in 2020. For k≥4k\ge4k≥4 the best known bounds (Leng–Sah–Sawhney) save only a power of log⁡log⁡N\log\log NloglogN in the exponent, far from what is needed. The Green–Tao method uses pseudorandom majorants specific to the primes and does not apply to arbitrary sets.

Formalization scope

All statements import the published definition file Erdos142Basic, which reproduces the formal-conjectures definitions IsAPOfLengthWith, IsAPOfLength and the counting function r k N (over {1,…,N}\{1,\dots,N\}{1,…,N}). The reciprocal-sum hypothesis is ¬ Summable (fun a : A ↦ 1 / (a : ℝ)); the element 000, if present, contributes 1/0=01/0=01/0=0. "Arbitrarily long" is written as ∃ᶠ k in atTop, which is equivalent to "every length" because sub-progressions of progressions are progressions. Bounds stated in the literature with ≪\ll≪ are written without a multiplicative constant and with "for all sufficiently large NNN"; the constant can be absorbed into the exponent. Contributions formalizing partial summation over sets of naturals and the equivalence of the "frequently" and "for every kkk" forms are welcome.

Selected references

  • P. Erdős and P. Turán, On some sequences of integers, J. London Math. Soc. 11 (1936).
  • K. F. Roth, On certain sets of integers, J. London Math. Soc. 28 (1953).
  • E. Szemerédi, On sets of integers containing no k elements in arithmetic progression, Acta Arith. 27 (1975).
  • W. T. Gowers, A new proof of Szemerédi's theorem, Geom. Funct. Anal. 11 (2001).
  • B. Green and T. Tao, The primes contain arbitrarily long arithmetic progressions, Ann. of Math. 167 (2008).
  • T. F. Bloom and O. Sisask, Breaking the logarithmic barrier in Roth's theorem on arithmetic progressions, arXiv:2007.03528 (2020).
  • Z. Kelley and R. Meka, Strong bounds for 3-progressions, FOCS 2023, arXiv:2302.05537.
  • J. Leng, A. Sah and M. Sawhney, Improved bounds for Szemerédi's theorem, arXiv:2402.17995 (2024).
  • T. F. Bloom, Erdős Problem #3, https://www.erdosproblems.com/3
16 thms3 active usersReviewed
Captain: Lucas

Erdős Problem 1210: reciprocal gaps of pairwise coprime setsOpen Problem

Motivation

A set AAA of positive integers is pairwise coprime if gcd⁡(a,b)=1\gcd(a,b)=1gcd(a,b)=1 for all distinct a,b∈Aa,b\in Aa,b∈A. The primes are the model example, and a recurring theme in Erdős's combinatorial number theory is that pairwise coprime sets cannot do much better than the primes on natural additive or harmonic statistics. Erdős Problem 1210 asks for a sharp version of this principle for the harmonic weight 1/(n−a)1/(n-a)1/(n−a), which measures how densely a coprime set can crowd the point nnn from below.

Timeline.

  • 1977. In [Er77c, p.64] Erdős posed a question about the primes q1<⋯<qkq_1<\dots<q_kq1​<⋯<qk​ in an interval (n,m](n,m](n,m]: is ∑i1/(qi−n)<∑p<m−n1/p+O(1)\sum_i 1/(q_i-n)<\sum_{p<m-n}1/p+O(1)∑i​1/(qi​−n)<∑p<m−n​1/p+O(1)?
  • 1980. In [Er80, p.112] he wrote that he had "not stated [this] quite correctly" in [Er77c] and posed the question for arbitrary pairwise coprime sets A⊆[1,n)A\subseteq[1,n)A⊆[1,n), which is the form recorded as Problem 1210.
  • 2026. On the erdosproblems.com forum, a reduction to a counting bound for A∩[n−x,n)A\cap[n-x,n)A∩[n−x,n) was suggested; it was then observed that this counting bound would itself imply an open inequality of the type π(x+y)≤π(x)+π(y)+O(y/(log⁡y)2)\pi(x+y)\le\pi(x)+\pi(y)+O(y/(\log y)^2)π(x+y)≤π(x)+π(y)+O(y/(logy)2) (compare Problem 855). The problem remains open.

Setting

Fix a natural number nnn. Consider finite sets AAA of integers with 1≤a<n1\le a<n1≤a<n for every a∈Aa\in Aa∈A, and with gcd⁡(a,b)=1\gcd(a,b)=1gcd(a,b)=1 for all distinct a,b∈Aa,b\in Aa,b∈A. For such a set define the reciprocal gap sum

Sn(A)=∑a∈A1n−a.S_n(A)=\sum_{a\in A}\frac{1}{n-a}.Sn​(A)=a∈A∑​n−a1​.

Every term is at most 111, and the element a=n−da=n-da=n−d contributes 1/d1/d1/d. Write ∑p<n1/p\sum_{p<n}1/p∑p<n​1/p for the sum of reciprocals of the primes below nnn; by Mertens' theorem it equals log⁡log⁡n+O(1)\log\log n+O(1)loglogn+O(1). Throughout, π(x)\pi(x)π(x) denotes the number of primes p≤xp\le xp≤x.

Target

The goal of the mission is the affirmative answer to Problem 1210: there is an absolute constant CCC such that

Sn(A)  ≤  ∑p<n1p+CS_n(A)\;\le\;\sum_{p<n}\frac1p+CSn​(A)≤p<n∑​p1​+C

for every nnn and every pairwise coprime A⊆[1,n)A\subseteq[1,n)A⊆[1,n). A negative answer is equally welcome and is recorded by disproving the goal statement.

The milestones are:

  1. Small prime factors. For pairwise coprime AAA, at most π(x)\pi(x)π(x) elements of A∩[n−x,n)A\cap[n-x,n)A∩[n−x,n) have a prime factor ≤x\le x≤x.
  2. Partial summation reduction. If ∣A∩[n−x,n)∣≤π(x)+O(x/(log⁡x)2)|A\cap[n-x,n)|\le\pi(x)+O(x/(\log x)^2)∣A∩[n−x,n)∣≤π(x)+O(x/(logx)2) uniformly, then the goal holds.
  3. The [Er77c] variant. For the primes qiq_iqi​ in (n,m](n,m](n,m], ∑i1/(qi−n)<∑p<m−n1/p+O(1)\sum_i 1/(q_i-n)<\sum_{p<m-n}1/p+O(1)∑i​1/(qi​−n)<∑p<m−n​1/p+O(1).

Significance

The result itself. An affirmative answer would say that, for the weight 1/(n−a)1/(n-a)1/(n−a), no pairwise coprime set beats the primes by more than a constant, and would give a quantitative form of the heuristic that coprime sets behave like sets of primes near a point. A negative answer would exhibit coprime sets that concentrate near nnn more efficiently than the primes do in the harmonic sense. The [Er77c] variant concerns only primes, and relates the distribution of primes just above nnn to the primes below the interval length m−nm-nm−n.

Formalizing it. Neither the goal nor the [Er77c] variant is known. The mission produces Lean statements checked against the source, a reduction (milestone 2), to be verified in Lean, that isolates exactly which counting estimate would suffice, and the elementary coprimality lemma (milestone 1). These pin down what a proof or disproof must supply.

Difficulty

The natural first attempt splits A∩[n−x,n)A\cap[n-x,n)A∩[n−x,n) into elements with a prime factor ≤x\le x≤x, of which there are at most π(x)\pi(x)π(x), and xxx-rough elements, and then hopes that sieve bounds make the rough part O(x/(log⁡x)2)O(x/(\log x)^2)O(x/(logx)2). The obstruction is that AAA may contain many primes in [n−x,n)[n-x,n)[n−x,n). Bounding the number of primes in a short interval [n−x,n)[n-x,n)[n−x,n) by π(x)+O(x/(log⁡x)2)\pi(x)+O(x/(\log x)^2)π(x)+O(x/(logx)2) is a form of the second Hardy–Littlewood conjecture π(x+y)≤π(x)+π(y)\pi(x+y)\le\pi(x)+\pi(y)π(x+y)≤π(x)+π(y), which is open and known to be incompatible, in its exact form, with the prime kkk-tuples conjecture. So the counting route in milestone 2 needs input on primes in short intervals beyond current knowledge, and any proof of the goal must either supply such input or avoid pointwise counting.

Formalization scope

All objects are elementary: AAA is a Finset ℕ, coprimality is Nat.Coprime, primes are Nat.Prime, π\piπ is Nat.primeCounting, and the sums are real-valued. The source's O(1)O(1)O(1) is encoded as an existentially quantified real constant CCC chosen before nnn and AAA. The source asks a yes/no question; each statement is posed in its affirmative form, and a disproof (a proof of the negation) settles the negative answer. The standing hypothesis 1≤a<n1\le a<n1≤a<n means every denominator n−an-an−a is at least 111, so no division-by-zero default can make the statement trivial. The window [n−x,n)[n-x,n)[n−x,n) is written as a≥n−xa\ge n-xa≥n−x with truncated natural subtraction together with a<na<na<n.

No new definitions are required. Useful reusable contributions include Mertens-type estimates for ∑p<n1/p\sum_{p<n}1/p∑p<n​1/p, partial summation lemmas for finite sums over N\mathbb NN, and upper bounds for primes in short intervals.

Selected references

  • P. Erdős, Problems and results on combinatorial number theory. III, Number Theory Day (Proc. Conf., Rockefeller Univ., New York, 1976), (1977), 43–72. [Er77c]
  • P. Erdős, A survey of problems in combinatorial number theory, Ann. Discrete Math. (1980), 89–115. [Er80]
  • T. F. Bloom, Erdős Problem #1210, https://www.erdosproblems.com/1210 , and discussion thread https://www.erdosproblems.com/forum/thread/1210
  • Formal Conjectures project, ErdosProblems/1210.lean, https://github.com/google-deepmind/formal-conjectures
4 thms3 active usersReviewed
Combinatorics·Captain: Lucas

Erdős Problem 52: the Erdős–Szemerédi sum–product conjectureOpen Problem

Motivation

Addition and multiplication interact in rigid ways: a finite set of numbers that is highly structured with respect to one operation (an arithmetic progression, say) tends to be unstructured with respect to the other (a geometric progression). The sum–product problem asks for the sharp quantitative form of this principle. It was posed by Erdős and Szemerédi in 1983 (Erdős Problem 52) and has since become a central question of additive combinatorics, with applications in incidence geometry, exponential sum estimates, expanders and randomness extraction.

Timeline.

  • 1983 — Erdős and Szemerédi show that max⁡(∣A+A∣,∣AA∣)≥c∣A∣1+δ\max(|A+A|,|AA|)\ge c|A|^{1+\delta}max(∣A+A∣,∣AA∣)≥c∣A∣1+δ for some absolute δ>0\delta>0δ>0 and every finite set of integers AAA, and conjecture exponent 2−ε2-\varepsilon2−ε.
  • 1997 — Nathanson obtains the explicit exponent 1+1311+\tfrac1{31}1+311​; Ford (1998) improves it to 1+1151+\tfrac1{15}1+151​.
  • 1997 — Elekes, using the Szemerédi–Trotter incidence theorem, proves ∣A+A∣ ∣AA∣≫∣A∣5/2|A+A|\,|AA|\gg|A|^{5/2}∣A+A∣∣AA∣≫∣A∣5/2 for finite sets of reals, hence exponent 5/45/45/4.
  • 2009 — Solymosi proves ∣A+A∣2∣AA∣≫∣A∣4/log⁡∣A∣|A+A|^2|AA|\gg |A|^4/\log|A|∣A+A∣2∣AA∣≫∣A∣4/log∣A∣ for finite sets of positive reals, hence exponent 4/34/34/3 up to a logarithmic factor.
  • 2015–2022 — Konyagin and Shkredov first break the 4/34/34/3 barrier (exponent 4/3+c4/3+c4/3+c for a small explicit c>0c>0c>0); after several improvements, Rudnev and Stevens reach 4/3+2/11674/3+2/11674/3+2/1167 up to logarithmic factors.

The conjecture itself remains open.

Setting

For a finite set A⊂ZA\subset\mathbb ZA⊂Z define the sumset and product set

A+A={a+b:a,b∈A},AA={ab:a,b∈A}.A+A=\{a+b : a,b\in A\},\qquad AA=\{ab : a,b\in A\}.A+A={a+b:a,b∈A},AA={ab:a,b∈A}.

For nonempty AAA both contain at least ∣A∣|A|∣A∣ elements and at most (∣A∣+12)\binom{|A|+1}{2}(2∣A∣+1​). The quantity of interest is max⁡(∣A+A∣,∣AA∣)\max(|A+A|,|AA|)max(∣A+A∣,∣AA∣) as a function of ∣A∣|A|∣A∣.

Formalization targets

Goal (Erdős–Szemerédi conjecture)

For every 0<ε<10<\varepsilon<10<ε<1 there is Cε>0C_\varepsilon>0Cε​>0 such that for every finite set A⊂ZA\subset\mathbb ZA⊂Z,

max⁡(∣A+A∣,∣AA∣) ≥ Cε ∣A∣2−ε.\max(|A+A|,|AA|)\ \ge\ C_\varepsilon\,|A|^{2-\varepsilon}.max(∣A+A∣,∣AA∣) ≥ Cε​∣A∣2−ε.

Known lower bounds (milestones, weakest to strongest)

∣A+A∣≥2∣A∣−1,max⁡(∣A+A∣,∣AA∣)≫∣A∣1+δ,≫∣A∣5/4,≫∣A∣4/3(log⁡∣A∣)1/3,≫ε∣A∣4/3+2/1167−ε.|A+A|\ge 2|A|-1,\qquad \max(|A+A|,|AA|)\gg|A|^{1+\delta},\qquad \gg|A|^{5/4},\qquad \gg\frac{|A|^{4/3}}{(\log|A|)^{1/3}},\qquad \gg_\varepsilon |A|^{4/3+2/1167-\varepsilon}.∣A+A∣≥2∣A∣−1,max(∣A+A∣,∣AA∣)≫∣A∣1+δ,≫∣A∣5/4,≫(log∣A∣)1/3∣A∣4/3​,≫ε​∣A∣4/3+2/1167−ε.

Sharpness

The ε\varepsilonε cannot be removed: there is no C>0C>0C>0 with max⁡(∣A+A∣,∣AA∣)≥C∣A∣2\max(|A+A|,|AA|)\ge C|A|^2max(∣A+A∣,∣AA∣)≥C∣A∣2 for all AAA (take A={1,…,n}A=\{1,\dots,n\}A={1,…,n}; the multiplication table ∣AA∣|AA|∣AA∣ is o(n2)o(n^2)o(n2) by Erdős).

Significance

A proof of the goal would settle the sharp form of the sum–product phenomenon over Z\mathbb ZZ. Sum–product estimates are an input to incidence bounds, to Bourgain–Katz–Tao-type results over finite fields, and to explicit constructions in theoretical computer science; improvements of the exponent over R\mathbb RR have come together with new incidence-geometric tools.

Formalization status: the milestones are published theorems (Erdős–Szemerédi, Elekes, Solymosi, Rudnev–Stevens, Erdős's multiplication table bound), but their proofs are not known to be formalized in Lean/Mathlib. The Szemerédi–Trotter theorem, multiplicative energy, and the Elekes and Solymosi arguments are reusable infrastructure. The goal itself is an open problem.

Difficulty

Incidence-geometric methods (Szemerédi–Trotter and its descendants) naturally produce exponents near 4/34/34/3, and passing beyond 4/34/34/3 has required intricate higher-energy arguments yielding only small gains. None of the existing approaches is known to reach exponents close to 222; even over Z\mathbb ZZ, where arithmetic structure is available, the best general bounds are the real-number ones.

Formalization scope

Sets are Finset ℤ; A+AA+AA+A and AAAAAA are Mathlib's pointwise sumset and product set (open scoped Pointwise), and cardinalities are cast to R\mathbb RR. Powers are real powers (Real.rpow). The empty set is allowed; the hypothesis ε<1\varepsilon<1ε<1 keeps every exponent positive, so the empty set contributes the trivial inequality 0≥00\ge 00≥0 rather than a junk value 00=10^0=100=1. The constant CCC may depend on ε\varepsilonε but not on AAA. The Solymosi milestone is stated for ∣A∣≥2|A|\ge2∣A∣≥2 so that log⁡∣A∣>0\log|A|>0log∣A∣>0.

The statements are for integer sets only; results proved over R\mathbb RR specialise to them. Contributions of general-purpose infrastructure (Szemerédi–Trotter over R\mathbb RR, multiplicative energy, bounds for the multiplication table) are welcome.

Selected references

  • P. Erdős, E. Szemerédi, On sums and products of integers, Studies in Pure Mathematics, Birkhäuser, 1983, 213–218.
  • M. B. Nathanson, On sums and products of integers, Proc. Amer. Math. Soc. 125 (1997), 9–16.
  • K. Ford, Sums and products from a finite set of real numbers, Ramanujan J. 2 (1998), 59–66.
  • G. Elekes, On the number of sums and products, Acta Arith. 81 (1997), 365–367.
  • J. Solymosi, Bounding multiplicative energy by the sumset, Adv. Math. 222 (2009), 402–408. https://arxiv.org/abs/0806.1040
  • S. V. Konyagin, I. D. Shkredov, On sum sets of sets having small product set, Proc. Steklov Inst. Math. 290 (2015), 288–299. https://arxiv.org/abs/1503.05771
  • M. Rudnev, S. Stevens, An update on the sum-product problem, Math. Proc. Cambridge Philos. Soc. 173 (2022), 411–430. https://arxiv.org/abs/2005.11145
  • Erdős Problem 52, https://www.erdosproblems.com/52
7 thms3 active usersReviewed
Captain: Lucas

Gilbreath's ConjectureOpen Problem

Motivation

Write the primes in increasing order, take the absolute differences of consecutive entries, take the absolute differences of the resulting row, and repeat. Every row produced this way appears to begin with 111:

2357111317…122424…10222…1200…120…\begin{array}{llllllll} 2 & 3 & 5 & 7 & 11 & 13 & 17 & \dots\\ 1 & 2 & 2 & 4 & 2 & 4 & \dots\\ 1 & 0 & 2 & 2 & 2 & \dots\\ 1 & 2 & 0 & 0 & \dots\\ 1 & 2 & 0 & \dots \end{array}21111​32022​52200​7420…​1122…​134…​17……

Gilbreath's conjecture asserts that this never fails. The observation is due to Norman L. Gilbreath (1958), who rediscovered a statement already published by François Proth in 1878 together with an argument that is not accepted as a proof. It is attractive because it is elementary to state and because it is one of the few statements about the primes whose difficulty is not visibly analytic: it concerns the combinatorics of iterated differences rather than the distribution of primes directly.

Timeline.

  • 1878 — Proth states the property and publishes a proof that is now regarded as erroneous.
  • 1958 — Gilbreath rediscovers the pattern; it circulates as a conjecture.
  • 1959 — Killgrove and Ralston verify the leading entry for the first 63,41863{,}41863,418 rows (MTAC 13 (1959), 121–122).
  • 1993 — Odlyzko reports a verification of the leading entry for all rows of index at most π(1013)≈3.4×1011\pi(10^{13}) \approx 3.4 \times 10^{11}π(1013)≈3.4×1011, using an argument that propagates a long block of entries lying in {0,2}\{0,2\}{0,2} downwards through the triangle (Math. Comp. 61 (1993), 373–380).

No proof is known.

Setting

Let p0=2<p1=3<p2=5<…p_0 = 2 < p_1 = 3 < p_2 = 5 < \dotsp0​=2<p1​=3<p2​=5<… be the increasing enumeration of the prime numbers, indexed from 000. Define the rows of the Gilbreath triangle by

d0(n)=pn,dk+1(n)=∣dk(n+1)−dk(n)∣(k,n≥0).d^0(n) = p_n, \qquad d^{k+1}(n) = \bigl| d^{k}(n+1) - d^{k}(n) \bigr| \quad (k, n \ge 0).d0(n)=pn​,dk+1(n)=​dk(n+1)−dk(n)​(k,n≥0).

Thus dkd^kdk is an infinite sequence of natural numbers for every kkk, row 000 is the sequence of primes, row 111 is the sequence of prime gaps pn+1−pnp_{n+1}-p_npn+1​−pn​, and each later row is the sequence of absolute differences of consecutive entries of the row above it. Only the leading entry dk(0)d^k(0)dk(0) of each row is at issue.

More generally, for an arbitrary sequence a:N→Na : \mathbb{N} \to \mathbb{N}a:N→N write (Δa)(n)=∣a(n+1)−a(n)∣(\Delta a)(n) = |a(n+1) - a(n)|(Δa)(n)=∣a(n+1)−a(n)∣ and Δja\Delta^j aΔja for the jjj-fold iterate, so that dk=Δkpd^k = \Delta^k pdk=Δkp.

Formalization targets

Goal

∀k≥1,dk(0)=1.\forall k \ge 1,\qquad d^{k}(0) = 1 .∀k≥1,dk(0)=1.

This is the conjecture in its standard form: every row after the row of primes begins with 111. It fixes no constants and no ranges, so no computational advance can invalidate it.

Milestones

The milestone list collects the statements that a proof, or a further computational verification, would be built from: the two low-level structural facts about the triangle (row 111 is the gap sequence; from row 111 on, the leading entry is odd and all later entries are even), a finite verification of the first rows, and the two statements underlying Odlyzko's method — the propagation lemma for an arbitrary sequence beginning 111 and continuing in {0,2}\{0,2\}{0,2}, and the reduction of the conjecture to the existence, for each row index, of an earlier row with a long enough block of entries in {0,2}\{0,2\}{0,2}.

Significance

The result itself. The conjecture is not known to imply other open statements about the primes, and its interest lies elsewhere: it is a test case for how much of the fine structure of the prime sequence is forced by coarse information. The propagation mechanism shows that the conjecture for a given row index follows from purely local data about an earlier row, and that mechanism is what every verification to date has relied on. A proof would have to show that such blocks of entries in {0,2}\{0,2\}{0,2} always appear early enough, which is a statement about the density of small prime gaps in disguise.

Formalizing it. Nothing here is currently formalized: Mathlib has the prime enumeration n↦pnn \mapsto p_nn↦pn​ (Nat.nth Nat.Prime) and the basic facts about it, but not the iterated-difference triangle nor any of its properties. This mission contributes the definition of the triangle, the structural facts about its rows, and a machine-checked version of the reduction step that all computational work on the problem uses. The goal theorem itself is open — the milestones are known mathematics, and each is provable with current tools, while the goal is not.

Difficulty

The obvious attack is induction on the row index: to see that dk+1(0)=1d^{k+1}(0) = 1dk+1(0)=1 it suffices to know that dk(0)=1d^{k}(0) = 1dk(0)=1 and dk(1)∈{0,2}d^{k}(1) \in \{0,2\}dk(1)∈{0,2}. But controlling dk(1)d^{k}(1)dk(1) requires controlling dk−1(1)d^{k-1}(1)dk−1(1) and dk−1(2)d^{k-1}(2)dk−1(2), and so on: the invariant that closes is not "the row begins with 111" but "the row begins with 111 and its next mmm entries lie in {0,2}\{0,2\}{0,2}", and each application of the difference operator consumes one entry of that block. So a finite block of good entries only carries the conclusion a finite number of rows further down, and the conjecture needs such blocks to keep reappearing forever, arbitrarily far down the triangle. Nothing is known that produces them.

A second warning, due to Hallard Croft: the property is not specific to the primes. Sequences that start with 222, continue with odd numbers, and have gaps that are not too large empirically exhibit the same behaviour, so any proof that uses only such coarse features would prove a much more general statement — and conversely, an argument exploiting deep properties of primes is likely to be proving the wrong thing.

Formalization scope

Rows are total functions N→N\mathbb{N} \to \mathbb{N}N→N, defined for every index, and the whole triangle is a single family indexed by the row number. Differences are taken as Int.natAbs of a difference computed in Z\mathbb{Z}Z, so truncated natural subtraction never occurs; the one place where N\mathbb{N}N-subtraction does appear is the milestone identifying row 111 with the gap sequence, where the subtraction is justified by monotonicity of n↦pnn \mapsto p_nn↦pn​.

Primes are indexed from 000 via Mathlib's Nat.nth Nat.Prime, so p0=2p_0 = 2p0​=2; rows are indexed with row 000 the primes, and the goal quantifies over all k≥1k \ge 1k≥1 in the form d (k + 1) 0 = 1, with no upper bound and no extra hypothesis, so no vacuous or finitely-truncated reading of the goal is available. The general difference operator is stated for arbitrary sequences N→N\mathbb{N} \to \mathbb{N}N→N, which is what makes the propagation lemma usable as a black box, and reusable beyond this mission.

A complete development needs no analytic input for the milestones: Mathlib's Nat.nth, Nat.prime_nth_prime, Nat.nth_prime_zero_eq_two and the strict monotonicity of the prime enumeration suffice. Contributions that would extend the mission beyond its current list: a formal version of a concrete computational verification (checking that the leading entries of the first NNN rows are 111 for an NNN well beyond the hand-checkable range), and formalizations of the general statement for non-prime sequences of the Croft type.

Selected references

  • N. L. Gilbreath, as reported in R. B. Killgrove and K. E. Ralston, On a conjecture concerning the primes, Mathematical Tables and Other Aids to Computation 13 (1959), 121–122. https://doi.org/10.1090/S0025-5718-1959-0105398-3
  • A. M. Odlyzko, Iterated absolute values of differences of consecutive primes, Mathematics of Computation 61 (1993), 373–380. https://doi.org/10.1090/S0025-5718-1993-1192979-9
  • Gilbreath's conjecture, Wikipedia. https://en.wikipedia.org/wiki/Gilbreath%27s_conjecture
23 thms3 active usersReviewed
Harmonic Analysis·Captain: Lucas

Connes: Weil positivity and the Riemann zeta functionResearch Paper

Motivation

The Riemann hypothesis (RH) asserts that every zero of the Riemann zeta function ζ\zetaζ in the strip 0<Re⁡s<10 < \operatorname{Re} s < 10<Res<1 has Re⁡s=12\operatorname{Re} s = \tfrac12Res=21​. One of the few reformulations that turns RH into a positivity statement, rather than a statement about the location of points, goes back to A. Weil (1952): the explicit formula expresses a sum over the zeros of ζ\zetaζ as a sum of local contributions over the places of Q\mathbb{Q}Q, and RH is equivalent to the resulting functional being positive on elements of the form g⋆g∗g \star g^{*}g⋆g∗.

Connes' 1999 programme paper Noncommutative geometry and the Riemann zeta function takes this reformulation as its endpoint. It builds a geometric framework — the adele class space X=A/k∗X = \mathbb{A}/k^{*}X=A/k∗ carrying an action of the idele class group CkC_kCk​ — in which the explicit formula appears as a Lefschetz formula, the zeros of LLL-functions appear spectrally, and the paper's concluding assertion (§3, p. 22) is that the validity of the global trace formula implies, and is in fact equivalent to, positivity of the Weil distribution, i.e. RH for all LLL-functions with Grössencharakter.

This mission formalizes the arithmetic core of that endpoint in its simplest instance: the global field k=Qk = \mathbb{Q}k=Q with trivial Grössencharakter, so that the LLL-function is ζ\zetaζ itself. Concretely it asks for (i) the Riemann–Weil explicit formula for ζ\zetaζ in the shape of Connes' equation (11), and (ii) both directions of the equivalence between positivity of the resulting Weil distribution and RH.

A rough timeline of the statements involved: Riemann (1859) gave the first explicit formula; von Mangoldt (1895) proved it rigorously; Weil (Sur les "formules explicites" de la théorie des nombres premiers, 1952) extended it to all global fields and isolated the positivity criterion; Bombieri (Remarks on Weil's quadratic functional in the theory of prime numbers, 2000) studied the associated quadratic functional in detail; Connes (1996–1999) gave the trace-formula interpretation formalized in part here.

Setting

All objects live on the group R+∗\mathbb{R}^{*}_{+}R+∗​, the module of the idele class group of Q\mathbb{Q}Q, written additively through u=etu = e^{t}u=et, d∗u=dtd^{*}u = dtd∗u=dt.

A test function is a map g:R→Cg : \mathbb{R} \to \mathbb{C}g:R→C that is C∞C^{\infty}C∞ and has compact support (IsTest).

Its transform is

g^(z)  =  ∫Rg(t) e(z−1/2) t dt,\widehat{g}(z) \;=\; \int_{\mathbb{R}} g(t)\, e^{(z - 1/2)\,t}\, dt ,g​(z)=∫R​g(t)e(z−1/2)tdt,

which is Connes' h^(z)=∫Ckh(u) ∣u∣z d∗u\widehat h(z) = \int_{C_k} h(u)\,|u|^{z}\,d^{*}uh(z)=∫Ck​​h(u)∣u∣zd∗u in the coordinate u=etu = e^{t}u=et, shifted by 12\tfrac1221​ so that zzz is the variable of ζ\zetaζ (mellinHat). On the critical line, g^(12+ir)=∫Rg(t)eirt dt\widehat{g}(\tfrac12 + i r) = \int_{\mathbb{R}} g(t) e^{irt}\,dtg​(21​+ir)=∫R​g(t)eirtdt is the ordinary Fourier transform.

The involution is g∗(t)=g(−t)‾g^{*}(t) = \overline{g(-t)}g∗(t)=g(−t)​, i.e. h∗(u)=h(u−1)‾h^{*}(u) = \overline{h(u^{-1})}h∗(u)=h(u−1)​ (starInv), and convolution is (g1⋆g2)(t)=∫Rg1(s) g2(t−s) ds(g_1 \star g_2)(t) = \int_{\mathbb{R}} g_1(s)\, g_2(t-s)\, ds(g1​⋆g2​)(t)=∫R​g1​(s)g2​(t−s)ds (conv).

The Weil distribution of a test function ggg collects the pole terms, the finite places and the archimedean place:

W(g)  =  g^(0)+g^(1)  −  ∑n≥2Λ(n)n(g(log⁡n)+g(−log⁡n))  +  12π∫Rg^(12+ir)(Re⁡ψ(14+ir2)−log⁡π)dr,W(g) \;=\; \widehat{g}(0) + \widehat{g}(1) \;-\; \sum_{n \ge 2} \frac{\Lambda(n)}{\sqrt{n}}\bigl(g(\log n) + g(-\log n)\bigr) \;+\; \frac{1}{2\pi}\int_{\mathbb{R}} \widehat{g}\left(\tfrac12 + ir\right)\Bigl(\operatorname{Re}\psi\left(\tfrac14 + \tfrac{ir}{2}\right) - \log \pi\Bigr) dr ,W(g)=g​(0)+g​(1)−n≥2∑​n​Λ(n)​(g(logn)+g(−logn))+2π1​∫R​g​(21​+ir)(Reψ(41​+2ir​)−logπ)dr,

where Λ\LambdaΛ is the von Mangoldt function and ψ=Γ′/Γ\psi = \Gamma'/\Gammaψ=Γ′/Γ (weilDistribution, with the three pieces named mellinHat, primeSum, archTerm). The middle sum is Weil's contribution of the finite places v=pv = pv=p, the last integral the contribution of the real place.

The spectral side is

Z(g)  =  ∑ρmρ g^(ρ),Z(g) \;=\; \sum_{\rho} m_{\rho}\, \widehat{g}(\rho),Z(g)=ρ∑​mρ​g​(ρ),

the sum over the zeros ρ\rhoρ of ζ\zetaζ with 0<Re⁡ρ<10 < \operatorname{Re}\rho < 10<Reρ<1, each counted with its multiplicity mρm_\rhomρ​ (zeroSum, IsCriticalZero, zeroMult).

Formalization targets

Goal — positivity of the Weil distribution implies RH

(∀ g test:  Re⁡W(g⋆g∗)≥0)  ⟹  (∀ρ, ζ(ρ)=0, 0<Re⁡ρ<1  ⇒  Re⁡ρ=12).\Bigl(\forall\, g \text{ test}:\; \operatorname{Re} W\bigl(g \star g^{*}\bigr) \ge 0\Bigr) \;\Longrightarrow\; \bigl(\forall \rho,\ \zeta(\rho) = 0,\ 0 < \operatorname{Re}\rho < 1 \;\Rightarrow\; \operatorname{Re}\rho = \tfrac12\bigr).(∀g test:ReW(g⋆g∗)≥0)⟹(∀ρ, ζ(ρ)=0, 0<Reρ<1⇒Reρ=21​).

This is the direction that yields RH, and it is the weakest form of the endpoint of the paper: it fixes no rate, no test-function normalization beyond Cc∞C_c^{\infty}Cc∞​, and no numerical constant.

Milestone — the explicit formula (Connes (11))

∑ρmρ g^(ρ)  =  W(g)for every test function g,\sum_{\rho} m_\rho\,\widehat g(\rho) \;=\; W(g) \qquad \text{for every test function } g,ρ∑​mρ​g​(ρ)=W(g)for every test function g,

with the sum over zeros asserted to be (unconditionally) summable.

Milestone — the converse direction

RH  ⟹  ∀ g test: Re⁡W(g⋆g∗)≥0.\text{RH} \;\Longrightarrow\; \forall\, g \text{ test}:\ \operatorname{Re} W\bigl(g \star g^{*}\bigr) \ge 0 .RH⟹∀g test: ReW(g⋆g∗)≥0.

Together with the goal this is the equivalence asserted on p. 22 of the paper, in the case k=Qk = \mathbb{Q}k=Q, trivial Grössencharakter.

Supporting statements

The ∗*∗-identity g⋆g∗^(12+ir)=∣g^(12+ir)∣2\widehat{g \star g^{*}}\left(\tfrac12 + ir\right) = \bigl|\widehat g\left(\tfrac12+ir\right)\bigr|^{2}g⋆g∗​(21​+ir)=​g​(21​+ir)​2 on the critical line, and the fact that g⋆g∗g \star g^{*}g⋆g∗ is again a test function.

Significance

Weil's positivity criterion is one of the standard equivalent forms of RH, and the only one in which the arithmetic input (the primes, through Λ\LambdaΛ) and the archimedean input (the Γ\GammaΓ-factor) enter as separate, explicitly computable local terms. Formalizing it produces a machine-checked bridge between the zeros of ζ\zetaζ and prime sums: the explicit formula milestone is the reusable object here, since essentially every analytic application of zeta zeros — zero-density estimates, prime-counting error terms, pair-correlation statistics — is an instance of it.

Status honesty: neither RH nor the positivity statement is known; the explicit formula and both implications relating positivity to RH are classical theorems, proved but not, as far as the catalog shows, formalized in Lean. Mathlib currently provides ζ\zetaζ, its functional equation, the von Mangoldt function and Γ\GammaΓ, but no explicit formula of any kind. What this mission adds on top of the paper is therefore the formal proof of known results, not new mathematics.

Difficulty

The obvious route to the explicit formula — integrate −ζ′/ζ(s)g^(s)-\zeta'/\zeta(s)\widehat g(s)−ζ′/ζ(s)g​(s) over a vertical line, move the contour to the reflected line, collect residues — fails to be routine at exactly two points. First, moving the contour requires control of ζ′/ζ\zeta'/\zetaζ′/ζ on horizontal segments between zeros, which is where the classical proof invests most of its work; Mathlib has bounds near Re⁡s=1\operatorname{Re} s = 1Res=1 but nothing of this shape inside the strip. Second, the sum over zeros must be shown to converge unconditionally, which needs a zero-counting bound of Riemann–von Mangoldt type (N(T)≪Tlog⁡TN(T) \ll T\log TN(T)≪TlogT) that is not in Mathlib either.

For the goal implication, the naive idea — pick a test function whose transform is supported near a hypothetical off-line zero — is unavailable: g^\widehat gg​ is entire whenever ggg has compact support, so it cannot be localized. The classical argument instead exploits the symmetry ρ↦1−ρˉ\rho \mapsto 1-\bar\rhoρ↦1−ρˉ​ of the zero set and makes the off-line quadruple contribute a negative amount in the limit along a family of test functions.

Formalization scope

Conventions the Lean statements commit to. Test functions are C\mathbb{C}C-valued on R\mathbb{R}R, ContDiff ℝ (⊤ : ℕ∞) (so C∞C^{\infty}C∞, not analytic) with HasCompactSupport; the multiplicative group R+∗\mathbb{R}^{*}_{+}R+∗​ is always written additively. The transform carries the 12\tfrac1221​-shift shown above, so the critical line is Re⁡z=12\operatorname{Re} z = \tfrac12Rez=21​ and g^(0),g^(1)\widehat{g}(0), \widehat{g}(1)g​(0),g​(1) are the two pole terms. Zeros are indexed by the subtype {s:0<Re⁡s<1, ζ(s)=0}\{s : 0 < \operatorname{Re} s < 1,\ \zeta(s) = 0\}{s:0<Res<1, ζ(s)=0} and weighted by (analyticOrderAt riemannZeta s).toNat; the trivial zeros are excluded. Integrals are Bochner integrals and sums are tsum, so both take the junk value 000 when the integrand is not integrable or the family is not summable — for that reason the explicit formula is stated as a HasSum, which carries summability, rather than as an equation between tsums. The archimedean term is written with Re⁡ψ\operatorname{Re}\psiReψ, ψ=\psi = ψ= logDeriv Complex.Gamma, rather than as a principal value, to avoid a second regularization convention.

The goal is not trivially satisfiable: its hypothesis quantifies over a nonempty class (smooth bump functions exist), and its conclusion is RH for ζ\zetaζ, so no vacuous reading is available.

Infrastructure a complete development needs, all reusable beyond this mission: growth bounds for ζ′/ζ\zeta'/\zetaζ′/ζ inside the critical strip, a Riemann–von Mangoldt zero-counting bound, Fourier analysis of Cc∞C^\infty_cCc∞​ functions (Paley–Wiener style decay of g^\widehat gg​), and the Hadamard product / functional equation package for the completed zeta function.

Out of scope, and deliberately so: Connes' operator-theoretic trace formula (equations (41) and (45) of the paper) and the spectral realization theorem of p. 17. Both are statements about traces of operators on Hilbert space, and Mathlib presently has no trace-class operator theory to state them faithfully. The mission therefore formalizes the arithmetic side of the paper's endpoint; contributions that build the missing operator theory, or that extend the statements from ζ\zetaζ to Dirichlet LLL-functions and Hecke LLL-functions with Grössencharakter, are welcome.

Selected references

  • A. Connes, Noncommutative geometry and the Riemann zeta function, in Mathematics: Frontiers and Perspectives, AMS (2000) — the source of this mission (§3, equations (11) and (45), and the concluding assertion on p. 22).
  • A. Connes, Trace formula in noncommutative geometry and the zeros of the Riemann zeta function, Selecta Math. (N.S.) 5 (1999) — reference [9] of the source. https://arxiv.org/abs/math/9811068
  • A. Weil, Sur les "formules explicites" de la théorie des nombres premiers, Comm. Sém. Math. Univ. Lund (1952) — reference [27] of the source.
  • E. Bombieri, Remarks on Weil's quadratic functional in the theory of prime numbers, I (2000).
  • H. Iwaniec and E. Kowalski, Analytic Number Theory, AMS Colloquium Publications 53 (2004), Chapter 5 (explicit formulas).
7 thms3 active usersReviewed
🏆Completed
Analysis·Captain: Lucas

The de Bruijn–Newman Constant is Non-negativeResearch Paper

Motivation

The Riemann hypothesis asserts that all nontrivial zeros of the Riemann zeta function lie on the critical line. A classical way to measure how far the hypothesis is from failing runs through a one-parameter deformation of the Riemann ξ\xiξ function by the backward heat flow. De Bruijn (1950) introduced a family of entire functions HtH_tHt​, t∈Rt \in \mathbb{R}t∈R, with H0H_0H0​ essentially the ξ\xiξ function, and showed that HtH_tHt​ has only real zeros for t≥1/2t \ge 1/2t≥1/2. Newman (1976) proved that there is a finite constant Λ\LambdaΛ, now called the de Bruijn–Newman constant, such that HtH_tHt​ has only real zeros precisely when t≥Λt \ge \Lambdat≥Λ. The Riemann hypothesis is exactly the statement Λ≤0\Lambda \le 0Λ≤0, and Newman conjectured the complementary bound Λ≥0\Lambda \ge 0Λ≥0 — in his phrase, that if the Riemann hypothesis is true, then it is only barely so.

Timeline of lower bounds on Λ\LambdaΛ, all obtained before 2018 by exhibiting Lehmer pairs, that is, pairs of adjacent zeros of ζ\zetaζ that are unusually close together: Λ>−∞\Lambda > -\inftyΛ>−∞ (Newman 1976), Λ≥−50\Lambda \ge -50Λ≥−50 (Csordas–Norfolk–Varga 1988), Λ≥−5\Lambda \ge -5Λ≥−5 (te Riele 1991), Λ≥−0.385\Lambda \ge -0.385Λ≥−0.385 (Norfolk–Ruttan–Varga 1992), Λ≥−0.0991\Lambda \ge -0.0991Λ≥−0.0991 (Csordas–Ruttan–Varga 1991), Λ≥−4.379×10−6\Lambda \ge -4.379 \times 10^{-6}Λ≥−4.379×10−6 (Csordas–Smith–Varga 1994), Λ≥−5.895×10−9\Lambda \ge -5.895 \times 10^{-9}Λ≥−5.895×10−9 (Csordas–Odlyzko–Smith–Varga 1993), Λ≥−2.63×10−9\Lambda \ge -2.63 \times 10^{-9}Λ≥−2.63×10−9 (Odlyzko 2000), Λ≥−1.15×10−11\Lambda \ge -1.15 \times 10^{-11}Λ≥−1.15×10−11 (Saouter–Gourdon–Demichel 2011). Rodgers and Tao closed the gap in 2020 by proving Λ≥0\Lambda \ge 0Λ≥0. In the other direction, de Bruijn's bound Λ≤1/2\Lambda \le 1/2Λ≤1/2 was sharpened to Λ<1/2\Lambda < 1/2Λ<1/2 by Ki–Kim–Lee (2009) and to Λ≤0.22\Lambda \le 0.22Λ≤0.22 by the Polymath 15 project (2019).

Setting

For a real number uuu put

Φ(u):=∑n=1∞(2π2n4e9u−3πn2e5u)exp⁡(−πn2e4u),\Phi(u) := \sum_{n=1}^{\infty}\bigl(2\pi^2 n^4 e^{9u} - 3\pi n^2 e^{5u}\bigr)\exp\bigl(-\pi n^2 e^{4u}\bigr),Φ(u):=n=1∑∞​(2π2n4e9u−3πn2e5u)exp(−πn2e4u),

a function that decays super-exponentially as ∣u∣→∞|u| \to \infty∣u∣→∞ and satisfies Φ(u)=Φ(−u)\Phi(u) = \Phi(-u)Φ(u)=Φ(−u). For each t∈Rt \in \mathbb{R}t∈R define the entire function

Ht(z):=∫0∞etu2 Φ(u) cos⁡(zu) du.H_t(z) := \int_0^{\infty} e^{t u^2}\,\Phi(u)\,\cos(z u)\,du .Ht​(z):=∫0∞​etu2Φ(u)cos(zu)du.

Each HtH_tHt​ is even and satisfies Ht(zˉ)=Ht(z)‾H_t(\bar z) = \overline{H_t(z)}Ht​(zˉ)=Ht​(z)​; the function H0H_0H0​ is 18ξ(12+iz2)\tfrac18 \xi\bigl(\tfrac12 + \tfrac{iz}{2}\bigr)81​ξ(21​+2iz​), so the Riemann hypothesis says exactly that every zero of H0H_0H0​ is real. Write

S:={ t∈R:every zero of Ht is real },Λ:=inf⁡S.S := \{\, t \in \mathbb{R} : \text{every zero of } H_t \text{ is real} \,\},\qquad \Lambda := \inf S .S:={t∈R:every zero of Ht​ is real},Λ:=infS.

By Pólya and Newman, SSS is the ray [Λ,∞)[\Lambda, \infty)[Λ,∞) with −∞<Λ≤1/2-\infty < \Lambda \le 1/2−∞<Λ≤1/2.

When Λ<t≤0\Lambda < t \le 0Λ<t≤0 the zeros of HtH_tHt​ are real, simple, symmetric about the origin and avoid the origin, so they can be listed as (xj(t))j∈Z∗(x_j(t))_{j \in \mathbb{Z}^*}(xj​(t))j∈Z∗​, indexed by the nonzero integers, with 0<x1(t)<x2(t)<⋯0 < x_1(t) < x_2(t) < \cdots0<x1​(t)<x2​(t)<⋯ and x−j(t)=−xj(t)x_{-j}(t) = -x_j(t)x−j​(t)=−xj​(t). The classical locations ξj\xi_jξj​ are defined for j≥1j \ge 1j≥1 by Ψ(ξj)=j\Psi(\xi_j) = jΨ(ξj​)=j with

Ψ(T):=T4πlog⁡T4π−T4π,\Psi(T) := \frac{T}{4\pi}\log\frac{T}{4\pi} - \frac{T}{4\pi},Ψ(T):=4πT​log4πT​−4πT​,

extended by ξ−j=−ξj\xi_{-j} = -\xi_jξ−j​=−ξj​; they are the positions the zeros would occupy if the Riemann–von Mangoldt counting formula were exact. Throughout, log⁡+x:=log⁡(2+∣x∣)\log_+ x := \log(2 + |x|)log+​x:=log(2+∣x∣).

Formalization targets

Goal — Newman's conjecture

Λ≥0,equivalentlyevery t with Ht having only real zeros satisfies t≥0.\Lambda \ge 0, \qquad\text{equivalently}\qquad \text{every } t \text{ with } H_t \text{ having only real zeros satisfies } t \ge 0 .Λ≥0,equivalentlyevery t with Ht​ having only real zeros satisfies t≥0.

The goal is stated in both forms simultaneously, so that it does not depend on any convention for the infimum of a set that might be empty or unbounded below.

Milestones

The milestone list follows the architecture of Rodgers–Tao, which is a proof by contradiction: every milestone is stated under the standing hypothesis Λ<0\Lambda < 0Λ<0 of that paper, in the time ranges the paper uses (Λ<t≤0\Lambda < t \le 0Λ<t≤0, then Λ/2≤t≤0\Lambda/2 \le t \le 0Λ/2≤t≤0, then Λ/4≤t≤0\Lambda/4 \le t \le 0Λ/4≤t≤0). In order: an upper bound for HtH_tHt​ near the real axis (Lemma 4); Riemann–von Mangoldt type counting formulae for the zeros of HtH_tHt​ (Theorem 9); the resulting macroscopic description of the zeros (Corollary 10); the equations of motion ∂txk=2∑j≠k(xk−xj)−1\partial_t x_k = 2\sum_{j \ne k} (x_k - x_j)^{-1}∂t​xk​=2∑j=k​(xk​−xj​)−1 (Theorem 11); a quantitative lower bound on gaps between zeros (Proposition 13); a bound on the time-integrated renormalized energy (Theorem 17); and a bound on that energy at time t=0t = 0t=0 (Proposition 26). The last of these says that at time zero the zeros are, on average, locally in the equilibrium configuration of an arithmetic progression, which contradicts known results on the local distribution of zeros of ζ\zetaζ.

Significance

Λ≥0\Lambda \ge 0Λ≥0 settles Newman's conjecture, and together with the Riemann hypothesis it would force Λ=0\Lambda = 0Λ=0. Unconditionally, it says that the zeros of ξ\xiξ are not in local equilibrium: infinitely often, gaps between consecutive zeros deviate from the mean spacing, which is what makes the pair correlation phenomenology of Montgomery and of Conrey–Ghosh–Goldston–Gonek–Heath-Brown incompatible with Λ<0\Lambda < 0Λ<0. Any proof of the Riemann hypothesis must therefore be compatible with the hypothesis being tight in this sense.

The theorem has a complete published proof (Rodgers–Tao, Forum of Mathematics, Pi, 2020); it is not an open problem. What is missing is a machine-checked proof. To the extent the material has been formalized at all, the underlying objects — the ξ\xiξ function, the heat flow HtH_tHt​, the counting function for zeros, the zero dynamics, the renormalized energies — are not available in Mathlib, so the mission produces reusable analytic infrastructure: bounds for a Fourier–Laplace type integral by the saddle point method, a Riemann–von Mangoldt counting argument via the argument principle, and a gradient-flow monotonicity framework for an infinite particle system with logarithmic interaction.

Difficulty

The obvious route to Λ≥0\Lambda \ge 0Λ≥0 is the one used for every previous lower bound: exhibit Lehmer pairs of ever higher quality, since if Λ\LambdaΛ were very negative the zeros of H0H_0H0​ would repel each other and unusually close pairs of zeta zeros could not exist. Producing an infinite sequence of Lehmer pairs of arbitrarily high quality is possible under the GUE hypothesis, but the known unconditional upper bounds for small gaps between zeta zeros are too weak, even assuming the Riemann hypothesis. The proof instead upgrades repulsion to relaxation to local equilibrium: it must control the zeros of HtH_tHt​ uniformly for Λ<t≤0\Lambda < t \le 0Λ<t≤0 at length scales as fine as log⁡T\log TlogT, with only the weaker counting formulae available for negative ttt (an error term O(log⁡+2T)O(\log_+^2 T)O(log+2​T) rather than O(log⁡+T)O(\log_+ T)O(log+​T)), and must make sense of a Hamiltonian and an energy that are given by divergent series, which requires truncation, renormalization, and careful control of all the resulting boundary terms.

Formalization scope

The Lean development commits to the following conventions. Φ\PhiΦ is a tsum over the positive integers and Ht(z)H_t(z)Ht​(z) is the Bochner integral over (0,∞)(0, \infty)(0,∞) of etu2Φ(u)cos⁡(zu)e^{tu^2}\Phi(u)\cos(zu)etu2Φ(u)cos(zu); no convergence or entireness statement is built into the definition. Λ\LambdaΛ is sInf of the set of admissible times, and the goal theorem also states the quantifier form "every admissible ttt is nonnegative", so it cannot be satisfied by a junk value of the infimum. The zero families (xj(t))(x_j(t))(xj​(t)) and the classical locations (ξj)(\xi_j)(ξj​) are not defined by choice functions: they enter the milestones as universally quantified functions Z→R\mathbb{Z} \to \mathbb{R}Z→R subject to explicit predicates saying exactly which sequences they are, so a milestone asserts something about every valid enumeration. Asymptotic notation is unfolded: O(⋅)O(\cdot)O(⋅) becomes an explicit existential constant, oT→∞(⋅)o_{T \to \infty}(\cdot)oT→∞​(⋅) an explicit ε\varepsilonε–T0T_0T0​ statement, and a principal value sum a limit of symmetric partial sums. Where a statement asserts the value of a time integral, absolute integrability is part of the conclusion, so the statement cannot be satisfied by the convention that a non-integrable function has integral zero.

One degeneracy is inherent to the source and is stated here explicitly: since the paper argues by contradiction, each milestone carries the hypothesis Λ<0\Lambda < 0Λ<0 (directly, or through a time range such as Λ<t≤0\Lambda < t \le 0Λ<t≤0). Once the goal theorem is proved, those hypotheses are unsatisfiable and the milestones become vacuously true. They are the intended attack path on the goal, not independent targets, and a solver who derives one of them from the goal theorem contributes nothing.

Contributions welcome: the analytic estimates for HtH_tHt​ (Lemma 4) and the counting formulae (Theorem 9) are independent of the dynamical part and are the natural entry points; Mathlib-level infrastructure on the argument principle, the saddle point method, and Stirling asymptotics for Γ\GammaΓ in vertical strips is reusable well beyond this mission.

Selected references

  • B. Rodgers and T. Tao, The de Bruijn–Newman constant is non-negative, Forum of Mathematics, Pi 8 (2020), e6. https://doi.org/10.1017/fmp.2020.6
  • N. G. de Bruijn, The roots of trigonometric integrals, Duke Math. J. 17 (1950), 197–226. https://doi.org/10.1215/S0012-7094-50-01720-0
  • C. M. Newman, Fourier transforms with only real zeros, Proc. Amer. Math. Soc. 61 (1976), 246–251. https://doi.org/10.1090/S0002-9939-1976-0434982-5
  • G. Csordas, W. Smith and R. S. Varga, Lehmer pairs of zeros, the de Bruijn–Newman constant Λ\LambdaΛ, and the Riemann hypothesis, Constr. Approx. 10 (1994), 107–129. https://doi.org/10.1007/BF01205170
  • H. L. Montgomery, The pair correlation of zeros of the zeta function, Proc. Sympos. Pure Math. XXIV (1973), 181–193. https://doi.org/10.1090/pspum/024
  • J. B. Conrey, A. Ghosh, D. Goldston, S. M. Gonek and D. R. Heath-Brown, On the distribution of gaps between zeros of the zeta-function, Q. J. Math. 36 (1985), 43–51. https://doi.org/10.1093/qmath/36.1.43
  • D. H. J. Polymath, Effective approximation of heat flow evolution of the Riemann ξ\xiξ function, and a new upper bound for the de Bruijn–Newman constant, Res. Math. Sci. 6 (2019), 31. https://doi.org/10.1007/s40687-019-0193-1
25 thms3 active usersReviewed
AlgebraGroup Theory·Captain: Lucas

The Inverse Galois ProblemOpen Problem

Motivation

Galois theory attaches to every finite Galois extension L/KL/KL/K a finite group Gal(L/K)\mathrm{Gal}(L/K)Gal(L/K), the group of field automorphisms of LLL fixing KKK pointwise, and the fundamental theorem of Galois theory turns the subfield structure of L/KL/KL/K into the subgroup structure of that group. The inverse Galois problem asks whether this correspondence is surjective over the rationals: given an arbitrary finite group GGG, is there a Galois extension L/QL/\mathbb{Q}L/Q with Gal(L/Q)≅G\mathrm{Gal}(L/\mathbb{Q}) \cong GGal(L/Q)≅G? The question was posed in the early nineteenth century and is unsolved.

What makes it a live research question rather than a curiosity is that the known positive results come from genuinely different sources, and none of them covers all finite groups.

  • Cyclic and, more generally, finite abelian groups are realizable over Q\mathbb{Q}Q by an explicit cyclotomic construction resting on Dirichlet's theorem on primes in arithmetic progressions.
  • Symmetric and alternating groups are realizable over Q\mathbb{Q}Q; this is due to Hilbert, who realized them first over the rational function field Q(t)\mathbb{Q}(t)Q(t) and then specialized ttt using his irreducibility theorem.
  • Every finite solvable group is realizable over Q\mathbb{Q}Q; this is Shafarevich's theorem (I. R. Shafarevich, The imbedding problem for splitting extensions, Dokl. Akad. Nauk SSSR 120 (1958), 1217–1219), obtained by solving embedding problems.
  • Over C(t)\mathbb{C}(t)C(t) — and over K(t)K(t)K(t) for any algebraically closed KKK of characteristic zero — every finite group is realizable, by the Riemann existence theorem. The obstruction to the goal is not the group theory; it is descending the field of constants to Q\mathbb{Q}Q.
  • Case-by-case work covers large finite lists: all transitive permutation groups of degree at most 232323, and every sporadic simple group, are known to be realizable over Q\mathbb{Q}Q.

Setting

Fix a field KKK and a group GGG. A Galois realization of GGG over KKK is a field LLL equipped with a KKK-algebra structure such that the extension L/KL/KL/K is Galois — normal and separable — together with a group isomorphism

G  ≅  Gal(L/K),G \;\cong\; \mathrm{Gal}(L/K),G≅Gal(L/K),

where Gal(L/K)\mathrm{Gal}(L/K)Gal(L/K) denotes the group of KKK-algebra automorphisms of LLL under composition. The group GGG is realizable over KKK, written IsRealizable K G, when at least one Galois realization of GGG over KKK exists. No finiteness of L/KL/KL/K is imposed in the definition; it is automatic once GGG is finite, because an infinite Galois extension has infinite automorphism group.

Two base fields beyond Q\mathbb{Q}Q appear throughout. K(t)K(t)K(t) denotes the field of rational functions in one variable over KKK, written RatFunc K; and for the statement that a group is realizable over some number field, the base field ranges over the intermediate fields of C/Q\mathbb{C}/\mathbb{Q}C/Q.

Formalization targets

Goal — the inverse Galois problem

for every finite group G,∃ L/Q Galois with Gal(L/Q)≅G.\text{for every finite group } G, \qquad \exists\, L/\mathbb{Q} \text{ Galois with } \mathrm{Gal}(L/\mathbb{Q}) \cong G.for every finite group G,∃L/Q Galois with Gal(L/Q)≅G.

The goal fixes no degree, no polynomial and no construction: it asserts only the shape of the truth, so no later refinement of the known constructions can invalidate it.

Milestones — the known partial results

G cyclic  ⟹  G realizable over Q,G abelian  ⟹  G realizable over Q,G \text{ cyclic} \;\Longrightarrow\; G \text{ realizable over } \mathbb{Q}, \qquad G \text{ abelian} \;\Longrightarrow\; G \text{ realizable over } \mathbb{Q},G cyclic⟹G realizable over Q,G abelian⟹G realizable over Q, Sym(S),  An realizable over Q,G solvable  ⟹  G realizable over Q,\mathrm{Sym}(S),\; A_n \text{ realizable over } \mathbb{Q}, \qquad G \text{ solvable} \;\Longrightarrow\; G \text{ realizable over } \mathbb{Q},Sym(S),An​ realizable over Q,G solvable⟹G realizable over Q, ∃ K, Q⊆K⊆C, G realizable over K,\exists\, K,\ \mathbb{Q} \subseteq K \subseteq \mathbb{C},\ G \text{ realizable over } K,∃K, Q⊆K⊆C, G realizable over K, G realizable over C(t),G realizable over K(t) (K algebraically closed, char 0),G \text{ realizable over } \mathbb{C}(t), \qquad G \text{ realizable over } K(t) \ (K \text{ algebraically closed, char } 0),G realizable over C(t),G realizable over K(t) (K algebraically closed, char 0), G realizable over Q(t)  ⟹  G realizable over Q.G \text{ realizable over } \mathbb{Q}(t) \;\Longrightarrow\; G \text{ realizable over } \mathbb{Q}.G realizable over Q(t)⟹G realizable over Q.

The last milestone is the Hilbert-irreducibility descent step; together with the geometric milestones it makes precise which half of the classical programme is missing.

Significance

The result itself would settle a two-century-old question and, with it, the surjectivity of the Galois correspondence over Q\mathbb{Q}Q: every abstract finite group would be known to arise from an explicit arithmetic object, a polynomial with rational coefficients. Its absence is felt in practice — constructing a single new Galois group over Q\mathbb{Q}Q is publishable work, as the recent additions of the degree-171717 group 17T717T717T7 (van Bommel–Costa–Elkies–Keller–Schiavone–Voight, 2024) and of the Mathieu group M23M_{23}M23​ show.

Formalizing it produces something available today independently of the goal: a machine-checked library of the known realizability results. Mathlib has the fundamental theorem of Galois theory, cyclotomic extensions, the Kronecker–Weber theorem, solvability of groups and symmetric/alternating group theory, but it does not have a predicate for "GGG is a Galois group over KKK", nor any of the milestones above. Every milestone here is a proved theorem of classical number theory and an unformalized one; the cyclic and abelian cases are within reach of current Mathlib, while the Shafarevich and Riemann-existence milestones are substantial formalization projects in their own right.

Difficulty

The obvious strategy fails at a well-understood point. Over C(t)\mathbb{C}(t)C(t) the problem is solved: by the Riemann existence theorem every finite group occurs as the deck-transformation group of a branched cover of the projective line. Hilbert's irreducibility theorem then descends realizability from Q(t)\mathbb{Q}(t)Q(t) to Q\mathbb{Q}Q. What is missing is the step in between: producing the cover over Q\mathbb{Q}Q rather than over C\mathbb{C}C, i.e. showing that the geometric solution can be chosen with rational field of constants. The rigidity method makes this work for many groups, but there is no known argument covering all of them; an approach that only produces realizability over some number field is not enough, and that weaker statement is included as a milestone precisely to mark the line.

A second, purely formal difficulty: the milestones are classical but their published proofs are long. Shafarevich's theorem rests on a delicate analysis of embedding problems, and the Riemann existence theorem is analytic input that Mathlib does not currently have in the required form.

Formalization scope

The mission fixes one definition file, published first, carrying the structure GaloisRealization and the one-field class IsRealizable. Conventions it commits to:

  • IsGalois K L is Mathlib's Galois condition (normal and separable); finiteness of the extension is not assumed.
  • The isomorphism is with the full automorphism group L≃alg[K]LL \simeq_{\mathrm{alg}[K]} LL≃alg[K]​L, not with a quotient or a subgroup of it.
  • The carrier LLL of a realization is required to live in the same universe as KKK. This costs no generality for the statements of the mission — for finite GGG a realization is a finite extension of KKK — and keeps every statement universe-monomorphic.
  • Sym(S)\mathrm{Sym}(S)Sym(S) is Equiv.Perm S for a finite type SSS, and AnA_nAn​ is alternatingGroup (Fin n); degenerate small cases are included rather than excluded.
  • Solvability is Group.IsSolvable.

The statements cannot be satisfied vacuously: IsRealizable K G asserts the existence of data, so a solver must exhibit an extension; and the hypotheses of the milestones (cyclic, abelian, solvable, or none at all) are all satisfiable, so no milestone is empty. The one conditional milestone, Hilbert descent, is stated with realizability over Q(t)\mathbb{Q}(t)Q(t) as an explicit hypothesis.

Infrastructure a complete development needs, most of it reusable well beyond this mission: transport of a Galois realization along an isomorphism of groups and along an isomorphism of base fields; the fixed-field construction and the fundamental theorem in the form "Gal(L/LH)≅H\mathrm{Gal}(L/L^H) \cong HGal(L/LH)≅H"; Galois groups of cyclotomic fields; Dirichlet's theorem on primes in arithmetic progressions (already in Mathlib); Hilbert's irreducibility theorem (not in Mathlib). Contributions of any of these as reusable platform definitions or lemmas are welcome, as are decompositions of the harder milestones into sketches.

Selected references

  • Inverse Galois problem, Wikipedia. https://en.wikipedia.org/wiki/Inverse_Galois_problem
  • I. R. Shafarevich, The imbedding problem for splitting extensions, Dokl. Akad. Nauk SSSR 120 (1958), 1217–1219.
  • C. U. Jensen, A. Ledet, N. Yui, Generic Polynomials: Constructive Aspects of the Inverse Galois Problem, MSRI Publications 45, Cambridge University Press, 2002. http://library.msri.org/books/Book45/files/book45.pdf
  • G. Malle, B. H. Matzat, Inverse Galois Theory, Springer Monographs in Mathematics, 1999.
  • R. van Bommel, E. Costa, N. D. Elkies, T. Keller, S. Schiavone, J. Voight, 17T7 is a Galois group over the rationals, arXiv:2411.07857, 2024. https://arxiv.org/abs/2411.07857
22 thms3 active usersReviewed
Combinatorics·Captain: Zexuan Liu

Erdős Problem 142: Asymptotics for Sets Free of k-Term Arithmetic ProgressionsOpen Problem

Motivation

Erdős asked, repeatedly and with a rising price tag, for an asymptotic formula for the largest subset of {1,…,N}\{1,\dots,N\}{1,…,N} that contains no arithmetic progression of a given length. He offered 1000 dollars for it in [Er97c] and 10000 dollars in [Er81, p.4], where he called the question "probably enormously difficult"; elsewhere he described it as "probably unattackable at present". Most of modern additive combinatorics — the density increment method, the triangle removal lemma, Gowers uniformity norms, the arithmetic regularity lemma — grew out of attempts on this single question, and the answer is still unknown, even in the first non-trivial case k=3k=3k=3.

Timeline.

  • 1936: Erdős and Turán conjecture that rk(N)=o(N)r_k(N)=o(N)rk​(N)=o(N) for every kkk.
  • 1946: Behrend constructs large progression-free sets, giving r3(N)≥Nexp⁡(−clog⁡N)r_3(N)\ge N\exp(-c\sqrt{\log N})r3​(N)≥Nexp(−clogN​).
  • 1953: Roth proves r3(N)=o(N)r_3(N)=o(N)r3​(N)=o(N), with the quantitative form r3(N)≪N/log⁡log⁡Nr_3(N)\ll N/\log\log Nr3​(N)≪N/loglogN.
  • 1961: Rankin generalises Behrend, giving rk(N)≥Nexp⁡(−ck(log⁡N)1/(k−1))r_k(N)\ge N\exp(-c_k(\log N)^{1/(k-1)})rk​(N)≥Nexp(−ck​(logN)1/(k−1)).
  • 1969, 1975: Szemerédi proves r4(N)=o(N)r_4(N)=o(N)r4​(N)=o(N) and then rk(N)=o(N)r_k(N)=o(N)rk​(N)=o(N) for all kkk, settling Erdős–Turán.
  • 1977: Furstenberg reproves Szemerédi's theorem ergodically, with no effective bound.
  • 1998, 2001: Gowers introduces uniformity norms and obtains rk(N)≪N(log⁡log⁡N)−ckr_k(N)\ll N(\log\log N)^{-c_k}rk​(N)≪N(loglogN)−ck​, the first effective bound for general kkk.
  • 2017: Green and Tao obtain r4(N)≪N(log⁡N)−cr_4(N)\ll N(\log N)^{-c}r4​(N)≪N(logN)−c.
  • 2020: Bloom and Sisask obtain r3(N)≪N(log⁡N)−1−cr_3(N)\ll N(\log N)^{-1-c}r3​(N)≪N(logN)−1−c, the first bound past the N/log⁡NN/\log NN/logN barrier.
  • 2023: Kelley and Meka obtain r3(N)≤Nexp⁡(−c(log⁡N)1/12)r_3(N)\le N\exp(-c(\log N)^{1/12})r3​(N)≤Nexp(−c(logN)1/12).
  • 2024: Leng, Sah and Sawhney obtain rk(N)≪Nexp⁡(−(log⁡log⁡N)ck)r_k(N)\ll N\exp(-(\log\log N)^{c_k})rk​(N)≪Nexp(−(loglogN)ck​) for k≥5k\ge5k≥5.

Every upper bound in this list is still astronomically far from Behrend's lower bound, and no candidate asymptotic formula has been proposed for any k≥3k\ge3k≥3.

Setting

Fix an integer kkk. A non-trivial kkk-term arithmetic progression is a list a, a+d, a+2d,…,a+(k−1)da,\,a+d,\,a+2d,\dots,a+(k-1)da,a+d,a+2d,…,a+(k−1)d of natural numbers with common difference d>0d>0d>0; the requirement d>0d>0d>0 is what "non-trivial" means, and it forces the kkk terms to be distinct. A finite set A⊆NA\subseteq\mathbb NA⊆N is kkk-AP-free if it contains no such progression. Write

rk(N)  =  max⁡{ ∣A∣  :  A⊆{1,…,N}, A is k-AP-free }.r_k(N)\;=\;\max\bigl\{\,|A| \;:\; A\subseteq\{1,\dots,N\},\ A\ \text{is}\ k\text{-AP-free}\,\bigr\}.rk​(N)=max{∣A∣:A⊆{1,…,N}, A is k-AP-free}.

The mission takes its formal definition of rkr_krk​ verbatim from the formal-conjectures entry for this problem, so that the goal below is literally the statement recorded there. In that development a set is called free of progressions of length lll when every subset of it that is an arithmetic progression of length lll forces l≤1l\le1l≤1. Progressions of length 000 and 111 count as trivial, so under that convention every set is free of them and r0(N)=r1(N)=Nr_0(N)=r_1(N)=Nr0​(N)=r1​(N)=N; the interesting range begins at k≥2k\ge2k≥2. Every statement in this mission that depends on the convention carries an explicit hypothesis on kkk.

On that range, rk(N)r_k(N)rk​(N) is non-decreasing in both NNN and kkk, satisfies rk(M+N)≤rk(M)+rk(N)r_k(M+N)\le r_k(M)+r_k(N)rk​(M+N)≤rk​(M)+rk​(N), and hence, by Fekete's subadditivity lemma, rk(N)/Nr_k(N)/Nrk​(N)/N converges. Szemerédi's theorem is the statement that the limit is 000; the whole difficulty of this mission lies in how fast it goes to 000.

Formalization targets

Goal

rk(N)  =  ok ⁣(Nlog⁡N)for every k>1.r_k(N)\;=\;o_k\!\left(\frac{N}{\log N}\right)\qquad\text{for every }k>1 .rk​(N)=ok​(logNN​)for every k>1.

This is erdos_142.variants.lower of the formal-conjectures file for Erdős 142, reproduced binder for binder, over that file's own definition of rkr_krk​.

The headline theorem in that file, erdos_142, states rk(N)=Θ(f)r_k(N)=\Theta(f)rk​(N)=Θ(f) with the comparison function left as an answer(sorry) placeholder, and the same is true of its variants.upper and variants.three. Those are not closed propositions and cannot serve as a mission goal: the literal request of Erdős Problem #142 — "prove an asymptotic formula for rk(N)r_k(N)rk​(N)" — has no known right-hand side for any k≥3k\ge3k≥3, which is exactly why the file leaves a hole there. variants.lower is the one formalizable target in the file, and it is also the strongest precisely-stated form the problem page attaches to #142: Erdős offered 5000 dollars for (essentially) exactly it, as recorded under Erdős Problem #3. It is known for k=3k=3k=3 — it follows from Bloom–Sisask 2020, and a fortiori from Kelley–Meka 2023 — trivial for k=2k=2k=2, where r2(N)=1r_2(N)=1r2​(N)=1, and open for every k≥4k\ge4k≥4. It fixes no constants, so no future improvement can invalidate it.

A weaker open question

rk(n)rk+1(n)⟶0for some k≥3.\frac{r_k(n)}{r_{k+1}(n)}\longrightarrow 0\qquad\text{for some }k\ge3 .rk+1​(n)rk​(n)​⟶0for some k≥3.

Erdős remarked in [Er80, p.92] that even this separation between consecutive progression lengths is not known. Here [Er80] is the erdosproblems.com bibliography key for Erdős's 1980 paper; it is a citation, not a pointer to Erdős Problem #80, which is an unrelated question about books in graphs. This statement has no counterpart in formal-conjectures: the file for #142 contains only the Θ\ThetaΘ, ooo and OOO variants above, and the only two files in that repository that mention rkr_krk​ at all are the ones for #142 and #139.

Significance

Proving rk(N)=ok(N/log⁡N)r_k(N)=o_k(N/\log N)rk​(N)=ok​(N/logN) for all kkk yields, by a standard summation argument, Erdős's conjecture that every A⊆NA\subseteq\mathbb NA⊆N with ∑a∈A1/a=∞\sum_{a\in A}1/a=\infty∑a∈A​1/a=∞ contains arbitrarily long arithmetic progressions — the 5000-dollar Erdős Problem #3, of which the Green–Tao theorem on primes is the best-known special case. Below that threshold, quantitative bounds on rkr_krk​ control the density at which progressions must appear in any concrete set, and are the input to results on progressions in the primes, in sumsets, and in sparse random subsets of the integers.

Formalization status is uneven, and this mission is designed around that gap. Mathlib already contains the k=3k=3k=3 theory in a usable form: the predicate ThreeAPFree, the Roth number rothNumberNat, its subadditivity, and a complete formalization of Behrend's construction (Behrend.roth_lower_bound). Mathlib does not contain Roth's theorem, Szemerédi's theorem, or any of the modern upper bounds; to the best of current knowledge none of Roth, Szemerédi, Gowers, Green–Tao, Kelley–Meka or Leng–Sah–Sawhney has a machine-checked proof anywhere. The formal-conjectures entry states the problem but proves nothing: every declaration in it is a sorry. Two of this mission's targets are taken from that repository — the goal from its file for #142, and the Szemerédi milestone from its file for #139, which uses the same rkr_krk​; those are the only two files there that mention rkr_krk​. The milestones therefore split cleanly: the first six are reachable now on top of Mathlib, and the last five are open formalization projects of independent value.

Difficulty

Every known upper bound for rkr_krk​ runs a density increment: if A⊆{1,…,N}A\subseteq\{1,\dots,N\}A⊆{1,…,N} of density δ\deltaδ has no kkk-term progression, find a long subprogression on which AAA has density δ(1+c(δ))\delta(1+c(\delta))δ(1+c(δ)), and iterate. The bound this produces is governed entirely by two quantities — how large the increment c(δ)c(\delta)c(δ) is, and how much of the interval survives one step. For k≥4k\ge4k≥4 the increment is extracted from an inverse theorem for the Gowers Uk−1U^{k-1}Uk−1-norm, and the best available correlation bounds there are quasipolynomial in δ\deltaδ; iterating a quasipolynomial increment cannot do better than Nexp⁡(−(log⁡log⁡N)c)N\exp(-(\log\log N)^{c})Nexp(−(loglogN)c), which is nowhere near N/log⁡NN/\log NN/logN. Reaching N/log⁡NN/\log NN/logN requires an increment with polynomial dependence on δ\deltaδ together with a subprogression of polynomial length, and that combination is currently available only for k=3k=3k=3, through the sifting and almost-periodicity machinery of Kelley–Meka. No soft or averaging argument can substitute: Behrend's construction shows the truth at k=3k=3k=3 is Nexp⁡(−Θ(log⁡N))N\exp(-\Theta(\sqrt{\log N}))Nexp(−Θ(logN​)), so the answer is not a power of log⁡N\log NlogN and cannot be produced by any argument whose output has that shape.

Formalization scope

The mission's definition file Erdos142Basic carries two layers, and every statement in the mission is written against them.

  1. The source definitions, ported verbatim. IsAPOfLengthWith, IsAPOfLength, IsAPOfLengthFree and r are the declarations of the formal-conjectures entry, transcribed unchanged into the mission's namespace: a set is an arithmetic progression of length lll with first term aaa and difference ddd when it has exactly lll elements and equals {a+nd:n<l}\{a+nd : n<l\}{a+nd:n<l}; it is free of length-lll progressions when every progression of length lll inside it forces l≤1l\le1l≤1; and rk(N)r_k(N)rk​(N) is the supremum of ∣S∣|S|∣S∣ over subsets S⊆{1,…,N}S\subseteq\{1,\dots,N\}S⊆{1,…,N} free of length-kkk progressions. The ground set is Finset.Icc 1 N, and the supremum is sSup over N\mathbb NN; the file proves the two facts that make it a genuine maximum (le_r and r_le).

  2. An elementary handle. HasAP k A is ∃ a d, 0 < d ∧ ∀ i < k, a + i * d ∈ A, and APFree k A its negation. This form carries no cardinality side condition in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞} and is what a solver actually wants to induct on. The first milestone is exactly the bridge between the two layers.

Two consequences of the source convention are worth stating plainly, because the prose is silent about them. Length-000 and length-111 progressions are trivial, so every set is free of them and r0(N)=r1(N)=Nr_0(N)=r_1(N)=Nr0​(N)=r1​(N)=N; monotonicity of rkr_krk​ in kkk therefore holds only from k≥2k\ge2k≥2 onward, and the corresponding milestone carries that hypothesis. Asymptotic statements use Asymptotics.IsLittleO and Filter.atTop over N\mathbb NN with real-valued casts, and real division is Lean's, so (N : ℝ) / Real.log N is 000 at N=1N=1N=1; this is invisible to atTop.

A trivializing formalization is ruled out by construction: one milestone asserts r3(N)=r_3(N)=r3​(N)= rothNumberNat N, pinning this development against Mathlib's independently written definition of the Roth number, so a vacuous or mis-quantified notion of progression-freeness cannot survive. That milestone, the bridge milestone above it, and the monotonicity milestone have all been checked to be provable before this proposal was drafted.

A full development needs: discrete Fourier analysis on Z/NZ\mathbb Z/N\mathbb ZZ/NZ, Bohr sets and their regularity, the arithmetic regularity lemma, Gowers uniformity norms and the inverse theorem for them, and — for the lower bounds — sphere-counting in high-dimensional boxes (already in Mathlib via Behrend). All of this is reusable well beyond this mission. Contributions of any kind are welcome, including partial results: quantitative bounds weaker than the cited ones, the k=3k=3k=3 case of a general-kkk milestone, and reusable Fourier-analytic infrastructure are all valuable even when they do not close a milestone.

Selected references

  • Erdős Problem #142. https://www.erdosproblems.com/142
  • Erdős Problem #3. https://www.erdosproblems.com/3
  • Erdős Problem #139 (Szemerédi's theorem in the rkr_krk​ formulation), linked from #142. https://www.erdosproblems.com/139
  • Google DeepMind, formal-conjectures, FormalConjectures/ErdosProblems/142.lean — the source of the goal statement and of the definition of rkr_krk​. https://github.com/google-deepmind/formal-conjectures/blob/main/FormalConjectures/ErdosProblems/142.lean
  • F. A. Behrend, On sets of integers which contain no three terms in arithmetical progression, Proc. Nat. Acad. Sci. USA 32 (1946), 331–332. https://doi.org/10.1073/pnas.32.12.331
  • K. F. Roth, On certain sets of integers, J. London Math. Soc. 28 (1953), 104–109. https://doi.org/10.1112/jlms/s1-28.1.104
  • R. A. Rankin, Sets of integers containing not more than a given number of terms in arithmetical progression, Proc. Roy. Soc. Edinburgh Sect. A 65 (1961), 332–344.
  • E. Szemerédi, On sets of integers containing no kkk elements in arithmetic progression, Acta Arith. 27 (1975), 199–245. https://doi.org/10.4064/aa-27-1-199-245
  • W. T. Gowers, A new proof of Szemerédi's theorem, Geom. Funct. Anal. 11 (2001), 465–588. https://doi.org/10.1007/s00039-001-0332-9
  • B. Green and T. Tao, New bounds for Szemerédi's theorem, III: A polylogarithmic bound for r4(N)r_4(N)r4​(N), Mathematika 63 (2017), 944–1040. https://arxiv.org/abs/1705.01703
  • T. F. Bloom and O. Sisask, Breaking the logarithmic barrier in Roth's theorem on arithmetic progressions, arXiv:2007.03528. https://arxiv.org/abs/2007.03528
  • Z. Kelley and R. Meka, Strong bounds for 3-progressions, arXiv:2302.05537. https://arxiv.org/abs/2302.05537
  • J. Leng, A. Sah and M. Sawhney, Improved bounds for Szemerédi's theorem, arXiv:2402.17995. https://arxiv.org/abs/2402.17995
37 thms3 active usersReviewed
Captain: xuanji

Lagarias criterion is equivalent to RHResearch Paper

Motivation

The Riemann hypothesis concerns the zeros of a complex analytic function, yet it admits an equivalent formulation involving only positive integers, finite sums, the real exponential, and the natural logarithm. Jeffrey C. Lagarias established this formulation in An Elementary Problem Equivalent to the Riemann Hypothesis (Theorem 1.1). It connects the distribution of divisors of an integer with the analytic behavior of the zeta function. For number theorists and formalizers, the interest lies in making that connection precise without confusing an elementary statement with an elementary proof.

The objective is the known equivalence selected by LeanEval v1, not a resolution of RH. Lagarias's paper builds on results of Guy Robin concerning large values of the divisor-sum function; those results remain substantial parts of the formalization workload (Lagarias, §3).

Setting

For a positive integer nnn, its divisor sum is

σ(n)=∑d∣nd,\sigma(n)=\sum_{d\mid n}d,σ(n)=d∣n∑​d,

where the sum runs over positive divisors, including 111 and nnn. Its harmonic number is

Hn=∑j=1n1j.H_n=\sum_{j=1}^{n}\frac1j.Hn​=j=1∑n​j1​.

All inequalities below are inequalities of real numbers. The symbols exp⁡\expexp and log⁡\loglog denote the real exponential and natural logarithm. The Euler–Mascheroni constant is γ=lim⁡n→∞(Hn−log⁡n)\gamma=\lim_{n\to\infty}(H_n-\log n)γ=limn→∞​(Hn​−logn).

The Riemann zeta function is obtained by analytic continuation of ∑m=1∞m−s\sum_{m=1}^{\infty}m^{-s}∑m=1∞​m−s from Re⁡(s)>1\operatorname{Re}(s)>1Re(s)>1. RH asserts that its nontrivial zeros have real part 1/21/21/2. The Lagarias elementary criterion in this mission is the assertion that σ(n)≤Hn+exp⁡(Hn)log⁡(Hn)\sigma(n)\le H_n+\exp(H_n)\log(H_n)σ(n)≤Hn​+exp(Hn​)log(Hn​) for every positive integer nnn. These conventions agree with the arithmetic quantities in Lagarias, Problem E, with the precise equality-clause distinction stated below.

Formalization targets

Main goal: the exact LeanEval equivalence

RH⟺∀n∈N,  n>0⟹σ(n)≤Hn+exp⁡(Hn)log⁡(Hn).\mathrm{RH}\quad\Longleftrightarrow\quad \forall n\in\mathbb N,\;n>0\Longrightarrow \sigma(n)\le H_n+\exp(H_n)\log(H_n).RH⟺∀n∈N,n>0⟹σ(n)≤Hn​+exp(Hn​)log(Hn​).

There are no hypotheses on the goal theorem. The quantifier ranges over all positive integers, not a bounded test set or an unspecified tail. The benchmark uses a non-strict inequality and does not include an equality characterization. Lagarias's Problem E additionally requires equality only at n=1n=1n=1; the proof of the reverse implication in Theorem 1.1, p. 8 uses the non-strict inequality alone. The stronger source formulation is therefore not silently substituted for the benchmark.

Supporting targets from the paper

The milestone list records the following source statements, with their thresholds unchanged:

  • Lemma 3.1: for n≥3n\ge3n≥3,
eγnlog⁡log⁡n≤exp⁡(Hn)log⁡(Hn).e^\gamma n\log\log n\le\exp(H_n)\log(H_n).eγnloglogn≤exp(Hn​)log(Hn​).
  • Lemma 3.2: for n≥20n\ge20n≥20,
Hn+exp⁡(Hn)log⁡(Hn)≤eγnlog⁡log⁡n+7nlog⁡n.H_n+\exp(H_n)\log(H_n)\le e^\gamma n\log\log n+\frac{7n}{\log n}.Hn​+exp(Hn​)log(Hn​)≤eγnloglogn+logn7n​.
  • The finite check in the proof of Theorem 1.1: the criterion holds for 1≤n≤50401\le n\le50401≤n≤5040, with equality exactly at n=1n=1n=1.
  • Proposition 3.1, attributed to Robin: assuming RH, for n≥5041n\ge5041n≥5041,
σ(n)≤eγnlog⁡log⁡n.\sigma(n)\le e^\gamma n\log\log n.σ(n)≤eγnloglogn.
  • Proposition 3.2, attributed to Robin: if RH is false, some fixed 0<β<1/20<\beta<1/20<β<1/2 and C>0C>0C>0 satisfy
σ(n)≥eγnlog⁡log⁡n+Cnlog⁡log⁡n(log⁡n)β\sigma(n)\ge e^\gamma n\log\log n+ \frac{Cn\log\log n}{(\log n)^\beta}σ(n)≥eγnloglogn+(logn)βCnloglogn​

for arbitrarily large integers nnn.

All five are taken from Lagarias, §3, pp. 6–8; the finite check is explicitly an unnumbered step, not a newly attributed lemma.

Significance

The result identifies an exact arithmetic reformulation of RH. It does not make either side unconditional. A proof of the equivalence gives a bridge between statements in different mathematical languages; it does not certify the universal inequality merely because many instances can be checked. This distinction is central to the interpretation of Lagarias's theorem.

The formalization would connect existing Mathlib definitions of the zeta function, divisor sums, harmonic numbers, and Euler's constant through a machine-checked argument. Reusable outputs include explicit harmonic/exponential comparisons, certified finite real inequalities, and formal versions of Robin's conditional and oscillation results. The mathematical results are known; this proposal supplies open formalization targets, not completed proofs. No accepted LeanEval result is claimed by creating or launching the mission.

Difficulty

The elementary appearance of the criterion hides its main analytic requirements. Bounding the divisor sum crudely, or checking any finite number of integers, cannot establish the universal equivalence. The conditional upper bound and especially the quantitative oscillation theorem connect zeta zeros with unusually large divisor sums. They must be proved, not packaged as definitions or presumed available because the paper cites them (Lagarias, Propositions 3.1–3.2).

The oscillation statement requires uniform positive constants and arbitrarily large indices. Replacing it with one counterexample loses essential information. The bounded computation also requires rigorous control of exponential and logarithmic values: an ordinary floating-point loop is not a Lean proof. Beyond the listed milestones, completion still requires standard growth comparisons, threshold bookkeeping, and assembly of the two implications. The short length of the source's final argument should not be read as an estimate of total formalization effort.

Formalization scope

The goal is LeanEval.NumberTheory.riemann_hypothesis_iff_lagarias_elementary_criterion, with type RiemannHypothesis ↔ LagariasElementaryCriterion. The criterion definition is copied from the benchmark. σ 1 n is natural-valued and cast to the reals; harmonic n is rational-valued and cast to the reals. RH remains Mathlib's predicate on riemannZeta, excluding negative even trivial zeros and the point s=1s=1s=1. No replacement axiom, hidden RH assumption, altered zeta function, or circular child restatement is permitted.

The goal excludes n=0n=0n=0 and includes n=1n=1n=1. Thresholds ensure positive logarithm arguments in the analytic milestones. The oscillation milestone expresses an infinite subset of the naturals as an unbounded set, retaining n≥3n\ge3n≥3; deletion of the finitely many smaller indices does not change the source's infinitude claim. Its real exponent is represented by Real.rpow, not natural exponentiation. Constants are chosen before the arbitrary cutoff.

The benchmark pins Lean 4.33.0 and Mathlib 6f1ef4e5dd604a435bddba4747b13970cd65d2a1. The proposal targets Prove2me's supported Lean 4.33.1 environment, Mathlib 0df444a360eaa60ab8c11dca51a86af692955474. These environments are distinct; eventual benchmark credit requires the benchmark's own validation. Contributions to analytic infrastructure, source-faithful supporting results, and certified finite inequalities are welcome. The five milestones are an initial source-backed structure, not a claim that all required infrastructure is already present.

Selected references

  • Jeffrey C. Lagarias, An Elementary Problem Equivalent to the Riemann Hypothesis, American Mathematical Monthly 109 (2002), 534–543. arXiv:math/0008177v2, posted 6 May 2001. The theorem, proposition, equation, and page numbers in this proposal refer to this nine-page arXiv version.
  • Guy Robin, Grandes valeurs de la fonction somme des diviseurs et hypothèse de Riemann, Journal de Mathématiques Pures et Appliquées 63 (1984), 187–213. Bibliography entry [18] in Lagarias. The milestone formulations are those explicitly reproduced and attributed in Lagarias's Propositions 3.1 and 3.2.
7 thms3 active usersReviewed
Differential GeometryGeometry & Topology·Captain: t4v1

Thurston's Question 23: rational relations among hyperbolic volumesOpen Problem

Motivation

In the last of the twenty-four questions that closed his 1982 survey Three-dimensional manifolds, Kleinian groups and hyperbolic geometry (Bull. Amer. Math. Soc. 6 (1982), 357–381), Thurston asked to "show that volumes of hyperbolic 333-manifolds are not all rationally related" (p. 380). Twenty-two of the twenty-four have since been answered — geometrization by Perelman, tameness by Agol and by Calegari–Gabai, the ending lamination conjecture by Brock–Canary–Minsky, virtual fibering by Agol — and this one is among the two that remain open.

Some rational relations are forced, and for a trivial reason: a degree nnn cover of a hyperbolic 333-manifold has nnn times its volume, so any two commensurable manifolds have rationally related volumes. The question, which remains open, is whether every rational relation arises that way — equivalently, whether some two hyperbolic 333-manifolds have irrational volume ratio. Remarkably, not a single such pair is known.

Setting

The bundle fixes the meaning of every term. Hyperbolic 333-space is the upper half-space {(x,y,z):z>0}\{(x,y,z) : z > 0\}{(x,y,z):z>0}. Its volume is Lebesgue measure with density z−3z^{-3}z−3 — the Riemannian volume of the metric (dx2+dy2+dz2)/z2(dx^2+dy^2+dz^2)/z^2(dx2+dy2+dz2)/z2 written out, so that no Riemannian machinery is required. The hyperbolic distance is given by its closed formula

cosh⁡d(p,q)  =  1+∣p−q∣22 p3 q3.\cosh d(p,q) \;=\; 1 + \frac{|p-q|^2}{2\,p_3\,q_3}.coshd(p,q)=1+2p3​q3​∣p−q∣2​.

A Kleinian action is a free, properly discontinuous action by hyperbolic isometries; the quotient is a complete hyperbolic 333-manifold, discreteness and torsion freeness being consequences rather than hypotheses. The volume of the quotient is the measure of a fundamental domain, in the sense of Mathlib's MeasureTheory.IsFundamentalDomain, and the set of volumes collects those that are finite and positive.

Two conventions are stated rather than derived, and are worth flagging. Isometries are not required to preserve orientation, so the set of volumes also contains those of non-orientable quotients; this enlarges the set but not its Q\mathbb{Q}Q-span, so neither goal is affected. And preservation of the hyperbolic volume is a field of the structure rather than a consequence of preserving the distance: it holds for every hyperbolic isometry, but deriving it amounts to classifying Isom(H3)\mathrm{Isom}(\mathbb{H}^3)Isom(H3), which is not the subject of this mission.

Formalization targets

The goal is that the volumes are not all rationally related: there are two of them, vvv and www, with v≠qwv \neq q wv=qw for every rational qqq.

Two milestones support it. The first is that passing to a subgroup of index nnn multiplies the volume by nnn, a fundamental domain for the subgroup being the union of nnn translates of one for the whole group; this is the source of every known rational relation, and it is why the question is phrased as it is. The second is that the set of volumes is nonempty — that some finite-volume hyperbolic 333-manifold exists at all — without which the goal would be vacuously false rather than open.

A stronger form of the question, that the Q\mathbb{Q}Q-span of the set of volumes is infinite dimensional, is also stated.

Significance

The question is a geometric statement whose difficulty is arithmetic. For the Bianchi groups of an imaginary quadratic field FFF, Humbert's formula gives the covolume as ∣δF∣3/2ζF(2)/4π2|\delta_F|^{3/2}\zeta_F(2)/4\pi^2∣δF​∣3/2ζF​(2)/4π2, so the ratio of two such volumes is, up to explicit algebraic factors, a ratio of Dedekind zeta values at 222; and Neumann and Yang showed that the Bloch invariant of a hyperbolic 333-manifold lies in a subgroup of finite Q\mathbb{Q}Q-rank determined by its invariant trace field, so that manifolds sharing an invariant trace field with a single complex place, such as an imaginary quadratic one, have rationally related volumes. Producing one irrational ratio therefore means separating two such transcendentals — a statement of the same order of difficulty as the irrationality of ζ(5)\zeta(5)ζ(5). The value of formalizing the question is not that it will be closed, but that its statement, and the elementary relations that make its naive form false, are pinned down exactly.

20 thms3 active usersReviewed
🏆Completed
Captain: tianyipeng

FLT-5: Fermats Last Theorem for n=5Textbook

A complete formal proof of Fermats Last Theorem for exponent 5: for all positive natural numbers a,b,c, a^5 + b^5 != c^5. The proof follows the classical Legendre-Dirichlet approach (1825-1830): Case 1 (5 does not divide a,b,c) is dispatched by congruences, and Case 2 (5 divides one of them) uses infinite descent through the ring Z[zeta_5]. The open hard leaf is the Z[zeta_5] PID step (flt5_zeta5_ring_witnesses).

63 thms3 active usersReviewed
Group Theory·Captain: Lucas

The Mathieu group M23 is a Galois group over QResearch Paper

Motivation

The inverse Galois problem asks whether every finite group GGG occurs as the Galois group of a finite Galois extension of Q\mathbb{Q}Q. For finite simple groups, a large part of the problem was settled by the rigidity method (Shih, Fried, Belyi, Matzat, Thompson) and its refinements such as the braid-group method. Between 1984 and 1989 this machinery realized 25 of the 26 sporadic simple groups as Galois groups over Q\mathbb{Q}Q, in fact as Galois groups of regular extensions of Q(t)\mathbb{Q}(t)Q(t). The Mathieu group M23M_{23}M23​ was the single exception.

Timeline.

  • 1985–1987: Hoyden-Siedersleben and Häfner obtained regular M23M_{23}M23​-extensions of k(t)k(t)k(t) for k=Q(−23)k=\mathbb{Q}(\sqrt{-23})k=Q(−23​) and k=Q(−7)k=\mathbb{Q}(\sqrt{-7})k=Q(−7​), by passing through M24M_{24}M24​.
  • 1996: Granboulan constructed a regular M23M_{23}M23​-extension of k(t)k(t)k(t) for every field kkk over which a certain conic has a point; that conic has no rational point.
  • 2013: Elkies computed the four complex polynomials PPP of degree 23 with Gal(P(x)−t/C(t))≅M23\mathrm{Gal}(P(x)-t/\mathbb{C}(t))\cong M_{23}Gal(P(x)−t/C(t))≅M23​; each is defined over a quartic number field.
  • 2026: Huang, Jackson, Lee, Poonen, Pries and Zhang (arXiv:2608.08538) produced an explicit regular M23M_{23}M23​-extension of Q(t)\mathbb{Q}(t)Q(t) and explicit degree-23 polynomials over Q\mathbb{Q}Q with Galois group M23M_{23}M23​, completing the program for the sporadic groups.

Setting

S23S_{23}S23​ denotes the group of permutations of the 23 points {1,…,23}\{1,\dots,23\}{1,…,23}, acting on the left, so (στ)(x)=σ(τ(x))(\sigma\tau)(x)=\sigma(\tau(x))(στ)(x)=σ(τ(x)). The paper fixes three explicit permutations

g1=(1,11)(2,23)(3,8)(4,16)(5,21)(7,20)(15,19)(18,22),g_1=(1,11)(2,23)(3,8)(4,16)(5,21)(7,20)(15,19)(18,22),g1​=(1,11)(2,23)(3,8)(4,16)(5,21)(7,20)(15,19)(18,22), g2=(1,2,11,10,16,9,6,3,23,19,20,14,21,17,4,8,22,5,18,15,13,7,12),g_2=(1,2,11,10,16,9,6,3,23,19,20,14,21,17,4,8,22,5,18,15,13,7,12),g2​=(1,2,11,10,16,9,6,3,23,19,20,14,21,17,4,8,22,5,18,15,13,7,12), g3=(1,2,3,4,10,11,12,7,19,18,8,6,9,16,17,21,22,5,14,20,13,15,23),g_3=(1,2,3,4,10,11,12,7,19,18,8,6,9,16,17,21,22,5,14,20,13,15,23),g3​=(1,2,3,4,10,11,12,7,19,18,8,6,9,16,17,21,22,5,14,20,13,15,23),

and the Mathieu group M23M_{23}M23​ is the subgroup of S23S_{23}S23​ generated by g1g_1g1​ and g2g_2g2​. A GGG-extension of a field kkk is a Galois extension L/kL/kL/k together with an isomorphism Gal(L/k)≅G\mathrm{Gal}(L/k)\cong GGal(L/k)≅G. A finite extension LLL of Q(t)\mathbb{Q}(t)Q(t) is regular if it contains no nontrivial algebraic extension of Q\mathbb{Q}Q.

For a triple of conjugacy classes (C1,C2,C3)(C_1,C_2,C_3)(C1​,C2​,C3​) of M23M_{23}M23​, the set Σc\Sigma_cΣc​ consists of triples (h1,h2,h3)∈C1×C2×C3(h_1,h_2,h_3)\in C_1\times C_2\times C_3(h1​,h2​,h3​)∈C1​×C2​×C3​ with h1h2h3=1h_1h_2h_3=1h1​h2​h3​=1 that generate M23M_{23}M23​, and the Nielsen class Nic\mathrm{Ni}_cNic​ is the set of orbits of Σc\Sigma_cΣc​ under simultaneous conjugation by M23M_{23}M23​. The paper works with the classes C1=2C_1=2C1​=2, C2=23AC_2=23AC2​=23A, C3=23BC_3=23BC3​=23B, represented by g1,g2,g3g_1,g_2,g_3g1​,g2​,g3​.

Formalization targets

Goal (Theorem 1.1)

∃ K/Q finite Galois with Gal(K/Q)≅M23.\exists\ K/\mathbb{Q}\ \text{finite Galois with}\ \mathrm{Gal}(K/\mathbb{Q})\cong M_{23}.∃ K/Q finite Galois with Gal(K/Q)≅M23​.

Stronger forms

  • Theorem 1.3: there is a finite Galois extension L/Q(t)L/\mathbb{Q}(t)L/Q(t), regular over Q\mathbb{Q}Q, with Gal(L/Q(t))≅M23\mathrm{Gal}(L/\mathbb{Q}(t))\cong M_{23}Gal(L/Q(t))≅M23​.
  • Examples 1.2 and 3.7: two explicit monic degree-23 polynomials in Z[x]\mathbb{Z}[x]Z[x] whose splitting fields are M23M_{23}M23​-extensions of Q\mathbb{Q}Q, unramified outside {2,3,23}\{2,3,23\}{2,3,23} and {2,7,23}\{2,7,23\}{2,7,23} respectively.

Supporting milestones

Facts about M23M_{23}M23​ stated in §3 (order 10,200,96010{,}200{,}96010,200,960, simplicity, 4-transitivity, 17 conjugacy classes), the membership (g1,g2,g3)∈Σc(g_1,g_2,g_3)\in\Sigma_c(g1​,g2​,g3​)∈Σc​, the count ∣Nic∣=7|\mathrm{Ni}_c|=7∣Nic​∣=7, the Riemann–Hurwitz count of Lemma 3.1, the hyperbolic triangle of Lemma 3.2, and the group-theoretic step of Corollary 3.6 (M23M_{23}M23​ is not normal in any strictly larger subgroup of S23S_{23}S23​).

Significance

Theorem 1.1 removes the last sporadic exception: combined with earlier work, every sporadic simple group is the Galois group of a regular extension of Q(t)\mathbb{Q}(t)Q(t) (Corollary 1.4), and hence occurs as a Galois group over every number field, in infinitely many mutually independent ways.

The paper's proof is computer-assisted: Belyi maps were computed numerically and the final claims were certified in Magma and PARI/GP. No machine-checked proof in a proof assistant is known. A formal proof of the explicit-polynomial statements (Examples 1.2 and 3.7) would give an independent, kernel-checked certificate of Theorem 1.1. The group-theoretic milestones (order, simplicity, transitivity, class counts, the Nielsen count) are reusable for any later work on M23M_{23}M23​ or the other Mathieu groups.

Difficulty

The rigidity method fails for M23M_{23}M23​: for every GQG_{\mathbb{Q}}GQ​-stable triple of conjugacy classes the Nielsen class has size different from 111, so no rational point of a Hurwitz space is forced. The smallest positive size, ∣Nic∣=7|\mathrm{Ni}_c|=7∣Nic​∣=7 for {2,23A,23B}\{2,23A,23B\}{2,23A,23B}, leaves seven covers, and the fact that one of them has field of moduli Q\mathbb{Q}Q was found by explicit computation, with no conceptual explanation. Certifying that a particular degree-23 polynomial has Galois group exactly M23M_{23}M23​ requires both a lower bound (the group contains M23M_{23}M23​, for instance via cycle types and the classification of transitive groups of degree 23) and an upper bound (the group is contained in a conjugate of M23M_{23}M23​, for instance via resolvents or reduction modulo primes). Neither bound is a finite check inside current Mathlib.

Formalization scope

S23S_{23}S23​ is Equiv.Perm (Fin 23). The paper's point kkk corresponds to k - 1 : Fin 23, and permutations compose as functions, which matches the paper's left-action convention. M23M_{23}M23​ is defined as Subgroup.closure {g₁, g₂}, so it is not described up to isomorphism. Galois groups are groups of field automorphisms, K ≃ₐ[ℚ] K, and a GGG-extension is recorded as a group isomorphism with M23M_{23}M23​. Q(t)\mathbb{Q}(t)Q(t) is RatFunc ℚ. "Unramified outside SSS" for a number field is encoded as "every prime dividing the absolute discriminant lies in SSS", which is equivalent by Dedekind's discriminant theorem. The Nielsen class is the set of orbits of Σc\Sigma_cΣc​ under simultaneous conjugation. Hyperbolic angles in the Poincaré disk are defined through the hyperbolic law of cosines.

The geometric statements Lemma 3.3 and Proposition 3.5 concern the numerically computed curve XCX_{\mathbb{C}}XC​ and the polynomial F(T,V)F(T,V)F(T,V), which the paper does not print, so they are not milestones. Lemma 3.1 enters only through its Riemann–Hurwitz count.

Contributions are welcome on decidable certificates for permutation-group facts in Lean, on Galois-group certification for explicit polynomials, and on Hilbert irreducibility.

Selected references

  • X. Huang, B. Jackson, K.-H. Lee, B. Poonen, R. Pries, S. Zhang, The Mathieu group M23M_{23}M23​ is a Galois group over Q\mathbb{Q}Q, 2026. https://arxiv.org/abs/2608.08538
  • G. Malle, B. H. Matzat, Inverse Galois Theory, 2nd ed., Springer, 2018.
  • J.-P. Serre, Topics in Galois Theory, Jones and Bartlett, 1992.
  • N. D. Elkies, The complex polynomials P(x)P(x)P(x) with Gal(P(x)−t)≅M23\mathrm{Gal}(P(x)-t)\cong M_{23}Gal(P(x)−t)≅M23​, ANTS X, Open Book Series 1, 2013.
  • M. D. Fried, H. Völklein, The inverse Galois problem and rational points on moduli spaces, Math. Ann. 290 (1991), 771–800.
18 thms2 active usersReviewed
🏆Completed
Captain: xuanji

Every Odd Number Greater Than 1 is the Sum of at Most 351 PrimesResearch Paper

Motivation

Schnirelmann showed around 1930, by elementary means, that some absolute constant kkk makes every integer n>1n > 1n>1 a sum of at most kkk primes. For odd nnn:

  • Schnirelmann (1930s): some finite kkk, by elementary methods.
  • Vinogradov (1937): every sufficiently large odd integer is a sum of three primes.
  • Ramaré (1995): every even integer is a sum of at most six primes, so every odd n>1n > 1n>1 is a sum of at most seven. (Ann. Sc. Norm. Super. Pisa, 1995)
  • Tao (2014): at most five primes. (arXiv:1201.6656)
  • Helfgott (2013): every odd n>5n > 5n>5 is a sum of three primes. (arXiv:1312.7748)

The campaign's first proved value, 100 001100\,001100001, came from Schnirelmann's method with every constant written out. This entry records a sharper value, 351351351, from the same elementary circle of ideas.

Setting

A representation of nnn as a sum of at most kkk primes is a finite multiset of primes summing to nnn with at most kkk elements counted with multiplicity. The Schnirelmann density of A⊆Z≥0A \subseteq \mathbb{Z}_{\ge 0}A⊆Z≥0​ is σ(A)=inf⁡N≥1∣A∩{1,…,N}∣/N\sigma(A) = \inf_{N \ge 1} |A \cap \{1, \dots, N\}|/Nσ(A)=infN≥1​∣A∩{1,…,N}∣/N (Mathlib: schnirelmannDensity).

Formalization target

Goal

∀n∈N,n odd, n>1  ⟹  ∃ s multiset of primes, ∣s∣≤351, ∑s=n.\forall n \in \mathbb{N},\quad n \text{ odd},\ n > 1 \implies \exists\, s \text{ multiset of primes},\ |s| \le 351,\ \textstyle\sum s = n.∀n∈N,n odd, n>1⟹∃s multiset of primes, ∣s∣≤351, ∑s=n.

This is the campaign template with the value 351351351 filled in. The source proves the stronger statement that every odd n≥703n \ge 703n≥703 is a sum of exactly 351351351 primes; the at-most form for all odd n>1n > 1n>1 follows.

How the bound arises

It uses the same density estimate as the companion 485485485 entry, σ(A)≥1/175\sigma(A) \ge 1/175σ(A)≥1/175 for A=B+BA = B + BA=B+B with B={(p−3)/2:p odd prime}B = \{(p-3)/2 : p \text{ odd prime}\}B={(p−3)/2:p odd prime} (explicit Selberg sieve, weighted first moment, eighth moment of the singular-series factor, Hölder). It then replaces Schnirelmann's sumset inequality by Mann's theorem, σ(D+E)≥min⁡{1,σ(D)+σ(E)}\sigma(D + E) \ge \min\{1, \sigma(D) + \sigma(E)\}σ(D+E)≥min{1,σ(D)+σ(E)} for sets containing 000:

  1. Mann's theorem gives 175A=Z≥0175A = \mathbb{Z}_{\ge 0}175A=Z≥0​, so 350B=Z≥0350B = \mathbb{Z}_{\ge 0}350B=Z≥0​.
  2. For odd n≥3K=1053n \ge 3K = 1053n≥3K=1053, write (n−3K)/2(n - 3K)/2(n−3K)/2 as a sum of 350350350 elements of BBB and add one more 333, giving K=351K = 351K=351 primes.
  3. For 703≤n<1053703 \le n < 1053703≤n<1053, use n−2Kn - 2Kn−2K threes and 3K−n3K - n3K−n twos.

Significance

The bound is far weaker than Tao's 555 or Helfgott's 333, but it rests on an elementary argument with no "sufficiently large" threshold and no prime number theorem, so it is a realistic target for a complete formalization and a large step down from 100 001100\,001100001. Reusable components:

  1. Explicit Chebyshev-type lower bound for π(y)\pi(y)π(y).
  2. Explicit Selberg upper-bound sieve for r(s)r(s)r(s) with threshold e100e^{100}e100.
  3. The eighth-moment bound ∑s≤xC(s)8≤800 000 x\sum_{s \le x} C(s)^8 \le 800\,000\,x∑s≤x​C(s)8≤800000x.
  4. Mann's theorem (αβ\alpha\betaαβ theorem) on Schnirelmann density.

Formalization scope

The Lean statement is the campaign template verbatim with 351351351 in place of the value. Mathlib already has schnirelmannDensity, the Λ² Selberg sieve setup (Mathlib/NumberTheory/SelbergSieve.lean) and central-binomial bounds.

Selected references

  • P. Pollack, Not Always Buried Deep, AMS, 2009, Chapter 6, §6. https://www.pollack-math.net/NABDofficial.pdf
  • K. S. Kedlaya, Notes on Analytic Number Theory, Chapter 13, "The Selberg sieve". https://kskedlaya.org/ant/chap-selberg.html
  • O. Ramaré, On Šnirel'man's constant, Ann. Sc. Norm. Super. Pisa (4) 22 (1995), 645–706.
  • T. Tao, Every odd number greater than 1 is the sum of at most five primes, Math. Comp. 83 (2014). https://arxiv.org/abs/1201.6656
  • H. A. Helfgott, The ternary Goldbach conjecture is true, 2013. https://arxiv.org/abs/1312.7748
  • Explicit improvement of the 100 001100\,001100001 constant (unpublished AI-assisted calculation, October 2026). Source of the constant 351351351; not peer reviewed.
6 thms2 active usersReviewed
🏆Completed
Captain: xuanji

Every Odd Number Greater Than 1 is the Sum of at Most 485 PrimesResearch Paper

Motivation

Schnirelmann showed around 1930, by elementary means, that some absolute constant kkk makes every integer n>1n > 1n>1 a sum of at most kkk primes. For odd nnn:

  • Schnirelmann (1930s): some finite kkk, by elementary methods.
  • Vinogradov (1937): every sufficiently large odd integer is a sum of three primes.
  • Ramaré (1995): every even integer is a sum of at most six primes, so every odd n>1n > 1n>1 is a sum of at most seven. (Ann. Sc. Norm. Super. Pisa, 1995)
  • Tao (2014): at most five primes. (arXiv:1201.6656)
  • Helfgott (2013): every odd n>5n > 5n>5 is a sum of three primes. (arXiv:1312.7748)

The campaign's first proved value, 100 001100\,001100001, came from Schnirelmann's method with every constant written out. This entry records a sharper value, 485485485, from the same elementary circle of ideas.

Setting

A representation of nnn as a sum of at most kkk primes is a finite multiset of primes summing to nnn with at most kkk elements counted with multiplicity. The Schnirelmann density of A⊆Z≥0A \subseteq \mathbb{Z}_{\ge 0}A⊆Z≥0​ is σ(A)=inf⁡N≥1∣A∩{1,…,N}∣/N\sigma(A) = \inf_{N \ge 1} |A \cap \{1, \dots, N\}|/Nσ(A)=infN≥1​∣A∩{1,…,N}∣/N (Mathlib: schnirelmannDensity).

Formalization target

Goal

∀n∈N,n odd, n>1  ⟹  ∃ s multiset of primes, ∣s∣≤485, ∑s=n.\forall n \in \mathbb{N},\quad n \text{ odd},\ n > 1 \implies \exists\, s \text{ multiset of primes},\ |s| \le 485,\ \textstyle\sum s = n.∀n∈N,n odd, n>1⟹∃s multiset of primes, ∣s∣≤485, ∑s=n.

This is the campaign template with the value 485485485 filled in. The source proves the stronger statement that every odd n≥971n \ge 971n≥971 is a sum of exactly 485485485 primes; the at-most form for all odd n>1n > 1n>1 follows.

How the bound arises

It keeps the explicit Selberg sieve and Schnirelmann's original sumset inequality from the 100 001100\,001100001 entry, and improves only the density estimate:

  1. Lower sieve threshold. With z=s/(log⁡s)2z = \sqrt{s}/(\log s)^2z=s​/(logs)2 the sieve gives r(s)≤13 C(s) s/(log⁡s)2r(s) \le 13\,C(s)\,s/(\log s)^2r(s)≤13C(s)s/(logs)2 for even s≥e100s \ge e^{100}s≥e100, where C(s)=∏p∣s(1+p/(p−1)2)C(s) = \prod_{p \mid s}\bigl(1 + p/(p-1)^2\bigr)C(s)=∏p∣s​(1+p/(p−1)2).
  2. Weighted first moment. Counting over the whole triangle p+q≤xp + q \le xp+q≤x and weighting by (log⁡s)2/s(\log s)^2/s(logs)2/s gives ∑e100<s≤xr(s)(log⁡s)2/s≥43100x\sum_{e^{100} < s \le x} r(s)(\log s)^2/s \ge \tfrac{43}{100}x∑e100<s≤x​r(s)(logs)2/s≥10043​x for x≥e200x \ge e^{200}x≥e200.
  3. Eighth moment of CCC. An Euler-product estimate (primes 3,5,73, 5, 73,5,7 handled individually, the tail bounded at once) gives ∑s≤x, 2∣sC(s)8≤800 000 x\sum_{s \le x,\, 2 \mid s} C(s)^8 \le 800\,000\,x∑s≤x,2∣s​C(s)8≤800000x.
  4. Hölder instead of Cauchy–Schwarz. This yields #{s≤x:r(s)>0}≥x/345\#\{s \le x : r(s) > 0\} \ge x/345#{s≤x:r(s)>0}≥x/345 for x≥e200x \ge e^{200}x≥e200, and with Chebyshev's bound for smaller scales, σ(A)≥1/175\sigma(A) \ge 1/175σ(A)≥1/175 for A=B+BA = B + BA=B+B, B={(p−3)/2}B = \{(p-3)/2\}B={(p−3)/2}.
  5. Schnirelmann's inequality with m=121m = 121m=121 (since (174/175)121<1/2(174/175)^{121} < 1/2(174/175)121<1/2) gives 242A=Z≥0242A = \mathbb{Z}_{\ge 0}242A=Z≥0​, hence K=4m+1=485K = 4m + 1 = 485K=4m+1=485.

Only Chebyshev-type prime bounds, the Selberg upper-bound sieve, Hölder's inequality and Schnirelmann's inequality are used.

Significance

The bound is far weaker than Tao's 555 or Helfgott's 333, but it rests on an elementary argument with no "sufficiently large" threshold and no prime number theorem, so it is a realistic target for a complete formalization and a large step down from 100 001100\,001100001. Reusable components:

  1. Explicit Chebyshev-type lower bound for π(y)\pi(y)π(y).
  2. Explicit Selberg upper-bound sieve for r(s)r(s)r(s) with threshold e100e^{100}e100.
  3. The eighth-moment bound ∑s≤xC(s)8≤800 000 x\sum_{s \le x} C(s)^8 \le 800\,000\,x∑s≤x​C(s)8≤800000x.
  4. Schnirelmann's inequality σ(D+E)≥σ(D)+σ(E)−σ(D)σ(E)\sigma(D+E) \ge \sigma(D)+\sigma(E)-\sigma(D)\sigma(E)σ(D+E)≥σ(D)+σ(E)−σ(D)σ(E).

Formalization scope

The Lean statement is the campaign template verbatim with 485485485 in place of the value. Mathlib already has schnirelmannDensity, the Λ² Selberg sieve setup (Mathlib/NumberTheory/SelbergSieve.lean) and central-binomial bounds.

Selected references

  • P. Pollack, Not Always Buried Deep, AMS, 2009, Chapter 6, §6. https://www.pollack-math.net/NABDofficial.pdf
  • K. S. Kedlaya, Notes on Analytic Number Theory, Chapter 13, "The Selberg sieve". https://kskedlaya.org/ant/chap-selberg.html
  • O. Ramaré, On Šnirel'man's constant, Ann. Sc. Norm. Super. Pisa (4) 22 (1995), 645–706.
  • T. Tao, Every odd number greater than 1 is the sum of at most five primes, Math. Comp. 83 (2014). https://arxiv.org/abs/1201.6656
  • H. A. Helfgott, The ternary Goldbach conjecture is true, 2013. https://arxiv.org/abs/1312.7748
  • Explicit improvement of the 100 001100\,001100001 constant (unpublished AI-assisted calculation, October 2026). Source of the constant 485485485; not peer reviewed.
6 thms2 active usersReviewed
🏆Completed
Captain: xuanji

Every Odd Number Greater Than 1 is the Sum of at Most 973 PrimesResearch Paper

Motivation

Schnirelmann showed around 1930, by elementary means, that some absolute constant kkk makes every integer n>1n > 1n>1 a sum of at most kkk primes. For odd nnn:

  • Schnirelmann (1930s): some finite kkk, by elementary methods.
  • Vinogradov (1937): every sufficiently large odd integer is a sum of three primes.
  • Ramaré (1995): every even integer is a sum of at most six primes, so every odd n>1n > 1n>1 is a sum of at most seven. (Ann. Sc. Norm. Super. Pisa, 1995)
  • Tao (2014): at most five primes. (arXiv:1201.6656)
  • Helfgott (2013): every odd n>5n > 5n>5 is a sum of three primes. (arXiv:1312.7748)

The campaign's first proved value, 100 001100\,001100001, came from Schnirelmann's method with every constant written out. This entry records a sharper value, 973973973, from the same elementary circle of ideas.

Setting

A representation of nnn as a sum of at most kkk primes is a finite multiset of primes summing to nnn with at most kkk elements counted with multiplicity. The Schnirelmann density of A⊆Z≥0A \subseteq \mathbb{Z}_{\ge 0}A⊆Z≥0​ is σ(A)=inf⁡N≥1∣A∩{1,…,N}∣/N\sigma(A) = \inf_{N \ge 1} |A \cap \{1, \dots, N\}|/Nσ(A)=infN≥1​∣A∩{1,…,N}∣/N (Mathlib: schnirelmannDensity).

Formalization target

Goal

∀n∈N,n odd, n>1  ⟹  ∃ s multiset of primes, ∣s∣≤973, ∑s=n.\forall n \in \mathbb{N},\quad n \text{ odd},\ n > 1 \implies \exists\, s \text{ multiset of primes},\ |s| \le 973,\ \textstyle\sum s = n.∀n∈N,n odd, n>1⟹∃s multiset of primes, ∣s∣≤973, ∑s=n.

This is the campaign template with the value 973973973 filled in. The source proves the stronger statement that every odd n≥1947n \ge 1947n≥1947 is a sum of exactly 973973973 primes; the at-most form for all odd n>1n > 1n>1 follows.

How the bound arises

It keeps the explicit Selberg sieve and Schnirelmann's original sumset inequality from the 100 001100\,001100001 entry, and improves the density estimate:

  1. Lower sieve threshold. With z=s/(log⁡s)2z = \sqrt{s}/(\log s)^2z=s​/(logs)2 the sieve gives r(s)≤13 C(s) s/(log⁡s)2r(s) \le 13\,C(s)\,s/(\log s)^2r(s)≤13C(s)s/(logs)2 for even s≥e100s \ge e^{100}s≥e100.
  2. Weighted first moment of at least 43100x\tfrac{43}{100}x10043​x for x≥e200x \ge e^{200}x≥e200.
  3. Fourth moment of CCC. An Euler-product estimate gives ∑s≤x, 2∣sC(s)4≤400 x\sum_{s \le x,\, 2\mid s} C(s)^4 \le 400\,x∑s≤x,2∣s​C(s)4≤400x.
  4. Hölder then yields σ(A)≥1/350\sigma(A) \ge 1/350σ(A)≥1/350 for A=B+BA = B + BA=B+B, B={(p−3)/2}B = \{(p-3)/2\}B={(p−3)/2}.
  5. Schnirelmann's inequality with m=243m = 243m=243 (the least mmm with (349/350)m<1/2(349/350)^m < 1/2(349/350)m<1/2) gives K=4m+1=973K = 4m + 1 = 973K=4m+1=973.

Significance

The bound is far weaker than Tao's 555 or Helfgott's 333, but it rests on an elementary argument with no "sufficiently large" threshold and no prime number theorem, so it is a realistic target for a complete formalization and a large step down from 100 001100\,001100001. Reusable components:

  1. Explicit Chebyshev-type lower bound for π(y)\pi(y)π(y).
  2. Explicit Selberg upper-bound sieve for r(s)r(s)r(s).
  3. Moment bounds for the singular-series factor C(s)C(s)C(s).
  4. Schnirelmann's inequality σ(D+E)≥σ(D)+σ(E)−σ(D)σ(E)\sigma(D+E) \ge \sigma(D)+\sigma(E)-\sigma(D)\sigma(E)σ(D+E)≥σ(D)+σ(E)−σ(D)σ(E).

Formalization scope

The Lean statement is the campaign template verbatim with 973973973 in place of the value. Mathlib already has schnirelmannDensity, the Λ² Selberg sieve setup (Mathlib/NumberTheory/SelbergSieve.lean) and central-binomial bounds.

Selected references

  • P. Pollack, Not Always Buried Deep, AMS, 2009, Chapter 6, §6. https://www.pollack-math.net/NABDofficial.pdf
  • K. S. Kedlaya, Notes on Analytic Number Theory, Chapter 13, "The Selberg sieve". https://kskedlaya.org/ant/chap-selberg.html
  • O. Ramaré, On Šnirel'man's constant, Ann. Sc. Norm. Super. Pisa (4) 22 (1995), 645–706.
  • T. Tao, Every odd number greater than 1 is the sum of at most five primes, Math. Comp. 83 (2014). https://arxiv.org/abs/1201.6656
  • H. A. Helfgott, The ternary Goldbach conjecture is true, 2013. https://arxiv.org/abs/1312.7748
  • Explicit improvement of the 100 001100\,001100001 constant (unpublished AI-assisted calculation, October 2026). Source of the constant 973973973; not peer reviewed.
6 thms2 active usersReviewed
Captain: xuanji

The irrationality measure of π is at most 7.103205334138 (Zeilberger–Zudilin 2020)Research Paper

Motivation

The irrationality measure μ(π)\mu(\pi)μ(π) is the supremum of the μ\muμ for which ∣π−p/q∣<q−μ|\pi - p/q| < q^{-\mu}∣π−p/q∣<q−μ has infinitely many rational solutions p/qp/qp/q. Every irrational number has μ≥2\mu \ge 2μ≥2 (Dirichlet), almost every real number has μ=2\mu = 2μ=2, and it is conjectured that μ(π)=2\mu(\pi) = 2μ(π)=2. Known upper bounds:

  • Mahler (1953): 424242, the first proof that π\piπ is not a Liouville number.
  • Mignotte (1974): 20.620.620.6.
  • Chudnovsky (1982): 19.8899944…19.8899944\ldots19.8899944…
  • Rhin–Viola (1993): 14.79707414.79707414.797074.
  • Hata (1993): 8.016045…8.016045\ldots8.016045…
  • Salikhov (2008): 7.606308…7.606308\ldots7.606308…
  • Zeilberger–Zudilin (2020): 7.103205334137…7.103205334137\ldots7.103205334137…, the current record.

The campaign's first proved value is Mahler's 424242. This entry records Zeilberger–Zudilin's bound.

Formalization target

The campaign template with the value 7.1032053341387.1032053341387.103205334138 filled in: PiIrrationality.UpperBound (7.103205334138 : ℝ), i.e. μ(π)≤7.103205334138\mu(\pi) \le 7.103205334138μ(π)≤7.103205334138.

Value. The paper states its bound as 7.103205334137…7.103205334137\ldots7.103205334137…, a truncation. This entry rounds the last digit up to 7.1032053341387.1032053341387.103205334138 so that the goal follows from the published proof.

How the bound arises

Zeilberger and Zudilin modify Salikhov's integrals and use the Almkvist–Zeilberger algorithm (creative telescoping) to find recurrences for them, then optimise the arithmetic and analytic estimates. This is the current record.

Significance

Each step down the list replaces Mahler's approximations with a sharper family. Formalizing 7.1032053341387.1032053341387.103205334138 would build reusable explicit machinery: integral constructions of rational approximations to π\piπ, bounds on their common denominators via prime-number estimates, and the standard lemma turning a sequence of good approximations into an irrationality-measure bound.

Selected references

  • D. Zeilberger, W. Zudilin, The irrationality measure of π\piπ is at most 7.103205334137…7.103205334137\ldots7.103205334137…, Moscow J. Combin. Number Theory 9 (2020), no. 4, 407–419. https://arxiv.org/abs/1912.06345
  • K. Mahler, On the approximation of π\piπ, Indag. Math. 15 (1953), 30–42.
  • F. Beukers, A rational approach to π\piπ, Nieuw Arch. Wiskd. (5) 1 (2000), 372–379.
  • Source table: https://teorth.github.io/optimizationproblems/constants/7a.html
2 thms2 active usersReviewed
🏆Completed
Captain: xuanji

The irrationality measure of π is at most 19.8899945 (Chudnovsky 1982)Research Paper

Motivation

The irrationality measure μ(π)\mu(\pi)μ(π) is the supremum of the μ\muμ for which ∣π−p/q∣<q−μ|\pi - p/q| < q^{-\mu}∣π−p/q∣<q−μ has infinitely many rational solutions p/qp/qp/q. Every irrational number has μ≥2\mu \ge 2μ≥2 (Dirichlet), almost every real number has μ=2\mu = 2μ=2, and it is conjectured that μ(π)=2\mu(\pi) = 2μ(π)=2. Known upper bounds:

  • Mahler (1953): 424242, the first proof that π\piπ is not a Liouville number.
  • Mignotte (1974): 20.620.620.6.
  • Chudnovsky (1982): 19.8899944…19.8899944\ldots19.8899944…
  • Rhin–Viola (1993): 14.79707414.79707414.797074.
  • Hata (1993): 8.016045…8.016045\ldots8.016045…
  • Salikhov (2008): 7.606308…7.606308\ldots7.606308…
  • Zeilberger–Zudilin (2020): 7.103205334137…7.103205334137\ldots7.103205334137…, the current record.

The campaign's first proved value is Mahler's 424242. This entry records Chudnovsky's bound.

Formalization target

The campaign template with the value 19.889994519.889994519.8899945 filled in: PiIrrationality.UpperBound (19.8899945 : ℝ), i.e. μ(π)≤19.8899945\mu(\pi) \le 19.8899945μ(π)≤19.8899945.

Value. The bound is quoted in the literature as 19.8899944…19.8899944\ldots19.8899944… (e.g. Hata 1993), a truncation. This entry rounds the last digit up to 19.889994519.889994519.8899945 so that the goal follows from the published constant.

How the bound arises

Chudnovsky determined the exact asymptotic behaviour of the Hermite-type contour integrals 12πi∮(n!z(z−1)⋯(z−n))kewz dz\frac{1}{2\pi i}\oint \left(\frac{n!}{z(z-1)\cdots(z-n)}\right)^k e^{wz}\,dz2πi1​∮(z(z−1)⋯(z−n)n!​)kewzdz behind Mahler's approximations, which sharpens the resulting exponent.

Significance

Each step down the list replaces Mahler's approximations with a sharper family. Formalizing 19.889994519.889994519.8899945 would build reusable explicit machinery: integral constructions of rational approximations to π\piπ, bounds on their common denominators via prime-number estimates, and the standard lemma turning a sequence of good approximations into an irrationality-measure bound.

Selected references

  • G. V. Chudnovsky, Hermite–Padé approximations to exponential functions and elementary estimates of the measure of irrationality of π\piπ, Lecture Notes in Math. 925, Springer (1982), 299–322.
  • K. Mahler, On the approximation of π\piπ, Indag. Math. 15 (1953), 30–42.
  • F. Beukers, A rational approach to π\piπ, Nieuw Arch. Wiskd. (5) 1 (2000), 372–379.
  • Source table: https://teorth.github.io/optimizationproblems/constants/7a.html
2 thms2 active usersReviewed
🏆Completed
Captain: xuanji

The irrationality measure of π is at most 20.6 (Mignotte 1974)Research Paper

Motivation

The irrationality measure μ(π)\mu(\pi)μ(π) is the supremum of the μ\muμ for which ∣π−p/q∣<q−μ|\pi - p/q| < q^{-\mu}∣π−p/q∣<q−μ has infinitely many rational solutions p/qp/qp/q. Every irrational number has μ≥2\mu \ge 2μ≥2 (Dirichlet), almost every real number has μ=2\mu = 2μ=2, and it is conjectured that μ(π)=2\mu(\pi) = 2μ(π)=2. Known upper bounds:

  • Mahler (1953): 424242, the first proof that π\piπ is not a Liouville number.
  • Mignotte (1974): 20.620.620.6.
  • Chudnovsky (1982): 19.8899944…19.8899944\ldots19.8899944…
  • Rhin–Viola (1993): 14.79707414.79707414.797074.
  • Hata (1993): 8.016045…8.016045\ldots8.016045…
  • Salikhov (2008): 7.606308…7.606308\ldots7.606308…
  • Zeilberger–Zudilin (2020): 7.103205334137…7.103205334137\ldots7.103205334137…, the current record.

The campaign's first proved value is Mahler's 424242. This entry records Mignotte's bound.

Formalization target

The campaign template with the value 20.620.620.6 filled in: PiIrrationality.UpperBound (20.6 : ℝ), i.e. μ(π)≤20.6\mu(\pi) \le 20.6μ(π)≤20.6.

Value. The paper's abstract states ∣π−p/q∣>q−20.6|\pi - p/q| > q^{-20.6}∣π−p/q∣>q−20.6 for all q≥2q \ge 2q≥2, which gives μ(π)≤20.6\mu(\pi) \le 20.6μ(π)≤20.6 exactly as stated. The paper also proves ∣π−p/q∣>q−20|\pi - p/q| > q^{-20}∣π−p/q∣>q−20 for q≥q0q \ge q_0q≥q0​ (explicit), so μ(π)≤20\mu(\pi) \le 20μ(π)≤20 follows from the same source; this entry uses the table value 20.620.620.6.

How the bound arises

Mignotte refined Mahler's method of explicit rational approximations to π\piπ (Hermite's approximation formulae for the exponential and logarithm) and sharpened the estimates that turn their size and denominators into an irrationality measure.

Significance

Each step down the list replaces Mahler's approximations with a sharper family. Formalizing 20.620.620.6 would build reusable explicit machinery: integral constructions of rational approximations to π\piπ, bounds on their common denominators via prime-number estimates, and the standard lemma turning a sequence of good approximations into an irrationality-measure bound.

Selected references

  • M. Mignotte, Approximations rationnelles de π\piπ et quelques autres nombres, Mém. Soc. Math. France 37 (1974), 121–132. https://doi.org/10.24033/msmf.139
  • K. Mahler, On the approximation of π\piπ, Indag. Math. 15 (1953), 30–42.
  • F. Beukers, A rational approach to π\piπ, Nieuw Arch. Wiskd. (5) 1 (2000), 372–379.
  • Source table: https://teorth.github.io/optimizationproblems/constants/7a.html
2 thms2 active usersReviewed
PreviousPage 2 of 5Next

Get started

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

About Prove2Me

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

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me