OFFSET
1,1
COMMENTS
282 = 2 * 3 * 47.
All terms are multiples of a(1) = 46.
For these numbers numerator(B_{n}) mod denominator(B_{n}) = 41.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..1000
EXAMPLE
Bernoulli B_{46} is 596451111593912163277961/282, hence 46 is in the sequence.
MAPLE
with(numtheory): P:=proc(q, h) local n; for n from 2 by 2 to q do
if denom(bernoulli(n))=h then print(n); fi; od; end: P(10^6, 282);
MATHEMATICA
Select[46 Range@ 320, Denominator@ BernoulliB@ # == 282 &] (* Michael De Vlieger, Apr 29 2016 *)
PROG
(PARI) isok(n) = denominator(bernfrac(n)) == 282; \\ Michel Marcus, Apr 22 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Apr 22 2016
EXTENSIONS
a(13)-a(28) from Michel Marcus, Apr 22 2016
a(29)-a(41) from Altug Alkan, Apr 22 2016
STATUS
approved