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