A141166 Primes of the form x^2+15*x*y-y^2.
37, 53, 173, 193, 229, 241, 347, 359, 383, 439, 443, 449, 461, 503, 509, 541, 593, 607, 617, 619, 643, 691, 907, 967, 977, 1019, 1051, 1063, 1097, 1109, 1249, 1277, 1291, 1303, 1321, 1399, 1429, 1583, 1667, 1741, 1783, 1993, 1997, 2003, 2087, 2137, 2143, 2333, 2347, 2351
Offset: 1
Keywords
Examples
a(2)=53 because we can write 53= 3^2+15*3*1-1^2
References
- Z. I. Borevich and I. R. Shafarevich, Number Theory
Links
- Juan Arias-de-Reyna, Table of n, a(n) for n = 1..10000
- N. J. A. Sloane et al., Binary Quadratic Forms and OEIS: Index to related sequences, programs, references. OEIS wiki, June 2014.
- D. B. Zagier, Zetafunktionen und quadratische Körper, Springer, 1981.
Crossrefs
Programs
-
Mathematica
lim = 100; Rest@ Union@ Abs@ Flatten@ Table[x^2 + 15 x y - y^2, {x, lim}, {y, lim}] /. n_ /; CompositeQ@ n -> Nothing (* Michael De Vlieger, Jan 27 2016 *)
-
PARI
is_A141166(p)=qfbsolve(Qfb(1,15,-1),p) \\ Returns nonzero (actually, a solution [x,y]) iff p is a member of the sequence. For efficiency it is assumed that p is prime. Example usage: select(is_A141166,primes(500)) - M. F. Hasler, Jan 27 2016
Comments