|
%I #19 Mar 29 2021 08:02:12
%S 16,34,37,52,70,74,88,106,111,124,142,148,160,178,185,196,214,222,232,
%T 250,259,268,286,296,304,322,333,340,358,370,376,394,407,412,430,444,
%U 448,466,481,484,502,518,520,538,555,556,574,592,610,628
%N Numbers k such that numerator of Bernoulli(2*k) is divisible by 37, the first irregular prime.
%H Amiram Eldar, <a href="/A091216/b091216.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..100 from Seiichi Manyama)
%t Select[ Range[ 681], Mod[ Numerator[ BernoulliB[2# ]], 37] == 0 &] (* _Robert G. Wilson v_, Feb 24 2004 *)
%o (PARI) for(j=1,260, if (! (numerator(bernfrac(2*j))%37), print1(j, ", ")))
%Y Cf. A000928, A092221, A092222, A092223, A092224, A092225, A092226, A092227, A092228, A092229, A092230, A092231.
%K nonn
%O 1,1
%A _N. J. A. Sloane_, Feb 24 2004
%E More terms from _Robert G. Wilson v_, Feb 24 2004
|
