cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A139827 Primes of the form 2x^2 + 2xy + 17y^2.

Original entry on oeis.org

2, 17, 29, 41, 101, 149, 173, 197, 233, 281, 293, 461, 557, 569, 593, 677, 701, 761, 809, 821, 857, 941, 953, 1097, 1217, 1229, 1289, 1361, 1481, 1493, 1553, 1601, 1613, 1733, 1877, 1889, 1913, 1949, 1997, 2081, 2129, 2141, 2153, 2213, 2273, 2309, 2393, 2417
Offset: 1

Views

Author

T. D. Noe, May 02 2008, May 07 2008

Keywords

Comments

Discriminant = -132.
Consider the quadratic form f(x,y) = ax^2 + bxy + cy^2. When the discriminant d=b^2-4ac is -4 times an idoneal number (A000926), there is exactly one class for each genus. As a result, the primes generated by f(x,y) are the same as the primes congruent to S (mod -d), where S is a set of numbers less than -d. The table on page 60 of Cox shows that there are exactly 331 quadratic forms having this property. The 217 sequences starting with this one complete the collection in the OEIS.
When a=1 and b=0, f(x,y) is a quadratic form whose congruences are discussed in A139642. Let N be an idoneal number. Then there are 2^r reduced quadratic forms whose discriminant is -4N, where r=1,2,3, or 4. By collecting the residuals p (mod 4N) for primes p generated by the i-th reduced quadratic form, we can empirically find a set Si. To show that the 2^r sets Si are complete, we only need to show that the union of the Si is equal to the set of numbers k such that the Jacobi symbol (-k/4N)=1.

References

  • David A. Cox, Primes of the Form x^2 + n y^2, Wiley, 1989.

Crossrefs

Cf. also A139653, A139904-A139906 (d=-1012), A139654, A139907-A139913 (d=-1092), A139655, A139914-A139920 (d=-1120), A139656, A139921-A139927 (d=-1248), A139657, A139928-A139934 (d=-1320), A139658, A139935-A139941 (d=-1380), A139659, A139942-A139948 (d=-1428), A139660, A139949-A139955 (d=-1540), A139661, A139956-A139962 (d=-1632), A139662, A139963-A139969 (d=-1848), A139663, A139970-A139976 (d=-2080), A139664, A139977-A139983 (d=-3040), A139665, A139984-A139998 (d=-3360), A139666, A139999-A140013 (d=-5280), A139667, A140014-A140028 (d=-5460), A139668, A140029-A140043 (d=-7392).
For a more complete list of sequences giving numbers and/or primes represented by binary quadratic forms, see the "Binary Quadratic Forms and OEIS" link.

Programs

  • Magma
    [ p: p in PrimesUpTo(2500) | p mod 132 in {2, 17, 29, 41, 65, 101}]; // Vincenzo Librandi, Jul 29 2012
    
  • Mathematica
    QuadPrimes2[2, -2, 17, 2500] (* see A106856 *)
    t = Table[{2, 17, 29, 41, 65, 101} + 132*n, {n, 0, 50}]; Select[Flatten[t], PrimeQ] (* T. D. Noe, Jun 21 2012 *)
  • PARI
    v=[2, 17, 29, 41, 65, 101]; select(p->setsearch(v,p%132),primes(100)) \\ Charles R Greathouse IV, Jan 08 2013

Formula

The primes are congruent to {2, 17, 29, 41, 65, 101} (mod 132).

A000926 Euler's "numerus idoneus" (or "numeri idonei", or idoneal, or suitable, or convenient numbers).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 15, 16, 18, 21, 22, 24, 25, 28, 30, 33, 37, 40, 42, 45, 48, 57, 58, 60, 70, 72, 78, 85, 88, 93, 102, 105, 112, 120, 130, 133, 165, 168, 177, 190, 210, 232, 240, 253, 273, 280, 312, 330, 345, 357, 385, 408, 462, 520, 760, 840, 1320, 1365, 1848
Offset: 1

Views

Author

Keywords

Comments

There are many equivalent definitions of these numbers. Based on Cox, Theorem 3.22 and Proposition 3.24 and a comment by Eric Rains (rains(AT)caltech.edu), we can say that a positive number n belongs to this sequence if and only if any of the following equivalent statements is true:
(1) Let m > 1 be an odd number relatively prime to n which can be written in the form x^2 + n*y^2 with x, y relatively prime. If the equation m = x^2 + n*y^2 has only one solution with x, y >= 0, then m is a prime number. [Euler]
(2) Every genus of quadratic forms of discriminant -4n consists of a single class. [Gauss]
(3) If a*x^2 + b*x*y + c*y^2 is a reduced quadratic form of discriminant -4n, then either b=0, a=b or a = c. [Cox]
(4) Two quadratic forms of discriminant -4n are equivalent if and only if they are properly equivalent. [Cox]
(5) The class group C(-4n) is isomorphic to (Z/2Z)^m for some integer m. [Cox]
(6) n is not of the form ab+ac+bc with 0 < a < b < c. (See proof in link below.) [Rains]
It is conjectured that the list given here is complete. Chowla showed that the list is finite and Weinberger showed that there is at most one further term.
If an additional term exists it is > 100000000. - Jud McCranie, Jun 27 2005
The terms shown are the union of {1,2,3,4,7}, A033266, A033267, A033268 and A033269 (corresponding to class numbers 1, 2, 4, 8 and 16 respectively).
Note that for n in this sequence, n+1 is either a prime, twice a prime, the square of a prime, 8 or 16. - T. D. Noe, Apr 08 2004. [This is a general theorem that is not hard to prove using genus theory. The "32" in the original comment was an error. - Tom Hagedorn (hagedorn(AT)tcnj.edu), Dec 29 2008]
Also numbers n such that for all primes p such that p is a quadratic residue (mod 4*n) and p-n is a quadratic residue (mod 4*n), p can be uniquely written into the form as x^2+n*y^2. - V. Raman, Nov 25 2013

References

  • Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See Table 97 at p. 272.
  • Z. I. Borevich and I. R. Shafarevich, Number Theory. Academic Press, NY, 1966, pp. 425-430.
  • David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, Section 3.
  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 1848, p. 146, Ellipses, Paris 2008.
  • C. F. Gauss, Disquisitiones Arithmeticae, 1801. English translation: Yale University Press, New Haven, CT, 1966, Sections 329-334.
  • G. B. Mathews, Theory of Numbers, Chelsea, no date, p. 263.
  • Paulo Ribenboim, My Numbers, My Friends, Chapter 11, Springer-Verlag, NY, 2000.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 142-143.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 103.
  • A. Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see pp. 188, 219-226.

Crossrefs

Sequence A025052 is a subsequence.
Cf. A139642 (congruences for idoneal quadratic forms).

Programs

  • Mathematica
    noSol={}; Do[lim=Ceiling[(n-2)/3]; found=False; Do[If[n>a*b && Mod[n-a*b, a+b]==0 && Quotient[n-a*b, a+b]>b, found=True; Break[]], {a, 1, lim-1}, {b, a+1, lim}]; If[ !found, AppendTo[noSol, n]], {n, 10000}]; noSol (* T. D. Noe, Apr 08 2004 *)
  • PARI
    A000926(Nmax=1e9)={for(n=1,Nmax,for(a=1,sqrtint(n\3),for(b=a+1,(n-a)\(3*a+2),n-a<(2*a+1+b)*b & break;(n-a*b)%(a+b)==0 & next(3)));print1(n", "))} \\ M. F. Hasler, Dec 04 2007
    
  • PARI
    ok(n)=!#select(k->k<>2, quadclassunit(-4*n).cyc) \\ Andrew Howroyd, Jun 08 2018

Extensions

Edited by N. J. A. Sloane, Dec 07 2007

A139643 Primes of the form x^2+Ny^2, with N=102.

Original entry on oeis.org

103, 127, 151, 223, 271, 409, 433, 457, 463, 577, 631, 727, 769, 919, 937, 967, 1033, 1039, 1063, 1087, 1249, 1279, 1327, 1447, 1471, 1543, 1657, 1753, 1759, 1777, 1783, 1801, 1879, 1951, 1993, 2089, 2143, 2161, 2287, 2311, 2473, 2503, 2551
Offset: 1

Views

Author

T. D. Noe, Apr 29 2008

Keywords

Comments

Discriminant=-408. N is an idoneal number (A000926), which means that the quadratic form's genus consists of a single class, which means that the primes of this form are identical to the primes that are congruent to c (mod 4N), where c is a set of numbers less than 4N. The sequence A139642 lists the set c for each idoneal number. That sequence also cross references the sequences for the quadratic forms with N equal to the first 36 idoneal numbers. The remaining quadratic forms are this sequence and the 28 listed in order below. Note that the sequences for N=120 and 240 are the same.
The primes are congruent to {1, 25, 49, 55, 103, 121, 127, 145, 151, 169, 217, 223, 247, 271, 319, 361} (mod 408).

References

  • David A. Cox, Primes of the Form x^2 + n y^2, Wiley, 1989.
  • L. E. Dickson, History of the Theory of Numbers, Vol 3, Chelsea, 1923.

Crossrefs

Programs

  • Magma
    [ p: p in PrimesUpTo(3000) | p mod 408 in {1, 25, 49, 55, 103, 121, 127, 145, 151, 169, 217, 223, 247, 271, 319, 361}]; // Vincenzo Librandi, Jul 28 2012
    
  • Magma
    k:=102; [p: p in PrimesUpTo(3000) | NormEquation(k, p) eq true]; // Bruno Berselli, Jun 01 2016
  • Maple
    C:= [1, 25, 49, 55, 103, 121, 127, 145, 151, 169, 217, 223, 247, 271, 319, 361]:
    select(isprime, [seq(seq(408*i+j,j=C),i=0..100)]); # Robert Israel, Jul 03 2016
  • Mathematica
    nn=102; pMax=10000; Union[Reap[Do[p=x^2+nn*y^2; If[p<=pMax && PrimeQ[p], Sow[p]], {x,Sqrt[pMax]}, {y, Sqrt[pMax/nn]}]][[2,1]]] (* T. D. Noe, Aug 02 2009 *)
    QuadPrimes2[1, 0, 102, 10000] (* see A106856 *)

A107008 Primes of the form x^2 + 24*y^2.

Original entry on oeis.org

73, 97, 193, 241, 313, 337, 409, 433, 457, 577, 601, 673, 769, 937, 1009, 1033, 1129, 1153, 1201, 1249, 1297, 1321, 1489, 1609, 1657, 1753, 1777, 1801, 1873, 1993, 2017, 2089, 2113, 2137, 2161, 2281, 2377, 2473, 2521, 2593, 2617, 2689, 2713
Offset: 1

Views

Author

T. D. Noe, May 09 2005

Keywords

Comments

Presumably this is the same as primes congruent to 1 mod 24, so a(n) = 24*A111174(n) + 1. - N. J. A. Sloane, Jul 11 2008. Checked for all terms up to 2 million. - Vladimir Joseph Stephan Orlovsky, May 18 2011.
Discriminant = -96.
Also primes of the forms x^2 + 48*y^2 and x^2 + 72*y^2. See A140633. - T. D. Noe, May 19 2008
Primes of the quadratic form are a subset of the primes congruent to 1 (mod 24). [Proof. For 0 <= x, y <= 23, the only values mod 24 that x^2 + 24*y^2 can take are 0, 1, 4, 9, 12 or 16. All of these r except 1 have gcd(r, 24) > 1 so if x^2 + 24*y^2 is prime its remainder mod 24 must be 1.] - David A. Corneth, Jun 08 2020
More advanced mathematics seems to be needed to determine whether this sequence lists all primes congruent to 1 (mod 24). Note the significance of 24 being a convenient number, as described in A000926. See also Sloane et al., Binary Quadratic Forms and OEIS, which explains how the table in A139642 may be used for this determination. - Peter Munn, Jun 21 2020
Primes == 1 (mod 2^3*3) are the intersection of the primes == 1 (mod 2^3) in A007519 and the primes == 1 (mod 3) in A002476, by the Chinese remainder theorem. - R. J. Mathar, Jun 11 2020

Crossrefs

Subset of A033199 (2y here = y there).
Is this the same as A141375?
See also the cross-references in A140633.

Programs

  • Mathematica
    QuadPrimes[1, 0, 24, 10000] (* see A106856 *)
  • PARI
    is(n) = isprime(n) && #qfbsolve(Qfb(1, 0, 24), n) == 2 \\ David A. Corneth, Jun 21 2020

Extensions

Recomputed b-file, deleted incorrect Mma program. - N. J. A. Sloane, Jun 08 2014
Showing 1-4 of 4 results.