OFFSET
1,2
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..5000 (first 2000 terms from Paolo P. Lava)
MATHEMATICA
f[x_] := Total@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger[x]]; s = 3; {1, 2}~Join~Reap[Do[s += Sow@ f[s], {n, 3, 1000}] ][[-1, 1]] (* Michael De Vlieger, Jun 24 2025 *)
PROG
(PARI) lista(nn) = {va = vector(nn); print1(va[1] = 1, ", "); print1(va[2] = 2, ", "); sp = vecsum(va); for (k=3, nn, f = factor(sp); va[k] = sum(j=1, #f~, f[j, 1]*f[j, 2]); print1(va[k], ", "); sp += va[k]; ); } \\ Michel Marcus, Feb 19 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
David James Sycamore, Feb 17 2016
STATUS
approved
