|
%I #7 Mar 29 2015 19:24:53
%S 1,2,3,4,7,15,18,178
%N Record lengths of transient part of aliquot sequences corresponding to entries in A098009.
%e 138 has a transient of length 177 (see Guy's book).
%t g[n_] := If[n > 0, DivisorSigma[1, n] - n, 0]; f[n_] := NestWhileList[g, n, UnsameQ, All]; a = -1; Do[b = Length[ f[n]] - 2; If[b > a, a = b; Print[b]], {n, 275}] (* _Robert G. Wilson v_, Sep 10 2004 *)
%Y See A098009 for discussion, see A098008 for references.
%K nonn,more
%O 1,2
%A _N. J. A. Sloane_, Sep 10 2004
%E 18 and 178 from _Robert G. Wilson v_, Sep 10 2004
|
