OFFSET
1,1
COMMENTS
A value of 1 indicates a Wolstenholme prime.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
MATHEMATICA
Table[Mod[Binomial[2p-1, p-1], p^4], {p, Prime[Range[30]]}] (* Harvey P. Dale, Jun 26 2017 *)
PROG
(PARI) forprime(n=2, 10^2, m=(binomial(2*n-1, n-1)%n^4); print1(m, ", "));
(Python)
from __future__ import division
from sympy import isprime
A242473_list, b = [], 1
for n in range(1, 10**4):
if isprime(n):
A242473_list.append(b % n**4)
b = b*2*(2*n+1)//(n+1) # Chai Wah Wu, Jan 26 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, May 26 2014
STATUS
approved