OFFSET
1,1
MATHEMATICA
a = 0; b = 0; Do[a = b; b = Apply[ Plus, Transpose[ FactorInteger[n + 1]] [[1]]]; If[ b - a == DivisorSigma[1, n] - n, Print[n]], {n, 1, 2*10^7} ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Sep 05 2001
EXTENSIONS
More terms from Robert G. Wilson v, Oct 05 2001
a(14)-a(27) from Donovan Johnson, Oct 13 2009
STATUS
approved