OFFSET
0,1
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
LINKS
Sean A. Irvine, Table of n, a(n) for n = 0..515 [Terms through a(250) from Hugo Pfoertner, terms a(251)-a(356) from N. J. A. Sloane]
Sean A. Irvine, Illustration of A291789 showing a(n+1)/a(n) (red), cumulative mean of a(n+1)/a(n) (green), and power of 2 in a(n) (blue)
N. J. A. Sloane, Three (No, 8) Lovely Problems from the OEIS, Experimental Mathematics Seminar, Rutgers University, Oct 05 2017, Part I, Part 2, Slides. (Mentions this sequence)
N. J. A. Sloane, "A Handbook of Integer Sequences" Fifty Years Later, arXiv:2301.03149 [math.NT], 2023, p. 14.
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 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 03 2017
STATUS
approved