A107203 Primes of the form 9x^2 + 10y^2.

19, 241, 331, 499, 571, 601, 691, 739, 769, 1009, 1129, 1249, 1291, 1531, 1579, 1699, 2011, 2089, 2131, 2161, 2521, 2689, 2731, 2851, 2971, 3001, 3049, 3259, 3331, 3499, 3691, 3739, 3889, 4051, 4129, 4219, 4339, 4441, 4561, 4729, 4801

Offset

1,1

Comments

Discriminant = -360. 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[9, 0, 10, 10000] (* see A106856 *)

Prog

(PARI) list(lim)=my(v=List(), w, t); for(x=1, sqrtint(lim\9), w=9*x^2; for(y=1, sqrtint((lim-w)\10), if(isprime(t=w+10*y^2), listput(v, t)))); Set(v) \\ Charles R Greathouse IV, Feb 10 2017

Crossrefs

Sequence in context: A021814 A299864 A025933 * A021544 A021772 A019783

Adjacent sequences: A107200 A107201 A107202 * A107204 A107205 A107206

Keyword

nonn,easy

Author

T. D. Noe, May 13 2005

Status

approved