A374540 a(1) = 0; for n >= 2, a(n) is the number of iterations needed for the map x -> x/A000005(x) to reach a least integer, when starting from x = A033950(n).

0, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1

Offset

1,3

Comments

The refactorability "depth" for refactorable numbers. Numbers from A159973 have the refactorability "depth" 0. Records reached for A033950(A360806(n)), i.e. the growth of the sequence is very slow.

Links

Table of n, a(n) for n=1..97.

Example

n = 2: A033950(2) = 2, 2/A000005(2) = 1, thus a(2) = 1.
n = 3: A033950(3) = 8, 8/A000005(8) = 2 --> 2/A000005(2) = 1, thus a(3) = 2.
n = 13: A033950(13) = 80, 80/A000005(80) = 8 --> 8/A000005(8) = 2 --> 2/A000005(2) = 1, thus a(13) = 3.

Mathematica

f[n_] := Module[{v = NestWhileList[# / DivisorSigma[0, #] &, n, IntegerQ[#] && # > 1 &], len}, len = Length[v]; If[IntegerQ[v[[2]]], If[v[[-1]] == 1, len - 1, len - 2], Nothing]]; f[1] = 0; Array[f, 1200] (* Amiram Eldar, Jul 11 2024 *)

Crossrefs

Cf. A000005, A033950, A036762, A159973, A360806.

Sequence in context: A030600 A079700 A357722 * A135362 A373387 A309913

Adjacent sequences: A374537 A374538 A374539 * A374541 A374542 A374543

Keyword

nonn

Author

Ctibor O. Zizka, Jul 11 2024

Status

approved