OFFSET
1,3
COMMENTS
First differs from A049865 at n = 24.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = a(A384247(n)) + 1 for n >= 2.
EXAMPLE
n | a(n) | iterations
--+------+----------------------
2 | 1 | 2 -> 1
3 | 2 | 3 -> 2 -> 1
4 | 3 | 4 -> 3 -> 2 -> 1
5 | 4 | 5 -> 4 -> 3 -> 2 -> 1
6 | 2 | 6 -> 2 -> 1
MATHEMATICA
f[p_, e_] := p^e*(1 - 1/p^(2^(IntegerExponent[e, 2]))); iphi[1] = 1; iphi[n_] := iphi[n] = Times @@ f @@@ FactorInteger[n];
a[n_] := Length @ NestWhileList[iphi, n, # != 1 &] - 1; Array[a, 100]
PROG
(PARI) iphi(n) = {my(f = factor(n)); n * prod(i = 1, #f~, (1 - 1/f[i, 1]^(1 << valuation(f[i, 2], 2)))); }
a(n) = if(n == 1, 0, 1 + a(iphi(n)));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Jul 08 2025
STATUS
approved
