(MAGMAMagma) [5] cat [ p: p in PrimesUpTo(3000) | p mod 40 in {13, 37} ]; // Vincenzo Librandi, Jul 24 2012
Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
(MAGMAMagma) [5] cat [ p: p in PrimesUpTo(3000) | p mod 40 in {13, 37} ]; // Vincenzo Librandi, Jul 24 2012
editing
approved
Primes of the form 5x^2 + 8y^2.
(PARI) list(lim)=my(v=List([5]), t); forprime(p=13, lim, t=p%40; if(t==13||t==37, listput(v, p))); Vec(v) \\ Charles R Greathouse IV, Feb 09 2017
approved
editing
Vincenzo Librandi and Ray Chandler, <a href="/A107144/b107144_1.txt">Table of n, a(n) for n = 1..10000</a> [First 1000 terms from Vincenzo Librandi]
editing
approved
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)
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)
QuadPrimesQuadPrimes2[5, 0, 8, 10000] (* see A106856 *)
QuadPrimes[5, 0, 8, 10000] (* see A106856 *) Note: the original QuadPrimes had a bug which could sometimes give wrong answers. This sequence should be checked (unless the coefficient of xy in the quadratic form is zero, in which case QuadPrimes gives correct answers). - N. J. A. Sloane, Jun 04 2014
QuadPrimes[5, 0, 8, 10000] (* see A106856 *)
QuadPrimes[5, 0, 8, 10000] (* see A106856 *) Note: the original QuadPrimes had a bug which could sometimes give wrong answers. This sequence should be checked (unless the coefficient of xy in the quadratic form is zero, in which case QuadPrimes gives correct answers). - N. J. A. Sloane, Jun 04 2014