A366573 a(n) = denominator(Bernoulli'(n, x)) / denominator(Bernoulli''(n, x)).

1, 1, 1, 2, 1, 6, 1, 6, 1, 10, 1, 6, 1, 210, 1, 2, 3, 10, 1, 42, 1, 110, 3, 2, 1, 546, 1, 2, 1, 30, 1, 462, 1, 170, 3, 2, 1, 51870, 1, 2, 3, 110, 1, 42, 1, 46, 15, 2, 1, 1326, 1, 22, 3, 10, 1, 798, 1, 290, 3, 2, 1, 930930, 1, 2, 3, 34, 5, 966, 1, 2, 3, 110, 1

Offset: 0

Author

Peter Luschny, Oct 13 2023

Keywords

• nonn

Crossrefs

Cf. A324370/A366168, A366572, A144845/A366168.

Programs

Maple

seq(denom(diff(bernoulli(n, x), x))/denom(diff(diff(bernoulli(n, x), x),x)), n = 0..100);

PARI

a(n) = lcm(apply(denominator, Vec(deriv(bernpol(n)))))/ lcm(apply(denominator, Vec(deriv(deriv(bernpol(n)))))); \\ Michel Marcus, Oct 14 2023

Formula

a(n) = A324370(n) / A366168(n).