OFFSET
0,2
COMMENTS
See, e.g., A157871 for details on B[d,a](n) with gcd(d,a) = 1.
LINKS
Antti Karttunen, Table of n, a(n) for n = 0..10101
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.
MATHEMATICA
Table[Denominator[BernoulliB[n, 1/5]]/5^n, {n, 0, 70}] (* Jean-François Alcover, Sep 24 2018, from PARI *)
PROG
(PARI) a(n)=denominator(subst(bernpol(n, x), x, 1/5))/5^n; \\ Michel Marcus, Jul 06 2017
(Python)
from sympy import bernoulli
def a(n): return bernoulli(n, 1/Integer(5)).denominator()//(5**n)
print([a(n) for n in range(41)]) # Indranil Ghosh, Jul 06 2017
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Wolfdieter Lang, Jul 05 2017
STATUS
approved