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