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: fully numerical thresholds for unequal Padé blocks

P699 research · Unequal Padé blocks


The twenty pair alternatives in the index-twelve ramp certificate have fully numerical starting heights. Combining the uniform coefficient-gcd estimates with a finite prime-interval bootstrap gives

n<103 071 170(i=12).(1)\boxed{n<10^{3\,071\,170}\qquad(i=12).} \tag{1}

The completed index-twelve proof excludes the entire finite range below this stopping bound. Every constant used below is evaluated or given by an explicit rational formula; no unevaluated asymptotic constant remains.

Finite factorial and integral bounds

Fix positive directions (λ,μ,κ)(\lambda,\mu,\kappa) and write H=λ+μ+κH=\lambda+\mu+\kappa, L0=λ+μL_0=\lambda+\mu, and L1=μ+κL_1=\mu+\kappa. The two adjacent rectangles are

A=λm+ε,B=μm−1−ε,C=κm+ε,ε∈{0,1}.A=\lambda m+\varepsilon,\quad B=\mu m-1-\varepsilon,\quad C=\kappa m+\varepsilon, \qquad\varepsilon\in\{0,1\}.

Assume m≥2m\ge2. Put

S=Hlog⁡H−λlog⁡λ−μlog⁡μ−κlog⁡κ,S=H\log H-\lambda\log\lambda-\mu\log\mu-\kappa\log\kappa, F=max⁡{μ,(H+1)μ2λκ}.F=\max\left\{\mu,\frac{(H+1)\mu^2}{\lambda\kappa}\right\}.

The multinomial theorem gives

(Hm)!(λm)!(μm)!(κm)!≤eSm.\frac{(Hm)!}{(\lambda m)!(\mu m)!(\kappa m)!}\le e^{Sm}.

For ε=0\varepsilon=0, the integral prefactor is exactly μm\mu m times this multinomial coefficient. For ε=1\varepsilon=1, its ratio to the multinomial coefficient is

(Hm+1)(μm)(μm−1)(λm+1)(κm+1)≤Fm.\frac{(Hm+1)(\mu m)(\mu m-1)}{(\lambda m+1)(\kappa m+1)}\le Fm.

Thus both prefactors are at most FmeSmFm e^{Sm}.

For the seed aP−bQ=D0>0aP-bQ=D_0>0, put T=aPT=aP and z=D0/Tz=D_0/T. Define

MQ=max⁡0<u<1uμ(1−u)κ(1−(1−z)u)λ,M_Q=\max_{0<u<1}u^\mu(1-u)^\kappa(1-(1-z)u)^\lambda, ME=max⁡0<u<1uλ(1−u)κ(1−zu)μ,M_E=\max_{0<u<1}u^\lambda(1-u)^\kappa(1-zu)^\mu,

and q=S+log⁡MQq=S+\log M_Q, e=S+log⁡MEe=S+\log M_E. These are precisely the two saddle-point quantities in the unequal-block lemma.

In either integral, remove the (m−2)(m-2)-th power of the corresponding base integrand. Every remaining exponent is nonnegative, even when μ=1\mu=1 and ε=1\varepsilon=1. Its remaining integrand is consequently at most one. Therefore the explicit bounds

∣QA,B,C(z)∣≤FmMQ−2eqm,∣EA,B,C(z)∣≤FmME−2eem(2)|Q_{A,B,C}(z)|\le FmM_Q^{-2}e^{qm},\qquad |E_{A,B,C}(z)|\le FmM_E^{-2}e^{em} \tag{2}

hold for both rectangles. No Stirling remainder or limiting argument is needed in (2).

A numerical coefficient-gcd starting index

The uniform gcd checker uses common prime intervals for both adjacent rectangles. For each eligible tt, put kX=⌊Xt/H⌋k_X=\lfloor Xt/H\rfloor for X∈{λ,μ,κ}X\in\{\lambda,\mu,\kappa\}. After its exact floor-stability bound, the interval

(mmax⁡{λ+1/m0kλ+1,μkμ+1,κ+1/m0kκ+1},Hm−1t]\left( m\max\left\{\frac{\lambda+1/m_0}{k_\lambda+1}, \frac\mu{k_\mu+1}, \frac{\kappa+1/m_0}{k_\kappa+1}\right\}, \frac{Hm-1}{t}\right]

is contained in both coefficient-gcd intervals for every m≥m0m\ge m_0. The intervals for distinct tt are disjoint. The retained explicit prime bounds give a rational lower bound for their total logarithmic prime product, divided by mm.

With 4,000 terms, all twenty original logarithmic gcd rates ℓ\ell are valid numerically: fifteen sources have m0=108m_0=10^8 and five have m0=109m_0=10^9. The gcd certificate records every exact rate comparison and names the source and ramp files.

The finite bootstrap then reduces all twenty starting indices to at most 1879518795, while preserving every original rate ℓ\ell. Here is its uniform interval argument. For m∈[L,U]m\in[L,U] and a tt whose floors have stabilized by LL, use the fixed common prime interval

(max⁡{λU+1kλ+1,μUkμ+1,κU+1kκ+1},HL−1t].\left( \max\left\{\frac{\lambda U+1}{k_\lambda+1}, \frac{\mu U}{k_\mu+1}, \frac{\kappa U+1}{k_\kappa+1}\right\}, \frac{HL-1}{t}\right].

Every prime in this interval divides both coefficient gcds for every integer m∈[L,U]m\in[L,U]. Summing directed lower bounds for their logarithms and comparing with ℓU\ell U certifies the entire interval at once. An exact sieve and directed fixed-point logarithm table cover endpoints through 10810^8; the retained BFT prime bounds cover all larger endpoints, which remain below 101110^{11}.

The bootstrap works backwards from each analytic starting index. Its certified intervals form a contiguous partition, checked at every step. Stopping immediately above the first unsupported singleton leaves no gap before the analytic tail. In total, the certificate records 1 356 6731\,356\,673 certified intervals. The resulting m0m_0 values, in the ramp's binding order, are

14192,17247,6201,8976,14366,9130,13835,16349,10921,18795,8140,9636,14996,11952,16697,4571,5683,7127,13737,15235.14192,17247,6201,8976,14366,9130,13835,16349,10921,18795, 8140,9636,14996,11952,16697,4571,5683,7127,13737,15235.

These are prime-interval certificates, not a scan of binomial rows.

Numerical cofactor alternatives

Let X=pkuX=p^ku, Y=qlvY=q^lv lie in [h,h+R][h,h+R], with ∣X−Y∣≤D|X-Y|\le D. The seed sides are powers of p,qp,q, respectively. Suppose u≤hepu\le h^{e_p} and v≤heqv\le h^{e_q}. Define

ρ=min⁡{1−epL0log⁡P,1−eqL1log⁡Q}.\rho=\min\left\{\frac{1-e_p}{L_0\log P}, \frac{1-e_q}{L_1\log Q}\right\}.

Choose a positive rational t<ρt<\rho such that tK0>1tK_0>1, and set m=⌊tlog⁡h⌋m=\lfloor t\log h\rfloor. Then both residual prime-power blocks are integral. The determinant supplies a nonzero integer, and (2), the gcd bound G≥eℓmG\ge e^{\ell m}, and the exact normalizations give

1≤m(A1e−C0m+YA2e−K0m),(3)1\le m\left(A_1e^{-C_0m}+Y A_2e^{-K_0m}\right),\tag{3}

where C0,K0C_0,K_0 are the stated unequal-block rates and

A1=DTFMQ−2,A2=D03TFME−2.(4)A_1=DTF M_Q^{-2},\qquad A_2=\frac{D_0^3}{T}F M_E^{-2}. \tag{4}

The factors TT and D03/TD_0^3/T cover both values of ε\varepsilon. For h>Rh>R, put x=log⁡hx=\log h and

η=min⁡{tC0,tK0−1}>0,\eta=\min\{tC_0,tK_0-1\}>0, L=log⁡(2t)+max⁡{log⁡A1+C0,log⁡(2A2)+K0}.L=\log(2t)+\max\{\log A_1+C_0,\log(2A_2)+K_0\}.

Since m≤txm\le tx, m≥tx−1m\ge tx-1, and Y<2hY<2h, (3) implies

1≤eLxe−ηx.(5)1\le e^L x e^{-\eta x}.\tag{5}

The checker finds an integer KK for which, at x=Klog⁡10x=K\log10,

tx≥m0+1,x≥log⁡(R+1),ηx>1,L+log⁡x−ηx<0.(6)tx\ge m_0+1,\qquad x\ge\log(R+1),\qquad \eta x>1,\qquad L+\log x-\eta x<0.\tag{6}

The last expression decreases for x>1/ηx>1/\eta. Thus the contradiction holds for every h≥10Kh\ge10^K, not just at the tested endpoint.

Index twelve

For each of the twenty exact edge bindings in the ramp data, use h=n−11h=n-11 and D=R=11D=R=11. The finite-constant checker certifies (6) for every source. Its numerical certificate records all prefactors, rational block parameters, directed decay margins, and twenty individual thresholds.

Their largest exponent is 30711693071169, from the (3,5)(3,5) alternative with targets (.40,.10)(.40,.10) and gcd starting index 1383513835. The existing F8F_8 ramp has gain 1/101/10 and algebraic bound 1084510^{845}. Adding one decimal power absorbs the row shift of eleven and proves (1).

For comparison, omitting the finite bootstrap gives the weaker but still fully numerical cap 1015893357621410^{158933576214}. The source certificate records both versions.

sh
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_index_12_ramp.py
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_off_diagonal_gcd_numeric.py
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_off_diagonal_gcd_bootstrap.py
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_off_diagonal_explicit_pairs.py