OFFSET
1,1
COMMENTS
Credit is given for Zak Seidov for doing the calculations.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..500
MAPLE
N:= 10^5: # to get all terms before the first term > N
A007504:= ListTools:-PartialSums([seq(ithprime(i), i=1..N)]):
found:= true:
for n from 1 while found do
found:= false;
for k from n+1 to N do
found:= true;
A[n]:= k;
break
fi
od
od:
seq(A[i], i=1..n-2); # Robert Israel, May 02 2018
MATHEMATICA
a[n_] := For[tpn = Prime[Range[n]] // Total; tpk = 0; k = n+1, True, k++, tpk += Prime[k]; If[Divisible[tpk, tpn], Return[k]]]; Array[a, 50] (* Jean-François Alcover, Mar 22 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot, Apr 27 2018
EXTENSIONS
More terms from Alois P. Heinz, Apr 27 2018
STATUS
approved