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.

A379146 Numbers k that are not in A378930 (i.e., are never the value of f(n) = n * d(n) / gcd(n, d(n))^2, where d = A000005).

Original entry on oeis.org

18, 27, 45, 63, 64, 72, 99, 112, 117, 144, 153, 160, 171, 207, 225, 243, 252, 261, 279, 288, 320, 333, 336, 352, 360, 369, 387, 396, 416, 423, 441, 468, 477, 504, 531, 544, 549, 567, 576, 603, 608, 612, 616, 625, 639, 657, 684, 711, 728, 736, 747, 792, 801, 828, 873, 880, 891, 909, 927, 928, 936, 952, 963, 981, 992
Offset: 1

Views

Author

Viliam Furík, Dec 16 2024

Keywords

Comments

Verified using the known values of f(n) up to a limit determined by the upper bound for divisor function 2*sqrt(n). The lower bound for f(n) is n/d(n), which can be combined with d(n) <= 2*sqrt(n) to yield f(n) >= sqrt(n)/2, and n <= 4*f(n)^2. E.g. for f(n) = 18, this means that checking f(n) for n <= 1296 is sufficient to verify it's never the value of f(n).
It appears that all numbers 9*p, p prime, are in this sequence.

Crossrefs

Cf. A000005. Complement of A378930.

A378929 Number of steps it takes for the chain of f(n) = n * d(n) / gcd(n, d(n))^2 to reach a cycle from a starting number n.

Original entry on oeis.org

0, 1, 1, 3, 1, 0, 1, 2, 2, 0, 1, 2, 1, 0, 3, 4, 1, 2, 1, 4, 3, 0, 1, 2, 1, 0, 4, 4, 1, 3, 1, 6, 3, 0, 1, 4, 1, 0, 3, 2, 1, 3, 1, 4, 4, 0, 1, 5, 1, 0, 3, 4, 1, 4, 1, 2, 3, 0, 1, 2, 1, 0, 4, 8, 1, 3, 1, 4, 3, 0, 1, 1, 1, 0, 0, 4, 1, 3, 1, 3, 1, 0, 1, 2, 1, 0, 3, 2
Offset: 1

Views

Author

Viliam Furík, Dec 11 2024

Keywords

Comments

It is conjectured that every number n eventually ends up in a cycle.

Crossrefs

Programs

  • Mathematica
    f[n_] := Module[{d = DivisorSigma[0, n]}, n*d/GCD[n, d]^2]; a[n_] := Module[{s = NestWhileList[f, n, UnsameQ, All]}, FirstPosition[s, s[[-1]]][[1]] - 1]; Array[a, 100] (* Amiram Eldar, Dec 16 2024 *)
Showing 1-2 of 2 results.