OFFSET
0,15
MAPLE
seq(numer(denom(diff(diff(bernoulli(n, x), x), x))/denom(bernoulli(n, 1))), n = 0..75);
PROG
(PARI) a(n) = numerator(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