A033470 Numerator of Bernoulli(2n,1/3).
1, -1, 13, -121, 1093, -49205, 61203943, -5580127, 25949996501, -2832495743227, 101471818419863, -4469253897850313, 11126033443528968583, -3623473373642099263, 90551915146257709805297, -295654219044074800178128205
Offset: 0
Crossrefs
Cf. A033471 (denominators).
Programs
-
Maple
with(numtheory): seq(numer(bernoulli(2*n, 1/3)), n=0..20);
-
Mathematica
Table[Numerator @ BernoulliB[2*n, 1/3], {n, 0, 15}] (* Amiram Eldar, Sep 16 2020 *)
-
PARI
a(n) = numerator(subst(bernpol(2*n), x, 1/3)); \\ Michel Marcus, Sep 16 2020