OFFSET
1,1
COMMENTS
This sequence is complete only as far as the last term given, for the eventual fate of the aliquot sequence generated by 276 is not (yet) known
REFERENCES
Benito, Manuel; Creyaufmueller, Wolfgang; Varona, Juan Luis; and Zimmermann, Paul; Aliquot Sequence 3630 Ends After Reaching 100 Digits; Experimental Mathematics, Vol. 11, No. 2, Natick, MA, 2002, pp. 201-206.
LINKS
Manuel Benito and Juan L. Varona, Advances In Aliquot Sequences, Mathematics of Computation, Vol. 68, No. 225, (1999), pp. 389-393.
Wolfgang Creyaufmueller, Aliquot sequences.
FORMULA
Define s(i)=sigma(i)-i=A000203(i)-i. Then if the aliquot sequence obtained by repeatedly applying the mapping i->s(i) terminates by encountering the prime 3 as a member of its trajectory, i is included in this sequence
EXAMPLE
a(5)=15 because the fifth integer whose aliquot sequence terminates by encountering the prime 3 as a member of its trajectory is 15. The complete aliquot sequence generated by iterating the proper divisors of 15 is 15->9->4->3->1->0
MATHEMATICA
s[n_] := DivisorSigma[1, n] - n; g[n_] := If[n > 0, s[n], 0]; Trajectory[n_] := Most[NestWhileList[g, n, UnsameQ, All]]; Select[Range[275], MemberQ[Trajectory[ # ], 3] &]
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
Ant King, Jan 07 2007
STATUS
approved