Plectis

Long paper

Factorial Linear Forms and Denominators: Detailed Proofs and Rationality Criteria

Erdős #68 By 1 September 2026; revised 7 October 2026 Authorship, AI use and citation

Summary. The record proves liminf (log LN)/(N3/2 log N) ≥22/3 for L_N=lcm(2!−1,…,N!−1), by comparing a terminal block of denominators with pairwise gcds rather than an external multiplicity theorem. Specialising a Louwsma–Martino valuation, it exhibits primes 139 and 2593 absent from all later reduced denominators. Coefficient constructions give integer forms MS+k; at fixed moment the remainder changes only by an integer. Rationality criteria use factorial digits and the next integer above a scaled partial sum; the digit argument also recovers the irrationality of e. Two exact calculations give q∤299999! and q≥239990>1012038. These exclusions do not prove irrationality.

This paper covers the complete problem-specific reasoning surface for Erdős #68, including all registered result families and their boundaries.

It does not establish the validity of claims tagged Lean, which belongs to the cited kernel-checked source, or a solution to Erdős #68, which remains open.

In this paper

Problem 0.1 (Erdős #68). Is S=∑n≥2(n!−1)−1 irrational?

Erdős posed this question in 1988, together with the expectation that the corresponding series with denominators n!+t should be transcendental for every integer t [1]. For negative t, the sum begins after all zero or negative denominators. Any further change of starting index adds a rational number. The expectation for all integer shifts remains unproved here. The problem number follows Bloom’s catalogue [2].

For finitely supported integer vectors on n≥2, cancellation forms MS+k give an explicitly computable ideal of multipliers. Correcting adjacent factorial differences at proper divisors gives an integral basis. Cancellation fixes its lower coordinates; excluding index one leaves a single Diophantine equation. Its gcd, computed below 2D2, gives the least positive multiplier and an attaining vector. At depth four, the support restriction strengthens 115∣M to 1380∣M. This classification is independent of nonintegrality.

The accompanying short paper proves this classification and bounds one resulting remainder between consecutive integers. Here we develop the support and coefficient-norm refinements and a primitive construction on arithmetic progressions. These solve different problems: at fixed M, changing the vector changes MS+k only by an integer, so it cannot alter nonintegrality. The progression construction makes M divisible by every fixed denominator eventually, but still needs a gap estimate for those same forms.

The later sections study that arithmetic difficulty from two directions. Common-denominator growth measures the cost of clearing terms; prime-power calculations determine what survives reduction. Carries and factorial digits instead characterize rationality itself. The resulting equivalences identify what an irrationality proof must show, without proving that the required departures occur arbitrarily far out.

Throughout, write

dn=n!−1,HN=∑n=2N1dn,LN=lcm⁡(d2,…,dN),AN=∑n=2NLNdn.

Thus HN=AN/LN. We write den⁡(x) for the positive reduced denominator of a rational number x. Under S=a/q, every form MS+k with M,k∈Z and q∣M is integral. For each positive integer q, an irrationality argument of this kind must produce a nonintegral form with q∣M. Both properties must hold for the same vector; cancellation and growth of M alone do not establish this.

The proof of the classification is in §1.1; the following table locates the additional arguments. In particular, the lower bound for LN concerns the unreduced common denominator, not automatically den⁡(HN). The appendices retain the supplementary identities, limits of the attempted arguments, and formal-source concordance.

Question

Where it is treated
Which M and coefficient vectors are possible? §1.1
What do support and norm minimisation change? §1.2
Which prime powers survive reduction? §2
How fast does the common denominator grow? §3
How does growth constrain the support parameter? §3.1
Which eventual patterns characterize rationality? §4
Which denominators are excluded by finite computations? §6
Where are the formal declarations? Appendix A.1

Integer coefficient vectors and cancellation

We need weights that preserve the fractional part of i!/(d!−1) and whose scaled adjacent differences have a simple form. Both requirements are met by removing one factor d! for each complete block of d indices. For a finitely supported integer vector c=(ci) on i≥2, put

M(c)=∑icii!,Wd,i=i!(d!)⌊i/d⌋,Vd(c)=∑iciWd,i(d≥2).

We allow either sign for M(c) and for Vd(c). The associated remainder is

R(c)=∑d≥2Vd(c)d!−1.

The index i selects a coefficient, whereas d selects a summand of the remainder. If c is supported on i≤N, then for d>N all floor exponents vanish and Vd(c)=M(c). The series therefore converges absolutely. Cancellation of the first weighted sums leaves a finite signed sum followed by M(c) times the original tail. We show below that the remainder differs from M(c)S by an integer.

The chosen exponent is a floor, rather than the largest possible power of d! dividing i!. For example, at d=2, i=6, it is 3 although 24∣6!. What matters is that it changes between i−1 and i exactly when d∣i. This makes adjacent differences amenable to divisor elimination in §1.1. Each Wd,i is an integer: writing i=kd+r with 0≤r<d, the quotient i!/((d!)kr!) is a multinomial coefficient. Since i!=(d!)⌊i/d⌋Wd,i and d!≡1(modd!−1),

(1)Vd(c)≡M(c)(modd!−1).

Theorem 1.1 (divisibility of the difference). For every finite integer support and every d≥2 there is an integer k with Vd(c)=M(c)+(d!−1)k.

Proof. By (1), d!−1 divides Vd(c)−M(c). Thus k=(Vd(c)−M(c))/(d!−1) is an integer and gives the identity. The modulus is positive also at d=2, when it equals one. ◻

For any N≥2 at least as large as every supported index, the terms with d>N equal M(c)/(d!−1). Therefore

R(c)−M(c)S=∑d=2NVd(c)−M(c)d!−1∈Z.

The sum is finite and each term is integral by (1). This argument also holds when index 1 is temporarily allowed in the basis calculation below. In particular, M=0 gives an integral remainder, and fixing M(c) fixes its fractional part. These conclusions hold for every finite vector, before imposing cancellation.

By (1), a vanishing d-th weighted sum forces (d!−1)∣M(c), and annihilating every weighted sum 2≤d≤D forces LD∣M(c), where LD=lcm2≤d≤D⁡(d!−1) is the same quantity as in §3. With M fixed, the same congruence fixes the fractional part of each quotient Vd(c)/(d!−1) separately.

An integral basis and the possible values of M

Fix D≥2. We temporarily allow a coefficient at index one so that the adjacent-difference vectors form a triangular basis. After solving V2=⋯=VD=0, we impose one further equation to remove that coordinate. Only the coefficient space is enlarged; the series S continues to start at d=2, with no term 1/(1!−1).

Let en be the unit vector at index n. For n≥2 put

Tn=nen−1−en.

We have M(Tn)=0. If d∤n, the two floor exponents in its weighted sum agree, and nWd,n−1=Wd,n. If d∣n, the exponent increases by one, so nWd,n−1=d!Wd,n. Hence

Vd(Tn)=(d!−1)Wd,n1d∣n.

For example, T4=4e3−e4 has weighted sums V2=6 and V4=23, and all its other weighted sums vanish. Since T2=2e1−e2 has only V2=1 nonzero, T4−6T2 isolates the fourth weighted sum.

We therefore define, in increasing order of n,

Un=Tn−∑d∣n2≤d<nWd,nUd.

Assume the construction has isolated the weighted sums below n. Subtracting Wd,nUd removes the contribution at a proper divisor d and changes no other weighted sum. At d=n we have Wn,n=1. Induction therefore gives M(Un)=0 and Vd(Un)=(d!−1)1d=n for every d≥2.

For the basis claim, restrict to indices 1,…,N. The matrix with columns e1,T2,…,TN is triangular over Z with diagonal 1,−1,…,−1. The change from Tn to Un is triangular with diagonal one. Both determinants are units in Z, so both inverses have integer entries. In back-substitution the only divisors are 1 and −1; no rational coefficients are needed. Thus e1,U2,…,UN represents every integer vector uniquely. If this expansion is be1+∑zdUd, then M(c)=b and Vd(c)=b+(d!−1)zd. These identities determine the coordinates:

c=M(c)e1+∑d≥2Vd(c)−M(c)d!−1Ud.

The congruence (1) also verifies the integrality of the displayed coefficients. They vanish for d>N, since then Vd(c)=M(c), so the expansion is finite. These same coefficients occur in the integer-difference identity at the start of §1: their sum is R(c)−M(c)S. Thus the coordinates also recover the integer term in the linear form.

Cancellation through D forces LD∣M(c). To obtain M(c)=LD in the enlarged space, start with LDe1, whose weighted sums all equal LD. The vector Ud changes only the dth sum, by d!−1. Hence define

KD=LDe1−∑d=2DLDd!−1Ud.

This vector satisfies M=LD and V2=⋯=VD=0. All solutions in the enlarged coefficient space consequently have the unique form

(2)c=tKD+∑n>DznUn,M(c)=tLD,

where t,zn are integers and only finitely many zn are nonzero. The only forbidden coefficient is now the one at index 1. Write aD=(KD)1 and un=(Un)1. Removing this coefficient is equivalent to

(3)taD+∑n>Dznun=0.

The free higher basis coefficients can change the first coefficient by every multiple of gcd{un:n>D}, and only by such a multiple. Thus the gcd gives the exact solvability condition. We next show how to compute it using a bounded list of coefficients. Here u2=2, whereas for n>2 the recurrence reads

un=−∑d∣n2≤d<nWd,nud.

It gives un=0 for odd n>1. If p is prime, the only proper divisor of 2p with a nonzero contribution is 2, and therefore u2p=−(2p)!/2p−1, including p=2.

Choose a prime p with D/2<p≤D. Such a prime exists by Bertrand’s postulate when D≥3, and p=2 works when D=2. The coefficient u2p lies in the tail and divides (2p)!. Once a finite initial interval includes that coefficient, its gcd therefore divides this factorial. To control the small-divisor weights beyond the interval, we shall require n/d≥2p whenever d≤D. This leads to the cutoff H=D(2p−1). We claim that

(4)gcd{un:n>D}=gcd(uD+1,…,uH)>0,H<2D2.

Let g be the finite gcd on the right. Its value could in principle decrease if we included later coefficients. To rule this out, we prove g∣un for every n>D by strong induction. The index 2p lies in the finite interval and has a nonzero coefficient, so g>0 and g∣(2p)!. The definition proves the induction claim through H.

For n>H, each term of the recurrence is Wd,nud with d<n. If d>D, divisibility comes from ud by induction. If d≤D, we cannot use that hypothesis and instead use the weight. Put k=n/d. The cutoff was chosen so that n>D(2p−1) forces k≥2p for every such divisor. Also,

Wd,nk!=n!(d!)kk!∈Z,

because the quotient counts partitions of n labelled objects into k unordered blocks of size d. Hence g∣(2p)!∣k!∣Wd,n. Thus every product in the recurrence is divisible by g, completing the induction. The finite gcd divides the entire tail, and inclusion of the finite interval gives the reverse divisibility. Finally, H≤D(2D−1)<2D2.

Theorem 1.2 (the set of attainable values of M). Fix D≥2 and a prime p with D/2<p≤D. Put

H=D(2p−1),GD=gcd{|un|:D<n≤H},μD=LDGDgcd(GD,aD).

For finite integer vectors supported on n≥2 and cancelling all weighted sums 2,…,D, the possible values of M are exactly μDZ. The integer μD is positive, is independent of the eligible prime p, and is attained by a vector of coefficients with gcd one, that is, by a primitive vector.

Proof. Equation (4) identifies GD with the gcd of all un for n>D. Their finite integer combinations are exactly GDZ. Hence (3) has a solution if and only if GD∣taD, or equivalently GD/gcd(GD,aD)∣t. Substitution in (2) gives the claimed set of values of M. Bézout’s identity applied to the finite list through H supplies a vector attaining its least positive element. A common coefficient divisor larger than one could be divided out, giving a smaller positive attainable value of M. The attaining vector is therefore primitive. Finally, the set of values of M is intrinsic to the cancellation equations, so its positive generator is independent of the eligible prime p. ◻

Fix M∈μDZ and put t=M/LD. The support equation (3) determines which basis coefficients are admissible. The following formula records their integer contribution and shows why changing the vector at that value of M cannot alter nonintegrality.

Theorem 1.3 (how the coefficient choices change the remainder). For the vector in (2),

R(tKD+∑n>DznUn)=tLD(S−HD)+∑n>Dzn.

The remainder series converges for every finite vector supported away from index zero. A vector with M=0 has integral remainder, and any two finite vectors with the same value of M have remainders differing by an integer.

Proof. For d>D, we have Vd(KD)=LD and Vd(Un)/(d!−1)=1d=n. Summing these identities gives tLD(S−HD)+∑n>Dzn, as claimed. The first series converges and the second sum is finite.

Absolute convergence and the general integer-difference identity were proved at the start of §1, including for the auxiliary index 1. Taking M=0 proves integrality; subtracting the identities for two vectors with equal M proves the last assertion. ◻

Here M=tLD, so every denominator in HD divides M. Thus MHD∈Z, and the integer in the original form MS+k is k=−MHD+∑n>Dzn.

Examples and support refinements

For D=4 we need a common divisor of the entire coefficient tail, as well as a finite combination attaining it. Put Λn=lcm⁡(1,…,n). We first prove Λn∣un. For d∣n and a prime r, Legendre’s formula gives

vr(Wd,n)=∑j≥1(⌊nrj⌋−nd⌊drj⌋).

Every summand is nonnegative. If d<rj≤n, the second floor is zero and the first is at least one. Counting these powers gives vr(Wd,n)≥vr(Λn)−vr(Λd), and therefore Λn∣Wd,nΛd. Starting from u2=2, induction in the recurrence proves Λn∣un.

Thus 60∣un for all n>4. Conversely, the finite calculation u6=−180, u8=−4200 gives 23u6−u8=60. The gcd is consequently G4=60. Since L4=115 and a4=−55, we obtain

μ4=11560gcd(60,55)=1380.

The vector 12K4+253U6−11U8 attains this value: its coefficient at index 1 is 12(−55)+253(−180)−11(−4200)=0.

The five-term vector used in the short paper is obtained from the same coordinates. The equality 12(−55)−27(−180)−4200=0 first gives 12K4−27U6+U8. Its coefficients at indices four and seven are removed by −26U5+8U7, respectively. Both corrections leave the first coordinate zero. Consequently

(5)12K4−26U5−27U6+8U7+U8=1482e2−784e3−136e5+83e6−e8.

The coefficients of the Un in this expansion sum to −44. Since H4=143/115, its remainder is 1380(S−H4)−44=1380S−1760. The exact interval calculation appears in §B.6. Thus the coordinates give both the original linear form and the cancellations, independently of a search over supports.

If support is additionally restricted to n≤6, only u5=0 and u6=−180 remain available in (3). The least positive value of M is then 115⋅180/gcd(180,55)=4140. Thus a restriction on the largest supported index can increase the least attainable value of M. The following refinements specify the support bound in each case.

There is also a direct congruence proof using only the second and fourth weighted sums. For every n≥2,

(6)1380∣11n!−46W2,n+12W4,n.

The expression is zero for n=2,3,4. For n≥4, it is divisible by 12: 12∣n! and 6∣W2,n, as follows from W2,2r=r!∏j=1r(2j−1) for r≥2 and the corresponding odd-index recurrence. For n≥5, it is divisible by 5 because 5∣n! and the denominators removed to form W2,n and W4,n are coprime to 5. Modulo 23, n!≡W4,n and 46W2,n≡0. This proves the congruence. Consequently V2=V4=0 forces 1380∣M, without a restriction on the largest supported index. Conversely, 1482e2−784e3−136e5+83e6−e8 satisfies M=1380 and V2=V3=V4=0. Thus imposing the additional equation V3=0 does not change the possible values of M, although 3e2−e3 shows that it changes the set of admissible vectors.

Among integer vectors supported in {2,…,8} with M=1380 and V2=V3=V4=0, this vector has the smallest support. Every such vector has the following coefficients for some t,x,y∈Z. The coefficient of K4 remains M/L4=12; here t is a new free parameter.

(7)λ2=1482−4200t,λ3=−784+2240t,λ4=−680+1890t−5x,λ5=x,λ6=83−238t−7y,λ7=y,λ8=−1+3t.

Indeed, the values of the expression in (6) at indices 2,…,8 are

(0,0,0,0,4140,28980,328440).

Therefore the dual identity and M=1380 give 11=3λ6+21λ7+238λ8. Modulo 3 this yields λ8=−1+3t, then λ6+7λ7=83−238t. Set x=λ5 and y=λ7; solving V2=V3=V4=0 for the coefficients at indices 2,3,4 gives (7). The coefficient matrix is (136214261), with determinant 55, and substitution verifies the converse. Indices 2,3,8 are nonzero for every integral t; neither the pair (λ4,λ5) nor (λ6,λ7) can vanish together. Thus five nonzero entries are necessary at this support bound.

The smallest coefficient ℓ1 norm for support in {2,…,8} is 2416, attained by (λ2,…,λ8)=(1482,−784,0,−136,−1,12,−1). For t≥1, the absolute values at indices 2,3,8 alone sum to 6443t−2267≥4176; for t≤−1 they sum to 2267−6443t≥8710. Thus a minimizer has t=0. The remaining independent minima are minx(|−680−5x|+|x|)=136 at x=−136 and miny(|83−7y|+|y|)=13 at y=12. For the first, 680≤|−680−5x|+5|x| is at most five times the objective. For the second, the objective is at least 83 when y≤0, at least 17 when 0≤y≤11, and equals 8y−83≥13 when y≥12. Adding the fixed contribution 1482+784+1=2267 gives 2416. Equality forces t=0, x=−136 and y=12, so this norm minimizer is unique. Among the five-term vectors, t=0 is again necessary for the least norm; the pair (λ6,λ7) can then have one nonzero entry only when y=0, since 7∤83. The pair (λ4,λ5) has one nonzero entry only at x=0 or x=−136, and the latter is smaller. Hence the displayed five-term vector is the unique least-norm vector among those of minimum support, with norm 2486.

The additional equation V3=0 changes the minimum support size and coefficient norm, while leaving the possible values of M unchanged. With only V2=V4=0, group the seven coefficients by a=λ2+3λ3, b=λ4+5λ5, c=λ6+7λ7 and d=λ8. Solving V2=V4=0 and M=1380 gives

(a,b,c,d)=(−870+2520t,−680+1890t,83−238t,−1+3t),t∈Z.

Each coordinate is nonzero for integral t, so four nonzero coefficients are necessary and sufficient. One four-term witness is −290e3−136e5+83e6−e8, with norm 510. For the coefficient norm on the same set of indices, |λ2|+|λ3|≥|a|/3 and |λ4|+|λ5|≥|b|/5; when t≠0, these bounds together with |d| already sum to at least 794. At t=0, the four groups have unique norm minima 290, 136, 13 and 1, respectively. Thus the unique norm minimizer under only V2=V4=0 is −290e3−136e5−e6+12e7−e8, with norm 440. For comparison, with M=1380 and support contained in {2,…,8}, the two optimization problems give:

Cancellation

Minimum support size Minimum ℓ1 norm
V2=V4=0 4 440
V2=V3=V4=0 5 2416

In both rows the norm minimizer has one more nonzero coefficient than a support minimizer. Thus even with M and the allowed indices fixed, least support and least norm need not coincide. These are finite optimization statements, not irrationality conclusions.

There is a uniform restriction 12LD∣M on support n≥2. For n=2,3 one has n!=2W2,n. For n≥4, both n! and 2W2,n are divisible by 12: if n=2k or 2k+1, then W2,2k=k!∏j=1k(2j−1) is divisible by 6 for k≥2. Thus 12∣M(c)−2V2(c). Every d!−1 is coprime to 6, so V2=⋯=VD=0 implies 12LD∣M(c). This factor is attained at D=2 and D=3 by −6e2+e4 and −6e2−8e3+5e4, respectively. The next example shows that this need not be the whole restriction.

At D=6, the same recurrence gives

L6=9839515,a6=−2242555,G6=840.

To verify the last value, u7=0, while 840=lcm⁡(1,…,8) divides every un for n≥8 by the preceding divisibility. Conversely, u8=−4200, u12=−14386680 and gcd(u8,u12)=840. Thus gcd(G6,a6)=35 and the least positive value of M is μ6=24L6, not 12L6. It is attained by 24K6−16240U8+U12: the coefficient at index 1 is zero and the coefficient at 12 is −1. The finite-gcd formula thus detects the extra factor at depth six, which the uniform divisibility by 12LD misses.

Support bounds and nonintegrality

There are three separate choices: the vector, its multiplier M, and the cutoff used to estimate its remainder. The basis and support equations classify the first two. Once M is fixed, R(c)≡M(c)S(modZ) fixes the fractional part: integer translation may put a representative near zero but cannot change its distance from an integer.

To exclude every rational denominator, each q>0 must divide the multiplier of some nonintegral form. Growth of the compulsory divisor 12LD alone is insufficient: its 3-adic valuation is exactly one for every D, so it does not even force 9∣M. By contrast, extending the cutoff changes neither M nor the remainder, only the finite enclosure used to test it.

For c=−6e2+e4, we have M(c)=12 and R(c)=12S−17. At cutoff 4, its finite part is −239/115: the gap 9/115 is smaller than the tail bound 24/119. At cutoff 5, the finite part is −27061/13685, whose gap 13376/13685 exceeds the new bound 24/719. The finite part has crossed −2, so the next integer is now −1 rather than −2. This explains the larger gap. The second cutoff certifies nonintegrality of the same remainder; failure at the first cutoff was inconclusive. The successful test excludes denominators dividing the unchanged value M=12, already covered by the finite exclusions.

For fixed c with M(c)>0, the finite part increases to R(c) with omitted tail M(c)(S−HN). If R(c)∉Z, its gap to the next integer is positive, whereas 2M(c)/((N+1)!−1) tends to zero; the test therefore succeeds at every sufficiently large cutoff. If R(c)∈Z, the eventual gap equals the omitted tail and is smaller than this upper bound. Eventual success is consequently equivalent to nonintegrality of this fixed form, not an independent reason for it. The short paper’s remainder comparison gives the inequality explicitly. For each q>0 we still need a vector with q∣M(c) for which it succeeds; the successful cutoff need not be the largest supported index.

The support itself has a separate geometric restriction. Cancellation with M≠0 must cross a change in the relevant floor exponent; otherwise the weighted sum is proportional to M.

Theorem 1.4 (constant values of the floor in the weights). Let d≥2 and k≥0, and suppose every supported index i satisfies kd≤i<(k+1)d. Then M(c)=(d!)kVd(c). In particular, cancellation on the interval d≤i<2d forces M(c)=0; and if every supported index is at least d while M(c)≠0 and Vd(c)=0, then some supported index is at least 2d.

Proof. On this interval the quotient ⌊i/d⌋ is the constant k, so i!=(d!)kWd,i for every supported index and the factor (d!)k comes out of the sum. The two consequences follow by taking k=1 and by contraposition. ◻

The constant-floor condition is local. A quantitative bound also needs control of the multiplier: cancellation alone does not give 0<M<(R+1)!−1. After proving the common-denominator estimates, §3.1 derives the resulting bounds on R.

At a prime index, the two-term vector Tp already changes only one weighted sum; there are no proper divisors to eliminate.

Theorem 1.5 (changing just one weighted sum). Let p≥3 be prime and let cp−1=p, cp=−1, with every other coefficient zero. Then M(c)=0, Vp(c)=p!−1, and Vd(c)=0 for every d≥2 with d≠p.

Proof. Since p≥3, both indices p−1 and p lie in the prescribed coefficient domain i≥2. Here M=p(p−1)!−p!=0. For 2≤d<p the quotient identity ⌊(p−1)/d⌋=⌊p/d⌋ holds, since d∤p, so the two factorial weights carry the same power of d! and the weighted sum evaluates to p(p−1)!/(d!)⌊p/d⌋−p!/(d!)⌊p/d⌋=0. For d>p both indices lie below d, so both weights are the plain factorials and the same cancellation occurs. At d=p the weights are (p−1)! and p!/p!=1, giving p(p−1)!−1=p!−1. ◻

For p=2, the same auxiliary identities would require c1=2. They do not define an admissible vector on the domain i≥2, since every vector supported on that domain has c1=0.

Adding an integer multiple of this vector changes Vp by that multiple of p!−1, without changing M or any other Vd. Suppose M(c)≠0, V2(c)=⋯=VD(c)=0, and every supported index exceeds a prescribed integer B. Choose a prime p>max(D,B+1,2) and replace c by

c−⌊R(c)+1/2⌋Tp.

Both new indices, p−1 and p, exceed B. Since p>D, the prescribed weighted sums remain zero, and M is unchanged. The new remainder is R(c)−⌊R(c)+1/2⌋∈[−1/2,1/2). The progression construction below supplies such an initial vector: choose its first index r>B.

An integral remainder rounds to zero; a nonintegral one rounds to a nonzero value. Thus rounding controls size but does not prove nonintegrality or exclude further denominators at fixed M. Moreover, ⌊R(c)+1/2⌋ depends on the unknown remainder: this is an existence formula for a representative, not an independent certificate of its nonvanishing.

A primitive solution on an arithmetic progression

We now separate two tasks that the basis classification alone does not combine: prescribe where support begins and make M absorb every fixed denominator as the support grows. Equal spacing turns each cancellation equation into a polynomial root condition. Scaling then makes the coefficients integral and the last coefficient equal to one. Let D,r≥2 and let ℓ be a positive multiple of 2,…,D. The roots will be (d!)−ℓ/d. Put

ij=r+jℓ (0≤j<D),N=r+(D−1)ℓ,αd=(d!)ℓ/d,
∏d=2D(αdX−1)=∑j=0D−1hjXj,A=∏d=2Dαd.

At index ij put the coefficient

cij=N!hjAij!,

and put every other coefficient equal to zero. We verify that these coefficients are integers and that V2=⋯=VD=0. For 2≤d≤D, divisibility of ℓ by d gives ⌊ij/d⌋=⌊r/d⌋+jℓ/d. Hence

Vd(c)=N!A(d!)⌊r/d⌋∑j=0D−1hjαd−j=0,

because αd−1 is a root of the displayed polynomial. Evaluating that polynomial at 1 gives M

M=N!∏d=2D(1−1/αd)>0.

For integrality, each term of hj/A is, up to sign, the reciprocal of a product of D−1−j distinct factors αd. In that denominator, each selected d contributes ℓ/d blocks of size d. The total size is (D−1−j)ℓ=N−ij, so the multinomial coefficient shows that this denominator divides N!/ij!. Thus every cij is an integer. Since hD−1=A, the coefficient at N is 1, and the coefficients have gcd one. This does not establish minimal support. The smallest permitted step is ℓ=lcm⁡(2,…,D); taking D=r=ℓ=2 recovers −6e2+e4.

The stronger divisibility

r!∏d=2D(ℓ/d)!∏d=2D(αd−1)∣M

follows because

N!r!A∏d=2D(ℓ/d)!

is an integer: it counts partitions into one distinguished block of size r and, for each d, an unordered collection of ℓ/d blocks of size d. Their sizes sum to N. Using M=(N!/A)∏d=2D(αd−1) now gives the asserted divisibility. Put v=max(r,ℓ/2), an integer since 2∣ℓ. The displayed factor divides M, so v!∣M. Since D≤ℓ≤2v and r≤v,

N=r+(D−1)ℓ≤4v2−v<4v2,v≥⌊N/2⌋+1.

Consequently (⌊N/2⌋+1)!∣M. Every positive integer q≤⌊N/2⌋+1 therefore divides M, as does every divisor of that factorial. In particular, N≥4(q−1)2 suffices for q∣M. This includes denominators growing with N.

The same vector also supplies the size estimate needed for a tail test. Every factor 1−1/αd lies strictly between 0 and 1, so 0<M<N! and

0<M(S−HN)<2M(N+1)!−1<2N,

where (N+1)!−1>NN! for N≥2. Thus the family supplies both eventual divisibility by each fixed q and the size condition M<(N+1)!−1, without an extra hypothesis. It does not yet compare this tail with the finite part’s gap to the next integer. Both may shrink together. Before studying that arithmetic question, we show why the cofactor formula gives the same primitive vectors.

Cofactors give the same primitive vector

A determinant formula may appear to offer a different family. For positive starting index, however, its primitive normalization gives exactly the preceding progression vector, up to sign. We first verify the determinant construction, then identify the two formulas. For integers n,t≥0, put sn=((n+2)!)2 and ij=(t+j)sn for 0≤j≤n+1. This step size matches the recorded construction; the determinant argument only needs a positive step divisible by every d=2,…,n+2. Here the factorial-weight formula is also allowed at index 0; the case t=0 will show why a lower support condition is necessary. Form the integer matrix A with first row (ij!)j and weighted-sum rows (Wd,ij)j for 2≤d≤n+2.

To see that detA≠0, divide column j by ij!. Since d∣sn, the dth weighted-sum row becomes ((d!)−(t+j)sn/d)j. Removing the nonzero factor (d!)−tsn/d from that row leaves a Vandermonde matrix with nodes

1, (2!)−sn/2, …, ((n+2)!)−sn/(n+2).

These nodes are distinct: (d!)1/d is strictly increasing, since d!<(d+1)d. The Vandermonde determinant is the product of their pairwise differences [32], and is therefore nonzero. None of the row or column factors removed is zero, so the original determinant is nonzero too.

Let c be the cofactor vector of the first row, and let Nd be the determinant obtained by replacing that row with (Wd,ij)j. Cofactor expansion gives M(c)=detA≠0 and Vd(c)=Nd. For 2≤d≤n+2 the replacement duplicates a row, so Vd(c)=0. By (1) each (Vd(c)−M(c))/(d!−1) is an integer and vanishes for d>maxjij, so

Rn,t:=∑d>n+2Ndd!−1,Rn,t−det(A)S∈Z.

A family with minjij=tsn→∞ and Rn,t∉Z would prove irrationality. Indeed, for a fixed denominator q, every entry of the first row is divisible by q once the least support index is at least q, so q∣detA. If S=a/q, the integer-difference identity would then force Rn,t∈Z. Conversely, if S is irrational, that same identity and detA≠0 make every Rn,t irrational. Taking, for example, n=0 and t→∞ gives the required family. Thus nonintegrality along such a family is equivalent to the original irrationality question, not a weaker conjecture obtained from the construction.

The condition tsn→∞ cannot be replaced by unboundedness of (n,t): at t=0, the least supported index stays zero. The coefficient at i0 really is nonzero: its cofactor uses columns 1,…,n+1, and after nonzero row and column scalings it is a Vandermonde determinant on the weighted-sum nodes. The eventual identity Nd=detA determines the tail, but not the finite intermediate sum, whose terms may have either sign. Thus the determinant proves cancellation, but its nonzero value does not prove nonintegrality of the remainder.

For t≥1, primitive normalisation identifies this vector exactly. Take D=n+2, ℓ=sn and r=tsn in the progression construction of §1.4, and put N=r+(D−1)ℓ. For any coefficient vector supported on these D indices, the equations V2=⋯=VD=0 say that the polynomial

∑j=0D−1cijij!Xj

vanishes at the D−1 distinct points αd−1=(d!)−ℓ/d, 2≤d≤D. Its degree is at most D−1, so it is a scalar multiple of ∏d=2D(αdX−1). Thus the rational solution space is one-dimensional. The progression vector has coefficient 1 at N, so every integer solution is its integer multiple, with multiplier cN. In particular, the cofactor gcd is |cN|, and the primitive cofactor vector is the progression vector up to sign. Choosing the sign so that M>0 gives

N!∏d=2D(1−(d!)−ℓ/d).

For example, n=0,t=1 gives support {4,8}, cofactor vector 2520e4−6e8, and primitive vector −420e4+e8 with M>0.

The factorial divisibility proved for the progression coefficient M therefore also applies to this primitive vector. Normalisation still leaves the nonintegrality condition equivalent to irrationality along families whose least support index tends to infinity: M is nonzero, every factorial below that index divides it, and (1) gives the integer-difference identity. This identification does not estimate the distance from an integer. The case t=0 is excluded from the identification with the stated progression construction, whose indices are at least 2.

Prime powers in reduced denominators

Let p be a prime, and let vp(a) denote the exponent of p in a positive integer a. We determine which prime powers of LM survive in den⁡(HM). At the largest exponent, this follows by specialising Louwsma and Martino’s valuation formula for an elementary symmetric sum [4]. Indeed, ∑i1/xi=(∑i∏j≠ixj)/∏ixi for positive integers xi. Their lemma evaluates the numerator’s valuation; subtracting the product’s valuation gives the reciprocal-sum form. The proof below uses the least common multiple instead of the product.

The distinction is already visible in 1/3+1/15=2/5: the common denominator contains 3, but the reduced denominator does not. By contrast, 1/9+1/45=2/15 only lowers the exponent of 3 from 2 to 1. The theorem below tests whether the largest exponent is preserved; a failure of that test need not remove the prime completely. The formula is standard. Its mechanism is simple: after clearing by the lcm and reducing modulo p, only denominators carrying the largest power of p contribute. The factorial examples and later tail comparison concern this particular series.

Theorem 2.1 (maximal prime-power survival). Let M≥2, let p be a prime dividing LM, and put e=vp(LM). Let J={n:2≤n≤M, vp(dn)=e} and write dn=peun for n∈J. Then, with inverses in Fp,

(8)vp(den⁡(HM))=e⟺∑n∈Jun−1≠0in Fp.

Proof. Write LM=peW with p∤W. If vp(dn)<e, then p∣LM/dn. If n∈J, then un(LM/dn)=W, so LM/dn≡Wun−1(modp). Summing over 2≤n≤M gives

(9)AM≡LMpe∑n∈Jun−1(modp).

Finally den⁡(HM)=LM/gcd(AM,LM), whose p-valuation is e exactly when p∤AM. Since W is invertible modulo p, that is exactly the stated condition. ◻

The complete valuation identity is

vp(den⁡(HM))=max{0,e−vp(AM)}.

A zero residue in (8) therefore lowers the exponent, whereas complete cancellation requires vp(AM)≥e. For e>1, reduction modulo p alone cannot distinguish partial from complete cancellation. The proof uses no property of factorials beyond the displayed denominators, so it applies to any finite sum of reciprocals of positive integers. For dn=n!−1, the examples below show that complete cancellation can actually occur; the general valuation formula alone does not predict it.

Two complete cancellations.

A unique maximal exponent makes the residue sum a single nonzero term, so the condition holds automatically. With several maximal terms it can fail. The following examples show cancellation in the factorial sequence itself, not just in arbitrary rational sums. Take M=p−1. The table lists every index with p∣dn, and every one has valuation exactly one.

p

indices n dn/p(modp) inverses modulo p
139 69, 122, 137 6, 49, 73 116, 122, 40
2593 349, 2243, 2591 1508, 1566, 1678 1367, 356, 870

The inverse sums are 278=2⋅139 and 2593. By Theorem 2.1, 139∤den⁡(H138) and 2593∤den⁡(H2592): each prime divides the common denominator of its partial sum and vanishes from its reduced denominator. The listed indices and their valuations are checked by the recurrence

r1=1,rn≡nrn−1(modp2),0≤rn<p2,

run through n=p−1. At an index with rn≡1(modp), the cofactor is (rn−1)/p modulo p. Valuation at least two would give rn=1, which occurs for neither prime among the tested indices 2≤n≤p−1. That enumeration is a finite calculation separate from the theorem.

The recurrence is run modulo p2, not merely modulo p: reduction modulo p would find the indices with p∣n!−1, but would not determine whether p2∣n!−1. This is why the zero inverse sums prove complete cancellation in these examples.

For each prime p≥3, the exponent vp(den⁡(HM)) is constant for M≥max(2,p−2). Wilson’s theorem gives (p−1)!−1≡−2(modp), and n!−1≡−1(modp) for n≥p. Hence every summand after index p−2 has denominator coprime to p. If a reduced fraction a/b is followed by a summand 1/d with p∤d, their sum has numerator ad+b over bd. When p∣b, this numerator is nonzero modulo p, so the exponent of p in the reduced denominator stays the same. When p∤b, the new denominator remains coprime to p. Consequently, the two cancellations persist in every later partial sum:

139∤den⁡(HM)(M≥137),2593∤den⁡(HM)(M≥2591).

The primes still divide the corresponding common denominators LM. These are assertions about the partial sums; they do not constrain a hypothetical denominator of their real limit S without a further approximation estimate.

We next ask where a prime first divides a denominator dm. This means it divides none of d2,…,dm−1; it may still divide a later denominator, so uniqueness at its first occurrence gives no uniqueness in every larger block.

Proposition 2.2 (cofinal first prime occurrences). For every integer B≥0 there are a prime q and an integer m>B with m<q, q∣m!−1 and gcd(q,k!−1)=1 for every k with 2≤k<m.

Proof. Choose a prime q≥B!+5. Wilson’s theorem gives (q−2)!≡1(modq), so there is a least index m≥2 with q∣m!−1, and m≤q−2. If m≤B, then q≤m!−1≤B!−1, contradicting the choice of q. Hence m>B, and minimality gives the coprimality assertions. ◻

The construction proves cofinality of first occurrences, but not the inequality involving q needed for the additional two-factor test in §5. The sufficient tail inequality there does not require choosing a prime q. Wilson reflection limits what can be inferred from a large prime factor alone: for odd 2≤n<q with q∣n!−1, the identity (q−1−n)!n!≡(−1)n+1(modq) gives the reflected congruence at q−n−1. It gives a second summand denominator only when q≠2n+1 and n≤q−3. At the Wilson endpoint n=q−2, it instead returns index 1, where 1!−1=0 is not a denominator of this series. The reflection identity is classical; see Stewart [8].

Growth of the common denominator

A large product of denominators gives a large least common multiple only if their overlap is controlled. We use the last k denominators: their logarithmic product contributes about kNlog⁡N, while the pairwise-gcd estimate costs about (k3/6)log⁡N. With k∼αN, the leading coefficient is α−α3/6, maximized at α=2. The lemmas below make this balance rigorous and yield lim inflog⁡LN/(N3/2log⁡N)≥22/3. Here LN is the common denominator before reduction, not den⁡(HN).

Lemma 3.1 (product, least common multiple, pairwise gcd). For positive integers x1,…,xk,

∏i=1kxi | lcm⁡(x1,…,xk)∏i<jgcd(xi,xj).

Proof. For each prime r, arrange ai=vr(xi) in increasing order. The valuation of the right-hand side is ak+∑i=1k−1(k−i)ai, whereas that of the left-hand side is ∑i=1kai. Their difference is ∑i=1k−1(k−i−1)ai≥0. This proves divisibility at every prime. ◻

The next divisibility is the case P=−1 of the relation gcd(i!+P(i),j!+P(j))∣(j!/i!)P(i)−P(j) for P∈Z[X], obtained by multiplying i!+P(i) by j!/i! and subtracting j!+P(j). Luca and Shparlinski use this relation to bound a common divisor of n!+P(n) and (n+h)!+P(n+h) in the proof of their Lemma 5 [13], and Lai uses it for a common prime-power divisor [14].

The earlier spacing method of Erdős and Stewart [18] is related background. We now state only the subtraction needed for the lcm argument, rather than importing the prime-factor estimates of those papers.

Lemma 3.2 (factorial-gap gcd). For 2≤i<j, the integer g=gcd(i!−1,j!−1) divides j!/i!−1, and g≤j!/i!−1<jj−i.

Proof. The identity (j!/i!)(i!−1)−(j!−1)=1−j!/i! shows that g divides j!/i!−1. This integer is positive. Moreover, j!/i! is a product of j−i integers at most j, so g≤j!/i!−1<jj−i. ◻

Lemma 3.3 (segment inequality). For 2≤k≤N−1,

(10)∑n=N−k+1Nlog⁡(n!−1) ≤ log⁡LN+(k+13)log⁡N.

Proof. Apply Lemma 3.1 to the k denominators with indices N−k+1,…,N. Their lcm divides LN. Lemma 3.2 bounds the logarithm of each pairwise gcd by (j−i)log⁡N. Since

∑N−k+1≤i<j≤N(j−i)=∑h=1k−1h(k−h)=(k+13),

taking logarithms gives (10). ◻

Theorem 3.4 (common-denominator growth).

lim infN→∞log⁡LNN3/2log⁡N ≥ 223.

Proof. The sum of the logarithms on a terminal block of length k is of order kNlog⁡N, while the gcd estimate costs order k3log⁡N. We therefore fix α>0 and take k=⌊αN⌋. For large N this is an admissible block length. If u=N−k+1, then n!−1≥n!/2 and log⁡n!≥nlog⁡n−n give

∑n=uNlog⁡(n!−1)≥k(ulog⁡u−u)−klog⁡2.

Here k=αN+O(1) and u=N+O(N), so kulog⁡u=αN3/2log⁡N+O(Nlog⁡N). The remaining two terms are O(N3/2). Also

(k+13)log⁡N=α36N3/2log⁡N+o(N3/2log⁡N).

After division by N3/2log⁡N, Lemma 3.3 therefore yields a lower limit of at least α−α3/6. This expression is maximised at α=2, giving 22/3. ◻

The displayed proof of the lcm bound uses only Lemmas 3.1 and 3.2. Theorem 12 of Garaev, Luca and Shparlinski [7] supplies a different input: a uniform O(N2/3) bound for the number of solutions of n!≡a(modp), for fixed a≢0(modp), in an interval of length N contained in 1≤n<p. The weaker lcm deduction from it is retained in §B.3, not used here. No exhaustive priority claim is made for the lcm estimate; its precise elementary derivation and the inherited subtraction are stated so that the claim can be assessed without relying on a search conclusion.

Remark (polynomial shifts). Fix P∈Z[X]∖{0}. By Lemma 3 of Luca and Shparlinski [13], in the form given by Lai [14], there is n0≥2, depending only on P, such that n!+P(n)>1 for all n≥n0 and P(n)(n+1)⋯(n+h)≠P(n+h) for all n≥n0 and h≥1. The nonvanishing has a short direct proof. Beyond a cutoff, P has constant nonzero sign and P(n+1)/P(n)→1, so 0<P(n+1)/P(n)<n+1 for every sufficiently large n. Multiplying these inequalities gives

0<P(n+h)P(n)<∏j=1h(n+j)(h≥1).

Thus the required subtraction is nonzero for every gap h, not merely for each fixed h after a separate cutoff. Increasing n0 also makes n!+P(n)>1. The proof of Theorem 3.4 then gives

lim infN→∞log⁡lcm⁡{n!+P(n):n0≤n≤N}N3/2log⁡N ≥ 223.

For n0≤i<j≤N, multiplying i!+P(i) by j!/i! and subtracting j!+P(j) gives

gcd(i!+P(i),j!+P(j)) | j!i!P(i)−P(j),

which replaces Lemma 3.2. The right-hand side is nonzero by the choice of n0, and its absolute value is at most CPNj−i+deg⁡P for a constant CP. Over a block of k indices the pairwise loss therefore grows by OP(k2log⁡N), which is OP(Nlog⁡N) when k=⌊αN⌋. Since n!+P(n)≥n!/2 for large n, the lower estimate for the sum of logarithms is unchanged. The cutoff is needed, since for P=−2 the factor at n=2 vanishes. The hypothesis P≠0 is needed as well, since for P=0 the least common multiple is N! whatever the cutoff.

A non-polynomial comparison.

For n!+2n−1, Luca and Shparlinski [19] eliminate the exponential term using three indices satisfying the congruence; multiplicative order enters the counting. Their valuation-layer identity [19] records repeated prime-power divisibility. These are related tools, not inputs to Theorem 3.4, and no extension to general non-polynomial perturbations is asserted.

Within the terminal-block estimate just proved, maximising α−α3/6 gives 22/3. This is an optimisation of that one-parameter bound, not an optimality theorem for all block selections or all gcd arguments. The discarded factorial term has size N3/2, so a finite normalised ratio need not be close to the limiting lower bound.

Theorem 3.4 concerns LN, not den⁡(HN). Theorem 2.1 and the cancellations at 139 and 2593 show why these cannot be interchanged: a prime present in the common denominator may be absent after reduction. The two examples establish this possibility, not an asymptotic estimate for cancellation.

Even two consecutive denominators obstruct clearing by the full lcm. For N≥3, the gcd of (N−1)!−1 and N!−1 divides N−1, whereas (N−1)!−1 is coprime to N−1. Thus those denominators are coprime, and

LN(S−HN)>((N−1)!−1)(N!−1)(N+1)!−1⟶∞.

This argument does not use the stronger asymptotic growth theorem.

Reduction could still help, but only in the opposite direction: a clearing argument needs a sufficiently small reduced denominator. Write HN=A^N/QN in lowest terms, reserving AN for the numerator over the full lcm. Under S=a/q, with a∈Z and q≥1, positivity gives

qQN(S−HN)=aQN−qA^N∈Z>0.

A contradiction requires qQN(S−HN)<1 at some N. A lower bound on QN cannot supply this upper estimate. For comparison, the rational limit 1 and the approximants 1−1/QN have arbitrarily large reduced denominators without any irrationality. The inequality dN(S−HN)<1 alone would not repair the argument, because dN generally does not clear the earlier summands of HN. The prime-power analysis describes which factors survive reduction, but it does not yet supply the required upper control on QN.

For linear forms in p-adic zeta values, with p≥5 a fixed prime and n the index of the form, Lai, Lupu and Sprang quantify a saving of this kind. They bound the denominators of the coefficients by powers of lcm⁡(1,…,n), show that the bound may be divided by a product Φn of prime powers, and compute the growth rate of Φn [15]. That saving enters the inequality [15] under which the rescaled forms, which have integer coefficients, meet the irrationality criterion [15] that they quote from Lai [16]. Along one unbounded subsequence, the forms must be nonzero and their p-adic absolute values, multiplied by the largest ordinary absolute value of an integer coefficient, must tend to zero. The coefficient bound is archimedean; the bound on the value is p-adic. The conclusion is that at least one of the finitely many p-adic numbers in the forms is irrational, not that each is irrational. The corresponding saving for HN is the ratio LN/QN.

This ratio is an integer. Merely naming it supplies neither its growth nor a small linear form; the analogy is about what an eventual estimate would have to accomplish. Theorem 2.1 decides, for each prime p∣LN, whether p divides this ratio; no asymptotic lower bound for the ratio is established here. This paragraph compares the two methods; no p-adic theorem is applied to S.

Consequences for the support parameter

The lcm estimates also constrain the size of a cancellation vector when its multiplier is bounded by a factorial. At depth D=2t2, cancellation gives LD∣M; the next theorem additionally assumes 0<M<(R+1)!−1. Its R is an index, not the remainder or the modulus Rp of §5. The conclusion bounds the largest supported index only when R denotes that index. Taking R to be a later cutoff gives a bound on the cutoff instead.

Theorem 3.5 (a lower bound for the support parameter). Let t,M,R∈N satisfy

t≥232,M>0,L2t2∣M,M<(R+1)!−1.

Then 3t3<2(R+1). Consequently no family satisfying these hypotheses for all sufficiently large t has (R(t)+1)/t3≤3/2 eventually, and none has R(t)=o(t3).

Proof. Put D=2t2, k=2t and u=D−k+1=2t2−2t+1. Lemma 3.3 at N=D, with n!−1≥n!/2 and log⁡n!≥nlog⁡n−n, and with LD≤M<(R+1)!, gives

(11)2t(ulog⁡u−u−log⁡2)<(R+1)log⁡(R+1)+(2t+13)log⁡(2t2).

Suppose R+1≤32t3. Since t≥16 one has 158t2≤u≤2t2 and log⁡u≥2log⁡t, so the left side of (11) is at least 152t3log⁡t−4t3−2tlog⁡2. Using (2t+13)<43t3 and log⁡32<log⁡2, the right side is at most

32t3(log⁡2+3log⁡t)+43t3(log⁡2+2log⁡t).

Dividing by t3 and rearranging, (11) would require

13log⁡t<4+2log⁡2t2+176log⁡2.

But log⁡t≥32log⁡2 and 23<log⁡2<1 give 13log⁡t−176log⁡2≥476log⁡2>5, while 4+2log⁡2/t2<5. This contradiction proves the finite inequality, and the two sequence conclusions follow. ◻

The hypotheses are compatible: take M=L2t2 and then choose R large enough. The conclusion rules out subcubic R under the stated size condition; it neither constructs a coefficient vector nor proves irrationality.

Theorem 3.4 raises the asymptotic constant in the same estimate.

Corollary 3.6 (the asymptotic lower bound). Let M(t),R(t) satisfy M(t)>0, L2t2∣M(t) and M(t)<(R(t)+1)!−1 for all sufficiently large t. Then

lim inft→∞R(t)+1t3 ≥ 169.

Proof. Put r=R(t)+1, so L2t2≤M(t)<r! and log⁡L2t2<rlog⁡r. Since

(2t2)3/2log⁡(2t2)t3log⁡t⟶42,

Theorem 3.4 gives, for every ε>0 and all sufficiently large t,

log⁡L2t2≥(163−ε)t3log⁡t.

Suppose r≤ct3 for arbitrarily large t, with 0<c<16/9 fixed. On that subsequence rlog⁡r≤ct3(3log⁡t+log⁡c)=(3c+o(1))t3log⁡t. Choosing 0<ε<16/3−3c contradicts the lower bound for sufficiently large t. ◻

This corollary does not assert the strict inequality R+1>169t3 eventually. Indeed, the finite logarithmic constraint (11) allows equality at that constant. To check this directly, set R+1=16t3/9 with t≥4. Since u=2t2−2t+1≤2t2 and xlog⁡x−x is increasing for x≥1, the left side of (11) is at most 4t3(log⁡(2t2)−1)−2tlog⁡2. Subtracting this bound from the right side gives

t3(4+169log⁡169−83log⁡2)+t(53log⁡2−23log⁡t)>43t3−23t2>0,

where log⁡2<1 and log⁡t≤t suffice for the last comparison. Integer pairs with 9(R+1)=16t3 exist for every multiple t≥6 of 3. This verifies the numerical constraint, not the existence of a coefficient vector with that support. The asymptotic assertion in Corollary 3.6 instead follows from Theorem 3.4.

Carries and factorial digits

Instead of clearing the partial sum, we now track the only integer that could equal m!S. Under rationality, m!S is integral for all large m; since the positive tail is less than one, it must be the strict successor of m!Hm. The recurrence for these successors therefore tests rationality without assuming that m! clears Hm. Put

Zm=⌊m!Hm⌋+1,Δm=Zm−1−(m−1)!Hm−1(m≥3).

Thus 0<Δm≤1, with Δm=1 when the preceding scaled partial sum is integral. Define the integer carry bm by

(12)Zm=mZm−1+1−bm.

For example, 2!H2=2 and 3!H3=36/5, so Z2=3 and Z3=8, giving b3=3⋅3+1−8=2. We use ⌊x⌋+1, not ⌈x⌉: at the integral value 2!H2=2, the required strict successor is 3.

Call bm=1 a unit carry. Put Em=m!(S−Hm) and εm=1/(m!−1). Since dn+1>(n+1)dn, a geometric majorant gives the tail estimate

(13)0<Em<2m!(m+1)!−1<2m≤1(m≥2).

Compare the factorial scaling in Hančl and Tijdeman’s tail-integrality lemma for factorial series with integer coefficients [11]. Their scaled partial sums are integers; ours need not be, since already 3!H3=36/5. The next proof instead uses the short positive tail to identify the least integer above a scaled partial sum under a rationality assumption. The carry-defect expansion in §B.2 gives a separate, direct application of factorial-series rationality criteria.

Theorem 4.1 (an exact criterion from successive partial sums). For m≥3,

(14)bm=1⟺m∣Zm⟺1+εm<mΔm≤2+εm.

Moreover

(15)S∈Q⟺bm=1 for all sufficiently large m,(16)S∉Q⟺∀B ∃m>B: m∤Zm.

If S=a/q with a∈Z, q≥1 and bm≠1, then q∤(m−1)! and q≥m.

Proof. From m!Hm=mZm−1−mΔm+1+εm one gets bm=⌈mΔm−1−εm⌉ with −1≤bm≤m−1, which gives both equivalences in (14), endpoints included.

Suppose S=a/q and q∣(m−1)!. For j=m−1 and j=m the number j!S is an integer, and j!Hj=j!S−Ej with 0<Ej<1 by (13), so Zj=j!S. Substituting into (12) forces bm=1. Every fixed q divides (m−1)! eventually, which gives one direction of (15) and, by contraposition, the final assertion; q<m would imply q∣(m−1)!. Conversely, an eventual unit-carry tail makes Zm/m! eventually constant, and 0<Zm/m!−Hm≤1/m! with Hm→S identifies that constant as S, which is then rational. Negating (15) gives (16). ◻

The conclusion q∤(m−1)! is about prime-power multiplicities, not just prime factors. It permits all prime factors of q to be at most m−1 if one occurs to a larger exponent than in (m−1)!.

The adjacent-unit window condition

Two successive unit carries imply a prime-power size condition after clearing denominators in their defining inequalities. The calculation below shows that the condition follows already from those carries. An additional arithmetic estimate is therefore needed to turn it into a contradiction.

Write Δk=PkΔ/QkΔ in lowest terms, with PkΔ,QkΔ>0. In one recurrence step the unreduced denominator is QkΔ(k!−1); let Gk=QkΔ(k!−1)/Qk+1Δ be the positive integer cancelled in reduction. For m≥3, put

Dm=QmΔ(m!−1)((m+1)!−1),Ωm=Dm(m(m+1)Δm−m−2−m+1m!−1−1(m+1)!−1).

Both are integers. The two unit carries are equivalent to 0<Ωm≤Dm. Independently of the carries, cancellation in two successive steps gives

Dm=Qm+2ΔGm+1Gm=QmΔ(m!−1)((m+1)!−1),

and under the pair assumption the offset factors to match, Ωm=Pm+2ΔGm+1Gm. Dividing by Gm+1Gm>0 reduces the proposed bound to

0<Pm+2Δ≤Qm+2Δ,

which every reduced positive gap in (0,1] already satisfies. Thus the resulting gap bound is not an additional restriction on a pair of unit carries. This calculation does not exclude an independent arithmetic obstruction involving the transition numerators or their valuations. Such an obstruction would need information not already implied by the pair assumption; it need not take the form of a bound on the preceding gap Δm.

Factorial digits of S−e+2

The small perturbation from the reciprocal-factorial series can also be isolated term by term. Subtracting e−2 leaves the absolutely convergent correction C=∑n≥2(n!(n!−1))−1. The identity 1/(n!−1)=1/n!+1/(n!(n!−1)) gives

(17)S=C+e−2.

The canonical factorial digits of a real x are am(x)=⌊m!x⌋−m⌊(m−1)!x⌋∈{0,…,m−1} for m≥2. At a rational endpoint, the floor convention selects the terminating representation. For example, 1/2=1/2! also equals ∑m≥3(m−1)/m!, but its canonical digits are a2=1 and am=0 for m≥3, not an eventually maximal tail. With this convention, x is rational exactly when its digits vanish from some index on. The criterion goes back to Cantor [5]; Galambos treats the rationality of Cantor series in [6], and Koepf and Schmersau prove the irrationality direction for nonterminating factorial expansions whose digits are not eventually maximal [9].

Proposition 4.2 (rationality and factorial residues).

S∈Q⟺⌊m!C⌋≡−2(modm)for all sufficiently large m.

Proof. Put Jm=∑k=0mm!/k!. Then Jm is an integer, every summand with k<m is divisible by m and the summand at k=m is 1, so Jm≡1(modm); and 0<m!e−Jm<1 for m≥2.

Suppose S is rational. For all large m both (m−1)!S and m!S are integers, so m∣m!S. Multiplying (17) by m! and using ⌊−m!e⌋=−Jm−1 gives ⌊m!C⌋=m!S+2m!−Jm−1≡−2(modm).

Conversely, assume the congruence from some index on. Because 0≤am(C)<m, for m≥3 the congruence is equivalent to am(C)=m−2. Choose N beyond the exceptional indices. The canonical expansion gives

C=⌊N!C⌋N!+∑m>Nm−2m!,

and the telescope ∑m>N(m−1)/m!=1/N! turns this into

S=C+e−2=⌊N!C⌋+1N!+∑m=2N1m!,

which is rational. ◻

The proof applies to every real C, not just this series: eventual digits am(C)=m−2 characterize rationality of C+e−2. Thus C=2−e satisfies the eventual pattern, whereas C=0 does not. For our particular positive C, neither possibility is decided by a finite digit list.

Negating the eventual congruence gives the equivalent formulation

(18)S∉Q⟺(∀B)(∃m>B) (⌊m!C⌋+2)modm≠0.

It remains to prove failures of this congruence at arbitrarily large indices for this particular C. A finite list of failures does not settle this question, and the equivalence gives no distribution estimate.

Recovering the finite carries.

To compare canonical digits with the carries, we must account for the tail removed before taking a floor. Set

δm=m!∑n>m1n!(n!−1),σm=1{{m!C}<δm}.

Here {x}=x−⌊x⌋. The factorial-tail bound gives 0<δm<1/((m+1)!−1)<1. Subtracting δm from m!C lowers its floor by one exactly when {m!C}<δm. Therefore

Zm=m!∑n=2m1n!+⌊m!C⌋+1−σm.

Substitution in (12) gives

(19)bm=m−1−am(C)+σm−mσm−1(m≥3).

When {m!C}=δm, the scaled partial sum is integral and σm=0. For a smaller fractional part, subtracting even this small tail crosses an integer. Tail size alone therefore permits neither correction term to be omitted.

Escape from a smaller interval

Write θm={m!S}. The next criterion compares a fractional part with the actual tail: θm−1 must leave [0,Em/m) at arbitrarily large indices. The interval has width less than 2/m2. Since Em contains S, this equivalence is not yet a finite test; the second assertion supplies a sufficient test using the finite quantities Δm and εm.

Proposition 4.3 (lower-interval criterion).

(20)S∉Q⟺∀B ∃m>B: Em≤mθm−1.

For m≥3 the finite condition

(21)mΔm≤1+εmor1+εm+2m≤mΔm

implies the escape inequality in (20). Cofinally many instances of (21) therefore imply S∉Q.

Proof. If S is rational, the remainders θm−1 vanish eventually while Em>0, so the escape inequality fails eventually. Conversely, suppose mθm−1<Em for every sufficiently large m. By (13) this gives mθm−1<1, so am(S)=0 and θm=mθm−1. Fix N after this recurrence begins. Then θN+k=((N+k)!/N!)θN<1 for every k≥0, forcing θN=0. Hence N!S is an integer and S is rational.

For the finite implication, suppose mθm−1<Em. The tail recurrence is mEm−1=1+εm+Em, so Em−1>Em/m>θm−1≥0, and with 0<Em−1<1 the strict-successor identity gives Δm=Em−1−θm−1. Hence

1+εm<mΔm=1+εm+Em−mθm−1≤1+εm+Em<1+εm+2m,

which contradicts (21). ◻

There are two comparisons to keep separate: a finite sufficient test for escape at one index, and an equivalence when indices are unbounded. The finite test excludes only the open interval 1+εm<mΔm<1+εm+2/m; it accepts either endpoint and all values outside. The width 2/m refers to mΔm; the interval for Δm has width 2/m2. At one index this is a sufficient test for escape, obtained by replacing Em with the upper bound 2/m.

At arbitrarily large integer indices, however, the finite test itself characterises irrationality:

S∉Q⟺(21) holds for arbitrarily large integers m.

Only the forward implication remains to be shown. If S is irrational, Theorem 4.1 gives bm≠1 at arbitrarily large indices. By (14), each such index satisfies mΔm≤1+εm or mΔm>2+εm. Since 2/m<1 for m≥3, either alternative implies (21). The converse is Proposition 4.3. This combines the two criteria; it does not prove either unbounded occurrence for S. The equivalence is over all integer indices. A restriction to primes gives only a sufficient condition here.

For comparison at a single index, the exact classification is

(22)Em≤mθm−1⟺bm≠1  or  am(S)=m−1(m≥3).

To verify the classification, use 0<Ej<1 to write

Zj=⌊j!S⌋+1{θj≥Ej}.

Substitution in the carry recurrence gives

bm=1−am(S)+m1{θm−1≥Em−1}−1{θm≥Em}.

Since 0≤am(S)≤m−1, a unit carry has exactly two possibilities: both indicators are zero and am(S)=0, or both are one and am(S)=m−1. In the first case mθm−1=θm<Em; in the second, mθm−1≥m−1>Em. Conversely, mθm−1<Em<1 forces am(S)=0 and both indicators to be zero, using mEm−1=1+εm+Em. It therefore forces a unit carry with zero digit. This proves (22), including the equality case in its left-hand inequality.

Thus escape can occur even when bm=1, provided the digit is m−1. This happens at m=52: exact rational arithmetic gives b52=1 together with 52Δ52>1+ε52+2/52, the margin being about 0.5689674908. This example separates the two events at a single index: the finite test can hold at a unit carry. Their occurrence at arbitrarily large indices nevertheless gives equivalent rationality criteria. Since 52 is composite, the example supplies no prime-index instance.

The series with denominators n!+t

The digit argument also applies to the other shifts mentioned by Erdős. For an integer t≥−1, put St=∑n≥21/(n!+t) and Ct=∑n≥21/(n!(n!+t)), so that St=−tCt+(e−2).

The restriction t≥−1 keeps every denominator positive for n≥2; t=−2 is excluded because its first denominator is zero. All series in this identity converge absolutely. The familiar case t=0 is e−2; this example will make the residue condition explicit.

Theorem 4.4 (a criterion for the shifts t≥−1). For every integer t≥−1, the series St is rational exactly when

⌈tm!Ct⌉≡2(modm)

for all sufficiently large m, and irrational exactly when that residue is missed cofinally. The member t=−1 is S, and the member t=0 is e−2, for which the scaled correction is always 0 and hence misses the residue class at every m≥3, proving the irrationality of e.

Proof. Put Y=−tCt, so that St=Y+e−2. The proof of Proposition 4.2 applies to any real Y: rationality of Y+e−2 forces ⌊m!Y⌋≡−2(modm) eventually. Conversely, that congruence forces the canonical digits of Y to equal m−2 eventually; adding the factorial series of e−2 leaves an eventually telescoping tail, hence a rational sum. Finally ⌈tm!Ct⌉=−⌊m!Y⌋ translates the congruence. Negation gives the cofinal statement. For t=0, the left side is zero, which is not congruent to 2 for m≥3. ◻

The t=0 specialisation recovers the classical irrationality of e. For t≠0, no estimate proving failure of the displayed congruence at arbitrarily large indices is supplied here. This is a statement about the scope of the present argument, not an assertion that every other member has the same current literature status.

A complementary-denominator criterion

We seek a multiplier with two properties: it absorbs each fixed rational denominator eventually, and the scaled partial sum has a computable gap to the next integer. Include a factorial to ensure the first property. For the second, the multiplier also absorbs every prime-power level shared by two summand denominators. Any exponent left in the denominator then comes from a unique maximal occurrence, so it cannot cancel on reduction. This multiplier is sufficient, not necessarily minimal. For an integer p≥3, not necessarily prime, put Ip={2,…,2p−1} and Fp=(p−1)!, and define

Dp=lcmi,j∈Ipi<j⁡gcd(di,dj),Cp=lcm⁡(Fp,Dp),Lpblk=lcm⁡(Fp,d2,…,d2p−1),Rp=Lpblk/Cp.

Prime by prime, Dp records the second-largest exponent, Cp raises it to at least the exponent in Fp, and Rp records what remains below the largest exponent. The factors Cp and Rp need not be coprime: a prime absent from Fp with largest exponent 5 and second-largest exponent 2 occurs to powers 2 and 3 in the two factors. This is an illustration of the valuation calculation, not a claimed pattern in a particular factorial block. It will be useful to retain the ratio C~p=Cp/Fp, so that

(23)Lpblk=CpRp=FpC~pRp.

Including Fp in both Cp and Lpblk ensures that a fixed rational denominator divides Cp for all large p. The ratio C~p removes precisely that factorial factor when we estimate the shared part. Put

Tp=∑i∈IpLpblkdi,ρp=(−Tp)modRp,Kp=2p2(2p−1)!,

with ρp the least nonnegative representative. All of these depend only on a finite partial sum.

For each prime, Cp contains the larger of its exponent in Fp and its second-largest exponent among the denominators. A prime dividing Rp therefore has a unique denominator whose exponent is maximal and exceeds the exponent in Fp, so the argument of Theorem 2.1 applies to it and gives gcd(Tp,Rp)=1. Thus Rp=den⁡(CpH2p−1), the reduced denominator of the scaled partial sum, not generally of H2p−1. When Rp>1, this ratio is nonintegral and

(24)ρpRp=⌈CpH2p−1⌉−CpH2p−1∈(0,1),

the distance to the next integer. When Rp=1, ρp=0 but the scaled partial sum is integral, so the least strictly larger integer is at distance 1. The theorem below uses only Rp>1.

At p=3, the denominators 1,5,23,119 are pairwise coprime. Thus D3=1, C3=2, R3=13685, and

C3H5=3426413685,ρ3R3=3−3426413685=679113685.

The tail estimate in the proof below gives C3(S−H5)<7/1080<6791/13685. This verifies the comparison for one partial sum. The theorem needs this comparison at arbitrarily large p, so that Fp can absorb any fixed rational denominator.

Theorem 5.1 (a sufficient tail inequality). Suppose that for every B there is a natural parameter p≥3 with p>B, Rp>1, and

(25)(2p+1)Lpblk<Kpρp.

Then S is irrational.

Proof. Suppose S=a/q with q≥1, and choose a parameter in the hypothesis with p>q. Since q∣Fp∣Cp, the number CpS is an integer strictly above CpH2p−1. The gap identity (24) gives

Cp(S−H2p−1)≥ρpRp.

On the other hand, 1/(n!−1)<2/n! for n≥2p, and comparison with a geometric series of ratio 1/(2p+1) gives

S−H2p−1<2(2p)!2p+12p=2p+1Kp.

Multiplication by Cp and (25), with Lpblk=CpRp, now give Cp(S−H2p−1)<ρp/Rp, a contradiction. ◻

Cancelling Rp in (25) through (23) displays the quantitative issue directly:

(26)(2p+1)Cp<KpρpRp.

The factor Rp has cancelled. Thus new prime divisors, even large ones, do not by themselves bound Cp or keep ρp/Rp away from zero. Proposition 2.2 supplies the primes but not these two estimates.

Under Rp>1, we have ρp>0, and taking logarithms gives

log⁡C~p−log⁡ρpRp<log⁡2p2(2p−1)!(2p+1)(p−1)!=plog⁡p+(2log⁡2−1)p+O(log⁡p).

Both terms on the left must be controlled at the same p. The elementary bound ρp≥1 gives only ρp/Rp≥1/Rp, which may be too small. For a concrete comparison, a gap of at least 1/2 would reduce the sufficient test to (2p+1)Cp<Kp/2. At the opposite extreme, the gap 1/Rp requires the stronger bound (2p+1)Cp<Kp/Rp. Neither type of favourable behaviour is proved here for an unbounded family of these particular prefixes.

A more restrictive test also uses a prime q that first divides m!−1 at an index m≥4. Choose p=⌊m/2⌋+1, which need not be prime. This choice puts m in the block and gives q∣Rp. Indeed, q>m≥p, so q∤Fp, and no earlier denominator contains q. The only possible later index in the block is m+1; when it occurs, (m+1)!−1≡m≢0(modq). Thus q occurs in just one block denominator and is absent from Cp.

The additional comparison uses the integer Rp/q. Requiring the tail bound to be smaller than both ρp and Rp/q gives min(ρp,Rp/q) and is equivalent to the following two inequalities:

(27)(2p+1)Lpblk<Kpmin(ρp,Rp/q)⟺{(2p+1)Lpblk<Kpρp,(2p+1)Cpq<Kp.

The equivalence follows by testing each entry of the minimum and using (23). Both inequalities must hold for the same arbitrarily large parameters. The first alone is already sufficient by Theorem 5.1; the second is an extra restriction, not an equivalent formulation of irrationality.

A fixed small index cannot supply new factors indefinitely. If n is fixed and p>n!−1, then n!−1 divides (p−1)!, so none of its prime powers remains outside the factorial base. For example, the factor 719 of 6!−1 is useful only over a finite range of p. Any argument using such factors at arbitrarily large p must let their indices increase.

Finite denominator exclusions

Two recorded computations give different restrictions on a hypothetical rational denominator. The carry calculation concerns divisibility, and the continued-fraction enclosure concerns size:

(28)q∤299999!,q≥239990>1012038

These restrictions hold simultaneously whenever S=a/q, with a∈Z and q>0; they are not two bounds to compare by taking a maximum. They also hold for an unreduced representation. The size bound passes to multiples of the reduced denominator, as does failure to divide a fixed factorial. The computations were performed outside Lean. The candidate finite-size certificate records that its computation has not been kernel checked, and the full carry certificate has the same outstanding formal obligation.

The two restrictions concern different properties of q. The prime 300007 does not divide 299999!, but lies below the size bound; 299999! exceeds that bound, but fails the nondivisibility condition. A sufficiently high power of a small prime need not divide 299999! either. Thus the carry exclusion is not an exclusion of all denominators with small prime factors. In Sondow’s terminology, it says that the least r for which q∣r! is at least 300000 [21]. His approximation estimates for e are not being applied to S.

Carry intervals.

The GMP calculation encloses each scaled sum by outward-rounded rational intervals, using scale 25025679 and 96 guard bits. It uses integer arithmetic throughout. A carry is entered only after the entire interval has been located in one half-open unit cell. All 299998 cells for 3≤m≤300000 were certified in the recorded run. The indices with bm=1 were exactly

52,591,1030,1407,1438,2164,4258,10991,21236.

In particular b300000≠1, which gives the first restriction in (28) by Theorem 4.1, and hence also q≥300000. The independently implemented GMP census is separate from the earlier calculation through 4000. That calculation gives the first six entries of this list and a non-unit carry at its endpoint, so it alone gives q≥4000. The implication from a non-unit carry to factorial nondivisibility is the recorded kernel-checked result; the full census is not thereby a kernel-checked computation.

A rational enclosure and continued fractions.

Set D=280000 and let N be the first integer satisfying N!−1>D. Put

ℓ=∑n=2N−1⌊Dn!−1⌋,u=ℓ+(N−2)+⌊2DN!−1⌋+1.

Each of the N−2 rounded terms loses less than one. Moreover, (n+1)!−1>(n+1)(n!−1), so the positive tail from N is less than (N+1)/(N(N!−1))<2/(N!−1). These inequalities give the strict enclosure ℓ/D<S<u/D. The omitted positive tail makes the lower endpoint strict, and the final +1 makes the upper endpoint strict even when 2D/(N!−1) is an integer.

The recorded values are N=7054 and u−ℓ=7053. Apply Euclid’s continued-fraction algorithm simultaneously to the two endpoints, stopping as soon as their integer parts differ or a remainder vanishes. At an inversion, reverse the order of the endpoints. There are 23449 common partial quotients a0,…,a23448. Their convergent denominators satisfy

Q−2=1,Q−1=0,Qj=ajQj−1+Qj−2,Q23448≥239990.

Every rational in the open enclosure shares these quotients. Such a rational either equals the last convergent P23448/Q23448, or has a subsequent complete quotient A/B>1, with coprime positive integers A,B, and equals

P23448A+P23447BQ23448A+Q23447B.

The determinant identity PjQj−1−Pj−1Qj=±1 proves that this fraction is reduced and its denominator is at least Q23448 [17]. No assumption that the last convergent lies inside the enclosure is needed. The integer inequality 239990>1012038 then gives the decimal bound in (28). The integer program and its recorded output give the 80000-bit scale, width 7053, prefix length 23449, binary exponent 39990 and strict decimal exponent 12038.

Carry computation record.

The carry calculation uses the program scripts/check_erdos68_strict_successor.py with the GMP backend scripts/check_erdos68_strict_successor_gmp.cpp; its output verification/erdos68-strict-successor.json records the scale and guard bits above, the unit carries, an event-trace digest and the final enclosure. The SHA-256 digests of the program and the backend, and the payload digest recorded in the output, are

Program: f25b5bd8ffbfd66bdbb42ca6a07a936aab9b4eb7d3a57446725be55acb642191.
Backend: 43308ddc3902dd537b789ae6e0648054a0d752adf07e5d6e1d829811e9342083.
Payload: 044ef27a60f7294c39b8279b7f04e22e1234a942e2862e9e8c64eb66505af29b.

Both calculations cover finite data. Refining an enclosure can leave the common prefix unchanged, and a later unit carry gives no new carry exclusion. In particular, neither calculation excludes an eventual tail of unit carries. The continued-fraction identities are classical [17]; the numerical parameters in this section are the outputs of the stated computations.

Further arithmetic questions

It remains to prove m∤Zm at arbitrarily large indices. Section 4 gives equivalent formulations through the actual remainders in (20) and the finite condition (21), when all integer indices are allowed. The estimates below concern the remaining arithmetic, not further ways to solve the coefficient equations. None establishes the required unbounded set. Restricting indices to primes is an extra condition; selected congruences likewise do not identify the sequence defined by S without its floor relations.

Value sets do not determine the weighted residue sum.

Garaev, Luca and Shparlinski’s harmonic-sum estimate [20] concerns unconditioned harmonic sums, not sums restricted to the indices where n!≡1(modq). The value-set literature addresses different invariants. Banks, Luca, Shparlinski and Stichtenoth [25] provide historical missing-value context; Klurman and Munsch [26] distinguish unconditional averages from GRH-dependent improvements. Grebennikov, Sagdeev, Semchankau and Vasilevskii [27] study sizes of value and quotient sets, as do the factorial-residue and representation results of Hu [22] and Garaev and Pardo [23]. None is a noncancellation theorem for the reciprocal unit weights here. A lower bound for the number of distinct residues does give an upper bound on the multiplicity of residue 1, since every other value needs at least one index. But that bound does not determine the reciprocal unit weights or test whether they cancel modulo q. The two explicit cancellations in §2 illustrate why the weights matter. Hu’s September preprint [24] considers pairwise distinct maps Tθ(x)=bθ+cθ/x over Fp, with cθ≠0. An internal transition is a pair (x,θ) for which x∈A∖{0} and Tθ(x)∈A. With at least M distinct such pairs, at most |A| parameters, and a bound μ≥1 on the number of ordered parameter pairs representing each nonidentity map TηTθ−1, the theorem gives |A|≫min{M,p}8/15μ−4/15. Its incidence estimate comes from Stevens and de Zeeuw [28]. To apply this to the residue classes of indices satisfying n!≡1(modq), one would first have to specify the set A and the transition maps. The number of such indices is not automatically a lower bound for M: different indices may produce the same pair (x,θ), and their factorial values are all 1. One must prove that enough distinct transitions remain, as well as bound the ordered-pair multiplicity μ. The unconditioned value-set theorem does not provide these additional facts.

The shared denominator and the gap must be controlled together.

For the blocks selected from first prime occurrences, the second inequality in (27) is

(29)log⁡C~p<log⁡(2p22p+1⋅∏j=p2p−1jq).

The first inequality is (25). Both must hold at the same arbitrarily large parameters.

The order of quantifiers is part of the proposed result: for every bound, one parameter must satisfy both inequalities. Two unrelated infinite subsequences would not prove the combined criterion.

For a prime r and an integer e≥1, let hr,e(p)=#{i∈Ip:re∣i!−1}. The exponent of r in the shared part is at least e exactly when hr,e(p)>1. Removing the exponent already present in Fp therefore gives

vr(C~p)=#{e≥1:hr,e+vr(Fp)(p)>1}.

When r∤Fp, this counts the positive exponents for which two or more denominators are divisible by re. The spacing bound hr,e(p)(e+1)≤2p+e−2 follows from the gap-product argument: if i<j are two such indices, then j<r and 0<re≤j!/i!−1<rj−i, so j−i≥e+1. All indices lie between 2 and 2p−1, giving the bound when hr,e(p)>0. It also holds when hr,e(p)=0. To estimate log⁡C~p, these exponent counts must be multiplied by log⁡r and summed over primes. A count ignoring the sizes of the primes does not give (29). For a prime q, Wilson reflection of an odd index n<q contributes a repeated divisor only when the reflected index is distinct and in the same block, with both indices at least 2. The endpoint n=q−2 reflects to 1 and contributes no such pair, as shown in §2.

Even for two admissible indices, the reflection identity is only modulo q [8]. For example, exact multiplication gives

609!≡1(mod9712),361!≡736019=1+758⋅971(mod9712),

where 361=971−609−1. Both indices lie in I306, but their factorial denominators have 971-exponents at least 2 and exactly 1, respectively. Thus both indices satisfy the congruence modulo q, but need not satisfy it modulo qe for e>1. Each level in the valuation sum requires its own congruence; neither distinctness nor a higher exponent may be inferred from the prime-level identity alone. Theorem 12 of Garaev, Luca and Shparlinski [7] is a multiplicity bound and supplies no lower bound for the distance to the next integer in (25).

For the first inequality alone, a joint mean bound

1|IX|∑p∈IX(log⁡C~p−log⁡ρpRp)<(1−ε)Xlog⁡X

on nonempty sets IX⊆[X,2X]∩N, with Rp>1, a fixed 0<ε<1, and arbitrarily large X, would suffice: at least one summand is at most the average, while the required upper bound is log⁡(Kp/((2p+1)Fp))=plog⁡p+O(p), uniformly for p∈[X,2X]. Two separate estimates on unrelated sets of indices would not give this joint comparison.

What two gap-product tests do not show.

For integers h≥1 and x≥0, let Fh(X)=∏j=1h(X+j). Suppose x+h<q for a prime q. If x! and (x+h)! both equal 1 modulo q and the harmonic sums through these indices are equal modulo q, then Fh(x)=1 and Fh′(x)=Fh(x)∑j=1h(x+j)−1=0 in that field. An exact calculation gives

gcd(F12−1,F12′)=(X−4626)(X−7848)in F12487[X].

Thus a uniform bound of one repeated root is false. But 4626!≡442 and 7848!≡6300 modulo 12487, so neither root also satisfies x!≡1. Counting repeated roots of Fh−1 therefore does not count the indices required by the factorial problem. These are finite modular calculations, not a bound on the number of indices satisfying both congruences.

Stewart’s gap-product estimate [8] requires the two products to be unequal. He supplies that step separately: in the proof of his bound (14), a prime-distribution argument separates the products when the common divisor is large; the smaller-divisor case is immediate [8].

The non-power theorem gives a different argument when the gaps are adjacent. For integers 2≤i<j<k, equality j!/i!=k!/j! would make k!/i! a square, contrary to that theorem [31]. Selecting two short gaps with nearly equal lengths need not preserve adjacency, so this observation alone cannot replace Stewart’s nonvanishing argument. The wider equal-product problem has a separate history [29]; the doubled-length theorem discussed by Nguyen Xuan Tho [30] is another special case, not arbitrary near-equal-length nonvanishing. No stronger lcm exponent is inferred from these local observations.

Prime powers already seen at earlier indices.

Recall that Δn=PnΔ/QnΔ is the positive gap in lowest terms. Define

An=∏r prime, r∣n!−1r∣k!−1 for some 2≤k<nvr(QnΔ)<vr(n!−1)rvr(n!−1).

Thus the product includes the full power of a prime in n!−1 when that prime has appeared before but has smaller exponent in QnΔ. The recurrence for the gaps is

Δn+1=nΔn−bn−1n!−1.

For a prime r in the product, put e=vr(n!−1). The reduced denominator of nΔn−bn has r-exponent at most vr(QnΔ)<e, whereas 1/(n!−1) has exponent e. Clearing denominators leaves exactly one numerator term nonzero modulo r, so the reduced denominator of Δn+1 has exponent exactly e. Hence An∣Qn+1Δ and, when An>1, Pn+1Δ≢0(modAn). The earlier occurrence of r is not used in this divisibility argument; it restricts the product to the primes proposed for this approach. A sufficient quantitative estimate, at arbitrarily large indices n, with m=n+1 and d=An>1, is

(30)((m+2)m!−2)QmΔ≤m2(m!−1)(PmΔmodd).

It implies the upper alternative of (21). A nonzero least representative can equal one even when the modulus is large. Neither a lower bound for log⁡An nor infinitely many increases in a previously occurring prime’s exponent therefore proves (30): the least representative itself must satisfy the displayed inequality. Repeated prime squares can first occur before the Wilson index r−2: for example, 971 divides 361!−1 and its square first divides 609!−1, with 609<971−2. Such an example rules out confinement to that last index but supplies no unbounded family satisfying the required inequality.

The upper alternative in the finite test.

It would suffice for (21) to hold at arbitrarily large indices; an unbounded set of prime indices would suffice as well. The upper alternative is exactly ((m+2)m!−2)QmΔ≤m2(m!−1)PmΔ. Replacing PmΔ by its least nonnegative residue modulo a specified divisor of QmΔ gives a stronger sufficient test. No bound forcing this inequality at arbitrarily large indices is established here. The full recurrence for the actual scaled partial sums, including its initial value, is fixed; a rational example for a weakened recurrence does not satisfy that full specification (§B.7). Nor does a zero canonical digit assert membership in the specified real interval.

Conditions at twice a prime.

For an odd prime p, the carry range and Z2p=2pZ2p−1+1−b2p give

(31)p2∣Z2p⟺{b2p=1 and p∣Z2p−1, orb2p=1+p and p∣2Z2p−1−1.

Reduction modulo p forces b2p≡1(modp), and the two displayed values are the only possibilities in [−1,2p−1]; division by p gives the predecessor conditions. For rational S=a/q and p>q, the tail estimate gives Z2p=(2p)!S, which is divisible by p2. A proof through p2∤Z2p must therefore exclude both displayed alternatives at the same prime, for arbitrarily large primes. That is a requirement of this particular test, not of every irrationality argument. Rationality also forces b2p=1 and Z2p−1=(2p−1)!a/q, with p∣Z2p−1: here q∣(p−1)!, and the remaining factorial factors include p. Thus either b2p≠1 at arbitrarily large odd primes or p∤Z2p−1 at arbitrarily large odd primes would already contradict rationality. The first is a special case of Theorem 4.1; the second uses the same factorial-clearing argument. Neither unbounded occurrence is proved here.

Criteria from the literature that do not apply.

Duverney’s Theorem 3.1 assumes, among other conditions, quadratic growth cun2≤un+1≤c′un2 for positive constants c,c′. For un=n!−1, however, un+1/un2→0. His Corollary 3.2 assumes convergence of the signed series ∑n(un+1/un2−1) in (3.6). Here its terms tend to −1, so that hypothesis fails as well [10]. For a single denominator sequence the rapid-growth criterion goes back to Erdős [3]; Barreto, Kang, Kim, Kovač and Zhang treat products of consecutive denominators and weighted extensions [12]. For an=n!−1, log⁡(n!−1)=O(nlog⁡n) makes (n!−1)1/ψn→1 for every fixed ψ>1, so neither Erdős’s hypothesis lim supnan1/2n=∞ nor the growth hypotheses of their Theorems 2 and 3 hold. These criteria do not decide S; the useful content is the target these papers identify: a subsequence with clearing integers small enough for the corresponding remainders, or enough exact cancellation in those integers. Theorem 3.4 measures the cost of clearing all summands before reduction; it does not bound cancellation in their sum. Theorem 2.1 tests whether a maximal prime power survives reduction. Neither supplies the upper bound for a clearing integer, relative to its remainder, that these arguments need.

No irrationality conclusion is obtained here. The finite conclusion is that every rational representation S=a/q, q>0, must satisfy both q∤299999! and q≥239990>1012038. The criteria above require a condition at arbitrarily large indices. Extending either finite computation alone does not establish it.

Sources and evidence

The margin marks and the table below identify the recorded Lean statements and supporting formulas at their original revisions. A link to a nearby statement does not certify a further deduction in the prose. The two partial-sum cancellations, the index-52 example and the continued-fraction enclosure use the finite integer calculations described above. The carry computation through 300000 uses a separate exact-interval computation outside Lean, described in §6.

Statement

Lean statements of its steps
Statement Lean statements of its steps (continued)

Theorem 2.1

maximal-power survival, with the residue formula and the denominator form
Proposition 2.2 first prime occurrences at arbitrarily large indices; the reflected divisibility
Theorem 4.1 failures of divisibility for the next integer; denominator exclusions: the size bound and the factorial nondivisibility
Proposition 4.2 factorial residues and rationality
Proposition 4.3 the equivalent fractional-part inequality; the finite window, the irrationality implication, the classification (22), the radius bound
Theorem 4.4 the criterion for the shifted series; the escape form, the member t=−1, and the irrationality of e
Theorem 5.1 the prime-parameter specialization; the margin mark gives the natural-parameter statement; the coprimality (24)
Equation (27) the two simultaneous inequalities, the floor identity, the irrationality implication
Removal of a fixed denominator factor a fixed factor eventually divides the factorial
Weight identities used in Theorem 1.1 and (1) integrality of the weights, the denominator-times-weight identity
Theorem 3.5 finite radius bound; the sequence form, the little-o form
Changing one weighted sum Vp adjustment at sufficiently large indices, the residual identity
Equation (31) doubled-prime criterion
Canonical factorial digits termination equivalence
Amplification modulus divisibility, nonvanishing

Attribution. Wilson’s theorem and the Wilson reflection identity are classical, the latter recorded by Stewart [8]. The factorial-digit termination criterion goes back to Cantor [5]; Koepf and Schmersau prove its irrationality direction for digits that are not eventually maximal [9], and Galambos treats rationality criteria for Cantor series [6]. Here the identity C=S−e+2 gives the eventual digit value m−2 and extends to the shifted family. The comparison with Hančl and Tijdeman [11] concerns factorial scaling: their lemma makes a scaled tail integral, whereas the proof in §4 identifies the next integer above a generally nonintegral scaled partial sum. The multiplicity bound is Garaev, Luca and Shparlinski’s [7], and the lcm deduction from it is not theirs. The divisibility in Lemma 3.2 is the case P=−1 of the relation used in the proof of Lemma 5 of Luca and Shparlinski [13] and at display (2.5) of Lai [14]. The nonvanishing cutoff for polynomial shifts is Lemma 3 of [13], restated with the bound n!+P(n)>1 in [14]. The survival criterion is a specialisation of Louwsma and Martino’s valuation formula [4]. The single-denominator growth criterion is Erdős’s [3], and the continued-fraction identities are those of [17]. The deductions from Wilson’s theorem, the conditional finite support bound and Theorem 3.4 are proved above. No further priority claim is inferred from this comparison.

Guide to the formal sources

Lean source links (115)

This concordance lists the floor, support and divisibility identities used alongside the principal results. All links retain their original source revisions. The last subsection contains the historical links from the short paper, with their targets unchanged.

Source links from the short paper

These links collect the declarations formerly listed in the short paper’s source appendix, at the same historical revisions and paths. They supplement the concordance above. No fresh compilation or independent review is asserted by a link, and a theorem about a conditional family leaves its hypotheses to be established.

The remaining links from the short paper use revision 99f4bf47422a; the links in the preceding concordance use this record’s own revision. The following numbers link to source lines within each named module, retaining both revisions where they differ.

AdjacentUnitCarryWindow.lean
ChannelBreakpointRigidity.lean
ChannelIntegralCongruence.lean
DivisorFactorialCentre.lean
EndpointWeightedPrivateSupport.lean

4766.

FactorialChannelCertificate.lean
FactorialZeroPlateau.lean
PrimeUnitTranslator.lean

1559.

PrimeZeroBranch.lean

6099.

The hypotheses about short intervals, least prime factors, and the nonzero real remainders are stated in the mathematical sections above. They remain hypotheses when a declaration link is followed.

Earlier links for the consolidated criteria.

The four criteria treated here retain the short paper’s earlier source and comparison links below. These identify historical declarations and recorded comparisons; they are not new executions.

Further deductions and limitations

The following sections give additional digit identities, denominator calculations and finite examples. For each unsuccessful argument we specify the missing estimate or the counterexample to its inference. The resulting limitations apply to those arguments, without ruling out other uses of the same identities.

Factorial digits and their remainders

We derive the factorial-digit identities used in §4 directly from their remainders. For a real x, start with θ0={x} and, for m≥1, put

am=⌊mθm−1⌋,θm=mθm−1−am,

so that θm={m!x} for m≥1. The recorded formal identities give the floor formula and the digit bounds 0≤am<m, the recurrence θm+1=(m+1)θm−am+1, the finite telescoping expansion

x=⌊x⌋+∑m=2Namm!+θNN!,

and the rule that a zero remainder at one index forces every later digit to vanish. The rational direction is also checked: if a∈Z, q∈N, and 0<q≤n, then ⌊n!a/q⌋=(n!/q)a and the canonical digit at radix n+1 vanishes, so every rational input has an eventually zero expansion. The converse holds for every real input, and gives Cantor’s termination criterion [5]

x∈Q⟺am(x)=0 for all large m.

If all digits after index N≥1 vanish, the recurrence gives θN+k≥(k+1)θN; since every remainder is below one, θN=0 and x=⌊N!x⌋/N!.

For x=S, a zero canonical digit means mθm−1<1. The stronger condition mθm−1<Em, where Em=m!(S−Hm), is the lower-interval event used in §4. The reported interval data contain zero digits at m=5 and m=23, but no occurrence of that stronger condition through m=100000.

There is also a useful identity for an abstract rational sequence Fm and integer sequence km. Suppose Fm=mFm−1+1+εm−km, and define

δm=⌈Fm⌉−Fm,qm=m⌈Fm−1⌉+1−km−⌈Fm⌉.

Then 0≤δm<1, and the exact identities are

qm=⌊mδm−1−εm⌋,δm=mδm−1−εm−qm.

The Lean sources prove these identities under the displayed recurrence. They do apply to the actual partial sums: take Fm=m!Hm, km=0 and εm=1/(m!−1) for m≥3. The recurrence then follows from Hm=Hm−1+1/(m!−1), starting with F2=2. The ceiling code qm need not equal the carry bm, because Zm is the least integer strictly above Fm, not always its ceiling. Precisely, Zm=⌈Fm⌉+1{Fm∈Z}, so

bm=qm+m1{Fm−1∈Z}−1{Fm∈Z}.

For example, F2=2 and F3=36/5 give q3=−1 but b3=2. The codes agree whenever both scaled prefixes are nonintegral. Once F2=2 and εm=1/(m!−1) are fixed, the recurrence with km=0 determines every Fm uniquely by induction. Thus matching the recurrence and its initial value is stronger than satisfying selected congruences or carry bounds. This specialization is an ordinary deduction from the partial-sum recurrence, not a new Lean check. It does not show that the rational gaps admit a finite-state description or establish the unbounded non-unit carries needed for irrationality.

How the classical criteria apply here

Let sn<a be partial sums tending to a. Koepf and Schmersau show that ⌊nsn⌋=⌊na⌋ for all sufficiently large n forces a to be irrational [9]. The strict inequality matters: if a=A/B were rational, then every sufficiently large multiple n of B would give ⌊nsn⌋<na=⌊na⌋. Their rational-term version uses a positive integer multiplier pn and obtains the floor equality from npnsn∈N0 and the strict tail bound a−sn<1/(npn) [9]. Indeed, integrality makes nsn a multiple of 1/pn, so its gap to the least strictly larger integer is at least 1/pn; the increase n(a−sn)<1/pn cannot cross that integer. For the choice they record after (2.1), the least common multiple of the reduced summand denominators, here pn=lcm⁡{k!−1:2≤k≤n}, the last two denominators already obstruct the tail bound. As proved in §3, they are coprime for n≥3, so

pn≥(n!−1)((n−1)!−1).

For n≥4 this gives npn>(n+1)!−1. Hence the first omitted summand 1/((n+1)!−1) already exceeds 1/(npn), and this pn cannot satisfy their tail hypothesis. This rules out the displayed choice of pn, not the criterion itself. In fact, the least positive integer pn for which npnHn is integral is exactly

pn=den⁡(nHn)=den⁡(Hn)gcd(n,den⁡(Hn)).

Indeed, writing Hn in lowest terms shows that its denominator must divide npn. With this least choice the required inequality is

0<S−Hn<gcd(n,den⁡(Hn))nden⁡(Hn).

Thus reduction of the partial sum, and the common factor of its reduced denominator with n, determine the best scale available in this criterion. The eventual inequality above remains unproved here; failure of a larger clearing multiplier does not establish its failure.

Duverney’s Theorem 3.1 includes the quadratic growth assumption cun2≤un+1≤c′un2 for positive constants c,c′, while un+1/un2→0 for un=n!−1 [10]. His Corollary 3.2, which allows signs an∈{−1,1}, assumes convergence of the signed series ∑n(un+1/un2−1) in (3.6) [10]. Its terms tend to −1 here, so the series does not converge. Lean checks the ratio limit and also the divergence of the absolute-value series; the signed condition fails already by the term test. To see the ratio limit directly, divide numerator and denominator by (n!)2: the numerator tends to zero and the denominator (1−1/n!)2 tends to one.

For a strictly increasing sequence (nk) of positive integers, Erdős proved irrationality of ∑k1/nk under the following conditions, with a fixed ε>0 [3]:

lim supk→∞nk1/2k=∞,nk>k1+εeventually.

Barreto, Kang, Kim, Kovač and Zhang treat products of consecutive denominators and weighted extensions [12]. For an=n!−1 and every fixed ψ>1,

0≤log⁡(n!−1)ψn≤n2ψn⟶0,(n!−1)1/ψn⟶1.

Thus neither Erdős’s limsup hypothesis, which uses ψ=2, nor the growth hypotheses of their Theorems 2 and 3 hold. The Lean proof of this limit treats the special case ψ=2; the estimate above also covers the other bases greater than one used in those theorems. Their proof of Theorem 3 uses the classical criterion, which they trace to Fourier’s proof that e is irrational, that a rational sum of nonnegative rationals with infinitely many positive terms admits no integers DN clearing the denominators of the partial sums with lim infNDNrN=0, where rN is the tail after N terms [12]; their Proposition 12 produces such integers under the hypotheses of that theorem [12]. For the present series, even the full summand lcm LN makes LN(S−HN) tend to infinity, as shown in §3. A smaller integer cannot clear every summand. It can nevertheless clear the partial sum after addition: the least positive such integer is den⁡(HN), and every other one is a multiple of it. Thus this criterion would require information about cancellation in the reduced partial sum, not merely an improvement from the product to the lcm. More precisely, if S=a/q, then

qden⁡(HN)(S−HN)∈Z>0,den⁡(HN)(S−HN)≥1q.

Thus lim infNden⁡(HN)(S−HN)=0 is a sufficient target for this rationality criterion; the preceding common-denominator estimates do not establish it.

Dividing the recurrence Zm=mZm−1+1−bm by m! and telescoping gives the exact finite identity

ZMM!=Z22!+∑m=3M1−bmm!,

so the carry defects 1−bm are integer coefficients of a factorial series. Hančl and Tijdeman classify the rational sums with polynomial coefficients [11]; their denominator is the cumulative linear product ∏n≤N(an+b), and the individual number N!−1 does not occur. In the factorial case a=1, b=0, the criterion of Oppenheim that they reproduce [11] says the following: for integer coefficients cn with |cn|<n eventually and lim infn|cn|/n=0, the sum ∑ncn/n! is rational exactly when cn=0 eventually. This formulation concerns integer coefficients and factorial denominators, not reciprocal terms with denominator n!−1. However, their introduction records a version of Oppenheim’s result without the liminf assumption: under |cm|<m−1 eventually, rationality is equivalent to eventual vanishing [11]. Changing finitely many terms adds a rational number, so eventual hypotheses suffice. Here 2−m≤cm=1−bm≤2 gives |cm|<m−1 for m≥4; the classical result therefore applies directly.

For completeness, the following elementary proof specialises that result to these one-sided bounds. For N≥2, telescoping ∑m>N(m−1)/m!=1/N! gives

−1<N!∑m>Ncmm!≤2N!∑m>N1m!<2N.

The lower bound follows from cm≥2−m and N!∑m>N(m−2)/m!=1−N!∑m>N1/m!<1. The upper bound follows by comparison with ∑j≥1(N+1)−j=1/N. If ∑m≥3cm/m! is rational, the scaled tails are integers for all sufficiently large N. They must then be zero, and subtracting consecutive tails gives cm=0 eventually. The converse is immediate. Since ZM/M!→S, the displayed finite identity recovers Theorem 4.1. This is a self-contained specialisation of the classical rationality criterion, not a strengthening of it. It proves an equivalence, not the nonvanishing needed for irrationality.

A superseded deduction

Factorial-congruence multiplicities give another lower bound for LN. Its exponent is weaker than that of the terminal-block argument in §3. We include the deduction to specify what follows from this external input. For an odd prime r, let mr count the indices 2≤n≤N with r∣n!−1; such an index satisfies n<r. The factorial-congruence multiplicity estimate of Garaev, Luca and Shparlinski [7], applied to the residue a=1 on 1≤n≤min(N,r−1), gives mr≪N2/3; the prime 2 divides none of these factors. The sum of the exponents at a prime is at most the number of occurrences times the largest exponent, which is its exponent in LN. Thus

∑n=2Nlog⁡(n!−1)=∑r∑n=2Nvr(n!−1)log⁡r≤(maxrmr)log⁡LN≪N2/3log⁡LN,

and Stirling’s formula makes the left side ≍N2log⁡N, whence log⁡LN≫N4/3log⁡N. Theorem 3.4 supersedes this: its exponent is 3/2, its constant is explicit, and it does not use the external multiplicity theorem. The deduction is retained because it is the only place a multiplicity bound enters the record.

Dividing all coefficients by their gcd does not evade the lcm restriction: the resulting integer vector still has V2=⋯=VD=0, so its value of M is still divisible by LD. For the nonzero cofactor vector, §1.5 justifies normalisation by identifying the primitive progression vector. It supplies no proof that the remainder is nonintegral. No formalised primitive cofactor construction is claimed here.

Rational grid points and first crossings

Fix q>0. To compare a partial sum H with a rational number of denominator q, consider the grid q−1Z. Its next point above H is (⌊qH⌋+1)/q. It is at most S exactly when ⌊qH⌋+1≤qS. Suppose H<G≤S, the number n!G is an integer, and n!(S−H)<1. Then

n!H<n!G≤n!S<n!H+1.

Thus n!G is both the least integer strictly above n!H and the floor of n!S. In particular, for n≥2, suppose a single rational level G satisfies Hn+1<G≤S and n!G∈Z. Since Hn<Hn+1, the tail bounds at both indices give

⌊(n+1)!S⌋=(n+1)!G=(n+1)⌊n!S⌋.

The next canonical digit is therefore zero. The same level must persist at both indices; choosing an unrelated rational level at each index would not imply this equality.

Set H1=0, the empty partial sum. Now let τ≥2 be a first crossing of a rational level G, so that Hτ−1<G≤Hτ, and write G−Hτ−1=a/v with positive integers a,v. The newly added summand gives

0<av≤1τ!−1,v≥a(τ!−1)≥τ!−1.

The scaled overshoot is

0≤τ!(Hτ−G)=1+1τ!−1−τ!av<2.

Consequently −⌊τ!(Hτ−G)⌋ is 0 or −1. Here the minus sign is outside the floor; taking the floor of the negative overshoot would give a different integer at nonintegral overshoots. The displayed integer is −1 precisely when the overshoot is at least 1, or equivalently when

av≤1τ!(τ!−1).

In that case v≥τ!(τ!−1). This extra small-gap condition is not asserted for every crossing. Nor is −⌊τ!(Hτ−G)⌋ automatically the carry bτ defined from Zτ: here G is an arbitrary rational level. Neither denominator bound requires a/v to be in lowest terms.

There is also a direct obstruction at prime indices. With Δm as defined in §4, the formal source proves for m≥3 the criterion

m∣Zm⟺1+1m!−1<mΔm≤2+1m!−1,

and if S=a/q with q>0, then p∣Zp for every prime p>q, so one exact missed prime p gives q≥p. Exact kernel reduction gives 11∤Z11, hence q≥11; exact rational normalisation gives 60∤Z60, 64∤Z64 and 67∤Z67, and the prime index 67 gives the checked bound

S=aq, q>0⟹q≥67.

The exact-interval census of the short note replaces 67 by 300000.

The finite geometric-series identity gives another exact decomposition. For a real x≠0,1 and an integer K≥0,

1x−1=∑j=1K1xj+1xK(x−1).

When x=k! and a chosen factorial scale is divisible by (k!)K, the scaled finite sum is integral and only the last term retains the factor k!−1 in its denominator. The identity isolates one remaining fraction before exact bounding, and it supplies no cofinal family of nonzero residuals.

The earlier record rejects a proposed divisibility strengthening of the first-crossing bound and reports examples at m=52 and m=591. Its notation for the rational quantity in that claim was not defined, so those reports are not used as a verified counterexample here. The established conclusion remains the lower bound on the denominator’s size, not a specified factor dividing it.

There is a second, more arithmetic mechanism at doubled prime indices, stated in (31). The formal theorem is not restricted to individually computed indices; the divisibility criterion holds for the actual partial sums at every odd prime.

Exact identities for shared prime powers

Let I⊆{2,3,…} be finite. We first compute the prime-power levels common to two summand denominators. We then compare those exponents with the ones already supplied by a factorial factor.

Write D(I)=lcmi<j, i,j∈I⁡gcd(di,dj), with D(I)=1 for an empty or singleton family. For each prime r, its valuation is the second-largest valuation among the di, counting missing values as zero; the denominator lcm has the largest valuation. This proves D(I)lcmi∈I⁡di∣∏i∈Idi prime by prime.

The integer D(I) can be updated exactly when one index is added. For the positive factorial-gap denominators, adjoining a new index a≥2, a∉I, gives

D(I∪{a})=lcm(D(I), gcd(da,lcmj∈I⁡dj)),

since finite-family gcd and lcm distributivity collapses the lcm of all pairwise gcds against da to a single gcd. For a fixed positive integer F, take the lcm with F on both sides:

lcm⁡(F,D(I∪{a}))=lcm⁡(F,D(I),gcd(da,lcmj∈I⁡dj)).

Thus one need only retain the denominator lcm and lcm⁡(F,D(I)) at each step. The prescribed factor F is not an extra summand denominator.

There is also an exact bound by the product of the denominators divided by their least common multiple. For a positive integer F, define C~(I)=lcm⁡(F,D(I))/F=D(I)/gcd(F,D(I)). Since C~(I) divides D(I), the preceding divisibility gives

C~(I)lcmj∈I⁡dj∣∏j∈Idj,C~(I)≤∏j∈Idjlcmj∈I⁡dj.

For the factorial block this specialises to

C~p≤∏n∈Ip(n!−1)lcmn∈Ip⁡(n!−1),

an upper bound for the shared factor in terms of the denominator product and lcm. The bound alone does not establish (29).

Dividing by Fp subtracts its prime exponents; it does not remove every prime that occurs in Fp. With Fp=(p−1)!, Dp=D(Ip) and Cp=lcm⁡(Fp,Dp), as in the main record,

C~p=CpFp=Dpgcd(Fp,Dp),vr(C~p)=max{0,vr(Dp)−vr(Fp)}.

For e>0, re∣C~p exactly when Dp is divisible by re+vr(Fp), which in the block forces two distinct denominators to be divisible by that higher power, by the second-largest-valuation formula. If i<j are two such indices and f=vr(Fp), then re+f∣j!/i!−1. Since r∣j!−1 forces j<r, 0<j!/i!−1<rj−i, hence e+f<j−i. For a prime r∣C~p, this bounds the surviving valuation by vr(C~p)+vr(Fp)<r. Since vr(C~p)≥1 and vr(Fp)≥⌊(p−1)/r⌋, we obtain ⌊(p−1)/r⌋≤r−2, hence p−1<r(r−1)<r2. Consequently C~p is coprime to k! whenever k(k−1)≤p−1: a prime r≤k dividing both would give p−1<r(r−1)≤k(k−1)≤p−1. This excludes the small primes but does not bound the product of the remaining prime powers.

If a prime r divides a denominator n!−1 with n≥p, then r∤Fp and re∣C~p⟺re∣Cp for every e>0. This has an exact incidence-count form,

re∣C~p⟺1<#{i∈Ip:re∣i!−1},

so a bound of at most one such index implies vr(C~p)<e. Counting the exponents gives

vr(C~p)=#{e∈[1,r−1]:1<#{i∈Ip:re∣i!−1}},

and for a prime r>2p−1 the range can be restricted to e∈[1,2p−4]. The same equivalence holds under the exact condition r∤Fp in place of r>2p−1, hence in particular for every prime r≥p. At e=2, an upper bound of one on the number of indices gives vr(C~p)≤1.

A surviving prime power also forces two indices to be far apart. If r is prime, e>0 and re∣C~p, there are i<j in Ip with re+vr(Fp)∣i!−1, re+vr(Fp)∣j!−1 and re+vr(Fp)≤jj−i; consequently (2p−1)d<re+vr(Fp) forces d<j−i. The required spacing follows directly: any two indices i<j with re∣i!−1 and re∣j!−1 satisfy e<j−i, because r∣j!−1 already forces j<r and the inequality for j!/i!−1 then forces the strict separation. More generally, for any two indices 2≤i<j and any prime power re with e>0,

re∣i!−1,re∣j!−1⟹re≤jj−i,

so successive solutions of i!≡1(modre) are at least e+1 apart, and

(e+1)#{i∈Ip:re∣i!−1} ≤ 2p+e−2.

The resulting bound on the exponent is

r∣C~p⟹vr(C~p)+vr(Fp)<2p−3

for every prime r and p≥2. Thus the exponent already present in Fp uses part of the same bound 2p−3. To apply (25), one still needs sufficiently strong bounds for the product of the shared prime powers and for the gap ρp/Rp.

The available squarefreeness evidence is finite. An exhaustive modular scan through r≤2,000,000 and n≤240 found four pairs (r,n) with r2∣n!−1 and no prime occurring in two such pairs. All 498,501 pairs 2≤a<b≤1000 have squarefree gcd(a!−1,b!−1). An aggregate scan through p=499 also reports a ratio below 0.374, but does not specify its logarithmic normalisation; that ratio is not used as a quantitative premise here. These are reported finite computations, not an asymptotic bound or a proof of squarefreeness at all indices.

A finite version of the argument for first prime occurrences compares the product of a chosen set of primes, each at least 5, with ∏2≤k≤B(k!−1); if the prime product is larger, at least one chosen prime q divides no k!−1 with 2≤k≤B, while Wilson’s theorem bounds its first occurrence by q−2. Wilson reflection limits what a linear-size divisor can be assumed to be: if n is odd, n<q, and q∣n!−1, then q∣(q−n−1)!−1. Suppose also that p≤n, both indices lie in Ip, and the reflected index is earlier, equivalently q<2n+1. The prime q then divides two denominators in the block. Since q>n≥p, it does not divide Fp=(p−1)!, so q∣C~p: removing the factorial factor does not remove this shared prime.

Stewart states that for every ε>0 there are infinitely many odd n whose least prime factor q of n!−1 satisfies

n<q<(145−18+ε)n,

estimate (9) of Theorem 1 being stated for n!+1 [8] and transferred to n!−1 in the text [8]. This controls q relative to the specified index n, not to its first occurrence. It does not by itself force a repeated divisor. For the minus sign, the bound is already met at all sufficiently large Wilson indices n=q−2. For every prime q≥5, Wilson gives q∣(q−2)!−1, and q is its least prime factor. Indeed, every prime divisor exceeds q−2, and q−1 is even; moreover q/(q−2)→1. The reflected index is then 1, outside the denominators of S.

For example, modulo 11 the only solution of k!≡1 with 2≤k≤9 is k=9, although 11/9<(145−1)/8. A small ratio q/n therefore supplies neither a second admissible index nor membership of a reflected index in the chosen block. The earlier reflection argument needs both block-membership hypotheses and distinct indices. Neither the required estimate at first occurrences nor a bound for the total shared part follows from the transferred least-prime-factor estimate alone.

For a selected prime q∣Rp, the factors 1 and q are coprime, and the associated moduli Rp and Rp/q have lcm Rp. Thus this specialisation requires no second selected prime.

The elementary fact behind the projection argument is as follows. Let Z,T,B be nonnegative integers, let R>0, and suppose Z≡T(modR). Every positive divisor Q of R with Z≤B<Q satisfies TmodQ=Z, since Z is already the least nonnegative representative. Consequently, unequal residues modulo two divisors exceeding B rule out such a Z. More generally, if TmodQ1≠TmodQ2, then min(Q1,Q2)≤T: otherwise both residues would equal T.

Finite vectors and exact numerical examples

The three depth-four vectors below all have M=1380. By Theorem 1.3, their remainders differ only by integers. A single finite calculation therefore tests all three:

255+45<242+1380∑d=581d!−1<255+56,2⋅13809!−1<1100.

The first remainder lies strictly between 255 and 256. Subtracting the corresponding integer offsets gives the other two intervals.

Vector ∑n>4zn k Interval for R
12K4+253U6−11U8 242 −1474 (255,256)
The five-term vector (5) −44 −1760 (−31,−30)
The norm minimizer on {2,…,8} −56 −1772 (−43,−42)

The second column is the offset from 1380(S−H4); the third is the integer k in R=1380S+k. Since 1380H4=1716, subtracting 1716 converts the second column into the third. The five-term remainder is 286 less than the first, and the norm-minimising remainder is another 12 less. All three exclude exactly the same divisors of 1380 by this argument. The coefficient and support improvements in §1.2 do not enlarge that set of exclusions; the larger finite computations in §6 already imply it. The exact offsets are reproduced by scripts/check_erdos68_channel_moment.py.

The finite-support vector λ=2e3−e4 has, by kernel check, V2(λ)=0, M=−12, V3(λ)=−2, V4(λ)=11, and Vd(λ)=−12 for every d≥5. Under the exact tail enclosure 1/119<∑d≥51/(d!−1)<1/50 its remainder lies strictly between −93/575 and −309/13685, so it is nonzero and has absolute value less than 1.

In the enlarged coefficient space of §1, exact integer computation verifies the vectors KD for every 2≤D≤12: weighted sums 2 through D vanish, M=LD, and the coefficients have gcd one. They lie outside the space of vectors supported on n≥2, since LD is odd whereas M is even for every such vector. At D=9,

L9=31540008254514077395,a9=[e1]K9=−3902884074990939115.

Since U11=T11 has u11=0, the vector K9−9553024718754U11 keeps the coordinate a9 and has coefficients with gcd one. By Theorem 1.3 its residual is L9(S−H9)−9553024718754, and exact rational arithmetic with the tail bound ∑n≥361/(n!−1)<2/(36!−1) places it strictly between 1353/100000 and 1354/100000. A different example is supported on n≥2: the vector c=(−40,55,−10,1) on the support (3,4,5,6) annihilates weighted sums 2 and 3, has M=600, and satisfies

0.09925341997208298<R(c)<0.09925341997208300,

which excludes denominators dividing 600.

A separate interval computation reports the stronger geometric statement that no lower-interval event mθm−1<Em occurs at any 3≤m≤100000; its executable and source digest are not available, so that classification remains external finite evidence. The carry computation through 300000 used in the short note is the exact GMP computation recorded in Section 6, with the separate calculation through 4000; it is computational evidence outside the Lean development.

Limits of the recurrences and clearing factors used here

Two distinctions are needed in using these identities: a recurrence may have solutions that violate the defining floor relations, and a fixed factor may eventually be absorbed by a growing factorial. The examples below keep these cases separate.

  • Without the defining floor relation, the carry recurrence and its range allow Zm/m! to be the constant 3/2. Taking Zm=3m!/2 for m≥2 and bm=1 for m≥3 gives Zm=mZm−1+1−bm and −1≤bm≤m−1, with the actual initial value Z2=3. But it gives Z3=9, whereas the actual partial sum H3=6/5 gives Z3=8. This example only disproves sufficiency of the stated recurrence, bounds and initial integer. It fails the floor definition. No claim is made that it satisfies the additional prime-index identities. The rational recurrence for Fm=m!Hm in §B.1 determines the actual partial sums uniquely.

  • The arguments considered with Wilson quotients, harmonic sums, p-adic gamma identities, and factorial residues still need a real gap estimate and the required modular divisibility at the same indices. No such unbounded family is obtained here. This records the missing step in these arguments, not an impossibility theorem for the use of those identities.

  • For the genus-zero product E(z)=∏n≥2(1−z/n!), local uniform convergence and logarithmic differentiation give −E′(1)/E(1)=S. Termwise clearing by ∏n=2N(n!−1) cannot give a positive remainder tending to zero: this product is at least LN, and LN(S−HN)→∞ by §3. The conclusion concerns this clearing factor alone; Hermite–Padé systems with other denominators are not ruled out.

  • Changing a coefficient vector without changing its value of M changes the remainder by an integer only. Also, cancelling the first D−1 weighted sums forces LD to divide the value of M; factorial divisibility does not remove that constraint.

  • A fixed pair of denominator indices cannot make the projection argument work at arbitrarily large parameters, because their factors eventually divide the factorial being removed.

An elementary criterion for any real number

The fractional-part condition in (20) is an instance of Cantor’s termination criterion. For any real x and positive sequence cm satisfying cm≤1/m eventually, we have

x∉Q⟺{(m−1)!x}≥cm for arbitrarily large m.

For rational x the fractional parts vanish eventually. Conversely, if {(m−1)!x}<cm≤1/m for all large m, then m{(m−1)!x}<1, so every sufficiently late canonical digit is zero. Cantor’s criterion then makes x rational. Equality cm=1/m is allowed: it is failure of the displayed weak inequality that makes m{(m−1)!x} strictly less than one. For x=S, the choice cm=Em/m gives exactly (20); indeed 0<Em<1 ensures the required threshold bound. The finite test (21) is a separate sufficient condition at one index. Its cofinal converse uses the carry theorem, as shown in §4, not substitution into this general criterion. Neither argument proves that S meets the test at arbitrarily large indices.

The threshold 1/m cannot be replaced by λ/m for a fixed λ>1. For x=e and m≥3, the factorial tail gives

m{(m−1)!e}=1+m!∑n>m1n!<1+1m<λeventually.

Thus the irrational number e would never meet the enlarged threshold at sufficiently large indices. Positivity is also essential: with cm=0, rational numbers would meet the threshold eventually. These examples explain the two assumptions without imposing any equidistribution hypothesis.

Relations among the criteria

The classification determines all cancellation vectors and attainable multipliers. The progression construction adds eventual divisibility by each fixed q and a vanishing tail bound for the same vector. Its primitive cofactor construction gives that vector up to sign, not another family. At fixed M, coefficient changes still alter only the integer in R(c)=M(c)S+k.

Common-denominator growth measures the cost of clearing, while prime-power survival describes reduction. Neither supplies the missing gap comparison. The carry, digit and interval criteria characterize irrationality but leave their required unbounded departures unproved. The finite exclusions restrict divisibility and size, respectively; they establish no unbounded assertion.

The remaining target for this method is therefore precise: for every q>0, find a nonintegral remainder whose own multiplier satisfies q∣M. The algebraic classification supplies choices, not this arithmetic conclusion.

Acknowledgements

An AI research pass supplied by Will Cook derived the depth-four dual congruence and shorter vectors from the weighted linear forms and the classification of their possible coefficients M developed earlier here. OpenAI Codex checked and integrated them; a separate AI pass reviewed the proof.

References

  1. Paul Erdős. On the irrationality of certain series: problems and results. In Alan Baker (ed.), New Advances in Transcendence Theory, Cambridge University Press (1988), pp. 102–109.

  2. Thomas F. Bloom. Erdős Problem #68. Online resource (2026). Historical access: 28 July 2026; present-page status not reverified.

  3. Paul Erdős. Some problems and results on the irrationality of the sum of infinite series. Journal of Mathematical Sciences 10 (1975), 1–7.

  4. Joel Louwsma and Joseph Martino. Rational numbers with odd greedy expansion of fixed length. Preprint (2023). arXiv:2309.07280.

  5. Georg Cantor. Über die einfachen Zahlensysteme. Zeitschrift für Mathematik und Physik 14 (1869), 121–128.

  6. János Galambos. Representations of Real Numbers by Infinite Series. Lecture Notes in Mathematics 502, Springer (1976).

  7. Moubariz Z. Garaev, Florian Luca and Igor E. Shparlinski. Character sums and congruences with n!. Transactions of the American Mathematical Society 356 (12) (2004), 5089–5102. arXiv:math/0403422.

  8. Cameron L. Stewart. On the greatest and least prime factors of n!+1, II. Publicationes Mathematicae Debrecen 65 (3–4) (2004), 461–480.

  9. Wolfram Koepf and Dieter Schmersau. Irrationality of certain infinite series II. Analysis 31 (2011), 117–124.

  10. Daniel Duverney. Irrationality of fast converging series of rational numbers. Journal of Mathematical Sciences, the University of Tokyo 8 (2001), 275–316.

  11. Jaroslav Hančl and Robert Tijdeman. On the irrationality of factorial series. Acta Arithmetica 118 (4) (2005), 383–401.

  12. Kevin Barreto, Jiwon Kang, Sang-hyun Kim, Vjekoslav Kovač and Shengtong Zhang. Irrationality of rapidly converging series: a problem of Erdős and Graham. Preprint (2026). arXiv:2601.21442.

  13. Florian Luca and Igor E. Shparlinski. Prime divisors of shifted factorials. Bulletin of the London Mathematical Society 37 (6) (2005), 809–817.

  14. Li Lai. On the largest prime divisor of n!+1. Bulletin of the Australian Mathematical Society 113 (3) (2026), 390–403. arXiv:2103.14894.

  15. Li Lai, Cezar Lupu and Johannes Sprang. On the irrationality of certain p-adic zeta values. Research in the Mathematical Sciences 12 (4) (2025), article 77. arXiv:2505.23088.

  16. Li Lai. On the irrationality of certain 2-adic zeta values. International Journal of Number Theory 21 (1) (2025), 207–235. arXiv:2304.00816.

  17. NIST Digital Library of Mathematical Functions. Continued fractions: convergents. Online resource (2026).

  18. Paul Erdős and Cameron L. Stewart. On the greatest and least prime factors of n!+1. Journal of the London Mathematical Society 13 (3) (1976), 513–519.

  19. Florian Luca and Igor E. Shparlinski. On the largest prime factor of n!+2n−1. Journal de Théorie des Nombres de Bordeaux 17 (3) (2005), 859–870.

  20. Moubariz Z. Garaev, Florian Luca and Igor E. Shparlinski. Distribution of harmonic sums and Bernoulli polynomials modulo a prime. Mathematische Zeitschrift 253 (2006), 855–865.

  21. Jonathan Sondow. A geometric proof that e is irrational and a new measure of its irrationality. American Mathematical Monthly 113 (7) (2006), 637–641. arXiv:0704.1282. Addendum: Amer. Math. Monthly 114 (2007), 659; reading copy v2.

  22. Xiyu Hu. Factorial residues modulo a prime: beyond the square-root bound. Preprint (2026). arXiv:2608.01781. Version 1, 3 August 2026.

  23. Moubariz Z. Garaev and Julio C. Pardo. Additive congruences with factorials modulo a prime. Proceedings of the American Mathematical Society 154 (10) (2026), 4179–4190. arXiv:2508.12127. Published online 1 July 2026; assigned October 2026 issue.

  24. Xiyu Hu. Lower bounds for some value sets over finite fields: incidence geometry and Bourgain’s group expansion theorem. Preprint (2026). arXiv:2609.05652. Version 1, 4 September 2026.

  25. William D. Banks, Florian Luca, Igor E. Shparlinski and Henning Stichtenoth. On the Value Set of n! Modulo a Prime. Turkish Journal of Mathematics 29 (2) (2005), 169–174.

  26. Oleksiy Klurman and Marc Munsch. Distribution of factorials modulo p. Journal de Théorie des Nombres de Bordeaux 29 (1) (2017), 169–177.

  27. Alexandr Grebennikov, Arsenii Sagdeev, Aliaksei Semchankau and Aliaksei Vasilevskii. On the sequence n! mod p. Revista Matemática Iberoamericana 40 (2) (2024), 637–648.

  28. Sophie Stevens and Frank de Zeeuw. An improved point-line incidence bound over arbitrary fields. Bulletin of the London Mathematical Society 49 (5) (2017), 842–858. arXiv:1609.06284.

  29. R. A. Macleod and I. Barrodale. On Equal Products of Consecutive Integers. Canadian Mathematical Bulletin 13 (2) (1970), 255–259.

  30. Nguyen Xuan Tho. On equal products of consecutive integers. Elemente der Mathematik 81 (3) (2026), 119–124. First published online 10 July 2025.

  31. Paul Erdős and John L. Selfridge. The product of consecutive integers is never a power. Illinois Journal of Mathematics 19 (2) (1975), 292–301.

  32. NIST Digital Library of Mathematical Functions. Vandermonde determinant, §1.3(ii), formula (1.3.13). Online resource, consulted 18 September 2026.

About this paper

Authorship and AI use. Will Cook built and directed the research infrastructure and maintains the public release. He reviewed claims when he could. AI agents did most of the research and drafting. Cook did not independently verify every claim.

Cite and contact. Cite this paper by its title, author and date above, with its PDF; cite the earlier sources it uses for a mathematical result. For software, use the release citation and give the commit used. Contact Will with questions or corrections.

Prefer the manuscript?

Equations typeset from the exact TeX. Rendered from the LaTeX manuscript source with file fingerprint sha256:0af0d85ec1beab46. Its file fingerprint matched the published source record, so the PDF above, the LaTeX, and this page are one manuscript. Redeploying the site regenerates this page from whatever the public repository holds at that moment.