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-10 of 23 results. Next

A106856 Primes of the form x^2 + xy + 2y^2, with x and y nonnegative.

Original entry on oeis.org

2, 11, 23, 37, 43, 53, 71, 79, 107, 109, 127, 137, 149, 151, 163, 193, 197, 211, 233, 239, 263, 281, 317, 331, 337, 373, 389, 401, 421, 431, 443, 463, 487, 491, 499, 541, 547, 557, 569, 599, 613, 617, 641, 653, 659, 673, 683, 739, 743, 751, 757, 809, 821
Offset: 1

Views

Author

T. D. Noe, May 09 2005, Apr 28 2008

Keywords

Comments

Discriminant=-7. Binary quadratic forms ax^2 + bxy + cy^2 have discriminant d = b^2 - 4ac.
Consider sequences of primes produced by forms with -100
The Mathematica function QuadPrimes2 is useful for finding the primes less than "lim" represented by the positive definite quadratic form ax^2 + bxy + cy^2 for any a, b and c satisfying a>0, c>0, and discriminant d<0. It does this by examining all x>=0 and y>=0 in the ellipse ax^2 + bxy + cy^2 <= lim. To find the primes generated by positive and negative x and y, compute the union of QuadPrimes2[a,b,c,lim] and QuadPrimes2[a,-b,c,lim]. - T. D. Noe, Sep 01 2009
For other programs see the "Binary Quadratic Forms and OEIS" link.

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.

Links

Crossrefs

Discriminants in the range -3 to -100: A007645 (d=-3), A002313 (d=-4), A045373, A106856 (d=-7), A033203 (d=-8), A056874, A106857 (d=-11), A002476 (d=-12), A033212, A106858-A106861 (d=-15), A002144, A002313 (d=-16), A106862-A106863 (d=-19), A033205, A106864-A106865 (d=-20), A106866-A106869 (d=-23), A033199, A084865 (d=-24), A002476, A106870 (d=-27), A033207 (d=-28), A033221, A106871-A106874 (d=-31), A007519, A007520, A106875-A106876 (d=-32), A106877-A106881 (d=-35), A040117, A068228, A106882 (d=-36), A033227, A106883-A106888 (d=-39), A033201, A106889 (d=-40), A106890-A106891 (d=-43), A033209, A106282, A106892-A106893 (d=-44), A033232, A106894-A106900 (d=-47), A068229 (d=-48), A106901-A106904 (d=-51), A033210, A106905-A106906 (d=-52), A033235, A106907-A106913 (d=-55), A033211, A106914-A106917 (d=-56), A106918-A106922 (d=-59), A033212, A106859 (d=-60), A106923-A106930 (d=-63), A007521, A106931 (d=-64), A106932-A106933 (d=-67), A033213, A106934-A106938 (d=-68), A033246, A106939-A106948 (d=-71), A106949-A106950 (d=-72), A033212, A106951-A106952 (d=-75), A033214, A106953-A106955 (d=-76), A033251, A106956-A106962 (d=-79), A047650, A106963-A106965 (d=-80), A106966-A106970 (d=-83), A033215, A102271, A102273, A106971-A106974 (d=-84), A033256, A106975-A106983 (d=-87), A033216, A106984 (d=-88), A106985-A106989 (d=-91), A033217 (d=-92), A033206, A106990-A107001 (d=-95), A107002-A107008 (d=-96), A107009-A107013 (d=-99).
Other collections of quadratic forms: A139643, A139827.
For a more comprehensive list of sequences giving numbers and/or primes represented by binary quadratic forms, see the "Binary Quadratic Forms and OEIS" link.
Cf. also A242660.

Programs

  • Mathematica
    QuadPrimes2[a_, b_, c_, lmt_] := Module[{p, d, lst = {}, xMax, yMax}, d = b^2 - 4a*c; If[a > 0 && c > 0 && d < 0, xMax = Sqrt[lmt/a]*(1+Abs[b]/Floor[Sqrt[-d]])]; Do[ If[ 4c*lmt + d*x^2 >= 0, yMax = ((-b)*x + Sqrt[4c*lmt + d*x^2])/(2c), yMax = 0 ]; Do[p = a*x^2 + b*x*y + c*y^2; If[ PrimeQ[ p]  && p <= lmt && !MemberQ[ lst, p], AppendTo[ lst, p]], {y, 0, yMax}], {x, 0, xMax}]; Sort[ lst]];
QuadPrimes2[1, 1, 2, 1000]
(This is a corrected version of the old, incorrect, program QuadPrimes. - N. J. A. Sloane, Jun 15 2014)
max = 1000; Table[yy = {y, 1, Floor[Sqrt[8 max - 7 x^2]/4 - x/4]}; Table[ x^2 + x y + 2 y^2, yy // Evaluate], {x, 0, Floor[Sqrt[max]]}] // Flatten // Union // Select[#, PrimeQ]& (* Jean-François Alcover, Oct 04 2018 *)
  • PARI
    list(lim)=my(q=Qfb(1,1,2), v=List([2])); forprime(p=2, lim, if(vecmin(qfbsolve(q, p))>0, listput(v,p))); Vec(v) \\ Charles R Greathouse IV, Aug 05 2016

Extensions

Removed old Mathematica programs - T. D. Noe, Sep 09 2009
Edited (pointed out error in QuadPrimes, added new version of program, checked and extended b-file). - N. J. A. Sloane, Jun 06 2014

A139490 Numbers n such that the quadratic form x^2 + n*x*y + y^2 represents exactly the same primes as the quadratic form x^2 + m*y^2 for some m.

Original entry on oeis.org

1, 4, 6, 7, 8, 10, 14, 16, 18, 22, 26, 38, 58, 82, 86
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008, Apr 26 2008, Apr 27 2008

Keywords

Comments

For the numbers m see A139491.
Conjecture: This sequence is finite and complete (checked for range n<=200 and m<=500).
Three more terms were found by searching n <= 1000 and m <= 4000. The corresponding m are 840, 840, and 1848, which are idoneal numbers A000926. The sequence is probably complete now. [T. D. Noe, Apr 27 2009]

Examples

			a(1)=1 because the primes represented by x^2+xy+y^2 are the same as the primes represented by x^2 + 3*y^2 (see A007645).
The known pairs (n,m) are the following (checked for range n<=200 and m<=500):
n={1, 4, 4, 6, 6, 7, 8, 8, 10, 10, 10, 14, 14, 14, 16, 18, 22, 22, 26, 38, 38}
m={3, 9, 12, 8, 16, 15, 45, 60, 24, 48, 72, 24, 48, 72, 21, 40, 120, 240, 168, 120, 240}.

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139491.

Programs

  • Mathematica
    f = 200; g = 300; h = 30; j = 100; b = {}; Do[a = {}; Do[Do[If[PrimeQ[x^2 + n y^2], AppendTo[a, x^2 + n y^2]], {x, 0, g}], {y, 1, g}]; AppendTo[b, Take[Union[a], h]], {n, 1, f}]; Print[b]; c = {}; Do[a = {}; Do[Do[If[PrimeQ[n^2 + w*n*m + m^2], AppendTo[a, n^2 + w*n*m + m^2]], {n, m, g}], {m, 1, g}]; AppendTo[c, Take[Union[a], h]], {w, 1, j}]; Print[c]; bb = {}; cc = {}; Do[Do[If[b[[p]] == c[[q]], AppendTo[bb, p]; AppendTo[cc, q]], {p, 1, f}], {q, 1, j}]; Union[cc] (*Artur Jasinski*)

Extensions

Edited by N. J. A. Sloane, Apr 25 2008
Extended by T. D. Noe, Apr 27 2009
Typo fixed by Charles R Greathouse IV, Oct 28 2009

A139492 Primes of the form x^2 + 5x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

7, 37, 43, 67, 79, 109, 127, 151, 163, 193, 211, 277, 331, 337, 373, 379, 421, 457, 463, 487, 499, 541, 547, 571, 613, 631, 673, 709, 739, 751, 757, 823, 877, 883, 907, 919, 967, 991, 1009, 1033, 1051, 1087, 1093, 1117, 1129, 1171, 1201, 1213, 1297, 1303
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Comments

Reduced form is [1, 3, -3]. Discriminant = 21. Class number = 2.
Values of the quadratic form are {0, 1, 3, 4} mod 6, so this is a subsequence of A002476. - R. J. Mathar, Jul 30 2008
It can be checked that the primes p of the form x^2 + n*x*y + y^2, n >= 3, where x and y are nonnegative, depend on n mod 6 as follows: n mod 6 = 0 => p mod 12 = {1,5}; n mod 6 = 1 => p mod 12 = {1,7}; n mod 6 = 2 => p mod 12 = {1}; n mod 6 = 3 => p mod 12 = {1,5,7,11}; n mod 6 = 4 => p mod 12 = {1}; n mod 6 = 5 => p mod 12 = {1,7}. - Walter Kehowski, Jun 01 2008

Examples

			a(1) = 7 because we can write 7 = 1^2 + 5*1*1 + 1^2.

References

  • Z. I. Borevich and I. R. Shafarevich, Number Theory.
  • David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.

Links

Crossrefs

Primes in A243172.
Cf. A002476, A007645, A007519, A033212, A033215, A068228, A107008, A107145, A107152, A139490, A139489, A139491, A139492, A141159.
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
    a = {}; w = 5; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a]
  • Sage
    # uses[binaryQF]
# The function binaryQF is defined in the link 'Binary Quadratic Forms'.
Q = binaryQF([1, 5, 1])
print(Q.represented_positives(1303, 'prime')) # Peter Luschny, May 12 2021

A139502 Primes of the form x^2 + 22x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

241, 409, 601, 769, 1009, 1129, 1201, 1249, 1321, 1489, 1609, 1801, 2089, 2161, 2281, 2521, 2689, 3001, 3049, 3121, 3169, 3361, 3529, 3769, 3889, 4129, 4201, 4441, 4561, 4729, 4801, 4969, 5209, 5281, 5449, 5521, 5569, 5641, 5689, 5881, 6121, 6361, 6481
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Comments

Also primes of the form x^2 + 120y^2. - T. D. Noe, Apr 29 2008
Also primes of the form x^2+240y^2. See A140633. - T. D. Noe, May 19 2008
In base 12, the sequence is 181, 2X1, 421, 541, 701, 7X1, 841, 881, 921, X41, E21, 1061, 1261, 1301, 13X1, 1561, 1681, 18X1, 1921, 1981, 1X01, 1E41, 2061, 2221, 2301, 2481, 2521, 26X1, 2781, 28X1, 2941, 2X61, 3021, 3081, 31X1, 3241, 3281, 3321, 3361, 34X1, 3661, 3821, 3901, where X is 10 and E is 11. Moreover, the discriminant is 340. - Walter Kehowski, Jun 01 2008

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.
Cf. A139643.

Programs

  • Magma
    [ p: p in PrimesUpTo(7000) | p mod 120 in {1, 49}]; // Vincenzo Librandi, Jul 28 2012
  • Mathematica
    QuadPrimes2[1, 0, 120, 10000] (* see A106856 *)

Formula

The primes are congruent to {1, 49} (mod 120). - T. D. Noe, Apr 29 2008

A139494 Primes of the form x^2 + 11x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

13, 43, 61, 79, 103, 127, 139, 157, 181, 199, 211, 277, 283, 313, 337, 367, 373, 433, 439, 523, 547, 571, 601, 607, 673, 727, 751, 757, 823, 829, 859, 883, 907, 919, 937, 991, 997, 1039, 1063, 1069, 1093, 1117, 1153, 1171, 1213, 1231, 1249, 1291, 1297
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.

Programs

  • Mathematica
    a = {}; w = 11; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)

A139506 Primes of the form x^2 + 26x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

193, 337, 457, 673, 1009, 1033, 1129, 1201, 1297, 1801, 1873, 2017, 2137, 2377, 2473, 2521, 2689, 2713, 2857, 3049, 3217, 3313, 3361, 3529, 3697, 3889, 4057, 4153, 4201, 4561, 4657, 4729, 4993, 5209, 5233, 5569, 5737, 5881, 6073, 6217, 6337, 6553, 6577
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Comments

Also primes of the form x^2 + 168y^2. - T. D. Noe, Apr 29 2008
In base 12, the sequence is 141, 241, 321, 481, 701, 721, 7X1, 841, 901, 1061, 1101, 1201, 12X1, 1461, 1521, 1561, 1681, 16X1, 17X1, 1921, 1X41, 1E01, 1E41, 2061, 2181, 2301, 2421, 24X1, 2521, 2781, 2841, 28X1, 2X81, 3021, 3041, 3281, 33X1, 34X1, 3621, 3721, 3801, 3961, 3981, where X is 10 and E is 11. Moreover, the discriminant is 480. - Walter Kehowski, Jun 01 2008

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A107152, A107008, A033215, A107145, A139490, A139491.
Cf. A139643.

Programs

  • Mathematica
    a = {}; w = 26; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a]

Formula

The primes are congruent to {1, 25, 121} (mod 168). - T. D. Noe, Apr 29 2008

A139512 Primes of the form x^2 + 32*x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

229, 349, 409, 421, 661, 769, 829, 1021, 1069, 1249, 1381, 1429, 1549, 1789, 1801, 1861, 2089, 2161, 2269, 2389, 3001, 3061, 3109, 3181, 3229, 3469, 3889, 4021, 4129, 4201, 4441, 4861, 4909, 5101, 5449, 5521, 5869, 5881, 6121, 6469, 6481, 6529, 6781
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Comments

Are all terms == 1 mod 12? - Zak Seidov, Apr 25 2008
Yes: (i) all terms == 1 mod 3 because the quadratic form has terms == {0,1} mod 3 and the values ==0 mod 3 are not primes. (ii) all terms == 1 mod 4 because the quadratic form has terms == {0,1,2} mod 4 and the values = {0,2} mod 4 are not primes. By the Chinese remainder constructions for coprime 3 and 4 all prime terms are == 1 mod 12. - R. J. Mathar, Jun 10 2020

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.

Programs

  • Mathematica
    a = {}; w = 32; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)

A139505 Primes of the form x^2 + 25x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

151, 163, 307, 397, 409, 541, 547, 601, 673, 811, 823, 859, 967, 997, 1153, 1231, 1237, 1327, 1567, 1669, 1741, 1879, 2083, 2143, 2281, 2293, 2557, 2677, 2707, 2833, 2971, 3037, 3259, 3313, 3433, 3877, 4003, 4129, 4153, 4603, 4639, 4861, 4957, 5101, 5227
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.

Programs

  • Mathematica
    a = {}; w = 25; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)
With[{nn=80},Select[Union[#[[1]]^2+25#[[1]]#[[2]]+#[[2]]^2&/@Tuples[ Range[ 0,nn],2]],PrimeQ[#]&&#Harvey P. Dale, Feb 10 2020 *)

A139493 Primes of the form x^2 + 9x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

11, 23, 37, 53, 67, 71, 113, 137, 163, 179, 191, 317, 331, 379, 389, 401, 421, 443, 449, 463, 487, 499, 599, 617, 631, 641, 653, 683, 709, 751, 757, 823, 863, 883, 907, 911, 947, 977, 991, 1061, 1087, 1093, 1103, 1171, 1213, 1303, 1367, 1373, 1409, 1423
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Comments

This is a member of the family of sequences of primes of the forms x^2 + kxy + y^2.
See for k=1 A007645 = x^2+3y^2, k=2 squares no primes, k=3 A038872, k=4 A068228 = x^2+9y^2, k=5 A139492, k=6 A007519 = x^2+8y^2, k=7 A033212 = x^2+15y^2, k=8 A107152 = x^2+45y^2, k=9 A139493, k=10 A107008 = x^2+24y^2, k=11 A139494, k=12 A139495, k=13 A139496, k=14* = 10 A107008 = x^2+24y^2, k=15 A139497, k=16 A033215 = x^2+21y^2, k=17 A139498, k=18 A107145 = x^2+40y^2, k=19 A139499, k=20 A139500, k=21 A139501, k=22 A139502, k=23 A139503, k=24 A139504, k=25 A139505, k=26,A139506, k=27 A139507, k=28 A139508, k=29 A139509, k=30 A139510, k=31 A139511, k=32 A139512

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.

Programs

  • Mathematica
    a = {}; w = 9; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)

A139495 Primes of the form x^2 + 12x*y + y^2 for x and y nonnegative.

Original entry on oeis.org

29, 109, 149, 281, 389, 401, 421, 449, 541, 569, 641, 701, 709, 809, 821, 1009, 1061, 1129, 1201, 1229, 1289, 1381, 1409, 1429, 1481, 1549, 1621, 1709, 1789, 1801, 1901, 2069, 2081, 2129, 2221, 2269, 2381, 2389, 2521, 2549, 2689, 2741, 2801, 2909, 2969
Offset: 1

Views

Author

Artur Jasinski, Apr 24 2008

Keywords

Links

Crossrefs

Cf. A139489, A007645, A068228, A007519, A033212, A033212, A107152, A107008, A033215, A107145, A139490, A139491.

Programs

  • Mathematica
    a = {}; w = 12; k = 1; Do[Do[If[PrimeQ[n^2 + w*n*m + k*m^2], AppendTo[a, n^2 + w*n*m + k*m^2]], {n, m, 400}], {m, 1, 400}]; Union[a] (*Artur Jasinski*)
With[{nn=50},Take[Union[Select[#[[1]]^2+12#[[1]]#[[2]]+#[[2]]^2&/@ Tuples[ Range[ nn],2],PrimeQ]],nn]] (* Harvey P. Dale, Dec 18 2015 *)
Showing 1-10 of 23 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.