OFFSET
1,2
MATHEMATICA
With[{max = 60}, sp = Accumulate[Array[# * Boole@PrimeQ[#]&, max]]; a[n_] := DivisorSum[n, sp[[#]] &]; Array[a, max]] (* Amiram Eldar, Aug 22 2023 *)
PROG
(PARI) vector(80, n, sumdiv(n, d, sum(i=1, primepi(d), prime(i))))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Feb 21 2009
STATUS
approved