|
%I #16 Dec 15 2017 17:35:24
%S 1,4,9,16,25,49,64,81,121,169,225,256,289,361,441,484,529,625,676,729,
%T 841,961,1024,1089,1156,1225,1369,1444,1521,1681,1849,2025,2116,2209,
%U 2401,2601,2809,3025,3249,3364,3481,3721,3844,3969,4096,4225,4489,4761
%N Square deficient numbers.
%H Harry J. Smith, <a href="/A063735/b063735.txt">Table of n, a(n) for n=1..1000</a>
%t Select[#^2 & /@ Range[70], DivisorSigma[1, #] < 2 # &]
%o (PARI) j=[]; for(n=1,8000, if(sigma(n)<(2*n),a=n; if(issquare(a),j=concat(j,a)))); j
%o (PARI) { n=0; for (m=1, 10^9, if (issquare(m) && sigma(m)<(2*m), write("b063735.txt", n++, " ", m); if (n==1000, break)) ) } \\ _Harry J. Smith_, Aug 29 2009
%o (PARI) { n=0; for (m=1, 10^9, s=m^2; if (sigma(s)<(2*s), write("b063735.txt", n++, " ", s); if (n==1000, break)) ) } \\ _Harry J. Smith_, Aug 29 2009]
%Y Intersection of A000290 and A005100.
%K nonn
%O 1,2
%A _Jason Earls_, Aug 13 2001
|
