OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(9) = 22; the primes dividing 10 are 2 and 5, but 2 also divides 22, so a(10) = 22 + 5 = 27.
MAPLE
R:=1: v:= 1:
for n from 2 to 100 do
v:= v + convert(select(p -> v mod p <> 0, numtheory:-factorset(n)), `+`);
R:= R, v
od:
R;
MATHEMATICA
n = 1; NestList[Function[k, n++; k + Total@ Select[FactorInteger[n][[All, 1]], CoprimeQ[#, k] &]], 1, 63] (* Michael De Vlieger, Jun 28 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Jun 27 2024
STATUS
approved