A366152 a(n) = denominator(denominator(Bernoulli'(n, x)) / denominator(Bernoulli(n, 1))).
1, 2, 6, 1, 30, 1, 42, 1, 10, 1, 66, 1, 2730, 1, 2, 1, 170, 1, 798, 1, 110, 1, 46, 1, 546, 1, 2, 1, 870, 1, 14322, 1, 170, 1, 2, 1, 1919190, 1, 2, 1, 4510, 1, 1806, 1, 46, 1, 94, 1, 1326, 1, 22, 1, 530, 1, 798, 1, 290, 1, 118, 1, 56786730, 1, 2, 1, 34, 1, 64722
Offset: 0
Programs
-
Maple
seq(denom(denom(diff(bernoulli(n, x), x))/denom(bernoulli(n, 1))), n = 0..66);
-
PARI
a(n) = denominator(lcm(apply(denominator, Vec(deriv(bernpol(n)))))/denominator(subst(bernpol(n, x), x, 1))); \\ Michel Marcus, Oct 14 2023