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

13, 53, 109, 149, 157, 229, 317, 373, 397, 421, 461, 557, 709, 733, 797, 829, 853, 997, 1013, 1021, 1061, 1213, 1229, 1277, 1381, 1429, 1453, 1493, 1549, 1621, 1637, 1669, 1709, 1861, 1901, 1949, 1973, 2213, 2237, 2293, 2309, 2333, 2341, 2381, 2477, 2549

Offset

1,1

Links

Robert Israel, Table of n, a(n) for n = 1..10000

Maple

N:= 10000:
A:= NULL:
for x from 3 by 4 while x^2 < N do
  for y from 2 by 4 while x^2 + y^2 < N do
    v:= x^2 + y^2;
    if isprime(v) then A:= A, v fi
od od:
sort(convert({A}, list)); # Robert Israel, Aug 26 2018

Mathematica

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 *)

Crossrefs

Sequence in context: A307749 A031905 A214523 * A228872 A041322 A022284

Adjacent sequences: A087877 A087878 A087879 * A087881 A087882 A087883

Keyword

nonn

Author

Cino Hilliard, Oct 11 2003

Status

approved