A107213 Primes of the form x^2 + 96y^2.
97, 409, 433, 457, 673, 937, 1009, 1033, 1153, 1321, 1489, 1609, 1657, 1753, 1777, 2089, 2161, 2377, 2521, 2593, 2689, 2713, 3121, 3217, 3361, 3457, 3673, 3769, 3889, 4297, 4561, 4729, 4801, 4993, 5209, 5233, 5857, 5881, 6073, 6121, 6337
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, 96, 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)\96), if(isprime(t=w+96*y^2), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Feb 10 2017
Comments