A360806 -id:A360806 - cp's OEIS frontend

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

Ctibor O. Zizka, Jul 11 2024

nonn

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.

			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.

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

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 *)

A375147 a(1) = 0; for n >= 2, a(n) is the number of iterations needed for the map: x -> x / A000005(x) if x is divisible by A000005(x), x -> x + 1 otherwise, to reach 1.

Original entry on oeis.org

0, 1, 7, 6, 5, 4, 3, 2, 8, 4, 3, 2, 13, 12, 11, 10, 9, 8, 13, 12, 11, 10, 9, 8, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 9, 8, 7, 6, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 9, 8, 7, 6, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 10, 9, 8, 7, 6, 5, 4, 3, 7, 6, 5, 4

Offset: 1

Ctibor O. Zizka, Aug 01 2024

nonn

The trajectory length is a repeated sum of steps up to the next refactorable number (A360778) and its refactoring "depth" (A374540). The sequence always reach 1 as soon as an iterate reaches the value x from A330816. Assuming A330816 to be finite (conjectured by David A. Corneth) and A360806 to be infinite, may there be a set of numbers n > 10^42, which is not reaching 1 ?

			x = 3:  the trajectory is 3 --> 4 --> 5 --> 6 --> 7 --> 8 --> 2 --> 1, number of steps needed to reach 1 is 7, thus a(3) = 7.
x = 81: the trajectory is 81 --> 82 --> 83 --> 84 --> 7 --> 8 --> 2 --> 1, number of steps needed to reach 1 is 7, thus a(81) = 7.

Cf. A000005, A033950, A330816, A360778, A360806, A374540.

a[n_] := -1 + Length[NestWhileList[If[IntegerQ[(r = #/DivisorSigma[0, #])], r, # + 1] &, n, # > 1 &]]; Array[a, 100] (* Amiram Eldar, Aug 01 2024 *)

a(A360806(n)) = n.

cp's OEIS Frontend

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).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica