|
%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
|
