login
A372918
a(n) = Sum_{k=1..n} gcd(k^3,n^2).
1
1, 5, 11, 26, 29, 55, 55, 148, 141, 145, 131, 286, 181, 275, 319, 680, 305, 705, 379, 754, 605, 655, 551, 1628, 1145, 905, 2367, 1430, 869, 1595, 991, 3408, 1441, 1525, 1595, 3666, 1405, 1895, 1991, 4292, 1721, 3025, 1891, 3406, 4089, 2755, 2255, 7480, 4501, 5725, 3355
OFFSET
1,2
LINKS
FORMULA
Multiplicative with a(p^e) = p^e * (1 + ((p-1)/p) * Sum_{i=1..2*e} p^(floor(2*i/3))). - Amiram Eldar, May 24 2024
MATHEMATICA
a[n_] := Sum[GCD[k^3, n^2], {k, 1, n}]; Array[a, 100] (* Amiram Eldar, May 24 2024 *)
PROG
(PARI) a(n) = sum(k=1, n, gcd(k^3, n^2));
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Seiichi Manyama, May 16 2024
STATUS
approved