A366572 - OEIS

%I #15 Oct 14 2023 13:13:24

%S 1,2,6,2,30,6,42,6,10,10,66,6,2730,210,2,6,510,10,798,42,110,330,138,

%T 2,546,546,2,2,870,30,14322,462,170,510,6,2,1919190,51870,2,6,13530,

%U 110,1806,42,46,690,1410,2,1326,1326,22,66,1590,10,798,798,290,870

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

%F a(n) = A144845(n) / A366168(n).

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

%o (PARI) a(n) = lcm(apply(denominator, Vec(bernpol(n))))/lcm(apply(denominator, Vec(deriv(deriv(bernpol(n)))))); \\ _Michel Marcus_, Oct 14 2023

%Y Cf. A144845/A366168, A366571, A144845/A324370 (1st derivative).

%K nonn

%O 0,2

%A _Peter Luschny_, Oct 13 2023