OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of A124652(n), n = 1..10^5, showing A124652(a(n)) in green, but A124652(a(n)) that are prime in red.
EXAMPLE
MATHEMATICA
nn = 120; c[_] := False;
rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
f[x_] := Select[Range[x], Divisible[x, rad[#]] &];
Array[Set[{a[#], c[#]}, {#, True}] &, 2]; s = a[1] + a[2];
Reap[Do[r = f[s]; k = SelectFirst[r, ! c[#] &];
If[Divisible[s, k], Sow[i]]; c[k] = True;
s += k, {i, 3, nn}] ][[-1, 1]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, May 05 2024
STATUS
approved