OFFSET
1,1
COMMENTS
1410 = 2 * 3 * 5 * 47.
All terms are multiple of a(1) = 92.
For these numbers numerator(B_{n}) mod denominator(B_{n}) = 1333.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
EXAMPLE
Bernoulli B_{92} is -1295585948207537527989427828538576749659341483719435143023316326829946247/1410, hence 92 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, 1410);
MATHEMATICA
Select[92 Range@ 360, Denominator@ BernoulliB@ # == 1410 &] (* Michael De Vlieger, Apr 28 2016 *)
PROG
(PARI) lista(nn) = for(n=1, nn, if(denominator(bernfrac(n)) == 1410, print1(n, ", "))); \\ Altug Alkan, Apr 28 2016
(Python)
from sympy import divisors, isprime
A272369_list = []
for i in range(92, 10**6, 92):
for d in divisors(i):
if d not in (1, 2, 4, 46) and isprime(d+1):
break
else:
A272369_list.append(i) # Chai Wah Wu, May 02 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Apr 28 2016
EXTENSIONS
a(18)-a(29) from Altug Alkan, Apr 28 2016
More terms from Michael De Vlieger, Apr 28 2016
STATUS
approved