OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)
MATHEMATICA
nusigma[ n_ ] := DivisorSigma[ 1, n ]-Times@@(1+Power@@#&/@FactorInteger[ n ]); For[ n=1, True, n++, If[ nusigma[ n ]<n, Print[ n ] ] ]
PROG
(PARI) usigma(n)= { local(f, s=1); f=factor(n); for(i=1, matsize(f)[1], s*=1 + f[i, 1]^f[i, 2]); return(s) } { n=0; for (m=1, 10^9, if (sigma(m) - usigma(m) < m, write("b064598.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Sep 19 2009
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Dean Hickerson, Sep 25 2001
STATUS
approved