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,1

Comments

Discriminant = -384. See A107132 for more information.

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)

Mathematica

QuadPrimes2[1, 0, 96, 10000] (* see A106856 *)

Prog

(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

Crossrefs

Sequence in context: A174338 A050666 A160440 * A157331 A142834 A361678

Adjacent sequences: A107210 A107211 A107212 * A107214 A107215 A107216

Keyword

nonn,easy

Author

T. D. Noe, May 13 2005

Status

approved