|
%I #110 Nov 18 2023 08:36:34
%S 945,1575,2205,2835,3465,4095,4725,5355,5775,5985,6435,6615,6825,7245,
%T 7425,7875,8085,8415,8505,8925,9135,9555,9765,10395,11025,11655,12285,
%U 12705,12915,13545,14175,14805,15015,15435,16065,16695,17325,17955
%N Odd abundant numbers (odd numbers m whose sum of divisors exceeds 2m).
%C While the first even abundant number is 12 = 2^2*3, the first odd abundant is 945 = 3^3*5*7, the 232nd abundant number.
%C Schiffman notes that 945+630k is in this sequence for all k < 52. Most of the first initial terms are of the form. Among the 1996 terms below 10^6, 1164 terms are of that form, and only 26 terms are not divisible by 5, cf. A064001. - _M. F. Hasler_, Jul 16 2016
%C From _M. F. Hasler_, Jul 28 2016: (Start)
%C Any multiple of an abundant number is again abundant, see A006038 for primitive terms, i.e., those which are not a multiple of an earlier term.
%C An odd abundant number must have at least 3 distinct prime factors, and 5 prime factors when counted with multiplicity (A001222), whence a(1) = 3^3*5*7. To see this, write the relative abundancy A(N) = sigma(N)/N = sigma[-1](N) as A(Product p_i^e_i) = Product (p_i-1/p_i^e_i)/(p_i-1) < Product p_i/(p_i-1).
%C See A115414 for terms not divisible by 3, A064001 for terms not divisible by 5, A112640 for terms coprime to 5*7, and A047802 for other generalizations.
%C As of today, we don't know of a difference between this set S of odd abundant numbers and the set S' of odd semiperfect numbers: Elements of S' \ S would be perfect (A000396), and elements of S \ S' would be weird (A006037), but no odd weird or perfect number is known. (End)
%C For any term m in this sequence, A064989(m) is also an abundant number (in A005101), and for any term x in A115414, A064989(x) is in this sequence. If there are no odd perfect numbers, then applying A064989 to these terms and sorting into ascending order gives A337386. - _Antti Karttunen_, Aug 28 2020
%C There exist infinitely many terms m such that 2*m+1 is also a term. An example of such a term is given by m = 985571808130707987847768908867571007187. - _Max Alekseyev_, Nov 16 2023
%D W. Dunham, Euler: The Master of Us All, The Mathematical Association of America Inc., Washington, D.C., 1999, p. 13.
%D R. K. Guy, Unsolved Problems in Number Theory, B2.
%H Metin Sariyar, <a href="/A005231/b005231.txt">Table of n, a(n) for n = 1..32000</a> (terms 1..1000 from T. D. Noe)
%H Jill Britton, <a href="http://britton.disted.camosun.bc.ca/perfect/jbperfect.htm">Perfect Number Analyzer</a>.
%H L. E. Dickson, <a href="http://www.jstor.org/stable/2370405">Finiteness of the odd perfect and primitive abundant numbers with n distinct prime factors</a>, American Journal of Mathematics 35 (1913), pp. 413-422.
%H Mitsuo Kobayashi, Paul Pollack and Carl Pomerance, <a href="https://doi.org/10.1016/j.jnt.2008.10.011">On the distribution of sociable numbers</a>, Journal of Number Theory, Vol. 129, No. 8 (2009), pp. 1990-2009. See Theorem 10 on p. 2007.
%H Victor Meally, <a href="/A006556/a006556.pdf">Letter to N. J. A. Sloane</a>, no date.
%H Walter Nissen, <a href="http://upforthecount.com/math/abundance.html">Abundancy : Some Resources</a>.
%H Jay L. Schiffman, <a href="http://www.appliedprobability.org/data/files/MS%20issues/Vol37_No2.pdf">Odd Abundant Numbers</a>, Mathematical Spectrum, Volume 37, Number 2 (January 2005), pp 73-75.
%H Jay L. Schiffman and Christopher S. Simons, <a href="http://users.rowan.edu/%7Esimons/abundant.pdf">More Odd Abundant Sequences</a>, Volume 38, Number 1 (September 2005), pp. 7-8.
%F a(n) ~ k*n for some constant k (perhaps around 500). - _Charles R Greathouse IV_, Apr 21 2022
%F 482.8 < k < 489.8 (based on density bounds by Kobayashi et al., 2009). - _Amiram Eldar_, Jul 17 2022
%p A005231 := proc(n) option remember ; local a ; if n = 1 then 945 ; else for a from procname(n-1)+2 by 2 do if numtheory[sigma](a) > 2*a then return a; end if; end do: end if; end proc: # _R. J. Mathar_, Mar 20 2011
%t fQ[n_] := DivisorSigma[1, n] > 2n; Select[1 + 2Range@ 9000, fQ] (* _Robert G. Wilson v_, Mar 20 2011 *)
%o (PARI) je=[]; forstep(n=1,15000,2, if(sigma(n)>2*n, je=concat(je,n))); je
%o (PARI) is_A005231(n)={bittest(n,0)&&sigma(n)>2*n} \\ _M. F. Hasler_, Jul 28 2016
%o (PARI) list(lim)=my(v=List()); forfactored(n=945,lim\1, if(n[2][1,1]>2 && sigma(n,-1)>2, listput(v,n[1]))); Vec(v) \\ _Charles R Greathouse IV_, Apr 21 2022
%Y Cf. A000203, A005835, A006038, A115414, A064001, A112640, A122036, A136446, A005101, A173490, A039725, A064989, A337386.
%K nonn
%O 1,1
%A _N. J. A. Sloane_
%E More terms from _James A. Sellers_
|
