OFFSET
1,1
COMMENTS
The terms are ordered according to their smaller counterparts (A348602).
MATHEMATICA
esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s[n_] := DivisorSigma[1, n] - esigma[n]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, m]], {n, 1, 1.7*10^6}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 25 2021
STATUS
approved