OFFSET
1,1
COMMENTS
EXAMPLE
(23, 37) are the first pair since sigma(sigma(23)) = sigma(sigma(37)) = 60 = 23 + 37.
MATHEMATICA
seq={}; s[n_]:=DivisorSigma[1, DivisorSigma[1, n]]-n; Do[m=s[n]; If[m>n && s[m]==n, AppendTo[seq, m]], {n, 1, 60000}]; seq
PROG
(PARI) f(n) = sigma(sigma(n)) - n;
lista(nn) = {for (n=1, nn, my(fn = f(n)); if ((fn > n) && (f(fn) == n), print1(fn, ", ")); ); } \\ Michel Marcus, Feb 20 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 19 2019
STATUS
approved