OFFSET
1,2
FORMULA
a(n) = n * Sum_{d|n, gcd(d,n/d)=1} mu(d)^2 / d.
MATHEMATICA
a[n_] := DivisorSum[n, n/# &, CoprimeQ[#, n/#] && SquareFreeQ[#] &]; Array[a, 100] (* Amiram Eldar, May 26 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, if ((gcd(d, n/d) == 1) && issquarefree(d), n/d)); \\ Michel Marcus, May 26 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wesley Ivan Hurt, May 26 2023
STATUS
approved