OFFSET
1,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = (p^(2*e+1) - (p + 1) * p^(2*e-1) + 1) / (p + 1).
Sum_{k=1..n} a(k) ~ c * n^3, where c = (1/3) * Product_{p prime} (1 - 1/p^2 - 1/(1 + p + p^2)) = 0.1381393084... . - Amiram Eldar, Nov 18 2022
MATHEMATICA
Table[Sum[If[GCD[n, k] == 1, n/GCD[n, k - 1], 0], {k, n}], {n, 60}]
f[p_, e_] := (p^(2 e + 1) - (p + 1) p^(2 e - 1) + 1)/(p + 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 60]
PROG
(PARI) a(n) = sum(k=1, n, if (gcd(n, k)==1, n/gcd(n, k-1))); \\ Michel Marcus, Sep 27 2021
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Ilya Gutkovskiy, Sep 26 2021
STATUS
approved