A091216 Numbers k such that numerator of Bernoulli(2*k) is divisible by 37, the first irregular prime.

16, 34, 37, 52, 70, 74, 88, 106, 111, 124, 142, 148, 160, 178, 185, 196, 214, 222, 232, 250, 259, 268, 286, 296, 304, 322, 333, 340, 358, 370, 376, 394, 407, 412, 430, 444, 448, 466, 481, 484, 502, 518, 520, 538, 555, 556, 574, 592, 610, 628

Offset

1,1

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..100 from Seiichi Manyama)

Mathematica

Select[ Range[ 681], Mod[ Numerator[ BernoulliB[2# ]], 37] == 0 &] (* Robert G. Wilson v, Feb 24 2004 *)

Prog

(PARI) for(j=1, 260, if (! (numerator(bernfrac(2*j))%37), print1(j, ", ")))

Crossrefs

Cf. A000928, A092221, A092222, A092223, A092224, A092225, A092226, A092227, A092228, A092229, A092230, A092231.

Sequence in context: A132370 A185467 A245664 * A350522 A132760 A209377

Adjacent sequences: A091213 A091214 A091215 * A091217 A091218 A091219

Keyword

nonn

Author

N. J. A. Sloane, Feb 24 2004

Extensions

More terms from Robert G. Wilson v, Feb 24 2004

Status

approved