cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A088164 Wolstenholme primes: primes p such that binomial(2p-1,p-1) == 1 (mod p^4).

Original entry on oeis.org

16843, 2124679
Offset: 1

Views

Author

Christian Schroeder, Sep 21 2003

Keywords

Comments

McIntosh and Roettger showed that the next term, if it exists, must be larger than 10^9. - Felix Fröhlich, Aug 23 2014
When cb(m)=binomial(2m,m) denotes m-th central binomial coefficient then, obviously, cb(a(n))=2 mod a(n)^4. I have verified that among all naturals 1A246134). One might therefore wonder whether this is true in general. - Stanislav Sykora, Aug 26 2014
Romeo Mestrovic, Congruences for Wolstenholme Primes, Lemma 2.3, shows that the criterion for p to be a Wolstenholme prime is equivalent to p dividing A027641(p-3). In 1847 Cauchy proved that this was a necessary condition for the failure of the first case of Fermat's Last Theorem for the exponent p (see Ribenboim, 13 Lectures, p. 29). - John Blythe Dobson, May 01 2015
Primes p such that p^3 divides A001008(p-1) (Zhao, 2007, p. 18). Also: Primes p such that (p, p-3) is an irregular pair (cf. Buhler, Crandall, Ernvall, Metsänkylä, 1993, p. 152). Keith Conrad observes that for the two known (as of 2015) terms ord_p(H_p-1) = 3 is satisfied, where ord_p(H_p-1) gives the p-adic valuation of H_p-1 (cf. Conrad, p. 5). Romeo Mestrovic conjectures that p is a Wolstenholme prime if and only if S_(p-2)(p) == 0 (mod p^3), where S_k(i) denotes the sum of the k-th powers of the positive integers up to and including (i-1) (cf. Mestrovic, 2012, conjecture 2.10). - Felix Fröhlich, May 20 2015
Primes p that divide the Wolstenholme quotient W_p (A034602). Also, primes p such that p^2 divides the Babbage quotient b_p (A263882). - Jonathan Sondow, Nov 24 2015
The only known composite numbers n such that binomial(2n-1, n-1) is congruent to 1 mod n^2 are the numbers n = p^2 where p is a Wolstenholme prime: see A267824. - Jonathan Sondow, Jan 27 2016
The converse of Wolstenholme's theorem implies that if an integer n satisfies the congruence binomial(2*n-1, n-1) == 1 (mod n^4), then n is a term of this sequence, i.e., then n is necessarily prime, or, equivalently, A298946(i) > 1 for all i > 0. Whether this is true for all such n is an open problem. - Felix Fröhlich, Feb 21 2018
Primes p such that binomial(2*p-1, p-1) == 1-2*p*Sum_{k=1..p-1} 1/k - 2*p^2*Sum_{k=1..p-1} 1/k^2 (mod p^7) (cf. Mestrovic, 2011, Corollary 4). - Felix Fröhlich, Feb 21 2018
These are primes p such that p^2 divides A007406(p-1) (Mestrovic, 2015, p. 241, Lemma 2.3). - Amiram Eldar and Thomas Ordowski, Jul 29 2019
If a third Wolstenholme prime exists it is larger than 6*10^10 (cf. Hathi, Mossinghoff, Trudgian, 2021). - Felix Fröhlich, Apr 27 2021
Named after the English mathematician Joseph Wolstenholme (1829-1891). - Amiram Eldar, Jun 10 2021
Primes p such that tanh(Sum_{k=1..p-1} artanh(k/p)) == 0 (mod p^4). - Thomas Ordowski, Apr 17 2025

References

  • Richard K. Guy, Unsolved Problems in Number Theory, Sect. B31.
  • Paulo Ribenboim, 13 Lectures on Fermat's Last Theorem (Springer, 1979).
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 23.

Crossrefs

Programs

  • Magma
    [p: p in PrimesUpTo(2*10^4)| (Binomial(2*p-1,p-1) mod (p^4)eq 1)]; // Vincenzo Librandi, May 02 2015
  • Mathematica
    For[i = 2, i <= 20000, i++, {If[PrimeQ[i] && Mod[Binomial[2*i - 1, i - 1], i^4] == 1, Print[i]]}] (* Dylan Delgado, Mar 02 2021 *)
  • PARI
    forprime(n=2, 10^9, if(Mod(binomial(2*n-1, n-1), n^4)==1, print1(n, ", "))); \\ Felix Fröhlich, May 18 2014
    

Formula

A000984(a(n)) = 2 mod a(n)^4. - Stanislav Sykora, Aug 26 2014
A099908(a(n)) == 1 mod a(n)^4. - Jonathan Sondow, Nov 24 2015
A034602(PrimePi(a(n))) == 0 mod a(n) and A263882(PrimePi(a(n))) == 0 mod a(n)^2. - Jonathan Sondow, Dec 03 2015

A246130 Binomial(2n,n)-2 mod n.

Original entry on oeis.org

0, 0, 0, 0, 0, 4, 0, 4, 0, 4, 0, 2, 0, 4, 13, 4, 0, 4, 0, 18, 4, 4, 0, 10, 0, 4, 18, 26, 0, 2, 0, 4, 7, 4, 5, 14, 0, 4, 18, 18, 0, 40, 0, 2, 43, 4, 0, 10, 0, 4, 1, 42, 0, 4, 30, 30, 37, 4, 0, 34, 0, 4, 10, 4, 3, 64, 0, 34, 64, 38, 0, 34, 0, 4, 43, 30, 75, 64, 0, 18, 18, 4, 0, 26, 63, 4, 76, 86, 0, 38, 89, 22, 18, 4, 3, 58, 0
Offset: 1

Views

Author

Stanislav Sykora, Aug 16 2014

Keywords

Comments

By Wolstenholme's theorem, when n>3 is prime and cb(n) is the central binomial coefficient A000984(n), then cb(n)-2 is divisible by n^3. This implies that it is also divisible by n^e for e=1,2 and 3, but not necessarily for e=4. It follows also that cn(n)-2, with cn(n)=cb(n)/(n+1) being the n-th Catalan number A000108(n), is divisible by any prime n. In fact, for any n>0, cn(n)-2 = (n+1)cb(n)-2 implies (cn(n)-2) mod n = (cb(n)-2) mod n = a(n). The sequence a(n) is of interest as a prime-testing sequence similar to Fermat's, albeit not a practical one until/unless an efficient algorithm to compute moduli of binomial coefficients is found. For more info, see A246131 through A246134.

Examples

			a(7)=0 because cb(7)-2 = binomial(14,7) -2 = 3432-2 = 490*7. Check also that cn(7) = 3432/8 = 429 and 429-2 = 61*7 so that (cn(7)-2) mod 7 = 0.
		

Crossrefs

Cf. A000108, A000984, A128311, A246131 (pseudoprimes of a(n)), A246132 (e=2), A246133 (e=3), A246134 (e=4).

Programs

  • PARI
    a(n) = (binomial(2*n,n)-2)%n

Formula

For any prime p, a(p)=0.

A246133 a(n) = (binomial(2n, n) - 2) mod n^3.

Original entry on oeis.org

0, 4, 18, 4, 0, 58, 0, 68, 504, 754, 0, 1562, 0, 2062, 2518, 580, 0, 922, 0, 818, 6535, 7990, 0, 12058, 250, 4398, 2691, 10358, 0, 12422, 0, 16964, 10666, 29482, 3680, 42818, 0, 41158, 19791, 13618, 0, 54430, 0, 71942, 40993, 73006, 0, 12058, 3430, 122254, 98278, 127494, 0
Offset: 1

Views

Author

Stanislav Sykora, Aug 16 2014

Keywords

Comments

When e=3, the numbers binomial(2n, n) - 2 mod n^e are 0 whenever n is a prime greater than 3 (Wolstenholme's theorem; see A246130 for introductory comments). No composite number n for which a(n)=0 was found up to n=431500 (conjecture: there are none, and a(n)=0 for n>3 is a deterministic primality test).

Examples

			a(7)= (binomial(14,7)-2) mod 7^3 = (3432-2) mod 343 = 10*343 mod 343 = 0.
		

Crossrefs

Cf. A000984, A246130 (e=1), A246132 (e=2), A246134 (e=4).

Programs

  • Maple
    seq(binomial(2*n,n)-2 mod n^3, n=1..100); # Robert Israel, Aug 17 2014
  • Mathematica
    Table[Mod[Binomial[2 n, n] - 2, n^3], {n, 60}] (* Wesley Ivan Hurt, May 25 2024 *)
  • PARI
    a(n) = (binomial(2*n,n)-2)%n^3

Formula

For any prime p>3, a(p)=0.

A246134 Binomial(2n, n) - 2 mod n^4.

Original entry on oeis.org

0, 4, 18, 68, 250, 922, 1029, 580, 2691, 4754, 2662, 8474, 4394, 10294, 2518, 49732, 29478, 65074, 123462, 128818, 6535, 93174, 36501, 12058, 187750, 162582, 297936, 273782, 536558, 741422, 59582, 16964, 118477, 540434, 132305, 136130, 1114366, 1138598, 2214594, 2381618, 1860867, 2795686, 1828661, 1775622, 2683618, 1435710, 1557345, 3882778
Offset: 1

Views

Author

Stanislav Sykora, Aug 16 2014

Keywords

Comments

For e > 3, unlike the cases e=1,2,3, the numbers binomial(2n, n) - 2 mod n^e are not necessarily 0 for any n>1, be it prime or composite (see A246130 for introductory comments). Testing up to n=278000, the only number n>1 for which a(n)=0 is the first Wolstenholme prime 16843 (A088164), but no composite.

Examples

			a(7) = (binomial(14,7)-2) mod 7^4 = (3432-2) mod 2401 = 1029.
		

Crossrefs

Cf. A000984, A088164, A246130 (e=1), A246132 (e=2), A246133 (e=3).

Programs

  • PARI
    a(n) = (binomial(2*n,n)-2)%n^4

A267824 Composite numbers n such that binomial(2n-1, n-1) == 1 (mod n^2).

Original entry on oeis.org

283686649, 4514260853041
Offset: 1

Views

Author

Jonathan Sondow, Jan 25 2016

Keywords

Comments

Babbage proved the congruence holds if n > 2 is prime.
See A088164 and A263882 for references, links, and additional comments.
Conjecture: n is a term if and only if n = A088164(i)^2 for some i >= 1 (cf. McIntosh, 1995, p. 385). - Felix Fröhlich, Jan 27 2016
The "if" part of the conjecture is true: see the McIntosh reference. - Jonathan Sondow, Jan 28 2016
The above conjecture implies that this sequence and A228562 are disjoint. - Felix Fröhlich, Jan 27 2016
Composites c such that A281302(c) > 1. - Felix Fröhlich, Feb 21 2018

Examples

			a(1) = 16843^2 and a(2) = 2124679^2 are squares of Wolstenholme primes A088164.
		

Crossrefs

Showing 1-5 of 5 results.