A141181 Primes of the form 2*x^2+3*x*y-4*y^2 (as well as of the form 2*x^2+7*x*y+y^2).
2, 5, 23, 31, 37, 41, 43, 59, 61, 73, 83, 103, 107, 113, 127, 131, 139, 163, 173, 197, 223, 241, 251, 269, 271, 277, 283, 307, 337, 349, 353, 359, 367, 373, 379, 389, 401, 409, 419, 431, 433, 443, 449, 461, 467, 487, 491, 523, 541, 569, 599, 607, 613, 617, 619
Offset: 1
Keywords
Examples
a(3) = 23 because we can write 23 = 2*3^2+3*3*1-4*1^2 (or 23 = 2*2^2+7*2*1+1^2).
References
- Z. I. Borevich and I. R. Shafarevich, Number Theory. Academic Press, NY, 1966.
- D. B. Zagier, Zetafunktionen und quadratische Körper, Springer, 1981.
Links
- N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
Crossrefs
Cf. A038872 (d=5). A038873 (d=8). A068228, A141123 (d=12). A038883 (d=13). A038889 (d=17). A141111, A141112 (d=65).
Primes in A035269.
A subsequence of (and may possibly coincide with) A038919. - R. J. Mathar, Jul 22 2008
For a list of sequences giving numbers and/or primes represented by binary quadratic forms, see the "Binary Quadratic Forms and OEIS" link.
Programs
-
Mathematica
Reap[For[p = 2, p < 1000, p = NextPrime[p], If[FindInstance[p == 2*x^2 + 3*x*y - 4*y^2, {x, y}, Integers, 1] =!= {}, Print[p]; Sow[p]]]][[2, 1]] (* Jean-François Alcover, Oct 25 2016 *)
-
PARI
select(p->isprime(p)&&qfbsolve(Qfb(1,7,2),p),[1..1500]) \\ This is to provide a generic characteristic function ("is_A141181") as 1st arg of select(), there are other ways to produce the sequence more efficiently. - M. F. Hasler, Jan 15 2016
Comments