%I #10 Mar 15 2023 12:10:28
%S 61,181,269,277,349,389,509,541,613,653,661,701,757,773,821,877,941,
%T 1069,1109,1117,1181,1237,1301,1597,1613,1693,1733,1741,1789,1877,
%U 1933,1997,2053,2069,2141,2221,2269,2357,2389,2437,2557,2741,2789,2837,2909
%N Primes of the form (4n+1)^2 + (4m+2)^2 where n=1,2,3.., m=1,2,3..
%H Vincenzo Librandi, <a href="/A087870/b087870.txt">Table of n, a(n) for n = 1..10000</a>
%t Union[Select[Flatten[Table[(4n+1)^2+(4m+2)^2, {n,100}, {m,100}]], PrimeQ]] (* _T. D. Noe_, Nov 02 2006 *)
%o (PARI) fourmp1a(n) = { for(y=1,n, for(x=1,n, z=(4*y+1)^2+(4*x+2)^2; if(isprime(z),print1(z","); break) ); ) }
%K nonn
%O 1,1
%A _Cino Hilliard_, Oct 10 2003
%E Corrected by _T. D. Noe_, Nov 02 2006
|