OFFSET
1,2
MAPLE
a:= n-> add(n^(k*add(i[2]/i[1], i=ifactors(k)[2])), k=1..n):
seq(a(n), n=1..18); # Alois P. Heinz, Oct 15 2021
MATHEMATICA
d[0] = d[1] = 0; d[n_] := n * Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); a[n_] := Sum[n^d[k], {k, 1, n}]; Array[a, 20] (* Amiram Eldar, Oct 16 2021 *)
PROG
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
a(n) = sum(k=1, n, n^ad(k)); \\ Michel Marcus, Oct 18 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Oct 15 2021
STATUS
approved