%I #20 Sep 08 2022 08:45:33
%S 181,193,241,277,373,433,577,661,709,733,757,829,853,877,997,1069,
%T 1201,1237,1549,1597,1609,1621,1657,1669,1693,1777,1789,1933,1993,
%U 2113,2269,2293,2377,2389,2557,2617,2677,2749,2833,2857,2917,2953
%N Primes of the form x^2 + 177*y^2.
%C Discriminant = -708.
%C The primes are congruent to {1, 25, 49, 85, 121, 133, 145, 169, 181, 193, 205, 241, 253, 265, 277, 289, 361, 373, 433, 481, 493, 517, 529, 553, 577, 625, 661, 685, 697} (mod 708).
%H Vincenzo Librandi and Ray Chandler, <a href="/A139649/b139649.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, 177, 10000] (* see A106856 *)
%o (Magma) [ p: p in PrimesUpTo(4000) | p mod 708 in {1, 25, 49, 85, 121, 133, 145, 169, 181, 193, 205, 241, 253, 265, 277, 289, 361, 373, 433, 481, 493, 517, 529, 553, 577, 625, 661, 685, 697}]; // _Vincenzo Librandi_, Jul 28 2012
%o (Magma) k:=177; [p: p in PrimesUpTo(3000) | NormEquation(k, p) eq true]; // _Bruno Berselli_, Jun 01 2016
%K nonn,easy
%O 1,1
%A _T. D. Noe_, Apr 29 2008
|