A107203 Primes of the form 9x^2 + 10y^2.
19, 241, 331, 499, 571, 601, 691, 739, 769, 1009, 1129, 1249, 1291, 1531, 1579, 1699, 2011, 2089, 2131, 2161, 2521, 2689, 2731, 2851, 2971, 3001, 3049, 3259, 3331, 3499, 3691, 3739, 3889, 4051, 4129, 4219, 4339, 4441, 4561, 4729, 4801
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[9, 0, 10, 10000] (* see A106856 *)
-
PARI
list(lim)=my(v=List(),w,t); for(x=1, sqrtint(lim\9), w=9*x^2; for(y=1, sqrtint((lim-w)\10), if(isprime(t=w+10*y^2), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Feb 10 2017
Comments