login
Product of the first n semiprimes.
22

%I #30 Aug 31 2021 21:21:37

%S 4,24,216,2160,30240,453600,9525600,209563200,5239080000,136216080000,

%T 4495130640000,152834441760000,5349205461600000,203269807540800000,

%U 7927522494091200000,364666034728195200000,17868635701681564800000,911300420785759804800000

%N Product of the first n semiprimes.

%C Semiprime analog of primorial (A002110). Equivalent for product of what A062198 is for sum.

%H T. D. Noe, <a href="/A112141/b112141.txt">Table of n, a(n) for n = 1..100</a>

%F a(n) = Product_{i=1..n} A001358(i).

%F A001222(a(n)) = 2*n.

%e a(10) = 4*6*9*10*14*15*21*22*25*26 = 136216080000, the product of the first 10 semiprimes.

%e From _Gus Wiseman_, Dec 06 2020: (Start)

%e The sequence of terms together with their prime signatures begins:

%e 4: (2)

%e 24: (3,1)

%e 216: (3,3)

%e 2160: (4,3,1)

%e 30240: (5,3,1,1)

%e 453600: (5,4,2,1)

%e 9525600: (5,5,2,2)

%e 209563200: (6,5,2,2,1)

%e 5239080000: (6,5,4,2,1)

%e 136216080000: (7,5,4,2,1,1)

%e 4495130640000: (7,6,4,2,2,1)

%e 152834441760000: (8,6,4,2,2,1,1)

%e 5349205461600000: (8,6,5,3,2,1,1)

%e 203269807540800000: (9,6,5,3,2,1,1,1)

%e 7927522494091200000: (9,7,5,3,2,2,1,1)

%e 364666034728195200000: (10,7,5,3,2,2,1,1,1)

%e 17868635701681564800000: (10,7,5,5,2,2,1,1,1)

%e (End)

%p A112141 := proc(n)

%p mul(A001358(i),i=1..n) ;

%p end proc:

%p seq(A112141(n),n=1..10) ; # _R. J. Mathar_, Jun 30 2020

%t NextSemiPrime[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sp = n + sgn; While[c < Abs[k], While[ PrimeOmega[sp] != 2, If[sgn < 0, sp--, sp++]]; If[sgn < 0, sp--, sp++]; c++]; sp + If[sgn < 0, 1, -1]]; f[n_] := Times @@ NestList[ NextSemiPrime@# &, 2^2, n - 1]; Array[f, 18] (* _Robert G. Wilson v_, Jun 13 2013 *)

%t FoldList[Times,Select[Range[30],PrimeOmega[#]==2&]] (* _Gus Wiseman_, Dec 06 2020 *)

%o (PARI) a(n)=my(v=vector(n),i,k=3);while(i<n,if(bigomega(k++)==2,v[i++]=k)); prod(i=1,n,v[i]) \\ _Charles R Greathouse IV_, Apr 04 2013

%o (Python)

%o from sympy import factorint

%o def aupton(terms):

%o alst, k, p = [], 1, 1

%o while len(alst) < terms:

%o if sum(factorint(k).values()) == 2:

%o p *= k

%o alst.append(p)

%o k += 1

%o return alst

%o print(aupton(18)) # _Michael S. Branicky_, Aug 31 2021

%Y Partial sums of semiprimes are A062198.

%Y First differences of semiprimes are A065516.

%Y A000040 lists primes, with partial products A002110 (primorials).

%Y A000142 lists factorials, with partial products A000178 (superfactorials).

%Y A001358 lists semiprimes, with partial products A112141 (this sequence).

%Y A005117 lists squarefree numbers, with partial products A111059.

%Y A006881 lists squarefree semiprimes, with partial products A339191.

%Y A101048 counts partitions into semiprimes (restricted: A338902).

%Y A320655 counts factorizations into semiprimes.

%Y A338898/A338912/A338913 give the prime indices of semiprimes, with product/sum/difference A087794/A176504/A176506.

%Y A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.

%Y Cf. A001222, A084126, A084127, A115392, A168472, A320732, A320892.

%K easy,nonn

%O 1,1

%A _Jonathan Vos Post_, Nov 28 2005