|
%I #26 Mar 07 2024 01:30:17
%S 1800,2700,3528,4500,5292,8712,9800,12168,12348,13068,18252,20808,
%T 24200,24500,25992,31212,33075,33800,34300,38088,38988,47432,47916,
%U 55125,57132,57800,60500,60552,66248,69192,72200,77175,79092,81675,84500
%N Numbers of the form p^3*q^2*r^2 where p, q, and r are distinct primes.
%C Subsequence of A225228. - _Reinhard Zumkeller_, May 03 2013
%H T. D. Noe, <a href="/A179695/b179695.txt">Table of n, a(n) for n = 1..1000</a>
%H Will Nicholes, <a href="https://willnicholes.com/2010/06/06/list-of-prime-signatures">List of prime signatures</a>, 2010.
%H <a href="/index/Pri#prime_signature">Index to sequences related to prime signature</a>.
%F A050326(a(n)) = 5. - _Reinhard Zumkeller_, May 03 2013
%F Sum_{n>=1} 1/a(n) = P(2)^2*P(3)/2 - P(3)*P(4)/2 - P(2)*P(5) + P(7) = 0.0032578591481263202818..., where P is the prime zeta function. - _Amiram Eldar_, Mar 07 2024
%t f[n_]:=Sort[Last/@FactorInteger[n]]=={2,2,3}; Select[Range[10^5], f]
%t f[n_]:={Times@@(n^{2,2,3}),Times@@(n^{2,3,2}),Times@@(n^{3,2,2})}; Module[ {nn=20},Select[Flatten[f/@Subsets[Prime[Range[nn]],{3}]],#<= 72*Prime[ nn]^2&]]//Union (* _Harvey P. Dale_, Jul 05 2019 *)
%o (PARI) list(lim)=my(v=List(),t1,t2);forprime(p=2, (lim\36)^(1/3), t1=p^3;forprime(q=2, sqrt(lim\t1), if(p==q, next);t2=t1*q^2;forprime(r=q+1, sqrt(lim\t2), if(p==r,next);listput(v,t2*r^2)))); vecsort(Vec(v)) \\ _Charles R Greathouse IV_, Jul 24 2011
%Y Cf. A050326, A225228.
%Y Cf. A085548, A085541, A085964, A085965, A085967.
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Jul 24 2010
|
