Wiki
Wiki

Research notes on every problem, and a library of the papers behind them. Built from the open erdos repository.

Updated

Erdős 699: finite prime-interval bounds for the BFT coefficient gcd

P699 research · The asymmetric BFT argument


The six coefficient-gcd bounds needed by the index-ten cofactor argument hold at much smaller numerical thresholds. No target cofactor exponent is weakened.

Write G(c,d,dm−δ)G(c,d,dm-\delta) for the coefficient gcd in BFT, with δ∈{0,1}\delta\in\{0,1\}. The following bounds are proved for every integer m≥mstartm\ge m_{\mathrm{start}}:

G(c,d,dm−δ)≥exp⁡(dmℓ).(1)G(c,d,dm-\delta)\ge\exp(dm\ell). \tag{1}
Prime pairc,dc,dExact ℓ\ellPrevious analytic startNew mstartm_{\mathrm{start}}Certified finite intervals
2,32,329,2029,201702881759/40000000001702881759/400000000010000000010000000012744127444943649436
2,52,543,3543,35280222279237/1012280222279237/10^{12}30000003000000104610461325513255
2,72,747,3247,32444329464383/1012444329464383/10^{12}100000000100000000982598256016860168
3,53,555,3455,34477383377397/1012477383377397/10^{12}3000000030000000403840383914339143
3,73,728,1728,17120244724163/250000000000120244724163/25000000000010000000010000000013075130754992849928
5,75,73,23,27614497793/156250000007614497793/1562500000010000000010000000028602286021937619376

The 231306 certified intervals cover each finite range from its new start to one less than its previous analytic start. The existing analytic bound then covers the infinite tail. All interval margins are positive; the smallest is greater than 2.9⋅10−82.9\cdot10^{-8} after dividing by dmdm. The finite calculations concern primes and Padé parameters, not binomial rows or columns.

The checker, input parameters, and certificate retain all constants, exact margins, and partition hashes.

Prime intervals that divide the coefficient gcd

Put h=c+dh=c+d. For an integer t≥1t\ge1, write

dt=kh+r,0≤r<h.dt=kh+r,\qquad 0\le r<h.

Retain tt when 2r>h2r>h. For m≥2tm\ge2t, define

Ut(m)=hm−2t,Vt(m)=mmax⁡{dk+1,c−dt−2k−1}.(2)U_t(m)=\frac{hm-2}{t},\qquad V_t(m)=m\max\left\{\frac d{k+1},\frac{c-d}{t-2k-1}\right\}. \tag{2}

The denominators in (2) are positive. BFT Lemma 5.2 and equations (5.7)--(5.8) show that every prime

Vt(m)<p≤Ut(m)V_t(m)<p\le U_t(m)

divides both coefficient gcds in (1). The upper endpoint uses the smaller value for δ=0,1\delta=0,1, and the lower endpoint uses upper bounds valid for both. These prime intervals are disjoint for different tt.

Consequently, for any finite retained set of eligible t≤m/2t\le m/2,

log⁡G(c,d,dm−δ)≥∑tmax⁡{θ(Ut(m))−θ(Vt(m)),0}.(3)\log G(c,d,dm-\delta) \ge\sum_t\max\{\theta(U_t(m))-\theta(V_t(m)),0\}. \tag{3}

BFT describes this finite prime-interval method and its extension across ranges of mm in Section 5, following Lemma 5.4.

Certifying a whole parameter interval

For an integer interval L≤m≤RL\le m\le R, keep eligible t≤min⁡(4000,L/2)t\le\min(4000,L/2). Both endpoints in (2) increase with mm, so

S(L,R)=∑tmax⁡{θ(Ut(L))−θ(Vt(R)),0}≤log⁡G(c,d,dm−δ).(4)S(L,R)=\sum_t\max\{\theta(U_t(L))-\theta(V_t(R)),0\} \le\log G(c,d,dm-\delta). \tag{4}

The sufficient check is simply

S(L,R)≥dRℓ.(5)S(L,R)\ge dR\ell. \tag{5}

All arguments of θ\theta are rounded down to integers, exactly matching its definition. The program evaluates a rigorous integer-scaled lower bound for (4) and checks (5) by integer multiplication. It subdivides an interval when this lower bound is insufficient.

The traversal runs from the right endpoint downwards. Once it encounters an unsupported singleton, all larger integers have already been covered by certified intervals. The chosen new start is one greater than that singleton. This does not claim that (1) fails at the singleton: the particular lower bound (4) was merely insufficient there.

Rigorous logarithms and the sieve

The sieve retains every prime at most 10810^8, a total of 57614555761455 primes. No numerical approximation is used to decide primality.

For logarithmic weights, let Q=232Q=2^{32} and M=214M=2^{14}. Directed interval arithmetic gives integer lower and upper bounds for

Qlog⁡2,Qlog⁡(k/M),M≤k≤2M.Q\log2,\qquad Q\log(k/M),\qquad M\le k\le2M.

For each prime pp, calculate exactly

b=⌊log⁡2p⌋,k=⌊pM2b⌋.b=\lfloor\log_2p\rfloor,\qquad k=\left\lfloor\frac{pM}{2^b}\right\rfloor.

Then

blog⁡2+log⁡(k/M)≤log⁡p≤blog⁡2+log⁡((k+1)/M).(6)b\log2+\log(k/M)\le\log p \le b\log2+\log((k+1)/M). \tag{6}

The program obtains integer bounds for Qlog⁡pQ\log p by adding the corresponding table entries. Both prefix sums are integer arrays, with explicit bounds preventing overflow. When both endpoints in a prime interval lie below the sieve limit, subtracting two entries of the lower prefix sums the lower weights of exactly the primes in that interval. This avoids losing the errors from primes below its lower endpoint.

The sieve, both logarithm tables, and both prefix sums have retained SHA-256 hashes. The real arithmetic used to construct the logarithm table is directed interval arithmetic; every subsequent certificate comparison is an integer calculation.

Larger prime arguments and the infinite parameter tail

Above the sieve limit, the checker combines

x−2.072x<θ(x)<x(1≤x≤1011)x-2.072\sqrt x<\theta(x)<x\qquad(1\le x\le10^{11})

from BFT Lemma 5.4 with the following retained bounds from Dusart's Theorem 5.2:

∣θ(x)−x∣<ηx/log⁡kx.|\theta(x)-x|<\eta x/\log^k x.
kkη\etaValid for x≥x\ge
33118996780389967803
221/201/20122568683122568683
3339/5039/50158822621158822621
111/10001/1000908994923908994923
221/1001/10077131338537713133853

The retained Dusart reading copy records the theorem and its exact thresholds. Integer square roots and the lower logarithmic table bounds make these inequalities rigorous at every requested argument. All arguments used in the finite certificate are below 101110^{11}, where the simpler upper bound θ(x)<x\theta(x)<x applies.

Finally, the checker regenerates each original infinite-tail certificate using the rational Dusart gcd calculation and checks that it supplies exactly the displayed ℓ\ell. Thus the finite interval calculations and the existing analytic tail together prove (1) for every integer mm at or above the new start.

To reproduce all six finite certificates:

sh
mkdir -p tmp/erdos699
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_bft_prime_bootstrap.py \
  wiki/research/erdos_699/computations/erdos_699_bft_prime_bootstrap_parameters.json \
  --output tmp/erdos699/erdos699_prime_bootstrap.json

The data include elapsed times, which can differ between runs. The exact bounds, interval counts and partition hashes are deterministic.