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: off-diagonal Padé approximants and unequal prime-power blocks

P699 research · Joint rational-block estimates


The diagonal choice A=CA=C in BFT's application of Padé approximation is not necessary. The general polynomials in its Lemma 3.1 give two different prime-power block lengths. Retaining this difference improves some of the cofactor estimates. The conclusion below is effective; it does not assume an ineffective approximation theorem.

These are conditional source inequalities. Excluding one previously difficult cofactor vector does not establish an entire binomial index: every new alternative and every position assignment still needs coverage.

The exact identity and nonvanishing

Use the integral normalizations of BFT (3.1)--(3.4), so that

PA,B,C(z)−(1−z)B+C+1QA,B,C(z)=zA+C+1EA,B,C(z).(1)P_{A,B,C}(z)-(1-z)^{B+C+1}Q_{A,B,C}(z) =z^{A+C+1}E_{A,B,C}(z).\tag{1}

Their degrees are respectively C,A,BC,A,B. Let P′,Q′,E′P',Q',E' refer to (A+1,B−1,C+1)(A+1,B-1,C+1), where B≥1B\ge1. The exponent B+C+1B+C+1 in (1) is unchanged, and subtracting the two identities gives

PQ′−QP′=zA+C+1(EQ′−z2QE′).PQ'-QP'=z^{A+C+1}(EQ'-z^2QE').

The left side has degree at most A+C+1A+C+1. Therefore

PQ′−QP′=DA,B,CzA+C+1,DA,B,C≠0.(2)\boxed{PQ'-QP'=D_{A,B,C}z^{A+C+1},\qquad D_{A,B,C}\ne0.}\tag{2}

Indeed its coefficient of zA+C+1z^{A+C+1} is E(0)Q′(0)E(0)Q'(0), which is nonzero by the integral formulas, or by the binomial-coefficient expansions. This proves the required determinant statement for all rectangles, not only A=CA=C. The integral normalization fixes the signs consistently; only absolute values are used below.

Choose a seed aP0−bQ0=D0>0aP_0-bQ_0=D_0>0, and set z=D0/(aP0)z=D_0/(aP_0). Let GG be the gcd of the coefficients of QA,B,CQ_{A,B,C}. The identity shows that GG also divides every coefficient of PA,B,CP_{A,B,C} and EA,B,CE_{A,B,C}. Thus the integers

P=G−1(aP0)CP(z),Q=G−1(aP0)AQ(z),E=G−1(aP0)BD0A+C+1E(z)\mathcal P=G^{-1}(aP_0)^C P(z),\quad \mathcal Q=G^{-1}(aP_0)^A Q(z),\quad \mathcal E=G^{-1}(aP_0)^B D_0^{A+C+1}E(z)

satisfy

(aP0)A+B+1P−(bQ0)B+C+1Q=E.(3)\boxed{(aP_0)^{A+B+1}\mathcal P -(bQ_0)^{B+C+1}\mathcal Q=\mathcal E.}\tag{3}

Both block lengths A+B+1A+B+1 and B+C+1B+C+1 are preserved by the adjacent rectangle in (2). Hence the two choices still supply a nonzero integer linear combination after any residual denominators are cleared.

Using a common denominator power max⁡(A,C)\max(A,C) instead would lose this advantage. Formula (3) is the reason unequal rectangles help.

An effective coefficient-gcd lower bound

Fix positive integer directions λ,μ,κ\lambda,\mu,\kappa, and take

A=λm,B=μm−1,C=κm,H=λ+μ+κ.A=\lambda m,\qquad B=\mu m-1,\qquad C=\kappa m, \qquad H=\lambda+\mu+\kappa.

For a finite set of positive integers tt satisfying

⌊λtH⌋+⌊μtH⌋+⌊κtH⌋=t−2,\left\lfloor\frac{\lambda t}{H}\right\rfloor+ \left\lfloor\frac{\mu t}{H}\right\rfloor+ \left\lfloor\frac{\kappa t}{H}\right\rfloor=t-2,

put

LT=∑t[Ht−max⁡{λ⌊λt/H⌋+1,μ⌊μt/H⌋+1,κ⌊κt/H⌋+1}].(4)L_T=\sum_t\left[\frac Ht- \max\left\{ \frac\lambda{\lfloor\lambda t/H\rfloor+1}, \frac\mu{\lfloor\mu t/H\rfloor+1}, \frac\kappa{\lfloor\kappa t/H\rfloor+1} \right\}\right].\tag{4}

This is an explicit positive rational number. Every summand comes from a disjoint prime interval in BFT Lemma 5.2. For any rational ℓ<LT\ell<L_T,

G≥eℓm(5)G\ge e^{\ell m}\tag{5}

holds above an effectively computable m0m_0, for both adjacent rectangles.

Here is why the bounded shifts in A,B,CA,B,C cause no difficulty. For an eligible tt, the three fractional parts in (4) sum to 22. None can be zero, since the other two are strictly less than 22 in total. Consequently all the floors stabilize effectively under these shifts. The corresponding interval endpoints are mm times the displayed constants, up to bounded errors. The effective prime-number estimates already retained in the Dusart note give θ(x)/x⟶1\theta(x)/x\longrightarrow1 with a computable rate. Apply them to the finitely many interval endpoints and absorb the strict difference LT−ℓL_T-\ell. This proves (5) with an effective threshold. No estimate for an infinite tail of (4) is needed.

The two exponential rates

Define

S=Hlog⁡H−λlog⁡λ−μlog⁡μ−κlog⁡κ,S=H\log H-\lambda\log\lambda-\mu\log\mu-\kappa\log\kappa, q=S+max⁡0<u<1{μlog⁡u+κlog⁡(1−u)+λlog⁡(1−(1−z)u)},q=S+\max_{0<u<1} \{\mu\log u+\kappa\log(1-u) +\lambda\log(1-(1-z)u)\}, e=S+max⁡0<u<1{λlog⁡u+κlog⁡(1−u)+μlog⁡(1−zu)}.e=S+\max_{0<u<1} \{\lambda\log u+\kappa\log(1-u) +\mu\log(1-zu)\}.

The two maxima are at unique points in (0,1)(0,1), determined by quadratic equations. They can therefore be bounded with directed intervals. The integral formulas and effective factorial estimates give the growth rates q,eq,e for Q,EQ,E, with at most a polynomial factor in mm. Such a factor is bounded effectively by eϵme^{\epsilon m} for every fixed ϵ>0\epsilon>0.

After division by the coefficient gcd, the decay rates needed in the integer argument are

C0=μlog⁡P0−λlog⁡a−(μ+κ)log⁡b−q+ℓ,C_0=\mu\log P_0-\lambda\log a-(\mu+\kappa)\log b-q+\ell, K0=λlog⁡P0+(μ+κ)log⁡Q0−μlog⁡a−(λ+κ)log⁡D0−e+ℓ.(6)K_0=\lambda\log P_0+(\mu+\kappa)\log Q_0 -\mu\log a-(\lambda+\kappa)\log D_0-e+\ell.\tag{6}

The block lengths per unit mm are L0=λ+μL_0=\lambda+\mu and L1=μ+κL_1=\mu+\kappa.

The effective joint criterion

Write nearby factors as fp=pkpup∈[h,h+R]f_p=p^{k_p}u_p\in[h,h+R], for distinct primes, and form products X,YX,Y of equal degree dd with positive multiplicities bpb_p. Assume ∣X−Y∣≤D|X-Y|\le D, with a fixed positive bound DD. The primes occurring in each product must be precisely those in the corresponding seed side

P0=∏p∈I0pap,Q0=∏p∈I1pap.P_0=\prod_{p\in I_0}p^{a_p},\qquad Q_0=\prod_{p\in I_1}p^{a_p}.

For a proposed upper cofactor vector epe_p and a real t>0t>0, set

Δ(t)=∑p(tLs(p)aplog⁡p−bp(1−ep))+,\Delta(t)=\sum_p \bigl(tL_{s(p)}a_p\log p-b_p(1-e_p)\bigr)_+,

where s(p)s(p) is the seed side. If

tC0−Δ(t)>0,tK0−Δ(t)−d>0,(7)\boxed{tC_0-\Delta(t)>0,\qquad tK_0-\Delta(t)-d>0,}\tag{7}

then the simultaneous bounds up≤hepu_p\le h^{e_p} are impossible above an effectively computable height.

To prove this, use m=⌊tlog⁡h⌋m=\lfloor t\log h\rfloor and clear the residual denominators of X/P0L0mX/P_0^{L_0m} and Y/Q0L1mY/Q_0^{L_1m}. Their logarithmic cost is at most Δ(t)log⁡h\Delta(t)\log h. Equations (2), (3), and (6) yield, for a computable constant CϵC_\epsilon,

1≤Cϵeϵm(elog⁡Z−mC0+Yelog⁡Z−mK0).1\le C_\epsilon e^{\epsilon m} \left(e^{\log Z-mC_0}+Y e^{\log Z-mK_0}\right).

Since Y≤2dhdY\le2^d h^d for h>Rh>R, both terms tend effectively to zero when ϵ>0\epsilon>0 is chosen small enough relative to the strict margins in (7). This contradiction proves the criterion.

The actual product difference may be zero. Nonvanishing of the integer linear combination comes from the seed determinant (2), whose seed gap D0D_0 is positive; it does not come from X−YX-Y. If X=YX=Y, the term involving X−YX-Y simply vanishes. Thus the same joint constraints apply when two selected prime powers belong to the same numerator factor. Their valuation and factorization hypotheses must still be checked.

The integral-block alternative is an immediate special case. Provided C0>0C_0>0, every target strictly below

1−dLs(p)aplog⁡pbpK0(8)1-\frac{dL_{s(p)}a_p\log p}{b_pK_0}\tag{8}

is a valid coordinate in a cofactor disjunction. The joint criterion (7) can be stronger than these coordinatewise limits.

The inequalities can also be replaced by rational linear alternatives for exact subsequent computations. Fix rational t>0t>0, choose rational gp>tLs(p)aplog⁡p−bpg_p>tL_{s(p)}a_p\log p-b_p, and choose a positive rational R<min⁡(tC0,tK0−d)R<\min(tC_0,tK_0-d). Above an effective threshold, necessarily

∑p(bpep+gp)+≥R.\sum_p(b_pe_p+g_p)_+\ge R.

Otherwise both margins in (7) have a uniform positive lower bound. The last inequality is equivalent to the finite disjunction

∑p∈Jbpep≥R−∑p∈Jgpfor at least one nonempty subset J.\sum_{p\in J}b_pe_p\ge R-\sum_{p\in J}g_p \quad\text{for at least one nonempty subset }J.

The rational-clause checker rounds the intercepts upwards and the threshold downwards, then verifies every rounding and a strict source margin using directed intervals. Its example certificate contains only rational coefficients for the resulting linear alternatives. This separates the analytic source verification from exact polyhedral coverage calculations.

A four-prime example

For P0=847=7⋅112P_0=847=7\cdot11^2, Q0=845=5⋅132Q_0=845=5\cdot13^2, a=b=1a=b=1, use directions (λ,μ,κ)=(39,7,33)(\lambda,\mu,\kappa)=(39,7,33). The finite sum (4) through t=200t=200 exceeds 3.04213.0421, so take ℓ=3\ell=3. Directed evaluation gives

C0/40>.01316,K0/40>11.5489.C_0/40>.01316,\qquad K_0/40>11.5489.

Thus, for four distinct positions satisfying r7+r11=r5+r13r_7+r_{11}=r_5+r_{13}, the following alternative holds effectively:

u7>h.612oru11>h.044oru5>h.721oru13>h.111.u_7>h^{.612}\quad\text{or}\quad u_{11}>h^{.044} \quad\text{or}\quad u_5>h^{.721}\quad\text{or}\quad u_{13}>h^{.111}.

The product gap is a fixed nonzero difference of the position products. This strengthens the diagonal source for the same seed. Applying it in a whole-index proof still requires covering every branch of the alternative.

The checker and certificate record the finite rational gcd sum, the two saddle points, directed bounds for both rates, and a strict margin for each listed exponent. The same checker exports the general growth calculation for other seeds and rectangles.

sh
.venv/bin/python wiki/research/erdos_699/computations/erdos_699_off_diagonal_pade.py