OFFSET
1,2
FORMULA
a(n) = Sum_{d|n} A053820(d).
MATHEMATICA
Table[Sum[(k/GCD[n, k])^4, {k, 1, n}], {n, 1, 39}]
PROG
(PARI) a(n) = sum(k=1, n, (k/gcd(n, k))^4); \\ Michel Marcus, Apr 17 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 17 2021
STATUS
approved