OFFSET
1,2
COMMENTS
a(51) > 1.5*10^13. - Bruce Garner, Jun 02 2021
LINKS
Bruce Garner, Table of n, a(n) for n = 1..50 (first 42 terms from Robert Price)
EXAMPLE
6 is a term because 1 plus the sum of the first 6 primes^19 is 1523090798793695143992 which is divisible by 6.
MATHEMATICA
p = 2; k = 0; s = 1; lst = {}; While[k < 40000000000, s = s + p^19; If[Mod[s, ++k] == 0, AppendTo[lst, k]; Print[{k, p}]]; p = NextPrime@ p] (* derived from A128169 *)
Module[{nn=74*10^5, apr}, apr=Accumulate[Prime[Range[nn]]^19]; Select[Range[ nn], Divisible[1+apr[[#]], #]&]] (* The program generates the first 25 terms of the sequence. To generate more, increase the value of nn, but the program may take a long time to run. *) (* Harvey P. Dale, Oct 02 2021 *)
CROSSREFS
Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).
KEYWORD
nonn
AUTHOR
Robert Price, Dec 15 2013
STATUS
approved