OFFSET
0,2
MAPLE
seq(denom(denom(diff(diff(bernoulli(n, x), x), x))/denom(bernoulli(n, 1))), n = 0..65);
PROG
(PARI) a(n) = denominator(lcm(apply(denominator, Vec(deriv(deriv(bernpol(n))))))/denominator(subst(bernpol(n, x), x, 1))); \\ Michel Marcus, Oct 14 2023
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Peter Luschny, Oct 13 2023
STATUS
approved