login
Number of iterations of the number of divisors function (A000005) required to reach a fixed point (1 or 2) when started at n!.
3

%I #20 Aug 17 2024 09:01:26

%S 1,1,4,5,4,6,7,7,7,5,7,7,7,7,7,6,7,7,7,7,7,8,7,7,7,7,6,6,8,6,7,6,7,8,

%T 6,8,7,5,7,6,8,6,6,8,8,8,8,8,8,6,8,8,8,8,8,8,8,8,8,8,7,8,7,8,7,8,8,8,

%U 8,8,8,7,8,8,8,8,8,8,7,8,8,8,8,8,8,8,8,8,8,8,7,8,7,8,8,8,7,8,8,7,8,8,8,8,8

%N Number of iterations of the number of divisors function (A000005) required to reach a fixed point (1 or 2) when started at n!.

%H Amiram Eldar, <a href="/A053025/b053025.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A036459(A000142(n)) = A036459(n!).

%e For n = 108, a(108) = 9 because the sequence of iterates is {108!, 798687560466432000, 7920, 60, 12, 6, 4, 3, 2}, and its length is 9.

%t a[n_] := -1 + Length @ FixedPointList[DivisorSigma[0, #] &, n!]; Array[a, 100] (* _Amiram Eldar_, Aug 17 2024 *)

%o (PARI) a(n) = {my(f = n!, c = 1); while(f > 2, f = numdiv(f); c++); c;} \\ _Amiram Eldar_, Aug 17 2024

%Y Cf. A000005, A036459, A000142.

%K nonn

%O 1,3

%A _Labos Elemer_, Feb 24 2000