cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A334814 Least number that reaches 1 after n iterations of the map k -> sigma(k)/d(k) if d(k) | sigma(k), and k -> 1 otherwise, where d(k) is the number of divisors of k (A000005) and sigma(k) is their sum (A000203).

Original entry on oeis.org

1, 2, 3, 5, 11, 29, 107, 257, 941, 2017, 11261, 45039, 441073, 2151073, 8575873, 42884161, 220268161, 440536321
Offset: 0

Views

Author

Amiram Eldar, May 12 2020

Keywords

Comments

Apparently, most of the terms are primes. 45039 = 3 * 15013 is the first composite term.
a(18) > 2*10^10, if it exists.

Examples

			a(3) = 5 since sigma(5)/d(5) = 6/2 = 3, sigma(3)/d(3) = 4/2 = 2, and sigma(2)/d(2) = 3/2 is not an integer, hence there are 3 iterations: 5 -> 3 -> 2 -> 1, and 5 is the least number with 3 iterations.
		

Crossrefs

Programs

  • Mathematica
    rat[n_] := If[IntegerQ[r = DivisorSigma[1, n]/DivisorSigma[0, n]], r, 1]; f[n_] := Length @ FixedPointList[rat, n] - 1; max = 10; seq = Table[0, {max}]; c = 0; n = 1; While[c < max, i = f[n]; If[i <= max && seq[[i]] == 0, c++; seq[[i]] = n]; n++]; seq

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

Views

Author

Ctibor O. Zizka, Aug 01 2024

Keywords

Comments

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 ?

Examples

			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.
		

Crossrefs

Programs

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

Formula

a(A360806(n)) = n.
Showing 1-2 of 2 results.