OFFSET
1,1
COMMENTS
See A189683 for pairs (p,2n) for the primes p in this sequence.
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..154 (first 66 terms from T. D. Noe)
Hisanori Mishima, WIFC (World Integer Factorization Center)
MATHEMATICA
Flatten[Table[Select[First /@ FactorInteger[Abs[Numerator[BernoulliB[n]]]], # >= n+3 &], {n, 2, 70, 2}]] (* T. D. Noe, Apr 25 2011 *)
PROG
(Macsyma) for n do for p in map('first, factor_number(abs(num(bern(2*n))))) do if p>=2*n+3 then (?prin1(p), ?prin1('?\-));
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Definition modified by Jonathan Sondow, Apr 27 2011
STATUS
approved