Two lower bounds on a possible rational denominator
Here the enclosure is for Π
itself, rather than a tail difference. A finite prime sum gives the
lower endpoint; the polynomial prime bound gives the upper endpoint. Two
fractions of determinant one bracket this entire interval, so every
rational in it has denominator at least the sum of their denominators.
This last step is an integer identity, independent of how the fractions
were found. Appendix B prints the
smaller certificate; the linked Lean development checks the integer
comparisons for both. The bounds exclude finite denominator ranges, not
rationality.
Theorem 7.1 (a denominator bound checked by exact
integer comparisons). Let 𝑎 ∈ℤ and let 𝑏 be a positive integer. If Π =𝑎/𝑏 then 𝑏 ≥2589 >10177, and the same
floor holds for every rational equal to 𝑆.
Proof. Put 𝑐 =1229, the
number of primes below 104, and
𝐴 =∑𝑐−1𝑖=0𝑝𝑖2𝑐−𝑖−1,
so that 𝐴/2𝑐 is the 𝑐th partial sum of Π. All terms are positive, so 𝐴/2𝑐 ≤Π. For the upper bound, 𝑝𝑐+𝑗 ≤1250(𝑐 +𝑗 +1)4. Since 𝑐 ≥9, each successive ratio of (𝑐 +1 +𝑗)4 is at most (11/10)4 <3/2. Hence (𝑐 +1 +𝑗)4 ≤(𝑐 +1)4(3/2)𝑗 for all 𝑗 ≥0. The omitted dyadic terms are
therefore bounded by a geometric series with ratio 3/4, whose sum gives
𝐴2𝑐≤Π≤2𝐴+5000(𝑐+1)42𝑐+1.
The four positive integers 𝑢,𝑣,𝑢′,𝑣′ printed in
Appendix B satisfy
𝑢′𝑣−𝑢𝑣′=1,𝑢𝑣<𝐴2𝑐,2𝐴+5000(𝑐+1)42𝑐+1<𝑢′𝑣′,𝑣+𝑣′≥2589,
all four being integer comparisons after
clearing the positive denominators. If Π =𝑎/𝑏, the strict lower and upper
bounds make 𝑎𝑣 −𝑏𝑢 and 𝑏𝑢′ −𝑎𝑣′ positive integers. Each
is therefore at least 1. The
determinant identity expresses 𝑏 as
their positive integer combination of 𝑣′ and 𝑣:
𝑏=𝑏(𝑢′𝑣−𝑢𝑣′)=(𝑎𝑣−𝑏𝑢)𝑣′+(𝑏𝑢′−𝑎𝑣′)𝑣≥𝑣+𝑣′≥2589.
The only properties needed of the
bracketing fractions are their order, positive denominators and
determinant one. Their decimal expansions and the way they were found
play no part in this denominator bound. Since 589log102 >177, the floor exceeds
10177. If 𝑆 =𝑎/𝑏, apply the same argument to Π =(𝑎 +2𝑏)/𝑏. ◻
The Lean kernel evaluates the four inequalities itself, on the same
integers, after recomputing the primes below 104 by trial division; no compiled code
is trusted. Its analytic input is the same enclosure, derived from the
polynomial
bound 𝑝𝑛 ≤1250(𝑛 +1)4. The
four integers 𝑢,𝑣,𝑢′,𝑣′
have 177 and 178 digits.
Theorem 7.2 (certified continued-fraction exclusion).
Every rational equal to Π,
and hence every rational equal to 𝑆, has reduced denominator 𝑞 ≥239997, and therefore 𝑞 >1012040.
The same determinant argument gives the stronger floor 240062 below. The smaller constants in
the statement and the different continued-fraction pairs are accounted
for after the proof.
Proof. Put 𝑐 =80200, the
number of primes below 1 023 068,
and 𝐴 =∑𝑐−1𝑖=0𝑝𝑖2𝑐−𝑖−1. The
estimate in the proof of Theorem 7.1 holds
for every 𝑐 ≥9 and gives
𝐴2𝑐≤Π≤2𝐴+5000(𝑐+1)42𝑐+1.
Four positive integers 𝑢,𝑣,𝑢′,𝑣′, described after this
proof, satisfy
𝑢′𝑣−𝑢𝑣′=1,𝑢𝑣<𝐴2𝑐,2𝐴+5000(𝑐+1)42𝑐+1<𝑢′𝑣′,𝑣+𝑣′≥240062,
all four being integer comparisons after
clearing the positive denominators. If Π =𝑎/𝑏 with 𝑎 ∈ℤ and 𝑏 a positive integer, the determinant
calculation in the proof of Theorem 7.1 gives
𝑏≥𝑣+𝑣′≥240062>1012059.
The last inequality holds because 1012059 =(1059)204 1023 and
240062 =(2196)204 278,
with 1059 <2196 and 1023 <278. The reduced
denominator 𝑞 is one such 𝑏, so 𝑞 ≥240062 >239997 and 𝑞 >1012059 >1012040. If 𝑆 =𝑎/𝑏, apply the same argument to Π =(𝑎 +2𝑏)/𝑏. ◻
The integers 𝑢,𝑣,𝑢′,𝑣′
come from the continued fractions of the two endpoints 𝐴/2𝑐 and (2𝐴 +5000(𝑐 +1)4)/2𝑐+1, which agree
in their first 23407 partial
quotients. The fraction 𝑢/𝑣 is the
convergent formed by those partial quotients, and 𝑢′/𝑣′ is the mediant of 𝑢/𝑣 and the convergent before it. Only
the four comparisons enter the proof. The numerators 𝑢,𝑢′ have 12061 decimal digits and the denominators
𝑣,𝑣′ have 12060.
The Lean kernel checks every step of this proof and evaluates every
computation itself; no compiled code is trusted. It recomputes 𝑐 and 𝐴 from the definition of primality,
scanning the integers below 1 023 068 in 250 consecutive blocks and checking the
prime count and partial numerator at the end of each block. It decides
primality by trial division below 1012. On [1012,10132), which contains the rest of
the scan because 10132 =1 026 169, it uses the fact
that such an integer is prime exactly when it is coprime to the product
of the primes below 1012, and it
recomputes that product by trial division. Lean proves that the two
primality tests agree at every natural number. The kernel then decides
the four comparisons, together with 1059 <2196 and 1023 <278; the enclosure and the
determinant step are proved in general form. Lean proves the bound
𝑏 ≥240062 and 𝑏 >1012059 for both series and
derives the theorem as stated from it, using only Lean’s three standard
axioms.
An independent program
outside Lean recomputes 𝐴 by a
prime sieve and uses exact integer arithmetic. Its recorded
output has the same 23407
common partial quotients. It brackets the interval with the last two
consecutive convergents, whose denominator sum has bit length 40062, and checks the advertised bound
239997 >1012040. This pair
differs from the convergent–mediant pair in the Lean certificate above,
whose denominator sum has bit length 40063 and is at least 240062. The two bit lengths refer to
different valid pairs.
The constants 239997 and
1012040 in the statement come
from a separate computation with a bracket of width 2−80000, which used the
best-approximation property of convergents [10]; Short gives an accessible statement
and proof [19]. Its recorded
output lists the generating program and its digest and a last
convergent whose denominator has bit length 39998, hence is at least 239997; the program also checked its
arithmetic against the known expansions of 𝑒 and 𝜋. That bit length describes one
convergent only; the bound for every rational in the enclosure comes
from the certified sum 𝑣 +𝑣′
above, and 239997 >1012040
is a separate integer comparison.
These computations exclude finite denominator ranges. Irrationality
would require unbounded exclusions, not just a larger finite bound. Even
extending the prefix need not improve the bound at every step. No
irrationality-exponent estimate is asserted.