|
|
A092221
|
|
Numbers k such that numerator of Bernoulli(2*k) is divisible by 59, the second irregular prime.
|
|
14
|
|
|
22, 51, 59, 80, 109, 118, 138, 167, 177, 196, 225, 236, 254, 283, 295, 312, 341, 354, 370, 399, 413, 428, 457, 472, 486, 515, 531, 544, 573, 590, 602, 631, 649, 660, 689, 708, 718, 747, 767, 776, 805, 826, 834, 863, 885, 892, 921, 944, 950, 979, 1003, 1008
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Satisfies a(n) = 2*a(n-3) - a(n-6) for n < 67. - Chai Wah Wu, May 28 2016
|
|
LINKS
|
|
|
MATHEMATICA
|
Select[ Range[ 1036], Mod[ Numerator[ BernoulliB[2# ]], 59] == 0 &]
|
|
PROG
|
(Python)
from sympy import bernoulli
A092221_list = [n for n in range(10**3) if not bernoulli(2*n).p % 59] # Chai Wah Wu, May 28 2016
(PARI) for(n=0, 10^3, if( numerator(bernfrac(2*n))%59==0, print1(n, ", ") ) ); \\ Joerg Arndt, May 29 2016
|
|
CROSSREFS
|
Cf. A000928, A091216, A092222, A092223, A092224, A092225, A092226, A092227, A092228, A092229, A092231.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|