OFFSET
1,1
EXAMPLE
792736 + sigma(792736) + sigma(sigma(792736)) = 56^4.
PROG
(PARI) for(n=1, 10^9, s=sigma(n); if(ispower(n+s+sigma(s), 4), print(n))) /* Donovan Johnson, Feb 17 2013 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Resta, Feb 13 2006
EXTENSIONS
a(5)-a(26) from Donovan Johnson, Feb 17 2013
STATUS
approved