OFFSET
1,1
COMMENTS
If p is a Sophie Germain prime (A005384) then denominator(B(2p)) = 6*(2p+1).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
dbn[n_]:=Module[{d=Denominator[BernoulliB[2n]]}, (d-6)/(12n)]; dbn/@ Prime[ Range[100]] (* Harvey P. Dale, May 19 2012 *)
PROG
(PARI) a(n) = my(p=prime(n)); (denominator(bernfrac(2*p)) - 6)/(12*p); \\ Michel Marcus, Sep 02 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 09 2003
STATUS
approved
