login
A366426
a(n) = numerator(denominator(Bernoulli''(n, x)) / denominator(Bernoulli(n, 1))).
3
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,15
FORMULA
a(n) = numerator(A366168(n) / A027642(n)).
MAPLE
seq(numer(denom(diff(diff(bernoulli(n, x), x), x))/denom(bernoulli(n, 1))), n = 0..75);
PROG
(PARI) a(n) = numerator(lcm(apply(denominator, Vec(deriv(deriv(bernpol(n))))))/denominator(subst(bernpol(n, x), x, 1))); \\ Michel Marcus, Oct 14 2023
CROSSREFS
Cf. A366168/A027642, A366427 (denominator), A366570/A366152 (1st derivative).
Sequence in context: A359667 A010129 A073050 * A154740 A179261 A154567
KEYWORD
nonn,frac
AUTHOR
Peter Luschny, Oct 13 2023
STATUS
approved