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-3 of 3 results.

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

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

References

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

Crossrefs

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

Programs

  • Magma
    [(Binomial(2*p-1,p)-1) div p^3: p in PrimesInInterval(4,100)]; // Vincenzo Librandi, Nov 23 2015
  • Maple
    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
  • Mathematica
    Table[(Binomial[2 Prime[n] - 1, Prime[n] - 1] - 1)/Prime[n]^3, {n, 3, 20}] (* Vincenzo Librandi, Nov 23 2015 *)

Formula

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

Extensions

Edited by Max Alekseyev, May 14 2010
More terms from Vincenzo Librandi, Nov 23 2015

A092101 Harmonic primes.

Original entry on oeis.org

5, 13, 17, 23, 41, 67, 73, 79, 107, 113, 139, 149, 157, 179, 191, 193, 223, 239, 241, 251, 263, 277, 281, 293, 307, 311, 317, 331, 337, 349, 431, 443, 449, 461, 467, 479, 487, 491, 499, 503, 541, 547, 557, 563, 569, 593, 619, 653, 683, 691, 709, 757, 769, 787
Offset: 1

Views

Author

T. D. Noe, Feb 20 2004

Keywords

Comments

For p = prime(n), Boyd defines J_p to be the set of numbers k such that p divides A001008(k), the numerator of the harmonic number H(k). For harmonic primes, J_p contains only the three numbers p-1, (p-1)p and (p-1)(p+1). It has been conjectured that there are an infinite number of these primes and that their density in the primes is 1/e.
Prime p=A000040(n) is in this sequence iff neither H(k) == 0 (mod p), nor H(k) == -A177783(n) (mod p) have solutions for 1 <= k <= p-2. - Max Alekseyev, May 13 2010

Crossrefs

Cf. A092102 (non-harmonic primes), A092103 (size of J_p).

Programs

  • PARI
    is(p)=my(K=-Mod((binomial(2*p-1, p)-1)/2/p^3,p),H=Mod(0,p));for(k=1,p-2,H+=1/k;if(H==0||H==K,return(0)));1 \\ Charles R Greathouse IV, Mar 16 2014

Extensions

More terms from Max Alekseyev, May 13 2010

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-3 of 3 results.