OFFSET
3,1
COMMENTS
a(n)==1 mod (prime(n))
MATHEMATICA
Table[p=Prime[n]; 3/(2p) (2p+1)(3^p+1)BernoulliB[2p], {n, 3, 10}] (* Harvey P. Dale, Aug 21 2013 *)
PROG
(PARI) a(n)=3/2/prime(n)*(2*prime(n)+1)*(3^prime(n)+1)*bernfrac(2*prime(n))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 11 2004
STATUS
approved