OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 1000 terms from Harry J. Smith)
PROG
(PARI) u(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
j=[]; for(n=1, 20000, if(sigma(n)-u(n)>2*n, j=concat(j, n))); j
(PARI) u(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d))
{ n=0; for (m=1, 10^9, if(sigma(m) - u(m) > 2*m, write("b063846.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Sep 01 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Aug 25 2001
STATUS
approved