OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..400
PROG
(PARI) u(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
for(n=1, 10^9, if(sigma(n)-u(n)>3*n, print(n)))
(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) > 3*m, write("b063875.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Sep 01 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Aug 27 2001
EXTENSIONS
a(20)-a(25) from Harry J. Smith, Sep 01 2009
STATUS
approved