login
Products of 8 distinct primes (squarefree 8-almost primes).
8

%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