|
%I #23 Feb 28 2021 20:28:06
%S 6,8,9,12,14,15,16,18,20,21,24,27,28,30,32,35,36,40,42,44,45,48,49,50,
%T 52,54,56,60,63,64,66,68,70,72,75,76,78,80,81,84,88,90,92,96,98,99,
%U 100,102,104,105,108,110,112,114,117,120,124,125,126,128,130,132,135,136,138,140,144,147,148,150,152,153,154
%N Numbers k for which A003973(k) >= 2*sigma(k).
%C Note that A003973(n) >= sigma(n) for all n. See A336852.
%C Like the abundancy index (ratio A000203(n)/n), and ratio A003961(n)/n, the ratio A003973(n)/sigma(n) is also multiplicative and > 1 for all n > 1. Thus if the number has a proper divisor that is in this sequence, then the number itself is also. See A337543 for those terms included here, but which have no proper divisor in this sequence. - _Antti Karttunen_, Aug 31 2020
%C All terms are in A246282 because A341528(n) < A341529(n) for all n > 1. - _Antti Karttunen_, Feb 22 2021
%H Antti Karttunen, <a href="/A337381/b337381.txt">Table of n, a(n) for n = 1..20000</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%o (PARI)
%o A003973(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); sigma(factorback(f)); };
%o isA337381(n) = (A003973(n)>=2*sigma(n));
%Y Cf. A000203, A003961, A003973, A336852, A337541, A337542, A341528, A341529.
%Y Cf. A337382 (complement), A337383 (characteristic function).
%Y Subsequences: A337378, A337384, A337386, A337543 (primitive terms).
%Y Subsequence of A246282.
%K nonn
%O 1,1
%A _Antti Karttunen_, Aug 27 2020
|
