OFFSET
1,2
EXAMPLE
MATHEMATICA
esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s[1] = 1; s[n_] := DivisorSigma[1, n] - esigma[n]; Select[Range[500000], s[#] == s[# + 1] &]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Oct 26 2021
STATUS
approved