A366426 a(n) = numerator(denominator(Bernoulli''(n, x)) / denominator(Bernoulli(n, 1))).

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 3, 5, 5, 7, 1, 5, 15, 1, 1, 7, 7, 1, 1, 1, 1, 77, 7, 35, 3, 1, 1, 455, 35, 7, 21, 55, 55, 7, 7, 7, 105, 1, 5, 221, 13, 11, 33, 55, 1, 19, 1, 5, 15, 1, 1, 5005, 715, 143, 33, 17, 85, 161, 35, 1, 3, 11, 55, 95095

Offset: 0

Author

Peter Luschny, Oct 13 2023

Keywords

• nonn
• frac

Crossrefs

Cf. A366168/A027642, A366427 (denominator), A366570/A366152 (1st derivative).

Programs

Maple

seq(numer(denom(diff(diff(bernoulli(n, x), x),x))/denom(bernoulli(n, 1))), n = 0..75);

PARI

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

Formula

a(n) = numerator(A366168(n) / A027642(n)).