|
%I #41 Dec 21 2020 02:20:27
%S 6,10,12,14,20,22,24,26,28,34,38,40,44,46,48,52,56,58,62,68,74,76,80,
%T 82,86,88,92,94,96,104,106,112,116,118,122,124,134,136,142,146,148,
%U 152,158,160,164,166,172,176,178,184,188,192,194,202,206,208,212,214,218,224
%N Numbers of the form 2^k * p where k > 0 and p is an odd prime.
%C Even numbers with 2 distinct prime factors where the odd factor is prime.
%C A proper subset of A098202. E.g., 210 is not here, but it is there. Also differs from A100367: 36, 100, 108, 196, etc. are missing here. Different also from A036348 because 90 and 180 are not here.
%C A128691 is a subsequence; A078834(a(n)) = A006530(a(n)). - _Reinhard Zumkeller_, Sep 19 2011
%C Composite numbers k having the property that the number of divisors of 2k equals the number of divisors of k + 2. All primes satisfy this property. - _Gary Detlefs_, Jan 23 2019
%H Reinhard Zumkeller, <a href="/A100368/b100368.txt">Table of n, a(n) for n = 1..10000</a>
%F Numbers of the form 2^k*p where k > 0, p is an odd prime.
%F a(n) = 2*A038550(n). - _Amiram Eldar_, Dec 21 2020
%p N:= 1000: # to get all terms <= N
%p P:= select(isprime, [seq(i,i=3..N/2,2)]):
%p S:= {seq(seq(2^i*p,i=1..ilog2(N/p)),p=P)}:
%p sort(convert(S,list)); # _Robert Israel_, Jul 09 2017
%p with(numtheory): for n from 1 to 224 do if tau(2*n)=tau(n)+2 and not isprime(n) then print(n) fi od # _Gary Detlefs_, Jan 22 2019
%t <<NumberTheory`NumberTheoryFunctions` p2[x_] :=Part[PrimeFactorList[x], 2]; lf[x_] :=Length[FactorInteger[x]];ta={{0}}; Do[If[Equal[lf[n], 2]&&EvenQ[n]&&IntegerQ[Log[2, n/p2[n]]], ta=Append[ta, n];Print[n]], {n, 1, 256}];ta=Delete[ta, 1]
%t dpf2Q[n_]:=Module[{f=FactorInteger[n]},Length[f]==2&&f[[1,1]]==2&&f[[-1,2]]==1]; Select[Range[2,250,2],dpf2Q] (* _Harvey P. Dale_, Sep 03 2016 *)
%o (Haskell)
%o import Data.Set (singleton, deleteFindMin, insert)
%o a100368 n = a100368_list !! (n-1)
%o a100368_list = f (singleton 6) (tail a065091_list) where
%o f s ps'@(p:ps) | mod m 4 > 0 = m : f (insert (2*p) $ insert (2*m) s') ps
%o | otherwise = m : f (insert (2*m) s') ps'
%o where (m,s') = deleteFindMin s
%o -- _Reinhard Zumkeller_, Sep 19 2011
%o (PARI) is(n)=n%2==0 && isprime(n>>valuation(n,2)) \\ _Charles R Greathouse IV_, Jul 09 2017
%o (PARI) list(lim)=my(v=List()); for(k=1,logint(lim\3,2), forprime(p=3,lim>>k, listput(v,p<<k))); Set(v) \\ _Charles R Greathouse IV_, Jul 09 2017
%o (GAP) a:=Filtered([1..224],n->Tau(2*n)=Tau(n)+2 and not IsPrime(n));; Print(a); # _Muniru A Asiru_, Jan 22 2019
%Y Cf. A001221, A038550, A098902, A100367, A036348, A100484.
%K nonn,easy
%O 1,1
%A _Labos Elemer_, Nov 22 2004
%E Name edited by _Charles R Greathouse IV_, Jul 09 2017
|
