OFFSET
1,1
COMMENTS
The subsequence of irregular primes p is A046753.
LINKS
T. D. Noe, Table of n, a(n) for n = 1..132
Bernd C. Kellner, On irregular prime power divisors of the Bernoulli numbers, Math. Comp. 76 (2007), 405-441; arXiv:0409223 [math.NT], 2004.
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.
Wikipedia, Regular prime: Irregular pairs.
EXAMPLE
The first few irregular pairs are (691,12), (3617,16), (43867,18), (283,20), (617,20), (131,22), (593,22), ...
MATHEMATICA
Flatten[Table[p = Select[First /@ FactorInteger[Abs[Numerator[BernoulliB[n]]]], # >= n+3 &]; Transpose[{p, Table[n, {Length[p]}]}], {n, 2, 70, 2}]] (* T. D. Noe, Apr 25 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Sondow, Apr 25 2011
STATUS
approved