OFFSET
1,2
FORMULA
a(n) = Sum_{d|n} d*A002322(n/d).
MATHEMATICA
a[n_] := DivisorSum[n, # * CarmichaelLambda[n/#] &]; Array[a, 50] (* Amiram Eldar, Jan 13 2020 *)
PROG
(PARI) a(n) = sumdiv(n, d, d * lcm(znstar(n/d)[2])); \\ Daniel Suteu, Jan 13 2020
(Magma) [1] cat [&+[d*CarmichaelLambda(n div d):d in Set(Divisors(n)) diff {n}]+n:n in [2..60]]; // Marius A. Burtea, Jan 14 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Torlach Rush, Jan 13 2020
STATUS
approved