A107206 Primes of the form x^2 + 90y^2.
139, 211, 379, 409, 619, 811, 859, 1051, 1171, 1201, 1321, 1459, 1489, 1609, 1801, 2179, 2251, 2281, 2371, 2539, 2659, 3019, 3121, 3169, 3361, 3529, 3571, 3769, 3931, 4099, 4201, 4651, 5059, 5419, 5449, 5641, 5689, 5779, 5851, 5881, 6091
Offset: 1
Links
- Vincenzo Librandi and Ray Chandler, Table of n, a(n) for n = 1..10000 [First 1000 terms from Vincenzo Librandi]
- N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
Programs
-
Mathematica
QuadPrimes2[1, 0, 90, 10000] (* see A106856 *)
-
PARI
list(lim)=my(v=List(),w,t); for(x=1, sqrtint(lim\1), w=x^2; for(y=1, sqrtint((lim-w)\90), if(isprime(t=w+90*y^2), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Feb 10 2017
Comments