A079739 - OEIS

%I #17 Dec 22 2014 02:06:36

%S 2,3,7,11,19,31,43,47,67,83,103,127,139,151,199,223,227,263,271,283,

%T 307,367,379,443,463,479,487,523,547,571,587,607,619,631,643,659,691,

%U 727,787,811,823,859,883,907,911,967,983,1019,1039,1051,1063,1091,1231

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

%H Vincenzo Librandi, <a href="/A079739/b079739.txt">Table of n, a(n) for n = 1..1000</a>

%p S := {}: for x from 0 to 100 do for y from 0 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

%t f[upto_]:=Module[{max=Ceiling[Sqrt[upto]]},Select[Select[ Union[Total[#]+2&/@(Tuples[Range[0,max],{2}]^2)], PrimeQ], #<=upto&]]; f[1250] (* _Harvey P. Dale_, Mar 19 2011 *)

%Y Cf. A079545, A079740.

%K nonn,easy

%O 1,1

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

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