%I #9 Dec 22 2014 02:03:28

%S 7,19,31,43,47,67,103,127,139,151,199,223,227,263,271,283,307,367,379,

%T 463,479,487,523,547,571,587,607,619,631,643,659,691,727,787,811,823,

%U 859,883,907,911,967,983,1019,1039,1051,1063,1231,1279,1291,1303,1307

%N Primes of the form x^2 + y^2 + 2 (x,y positive).

%p S := {}: for x from 1 to 100 do for y from 1 to 100 do S := S union {x^2+y^2+2} od:od:S := sort(convert(S, list)): for i from 1 to 500 do if isprime(S[i]) then printf(`%d,`,S[i]) fi:od: # _James A. Sellers_, Feb 25 2003

%Y Cf. A079545, A079739.

%K nonn,easy

%O 1,1

%A _N. J. A. Sloane_, Feb 18 2003

%E More terms from _James A. Sellers_, Feb 25 2003