A120047 -id:A120047 - cp's OEIS frontend

A046306 Numbers that are divisible by exactly 6 primes with multiplicity.

64, 96, 144, 160, 216, 224, 240, 324, 336, 352, 360, 400, 416, 486, 504, 528, 540, 544, 560, 600, 608, 624, 729, 736, 756, 784, 792, 810, 816, 840, 880, 900, 912, 928, 936, 992, 1000, 1040, 1104, 1134, 1176, 1184, 1188, 1215, 1224, 1232, 1260, 1312, 1320

Offset: 1

Patrick De Geest, Jun 15 1998

Also called 6-almost primes. Products of exactly 6 primes (not necessarily distinct). Any 6-almost prime can be represented in several ways as a product of two 3-almost primes A014612 and in several ways as a product of three semiprimes A001358. - Jonathan Vos Post, Dec 11 2004

T. D. Noe, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Almost Prime

Cf. A046305, A120047 (number of 6-almost primes <= 10^n).
Cf. A101605, A101606.
Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), this sequence (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).

Select[Range[500], Plus @@ Last /@ FactorInteger[ # ] == 6 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[1400],PrimeOmega[#]==6&] (* Harvey P. Dale, May 21 2012 *)

is(n)=bigomega(n)==6 \\ Charles R Greathouse IV, Mar 21 2013

from math import isqrt, prod
from sympy import primepi, primerange, integer_nthroot
def A046306(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n-1+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,6)))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 23 2024

Product p_i^e_i with Sum e_i = 6.
a(n) ~ 120n log n / (log log n)^5. - Charles R Greathouse IV, May 06 2013
a(n) = A078840(6,n). - R. J. Mathar, Jan 30 2019

A116430 The number of n-almost primes less than or equal to 10^n, starting with a(0)=1.

Original entry on oeis.org

1, 4, 34, 247, 1712, 11185, 68963, 409849, 2367507, 13377156, 74342563, 407818620, 2214357712, 11926066887, 63809981451, 339576381990, 1799025041767, 9494920297227, 49950199374227, 262036734664892

Offset: 0

Robert G. Wilson v, Feb 10 2006, Jun 01 2006

nonn

If instead we asked for those less than or equal to 2^n, then the sequence is A000012.

Cf. A036352, A114106, A114453.
Cf. A078840, A078841, A078842, A116432, A078843, A116426, A078844, A116427, A078845, A116428, A116429, A116430, A078846, A116431.
Cf. A006880, A036352, A109251, A114106, A114453, A120047 - A120053.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[ AlmostPrimePi[n, 10^n], {n, 0, 13}]

almost_prime_count(N, k) = if(k==1, return(primepi(N))); (f(m, p, k, j=0) = my(c=0, s=sqrtnint(N\m, k)); if(k==2, forprime(q=p, s, c += primepi(N\(m*q))-j; j += 1), forprime(q=p, s, c += f(m*q, q, k-1, j); j += 1)); c); f(1, 2, k);
a(n) = if(n == 0, 1, almost_prime_count(10^n, n)); \\ Daniel Suteu, Jul 10 2023

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A116430(n):
    if n<=1: return 3*n+1
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    return int(sum(primepi(10**n//prod(c[1] for c in a))-a[-1][0] for a in g(10**n,0,1,1,n))) # Chai Wah Wu, Aug 23 2024

Edited by N. J. A. Sloane, Aug 08 2008 at the suggestion of R. J. Mathar
a(15)-a(16) from Donovan Johnson, Oct 01 2010
a(17)-a(19) from Daniel Suteu, Jul 10 2023

A120049 Number of 8-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 7, 105, 1418, 17572, 207207, 2367507, 26483012, 291646797, 3173159326, 34192782745, 365561221293, 3882841742380, 41015564702074, 431227959019552, 4515480975731045, 47115876816676830

Offset: 0

Robert G. Wilson v, Feb 07 2006

nonn

			There are 7 eight-almost primes up to 1000: 256, 384, 576, 640, 864, 896 & 960.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[8, 10^n], {n, 12}]

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A120049(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    return int(sum(primepi(10**n//prod(c[1] for c in a))-a[-1][0] for a in g(10**n,0,1,1,8))) # Chai Wah Wu, Aug 23 2024

a(13)-a(14) from Robert G. Wilson v, Jan 07 2007
Example corrected by Harvey P. Dale, Aug 13 2018
a(15)-a(18) from Henri Lifchitz, Mar 18 2025

A120048 Number of 7-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 14, 231, 2973, 35585, 409849, 4600247, 50678212, 550454756, 5913771637, 62981797962, 665997804082, 7001087934965, 73232029374751, 762783057783010, 7916319351632036, 81898808371556517

Offset: 0

Robert G. Wilson v, Feb 07 2006

nonn

			There are 14 seven-almost primes up to 1000: 128, 192, 288, 320, 432, 448, 480, 648, 672, 704, 720, 800, 832 & 972.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[7, 10^n], {n, 11}]

More terms from Robert G. Wilson v, Jan 07 2007
Example corrected by Harvey P. Dale, Jan 25 2013
a(15)-a(18) from Henri Lifchitz, Mar 18 2025

A120050 Number of 9-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 2, 47, 671, 8491, 101787, 1180751, 13377156, 148930536, 1636170477, 17787688377, 191742524399, 2052389350029, 21838745177567, 231206458686127, 2437121982958248, 25591920108631224, 267840642082525459

Offset: 0

Robert G. Wilson v, Feb 07 2006

nonn

			There are 2 nine-almost primes up to 1000: 512 & 768.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[9, 10^n], {n, 12}]

a(13) and a(14) from Robert G. Wilson v, Jan 07 2007

a(15)-a(19) from Henri Lifchitz, Mar 18 2025

A120051 Number of 10-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 0, 22, 306, 4016, 49163, 578154, 6618221, 74342563, 823164388, 9011965866, 97765974368, 1052666075366, 11263041623194, 119864659464824, 1269754732725522, 13396817167474205, 140847445420555406

Offset: 0

Robert G. Wilson v, Feb 07 2006

			There are 22 ten-almost primes up to 10000: 1024, 1536, 2304, 2560, 3456, 3584, 3840, 5184, 5376, 5632, 5760, 6400, 6656, 7776, 8064, 8448, 8640, 8704, 8960, 9600, 9728, and 9984.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[10, 10^n], {n, 12}]

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A120051(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    return int(sum(primepi(10**n//prod(c[1] for c in a))-a[-1][0] for a in g(10**n,0,1,1,10))) # Chai Wah Wu, Nov 03 2024

More terms from Robert G. Wilson v, Jan 07 2007

a(15)-a(19) from Henri Lifchitz, Mar 20 2025

A120053 Number of 12-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 0, 3, 63, 865, 11068, 133862, 1563465, 17836903, 200051717, 2214357712, 24255601105, 263439785143, 2841076717752, 30457549169277, 324855769153426, 3449587218984911, 36489283363168885

Offset: 0

Robert G. Wilson v, Feb 07 2006

			There are 3 twelve-almost primes up to 10000: 4096, 6144, and 9216.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[12, 10^n], {n, 11}]

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A120053(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    return int(sum(primepi(10**n//prod(c[1] for c in a))-a[-1][0] for a in g(10**n,0,1,1,12))) # Chai Wah Wu, Aug 23 2024

a(13) and a(14) from Robert G. Wilson v, Jan 07 2007
a(15) from Chai Wah Wu, Aug 24 2024
a(16)-a(19) from Henri Lifchitz, Mar 18 2025

A120052 Number of 11-almost primes less than or equal to 10^n.

Original entry on oeis.org

0, 0, 0, 0, 7, 138, 1878, 23448, 279286, 3230577, 36585097, 407818620, 4490844534, 48972151631, 529781669333, 5693047157230, 60832290450373, 646862625625663, 6849459596884350, 72259172519243461

Offset: 0

Robert G. Wilson v, Feb 07 2006

			There are 7 eleven-almost primes up to 10000: 2048, 3072, 4608, 5120, 6912, 7168, and 7680.

Cf. A066265, A014613, A116382.

Cf. A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053, A116430.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[AlmostPrimePi[11, 10^n], {n, 12}]

a(14) from Robert G. Wilson v, Jan 07 2007

a(15)-a(19) from Henri Lifchitz, Mar 18 2025

A126280 Triangle read by rows: T(k,n) is number of numbers <= 10^n that are products of k primes.

Original entry on oeis.org

4, 4, 1, 25, 34, 22, 12, 4, 2, 168, 299, 247, 149, 76, 37, 14, 7, 2, 1229, 2625, 2569, 1712, 963, 485, 231, 105, 47, 22, 7, 3, 1, 9592, 23378, 25556, 18744, 11185, 5933, 2973, 1418, 671, 306, 138, 63, 25, 11, 4, 2, 78498, 210035, 250853, 198062, 124465, 68963

Offset: 1

Jonathan Vos Post & Robert G. Wilson v, Dec 22 2006

The n-th row's sum is 10^n - 1.

			4 4 1
25 34 22 12 4 2
168 299 247 149 76 37 14
7 2
1229 2625 2569 1712 963 485 231
105 47 22 7 3 1
9592 23378 25556 18744 11185 5933 2973
1418 671 306 138 63 25 11 4 2
78498 210035 250853 198062 124465 68963 35585 17572
8491 4016 1878 865 400 179 79 35 14 7 2
664579 1904324 2444359 2050696 1349779 774078 409849 207207
101787 49163 23448 11068 5210 2406 1124 510 233 102 45 21 7 3 1

Robert G. Wilson v, Table of n, a(n) for n = 1..342
Eric Weisstein's World of Mathematics, Semiprime.
Eric Weisstein's World of Mathematics, Almost Prime.

By columns: A006880, A036352, A109251, A114106, A114453, A120047, A120048, A120049, A120050, A120051, A120052, A120053.

The n-th row's sum: A002283 = 10^n -1, A116430, A126279: same array but for powers of two.

AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[ PrimePi[ n/Times @@ Prime[ Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[ Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]] ]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[ AlmostPrimePi[m, 10^n], {n, 6}, {m, Floor[Log[2, 10^n]] }] // Flatten

A124033 Number of n-digit numbers having exactly n prime factors (counted with multiplicity).

Original entry on oeis.org

4, 31, 225, 1563, 10222, 63030, 374264, 2160300, 12196405, 67724342, 371233523, 2014305995, 10841722966, 57974736592, 308361428628, 1632877406997

Offset: 1

J. M. Bergot, Apr 08 2011

Essentially the same as A036335.

What would be the ratio between a(n) and all possible numbers with n digits for each n?

			a(1) = A006880(1) = 4.
a(2) = A066265(2) - A066265(1) = 34 - 3 = 31.
a(3) = A109251(3) - A109251(2) = 247 - 22 = 225.
a(4) = A114106(4) - A114106(3) = 1712 - 149 = 1563.
a(5) = A114453(5) - A114453(4) = 11185 - 963 = 10222.
a(6) = A120047(6) - A120047(5) = 68963 - 5933 = 63030.
a(7) = A120048(7) - A120048(6) = 409849 - 35585 = 374264.
a(8) = A120049(8) - A120049(7) = 2367507 - 207207 = 2160300.
a(9) = A120050(9) - A120050(8) = 13377156 - 1180751 = 12196405.
a(10) = A120051(10) - A120051(9) = 74342563 - 6618221 = 67724342.
a(11) = A120052(11) - A120052(10) = 407818620 - 36585097 = 371233523.
a(12) = A120053(12) - A120053(11) = 2214357712 - 200051717 = 2014305995.

Table[Count[Range[10^(n-1),10^n-1],?(PrimeOmega[#]==n&)],{n,8}]  (* _Harvey P. Dale, Apr 22 2011 *)
AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
f[n_] := AlmostPrimePi[n, 10^n - 1] - AlmostPrimePi[n, 10^(n - 1) - 1]; Array[f, 12] (* Robert G. Wilson v, Jul 06 2012 *)

Corrected and extended by Ray Chandler, Apr 11 2011
a(9)-a(12) from Ray Chandler, Apr 12 2011
a(13)-a(16) from Robert G. Wilson v, Jul 06 2012

cp's OEIS Frontend

A046306 Numbers that are divisible by exactly 6 primes with multiplicity.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A116430 The number of n-almost primes less than or equal to 10^n, starting with a(0)=1.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

Python

Extensions

A120049 Number of 8-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Python

Extensions

A120048 Number of 7-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Extensions

A120050 Number of 9-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Extensions

A120051 Number of 10-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Python

Extensions

A120053 Number of 12-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Python

Extensions

A120052 Number of 11-almost primes less than or equal to 10^n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica