OFFSET
1,1
COMMENTS
Satisfies a(n) = 2*a(n-3) - a(n-6) for n < 67. - Chai Wah Wu, May 28 2016
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..6450 (terms 1..1000 from Joerg Arndt)
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
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Feb 25 2004
STATUS
approved