OFFSET
1,10
COMMENTS
a(n) is the number of primes p >= 2n+3 that divide the numerator of the Bernoulli number B_{2n}.
LINKS
B. Mazur, How can we construct abelian Galois extensions of basic number fields?, Bull. Amer. Math. Soc., 48 (2011), 155-209. See footnote 64 on p. 205.
Samual S. Wagstaff, Jr., Prime divisors of the Bernoulli and Euler numbers
Wikipedia, Irregular pairs
MATHEMATICA
Table[p = Select[First /@ FactorInteger[Abs[Numerator[BernoulliB[2n]]]], # >= 2*n+3 &]; Length[p], {n, 35}] (* T. D. Noe, Apr 25 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Sondow, Apr 25 2011
STATUS
approved