%I #18 Oct 28 2022 18:23:16
%S 1,1,1,1,5,16,7,17,28,20,26,34,25,47,53,44,79,79,88,82,88,86,104,146,
%T 113,152,145,164,152,158,160,206,164,190,193,250,205,262,259,341,257,
%U 278,331,301,323,362,295,356,349,355,346,359,380,475,457,449,415,464
%N Sum of digits of numerator of Bernoulli number B(2n).
%F a(n) = A007953(A000367(n)). - _Michel Marcus_, Feb 16 2021
%e Numerator(B(2*9))=43867 and 4+3+8+6+7=28 hence a(9)=28.
%t Total[IntegerDigits[Numerator[#]]]&/@BernoulliB[2*Range[60]] (* _Harvey P. Dale_, Oct 28 2022 *)
%o (PARI) a(n) = sumdigits(numerator(bernfrac(2*n))); \\ _Michel Marcus_, Feb 16 2021
%o (Python)
%o from sympy import bernoulli
%o def a(n): return sum(map(int, str(abs(bernoulli(2*n).numerator()))))
%o print([a(n) for n in range(1, 59)]) # _Michael S. Branicky_, Jun 03 2021
%Y Cf. A000367, A007953, A043296.
%K base,easy,nonn
%O 1,5
%A _Benoit Cloitre_, Mar 24 2002