A033470 - OEIS

%I #18 Sep 16 2020 07:37:40

%S 1,-1,13,-121,1093,-49205,61203943,-5580127,25949996501,

%T -2832495743227,101471818419863,-4469253897850313,

%U 11126033443528968583,-3623473373642099263,90551915146257709805297,-295654219044074800178128205

%N Numerator of Bernoulli(2n,1/3).

%H <a href="/index/Be#Bernoulli">Index entries for sequences related to Bernoulli numbers.</a>

%p with(numtheory): seq(numer(bernoulli(2*n, 1/3)), n=0..20);

%t Table[Numerator @ BernoulliB[2*n, 1/3], {n, 0, 15}] (* _Amiram Eldar_, Sep 16 2020 *)

%o (PARI) a(n) = numerator(subst(bernpol(2*n), x, 1/3)); \\ _Michel Marcus_, Sep 16 2020

%Y Cf. A033471 (denominators).

%K sign,frac

%O 0,3

%A _N. J. A. Sloane_.