%I A107145
%S A107145 41,89,241,281,401,409,449,521,569,601,641,761,769,809,881,929,1009,
%T A107145 1049,1129,1201,1249,1289,1321,1361,1409,1481,1489,1601,1609,1721,1801,
%U A107145 1889,2081,2089,2129,2161,2281,2441,2521,2609,2689,2729,2801
%N A107145 Primes of the form x^2+40y^2.
%C A107145 Discriminant=-160. See A107132 for more information.
%F A107145 The primes are congruent to {1, 9} (mod 40). - T. D. Noe (noe(AT)sspectra.com), 
               Apr 29 2008
%t A107145 f[x_,y_]:=x^2+40*y^2; lst={};Do[Do[p=f[x,y];If[PrimeQ[p]&&p>0,AppendTo[lst,
               p]],{y,-4!,3*4!}],{x,-4!,3*4!}];Take[Union[lst],90] [From Vladimir 
               Orlovsky (4vladimir(AT)gmail.com), Jun 30 2009]
%t A107145 QuadPrimes[1, 0, 40, 10000] (* see A106856 *)
%Y A107145 Cf. A139643.
%Y A107145 Sequence in context: A142411 A139924 A155572 this_sequence A087857 A139995 
               A044179
%Y A107145 Adjacent sequences: A107142 A107143 A107144 this_sequence A107146 A107147 
               A107148
%K A107145 nonn,easy
%O A107145 1,1
%A A107145 T. D. Noe (noe(AT)sspectra.com), May 13 2005