OFFSET
1,1
COMMENTS
EXAMPLE
25 = 5^2 is a term because sigma(sigma(sigma(5^2))) = sigma(2^5) = sigma(sigma(5^2)) + sigma(5^2).
MATHEMATICA
f[n_] := DivisorSigma[1, n]; fQ[n_] := f[f[f[n]]] == f[f[n]] + f[n]; Select[ Range@1000000, fQ] (* Robert G. Wilson v, Sep 23 2017 *)
PROG
(PARI) s(n) = sigma(n);
isok(n) = s(s(s(n)))==s(s(n))+s(n);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Sep 15 2017
EXTENSIONS
a(9) from Giovanni Resta, Sep 15 2017
STATUS
approved