OFFSET
0,2
COMMENTS
4.5*10^14 < a(4) <= 1081895452249824576000.
MATHEMATICA
nusigma[ n_ ] := DivisorSigma[ 1, n ]-Times@@(1+Power@@#&/@FactorInteger[ n ]); a[ n_ ] := For[ k=1, True, k++, If[ nusigma[ k ]==n k, Return[ k ] ] ]
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Dean Hickerson, Sep 25 2001
STATUS
approved