OFFSET
1,1
COMMENTS
272118 = 2*3*7*11*19*31.
All terms are multiples of a(1) = 90.
For these numbers numerator(B_{k}) mod denominator(B_{k}) = 230759.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..1000
EXAMPLE
272118 = 2*3*7*11*19*31.
Bernoulli B_{90} is 1179057279021082799884123351249215083775254949669647116231545215727922535/ 272118 hence 90 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, 272118);
# Alternative: according to Robert Israel code in A282773
with(numtheory): filter:= n ->
select(isprime, map(`+`, divisors(n), 1)) = {2, 3, 7, 11, 19, 31}:
select(filter, [seq(i, i=1..10^5)]);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Nov 24 2017
STATUS
approved
