|
%I #32 Aug 02 2015 16:44:47
%S 60,84,90,120,126,132,140,150,156,168,180,198,204,220,228,234,240,252,
%T 260,264,270,276,280,294,300,306,308,312,315,336,340,342,348,350,360,
%U 364,372,378,380,396,408,414,440,444,450,456,460,468,476,480,490,492
%N Nonsquarefree numbers divisible by exactly three distinct primes.
%H Reinhard Zumkeller, <a href="/A178212/b178212.txt">Table of n, a(n) for n = 1..10000</a>
%F A001221(a(n)) = 3; A001222(a(n)) > 3; A000005(n) >= 12;
%F a(n) = A123712(n) for n <= 52, possibly more.
%e 60 is in the sequence because it is not squarefree and it is divisible by three distinct primes: 2, 3, 5.
%e 72 is not in the sequence, because although it is not squarefree, it is divisible by only two distinct primes: 2 and 3.
%t nsD3Q[n_] := Block[{fi = FactorInteger@ n}, Length@ fi == 3 && Plus @@ Last /@ fi > 3]; Select[ Range@ 494, nsD3Q] (* _Robert G. Wilson v_, Feb 09 2012 *)
%t Select[Range[500], PrimeNu[#] == 3 && PrimeOmega[#] > 3 &] (* _Alonso del Arte_, Mar 23 2015, based on a comment from _Robert G. Wilson v_, Feb 09 2012; requires Mathematica 7.0+ *)
%o (PARI) is_A178212(n)={ omega(n)==3 & bigomega(n)>3 }
%o for(n=1,999,is_A178212(n) & print1(n",")) \\ _M. F. Hasler_, Feb 09 2012
%o (Haskell)
%o a178212 n = a178212_list !! (n-1)
%o a178212_list = filter f [1..] where
%o f x = length (a027748_row x) == 3 && any (> 1) (a124010_row x)
%o -- _Reinhard Zumkeller_, Apr 03 2015
%Y A subsequence of A033987.
%Y A085987 is a subsequence.
%Y Cf. A007304, A200511, A200521.
%Y Cf. A000005, A001221, A001222, A013929, A085987, A123712.
%Y Cf. A027748, A124010.
%K nonn
%O 1,1
%A _Reinhard Zumkeller_, May 24 2010
|
