OFFSET
1,84
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
MATHEMATICA
f[p_, e_] := p^e - 1; uphi[1] = 1; uphi[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := Plus @@ DeleteDuplicates[uphi /@ Select[Divisors[n], CoprimeQ[#, n/#] &]]; a[n_] := -2 + Length@ FixedPointList[s, n]; Array[a, 100]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 09 2021
STATUS
approved