login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A046315 Odd semiprimes: odd numbers divisible by exactly 2 primes (counted with multiplicity). 90

%I

%S 9,15,21,25,33,35,39,49,51,55,57,65,69,77,85,87,91,93,95,111,115,119,

%T 121,123,129,133,141,143,145,155,159,161,169,177,183,185,187,201,203,

%U 205,209,213,215,217,219,221,235,237,247,249,253,259,265,267,287,289

%N Odd semiprimes: odd numbers divisible by exactly 2 primes (counted with multiplicity).

%C In general, the prime factors, p, of a(n) are given by: p = sqrt(a(n) + (k/2)^2) +- (k/2) where k is the positive difference of the prime factors. Equivalently, p = (1/2)( sqrt(4a(n) + k^2) +- k ). - _Wesley Ivan Hurt_, Jun 28 2013

%H Zak Seidov and K. D. Bajpai, <a href="/A046315/b046315.txt">Table of n, a(n) for n = 1..10000</a> (first 1956 terms from Zak Seidov)

%F Sum_{n>=1} 1/a(n)^s = (1/2)*(P(s)^2 + P(2*s)) - P(s)/2^s, for s>1, where P is the prime zeta function. - _Amiram Eldar_, Nov 21 2020

%e From _K. D. Bajpai_, Jul 05 2014: (Start)

%e 15 is a term because it is odd number and also 15 = 3 * 5, which is semiprime.

%e 39 is a term because it is odd number and also 39 = 3 * 13, which is semiprime. (End)

%p A046315 := proc(n) option remember; local r;

%p if n = 0 then RETURN(0) fi;

%p for r from A046315(n - 1) + 1 do

%p if r mod 2 = 1 and numtheory[bigomega](r) = 2

%p then RETURN(r) fi

%p od end:

%p seq(A046315(n),n=1..56); # _Peter Luschny_, Feb 15 2011

%t Reap[Do[If[Total[FactorInteger[n]][[2]] == 2, Sow[n]], {n, 1, 400, 2}]][[2,1]] (* _Zak Seidov_ *)

%t fQ[n_] := Plus @@ Last /@ FactorInteger@ n == 2; Select[2 Range@ 150 - 1, fQ] (* _Robert G. Wilson v_, Feb 15 2011 *)

%t Select[Range[5,301,2],PrimeOmega[#]==2&] (* _Harvey P. Dale_, May 22 2015 *)

%o (PARI) list(lim)=my(u=primes(primepi(lim\3)),v=List(),t); for(i=2,#u, for(j=i,#u, t=u[i]*u[j];if(t>lim,break); listput(v,t))); vecsort(Vec(v)) \\ _Charles R Greathouse IV_, Jul 19 2011

%o (Haskell)

%o a046315 n = a046315_list !! (n-1)

%o a046315_list = filter odd a001358_list -- _Reinhard Zumkeller_, Jan 02 2014

%Y Odd members of A001358.

%Y A046388 is a subsequence.

%Y Cf. A085770 (number of odd semiprimes < 10^n). - _Robert G. Wilson v_, Aug 25 2011

%K nonn

%O 1,1

%A _Patrick De Geest_, Jun 15 1998

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 16 17:01 EDT 2021. Contains 343050 sequences. (Running on oeis4.)