OFFSET
1,2
PROG
(PARI) a(n) = sumdiv(n, d, lcm(znstar(d)[2]) * lcm(znstar(n/d)[2]));
CROSSREFS
KEYWORD
nonn
AUTHOR
Torlach Rush, Jan 14 2020
STATUS
approved