A062198 -id:A062198 - cp's OEIS frontend

A001358 Semiprimes (or biprimes): products of two primes.

4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, 34, 35, 38, 39, 46, 49, 51, 55, 57, 58, 62, 65, 69, 74, 77, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 121, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 169, 177, 178, 183, 185, 187

Offset: 1

N. J. A. Sloane, R. K. Guy

Numbers of the form p*q where p and q are primes, not necessarily distinct.
These numbers are sometimes called semiprimes or 2-almost primes.
Numbers n such that Omega(n) = 2 where Omega(n) = A001222(n) is the sum of the exponents in the prime decomposition of n.
Complement of A100959; A064911(a(n)) = 1. - Reinhard Zumkeller, Nov 22 2004
The graph of this sequence appears to be a straight line with slope 4. However, the asymptotic formula shows that the linearity is an illusion and in fact a(n)/n ~ log(n)/log(log(n)) goes to infinity. See also the graph of A066265 = number of semiprimes < 10^n.
For numbers between 33 and 15495, semiprimes are more plentiful than any other k-almost prime. See A125149.
Numbers that are divisible by exactly 2 prime powers (not including 1). - Jason Kimberley, Oct 02 2011
The (disjoint) union of A006881 and A001248. - Jason Kimberley, Nov 11 2015
An equivalent definition of this sequence is a'(n) = smallest composite number which is not divided by any smaller composite number a'(1),...,a'(n-1). - Meir-Simchah Panzer, Jun 22 2016
The above characterization can be simplified to "Composite numbers not divisible by a smaller term." This shows that this is the equivalent of primes computed via Eratosthenes's sieve, but starting with the set of composite numbers (i.e., complement of 1 union primes) instead of all positive integers > 1. It's easy to see that iterating the method (using Eratosthenes's sieve each time on the remaining numbers, complement of the previously computed set) yields numbers with bigomega = k for k = 0, 1, 2, 3, ..., i.e., {1}, A000040, this, A014612, etc. - M. F. Hasler, Apr 24 2019
For all n except n = 2, a(n) is a deficient number. - Amrit Awasthi, Sep 10 2024
It is reasonable to assume that the "comforting numbers" which John T. Williams found in Chapter 3 of Milne's book "The House at Pooh Corner" are these semiprimes. Winnie-the-Pooh wonders whether he has 14 or 15 honey pots and concludes: "It's sort of comforting." To arrange a semiprime number of honey pots in a rectangular way, let's say on a shelf, with the larger divisor parallel to the wall, there is only one solution and this is for a simple mind like Winnie-the-Pooh comforting. - Ruediger Jehn, Dec 12 2024

			From _Gus Wiseman_, May 27 2021: (Start)
The sequence of terms together with their prime factors begins:
   4 = 2*2     46 = 2*23     91 = 7*13    141 = 3*47
   6 = 2*3     49 = 7*7      93 = 3*31    142 = 2*71
   9 = 3*3     51 = 3*17     94 = 2*47    143 = 11*13
  10 = 2*5     55 = 5*11     95 = 5*19    145 = 5*29
  14 = 2*7     57 = 3*19    106 = 2*53    146 = 2*73
  15 = 3*5     58 = 2*29    111 = 3*37    155 = 5*31
  21 = 3*7     62 = 2*31    115 = 5*23    158 = 2*79
  22 = 2*11    65 = 5*13    118 = 2*59    159 = 3*53
  25 = 5*5     69 = 3*23    119 = 7*17    161 = 7*23
  26 = 2*13    74 = 2*37    121 = 11*11   166 = 2*83
  33 = 3*11    77 = 7*11    122 = 2*61    169 = 13*13
  34 = 2*17    82 = 2*41    123 = 3*41    177 = 3*59
  35 = 5*7     85 = 5*17    129 = 3*43    178 = 2*89
  38 = 2*19    86 = 2*43    133 = 7*19    183 = 3*61
  39 = 3*13    87 = 3*29    134 = 2*67    185 = 5*37
(End)

Archimedeans Problems Drive, Eureka, 17 (1954), 8.
Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; Chapter II, Problem 60.
Edmund Landau, Handbuch der Lehre von der Verteilung der Primzahlen, Vol. 1, Teubner, Leipzig; third edition: Chelsea, New York (1974). See p. 211.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
John T. Williams, Pooh and the Philosophers, Dutton Books, 1995.

N. J. A. Sloane, Table of n, a(n) for n = 1..20000 (first 10000 terms from T. D. Noe)
Dragos Crisan and Radek Erban, On the counting function of semiprimes, INTEGERS, Vol. 21 (2021), #A122.
Daniel A. Goldston, Sidney W. Graham, János Pintz and Cem Y. Yildirim, Small gaps between primes or almost primes, Transactions of the American Mathematical Society, Vol. 361, No. 10 (2009), pp. 5285-5330, arXiv preprint, arXiv:math/0506067 [math.NT], 2005.
Richard K. Guy, Letters to N. J. A. Sloane, June-August 1968
Sh. T. Ishmukhametov and F. F. Sharifullina, On distribution of semiprime numbers, Izvestiya Vysshikh Uchebnykh Zavedenii. Matematika, 2014, No. 8, pp. 53-59. English translation, Russian Mathematics, Vol. 58, No. 8 (2014), pp. 43-48, alternative link.
Donovan Johnson, Jonathan Vos Post, and Robert G. Wilson v, Selected n and a(n). (2.5 MB)
Dixon Jones, Quickie 593, Mathematics Magazine, Vol. 47, No. 3, May 1974, p. 167.
Edmund Landau, Handbuch der Lehre von der Verteilung der Primzahlen, vol. 1 and vol. 2, Leipzig, Berlin, B. G. Teubner, 1909. See Vol. 1, p. 211.
Xianmeng Meng, On sums of three integers with a fixed number of prime factors, Journal of Number Theory, Vol. 114, No. 1 (2005), pp. 37-65.
Michael Penn, What makes a number "good"?, YouTube video, 2022.
Carlos Rivera, Conjecture 108: On the counting function of semiprimes, The Prime Puzzles & Problems Connection.
Eric Weisstein's World of Mathematics, Semiprime.
Eric Weisstein's World of Mathematics, Almost Prime.
Wikipedia, Almost prime.
Robert G. Wilson v, Subsequences at various powers of 10.
Index to sequences related to sums of cubes
Index entries for "core" sequences

Cf. A064911 (characteristic function).
Cf. A048623, A048639, A000040 (primes), A014612 (products of 3 primes), A014613, A014614, A072000 ("pi" for semiprimes), A065516 (first differences).
Cf. A077554, A077555, A002024, A072966, A100592, A014673, A068318, A061299, A087718, A089994, A089995, A096916, A096932, A106550, A106554, A108541, A108542, A126663, A131284, A138510, A138511, A072931, A088183, A171963, A237040 (semiprimes of form n^3 + 1).
Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r=1), this sequence (r=2), A014612 (r=3), A014613 (r=4), A014614 (r=5), A046306 (r=6), A046308 (r=7), A046310 (r=8), A046312 (r=9), A046314 (r=10), A069272 (r=11), A069273 (r=12), A069274 (r=13), A069275 (r=14), A069276 (r=15), A069277 (r=16), A069278 (r=17), A069279 (r=18), A069280 (r=19), A069281 (r=20).
These are the Heinz numbers of length-2 partitions, counted by A004526.
The squarefree case is A006881 with odd/even terms A046388/A100484 (except 4).
Including primes gives A037143.
The odd/even terms are A046315/A100484.
Partial sums are A062198.
The prime factors are A084126/A084127.
Grouping by greater factor gives A087112.
The product/sum/difference of prime indices is A087794/A176504/A176506.
Positions of even/odd terms are A115392/A289182.
The terms with relatively prime/divisible prime indices are A300912/A318990.
Factorizations using these terms are counted by A320655.
The prime indices are A338898/A338912/A338913.
Grouping by weight (sum of prime indices) gives A338904, with row sums A024697.
The terms with even/odd weight are A338906/A338907.
The terms with odd/even prime indices are A338910/A338911.
The least/greatest term of weight n is A339114/A339115.
Cf. A014342, A098350, A112141, A320732, A332765, A339003/A339004, A339116.

a001358 n = a001358_list !! (n-1)
a001358_list = filter ((== 2) . a001222) [1..]

[n: n in [2..200] | &+[d[2]: d in Factorization(n)] eq 2]; // Bruno Berselli, Sep 09 2015

A001358 := proc(n) option remember; local a; 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:
seq(A001358(n), n=1..120) ; # R. J. Mathar, Aug 12 2010

Select[Range[200], Plus@@Last/@FactorInteger[#] == 2 &] (* Zak Seidov, Jun 14 2005 *)
Select[Range[200], PrimeOmega[#]==2&] (* Harvey P. Dale, Jul 17 2011 *)

select( isA001358(n)={bigomega(n)==2}, [1..199]) \\ M. F. Hasler, Apr 09 2008; added select() Apr 24 2019

list(lim)=my(v=List(),t);forprime(p=2, sqrt(lim), t=p;forprime(q=p, lim\t, listput(v,t*q))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Sep 11 2011

A1358=List(4); A001358(n)={while(#A1358M. F. Hasler, Apr 24 2019

from sympy import factorint
def ok(n): return sum(factorint(n).values()) == 2
print([k for k in range(1, 190) if ok(k)]) # Michael S. Branicky, Apr 30 2022

from math import isqrt
from sympy import primepi, prime
def A001358(n):
    def f(x): return int(n+x-sum(primepi(x//prime(k))-k+1 for k in range(1, primepi(isqrt(x))+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Jul 23 2024

a(n) ~ n*log(n)/log(log(n)) as n -> infinity [Landau, p. 211], [Ayoub].
Recurrence: a(1) = 4; for n > 1, a(n) = smallest composite number which is not a multiple of any of the previous terms. - Amarnath Murthy, Nov 10 2002
A174956(a(n)) = n. - Reinhard Zumkeller, Apr 03 2010
a(n) = A088707(n) - 1. - Reinhard Zumkeller, Feb 20 2012
Sum_{n>=1} 1/a(n)^s = (1/2)*(P(s)^2 + P(2*s)), where P is the prime zeta function. - Enrique Pérez Herrero, Jun 24 2012
sigma(a(n)) + phi(a(n)) - mu(a(n)) = 2*a(n) + 1. mu(a(n)) = ceiling(sqrt(a(n))) - floor(sqrt(a(n))). - Wesley Ivan Hurt, May 21 2013
mu(a(n)) = -Omega(a(n)) + omega(a(n)) + 1, where mu is the Moebius function (A008683), Omega is the count of prime factors with repetition, and omega is the count of distinct prime factors. - Alonso del Arte, May 09 2014
a(n) = A078840(2,n). - R. J. Mathar, Jan 30 2019
A100484 UNION A046315. - R. J. Mathar, Apr 19 2023
Conjecture: a(n)/n ~ (log(n)/log(log(n)))*(1-(M/log(log(n)))) as n -> oo, where M is the Mertens's constant (A077761). - Alain Rocchelli, Feb 02 2025

More terms from James Sellers, Aug 22 2000

A176506 Difference between the prime indices of the two factors of the n-th semiprime.

Original entry on oeis.org

0, 1, 0, 2, 3, 1, 2, 4, 0, 5, 3, 6, 1, 7, 4, 8, 0, 5, 2, 6, 9, 10, 3, 7, 11, 1, 12, 4, 13, 8, 2, 9, 14, 5, 15, 10, 6, 16, 3, 0, 17, 11, 12, 4, 18, 13, 19, 1, 7, 20, 8, 21, 14, 5, 22, 0, 15, 23, 16, 9, 2, 24, 17, 25, 6, 10, 26, 3, 18, 27, 11, 7, 28, 19, 1, 29, 12, 20, 2, 21, 4, 30, 8, 31, 13, 22

Offset: 1

Juri-Stepan Gerasimov, Apr 19 2010

Are there no adjacent equal terms? I have verified this up to n = 10^6. - Gus Wiseman, Dec 04 2020

			From _Gus Wiseman_, Dec 04 2020: (Start)
The sequence of semiprimes together with the corresponding differences begins:
   4: 1 - 1 = 0
   6: 2 - 1 = 1
   9: 2 - 2 = 0
  10: 3 - 1 = 2
  14: 4 - 1 = 3
  15: 3 - 2 = 1
  21: 4 - 2 = 2
  22: 5 - 1 = 4
  25: 3 - 3 = 0
  26: 6 - 1 = 5
  33: 5 - 2 = 3
(End)

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A109313.
Cf. A001358, A049084, A084126, A084127.
A087794 is product of the same indices.
A176504 is the sum of the same indices.
A115392 lists positions of first appearances.
A128301 lists positions of 0's.
A172348 lists positions of 1's.
A338898 has this sequence as row differences.
A338900 is the squarefree case.
A338912/A338913 give the two prime indices of semiprimes.
A006881 lists squarefree semiprimes.
A024697 is the sum of semiprimes of weight n.
A056239 gives sum of prime indices (Heinz weight).
A087112 groups semiprimes by greater factor.
A270650/A270652/A338899 give the prime indices of squarefree semiprimes.
A338904 groups semiprimes by weight.
A338907/A338906 list semiprimes of odd/even weight.
A339114/A339115 give the least/greatest semiprime of weight n.
Cf. A000040, A001222, A062198, A112141, A112798, A325699, A320655, A339005.

isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:
A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:
A084126 := proc(n) min(op(numtheory[factorset](A001358(n)))) ; end proc:
A084127 := proc(n) max(op(numtheory[factorset](A001358(n)))) ; end proc:
A176506 := proc(n) numtheory[pi](A084127(n)) - numtheory[pi](A084126(n)) ; end proc: seq(A176506(n),n=1..120) ; # R. J. Mathar, Apr 22 2010
# Alternative:
N:= 500: # to use the first N semiprimes
Primes:= select(isprime, [2,seq(i,i=3..N/2,2)]):
SP:= NULL:
for i from 1 to nops(Primes) do
  for j from 1 to i do
    sp:= Primes[i]*Primes[j];
    if sp > N then break fi;
    SP:= SP, [sp, i-j]
od od:
SP:= sort([SP],(s,t) -> s[1] t[2], SP); # Robert Israel, Jan 17 2019

M = 500; (* to use the first M semiprimes *)
primes = Select[Join[{2}, Range[3, M/2, 2]], PrimeQ];
SP = {};
For[i = 1, i <= Length[primes], i++,
  For[j = 1, j <= i, j++,
    sp = primes[[i]] primes[[j]];
    If[sp > M, Break []];
    AppendTo[SP, {sp, i - j}]
]];
SortBy[SP, First][[All, 2]] (* Jean-François Alcover, Jul 18 2020, after Robert Israel *)
Table[If[!SquareFreeQ[n],0,-Subtract@@PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

lista(nn) = {my(vsp = select(x->(bigomega(x)==2), [1..nn])); vector(#vsp, k, my(f=factor(vsp[k])[,1]); primepi(vecmax(f)) - primepi(vecmin(f)));} \\ Michel Marcus, Jul 18 2020

a(n) = A049084(A084127(n)) - A049084(A084126(n)). [corrected by R. J. Mathar, Apr 22 2010]

a(n) = A338913(n) - A338912(n). - Gus Wiseman, Dec 04 2020

a(51) and a(69) corrected by R. J. Mathar, Apr 22 2010

A338904 Irregular triangle read by rows where row n lists all semiprimes whose prime indices sum to n.

Original entry on oeis.org

4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, 35, 34, 39, 49, 55, 38, 51, 65, 77, 46, 57, 85, 91, 121, 58, 69, 95, 119, 143, 62, 87, 115, 133, 169, 187, 74, 93, 145, 161, 209, 221, 82, 111, 155, 203, 247, 253, 289, 86, 123, 185, 217, 299, 319, 323, 94, 129, 205

Offset: 2

Gus Wiseman, Nov 28 2020

A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

			Triangle begins:
   4
   6
   9  10
  14  15
  21  22  25
  26  33  35
  34  39  49  55
  38  51  65  77
  46  57  85  91 121
  58  69  95 119 143
  62  87 115 133 169 187
  74  93 145 161 209 221
  82 111 155 203 247 253 289
  86 123 185 217 299 319 323
  94 129 205 259 341 361 377 391

A004526 gives row lengths.
A024697 gives row sums.
A087112 is a different triangle of semiprimes.
A098350 has antidiagonals with the same distinct terms as these rows.
A338905 is the squarefree case, with row sums A025129.
A338907/A338906 are the union of odd/even rows.
A339114/A339115 are the row minima/maxima.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A014342 is the self-convolution of primes.
A037143 lists primes and semiprimes.
A056239 gives sum of prime indices (Heinz weight).
A062198 gives partial sums of semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A289182/A115392 list the positions of odd/even terms in A001358.
A332765 gives the greatest squarefree semiprime of weight n.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
Cf. A000040, A001221, A001222, A005117, A112798, A320732, A332877, A338908, A338910, A338911, A339116.

Table[Sort[Table[Prime[k]*Prime[n-k],{k,n/2}]],{n,2,10}]

A338907 Semiprimes whose prime indices sum to an odd number.

Original entry on oeis.org

6, 14, 15, 26, 33, 35, 38, 51, 58, 65, 69, 74, 77, 86, 93, 95, 106, 119, 122, 123, 141, 142, 143, 145, 158, 161, 177, 178, 185, 201, 202, 209, 214, 215, 217, 219, 221, 226, 249, 262, 265, 278, 287, 291, 299, 302, 305, 309, 319, 323, 326, 327, 329, 346, 355

Offset: 1

Gus Wiseman, Nov 28 2020

nonn

All terms are squarefree (A005117).
A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.
The semiprimes in A300063; the semiprimes in A332820. - Peter Munn, Dec 25 2020

			The sequence of terms together with their prime indices begins:
      6: {1,2}      95: {3,8}     202: {1,26}
     14: {1,4}     106: {1,16}    209: {5,8}
     15: {2,3}     119: {4,7}     214: {1,28}
     26: {1,6}     122: {1,18}    215: {3,14}
     33: {2,5}     123: {2,13}    217: {4,11}
     35: {3,4}     141: {2,15}    219: {2,21}
     38: {1,8}     142: {1,20}    221: {6,7}
     51: {2,7}     143: {5,6}     226: {1,30}
     58: {1,10}    145: {3,10}    249: {2,23}
     65: {3,6}     158: {1,22}    262: {1,32}
     69: {2,9}     161: {4,9}     265: {3,16}
     74: {1,12}    177: {2,17}    278: {1,34}
     77: {4,5}     178: {1,24}    287: {4,13}
     86: {1,14}    185: {3,12}    291: {2,25}
     93: {2,11}    201: {2,19}    299: {6,9}

A031368 looks at primes instead of semiprimes.
A098350 has this as union of odd-indexed antidiagonals.
A300063 looks at all numbers (not just semiprimes).
A338904 has this as union of odd-indexed rows.
A338906 is the even version.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A056239 gives the sum of prime indices (Heinz weight).
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A289182/A115392 list the positions of odd/even terms in A001358.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338908 lists squarefree semiprimes of even weight.
A339114/A339115 give the least/greatest semiprime of weight n.
Cf. A000040, A001222, A014342, A024697, A062198, A112798, A300061, A319242, A320655, A338910, A339003.
Subsequence of A332820.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],PrimeOmega[#]==2&&OddQ[Total[primeMS[#]]]&]

from math import isqrt
from sympy import primepi, primerange
def A338907(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum((primepi(x//p)-a>>1) for a,p in enumerate(primerange(isqrt(x)+1)))
    return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

Complement of A338906 in A001358.

A112141 Product of the first n semiprimes.

Original entry on oeis.org

4, 24, 216, 2160, 30240, 453600, 9525600, 209563200, 5239080000, 136216080000, 4495130640000, 152834441760000, 5349205461600000, 203269807540800000, 7927522494091200000, 364666034728195200000, 17868635701681564800000, 911300420785759804800000

Offset: 1

Jonathan Vos Post, Nov 28 2005

Semiprime analog of primorial (A002110). Equivalent for product of what A062198 is for sum.

			a(10) = 4*6*9*10*14*15*21*22*25*26 = 136216080000, the product of the first 10 semiprimes.
From _Gus Wiseman_, Dec 06 2020: (Start)
The sequence of terms together with their prime signatures begins:
                        4: (2)
                       24: (3,1)
                      216: (3,3)
                     2160: (4,3,1)
                    30240: (5,3,1,1)
                   453600: (5,4,2,1)
                  9525600: (5,5,2,2)
                209563200: (6,5,2,2,1)
               5239080000: (6,5,4,2,1)
             136216080000: (7,5,4,2,1,1)
            4495130640000: (7,6,4,2,2,1)
          152834441760000: (8,6,4,2,2,1,1)
         5349205461600000: (8,6,5,3,2,1,1)
       203269807540800000: (9,6,5,3,2,1,1,1)
      7927522494091200000: (9,7,5,3,2,2,1,1)
    364666034728195200000: (10,7,5,3,2,2,1,1,1)
  17868635701681564800000: (10,7,5,5,2,2,1,1,1)
(End)

T. D. Noe, Table of n, a(n) for n = 1..100

Partial sums of semiprimes are A062198.
First differences of semiprimes are A065516.
A000040 lists primes, with partial products A002110 (primorials).
A000142 lists factorials, with partial products A000178 (superfactorials).
A001358 lists semiprimes, with partial products A112141 (this sequence).
A005117 lists squarefree numbers, with partial products A111059.
A006881 lists squarefree semiprimes, with partial products A339191.
A101048 counts partitions into semiprimes (restricted: A338902).
A320655 counts factorizations into semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes, with product/sum/difference A087794/A176504/A176506.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.
Cf. A001222, A084126, A084127, A115392, A168472, A320732, A320892.

A112141 := proc(n)
    mul(A001358(i),i=1..n) ;
end proc:
seq(A112141(n),n=1..10) ; # R. J. Mathar, Jun 30 2020

NextSemiPrime[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sp = n + sgn; While[c < Abs[k], While[ PrimeOmega[sp] != 2, If[sgn < 0, sp--, sp++]]; If[sgn < 0, sp--, sp++]; c++]; sp + If[sgn < 0, 1, -1]]; f[n_] := Times @@ NestList[ NextSemiPrime@# &, 2^2, n - 1]; Array[f, 18] (* Robert G. Wilson v, Jun 13 2013 *)
FoldList[Times,Select[Range[30],PrimeOmega[#]==2&]] (* Gus Wiseman, Dec 06 2020 *)

a(n)=my(v=vector(n),i,k=3);while(iCharles R Greathouse IV, Apr 04 2013

from sympy import factorint
def aupton(terms):
    alst, k, p = [], 1, 1
    while len(alst) < terms:
        if sum(factorint(k).values()) == 2:
            p *= k
            alst.append(p)
        k += 1
    return alst
print(aupton(18)) # Michael S. Branicky, Aug 31 2021

a(n) = Product_{i=1..n} A001358(i).

A001222(a(n)) = 2*n.

A339114 Least semiprime whose prime indices sum to n.

Original entry on oeis.org

4, 6, 9, 14, 21, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, 146, 158, 166, 178, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 386, 394, 398, 422, 446, 454, 458, 466, 478, 482, 502, 514, 526

Offset: 2

Gus Wiseman, Nov 28 2020

nonn

Converges to A100484.
After a(4) = 9, also the least squarefree semiprime whose prime indices sum to n.
A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

			The sequence of terms together with their prime indices begins:
      4: {1,1}     106: {1,16}    254: {1,31}
      6: {1,2}     118: {1,17}    262: {1,32}
      9: {2,2}     122: {1,18}    274: {1,33}
     14: {1,4}     134: {1,19}    278: {1,34}
     21: {2,4}     142: {1,20}    298: {1,35}
     26: {1,6}     146: {1,21}    302: {1,36}
     34: {1,7}     158: {1,22}    314: {1,37}
     38: {1,8}     166: {1,23}    326: {1,38}
     46: {1,9}     178: {1,24}    334: {1,39}
     58: {1,10}    194: {1,25}    346: {1,40}
     62: {1,11}    202: {1,26}    358: {1,41}
     74: {1,12}    206: {1,27}    362: {1,42}
     82: {1,13}    214: {1,28}    382: {1,43}
     86: {1,14}    218: {1,29}    386: {1,44}
     94: {1,15}    226: {1,30}    394: {1,45}

A024697 is the sum of the same semiprimes.
A098350 has this sequence as antidiagonal minima.
A338904 has this sequence as row minima.
A339114 (this sequence) is the squarefree case for n > 4.
A339115 is the greatest among the same semiprimes.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A037143 lists primes and semiprimes.
A056239 gives the sum of prime indices of n.
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A320655 counts factorizations into semiprimes.
A332765/A332877 is the greatest squarefree semiprime of weight n.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338907/A338906 list semiprimes of odd/even weight.
A338907/A338908 list squarefree semiprimes of odd/even weight.
Cf. A000040, A001221, A001222, A014342, A025129, A062198, A112798, A338905, A339116.

Table[Min@@Table[Prime[k]*Prime[n-k],{k,n-1}],{n,2,30}]
Take[DeleteDuplicates[SortBy[{Times@@#,Total[PrimePi[#]]}&/@Tuples[ Prime[ Range[ 200]],2],{Last,First}],GreaterEqual[#1[[2]],#2[[2]]]&][[All,1]],60] (* Harvey P. Dale, Sep 06 2022 *)

a(n) = vecmin(vector(n-1, k, prime(k)*prime(n-k))); \\ Michel Marcus, Dec 03 2020

A025129 a(n) = p(1)p(n) + p(2)p(n-1) + ... + p(k)p(n-k+1), where k = [ n/2 ], p = A000040, the primes.

Original entry on oeis.org

0, 6, 10, 29, 43, 94, 128, 231, 279, 484, 584, 903, 1051, 1552, 1796, 2489, 2823, 3784, 4172, 5515, 6091, 7758, 8404, 10575, 11395, 14076, 15174, 18339, 19667, 23414, 24906, 29437, 31089, 36500, 38614, 44731, 47071, 54198, 56914, 65051, 68371, 77402, 81052, 91341

Offset: 1

Clark Kimberling

nonn

This is the sum of distinct squarefree semiprimes with prime indices summing to n + 1. A squarefree semiprime is a product of any two distinct prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798. - Gus Wiseman, Dec 05 2020

			From _Gus Wiseman_, Dec 05 2020: (Start)
The sequence of sums begins (n > 1):
    6 =  6
   10 = 10
   29 = 14 + 15
   43 = 22 + 21
   94 = 26 + 33 + 35
  128 = 34 + 39 + 55
  231 = 38 + 51 + 65 + 77
  279 = 46 + 57 + 85 + 91
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Gus Wiseman, Sum of prime(i) * prime(j) for i + j = n, i != j.

Cf. A000040, A258323.
The nonsquarefree version is A024697 (shifted right).
Row sums of A338905 (shifted right).
A332765 is the greatest among these squarefree semiprimes.
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A014342 is the self-convolution of the primes.
A056239 is the sum of prime indices of n.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
A339194 sums squarefree semiprimes grouped by greater prime factor.
Cf. A001221, A005117, A062198, A098350, A168472, A320656, A338900, A338901, A338904, A339114, A339116.

a025129 n = a025129_list !! (n-1)
a025129_list= f (tail a000040_list) [head a000040_list] 1 where
   f (p:ps) qs k = sum (take (div k 2) $ zipWith (*) qs $ reverse qs) :
                   f ps (p : qs) (k + 1)
-- Reinhard Zumkeller, Apr 07 2014

f[n_] := Block[{primeList = Prime@ Range@ n}, Total[ Take[ primeList, Floor[n/2]]*Reverse@ Take[ primeList, {Floor[(n + 3)/2], n}]]]; Array[f, 44] (* Robert G. Wilson v, Apr 07 2014 *)

A025129=n->sum(k=1,n\2,prime(k)*prime(n-k+1)) \\ M. F. Hasler, Apr 06 2014

a(n) = A024697(n) for even n. - M. F. Hasler, Apr 06 2014

Following suggestions by Robert Israel and N. J. A. Sloane, initial 0=a(1) added by M. F. Hasler, Apr 06 2014

A332765 Consider all permutations p_i of the first n primes; a(n) is the minimum over p_i of the maximal product of two adjacent primes in the permutation.

Original entry on oeis.org

6, 10, 15, 22, 35, 55, 77, 91, 143, 187, 221, 253, 323, 391, 493, 551, 667, 713, 899, 1073, 1189, 1271, 1517, 1591, 1763, 1961, 2183, 2419, 2537, 2773, 3127, 3233, 3599, 3953, 4189, 4331, 4757, 4897, 5293, 5723, 5963, 6499, 6887, 7171, 7663, 8051, 8633, 8989, 9797, 9991, 10403, 10807

Offset: 2

Bobby Jacobs, Apr 23 2020

nonn

The optimal permutation of n primes is {p_n, p_1, p_n-1, p_2, …, p_ceiling(n/2)}. - Ivan N. Ianakiev, Apr 28 2020

Also the greatest squarefree semiprime whose prime indices sum to n + 1. A squarefree semiprime (A006881) is a product of any two distinct prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798. - Gus Wiseman, Dec 06 2020

			Here are the ways (up to reversal) to order the first four primes:
  2, 3, 5, 7: Products: 6, 15, 35;  Largest product: 35
  2, 3, 7, 5: Products: 6, 21, 35;  Largest product: 35
  2, 5, 3, 7: Products: 10, 15, 21; Largest product: 21
  2, 5, 7, 3: Products: 10, 35, 21; Largest product: 35
  2, 7, 3, 5: Products: 14, 21, 15; Largest product: 21
  2, 7, 5, 3: Products: 14, 35, 15; Largest product: 35
  3, 2, 5, 7: Products: 6, 10, 35;  Largest product: 35
  3, 2, 7, 5: Products: 6, 14, 35;  Largest product: 35
  3, 5, 2, 7: Products: 15, 10, 14; Largest product: 15
  3, 7, 2, 5: Products: 21, 14, 10; Largest product: 21
  5, 2, 3, 7: Products: 10, 6, 21;  Largest product: 21
  5, 3, 2, 7: Products: 15, 6, 14;  Largest product: 15
The minimum largest product is 15, so a(4) = 15.
From _Gus Wiseman_, Dec 06 2020: (Start)
The sequence of terms together with their prime indices begins:
      6: {1,2}     551: {8,10}    3127: {16,17}
     10: {1,3}     667: {9,10}    3233: {16,18}
     15: {2,3}     713: {9,11}    3599: {17,18}
     22: {1,5}     899: {10,11}   3953: {17,19}
     35: {3,4}    1073: {10,12}   4189: {17,20}
     55: {3,5}    1189: {10,13}   4331: {18,20}
     77: {4,5}    1271: {11,13}   4757: {19,20}
     91: {4,6}    1517: {12,13}   4897: {17,23}
    143: {5,6}    1591: {12,14}   5293: {19,22}
    187: {5,7}    1763: {13,14}   5723: {17,25}
    221: {6,7}    1961: {12,16}   5963: {19,24}
    253: {5,9}    2183: {12,17}   6499: {19,25}
    323: {7,8}    2419: {13,17}   6887: {20,25}
    391: {7,9}    2537: {14,17}   7171: {20,26}
    493: {7,10}   2773: {15,17}   7663: {22,25}
(End)

Cf. A332877, A333747.
A338904 and A338905 have this sequence as row maxima.
A339115 is the not necessarily squarefree version.
A001358 lists semiprimes.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A025129 gives the sum of squarefree semiprimes of weight n.
A056239 (weight) gives the sum of prime indices of n.
A320656 counts factorizations into squarefree semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes, with product/sum/difference A087794/A176504/A176506.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.
A338907/A338908 list squarefree semiprimes of odd/even weight.
A339114 is the least (squarefree) semiprime of weight n.
A339116 groups squarefree semiprimes by greater prime factor.
Cf. A001221, A014342, A024697, A046388, A062198, A098350, A112798, A168472, A320655, A338901.

primes[n_]:=Reverse[Prime/@Range[n]]; partition[n_]:=Partition[primes[n],UpTo[Ceiling[n/2]]];
riffle[n_]:=Riffle[partition[n][[1]],Reverse[partition[n][[2]]]];
a[n_]:=Max[Table[riffle[n][[i]]*riffle[n][[i+1]],{i,1,n-1}]];a/@Range[2,53]
(* Ivan N. Ianakiev, Apr 28 2020 *)

It appears that a(n) = A332877(n - 1) for n > 5.

a(12)-a(13) from Jinyuan Wang, Apr 24 2020

More terms from Ivan N. Ianakiev, Apr 28 2020

A339115 Greatest semiprime whose prime indices sum to n.

Original entry on oeis.org

4, 6, 10, 15, 25, 35, 55, 77, 121, 143, 187, 221, 289, 323, 391, 493, 551, 667, 841, 899, 1073, 1189, 1369, 1517, 1681, 1763, 1961, 2183, 2419, 2537, 2809, 3127, 3481, 3599, 3953, 4189, 4489, 4757, 5041, 5293, 5723, 5963, 6499, 6887, 7171, 7663, 8051, 8633

Offset: 2

Gus Wiseman, Nov 28 2020

nonn

A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

			The sequence of terms together with their prime indices begins:
        4: {1,1}      493: {7,10}      2809: {16,16}
        6: {1,2}      551: {8,10}      3127: {16,17}
       10: {1,3}      667: {9,10}      3481: {17,17}
       15: {2,3}      841: {10,10}     3599: {17,18}
       25: {3,3}      899: {10,11}     3953: {17,19}
       35: {3,4}     1073: {10,12}     4189: {17,20}
       55: {3,5}     1189: {10,13}     4489: {19,19}
       77: {4,5}     1369: {12,12}     4757: {19,20}
      121: {5,5}     1517: {12,13}     5041: {20,20}
      143: {5,6}     1681: {13,13}     5293: {19,22}
      187: {5,7}     1763: {13,14}     5723: {17,25}
      221: {6,7}     1961: {12,16}     5963: {19,24}
      289: {7,7}     2183: {12,17}     6499: {19,25}
      323: {7,8}     2419: {13,17}     6887: {20,25}
      391: {7,9}     2537: {14,17}     7171: {20,26}

Robert Israel, Table of n, a(n) for n = 2..10000

A024697 is the sum of the same semiprimes.
A332765/A332877 is the squarefree case.
A338904 has this sequence as row maxima.
A339114 is the least among the same semiprimes.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A037143 lists primes and semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A320655 counts factorizations into semiprimes.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338907/A338906 list semiprimes of odd/even weight.
A338907/A338908 list squarefree semiprimes of odd/even weight.
Cf. A000040, A001221, A001222, A014342, A025129, A056239, A062198, A098350, A112798, A338905, A339116.

P:= [seq(ithprime(i),i=1..200)]:
[seq(max(seq(P[i]*P[j-i],i=1..j-1)),j=2..200)]; # Robert Israel, Dec 06 2020

Table[Max@@Table[Prime[k]*Prime[n-k],{k,n-1}],{n,2,30}]

cp's OEIS Frontend

A001358 Semiprimes (or biprimes): products of two primes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Formula

Extensions

A176506 Difference between the prime indices of the two factors of the n-th semiprime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A338904 Irregular triangle read by rows where row n lists all semiprimes whose prime indices sum to n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A338907 Semiprimes whose prime indices sum to an odd number.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

Formula

A112141 Product of the first n semiprimes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A339114 Least semiprime whose prime indices sum to n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI