|
%I #16 May 02 2016 17:35:26
%S 92,184,1564,1748,2116,3496,4232,4324,5428,5612,6532,8648,9476,9844,
%T 10028,10856,11224,12604,14444,15364,16652,18124,18952,19412,20056,
%U 20884,21068,23644,24932,26036,26588,28612,28796,28888,29164,30728,31004,31924,32108,32476,33304,34868,35236,35788,36248,36524
%N Numbers n such that Bernoulli number B_{n} has denominator 1410.
%C 1410 = 2 * 3 * 5 * 47.
%C All terms are multiple of a(1) = 92.
%C For these numbers numerator(B_{n}) mod denominator(B_{n}) = 1333.
%H Chai Wah Wu, <a href="/A272369/b272369.txt">Table of n, a(n) for n = 1..10000</a>
%e Bernoulli B_{92} is -1295585948207537527989427828538576749659341483719435143023316326829946247/1410, hence 92 is in the sequence.
%p with(numtheory): P:=proc(q,h) local n; for n from 2 by 2 to q do
%p if denom(bernoulli(n))=h then print(n); fi; od; end: P(10^6,1410);
%t Select[92 Range@ 360, Denominator@ BernoulliB@ # == 1410 &] (* _Michael De Vlieger_, Apr 28 2016 *)
%o (PARI) lista(nn) = for(n=1, nn, if(denominator(bernfrac(n)) == 1410, print1(n, ", "))); \\ _Altug Alkan_, Apr 28 2016
%o (Python)
%o from sympy import divisors, isprime
%o A272369_list = []
%o for i in range(92, 10**6, 92):
%o for d in divisors(i):
%o if d not in (1,2,4,46) and isprime(d+1):
%o break
%o else:
%o A272369_list.append(i) # _Chai Wah Wu_, May 02 2016
%Y Cf. A045979, A051222, A051225, A051226, A051227, A051228, A051229, A051230, A119456, A119480, A249134, A255684, A271634, A271635, A272138, A272139, A272140, A272183, A272184, A272185, A272186.
%K nonn,easy
%O 1,1
%A _Paolo P. Lava_, Apr 28 2016
%E a(18)-a(29) from _Altug Alkan_, Apr 28 2016
%E More terms from _Michael De Vlieger_, Apr 28 2016
|
