A046308 -id:A046308 - cp's OEIS frontend

A046310 Numbers that are divisible by exactly 8 primes counting multiplicity.

256, 384, 576, 640, 864, 896, 960, 1296, 1344, 1408, 1440, 1600, 1664, 1944, 2016, 2112, 2160, 2176, 2240, 2400, 2432, 2496, 2916, 2944, 3024, 3136, 3168, 3240, 3264, 3360, 3520, 3600, 3648, 3712, 3744, 3968, 4000, 4160, 4374, 4416, 4536, 4704, 4736

Offset: 1

Patrick De Geest, Jun 15 1998

nonn

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

Odd terms are in A046321; first odd term is a(64)=6561=3^8. - Zak Seidov, Feb 08 2016

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

Cf. A046309, A120049 (number of 8-almost primes <= 10^n).
Cf. A101637, A101638, 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), A046306 (r = 6), A046308 (r = 7), this sequence (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).
Cf. A046321.

A046310 := proc(n)
    option remember;
    if n = 1 then
        2^8 ;
    else
        for a from procname(n-1)+1 do
            if numtheory[bigomega](a) = 8 then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A046310(n),n=1..30) ; # R. J. Mathar, Dec 21 2018

Select[Range[1600], Plus @@ Last /@ FactorInteger[ # ] == 8 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[5000],PrimeOmega[#]==8&]  (* Harvey P. Dale, Apr 19 2011 *)

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

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A046310(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,8)))
    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 = 8.
a(n) ~ 5040n log n / (log log n)^7. - Charles R Greathouse IV, May 06 2013
a(n) = A078840(8,n). - R. J. Mathar, Jan 30 2019

A046314 Numbers that are divisible by exactly 10 primes with multiplicity.

Original entry on oeis.org

1024, 1536, 2304, 2560, 3456, 3584, 3840, 5184, 5376, 5632, 5760, 6400, 6656, 7776, 8064, 8448, 8640, 8704, 8960, 9600, 9728, 9984, 11664, 11776, 12096, 12544, 12672, 12960, 13056, 13440, 14080, 14400, 14592, 14848, 14976, 15872, 16000, 16640

Offset: 1

Patrick De Geest, Jun 15 1998

nonn

Also called 10-almost primes. Products of exactly 10 primes (not necessarily distinct). Any 10-almost prime can be represented in several ways as a product of two 5-almost primes A014614 and in several ways as a product of five 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, Reference

Cf. A046313, A120051 (number of 10-almost primes <= 10^n).
Cf. A101637, A101638, 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), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), this sequence (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). - Jason Kimberley, Oct 02 2011

Select[Range[5000], Plus @@ Last /@ FactorInteger[ # ] == 10 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[17000],PrimeOmega[#]==10&] (* Harvey P. Dale, Jun 23 2018 *)

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

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A046314(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,10)))
    return bisection(f,n,n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 10.

a(n) ~ 362880n log n / (log log n)^9. - Charles R Greathouse IV, May 06 2013

A046312 Numbers that are divisible by exactly 9 primes with multiplicity.

Original entry on oeis.org

512, 768, 1152, 1280, 1728, 1792, 1920, 2592, 2688, 2816, 2880, 3200, 3328, 3888, 4032, 4224, 4320, 4352, 4480, 4800, 4864, 4992, 5832, 5888, 6048, 6272, 6336, 6480, 6528, 6720, 7040, 7200, 7296, 7424, 7488, 7936, 8000, 8320, 8748, 8832, 9072, 9408

Offset: 1

Patrick De Geest, Jun 15 1998

nonn

Also called 9-almost primes. Products of exactly 9 primes (not necessarily distinct). - Jonathan Vos Post, Dec 11 2004

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

Cf. A046311, A120050 (number of 9-almost primes <= 10^n).
Cf. A101637, A101638, 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), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), this sequence (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). - Jason Kimberley, Oct 02 2011

Select[Range[2200], Plus @@ Last /@ FactorInteger[ # ] == 9 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[10000],PrimeOmega[#]==9&] (* Harvey P. Dale, Oct 24 2020 *)

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

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A046312(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,9)))
    return bisection(f,n,n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 9.

a(n) ~ 40320n log n / (log log n)^8. - Charles R Greathouse IV, May 06 2013

A078840 Table of n-almost-primes T(n,k) (n >= 0, k > 0), read by antidiagonals, starting at T(0,1)=1 followed by T(1,1)=2.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 7, 9, 12, 16, 11, 10, 18, 24, 32, 13, 14, 20, 36, 48, 64, 17, 15, 27, 40, 72, 96, 128, 19, 21, 28, 54, 80, 144, 192, 256, 23, 22, 30, 56, 108, 160, 288, 384, 512, 29, 25, 42, 60, 112, 216, 320, 576, 768, 1024, 31, 26, 44, 81, 120, 224, 432, 640, 1152

Offset: 0

Benoit Cloitre and Paul D. Hanna, Dec 10 2002

An n-almost-prime is a positive integer that has exactly n prime factors.
This sequence is a rearrangement of the natural numbers. - Robert G. Wilson v, Feb 11 2006
Each antidiagonal begins with the n-th prime and ends with 2^n.
From Eric Desbiaux, Jun 27 2009: (Start)
(A001222 gives this sequence)
A001221 gives another table:
  1
  -    2    3    4    5    7    8    9   11 ... A000961
  -    6   10   12   14   15   18   20   21 ... A007774
  -   30   42   60   66   70   78   84   90 ... A033992
  -  210  330  390  420  462  510  546  570 ... A033993
  - 2310 2730 3570 3990 4290 4620 4830 5460 ... A051270
Antidiagonals begin with A000961 and end with A002110.
Diagonal is A073329 which is last term in n-th row of A048692. (End)

			Table begins:
  1
  -  2  3   5   7  11  13  17  19  23  29 ...
  -  4  6   9  10  14  15  21  22  25  26 ...
  -  8 12  18  20  27  28  30  42  44  45 ...
  - 16 24  36  40  54  56  60  81  84  88 ...
  - 32 48  72  80 108 112 120 162 168 176 ...
  - 64 96 144 160 216 224 240 324 336 352 ...

Robert G. Wilson v, Table of n, a(n) for n = 0..10011 (corrected by Ivan Neretin).
Eric Weisstein's World of Mathematics, Almost Prime.

Cf. A078840, A078841, A078842, A078843, A078844, A078445, A078846, A109636.
T(1, k)=A000040(k), T(2, k)=A001358(k), T(3, k)=A014612(k), T(4, k)=A014613(k), T(5, k)=A014614(k), T(6, k)=A046306(k), T(7, k)=A046308(k), T(8, k)=A046310(k), T(9, k)=A046312(k), T(10, k)=A046314(k).
T(11, k)=A069272(k), T(12, k)=A069273(k), T(13, k)=A069274(k), T(14, k)=A069275(k), T(15, k)=A069276(k), T(16, k)=A069277(k), T(17, k)=A069278(k), T(18, k)=A069279(k), T(19, k)=A069280(k), T(20, k)=A069281(k).
T(k, 1)=A000079(k), T(k, 2)=A007283(k), T(k, 3)=A116453(k), T(k, k)=A101695(k), T(k, k+1)=A078841(k).
A091538 is this sequence with zeros inserted, making a square array.

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 *)
AlmostPrime[k_, n_] := Block[{e = Floor[Log[2, n]+k], a, b}, a = 2^e; Do[b = 2^p; While[ AlmostPrimePi[k, a] < n, a = a + b]; a = a - b/2, {p, e, 0, -1}]; a + b/2]; Table[ AlmostPrime[k, n - k + 1], {n, 11}, {k, n}] // Flatten (* Robert G. Wilson v *)
mx = 11; arr = NestList[Take[Union@Flatten@Outer[Times, #, primes], mx] &, primes = Prime@Range@mx, mx]; Prepend[Flatten@Table[arr[[k, n - k + 1]], {n, mx}, {k, n}], 1] (* Ivan Neretin, Apr 30 2016 *)
(* The next code skips the initial 1. *)
width = 15; (seq = Table[
  Rest[NestList[1 + NestWhile[# + 1 &, #, ! PrimeOmega[#] == z &] &,
  2^z, width - z + 1]] - 1, {z, width}]) // TableForm
Flatten[Map[Reverse[Diagonal[Reverse[seq], -width + #]] &, Range[width]]]
(* Peter J. C. Moses, Jun 05 2019 *)
Grid[Table[Select[Range[200], PrimeOmega[#] == n &], {n, 0, 7}]]
(* Clark Kimberling, Nov 17 2024 *)

T(n,k)=if(k<0,0,s=1; while(sum(i=1,s,if(bigomega(i)-n,0,1))

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi, prime
def A078840_T(n,k):
    if n == 1: return prime(k)
    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(k-1+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,n)))
    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

Edited by Robert G. Wilson v, Feb 11 2006

A069272 11-almost primes (generalization of semiprimes).

Original entry on oeis.org

2048, 3072, 4608, 5120, 6912, 7168, 7680, 10368, 10752, 11264, 11520, 12800, 13312, 15552, 16128, 16896, 17280, 17408, 17920, 19200, 19456, 19968, 23328, 23552, 24192, 25088, 25344, 25920, 26112, 26880, 28160, 28800, 29184, 29696

Offset: 1

Rick L. Shepherd, Mar 12 2002

nonn

Product of 11 not necessarily distinct primes.

Divisible by exactly 11 prime powers (not including 1).

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

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), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), this sequence (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). - Jason Kimberley, Oct 02 2011

Select[Range[9000], Plus @@ Last /@ FactorInteger[ # ] == 11 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[30000],PrimeOmega[#]==11&] (* Harvey P. Dale, Jul 13 2013 *)

k=11; start=2^k; finish=30000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

is(n)=bigomega(n)==11 \\ Charles R Greathouse IV, Oct 15 2015

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069272(n):
    def bisection(f, kmin=0, kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x, 0, 1, 1, 11)))
    return bisection(f, n, n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 11.

a(n) ~ 3628800n log n / (log log n)^10. - Charles R Greathouse IV, May 06 2013

A069273 12-almost primes (generalization of semiprimes).

Original entry on oeis.org

4096, 6144, 9216, 10240, 13824, 14336, 15360, 20736, 21504, 22528, 23040, 25600, 26624, 31104, 32256, 33792, 34560, 34816, 35840, 38400, 38912, 39936, 46656, 47104, 48384, 50176, 50688, 51840, 52224, 53760, 56320, 57600, 58368, 59392

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 12 not necessarily distinct primes.
Divisible by exactly 12 prime powers (not including 1).
Any 12-almost prime can be represented in at least one way as a product of two 6-almost primes A046306, three 4-almost primes A014613, four 3-almost primes A014612, or six semiprimes A001358. - Jonathan Vos Post, Dec 11 2004

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

Cf. A101637, A101638, 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), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), this sequence (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). - Jason Kimberley, Oct 02 2011

Select[Range[20000], Plus @@ Last /@ FactorInteger[ # ] == 12 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[60000],PrimeOmega[#]==12&] (* Harvey P. Dale, May 01 2019 *)

k=12; start=2^k; finish=70000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069273(n):
    def bisection(f, kmin=0, kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x, 0, 1, 1, 12)))
    return bisection(f, n, n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 12.

A069279 Products of exactly 18 primes (generalization of semiprimes).

Original entry on oeis.org

262144, 393216, 589824, 655360, 884736, 917504, 983040, 1327104, 1376256, 1441792, 1474560, 1638400, 1703936, 1990656, 2064384, 2162688, 2211840, 2228224, 2293760, 2457600, 2490368, 2555904, 2985984, 3014656, 3096576, 3211264, 3244032, 3317760, 3342336, 3440640

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 18 not necessarily distinct primes.

Divisible by exactly 18 prime powers (not including 1).

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

Cf. A101637, A101638, 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), 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), this sequence (r = 18), A069280 (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

Select[Range[31*10^5],PrimeOmega[#]==18&] (* Harvey P. Dale, Apr 05 2015 *)

k=18; start=2^k; finish=4000000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A069279(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,18)))
    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 = 18.

A069281 20-almost primes (generalization of semiprimes).

Original entry on oeis.org

1048576, 1572864, 2359296, 2621440, 3538944, 3670016, 3932160, 5308416, 5505024, 5767168, 5898240, 6553600, 6815744, 7962624, 8257536, 8650752, 8847360, 8912896, 9175040, 9830400, 9961472, 10223616, 11943936, 12058624

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 20 not necessarily distinct primes.
Divisible by exactly 20 prime powers (not including 1).
Any 20-almost prime can be represented in several ways as a product of two 10-almost primes A046314; in several ways as a product of four 5-almost primes A014614; in several ways as a product of five 4-almost primes A014613; and in several ways as a product of ten semiprimes A001358. - Jonathan Vos Post, Dec 12 2004

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

Cf. A101637, A101638, 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), 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), this sequence (r = 20). - Jason Kimberley, Oct 02 2011

Select[Range[2*9!,5*10! ],Plus@@Last/@FactorInteger[ # ]==20 &] (* Vladimir Joseph Stephan Orlovsky, Feb 26 2009 *)

k=20; start=2^k; finish=15000000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v \\ Depending upon the size of k and how many terms are needed, a much more efficient algorithm than the brute-force method above may be desirable. See additional comments in this section of A069280.

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069281(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,20)))
    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 = 20.

a(n) = A078840(20,n). - R. J. Mathar, Jan 30 2019

A069275 14-almost primes (generalization of semiprimes).

Original entry on oeis.org

16384, 24576, 36864, 40960, 55296, 57344, 61440, 82944, 86016, 90112, 92160, 102400, 106496, 124416, 129024, 135168, 138240, 139264, 143360, 153600, 155648, 159744, 186624, 188416, 193536, 200704, 202752, 207360, 208896, 215040, 225280

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 14 not necessarily distinct primes.
Divisible by exactly 14 prime powers (not including 1).
Any 14-almost prime can be represented in several ways as a product of two 7-almost primes A046308; and in several ways as a product of seven semiprimes A001358. - Jonathan Vos Post, Dec 11 2004

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

Cf. A101637, A101638, 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), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), this sequence(r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

Select[Range[50000], Plus @@ Last /@ FactorInteger[ # ] == 14 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)

k=14; start=2^k; finish=240000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069275(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,14)))
    return bisection(f,n,n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 14.

A069276 15-almost primes (generalization of semiprimes).

Original entry on oeis.org

32768, 49152, 73728, 81920, 110592, 114688, 122880, 165888, 172032, 180224, 184320, 204800, 212992, 248832, 258048, 270336, 276480, 278528, 286720, 307200, 311296, 319488, 373248, 376832, 387072, 401408, 405504, 414720, 417792, 430080

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 15 not necessarily distinct primes.
Divisible by exactly 15 prime powers (not including 1).
Any 15-almost prime can be represented in several ways as a product of three 5-almost primes A014614, and in several ways as a product of five 3-almost primes A014612. - Jonathan Vos Post, Dec 11 2004

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

Cf. A101637, A101638, 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), 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), this sequence (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

Select[Range[90000], Plus @@ Last /@ FactorInteger[ # ] == 15 &] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2008 *)
Select[Range[450000],PrimeOmega[#]==15&] (* Harvey P. Dale, Aug 14 2019 *)

k=15; start=2^k; finish=500000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069276(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    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+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,15)))
    return bisection(f,n,n) # Chai Wah Wu, Nov 03 2024

Product p_i^e_i with Sum e_i = 15.

cp's OEIS Frontend

A046310 Numbers that are divisible by exactly 8 primes counting multiplicity.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A046314 Numbers that are divisible by exactly 10 primes with multiplicity.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A046312 Numbers that are divisible by exactly 9 primes with multiplicity.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A078840 Table of n-almost-primes T(n,k) (n >= 0, k > 0), read by antidiagonals, starting at T(0,1)=1 followed by T(1,1)=2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Extensions

A069272 11-almost primes (generalization of semiprimes).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

A069273 12-almost primes (generalization of semiprimes).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A069279 Products of exactly 18 primes (generalization of semiprimes).

Views

Author

Keywords