|
%I M5486 #68 Nov 22 2020 20:24:08
%S 945,1575,2205,3465,4095,5355,5775,5985,6435,6825,7245,7425,8085,8415,
%T 8925,9135,9555,9765,11655,12705,12915,13545,14805,15015,16695,18585,
%U 19215,19635,21105,21945,22365,22995,23205,24885,25935,26145,26565,28035,28215
%N Odd primitive abundant numbers.
%C Dickson proves that there are only a finite number of odd primitive abundant numbers having n distinct prime factors. Sequence A188342 lists the smallest such numbers. - _T. D. Noe_, Mar 28 2011
%C Sequence A188439 sorts the numbers in this sequence by the number of distinct prime factors. Eight numbers have exactly three prime factors; 576 have exactly four prime factors. - _T. D. Noe_, Apr 04 2011
%C Any multiple of an abundant number (A005101) is again an abundant number. Primitive abundant numbers (A091191) are those not of this form, i.e., without an abundant proper divisor. We don't know any odd perfect number (A000396), so the (odd) terms here have only deficient proper divisors (A071395), and their prime factors p are less than sigma(n/p)/deficiency(n/p). See A005231 (odd abundant numbers) for an explanation why all terms have at least 3 distinct prime factors, and 5 prime factors when counted with multiplicity (A001222), whence a(1) = 3^3*5*7. All known terms are semiperfect (A005835, and thus in A006036): no odd weird number (A006037) is known, but if it exists, the smallest one is in this sequence. - _M. F. Hasler_, Jul 28 2016
%C So far, a(173) = 351351 is the only known term of A122036, i.e., which can't be written as sum of its proper divisors > 1. - _M. F. Hasler_, Jan 26 2020
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H T. D. Noe, <a href="/A006038/b006038.txt">Table of n, a(n) for n = 1..10000</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 R. K. Guy, <a href="/A001599/a001599_1.pdf">Letter to N. J. A. Sloane with attachment, Jun. 1991</a>
%H Jacob Liddy, <a href="https://ideaexchange.uakron.edu/honors_research_projects/728">An algorithm to determine all odd primitive abundant numbers with d prime divisors</a>, Honors Research Projects (2018), 728.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimitiveAbundantNumber.html">Primitive Abundant Number</a>
%p isA005101 := proc(n) is(numtheory[sigma](n) > 2*n ); end proc:
%p isA005100 := proc(n) is(numtheory[sigma](n) < 2*n ); end proc:
%p isA006038 := proc(n) local d; if type(n,'odd') and isA005101(n) then for d in numtheory[divisors](n) minus {1,n} do if not isA005100(d) then return false; end if; end do: return true;else false; end if; end proc:
%p n := 1 ; for a from 1 by 2 do if isA006038(a) then printf("%d %d\n",n,a) ; n := n+1 ; end if; end do: # _R. J. Mathar_, Mar 28 2011
%t t = {}; n = 1; While[Length[t] < 50, n = n + 2; If[DivisorSigma[1, n] > 2 n && Intersection[t, Divisors[n]] == {}, AppendTo[t, n]]]; t (* _T. D. Noe_, Mar 28 2011 *)
%o (PARI) is(n)=n%2 && sumdiv(n,d,sigma(d,-1)>2)==1 \\ _Charles R Greathouse IV_, Jun 10 2013
%o (PARI) is_A006038(n)=bittest(n,0) && sigma(n)>2*n && !for(i=1,#f=factor(n)[,1],sigma(n\f[i],-1)>2&&return) \\ More than 5 times faster. - _M. F. Hasler_, Jul 28 2016
%o (Haskell)
%o a006038 n = a006038_list !! (n-1)
%o a006038_list = filter f [1, 3 ..] where
%o f x = sum pdivs > x && all (<= 0) (map (\d -> a000203 d - 2 * d) pdivs)
%o where pdivs = a027751_row x
%o -- _Reinhard Zumkeller_, Jan 31 2014
%Y Cf. A005101, A005231. Subsequence of A091191.
%Y Cf. A000203, A027751.
%K nonn
%O 1,1
%A _N. J. A. Sloane_
|
