OFFSET
1,2
COMMENTS
Inverse Möbius transform of mu(n)^2 * (n mod 10). - Wesley Ivan Hurt, Jun 21 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = Sum_{d|n} mu(d)^2 * (d mod 10).
EXAMPLE
a(10) = 8. The squarefree divisors of 10 are 1, 2, 5, 10 and the sum of their final digits is 1 + 2 + 5 + 0 = 8.
MAPLE
f:= proc(n) local t; add(t mod 10, t = map(convert, combinat:-powerset(numtheory:-factorset(n)), `*`)) end proc:
map(f, [$1..100]); # Robert Israel, Nov 21 2023
MATHEMATICA
Table[DivisorSum[n, MoebiusMu[#]^2*Mod[#, 10] &], {n, 100}]
PROG
(PARI) a(n) = sumdiv(n, d, if (issquarefree(d), d%10)); \\ Michel Marcus, Nov 21 2023
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Wesley Ivan Hurt, Nov 20 2023
STATUS
approved