OFFSET
1,2
COMMENTS
a(8) > 3.4*10^10, if it exists. - Amiram Eldar, Apr 06 2024
EXAMPLE
abundance(84) = sigma(84) - 2(84) = 56 = abundance(224) = abundance(sigma(84)), so 84 is a term in the sequence.
MATHEMATICA
abund[n_] := DivisorSigma[1, n]-2n; l = {}; For[i = 1, i <= 10^6, i++, If[abund[i] == abund[DivisorSigma[1, i]], l = Append[l, i]]]; l
PROG
(PARI) is(k) = {my(s = sigma(k)); s - 2*k == sigma(s) - 2*s; } \\ Amiram Eldar, Apr 06 2024
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Joseph L. Pe, Mar 05 2008
EXTENSIONS
a(7) from Donovan Johnson, Apr 27 2008
STATUS
approved