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!)
A106350 Semiprimes indexed by primes. 16
6, 9, 14, 21, 33, 35, 49, 55, 65, 86, 91, 115, 122, 129, 142, 159, 183, 187, 206, 215, 218, 247, 259, 287, 303, 319, 323, 334, 339, 358, 403, 415, 446, 451, 482, 489, 511, 527, 537, 553, 573, 581, 626, 633, 655, 667, 698, 737, 753, 758, 771, 791, 794, 835, 851 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

This is the sequence of the n-th semiprime for n = {2,3,5,7,11,13,17,19,23,29...}. Not to be confused with A106349: Primes indexed by semiprimes. We seek to know what this sequence is asymptotically, as J. B. Rosser's result, subsequently modified, is that prime(n) ~ n*(log n + log log n - 1). hence semiprime(prime(n)) ~ semiprime(n)*(log semiprime(n) + log log semiprime(n) - 1). But what is, asymptotically, semiprime(n)?

Semiprime(n) ~ n log n / log log n, hence a(n) ~ n log^2 n / log log n. - Charles R Greathouse IV, Dec 28 2011

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..10000

J. B. Rosser, The n-th Prime is Greater than n log(n), Proc. London Math. Soc. 45, 21-44, 1939.

Eric Weisstein's World of Mathematics, Semiprime.

FORMULA

a(n) = semiprime(prime(n)). a(n) = A001358(A000040(n)).

a(n) ~ n log^2 n / log log n. - Charles R Greathouse IV, Dec 28 2011

EXAMPLE

a(1) = semiprime(prime(1)) = semiprime(2) = 6.

a(2) = semiprime(prime(2)) = semiprime(3) = 9.

MAPLE

A001358 := proc(n) if n = 1 then 4; else for a from procname(n-1)+1 do if numtheory[bigomega](a) = 2 then return a ; end if; end do ; end if ; end proc: A106350 := proc(n) A001358(ithprime(n)) ; end proc: seq(A106350(n), n=1..80) ; # R. J. Mathar, Dec 14 2009

MATHEMATICA

terms = 55;

semiPrimes = Select[Range[16 terms], PrimeOmega[#] == 2&];

(* NB If the index Prime[terms] exceeds the size of the table semiPrimes, then the coefficient 16 has to be increased according to the number of terms desired: for instance, for 1000 terms, replace 16 with 32. *)

a[n_] := semiPrimes[[Prime[n]]];

Array[a, terms] (* Jean-Fran├žois Alcover, Apr 13 2020 *)

CROSSREFS

Cf. A000040, A001358, A007097, A091022, A105997, A105998, A106349.

Sequence in context: A316013 A154778 A316014 * A217851 A218697 A316015

Adjacent sequences:  A106347 A106348 A106349 * A106351 A106352 A106353

KEYWORD

easy,nonn

AUTHOR

Jonathan Vos Post, Apr 30 2005

EXTENSIONS

All values after a(32) corrected by R. J. Mathar, Dec 14 2009

STATUS

approved

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 June 1 04:57 EDT 2020. Contains 334758 sequences. (Running on oeis4.)