A128673 -id:A128673 - cp's OEIS frontend

A034602 Wolstenholme quotient W_p = (binomial(2p-1,p) - 1)/p^3 for prime p=A000040(n).

1, 5, 265, 2367, 237493, 2576561, 338350897, 616410400171, 7811559753873, 17236200860123055, 3081677433937346539, 41741941495866750557, 7829195555633964779233, 21066131970056662377432067, 59296957594629000880904587621, 844326030443651782154010715715

Offset: 3

N. J. A. Sloane

nonn

Equivalently, (binomial(2p,p)-2)/(2*p^3) where p runs through the primes >=5.

The values of this sequence's terms are replicated by conjectured general formula, given in A223886 (and also added to the formula section here) for k=2, j=1 and n>=3. - Alexander R. Povolotsky, Apr 18 2013

			Binomial(10,5)-2 = 250; 5^3=125 hence a(5)=1.

R. K. Guy, Unsolved Problems in Number Theory, Sect. B31.

Robert Israel, Table of n, a(n) for n = 3..263
R. R. Aidagulov and M. A. Alekseyev. On p-adic approximation of sums of binomial coefficients. Journal of Mathematical Sciences 233:5 (2018), 626-634. doi:10.1007/s10958-018-3948-0 arXiv:1602.02632
R. J. McIntosh, On the converse of Wolstenholme's theorem, Acta Arithmetica 71:4 (1995), 381-389.
Romeo Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.
Jonathan Sondow, Extending Babbage's (non-)primality tests, in Combinatorial and Additive Number Theory II, Springer Proc. in Math. & Stat., Vol. 220, 269-277, CANT 2015 and 2016, New York, 2017; arXiv:1812.07650 [math.NT], 2018.

Cf. A177783 (alternative definition of Wolstenholme quotient), A072984, A092101, A092103, A092193, A128673, A217772, A223886, A263882.

Cf. A268512, A268589, A268590.

[(Binomial(2*p-1,p)-1) div p^3: p in PrimesInInterval(4,100)]; // Vincenzo Librandi, Nov 23 2015

f:= proc(n) local p;
p:= ithprime(n);
(binomial(2*p-1,p)-1)/p^3
end proc:
map(f, [$3..30]); # Robert Israel, Dec 19 2018

Table[(Binomial[2 Prime[n] - 1, Prime[n] - 1] - 1)/Prime[n]^3, {n, 3, 20}] (* Vincenzo Librandi, Nov 23 2015 *)

a(n) = (A088218(p)-1)/p^3 = (A001700(p-1)-1)/p^3 = (A000984(p)-2)/(2*p^3), where p=A000040(n).
a(n) = A087754(n)/2.
a(n) = (binomial(j*k*prime(n), j*prime(n)) - binomial(k*j, j)) / (k*prime(n)^3) for k=2, j=1, and n>=3. - Alexander R. Povolotsky, Apr 18 2013
a(n) = A263882(n)/prime(n) for n > 2. - Jonathan Sondow, Nov 23 2015
a(n) = numerator(tanh(Sum_{k=1..p-1} artanh(k/p)))/p^3, where p = prime(n) for n >= 3. - Thomas Ordowski, Apr 17 2025

Edited by Max Alekseyev, May 14 2010

More terms from Vincenzo Librandi, Nov 23 2015

A128672 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the m-th alternating generalized harmonic number H'(m,k), for k = 2.

Original entry on oeis.org

20, 42, 100, 110, 156, 272, 294, 342, 500, 506, 812, 930, 1210, 1332, 1640, 1806, 2028, 2058, 2162, 2500, 2756, 3422, 3660, 4422, 4624, 4970, 5256, 6162, 6498, 6806, 7832, 9312, 10100, 10506, 11026, 11342, 11638, 11772, 12500, 12656, 13310, 14406, 16002, 17030

Offset: 1

Alexander Adamchuk, Mar 20 2007

nonn

Generalized harmonic numbers are defined as H(m,k) = Sum_{j=1..m} 1/j^k. Alternating generalized harmonic numbers are defined as H'(m,k) = Sum_{j=1..m} (-1)^(j+1)/j^k.

Sequence contains all geometric progressions of the form (p-1)*p^k for k > 0 and some primes p > 3. Note the factorization of initial terms of {a(n)} = {4*5, 6*7, 4*5^2, 10*11, 12*13, 16*17, 6*7^2, 18*19, 4*5^3, 22*23, 28*29, 30*31, 10*11^2, 36*37, 40*41, 42*43, 12*13^2, 6*7^3, 46*47, 4*5^4, 52*53, 58*59, 60*61, 66*67, 16*17^2, 70*71, 72*73, 78*79, 18*19^2, 82*83, ...}. The smallest term that does not fit this pattern is 11026 = ((149-1)/2) * 149.

Eric Weisstein's World of Mathematics, Harmonic Number

Similar sequences for generalized harmonic numbers with different k: A125581 (k=1), A128673 (k=3), A128674 (k=4), A128675 (k=5); A128676 (k=6).
For the least numbers k > 0 such that k^n does not divide the denominator of H(k,n) nor the denominator of H'(k,n), see A128670. See also A128671(n) = A128670(prime(n)).
Cf. A001008, A002805, A058313, A058312, A007406, A007407, A119682, A125581, A126196, A126197, A128674, A128675, A128676, A128673, A128670, A129671.

k=2; f=0; g=0; Do[ f=f+1/n^k; g=g+(-1)^(n+1)*1/n^k; kf=Denominator[f]; kg=Denominator[g]; If[ !IntegerQ[kf/n^k] && !IntegerQ[kg/n^k], Print[n] ], {n,1,7000} ]

Edited and extended by Max Alekseyev, May 07 2010

A128676 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the m-th alternating generalized harmonic number H'(m,k), for k = 6.

Original entry on oeis.org

20, 100, 110, 156, 161, 272, 342, 345, 500, 506, 812, 930, 1210, 1332, 1640, 1806, 2028, 2162, 2500, 2756, 3051, 3422, 3660, 3703, 4422, 4624, 4970, 5256, 6162, 6498, 6806, 7832, 7935, 9312, 9605, 10100, 10506, 11342, 11638, 11772, 12500, 12656, 13310

Offset: 1

Alexander Adamchuk, Mar 20 2007

nonn

Generalized harmonic numbers are defined as H(m,k) = Sum_{j=1..m} 1/j^k. Alternating generalized harmonic numbers are defined as H'(m,k) = Sum_{j=1..m} (-1)^(j+1)/j^k.

Sequence contains all terms of geometric progressions of the form (p-1)*p^k, k > 0, for some primes p >= 5, such as 4*5^k, 7*23^k, 15*23^k, 27*113^k, etc. Note the factorization of initial terms of {a(n)} = {4*5, 4*5^2, 10*11, 12*13, 7*23, 16*17, 18*19, 15*23, 4*5^3, 22*23, 28*29, 30*31, 10*11^2, 36*37, 40*41, 42*43, 12*13^2, 46*47, 4*5^4, 52*53, 27*113, 58*59, 60*61, 7*23^2, ...}. See more details in Comments at A128672 and A125581.

Eric Weisstein's World of Mathematics, Harmonic Number

Cf. A001008, A002805, A058313, A058312.

Cf. A007406, A007407, A119682, A007410, A120296, A099828, A125581, A126196, A126197, A128672, A128673, A128676.

k=6; f=0; g=0; Do[ f=f+1/n^k; g=g+(-1)^(n+1)*1/n^k; kf=Denominator[f]; kg=Denominator[g]; If[ !IntegerQ[kf/n^k] && !IntegerQ[kg/n^k], Print[n] ], {n,1,3703} ]

Edited and extended by Max Alekseyev, May 08 2010

A128675 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the n-th alternating generalized harmonic number H'(m,k), for k = 5.

Original entry on oeis.org

444, 666, 888, 1332, 16428, 24642, 32856, 49284, 607836, 911754, 1215672, 1823508

Offset: 1

Alexander Adamchuk, Mar 20 2007

Generalized harmonic numbers are defined as H(m,k) = Sum_{j=1..m} 1/j^k. Alternating generalized harmonic numbers are defined as H'(m,k) = Sum_{j=1..m} (-1)^(j+1)/j^k.

Sequence contains all terms of geometric progressions 37^k*(37-1)/3, 37^k*(37-1)/2, 37^k*(37-1)*2/3, 37^k*(37-1) for k > 0. Note the factorization of initial terms of {a(n)} = {37*12, 37*18, 37*24, 37*36, ...}. See more details in Comments at A128672 and A125581.

Eric Weisstein's World of Mathematics, Harmonic Number

Cf. A001008, A002805, A058313, A058312, A007406, A007407, A119682, A007410, A120296, A099828, A125581, A126196, A126197, A128672, A128673, A128676.

k=5; f=0; g=0; Do[ f=f+1/n^k; g=g+(-1)^(n+1)*1/n^k; kf=Denominator[f]; kg=Denominator[g]; If[ !IntegerQ[kf/n^k] && !IntegerQ[kg/n^k], Print[n] ], {n,1,2000} ]

Eight more terms from Max Alekseyev, May 08 2010

A128670 Least number k > 0 such that k^n does not divide the denominator of generalized harmonic number H(k,n) nor the denominator of alternating generalized harmonic number H'(k,n).

Original entry on oeis.org

77, 20, 94556602, 42, 444, 20, 104, 42, 76, 20, 77, 110, 3504, 20, 903, 42, 1107, 20, 104, 42, 77, 20, 2948, 110, 136, 20, 76, 42, 903, 20, 77, 42, 268, 20, 7004, 110, 1752, 20, 19203, 42, 77, 20, 104, 42, 76, 20, 370, 110, 1107, 20, 77, 42, 12246, 20, 104, 42

Offset: 1

Alexander Adamchuk, Mar 24 2007

nonn

Generalized harmonic numbers are defined as H(m,k) = Sum_{j=1..m}1/j^k. Alternating generalized harmonic numbers are defined as H'(m,k) = Sum_{j=1..m} (-1)^(j+1)/j^k.
Some apparent periodicity in {a(n)} (not without exceptions): a(n) = 20 for n = 2 + 4m, a(n) = 42 for n = 4 + 12m and 8 + 12m, a(n) = 76 for n = 9 + 18m, a(n) = 77 for n = 1 + 10m, a(n) = 104 for n = 7 + 12m, a(n) = 110 for n = 12m, a(n) = 136 for n = 25 + 32m, etc.
See more details in Comments at A128672 and A125581.

Max Alekseyev, Table of n, a(n) for n=1..158.
Eric Weisstein's World of Mathematics, Harmonic Number

Cf. A001008, A002805, A058313, A058312, A007406, A007407, A119682, A007410, A120296, A125581, A126196, A126197, A128672, A128673, A128674, A128675, A128676, A128671, A128670.

More terms and b-file from Max Alekseyev, May 07 2010

A128671 Least number k > 0 such that k^p does not divide the denominator of generalized harmonic number H(k,p) nor the denominator of alternating generalized harmonic number H'(k,p), where p = prime(n).

Original entry on oeis.org

20, 94556602, 444, 104, 77, 3504, 1107, 104, 2948, 903, 77, 1752, 77, 104, 370

Offset: 1

Alexander Adamchuk, Mar 24 2007, Mar 26 2007

Generalized harmonic numbers are defined as H(m,k) = Sum_{i=1..m} 1/i^k. Alternating generalized harmonic numbers are defined as H'(m,k) = Sum_{i=1..m} (-1)^(i+1)*1/i^k.
a(18)..a(24) = {77,104,77,136,104,370,136}. a(26)..a(27) = {77,104}.
a(n) is currently unknown for n = {16,17,25,...}. See more details in Comments at A128672 and A125581.

			a(2) = A128673(1) = 94556602.

Eric Weisstein's World of Mathematics, Harmonic Number

Cf. A128670, A001008, A002805, A058313, A058312, A007406, A007407, A119682, A007410, A120296, A125581, A126196, A126197, A128672, A128673, A128674, A128675, A128676.

a(n) = A128670(prime(n)).

a(9) = 2948 and a(12) = 1752 from Max Alekseyev

Edited by Max Alekseyev, Feb 20 2019

A177783 Wolstenholme quotient of prime p=A000040(n), i.e., such integer m

Original entry on oeis.org

3, 6, 6, 7, 10, 14, 18, 20, 16, 24, 17, 38, 39, 19, 29, 28, 12, 53, 31, 19, 53, 58, 48, 42, 1, 33, 53, 37, 5, 81, 4, 17, 29, 13, 13, 72, 75, 70, 173, 159, 111, 150, 39, 178, 106, 163, 196, 163, 172, 30, 98, 24, 177, 261, 212, 223, 122, 147, 276, 17, 92, 111, 27, 209, 241

Offset: 3

Max Alekseyev, May 13 2010

nonn

a(n) = 0 iff A000040(n) is a Wolstenholme prime (given by A088164).

For n>2 and p=A000040(n), H(p^2-p) == H(p^2-1) == a(n)*p (mod p^2).

David W. Boyd, A p-adic study of the partial sums of the harmonic series, Experimental Math., Vol. 3 (1994), No. 4, 287-302.
R. Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.
R. Mestrovic, On a Congruence Modulo n^3 Involving Two Consecutive Sums of Powers, Journal of Integer Sequences, Vol. 17 (2014), 14.8.4.
Jianqiang Zhao, Bernoulli numbers, Wolstenholme's theorem, and p^5 variations of Lucas' theorem, Journal of Number Theory, Volume 123, Issue 1, March 2007, Pages 18-26.

Cf. A034602, A072984, A087754, A092101, A092103, A092193, A128673.

{ a(n) = my(p); p=prime(n); ((binomial(2*p-1,p)-1)/2/p^3)%p }

a(n) = H(p-1)/p^2 mod p = A001008(p-1)/A002805(p-1)/p^2 mod p = A034602(n)/2 mod p = (binomial(2*p-1,p)-1)/(2*p^3) mod p, where p = A000040(n).
a(n) = (-1/3)*B(p-3) mod p, with p=prime(n) and B(n) is the n-th Bernoulli number. - Michel Marcus, Feb 05 2016
a(n) = A087754(n)/4 mod A000040(n).

Edited by Max Alekseyev, May 16 2010

cp's OEIS Frontend

A034602 Wolstenholme quotient W_p = (binomial(2p-1,p) - 1)/p^3 for prime p=A000040(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

Extensions

A128672 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the m-th alternating generalized harmonic number H'(m,k), for k = 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Extensions

A128676 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the m-th alternating generalized harmonic number H'(m,k), for k = 6.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Extensions

A128675 Numbers m such that m^k does not divide the denominator of the m-th generalized harmonic number H(m,k) nor the denominator of the n-th alternating generalized harmonic number H'(m,k), for k = 5.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Extensions

A128670 Least number k > 0 such that k^n does not divide the denominator of generalized harmonic number H(k,n) nor the denominator of alternating generalized harmonic number H'(k,n).

Views

Author

Keywords

Comments

Links

Crossrefs

Extensions

A128671 Least number k > 0 such that k^p does not divide the denominator of generalized harmonic number H(k,p) nor the denominator of alternating generalized harmonic number H'(k,p), where p = prime(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A177783 Wolstenholme quotient of prime p=A000040(n), i.e., such integer m

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Extensions