login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A092221 Numbers n such that numerator of Bernoulli(2n) 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

Joerg Arndt, Table of n, a(n) for n = 1..1000

Eric Weisstein's World of Mathematics, Bernoulli Number..

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.

Sequence in context: A043168 A043948 A015226 * A191279 A200880 A111576

Adjacent sequences:  A092218 A092219 A092220 * A092222 A092223 A092224

KEYWORD

nonn

AUTHOR

Robert G. Wilson v, Feb 25 2004

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified February 18 05:48 EST 2018. Contains 299298 sequences. (Running on oeis4.)