A094376 -id:A094376 - cp's OEIS frontend

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

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

N. J. A. Sloane

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

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.

S. Chowla, An extension of Heilbronn's class number theorem, Quart. J. math., 5 (1934), 304-307.
K. S. Brown, Mathpages, Numeri Idonei
Günther Frei, Les nombres convenables de Leonhard Euler, Publications Université de Besançon, 1983-1984.
Günther Frei, Euler's convenient numbers, Math. Intell. Vol. 7 No. 3 (1985), 55-58 and 64.
E. Hertel, C. Richter, Tiling Convex Polygons with Congruent Equilateral Triangles, Discrete & Computational Geometry, 2014, DOI 10.1007/s00454-014-9576-7. Mentions this sequence. - _N. J. A. Sloane_, Mar 17 2014
O.-H. Keller, Über die "Numeri idonei" von Euler, Beitraege Algebra Geom., 16 (1983), 79-91. [Math. Rev. 85m:11019]
Robert Krzyzanowski, Euler's Convenient Numbers
David Masser, Alan Baker, arXiv:2010.10256 [math.HO], 2020. See p. 24.
Eric Rains, Comments on A000926
Paulo Ribenboim, Galimatias Arithmeticae, in Mathematics Magazine 71(5) 339 1998 MAA.
Rick L. Shepherd, Binary quadratic forms and genus theory, Master of Arts Thesis, University of North Carolina at Greensboro, 2013.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
J. Steinig, On Euler's ideoneal numbers, Elemente Math., 21 (1966), 73-88.
M. Waldschmidt, Open Diophantine problems, arXiv:math/0312440 [math.NT], 2003-2004
P. Weinberger, Exponents of the class groups of complex quadratic fields, Acta Arith., 22 (1973), 117-124.
Eric Weisstein's World of Mathematics, Idoneal Number

Sequence A025052 is a subsequence.
Cf. A014556, A026501, A093669, A094376, A094377, A094378.
Cf. A139642 (congruences for idoneal quadratic forms).

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 *)

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

ok(n)=!#select(k->k<>2, quadclassunit(-4*n).cyc) \\ Andrew Howroyd, Jun 08 2018

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

A094377 Greatest number having exactly n representations as ab+ac+bc with 0 < a < b < c.

Original entry on oeis.org

1848, 193, 1012, 862, 3040, 2062, 4048, 3217, 7392, 4162, 7837, 8002, 12397, 13297, 14722, 16417, 21253, 21058, 30493, 27358, 34357, 34318, 47338, 40177, 50317, 39502, 61462, 62302, 73117, 83218, 106177, 67138, 92698, 102958, 134773, 111577, 112942, 121522, 104938, 96958, 151237, 166798, 150382, 139393, 190513, 129838

Offset: 0

T. D. Noe, Apr 28 2004

nonn

Numbers up to 250,000 were checked. Note that the Mathematica program computes A094376, A094377 and A094378, but outputs only this sequence.

			a(1) = 193 because 193 is the largest number with a unique representation: (a,b,c) = (4,7,15).

See A025052.

Cf. A000926 (n having no representations), A093669 (n having one representation), A094376, A094378.

cntMax=10; nSol=Table[{0, 0, 0}, {cntMax+1}]; Do[lim=Ceiling[(n-2)/3]; cnt=0; Do[If[n>a*b && Mod[n-a*b, a+b]==0 && Quotient[n-a*b, a+b]>b, cnt++; If[cnt>cntMax, Break[]]], {a, 1, lim-1}, {b, a+1, lim}]; If[cnt<=cntMax, If[nSol[[cnt+1, 1]]==0, nSol[[cnt+1, 1]]=n]; nSol[[cnt+1, 2]]=n; nSol[[cnt+1, 3]]++;], {n, 10000}]; Table[nSol[[i, 2]], {i, cntMax+1}]

More terms (using limit 10^6) from Joerg Arndt, Oct 01 2017

A094378 Number of numbers having exactly n representations as ab+ac+bc with 0 < a < b < c.

Original entry on oeis.org

65, 23, 91, 40, 197, 39, 195, 56, 298, 87, 217, 60, 512, 97, 327, 77, 562, 125, 433, 88, 712, 125, 484, 115, 924, 121

Offset: 0

T. D. Noe, Apr 28 2004

nonn

Numbers up to 250,000 were checked. Note that there seem to be many more numbers having an even number of representations. Note that the Mathematica program computes A094376, A094377 and A094378, but outputs only this sequence.

			a(1) = 23 because there are 23 numbers (A093669) with unique representations.

See A025052

Cf. A000926 (n having no representations), A093669 (n having one representation), A094376, A094377.

cntMax=10; nSol=Table[{0, 0, 0}, {cntMax+1}]; Do[lim=Ceiling[(n-2)/3]; cnt=0; Do[If[n>a*b && Mod[n-a*b, a+b]==0 && Quotient[n-a*b, a+b]>b, cnt++; If[cnt>cntMax, Break[]]], {a, 1, lim-1}, {b, a+1, lim}]; If[cnt<=cntMax, If[nSol[[cnt+1, 1]]==0, nSol[[cnt+1, 1]]=n]; nSol[[cnt+1, 2]]=n; nSol[[cnt+1, 3]]++;], {n, 10000}]; Table[nSol[[i, 3]], {i, cntMax+1}]

A290870 a(n) is the number of ways to represent n as n = xy + yz + z*x where 0 < x < y < z.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 2, 0, 0, 2, 1, 0, 2, 0, 2, 1, 0, 1, 2, 1, 0, 2, 2, 0, 3, 0, 1, 3, 0, 1, 4, 0, 1, 2, 2, 1, 2, 2, 2, 3, 0, 0, 5, 0, 2, 3, 2, 1, 2, 2, 1, 4, 2, 0, 6, 0, 1, 4, 2, 3, 2, 0, 4, 3, 2, 1, 5, 2, 0, 4, 4, 0, 5, 2, 2, 4, 0, 3, 6, 2, 1, 3, 3, 1

Offset: 1

Joerg Arndt, Aug 13 2017

nonn

a(n) = 0 if and only if n is a term of A000926.

a(n) = 1 if and only if n is a term of A093669.

			For (x, y, z) = (1, 3, 5), we have x * y + y * z + z * x = 1 * 3 + 3 * 5 + 5 * 1 = 23 and similarily for (x, y, z) = (1, 2, 7), we have x * y + y * z + z * x = 23. Those 2 triples are all for n=23, so a(23) = 2. - _David A. Corneth_, Oct 01 2017

Joerg Arndt, Table of n, a(n) for n = 1..10000

Cf. A066955 (ways to represent n as n = x*y + y*z + z*x where 0 <= x <= y <= z).
Cf. A094377 (greatest number having exactly n representations).
Cf. A094376 (indices of records).

N=10^3; V=vector(N);
{ for (x=1, N,
    for (y=x+1, N, t=x*y; if( t > N, break() );
      for (z=y+1, N,
        tt = t + y*z + z*x;  if( tt > N, break() );
        V[tt]+=1;
); ); ); }
V \\ Joerg Arndt, Oct 01 2017

a(n) = {my(res = 0);
for(x = 1, sqrtint(n\3), for(y = x + 1, (n - x^2) \ (2 * x), z = (n - x*y) / (x + y); if(z > y && z == z\1, res++))); res} \\ David A. Corneth, Oct 01 2017

For the triples (x,y,z) we have x < sqrt(n / 3), y < (n - x^2) / (2 * x), z = (n - x*y) / (x + y) which must be integer. - David A. Corneth, Oct 01 2017

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Extensions

A094377 Greatest number having exactly n representations as ab+ac+bc with 0 < a < b < c.

Views

Author

Keywords

Comments

Examples

References

Crossrefs

Programs

Mathematica

Extensions

A094378 Number of numbers having exactly n representations as ab+ac+bc with 0 < a < b < c.

Views

Author

Keywords

Comments

Examples

References

Crossrefs

Programs

Mathematica

A290870 a(n) is the number of ways to represent n as n = x*y + y*z + z*x where 0 < x < y < z.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula

A290870 a(n) is the number of ways to represent n as n = xy + yz + z*x where 0 < x < y < z.