Introduction
Let
G(z)=exp(k≥1∑k(k+1)zk)=e(1−z)(1−z)/z=n≥0∑Enzn,
so that E0=1, E1=1/2, E2=7/24, E3=3/16, E4=743/5760. The second expression follows from ∑k≥1zk/(k(k+1))=1+(1−z)log(1−z)/z, obtained from 1/(k(k+1))=1/k−1/(k+1).
Putting z=1/(x+1) gives
e(1+x1)−x=n≥0∑(x+1)nEn.
Yang [4 Ya2] wrote (1+1/x)x=e(1−∑n≥1bn(x+1)−n) to refine Carleman’s inequality, so ∑nEnzn=1/(1−∑nbnzn); the En are the coefficients Sn of [NG], where the p-adic valuations of the bn are studied. Equivalently, 1/G(z)=(1−z)e−1(1−z)−1/z, where e−1(1+x)1/x has the coefficients tabulated in [OEIS]. We have found no previous result on the arithmetic of the En; see [1] for analytic properties of the bn.
Expanding the exponential, En is a sum over the partitions of n. For n=p the only terms whose denominator is divisible by p come from the partitions (p−1)+1, p and 1p, namely (p−1)!p2, p(p+1)1 and 2pp!1. Multiplied by p, they are congruent modulo p to −1/2, 1 and −1/2 (by Fermat and Wilson), so the poles cancel and Ep is p-integral. The point of the theorem is the next digit.
Theorem 1. For every odd prime p,
Ep=−2p(modp2). In particular vp(Ep)=1.
For p=2 the statement fails, since E2=7/24. For example E3=3/16 and E5=215/2304, with E3+3/2=27/16=33/24 and E5+5/2=5975/2304=52⋅239/(28⋅32).
The proof (Section 3) uses two expressions of Ep through the coefficients of Φ(u)=e(1+u)−1/u, the congruence (jp−1)≡(−1)j(1−pHj)(modp2), and the recurrence of the En obtained by differentiating G. Throughout, a≡b(modpj) means that a−b∈pjZ(p), for rationals a,b, and Hj=1+1/2+⋯+1/j, H0=0.
Identities
Put
Φ(u)=e(1+u)−1/u=exp(n≥1∑hnun),hn=n+1(−1)n+1,Φ(u)=k≥0∑φkuk.
the second expression coming from 1−log(1+u)/u=∑n≥1(−1)n+1un/(n+1). Thus φ0=1, φ1=1/2.
Lemma 2. As formal power series over Q,
(a)G(z)=1−zΦ(−z),(b)G(z)=Φ(1−zz),(c)(n+1)En+1=j=0∑nn+2−jEj(n≥0). Consequently, for n,k≥0 and every integer p≥1,
En=k=0∑n(−1)kφk,(−1)kφk=Ek−Ek−1 (E−1=0),Ep=k=1∑p(k−1p−1)φk. Proof. (a) Φ(−z)=e(1−z)1/z and G(z)=e(1−z)1/z−1. (b) With w=z/(1−z) we have 1+w=(1−z)−1 and −1/w=−(1−z)/z, so Φ(w)=e(1−z)(1−z)/z=G(z). (c) Differentiate logG=∑k≥1zk/(k(k+1)): G′=G∑m≥0zm/(m+2), and compare the coefficients of zn.
The first two identities in (1) follow from (a), since 1/(1−z)=∑j≥0zj. For the third,
[zp](1−zz)k=[zp−k](1−z)−k=(k−1p−1) for 1≤k≤p, and 0 for k=0 or k>p; apply (b). □
Fix an odd prime p. Let ψk be the coefficients of Ψ(u)=exp(∑n=1p−2hnun).
Lemma 3. (i) hn is p-integral for n≤p−2, and ψk is p-integral for k≤p−1. (ii) φk=ψk for k≤p−2; in particular φk and Ek are p-integral for k≤p−2. (iii) φp−1=ψp−1−1/p, and φp+1/p is p-integral. (iv) Ep=Ep−2+φp−1−φp; hence Ep is p-integral.
Proof. (i) The denominator of hn is n+1≤p−1. The coefficient ψk is a sum of terms ∏jhjmj/mj! over the partitions ∑jjmj=k with parts j≤p−2; for k≤p−1 every mj is less than $p.
(ii), (iii) Write Φ=Ψ⋅R with R=exp(∑n≥p−1hnun). Since 2(p−1)≥p+1, R≡1+hp−1up−1+hpup(modup+1), where hp−1=(−1)p/p=−1/p and hp=1/(p+1) is p-integral. Hence φk=ψk for k≤p−2, φp−1=ψp−1−1/p and φp=ψp+ψ1hp−1+hp=ψp−1/(2p)+hp.
In the expansion of ψp, all partitions of p into parts ≤p−2 have multiplicities less than p except 1p, which contributes h1p/p!=1/(2pp!). By Fermat and Wilson, p/(2pp!)=1/(2p(p−1)!)≡1/(2⋅(−1))=−1/2(modp), so ψp+1/(2p) is p-integral and so is φp+1/p. By (1), (i) and the first part of (ii), Ek=∑i≤k(−1)iφi is p-integral for k≤p−2.
(iv) By (1), Ep=Ep−2+(−1)p−1φp−1+(−1)pφp, and p is odd. The right-hand side is p-integral by (ii) and (iii). □
Proof of Theorem 1
Step A: reduction to a congruence modulo p. For 0≤j≤p−1, (jp−1)=∏i=1j(p−i)/i=(−1)j∏i=1j(1−p/i)≡(−1)j(1−pHj)(modp2). In the last expression of (1) the terms k≤p−2 have p-integral φk, so
Ep≡k=1∑p−2(−1)k−1φk−pH∗+(p−1)φp−1+φp(modp2),H∗=k=1∑p−2(−1)k−1Hk−1φk.
Here ∑k=1p−2(−1)k−1φk=1−Ep−2 and (p−1)φp−1+φp=pφp−1−(φp−1−φp). Adding Lemma 3(iv),
2Ep≡1+pφp−1−pH∗=p(φp−1−H∗)(modp2),
@ends_mid=1 by Lemma 3(iii). Since ψp−1 and H∗ are p-integral, this already gives Ep≡0(modp), and the theorem is equivalent to
ψp−1−H∗≡−1(modp).(1)
Step B: the residue ψp−1. From Ψ′=Ψ∑j=1p−2jhjuj−1, the coefficient of up−2 gives the exact identity (p−1)ψp−1=∑j=1p−2jhjψp−1−j, and ψp−1−j=φp−1−j by Lemma 3(ii). Put k=p−1−j, so 1≤k≤p−2. Modulo p, p−1≡−1, j≡−(k+1) and hj=(−1)p−k/(p−k)≡(−1)k/k, so jhj≡−(−1)k(1+1/k) and
ψp−1≡k=1∑p−2(−1)k(1+k1)φk(modp).
Since H∗=−∑k=1p−2(−1)kHk−1φk and Hk−1+1/k=Hk,
ψp−1−H∗≡k=1∑p−2(−1)k(1+Hk)φk=Ep−2−1+k=1∑p−2(−1)kHkφk(modp).
Hence (2) is equivalent to
k=1∑p−2(−1)kHkφk≡−Ep−2(modp).(2)
Step C: the recurrence. By (1) and summation by parts,
k=1∑p−2(−1)kHkφk=k=1∑p−2Hk(Ek−Ek−1)=Hp−2Ep−2−H1E0−k=1∑p−3(Hk+1−Hk)Ek=Hp−2Ep−2−1−k=1∑p−3k+1Ek.
Since Hp−1≡0(modp) (pair i with p−i), Hp−2=Hp−1−1/(p−1)≡1. So (3) is equivalent to
2Ep−2≡1+k=1∑p−3k+1Ek=j=0∑p−3j+1Ej(modp).
This is Lemma 2(c) with n=p−3, namely (p−2)Ep−2=∑j=0p−3Ej/(p−1−j), read modulo p: all Ej with j≤p−2 are p-integral, p−2≡−2 and 1/(p−1−j)≡−1/(j+1) for 0≤j≤p−3. This proves (3), hence (2) and the theorem. □
Remark 4. After Step A only E0,…,Ep−2 and φ1,…,φp−2 occur, and the congruence finally needed is the recurrence of the En at n=p−3. For p=3 all the sums in Steps B and C are empty or reduce to a single term, and the theorem reads E3=3/16≡−3/2(mod9).
Numerical verification and a remark
The ancillary file verify_Ep_congruence.py (Python 3, no dependency, exact rational arithmetic, about thirty seconds) computes En by the recurrence of Lemma 2(c) and φk from the exponential, and checks, for all odd primes p<200 and for p=251, 307, 401: Lemma 2(a) and (b) coefficientwise up to zp+1, and the third identity of (1); Lemma 3; each displayed congruence and equivalence of Steps A, B and C, and the summation by parts; and Theorem 1.
The next p-adic digit (Ep+p/2)/p2(modp) is 0,1,3,9,0,12,14,3 for p=3,5,…,23 and shows no pattern. In the tested range, Ep≡−p/2(modp3) holds only for p=3 and p=13, which is compatible with a heuristic probability 1/p for each p.
Acknowledgement
Acknowledgement. Technical assistance: Claude (Anthropic).
References
ORCID 0009-0009-1526-1178