OFFSET
1,2
PROG
(PARI) s(n)=sigma(n)-n; for(n=1, 10^8, x=sum(k=1, n, s(k)); if(Mod(x, n)==0, print(n)))
(PARI) x = 0; for(n=1, 10^10, x = x + sigma(n)-n; if(x%n == 0, print1(n, ", "))) \\ Ryan Propper, Jan 13 2008
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Aug 29 2001
EXTENSIONS
More terms from Michel ten Voorde Jun 20 2003
More terms from Ryan Propper, Jan 13 2008
a(28)-a(30) from Donovan Johnson, Dec 29 2008
a(31)-a(32) from Donovan Johnson, Aug 07 2011
Name clarified by Jon E. Schoenfield, Nov 17 2018
STATUS
approved