A033471 Denominator of Bernoulli(2n,1/3).
1, 18, 810, 10206, 65610, 1299078, 483611310, 9565938, 7317942570, 103053850074, 383546284110, 1443528742014, 257010878197710, 5083731656658, 6634269811938690, 982924264619854326, 315013432104812970
Offset: 0
Crossrefs
Cf. A033470 (numerators).
Programs
-
Maple
with(numtheory): seq(denom(bernoulli(2*n, 1/3)), n=0..20);
-
Mathematica
Denominator/@BernoulliB[2Range[0,20],1/3] (* Harvey P. Dale, May 11 2011 *)
-
PARI
a(n) = denominator(subst(bernpol(2*n), x, 1/3)); \\ Michel Marcus, Sep 16 2020