%I #20 Sep 08 2022 08:45:33
%S 181,229,421,661,709,829,1021,1321,1489,1549,1609,1621,1741,2029,2161,
%T 2269,2281,2341,2689,3001,3061,3169,3301,3469,3529,4129,4261,4621,
%U 4789,4801,4909,5281,5449,5569,5581,5641,5701,6121,6229,6301,6361
%N Primes of the form x^2 + 165*y^2.
%C Discriminant = -660.
%C The primes are congruent to {1, 49, 169, 181, 229, 289, 301, 361, 421, 529} (mod 660).
%H Vincenzo Librandi and Ray Chandler, <a href="/A139648/b139648.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from Vincenzo Librandi).
%H N. J. A. Sloane et al., <a href="https://oeis.org/wiki/Binary_Quadratic_Forms_and_OEIS">Binary Quadratic Forms and OEIS</a> (Index to related sequences, programs, references)
%t QuadPrimes2[1, 0, 165, 10000] (* see A106856 *)
%o (Magma) [ p: p in PrimesUpTo(7000) | p mod 660 in {1, 49, 169, 181, 229, 289, 301, 361, 421, 529}]; // _Vincenzo Librandi_, Jul 28 2012
%o (Magma) k:=165; [p: p in PrimesUpTo(7000) | NormEquation(k, p) eq true]; // _Bruno Berselli_, Jun 01 2016
%K nonn,easy
%O 1,1
%A _T. D. Noe_, Apr 29 2008
|