|
%I #15 Aug 13 2019 11:22:47
%S 1,840,900,1080,1225,1260,1440,1600,1680,1800,1848,1890,1980,2016,
%T 2100,2160,2340,2400,2520,2640,2700,2772,2800,2880,2970,3024,3080,
%U 3120,3136,3150,3240,3276,3300,3360,3465,3528,3600,3640,3696,3780,3900,3960,3969
%N Numbers n such that s(n) >= sigma(n), where s(n) = A095118(n) is the sum of the squares of the divisors of n which are <= sqrt(n) and sigma(n) = A000203(n) is the sum of the divisors of n.
%H Amiram Eldar, <a href="/A095119/b095119.txt">Table of n, a(n) for n = 1..10000</a>
%e 840 is in the sequence because s(840) = 3070 >= 2880 = sigma(840).
%t s[n_]:=Plus@@(Select[Divisors[n], #^2<=n&]^2); Select[Range[4000], s[ # ]>=DivisorSigma[1, # ]&]
%o (PARI) isok(n) = sumdiv(n, d, if (d^2 <= n, d^2)) >= sigma(n); \\ _Michel Marcus_, Aug 13 2019
%Y Cf. A000203, A095118.
%K nonn
%O 1,2
%A _Dean Hickerson_, following a suggestion of _Leroy Quet_, May 28 2004
|
