OFFSET
0,2
MATHEMATICA
r[n_] := Sum[BernoulliB[1 + k] / ((1 + k) * k^(1 + k)), {k, 1, 2*n}];
Table[r[n], {n, 0, 10}] // Denominator
PROG
(PARI) a(n) = denominator(sum(k=1, 2*n, bernfrac(1+k)/((1+k)*k^(1+k)))); \\ Michel Marcus, Aug 12 2021
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Peter Luschny, Aug 11 2021
STATUS
approved