A305702 a(n) is the denominator of Sum_{primes p < n} 1/(n-p).

1, 1, 1, 2, 6, 12, 20, 10, 84, 840, 72, 630, 1320, 2772, 1560, 90090, 42, 240240, 1904, 46410, 95760, 639540, 5040, 9699690, 637560, 14316120, 92400, 176125950, 308880, 20078358300, 475020, 67690350, 7447440, 116925953760, 110880, 1002802450650, 318240, 122453100, 906665760, 165723680430

Offset

1,4

Comments

The first n >= 2 for which a(n) <> A130492(n-1) is 32.

Links

Robert Israel, Table of n, a(n) for n = 1..3569

Math Overflow, Sum of reciprocals of integers minus primes.

Example

Sum_{primes p < 6} 1/(6-p) = 1/(6-2) + 1/(6-3) + 1/(6-5) = 19/12 so a(6) = 12.

Maple

N:= 100: # to get a(1)..a(N)
P:= select(isprime, [2, seq(i, i=3..N, 2)]):
seq(denom(add(1/(n-p), p=select(`<`, P, n))), n=1..N);

Prog

(PARI) a(n) = my(p=select(x->isprime(x), [1..n-1])); denominator(sum(k=1, #p, 1/(n-p[k]))); \\ Michel Marcus, Jun 09 2018

Crossrefs

Cf. A130492, A305700.

Sequence in context: A139084 A086958 A130492 * A152222 A139082 A061078

Adjacent sequences: A305699 A305700 A305701 * A305703 A305704 A305705

Keyword

nonn,frac

Author

Robert Israel, Jun 08 2018

Status

approved