OFFSET
0,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = Sum_{k=0..n} A258851(k).
MAPLE
with(numtheory):
d:= n-> n*add(i[2]*pi(i[1])/i[1], i=ifactors(n)[2]):
a:= proc(n) option remember; d(n)+`if`(n>0, a(n-1), 0) end:
seq(a(n), n=0..100);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jun 12 2015
STATUS
approved