OFFSET
0,2
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 *)
PROG
(PARI) a(n) = denominator(subst(bernpol(2*n), x, 1/3)); \\ Michel Marcus, Sep 16 2020
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
STATUS
approved