|
%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
|
