OFFSET
0,3
COMMENTS
From Wolfdieter Lang, Jul 05 2017: (Start)
a(n) gives also the numerators of the generalized Bernoulli numbers B[5,2](n) = 5^n*Bernoulli(n, 2/5) with the Bernoulli polynomials B(n, x) = Bernoulli(n, x) from A196838/A196839 or A053382/A053383. For the denominators see A288872(n) = A157867(n)/5^n. See a comment under A157866 for B[d,a](n).
(-1)^n*a(n) gives the numerators of the generalized Bernoulli numbers B[5,3](n); the denominators are A288872(n).
(End)
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..250
MATHEMATICA
Table[Numerator[BernoulliB[n, 2/5]], {n, 0, 50}] (* Vincenzo Librandi, Mar 16 2014 *)
PROG
(PARI) a(n) = numerator(subst(bernpol(n, x), x, 2/5)); \\ Michel Marcus, Jul 06 2017
(Python)
from sympy import bernoulli, Integer
def a(n): return bernoulli(n, Integer(2)/5).numerator()
print([a(n) for n in range(51)]) # Indranil Ghosh, Jul 06 2017
CROSSREFS
KEYWORD
sign,frac
AUTHOR
N. J. A. Sloane, Nov 08 2009
STATUS
approved