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

0, 0, 1, 3, 5, 19, 19, 17, 89, 673, 47, 979, 1297, 4883, 1771, 79613, 31, 393959, 2033, 85639, 116551, 616181, 4111, 16637083, 727403, 13117673, 72631, 122771983, 194803, 31691158757, 491951, 124085749, 9079549, 114103102711, 92671, 743246297281, 213649, 197986199, 972486919, 144015774883

Offset

1,4

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) = 19.

Maple

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

Mathematica

a[n_] := Sum[1/(n-p), {p, Prime[Range[PrimePi[n-1]]]}] // Numerator;
Array[a, 100] (* Jean-François Alcover, Apr 29 2019 *)

Prog

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

Crossrefs

Cf. A305702 (denominators).

Sequence in context: A270723 A242961 A025046 * A299073 A329797 A180931

Adjacent sequences: A305697 A305698 A305699 * A305701 A305702 A305703

Keyword

nonn,frac

Author

Robert Israel, Jun 08 2018

Status

approved