OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Range[600], DivisorSigma[1, #]^2>4*DivisorSigma[2, #]&] (* Harvey P. Dale, Dec 27 2015 *)
PROG
(PARI) for(n=1, 1000, if(sigma(n)^2>4*sumdiv(n, k, k^2), print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 04 2002
STATUS
approved