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