A221265 - OEIS

%I #11 Jan 26 2022 11:58:01

%S 5,10,13,17,20,25,26,29,34,37,40,41,50,52,53,58,61,65,68,73,74,80,82,

%T 85,89,97,100,101,104,106,109,113,116,122,125,130,136,137,145,146,148,

%U 149,157,160,164,169,170,173,178,181,185,193,194,195,197,200,202

%N Numbers having more distinct prime factors of form 4*k+1 than of 4*k+3.

%C A005089(a(n)) > A005091(a(n)); A005094(a(n)) > 0.

%H Reinhard Zumkeller, <a href="/A221265/b221265.txt">Table of n, a(n) for n = 1..10000</a>

%t mdpf1Q[n_]:=Module[{f=Transpose[FactorInteger[n]][[1]]},Count[f,_?(Mod[ #,4] == 1&)]>Count[f,_?(Mod[#,4]==3&)]]; Select[Range[2,250],mdpf1Q] (* _Harvey P. Dale_, Mar 03 2016 *)

%o (Haskell)

%o a221265 n = a221265_list !! (n-1)

%o a221265_list = filter ((> 0) . a005094) [1..]

%Y Cf. A078613, A221264.

%Y Cf. A005089, A005091, A005094.

%K nonn

%O 1,1

%A _Reinhard Zumkeller_, Jan 07 2013