|
%I #21 Jan 12 2024 06:15:40
%S 42,54,66,78,90,144,259,45,33,15,9,4,3,1,0
%N Aliquot sequence starting at 42.
%C The sum-of-divisor function A000203 and aliquot parts A001065 are defined only for positive integers, so the trajectory ends when 0 is reached, here at index 14. - _M. F. Hasler_, Feb 24 2018
%D R. K. Guy, Unsolved Problems in Number Theory, B6.
%H Christophe CLAVIER, <a href="http://christophe.clavier.free.fr/Aliquot/site/Aliquot.html">Aliquot Sequences</a>
%H <a href="/index/Al#ALIQUOT">Index entries for sequences related to aliquot parts</a>.
%F a(n+1) = A001065(a(n)). - _R. J. Mathar_, Oct 11 2017
%F a(n) = A008885(n+1). - _R. J. Mathar_, Jan 12 2024
%p f := proc(n) option remember; if n = 0 then 42; else sigma(f(n-1))-f(n-1); fi; end:
%t Join[NestList[DivisorSigma[1,#]-#&,42,14],PadRight[{},60,0]] (* _Harvey P. Dale_, Apr 28 2014 *)
%o (PARI) a(n,a=42)={for(i=1,n,a=sigma(a)-a);a} \\ _M. F. Hasler_, Feb 24 2018
%Y Cf. A098007 (length of aliquot sequences); some other examples: A008885 (starting at 30) .. A008892 (starting at 276), A014360 (starting at 552) .. A014365 (starting at 1134), see link to index for a more complete list.
%K nonn,fini,full
%O 0,1
%A _N. J. A. Sloane_.
%E Edited by _M. F. Hasler_, Feb 24 2018
|
