A087880 - OEIS

%I #9 Aug 26 2018 20:15:12

%S 13,53,109,149,157,229,317,373,397,421,461,557,709,733,797,829,853,

%T 997,1013,1021,1061,1213,1229,1277,1381,1429,1453,1493,1549,1621,1637,

%U 1669,1709,1861,1901,1949,1973,2213,2237,2293,2309,2333,2341,2381,2477,2549

%N Primes of the form (4n+3)^2 + (4m+2)^2, m,n = 0,1,2..

%H Robert Israel, <a href="/A087880/b087880.txt">Table of n, a(n) for n = 1..10000</a>

%p N:= 10000:

%p A:= NULL:

%p for x from 3 by 4 while x^2 < N do

%p   for y from 2 by 4 while x^2 + y^2 < N do

%p     v:= x^2 + y^2;

%p     if isprime(v) then A:= A,v fi

%p od od:

%p sort(convert({A},list)); # _Robert Israel_, Aug 26 2018

%t With[{nn=20},Select[Flatten[Table[(4n+3)^2+(4m+2)^2,{n,0,nn},{m,0,nn}],1],PrimeQ]]//Union (* _Harvey P. Dale_, Mar 04 2018 *)

%K nonn

%O 1,1

%A _Cino Hilliard_, Oct 11 2003