OFFSET
1,1
COMMENTS
Quote from the abstract of the article by te Riele: "In this note, the existence of an aliquot sequence with more than 5092 monotonically increasing even terms is proved". The author uses the perfect number corresponding to the Mersenne prime 2^p-1 with p=19937 (whereas the script below only uses p=521).
LINKS
H. J. J. te Riele, A note on the Catalan-Dickson conjecture, Math. Comp. 27 (1973), 189-192.
PROG
(PARI) lista(p=521, nb) = {perf = 2^(p-1)*(2^p-1); a = 27*perf; print1(a/perf, ", "); for (i=1, nb, a = sigma(a) - a; print1(a/perf, ", "); if (gcd(a/perf, p) != 1, return()); ); } \\ Michel Marcus, Mar 13 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Mar 13 2013
STATUS
approved