OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..387
László Tóth, On the number of cyclic subgroups of a finite abelian group, arXiv: 1203.6201 [math.GR], 2012.
FORMULA
a(n) = Sum_{x_1|n, x_2|n, ... , x_n|n} phi(x_1)*phi(x_2)* ... *phi(x_n)/phi(lcm(x_1, x_2, ... , x_n)).
a(n) = Sum_{d|n} b(d, n), where b(n, k) = ( Sum_{d|n} mu(n/d) * d^k )/phi(n).
If p is prime, a(p) = 1 + (p^p - 1)/(p - 1).
MATHEMATICA
b[n_, k_] := DivisorSum[n, MoebiusMu[n/#] * #^k &] / EulerPhi[n]; a[n_] := DivisorSum[n, b[#, n] &]; Array[a, 20] (* Amiram Eldar, Oct 04 2023 *)
PROG
(PARI) b(n, k) = sumdiv(n, d, moebius(n/d)*d^k)/eulerphi(n);
a(n) = sumdiv(n, d, b(d, n));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 12 2021
STATUS
approved