OFFSET
1,1
COMMENTS
6 + 345*k and 100 + 345*k are terms for k >= 0.
LINKS
Bernd C. Kellner, The Bernoulli Number Page.
Eric Weisstein's World of Mathematics, Bernoulli Number
Wikipedia, Kummer's congruences
FORMULA
a(n) = A119864(n)/2.
EXAMPLE
Bernoulli(2*6) = -691/2730. So 6 is a term.
MATHEMATICA
Select[Range[4930], Mod[Numerator[BernoulliB[2#]], 691] == 0 &] (* Indranil Ghosh, Mar 11 2017 *)
PROG
(PARI) is(n) = Mod(numerator(bernfrac(2*n)), 691)==0 \\ Felix Fröhlich, Jan 23 2017
(Python)
from itertools import count, islice
from sympy import bernoulli
def A281502gen(): return filter(lambda n:not bernoulli(2*n).p % 691, count(0))
A281502_list = list(islice(A281502gen(), 20)) # Chai Wah Wu, Dec 21 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 23 2017
EXTENSIONS
a(12) - a(36) from Seiichi Manyama, Jan 24 2017
More terms from Indranil Ghosh, Mar 11 2017
STATUS
approved