|
|
A272383
|
|
Numbers n such that Bernoulli number B_{n} has denominator 3318.
|
|
1
|
|
|
78, 1014, 2418, 3354, 7566, 8502, 10842, 11622, 12246, 12714, 13026, 15054, 15366, 15522, 16458, 17394, 23946, 26286, 27222, 27534, 29562, 29874, 30342, 31434, 31902, 33774, 34242, 35646, 36114, 40794, 42198, 43602, 44538, 47814, 48126, 48282, 49218, 50154, 52494, 55302, 57174, 57642, 59046, 59982
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
3318 = 2 * 3 * 7 * 79.
All terms are multiples of a(1) = 78.
For these numbers numerator(B_{n}) mod denominator(B_{n}) = 37.
|
|
LINKS
|
|
|
EXAMPLE
|
Bernoulli B_{78} is 414846365575400828295179035549542073492199375372400483487/3318, hence 78 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, 3318);
|
|
MATHEMATICA
|
Select[78 Range@ 800, Denominator@ BernoulliB@ # == 3318 &] (* Michael De Vlieger, Apr 28 2016 *)
|
|
PROG
|
(PARI) lista(nn) = for(n=1, nn, if(denominator(bernfrac(n)) == 3318, print1(n, ", "))); \\ Altug Alkan, Apr 28 2016
(Python)
from sympy import divisors, isprime
for i in range(78, 10**6, 78):
for d in divisors(i):
if d not in (1, 2, 6, 78) and isprime(d+1):
break
else:
|
|
CROSSREFS
|
Cf. A045979, A051222, A051225, A051226, A051227, A051228, A051229, A051230, A119456, A119480, A249134, A255684, A271634, A271635, A272138, A272139, A272140, A272183, A272184, A272185, A272186.
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|