OFFSET
1,1
MATHEMATICA
nusigma[ n_ ] := DivisorSigma[ 1, n ]-Times@@(1+Power@@#&/@FactorInteger[ n ]); For[ n=1, True, n++, v=nusigma[ n ]; If[ v>0&&Mod[ v, n ]==0, Print[ n ] ] ]
CROSSREFS
KEYWORD
nonn,hard
AUTHOR
Dean Hickerson, Sep 25 2001
STATUS
approved