Introduction
Erdős asked whether replacing 𝑛!
by 𝑛! −1 preserves the irrationality
of the reciprocal sum [1]. The question appears as Problem 68 in
Bloom’s catalogue [2].
Problem 1.1 (Erdős Problem #68). Is
𝑆=∑𝑛≥21𝑛!−1
irrational?
The elementary proof that 𝑒 is
irrational clears the partial sum by a factorial and bounds the
remaining positive tail. Here the first step fails: for 𝐻𝑚 =∑𝑚𝑛=2(𝑛! −1)−1, already
3!𝐻3 =36/5. We instead construct
forms 𝑀𝑆 +𝑘, with 𝑀,𝑘 ∈ℤ, by cancelling initial
weighted sums. If 𝑆 =𝑎/𝑞 and 𝑞 ∣𝑀, then 𝑀𝑆 +𝑘 is integral; proving it nonintegral
excludes that denominator.
At depth four, the weight congruences force 115 ∣𝑀. Restricting the coefficient
support to 𝑛 ≥2 strengthens this
to 1380 ∣𝑀, and every multiple
of 1380 occurs. We determine the
exact factor at every depth 𝐷 ≥2
and construct a vector attaining the least positive multiplier.
First allow a coefficient at index one. Theorem 2.1
constructs an integral basis by correcting adjacent factorial
differences at proper divisors to isolate each weighted sum. At fixed
𝑀, cancellation determines the
lower basis coefficients; removing index one leaves a single equation in
the higher coefficients. Theorem 3.1
reduces its infinite gcd to a calculation below 2𝐷2. Bézout coefficients then construct
an attaining vector. No assumption on the rationality of 𝑆 enters this classification.
Hančl and Tijdeman’s tail-integrality lemma treats factorial series
with integer coefficients [4]. Their factorial-scaled partial sums
are integers. Here the scaled partial sum need not be integral. Instead,
the weight congruences give R −𝑀𝑆 ∈ℤ (1); the
assumption 𝑆 =𝑎/𝑞 and 𝑞 ∣𝑀 then makes R integral. Koepf and
Schmersau’s factorial-digit criterion [3] is used in the companion’s
digit argument to give an equivalent eventual-pattern condition for
𝑆 −𝑒 +2. Failure of that pattern at
arbitrarily large indices also remains unproved.
Section 4
tests an attaining form; Section 5 gives the general
nonintegrality test.
An integral basis
We need integer weights congruent to 𝑛! modulo 𝑑! −1, with 𝑛𝑊𝑑,𝑛−1 =𝑊𝑑,𝑛 when 𝑛 ≥2 and 𝑑 ∤𝑛. Removing powers of 𝑑! preserves the residue; the floor
exponent gives the required jumps. Put
𝑊𝑑,𝑛=𝑛!(𝑑!)⌊𝑛/𝑑⌋(𝑑≥2, 𝑛≥1).
Writing 𝑛 =𝑘𝑑 +𝑟, 0 ≤𝑟 <𝑑, shows that 𝑊𝑑,𝑛
is 𝑟! times the multinomial
coefficient 𝑛!/((𝑑!)𝑘𝑟!). Thus it
is an integer, and 𝑛! =(𝑑!)𝑘𝑊𝑑,𝑛
gives 𝑊𝑑,𝑛 ≡𝑛!(mod𝑑! −1).
For example, 𝑊3,4 =4 and (4! −𝑊3,4)/(3! −1) =4.
For a finitely supported integer vector 𝜆 =(𝜆𝑛)𝑛≥1, define the
integer linear forms
𝑀(𝜆)=∑𝑛𝜆𝑛𝑛!,𝑉𝑑(𝜆)=∑𝑛𝜆𝑛𝑊𝑑,𝑛,
and let
R(𝜆)=∑𝑑≥2𝑉𝑑(𝜆)𝑑!−1.
Here 𝑛 indexes coefficients and 𝑑 indexes remainder summands. If the
support lies in 𝑛 ≤𝑁, 𝑁 ≥2, then 𝑉𝑑 =𝑀 for 𝑑 >𝑁. Thus the remainder converges
absolutely, and subtracting 𝑀𝑆
leaves only finitely many terms. The weight congruence makes each an
integer:
R(𝜆)−𝑀(𝜆)𝑆=𝑁∑𝑑=2𝑉𝑑(𝜆)−𝑀(𝜆)𝑑!−1∈ℤ.(1)
We solve 𝑉2 =⋯ =𝑉𝐷 =0, allowing either sign of
𝑀. Write 𝑒𝑛 for the unit vector at index 𝑛. Index one is auxiliary: it makes the
basis triangular but introduces no term 1/(1! −1), since the remainder still
starts at 𝑑 =2. We impose 𝜆1 =0 in Section 3.
The identity 𝑛(𝑛 −1)! −𝑛! =0
suggests starting with 𝑇𝑛 =𝑛𝑒𝑛−1 −𝑒𝑛. At 𝑛 =2, the vector 𝑈2 =2𝑒1 −𝑒2 has 𝑉2 =1 and every other weighted sum zero.
At 𝑛 =4, the difference 𝑇4 =4𝑒3 −𝑒4 has two nonzero weighted
sums. One correction suffices:
| Vector |
𝑀 |
𝑉2 |
𝑉3 |
𝑉4 |
| 𝑇4 |
0 |
6 |
0 |
23 |
| 6𝑈2 |
0 |
6 |
0 |
0 |
| 𝑇4 −6𝑈2 |
0 |
0 |
0 |
23 |
For 𝑑 >4 these weighted sums
equal 𝑀 =0. Thus 𝑈4 =𝑇4 −6𝑈2 isolates 𝑉4 =4! −1. In general, the unwanted sums
occur at proper divisors of 𝑛,
whose correcting vectors are already available in an inductive
construction.
Theorem 2.1 (an integer basis with prescribed weighted
sums). Set
𝑇𝑛=𝑛𝑒𝑛−1−𝑒𝑛,𝑈𝑛=𝑇𝑛−∑𝑑∣𝑛2≤𝑑<𝑛𝑊𝑑,𝑛𝑈𝑑(𝑛≥2).
Then
𝑀(𝑈𝑛)=0,𝑉𝑑(𝑈𝑛)=(𝑑!−1)𝟏𝑑=𝑛.
The vectors 𝑒1,𝑈2,𝑈3,… form an integral
basis. Every finite vector has the unique finite expansion
𝜆=𝑀(𝜆)𝑒1+∑𝑑≥2𝑉𝑑(𝜆)−𝑀(𝜆)𝑑!−1𝑈𝑑.(2)
Proof. We have 𝑀(𝑇𝑛) =0. If 𝑑 ∤𝑛, the floor exponents at 𝑛 −1 and 𝑛 agree, so 𝑛𝑊𝑑,𝑛−1 =𝑊𝑑,𝑛. If 𝑑 ∣𝑛, the exponent increases by one,
so 𝑛𝑊𝑑,𝑛−1 =𝑑!𝑊𝑑,𝑛. Hence
𝑉𝑑(𝑇𝑛)=(𝑑!−1)𝑊𝑑,𝑛𝟏𝑑∣𝑛.
Suppose the asserted identities hold below 𝑛. For each proper divisor 𝑑, subtracting 𝑊𝑑,𝑛𝑈𝑑 cancels the sum at 𝑑 and changes none of the others. The sum
at 𝑛 remains 𝑛! −1, since 𝑊𝑛,𝑛 =1. This proves the identities by
induction.
For completeness, restrict to indices 1,…,𝑁. The columns 𝑒1,𝑇2,…,𝑇𝑁 form a triangular
integer matrix with diagonal 1, −1,…, −1. Replacing each 𝑇𝑛 by 𝑈𝑛 subtracts an integer combination of
earlier columns, so the second change of basis is triangular with
diagonal one. Both matrices have integer inverses: back-substitution
divides only by 1 or −1. Hence 𝑒1,𝑈2,…,𝑈𝑁 is a basis over ℤ.
In the resulting expansion, applying 𝑀 gives the coefficient of 𝑒1. If the coefficient of 𝑈𝑑 is 𝑐𝑑, applying 𝑉𝑑 gives 𝑉𝑑 =𝑀 +(𝑑! −1)𝑐𝑑, which proves (2). For
𝑑 >𝑁 that coefficient is zero
because 𝑉𝑑 =𝑀. ◻
The 𝑈𝑑-coordinates in (2) are
exactly the integer summands in (1). Their sum
is therefore 𝑘 in R =𝑀𝑆 +𝑘.
For a prime 𝑝 ≥3 the divisor
sum is empty, and 𝑈𝑝 =𝑝𝑒𝑝−1 −𝑒𝑝.
At a composite index the correction terms remove the extra weighted
sums. For example,
𝑈9=9𝑒8−𝑒9−5040𝑒2+1680𝑒3.(3)
Adding 𝑈𝑛
changes only 𝑉𝑛 among the weighted
sums, preserves 𝑀, and adds one to
R. Integer multiples
therefore preserve the fractional part.
Possible values of 𝑀
Fix 𝐷 ≥2 and solve 𝑉2 =⋯ =𝑉𝐷 =0, allowing index one
temporarily. The coordinates in (2) give
𝑉𝑑(𝜆)≡𝑀(𝜆)(mod𝑑!−1).(4)
It follows that 𝑀 must be divisible by
𝐿𝐷=lcm2≤𝑑≤𝐷(𝑑!−1).
This necessary condition is sufficient while index
one is allowed: 𝐿𝐷𝑒1 has every
weighted sum equal to 𝐿𝐷, and
subtracting 𝐿𝐷/(𝑑! −1) copies of
𝑈𝑑 sets the 𝑑th sum to zero without changing the
others. Thus put
𝐾𝐷=𝐿𝐷𝑒1−𝐷∑𝑑=2𝐿𝐷𝑑!−1𝑈𝑑.
Uniqueness of the coordinates now gives the complete solution
𝑉2=⋯=𝑉𝐷=0⟺𝜆=𝑡𝐾𝐷+∑𝑛>𝐷𝑧𝑛𝑈𝑛,𝑀=𝑡𝐿𝐷,(5)
with 𝑡 ∈ℤ and finitely many nonzero integers 𝑧𝑛. The coefficient of 𝑒1 in (2) is
𝑀; the original coordinate 𝜆1 also receives contributions
from the 𝑈𝑛. Write 𝑎𝐷 =(𝐾𝐷)1 and 𝑢𝑛 =(𝑈𝑛)1. The support restriction is
exactly
𝑡𝑎𝐷+∑𝑛>𝐷𝑧𝑛𝑢𝑛=0.(6)
The higher corrections leave cancellation
untouched and change the first coordinate by exactly 𝑔𝐷ℤ, where 𝑔𝐷 =gcd{𝑢𝑛 :𝑛 >𝐷}. They can
therefore remove 𝑡𝑎𝐷 precisely
when 𝑔𝐷 ∣𝑡𝑎𝐷. For example, at
depth four we shall find 𝑎4 = −55
and 𝑔4 =60. Removing index one then
requires 60 ∣55𝑡, or 12 ∣𝑡; hence 𝑀 =115𝑡 must be a multiple of 1380.
In general, the possible values of 𝑀 form the ideal
𝜇𝐷ℤ,𝜇𝐷=𝐿𝐷𝑔𝐷gcd(𝑔𝐷,𝑎𝐷).(7)
Indeed, division by gcd(𝑔𝐷,𝑎𝐷) leaves two coprime
integers, so 𝑔𝐷 ∣𝑡𝑎𝐷 holds
exactly when 𝑔𝐷/gcd(𝑔𝐷,𝑎𝐷) ∣𝑡. Bézout’s identity supplies the required finite combination;
the next theorem bounds the indices needed. Every vector attaining 𝜇𝐷 has coefficient gcd one, since
division by a common factor would give a smaller positive
multiplier.
A finite gcd calculation
The recursion for the first coordinate is
𝑢2=2,𝑢𝑛=−∑𝑑∣𝑛2≤𝑑<𝑛𝑊𝑑,𝑛𝑢𝑑(𝑛>2).
Induction gives 𝑢𝑛 =0 at odd indices and 𝑢2𝑝 = −(2𝑝)!/2𝑝−1 ≠0 at twice a
prime, including 𝑝 =2. Choose a
prime ℓ with 𝐷/2 <ℓ ≤𝐷, using Bertrand’s
postulate for 𝐷 ≥3 and ℓ =2 for 𝐷 =2. A finite gcd containing 𝑢2ℓ is positive and divides (2ℓ)!.
The remaining issue is whether later coefficients reduce this gcd. In
each recurrence term 𝑊𝑑,𝑛𝑢𝑑,
induction controls 𝑢𝑑 for 𝑑 >𝐷. For 𝑑 ≤𝐷, unordered block counting gives
(𝑛/𝑑)! ∣𝑊𝑑,𝑛. The cutoff
ensures 𝑛/𝑑 ≥2ℓ, so the weight
supplies the factor (2ℓ)!
instead.
Theorem 3.1 (a finite formula for the gcd).
Choose a prime ℓ with 𝐷/2 <ℓ ≤𝐷 and put 𝐻 =𝐷(2ℓ −1). Then
𝑔𝐷=gcd(𝑢𝐷+1,…,𝑢𝐻),𝐻<2𝐷2.
Proof. Let 𝑔 =gcd(𝑢𝐷+1,…,𝑢𝐻). The index
2ℓ belongs to this interval, and
𝑢2ℓ = −(2ℓ)!/2ℓ−1 ≠0.
Consequently 𝑔 >0 and 𝑔 ∣(2ℓ)!.
We prove 𝑔 ∣𝑢𝑛 for every
𝑛 >𝐷 by strong induction. This is
true up to 𝐻 by definition. For
𝑛 >𝐻, consider a proper divisor
𝑑 ≥2 in the recurrence for 𝑢𝑛. If 𝑑 >𝐷, the induction hypothesis gives
𝑔 ∣𝑢𝑑. If 𝑑 ≤𝐷, put 𝑘 =𝑛/𝑑. The choice 𝐻 =𝐷(2ℓ −1) ensures 𝑘 ≥2ℓ. Moreover,
𝑊𝑑,𝑛𝑘!=𝑛!(𝑑!)𝑘𝑘!∈ℤ:
this quotient counts partitions of 𝑛 labelled objects into 𝑘 unordered blocks of size 𝑑. Hence (2ℓ)! ∣𝑘! ∣𝑊𝑑,𝑛, so 𝑔 ∣𝑊𝑑,𝑛 in the small-divisor case
as well. Every product 𝑊𝑑,𝑛𝑢𝑑
in the recurrence is therefore divisible by 𝑔.
Thus 𝑔 divides the whole tail.
Conversely, every common divisor of the tail divides its finite
subfamily 𝑢𝐷+1,…,𝑢𝐻. The
two gcds are equal, and 𝐻 ≤𝐷(2𝐷 −1) <2𝐷2. ◻
The coefficient at 2ℓ ≤2𝐷
supplies positivity; the larger cutoff controls the small-divisor
weights. The argument does not justify truncation at 2𝐷.
To attain the least positive multiplier, compute integers 𝑏𝑛 with ∑𝐻𝑛=𝐷+1𝑏𝑛𝑢𝑛 =𝑔𝐷 by the extended
Euclidean algorithm, and put ℎ =gcd(𝑔𝐷,𝑎𝐷). Taking
𝑡=𝑔𝐷/ℎ,𝑧𝑛=−(𝑎𝐷/ℎ)𝑏𝑛(𝐷<𝑛≤𝐻)
makes 𝑡𝑎𝐷 +∑𝑧𝑛𝑢𝑛 =0. Substitution in (5)
therefore gives 𝑀 =𝜇𝐷 with
support in 2,…,𝐻.
The associated remainder
Fix 𝑀 ∈𝜇𝐷ℤ and put
𝑡 =𝑀/𝐿𝐷. For choices of 𝑧𝑛 satisfying (6), the
identities R(𝑒1) =𝑆 and
R(𝑈𝑛) =1 give
R(𝑡𝐾𝐷+∑𝑛>𝐷𝑧𝑛𝑈𝑛)=𝑡𝐿𝐷(𝑆−𝐻𝐷)+∑𝑛>𝐷𝑧𝑛,𝐻𝐷=𝐷∑𝑑=21𝑑!−1.(8)
Since 𝑀𝐻𝐷 ∈ℤ, the offset from 𝑀𝑆 is
𝑘 = −𝑀𝐻𝐷 +∑𝑛>𝐷𝑧𝑛; ∑𝑧𝑛 alone is the offset from 𝑀(𝑆 −𝐻𝐷). All choices at fixed 𝑀 therefore have the same fractional
part.
The least multiplier at depth four
At depth four the recursion gives
𝐿4=115,𝐾4=−55𝑒1+16𝑒2+3𝑒3+5𝑒4,𝑎4=−55.
Take ℓ =3 in Theorem 3.1,
so that 𝐻 =20. Since 𝑢6 = −180, 𝑢8 = −4200 and 23𝑢6 −𝑢8 =60, we have 𝑔4 ∣60. For the converse,
lcm(1,…,𝑛)∣𝑢𝑛(9)
gives 60 ∣𝑢𝑛 for every 𝑛 >4,
hence 60 ∣𝑔4. To prove (9), put
Λ𝑛 =lcm(1,…,𝑛).
For a prime 𝑟, write 𝑣𝑟 for its valuation. If 𝑑 ∣𝑛, Legendre’s formula gives
𝑣𝑟(𝑊𝑑,𝑛)=∑𝑗≥1(⌊𝑛𝑟𝑗⌋−𝑛𝑑⌊𝑑𝑟𝑗⌋)≥𝑣𝑟(Λ𝑛)−𝑣𝑟(Λ𝑑).
Each summand is nonnegative,
and each power 𝑑 <𝑟𝑗 ≤𝑛
contributes at least one. Hence Λ𝑛 ∣𝑊𝑑,𝑛Λ𝑑, and
induction in the recurrence, starting from 𝑢2 =2, proves the assertion. Thus 𝑔4 =60, and (7) gives
𝜇4=11560gcd(60,55)=1380.
To attain 𝑀 =1380, take 𝑡 =12 in (6). The choice
𝑧6 = −27, 𝑧8 =1 works because
12(−55)−27(−180)−4200=0.
The vector
12𝐾4 −27𝑈6 +𝑈8 therefore has no
index-one coefficient and cancels 𝑉2,𝑉3,𝑉4. To remove its coefficients
at indices four and seven, add −26𝑈5 +8𝑈7. These corrections still
cancel through four and do not restore index one, since 𝑢5 =𝑢7 =0. Thus
𝜆=12𝐾4−26𝑈5−27𝑈6+8𝑈7+𝑈8,
which expands to
𝜆=1482𝑒2−784𝑒3−136𝑒5+83𝑒6−𝑒8,(𝑀,𝑉2,𝑉3,𝑉4)=(1380,0,0,0).(10)
Table 1 verifies all four
sums directly.
The depth-four checks: multiply each of the last four columns
by 𝜆𝑛 and sum to obtain the
last row.
| 𝑛 |
𝜆𝑛 |
𝑛! |
𝑊2,𝑛 |
𝑊3,𝑛 |
𝑊4,𝑛 |
| 2 |
1482 |
2 |
1 |
2 |
2 |
| 3 |
-784 |
6 |
3 |
1 |
6 |
| 5 |
-136 |
120 |
30 |
20 |
5 |
| 6 |
83 |
720 |
90 |
20 |
30 |
| 8 |
-1 |
40320 |
2520 |
1120 |
70 |
| Weighted totals |
1380 |
0 |
0 |
0 |
The basis coordinates give the offset −26 −27 +8 +1 = −44 in (8). Since
𝐻4 =143/115, the original linear
form is
R(𝜆)=1380(𝑆−𝐻4)−44=1380𝑆−1760.
For the remainder
estimate, the inequality (𝑛 +1)! −1 >(𝑛 +1)(𝑛! −1) gives the
geometric bound
0<𝑆−𝐻𝑁<𝑁+2(𝑁+1)((𝑁+1)!−1)<2(𝑁+1)!−1(𝑁≥2).(11)
The majorant starts with 1/((𝑁 +1)! −1) and has ratio 1/(𝑁 +2).
Exact rational arithmetic gives
−31+45<13808∑𝑑=51𝑑!−1−44<−31+56,2⋅13809!−1<1100.
The omitted
contribution to R is
positive and less than 1/100, so
−31 <R(𝜆) < −30.
If 𝑆 =𝑎/𝑞 with 𝑞 ∣1380, (1) would make
this remainder an integer. Every divisor of 1380 is therefore excluded.
The depth-six
calculation gives 𝜇6 =24𝐿6:
the extra factor need not equal 12.
The companion also gives a dual
congruence, support bounds and coefficient-norm refinements. Least
multiplier, least support size and least coefficient norm are distinct
problems; (7)
settles only the first.
Nonintegrality of the remainder
Fix a vector with 𝑉2 =⋯ =𝑉𝐷 =0 and 𝑀 >0, and choose any cutoff 𝑁 ≥𝐷 at or beyond its support. Separate
the signed finite sum from the positive tail:
𝐴𝑁=𝑁∑𝑑=𝐷+1𝑉𝑑(𝜆)𝑑!−1,R(𝜆)=𝐴𝑁+𝑀(𝑆−𝐻𝑁).
The remainder exceeds 𝐴𝑁. By (11), it lies
below the least integer strictly above 𝐴𝑁 whenever
2𝑀(𝑁+1)!−1<⌊𝐴𝑁⌋+1−𝐴𝑁.(12)
Then ⌊𝐴𝑁⌋ <R(𝜆) <⌊𝐴𝑁⌋ +1,
excluding every denominator dividing 𝑀. The strict successor matters: at an
integral 𝐴𝑁 the gap is one,
whereas just below an integer it can be arbitrarily small. Failure of
the comparison is inconclusive.
At depth two, fix 𝜆 = −6𝑒2 +𝑒4, with 𝑀 =12, and vary only the cutoff. At 𝑁 =4, its finite part 𝐴4 = −239/115 lies just below −2, and its gap 9/115 is smaller than the bound 24/119. At 𝑁 =5, the finite part 𝐴5 = −27061/13685 has crossed −2: the next integer is now −1, and the gap 13376/13685 exceeds 24/719. For any fixed vector with 𝑀 >0, 𝐴𝑁 increases to R. If R is nonintegral, the gap tends
to a positive number and the test eventually succeeds. If R is integral, the eventual gap
equals the omitted tail, so the strict upper-bound test fails. Changing
the cutoff detects nonintegrality; it neither creates it nor changes
𝑀. Coefficient changes at fixed
𝑀 preserve nonintegrality by (8).
To exclude every denominator, 𝑀
must vary. The companion’s
progression construction adds size control: polynomial root
conditions give primitive vectors (coefficient gcd one) with 0 <𝑀 <𝑁!, where 𝑁 is their largest supported index. Every
fixed 𝑞 divides their multipliers
for all sufficiently large 𝑁. The
tail at that endpoint is below 2/𝑁,
but the next-integer gap may shrink too. For each 𝑞 >0, one still needs a vector with
𝑞 ∣𝑀 and a cutoff satisfying
(12); the
cutoff may exceed its support endpoint. No minimal-support assertion is
made.
The companion treats the equivalent carry
condition and the complementary-denominator
tail inequality. It also records two larger finite computations: a
non-unit carry at 𝑚 =300000 excludes
every 𝑞 ∣299999!, and a
continued-fraction enclosure excludes every 𝑞 <239990. Both computations were
performed outside Lean; their algorithms,
exact outputs and scope are given there. The factorial-divisibility
exclusion does not exclude all denominators with only small prime
factors, whose exponents may be larger. Neither finite restriction
establishes the required condition at arbitrarily large indices.
Verification
The integral-basis theorem and finite gcd formula have recorded exact
Lean bindings. The companion
source concordance identifies their declarations and the theorem on
the possible values of 𝑀. Its large
finite computations are separate from these kernel proofs.
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 𝑀 developed earlier here. OpenAI Codex
checked and integrated them; a separate AI pass reviewed the proof.