%I #47 Aug 31 2024 18:02:56
%S 9699690,11741730,13123110,14804790,15825810,16546530,17160990,
%T 17687670,18888870,20030010,20281170,20930910,21111090,21411390,
%U 21637770,21951930,23130030,23393370,23993970,24534510,25555530,25571910
%N Products of 8 distinct primes (squarefree 8-almost primes).
%C Intersection of A005117 and A046310.
%H Robert G. Wilson v, <a href="/A123322/b123322.txt">Table of n, a(n) for n = 1..29422</a> (first 10000 terms from Rick Shepherd)
%H <a href="/index/Pri#prime_signature">Index to sequences related to prime signature</a>
%e a(1) = 9699690 = 2*3*5*7*11*13*17*19 = A002110(8).
%p N:= 3*10^7: # to get all terms <= N
%p pmax:= floor(N/mul(ithprime(i),i=1..7)):
%p Primes:= select(isprime,[2,seq(i,i=3..pmax,2)]):
%p sort(select(`<`,map(convert,combinat:-choose(Primes,8),`*`),N)); # _Robert Israel_, Dec 18 2018
%t f8Q[n_]:=Last/@FactorInteger[n]=={1, 1, 1, 1, 1, 1, 1, 1}; lst={};Do[If[f8Q[n], AppendTo[lst, n]], {n, 10!, 11!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Aug 26 2008 *)
%t Take[ Sort[ Times @@@ Subsets[ Prime@ Range@ 15, {8}]], 22] (* _Robert G. Wilson v_, Dec 18 2018 *)
%o (PARI) is(n)=issquarefree(n)&&omega(n)==8 \\ _Charles R Greathouse IV_, Feb 01 2017, corrected (following an observation from Zak Seidov) by _M. F. Hasler_, Dec 19 2018
%o (PARI) is(n) = my(f = factor(n)); omega(f) == 8 && bigomega(f) == 8 \\ _David A. Corneth_, Dec 18 2018
%o (Python)
%o from math import isqrt, prod
%o from sympy import primerange, integer_nthroot, primepi
%o def A123322(n):
%o def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b+1,isqrt(x//c)+1),a+1)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b+1,integer_nthroot(x//c,m)[0]+1),a+1) for d in g(x,a2,b2,c*b2,m-1)))
%o 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,8)))
%o def bisection(f,kmin=0,kmax=1):
%o while f(kmax) > kmax: kmax <<= 1
%o while kmax-kmin > 1:
%o kmid = kmax+kmin>>1
%o if f(kmid) <= kmid:
%o kmax = kmid
%o else:
%o kmin = kmid
%o return kmax
%o return bisection(f) # _Chai Wah Wu_, Aug 31 2024
%Y Cf. A001221, A001222, A005117, A046310, A048692, Squarefree k-almost primes: A000040 (k=1), A006881 (k=2), A007304 (k=3), A046386 (k=4), A046387 (k=5), A067885 (k=6), A123321 (k=7), A115343 (k=9).
%K nonn
%O 1,1
%A _Rick L. Shepherd_, Sep 25 2006
%E Edited by _Robert Israel_, Dec 18 2018