login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A072000 Number of semiprimes (A001358) <= n. 31

%I #51 Jul 01 2020 14:32:56

%S 0,0,0,1,1,2,2,2,3,4,4,4,4,5,6,6,6,6,6,6,7,8,8,8,9,10,10,10,10,10,10,

%T 10,11,12,13,13,13,14,15,15,15,15,15,15,15,16,16,16,17,17,18,18,18,18,

%U 19,19,20,21,21,21,21,22,22,22,23,23,23,23,24,24,24,24,24,25,25,25,26

%N Number of semiprimes (A001358) <= n.

%C Number of k <= n such that bigomega(k) = 2.

%D A. Hildebrand, On the number of prime factors of an integer, in Ramanujan Revisited (Urbana-Champaign, Ill., 1987), pp. 167-185, Academic Press, Boston, MA, 1988.

%D E. Landau, Handbuch der Lehre von der Verteilung der Primzahlen, vol. 1, Teubner, Leipzig, 1909; third edition : Chelsea, New York (1974).

%D G. Tenenbaum, Introduction à la théorie analytique et probabiliste des nombres, p. 203, Publications de l'Institut Cartan, 1990.

%H Daniel Forgues, <a href="/A072000/b072000.txt">Table of n, a(n) for n = 1..40882</a>

%H Dragos Crisan and Radek Erban, <a href="https://arxiv.org/abs/2006.16491">On the counting function of semiprimes</a>, arXiv:2006.16491 [math.NT], 2020.

%H E. Landau, Handbuch der Lehre von der Verteilung der Primzahlen, <a href="http://name.umdl.umich.edu/ABV2766.0001.001">vol. 1</a> and <a href="http://name.umdl.umich.edu/ABV2766.0002.001">vol. 2</a>, Leipzig, Berlin, B. G. Teubner, 1909.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Semiprime.html">Semiprime</a>

%F Let pi(x) denote the number of primes <= x (cf. A000720). Then 2*a(n) = Sum_{ primes p <= n/2 } Pi(n/p) + Pi(sqrt(n)). [Landau, p. 211]

%F Let pi(x) denote the number of primes <= x (cf. A000720). Then a(n) = Sum_{i=1..Pi(sqrt(n))} (Pi(n/prime(i)) - i + 1). - _Robert G. Wilson v_, Feb 07 2006

%F a(n) = card{ x <= n : bigomega(x) = 2 }.

%F Asymptotically a(n) ~ n*log(log(n))/log(n). [Landau, p. 211]

%F Let A be a positive integer. Then card{ x <= n : bigomega(x) = A } ~ (n/log(n))*log(log(n))^(A-1)/(A-1)! [Landau, p. 211]

%F a(n) = A072613(n)+A056811(n). - _R. J. Mathar_, Jun 10 2007

%F a(n) = sum_{i=1..n} A064911(i). - _Jonathan Vos Post_, Dec 30 2007

%F a(n)*A064911(n) = A174956(n). - _Reinhard Zumkeller_, Apr 03 2010

%p A072000 := proc(n) local sp,t ; sp := 0 ; for t from 1 to n do if numtheory[bigomega](t) = 2 then sp := sp+1 ; fi ; od ; sp ; end proc: # _R. J. Mathar_, Jun 10 2007

%t semiPrimePi[n_] := Sum[ PrimePi[n/Prime@i] -i + 1, {i, PrimePi@Sqrt@n}]; Array[semiPrimePi, 78] (* _Robert G. Wilson v_, Jan 03 2006 *)

%t (* If version >= 7 *) a[n_] := Select[Range[n], PrimeOmega[#] == 2 &] // Length; Table[a[n], {n, 1, 77}] (* _Jean-François Alcover_, Jun 29 2013 *)

%t Accumulate[Table[If[PrimeOmega[n]==2,1,0],{n,100}]] (* _Harvey P. Dale_, Jun 14 2014 *)

%o (PARI) for(n=1,100,print1(sum(i=1,n,if(bigomega(i)-2,0,1)),","))

%o (PARI) a(n)=my(s=0,i=0); forprime(p=2, sqrt(n), s+=primepi(n\p); i++); s - i * (i-1)/2 \\ _Charles R Greathouse IV_, Apr 21 2011

%Y Cf. A000720, A001358, A066265, A064911.

%K easy,nonn

%O 1,6

%A _Benoit Cloitre_, Jun 19 2002

%E Edited by _Robert G. Wilson v_, Feb 15 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:28 EDT 2024. Contains 371264 sequences. (Running on oeis4.)