A079733 Primes of the form x^2 + y^2 + 3 (x,y nonnegative).
3, 5, 7, 11, 13, 19, 23, 29, 37, 43, 53, 61, 67, 71, 83, 101, 103, 107, 109, 131, 139, 149, 151, 163, 167, 173, 181, 197, 199, 211, 229, 263, 277, 293, 317, 331, 349, 359, 373, 389, 397, 419, 439, 461, 467, 487, 491, 503, 523, 541, 547, 557, 587, 599, 613, 631, 643, 653
Offset: 1
Keywords
Links
- Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Programs
-
Mathematica
f[upto_]:=Module[{max=Ceiling[Sqrt[upto]]},Select[Select[Union[Total[#]+3&/@(Tuples[Range[0,max],{2}]^2)],PrimeQ],#<=upto&]];f[1200] (* Vincenzo Librandi, Jul 15 2012 *)