OFFSET
1,1
COMMENTS
Sequence naturally partitions into two sequences: all primes p with ord_p(-2) odd (A163183, the primes dividing 2^j +1 for some odd j) and certain primes p with ord_p(-2) even (A163185). - Christopher J. Smyth, Jul 23 2009
REFERENCES
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
FORMULA
MATHEMATICA
QuadPrimes2[1, 0, 2, 10000] (* see A106856 *)
Select[Prime[Range[200]], MemberQ[{1, 2, 3}, Mod[#, 8]]&] (* Harvey P. Dale, Mar 16 2013 *)
PROG
(Haskell)
a033203 n = a033203_list !! (n-1)
a033203_list = filter ((== 1) . a010051) a047476_list
-- Reinhard Zumkeller, Dec 29 2012, Jan 22 2012
(Magma) [p: p in PrimesUpTo(600) | p mod 8 in [1..3]]; // Vincenzo Librandi, Aug 11 2012
(Magma) [p: p in PrimesUpTo(800) | NormEquation(2, p) eq true]; // Bruno Berselli, Jul 03 2016
(PARI) is(n)=isprime(n) && issquare(Mod(-2, n)) \\ Charles R Greathouse IV, Nov 29 2016
CROSSREFS
KEYWORD
nonn,nice,easy
AUTHOR
STATUS
approved