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-3 of 3 results.

A291789 Trajectory of 270 under repeated application of k -> (phi(k)+sigma(k))/2.

Original entry on oeis.org

270, 396, 606, 712, 851, 852, 1148, 1416, 2032, 2488, 2960, 4110, 5512, 6918, 8076, 10780, 16044, 23784, 33720, 55240, 73230, 97672, 118470, 169840, 247224, 350260, 442848, 728448, 1213440, 2124864, 4080384, 8159616, 13515078, 15767596, 18626016, 29239504, 39012864, 62623600, 92580308
Offset: 0

Views

Author

N. J. A. Sloane, Sep 03 2017

Keywords

Comments

The ultimate fate of this trajectory is presently unknown. It may reach a fractional value (when it dies), it may reach a prime (which would be a fixed point), it may enter a cycle of length greater than 1, or it may be unbounded. - Hugo Pfoertner and N. J. A. Sloane, Sep 18 2017

Crossrefs

See A291914 and A292108 for the "big picture".

Programs

  • Maple
    orbit:= proc(n, m) uses numtheory;
      local V,k;
      V:= Vector(m);
      V[1]:= n;
      for k from 2 to m do V[k]:= (sigma(V[k-1])+ phi(V[k-1]))/2 od:
      convert(V,list)
    end proc:
    orbit(270, 200); # Robert Israel, Sep 07 2017
  • Mathematica
    NestWhileList[If[! IntegerQ@ #, -1/2, (DivisorSigma[1, #] + EulerPhi@ #)/2] &, 270, Nor[! IntegerQ@ #, SameQ@ ##] &, 2, 38] (* Michael De Vlieger, Sep 19 2017 *)

A292108 Iterate the map k -> (sigma(k) + phi(k))/2 starting at n; a(n) is the number of steps to reach either a fixed point or a fraction, or a(n) = -1 if neither of these two events occurs.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 3, 2, 1, 0, 1, 0, 2, 2, 1, 0, 4, 1, 2, 1, 4, 0, 2, 0, 1, 4, 3, 2, 1, 0, 3, 2, 1, 0, 9, 0, 2, 3, 1, 0, 7, 1, 1, 2, 1, 0, 8, 3, 2, 2, 1, 0, 3, 0, 8, 7, 1, 3, 2, 0, 1, 7, 6, 0, 1, 0, 3, 2, 4
Offset: 1

Views

Author

Hugo Pfoertner and N. J. A. Sloane, Sep 18 2017

Keywords

Comments

The first unknown value is a(270).
For an alternative version of this sequence, see A291914.
From Andrew R. Booker, Sep 19 2017 and Oct 03 2017: (Start)
Let f(n) = (sigma(n) + phi(n))/2. Then f(n) >= n, so the trajectory of n under f either terminates with a half-integer, reaches a fixed point, or increases monotonically. The fixed points of f are 1 and the prime numbers, and f(n) is fractional iff n>2 is a square or twice a square.
It seems likely that a(n) = -1 for all but o(x) numbers n <= x. See link for details of the argument. (End)

Examples

			Let f(k) = (sigma(k) + phi(k))/2. Under the action of f:
14 -> 15 -> 16 -> 39/2, taking 3 steps, so a(14) = 3.
21 -> 22 -> 23, a prime, in 2 steps, so a(21) = 2.
		

Crossrefs

Programs

  • Mathematica
    With[{i = 200}, Table[-1 + Length@ NestWhileList[If[! IntegerQ@ #, -1/2, (DivisorSigma[1, #] + EulerPhi@ #)/2] &, n, Nor[! IntegerQ@ #, SameQ@ ##] &, 2, i, -1] /. k_ /; k >= i - 1 -> -1, {n, 76}]] (* Michael De Vlieger, Sep 19 2017 *)

Formula

a(n) = 0 if n is 1 or a prime (these are fixed points).
a(n) = 1 if n>2 is a square or twice a square, since these reach a fraction in one step.

A291805 a(n) = phi(A291789(n)).

Original entry on oeis.org

72, 120, 200, 352, 792, 280, 480, 464, 1008, 1240, 1152, 1088, 2496, 2304, 2688, 3360, 4560, 7920, 8960, 22080, 19520, 47040, 28640, 61440, 82400, 137760, 126336, 207360, 319488, 552960, 1013760, 2716416, 4505024, 7193472
Offset: 1

Views

Author

N. J. A. Sloane, Sep 07 2017

Keywords

Comments

Created in the hope that it will help explain the long-term behavior of A291789, which is the mean of this and A291804.

Crossrefs

Showing 1-3 of 3 results.