This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A245214 #20 Jul 19 2024 14:27:56 %S A245214 144,192,216,240,288,336,360,384,432,480,504,540,576,600,648,672,720, %T A245214 768,792,840,864,900,936,960,1008,1056,1080,1152,1200,1248,1260,1296, %U A245214 1320,1344,1440,1512,1536,1560,1584,1620,1632,1680,1728,1800,1824,1848,1872,1920,1944,1980,2016,2040,2100,2112,2160,2240 %N A245214 Numbers k such that A245212(k) < 0. %C A245214 If d are divisors of k then values of sequence A245212(k) are by bending moments in point 0 of static forces of sizes tau(d) operating in places d on the cantilever as the nonnegative number axis of length k with bracket in point 0 by the schema: A245212(k) = (k * tau(k)) - Sum_{(d<k) | k} (d * tau(d)). %C A245214 Numbers k such that A038040(k) = k * tau(k) < A245211(k) = Sum_{(d<k) | k} (d * tau(d)). %C A245214 From _Amiram Eldar_, Jul 19 2024: (Start) %C A245214 Numbers whose divisors have a mean abundancy index that is larger than 2. %C A245214 The numbers of terms that do not exceed 10^k, for k = 3, 4, ..., are 24, 243, 2571, 25583, 254794, 2551559, 25514104, 255112225, ... . Apparently, the asymptotic density of this sequence exists and equals 0.02551... . %C A245214 The least odd term in this sequence is a(276918705) = 10854718875. (End) %H A245214 Jens Kruse Andersen, <a href="/A245214/b245214.txt">Table of n, a(n) for n = 1..10000</a> %e A245214 Number 144 is in sequence because 144 * tau(144) = 2160 < Sum_{(d<144) | 144} (d * tau(d)) = 2226. %t A245214 f[p_, e_] := ((e+1)*p^2 - (e+2)*p + p^(-e))/((e+1)*(p-1)^2); s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[2500], s[#] > 2 &] (* _Amiram Eldar_, Jul 19 2024 *) %o A245214 (Magma) [n: n in [1..100000] | (2*(n*(#[d: d in Divisors(n)]))-(&+[d*#([e: e in Divisors(d)]): d in Divisors(n)])) lt 0] %o A245214 (PARI) isok(n) = (n*numdiv(n) - sumdiv(n, d, (d<n)*d*numdiv(d))) < 0; \\ _Michel Marcus_, Aug 06 2014 %o A245214 (PARI) is(n) = {my(f = factor(n)); prod(i = 1, #f~, p=f[i,1]; e=f[i,2]; (-2*p - e*p + p^2 + e*p^2 + p^(-e))/((e + 1)*(p - 1)^2)) > 2;} \\ _Amiram Eldar_, Jul 19 2024 %Y A245214 Cf. A000005, A038040, A245211, A245212, A245213. %K A245214 nonn %O A245214 1,1 %A A245214 _Jaroslav Krizek_, Jul 23 2014