A065516 -id:A065516 - cp's OEIS frontend

A085809 Indices of semiprimes where largest gap occurs. Or, positions of records in A065516.

1, 2, 4, 6, 10, 34, 186, 422, 760, 1765, 4112, 4585, 8112, 8650, 8861, 75150, 223993, 327048, 712605, 1135940, 23958638, 42367759, 47848742, 54626559, 121984495, 157877985, 413509327, 798321315, 983679985, 1277946119, 2403158480

Offset: 1

Jason Earls, Jul 24 2003

nonn

A215231(n) = A065516(a(n)). - Reinhard Zumkeller, Mar 23 2014

Cf. A239674.

a085809 n = a085809_list !! (n-1)
-- See A215231 for definition of a085809_list.
-- Reinhard Zumkeller, Mar 23 2014

{sgr(m)=local(a,b,rec,c); c=0; a=0; b=4; rec=0; for(n=5,m,if(bigomega(n)==2,c++; a=n; if(a-b>rec,rec=a-b; print1(c","); b=a,b=a; )))}

a(19)-a(26) from Donovan Johnson, Jan 28 2009

a(27)-a(31) from Donovan Johnson, Apr 14 2010

A123375 Least k such that the difference between consecutive semiprimes A065516(k) equals n, or 0 if no such k exists.

Original entry on oeis.org

3, 1, 2, 4, 24, 6, 10, 56, 50, 78, 34, 320, 249, 186, 463, 762, 598, 1238, 422, 760, 3760, 3585, 9214, 1765, 4112, 13447, 6675, 4585, 68498, 8112, 10083, 8650, 86203, 49433, 35085, 20641, 458421, 8861, 366314, 157857, 169147, 487115, 277440, 563951, 511757, 920602, 75150

Offset: 1

Alexander Adamchuk, Nov 09 2006

a(n) equals least k such that A065516(k) = n, or 0 if no such k exists. Conjecture: a(n)>0 exists for all n.

			A065516(n) begins {2, 3, 1, 4, 1, 6, 1, 3, 1, 7, 1, 1, 3, 1, 7, 3, 2, 4, 2, 1, 4, 3, 4, 5, ...}.
Thus a(1) = 3, a(2) = 1, a(3) = 2, a(4) = 4, a(5) = 24.

David A. Corneth, Table of n, a(n) for n = 1..82
Eric Weisstein's World of Mathematics, Semiprime

Cf. A001358 (semiprimes), A065516 (differences between semiprimes), A131109.

nextSprime := proc(n) local res ; res := n+1 ; while numtheory[bigomega](res) <> 2 do res := res+1 ; od ; RETURN(res) ; end ; nmax := 500 ; kmax := 500000 ; a := array(1..nmax) ; for i from 1 to nmax do a[i] := 0 ; od : sp1 := 4 : sp2 := nextSprime(sp1) : n := sp2-sp1 : a[n] := 1 : for k from 2 to kmax do sp1 := sp2 ; sp2 := nextSprime(sp1) ; n := sp2-sp1 ; if a[n] = 0 then a[n] := k ; fi ; od : for i from 1 to nmax do if a[i] = 0 then break ; else printf("%d,",a[i]) ; fi ; od : # R. J. Mathar, Jan 13 2007

Table[k=6;While[FreeQ[b=Differences[Select[Range@k++,PrimeOmega[#]==2&]],n]];
Length@b,{n,11}] (* Giorgos Kalogeropoulos, Apr 02 2021 *)

Corrected and extended by R. J. Mathar, Jan 13 2007

More terms from David A. Corneth, Apr 02 2021

A174479 Triangle read by rows: T(n,k) = semiprime(n) mod A065516(k), where A065516 are differences between consecutive semiprimes.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 2, 0, 2, 0, 1, 0, 0, 3, 0, 3, 1, 0, 0, 1, 0, 3, 0, 0, 1, 0, 2, 0, 4, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 2, 0, 2, 0, 2, 0, 2, 0, 5, 1, 0, 0, 1, 0, 3, 0, 0, 0, 5, 0, 0, 1, 0, 2, 0, 4, 0, 1, 0, 6, 0, 0, 1, 2, 0, 3, 0, 5, 0, 2, 0, 0, 0, 0, 2

Offset: 1

Juri-Stepan Gerasimov, Nov 28 2010

			The triangle starts in row n=1 with columns 1 <= k <= n:
  0;
  0, 0;
  1, 0, 0;
  0, 1, 0, 2;
  0, 2, 0, 2, 0;
  1, 0, 0, 3, 0, 3;

Cf. A001358, A174433.

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:
A065516 := proc(n) A001358(n+1)-A001358(n) ;end proc:
A174479 := proc(n,k) A001358(n) mod A065516(k) ; end proc:

A356769 Semiprime gaps (A065516) in the order of first occurrences.

Original entry on oeis.org

2, 3, 1, 4, 6, 7, 5, 11, 9, 8, 10, 14, 13, 12, 19, 15, 17, 20, 16, 18, 24, 22, 21, 25, 28, 27, 30, 32, 38, 23, 31, 26, 36, 35, 34, 29, 47, 33, 40, 41, 54, 50, 43, 55, 39, 48, 37, 42, 45, 44, 53, 70, 46, 56, 74, 52, 62, 51, 66, 49, 58, 68, 59, 63, 67, 60, 57, 61, 72, 64, 65, 76, 69, 73, 75, 82, 85

Offset: 1

Zak Seidov, Aug 27 2022

nonn

Cf. A065516, A001358.

s={2};n=6;Do[m=n+1;While[PrimeOmega[m]!= 2,m++];If[FreeQ[s,m-n],AppendTo[s,m-n];Print[{n,m-n}]];n=m,{10^8}];s

A176464 The positions of primes in A065516.

Original entry on oeis.org

1, 2, 8, 10, 13, 15, 16, 17, 19, 22, 24, 25, 26, 27, 31, 34, 35, 37, 39, 45, 48, 51, 53, 54, 55, 58, 59, 60, 61, 62, 65, 67, 71, 74, 75, 77, 79, 82, 83, 84, 85, 86, 89, 90, 92, 94, 97, 99, 100, 101, 102, 103, 105, 106, 109, 110, 112, 113, 115, 121, 122, 125, 126, 129, 130

Offset: 1

Juri-Stepan Gerasimov, Apr 18 2010

nonn

Prime differences between products of 2 primes.

read("transforms3") : L := BFILETOLIST("b065516.txt") : for a from 1 to 200 do if isprime(op(a,L)) then printf("%d,", a) ; end if; end do: # R. J. Mathar, Apr 20 2010

A065516(a(n) ) =A000040(k).

Entries checked by R. J. Mathar, Apr 20 2010

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

Original entry on oeis.org

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

A338899 Concatenated sequence of prime indices of squarefree semiprimes (A006881).

Original entry on oeis.org

1, 2, 1, 3, 1, 4, 2, 3, 2, 4, 1, 5, 1, 6, 2, 5, 1, 7, 3, 4, 1, 8, 2, 6, 1, 9, 2, 7, 3, 5, 2, 8, 1, 10, 1, 11, 3, 6, 2, 9, 1, 12, 4, 5, 1, 13, 3, 7, 1, 14, 2, 10, 4, 6, 2, 11, 1, 15, 3, 8, 1, 16, 2, 12, 3, 9, 1, 17, 4, 7, 1, 18, 2, 13, 2, 14, 4, 8, 1, 19, 2, 15

Offset: 1

Gus Wiseman, Nov 16 2020

This is a triangle with two columns and strictly increasing rows, namely {A270650(n), A270652(n)}.

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.

			The sequence of terms together with their prime indices begins:
      6: {1,2}     57: {2,8}     106: {1,16}    155: {3,11}
     10: {1,3}     58: {1,10}    111: {2,12}    158: {1,22}
     14: {1,4}     62: {1,11}    115: {3,9}     159: {2,16}
     15: {2,3}     65: {3,6}     118: {1,17}    161: {4,9}
     21: {2,4}     69: {2,9}     119: {4,7}     166: {1,23}
     22: {1,5}     74: {1,12}    122: {1,18}    177: {2,17}
     26: {1,6}     77: {4,5}     123: {2,13}    178: {1,24}
     33: {2,5}     82: {1,13}    129: {2,14}    183: {2,18}
     34: {1,7}     85: {3,7}     133: {4,8}     185: {3,12}
     35: {3,4}     86: {1,14}    134: {1,19}    187: {5,7}
     38: {1,8}     87: {2,10}    141: {2,15}    194: {1,25}
     39: {2,6}     91: {4,6}     142: {1,20}    201: {2,19}
     46: {1,9}     93: {2,11}    143: {5,6}     202: {1,26}
     51: {2,7}     94: {1,15}    145: {3,10}    203: {4,10}
     55: {3,5}     95: {3,8}     146: {1,21}    205: {3,13}

A270650 is the first column.
A270652 is the second column.
A320656 counts multiset partitions using these rows, or factorizations into squarefree semiprimes.
A338898 is the version including squares, with columns A338912 and A338913.
A338900 gives row differences.
A338901 gives the row numbers for first appearances.
A001221 and A001222 count distinct/all prime indices.
A001358 lists semiprimes.
A004526 counts 2-part partitions, with strict case shifted right once.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A046315 and A100484 list odd and even semiprimes.
A046388 lists odd squarefree semiprimes.
A166237 gives first differences of squarefree semiprimes.
Cf. A030229, A056239, A065516, A112798, A115392, A167171, A176506, A320893, A338904, A338906, A338907, A338910, A338911.

Join@@Cases[Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&],k_:>PrimePi/@First/@FactorInteger[k]]

A338913 Greater prime index of the n-th semiprime.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 20 2020

nonn

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

After the first three terms, there appear to be no adjacent equal terms.

			The semiprimes are:
  2*2, 2*3, 3*3, 2*5, 2*7, 3*5, 3*7, 2*11, 5*5, 2*13, ...
so the greater prime factors are:
  2, 3, 3, 5, 7, 5, 7, 11, 5, 13, ...
with indices:
  1, 2, 2, 3, 4, 3, 4, 5, 3, 6, ...

Amiram Eldar, Table of n, a(n) for n = 1..10000

A115392 lists positions of first appearances of each positive integer.
A270652 is the squarefree case, with lesser part A270650.
A338898 has this as second column.
A338912 is the corresponding lesser prime index.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A087794/A176504/A176506 are product/sum/difference of semiprime indices.
A338910/A338911 list products of pairs of odd/even-indexed primes.
Cf. A037143, A056239, A065516, A084126, A084127, A112798, A128301, A289182, A320732, A320892, A338900, A338904, A338909.

Table[Max[PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}]

a(n) = A000720(A084127(n)).

A338898 Concatenated sequence of prime indices of semiprimes (A001358).

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 1, 3, 1, 4, 2, 3, 2, 4, 1, 5, 3, 3, 1, 6, 2, 5, 1, 7, 3, 4, 1, 8, 2, 6, 1, 9, 4, 4, 2, 7, 3, 5, 2, 8, 1, 10, 1, 11, 3, 6, 2, 9, 1, 12, 4, 5, 1, 13, 3, 7, 1, 14, 2, 10, 4, 6, 2, 11, 1, 15, 3, 8, 1, 16, 2, 12, 3, 9, 1, 17, 4, 7, 5, 5, 1, 18, 2

Offset: 1

Gus Wiseman, Nov 15 2020

This is a triangle with two columns and weakly increasing rows, namely {A338912(n), A338913(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 semiprimes together with their prime indices begins:
      4: {1,1}     46: {1,9}      91: {4,6}     141: {2,15}
      6: {1,2}     49: {4,4}      93: {2,11}    142: {1,20}
      9: {2,2}     51: {2,7}      94: {1,15}    143: {5,6}
     10: {1,3}     55: {3,5}      95: {3,8}     145: {3,10}
     14: {1,4}     57: {2,8}     106: {1,16}    146: {1,21}
     15: {2,3}     58: {1,10}    111: {2,12}    155: {3,11}
     21: {2,4}     62: {1,11}    115: {3,9}     158: {1,22}
     22: {1,5}     65: {3,6}     118: {1,17}    159: {2,16}
     25: {3,3}     69: {2,9}     119: {4,7}     161: {4,9}
     26: {1,6}     74: {1,12}    121: {5,5}     166: {1,23}
     33: {2,5}     77: {4,5}     122: {1,18}    169: {6,6}
     34: {1,7}     82: {1,13}    123: {2,13}    177: {2,17}
     35: {3,4}     85: {3,7}     129: {2,14}    178: {1,24}
     38: {1,8}     86: {1,14}    133: {4,8}     183: {2,18}
     39: {2,6}     87: {2,10}    134: {1,19}    185: {3,12}

A112798 restricted to rows of length 2 gives this triangle.
A115392 is the row number for the first appearance of each positive integer.
A176506 gives row differences.
A338899 is the squarefree version.
A338912 is column 1.
A338913 is column 2.
A001221 counts a number's distinct prime indices.
A001222 counts a number's prime indices.
A001358 lists semiprimes.
A004526 counts 2-part partitions.
A006881 lists squarefree semiprimes.
A037143 lists primes and semiprimes.
A046315 and A100484 list odd and even semiprimes.
A046388 and A100484 list odd and even squarefree semiprimes.
A065516 gives first differences of semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A270650 and A270652 give the prime indices of squarefree semiprimes.
A320655 counts factorizations into semiprimes.
Cf. A056239, A101048, A320892, A320912, A338900, A338901, A338904, A338906, A338907, A338910, A338911.

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

A338912 Lesser prime index of the n-th semiprime.

Original entry on oeis.org

1, 1, 2, 1, 1, 2, 2, 1, 3, 1, 2, 1, 3, 1, 2, 1, 4, 2, 3, 2, 1, 1, 3, 2, 1, 4, 1, 3, 1, 2, 4, 2, 1, 3, 1, 2, 3, 1, 4, 5, 1, 2, 2, 4, 1, 2, 1, 5, 3, 1, 3, 1, 2, 4, 1, 6, 2, 1, 2, 3, 5, 1, 2, 1, 4, 3, 1, 5, 2, 1, 3, 4, 1, 2, 6, 1, 3, 2, 6, 2, 5, 1, 4, 1, 3, 2, 1

Offset: 1

Gus Wiseman, Nov 20 2020

nonn

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

			The semiprimes are:
  2*2, 2*3, 3*3, 2*5, 2*7, 3*5, 3*7, 2*11, 5*5, 2*13, ...
so the lesser prime factors are:
  2, 2, 3, 2, 2, 3, 3, 2, 5, 2, ...
with indices:
  1, 1, 2, 1, 1, 2, 2, 1, 3, 1, ...

Amiram Eldar, Table of n, a(n) for n = 1..10000

A084126 is the lesser prime factor (not index).
A084127 is the greater factor, with index A338913.
A115392 lists positions of ones.
A128301 lists positions of first appearances of each positive integer.
A270650 is the squarefree case, with greater part A270652.
A338898 has this as first column.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes, with odds A046315 and evens A100484.
A006881 lists squarefree semiprimes, with odds A046388 and evens A100484.
A087794/A176504/A176506 are product/sum/difference of semiprime indices.
A338910/A338911 list products of pairs of odd/even-indexed primes.
Cf. A037143, A056239, A065516, A112798, A289182, A320732, A320892, A338900, A338904, A338909.

Table[Min[PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}]

a(n) = A000720(A084126(n)).

cp's OEIS Frontend

A085809 Indices of semiprimes where largest gap occurs. Or, positions of records in A065516.

Views

Author

Keywords

Comments

Crossrefs

Programs

Haskell

PARI

Extensions

A123375 Least k such that the difference between consecutive semiprimes A065516(k) equals n, or 0 if no such k exists.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A174479 Triangle read by rows: T(n,k) = semiprime(n) mod A065516(k), where A065516 are differences between consecutive semiprimes.

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

A356769 Semiprime gaps (A065516) in the order of first occurrences.

Views

Author

Keywords

Crossrefs

Programs

Mathematica

A176464 The positions of primes in A065516.

Views

Author

Keywords

Comments

Programs

Maple

Formula

Extensions

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

A338899 Concatenated sequence of prime indices of squarefree semiprimes (A006881).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A338913 Greater prime index of the n-th semiprime.

Views

Author