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: adjusting the Padé block while retaining the cofactor clause

P699 research · Numerical joint constants · Current seven and nine bounds


The numerical argument need not choose the Padé index as m=⌊tlog⁡h⌋m=\lfloor t\log h\rfloor, even when tt was used to define the retained cofactor clause. Choosing a smaller index decreases the residual denominator cost. Quantifying this decrease strengthens the numerical estimate without changing the clause or its polynomial proof forest.

The denominator saving

Put x=log⁡hx=\log h. Write a retained clause as

∑p(gp+bpep)+≥R,R>0,\sum_p(g_p+b_pe_p)_+\ge R,\qquad R>0,

where its outward-rounded intercepts satisfy

gp>tκp−bp,κp=Ls(p)aplog⁡p>0.g_p>t\kappa_p-b_p,\qquad \kappa_p=L_{s(p)}a_p\log p>0.

Here apa_p is the exponent of pp in its seed product, bpb_p is its factor multiplicity, and L0=A+BL_0=A+B, L1=B+CL_1=B+C are the two directional block lengths. Choose any positive rational

κ≤min⁡pκp.\kappa\le\min_p\kappa_p.

At integer Padé index m≤txm\le tx, the actual logarithmic residual denominator cost is bounded by

Dm(x)=∑p(mκp−bp(1−ep)x)+.D_m(x)=\sum_p\bigl(m\kappa_p-b_p(1-e_p)x\bigr)_+.

Lemma. If the retained clause fails, then

Dm(x)≤(Rx−κ(tx−m))+.(1)\boxed{D_m(x)\le\bigl(Rx-\kappa(tx-m)\bigr)_+.}\tag{1}

Indeed, set zp=(gp+bpep)xz_p=(g_p+b_pe_p)x. The expression inside each actual positive part is at most zp−κp(tx−m)z_p-\kappa_p(tx-m). For arbitrary real zpz_p and sp≥s≥0s_p\ge s\ge0,

∑p(zp−sp)+≤(∑p(zp)+−s)+.\sum_p(z_p-s_p)_+\le\left(\sum_p(z_p)_+-s\right)_+.

If no summand survives the subtraction, the assertion is immediate. Otherwise, summing over surviving indices subtracts at least ss and uses at most the original sum of positive parts. Apply this with s=κ(tx−m)s=\kappa(tx-m) and ∑p(zp)+<Rx\sum_p(z_p)_+<Rx to obtain (1).

The subtraction in (1) is made once, not once for every prime. No assumption about the active sign mask is needed.

A fixed index covering a height interval

Use the finite constants A1,A2A_1,A_2 from the linked numerical argument, and let dd be the product degree. Write c,kc,k for the two exponential decay rates before adding a coefficient-gcd contribution. If a common divisor of both adjacent Padé rectangles has logarithm at least GmG_m, the nonzero-integer argument gives

1≤mA1exp⁡(Dm(x)−cm−Gm)+2dmA2exp⁡(dx+Dm(x)−km−Gm).(2)1\le mA_1\exp\bigl(D_m(x)-cm-G_m\bigr) +2^dmA_2\exp\bigl(dx+D_m(x)-km-G_m\bigr).\tag{2}

Let y=tx≥my=tx\ge m. The upper bound in (1) becomes

Dm(x)≤max⁡{0,κm+(R/t−κ)y}.(3)D_m(x)\le\max\{0,\kappa m+(R/t-\kappa)y\}.\tag{3}

After taking logarithms, twice the first term in (2) is bounded by the maximum of two affine functions of yy. Twice the second term is bounded by the maximum of two affine functions with slopes

d/t,(R+d)/t−κ.d/t,\qquad (R+d)/t-\kappa.

These functions are convex. Checking strict negativity at both endpoints therefore certifies every height between them. If κt≥R\kappa t\ge R, the first term is nonincreasing with yy; only its left endpoint needs checking.

The left endpoint may exceed mm. Thus, at an already covered height x0x_0, one can choose m<tx0m<tx_0, check (2) there, and use the same index until the right endpoint fails. A sequence of overlapping certified intervals then joins a finite height to an existing uniform analytic tail. The usual condition h≥s+1h\ge s+1 ensures that the shifted factors are at most 2h2h; it is checked separately from the interval inequalities.

A smaller block for a uniform tail

Choose a rational 0<σ≤t0<\sigma\le t and set m=⌊σx⌋m=\lfloor\sigma x\rfloor. Formula (1) implies

Dm(x)≤Bσx,Bσ=max⁡{0,R−κ(t−σ)}.(4)D_m(x)\le B_\sigma x,\qquad B_\sigma=\max\{0,R-\kappa(t-\sigma)\}.\tag{4}

Suppose log⁡Gm≥ℓm\log G_m\ge\ell m for every m≥m0m\ge m_0. Put C=c+ℓC=c+\ell and K=k+ℓK=k+\ell. Sufficient positive decay margins are now

σC>Bσ,σK>Bσ+d.(5)\boxed{\sigma C>B_\sigma,\qquad \sigma K>B_\sigma+d.}\tag{5}

All floor and integral prefactors are handled exactly as in the original numerical argument, with block σ\sigma and denominator exponent BσB_\sigma. In particular σx≥m0+1\sigma x\ge m_0+1 ensures that the chosen index lies in the uniform gcd range. The two inequalities force C,K>0C,K>0, so replacing mm by σx−1\sigma x-1 introduces the valid factors eC,eKe^C,e^K.

When K>CK>C, the value

σ∗=dK−C\sigma_* = \frac{d}{K-C}

balances the two margins. It is a useful candidate when σ∗≤t\sigma_*\le t; a nearby rational value is checked with directed bounds. This is a choice to test, not an assertion that it is optimal in every parameter regime. The difference K−C=k−cK-C=k-c is independent of the chosen gcd lower rate.

Every whole-index bound must certify every source and join the resulting height to the polynomial and finite bounds. The complete proof at seven supplies these certificates. The complete proof at nine now does so as well.

Initial sources with explicit cofactor targets

For an initial alternative, its negation directly bounds each ep≤αp<1e_p\le\alpha_p<1. One can then use

Dm(x)≤∑p(mκp−bp(1−αp)x)+.(6)D_m(x)\le\sum_p\bigl(m\kappa_p-b_p(1-\alpha_p)x\bigr)_+.\tag{6}

This bound does not require m≤txm\le tx: the nominal block used to construct a previous clause is irrelevant. In numerical upper bounds for (6), each κp\kappa_p is rounded upward.

The right side is the maximum, over all subsets SS of the component primes, of

m∑p∈Sκp−x∑p∈Sbp(1−αp).m\sum_{p\in S}\kappa_p-x\sum_{p\in S}b_p(1-\alpha_p).

Every slope is nonpositive, so the first term of (2) need only be checked at the left endpoint. After adding dxdx, the second term is again the maximum of finitely many affine functions. Its endpoint checks certify the entire interval. For a pair there are only four subsets. The balanced choice m≈dx0/(K−C)m\approx d x_0/(K-C) can be tested even when it exceeds the old nominal value tx0tx_0.