A295591 Numbers k such that Bernoulli number B_{k} has denominator 61410.

88, 968, 5192, 5368, 13816, 15928, 19624, 19976, 22616, 23144, 23848, 24904, 27368, 27544, 27896, 29656, 31064, 33704, 34936, 38632, 40216, 40568, 40744, 45848, 46024, 48136, 49544, 50248, 51656, 53416, 56584, 56936, 57112, 59048, 60808, 61688, 67672, 68024, 71368

Offset

1,1

Comments

61410 = 2*3*5*23*89.

All terms are multiples of a(1) = 88.

For these numbers numerator(B_{k}) mod denominator(B_{k}) = 56003.

Links

Seiichi Manyama, Table of n, a(n) for n = 1..1000

Example

Bernoulli B_{88} is -1311426488674017507995511424019311843345750275572028644296919890574047/61410 hence 88 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, 61410);
# Alternative: # according to Robert Israel code in A282773
with(numtheory): filter:= n ->
select(isprime, map(`+`, divisors(n), 1)) = {2, 3, 5, 23, 89}:
select(filter, [seq(i, i=1..10^5)]);

Prog

(PARI) isok(n) = denominator(bernfrac(n)) == 61410; \\ Michel Marcus, Jan 07 2018

Crossrefs

Cf. A045979, A051222, A051225, A051226, A051227, A051228, A051229, A051230, A119456, A119480, A249134, A255684, A271634, A271635, A272138, A272139, A272140, A272183, A272184, A272185, A272186, A272369.

Sequence in context: A055749 A227479 A321061 * A182676 A228744 A249294

Adjacent sequences: A295588 A295589 A295590 * A295592 A295593 A295594

Keyword

nonn,easy

Author

Paolo P. Lava, Nov 24 2017

Status

approved