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

Author

N. J. A. Sloane, Feb 24 2004

Keywords

• nonn

Links

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

Crossrefs

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

Programs

Mathematica

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

PARI

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

Extensions

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