OFFSET
3,2
COMMENTS
p is a Wolstenholme prime (A088164) iff a(n) = 0. This holds for n = 1944 and n = 157504.
When performing a search for Wolstenholme primes, one can choose an integer constant c >= 0 and record all primes with a(n) <= c in order to get a larger data set.
The values here appear to have a nicer asymptotic growth behavior than those in A260209.
The formula only returns integers for primes greater than 3. - Robert G. Wilson v, Jul 29 2015
LINKS
Felix Fröhlich, Table of n, a(n) for n = 3..10000
FORMULA
MATHEMATICA
f[n_] := Block[{p = Prime@ n}, (Mod[ Binomial[2p - 1, p - 1], p^4] - 1)/p^3]; Array[f, 60, 3] (* Robert G. Wilson v, Jul 29 2015 *)
PROG
(PARI) a(n) = p=prime(n); lift(Mod(binomial(2*p-1, p-1)\p^3, p))
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Jul 19 2015
STATUS
approved