Home Quartic Polynomials with a Given Discriminant
Article Open Access

Quartic Polynomials with a Given Discriminant

  • Jiří Klaška EMAIL logo
Published/Copyright: February 16, 2022
Become an author with De Gruyter Brill

Abstract

Let 0 ≠ D ∈ ℤ and let QD be the set of all monic quartic polynomials x4 + ax3 + bx2 + cx + d ∈ ℤ[x] with the discriminant equal to D. In this paper we will devise a method for determining the set QD. Our method is strongly related to the theory of integral points on elliptic curves. The well-known Mordell’s equation plays an important role as well in our considerations. Finally, some new conjectures will be included inspired by extensive calculations on a computer.

2020 Mathematics Subject Classification: Primary 11D25; 11D45; 11Y50

1. Introduction

Let 0 ≠ D ∈ ℤ and let

(1.1)QD={f(x)=x4+ax3+bx2+cx+d[x];Df=D}

where

(1.2)Df=a2b2c24a2b3d4a3c3+18a3bcd27a4d24b3c2+16b4d+18abc380ab2cd6a2c2d+144a2bd227c4+144bc2d128b2d2192acd2+256d3

is the discriminant of f(x). In this paper, the set QD will be studied in detail. Most of the focus will be given to the problem of determining all polynomials in QD. Clearly, this is equivalent to finding all integer solutions of the Diophantine equation Df = D. In proving the main results, the following two known theorems will be needed.

THEOREM 1.1

(Mordell, 1920). For any given 0 ≠ k ∈ ℤ, the equation

(1.3)Y2=X3+k
has at most finitely many integer solutions.

Equation (1.3) is often called Mordell’s equation, in honour of the contribution Louis Joel Mordell [17] has made to this subject. An extension to Theorem 1.1 was later made by Carl Ludwig Siegel [18]. In its simplest form, Siegel’s result can be formulated as follows:

THEOREM 1.2

(Siegel, 1929). Let α, β ∈ ℤ be such that 4α3 + 27β2 ≠ 0. Then the equation

(1.4)η2=ξ3+αξ+β
has at most finitely many integer solutions.

There is a standard method for computing all integer solutions of (1.3) and (1.4) using David’s bounds and lattice reduction. This method can be found, for example, in [19]. At present, this method is implemented in several computer algebra packages, including Magma and Pari (Sage).

Remark 1.3.

Mordell’s equation has had a long history. First discoveries concerning (1.3) were given in Dickson [2: pp. 533–539] going back to the work of Bachet from 1621. Many interesting historical notes to (1.3) can be found in [1,5,7,16]. Perhaps the most extensive historical comments related to Mordell’s contribution to (1.3) can be found in the recent paper [6].

Throughout this paper, the following notation will be adopted. If A is a finite set, #A denotes the number of elements of A.

2. Equivalence on the set QD

Let f(x) = x4 + ax3 + bx2 + cx + d ∈ ℤ[x] and let Df be the discriminant of f(x).

Next, let rf(x) = f(xa/4). Then

(2.1)rf(x)=x4+Ax2+Bx+C[x]

where

(2.2)A=b3a28,B=cab2+a38,C=dac4+a2b163a4256.

Moreover, we have

(2.3)Drf=Df=16A4C4A3B227B4128A2C2+144AB2C+256C3.

From (2.2), it follows that there exist R, S, T ∈ ℤ such that

(2.4)A=R8,B=S8,C=T256,

where

(2.5)R=8b3a2,S=8c4ab+a3,T=256d64ac+16a2b3a4.

Hence, we can write (2.1) in the form

(2.6)rf(x)=x4+R8x2+S8x+T256[x]withR,S,T.

We start with a more general theorem.

THEOREM 2.1

Let n ∈ ℕ, n ≥ 2, 0 ≠ D ∈ ℤ and let

f(x)=xn+an1xn1++a1x+a0[x]
be an arbitrary polynomial with the discriminant equal to D. Further, for any w ∈ ℤ, let
(2.7)fw(x)=k=0nf(k)(w)k!xk,
where f(k)(w) denotes the k-th derivative of f(x) at w. Then fw(x) ∈ ℤ[x] and all polynomials in {fw(x); w ∈ ℤ} have the same discriminant equal to D.

Proof. First, by induction on k, it can be proved that k!|f(k)(w) for any k ∈ {0, 1, 2,… }. Hence, fw(x) ∈ ℤ[x]. Further, Taylor’s theorem yields

(2.8)f(x)=k=0nf(k)(w)k!(xw)kfor anyw.

Let α1,…,αn be the roots of f(x) in the set of complex numbers ℂ. Then

D=i=1n1j=i+1n(αjαi)2.

Next, by (2.8), for any α ∈ {α1,…, αn}, we have

(2.9)f(α)=k=0nf(k)(w)k!(αw)k=0.

Combining (2.7) with (2.9), we get fw(αw) = 0, and thus,

(2.10)β1=α1w,,βn=αnw

are the roots of fw(x) in ℂ. Using (2.10), we now get

Dfw=i=1n1j=i+1n(βjβi)2=i=1n1j=i+1n(αjw(αiw))2=i=1n1j=i+1n(αjαi)2=D,

as desired. □

Remark 2.2.

Observe that, in Theorem 2.1, f(x) = f0(x). Hence, f(x) ∈ {fw(x); w ∈ ℤ}.

COROLLARY 2.3.

Let 0 ≠ D ∈ ℤ and let f(x) = x4 + ax3 + bx2 + cx + dQD. Further, for any w ∈ ℤ, let

(2.11)fw(x)=x4+f(w)3!x3+f(w)2!x2+f(w)1!x+f(w).
Then (i) and (ii) hold:
  1. QD is an infinite set and {fw(x); w ∈ ℤ} ⊆ QD.

  2. For any w ∈ ℤ, we have rf(x) = rf(x) = x4 + Ax2 + Bx + C ∈ ℚ[x], where A, B, C satisfy (2.2).

Proof. Part (i) of Corollary 2.3 is a direct consequence of Theorem 2.1 for n = 4. Part (ii) can be verified by direct calculation. □

LEMMA 2.4.

Let 0 ≠ D ∈ ℤ and let f(x), g(x) ∈ QD. Then (i), (ii) and (iii) are equivalent:

  1. There exists w ∈ ℤ satisfying g(x) = f(x + w).

  2. There exists w ∈ ℤ satisfying g(x) = fw(x).

  3. rf(x) = rg(x).

Proof. Let f(x) = x4 + ax3 + bx2 + cx + d, g(x) = x4 + ax3 + bx2 + cx + dQD.

First we show that (i) is equivalent to (ii). Using Taylor’s theorem, we obtain

fw(x)=(xw)4+f(w)3!(xw)3+f(w)2!(xw)2+f(w)1!(xw)+f(w)

for any w ∈ ℤ. Therefore,

(2.12)f(x+w)=x4+f(w)3!x3+f(w)2!x2+f(w)1!x+f(w).

Combining (2.12) with (2.11), we get f(x + w) = fw(x). Hence, (i) and (ii) are equivalent.

Further we prove that (i) is equivalent to (iii). Assume that g(x) = f(x + w) for some w ∈ ℤ. Then (2.12) yields

g(x)=x4+(4w+a)x3+(6w2+3aw+b)x2+(4w3+3aw2+2bw+c)x+w4+aw3+bw2+cw+d.

Hence, rg(x) = g(x – (4w + a)/4) = g(xwa/4) = f(xwa/4 + w) = f(xa/4) = rf(x).

Finally, let rf(x) = rg(x). Then f(xa/4) = g(x/4). Hence, f(xa/4 + /4) = g(x/4 + /4) = g(x) and g(x) = f(x – (a)/4) follows. Put w = (a)/4. Clearly, if a (mod 4), then w ∈ ℤ. Suppose that a (mod 4). Using (2.5) we obtain R = 8b – 3a2 = 8 – 32, S = 8c – 4ab + a3 = 8 – 4a̅b̅ + 3, which implies a22 (mod 8) and a3a̅3 (mod 4). Therefore, a22 (mod 4), which yields, without loss of generality, that either a ≡ 0 (mod 4), ≡ 2 (mod 4) or a ≡ 1 (mod 4), ≡ 3 (mod 4). If a ≡ 0 (mod 4), ≡ 2 (mod 4), then a2 ≡ 0 (mod 8), 2 ≡ 4 (mod 8), which is in contradiction to a22 (mod 8). Similarly, if a ≡ 1 (mod 4), ≡ 3 (mod 4), then a3 ≡ 1 (mod 4), 3 ≡ 3 (mod 4), which is in contradiction to a33 (mod 4). □

Let 0 ≠ D ∈ ℤ and let QD ≠ ∅. For f(x), g(x) ∈ QD put

f(x)g(x)w:g(x)=f(x+w)=fw(x)rf(x)=rg(x).

It is evident that ∼ is an equivalence relation on the set QD. Moreover, QD/∼ has only finitely many equivalence classes. In Section 4, this fact will be proved using the results of Mordell and Siegel presented in Theorem 1.1 and Theorem 1.2. On the other hand, this claim also follows as a consequence of a more general theorem that has been proved by Kálmán Györy [8: p. 419]. See also [9: p. 475] or consult [3: p. 109].

3. Connection between Mordell’s equation Y2 = X3 – 21633D and the set QD

THEOREM 3.1

Let 0 ≠ D ∈ ℤ. If Mordell’s equation

(3.1)Y2=X3+kwithk=1769472D=21633D
has no integer solution, then QD = ϕ.

Proof. Let f(x) = x4 + ax3 + bx2 + cx + dQD and let rf(x) = x4 + Ax2 + Bx + C ∈ ℚ[x]. Direct calculation will verify that (2.3) can be written in the form

(3.2)Drf=427(A2+12C)3127(2A372AC+27B2)2.

Substituting (2.4) into (3.2), after short calculation, we obtain

(3.3)Drf=11769472((R2+3T)3(R39RT+108S2)2).

Put

(3.4)X=R2+3TandY=R39RT+108S2.

Then X, Y ∈ ℤ and (3.3) yields

Y2=X3+kwherek=1769472Drf=21633Drf.

Since Drf = Df = D, the proof is complete.

Remark 3.2.

If x4 + ax3 + bx2 + cx + dQD, then (2.2) yields A, B, C ∈ ℤ ⇔ 4|a. In this case, we can write (3.2) in the form V2 = 4U3 – 27Drf, where U = A2 + 12C and V = 2A3 –72AC+27B2.

Hence, we have

(3.5)(4V)2=(4U)3432Drf.

Since Drf = D, the substitutions X = 4U, Y = 4V reduce (3.5) to

(3.6)Y2=X3432D.

It is interesting that Mordell’s equation (3.6) plays a fundamental role also in the theory of cubic polynomials with the same discriminant D. Consult [10: p. 313].

The following notation will be useful. For an arbitrary 0 ≠ D ∈ ℤ, let MD denote the set of all [X0, Y0], where X0, Y0 ∈ ℤ and Y02=X0321633D.

LEMMA 3.3.

Let 0 ≠ D ∈ ℤ and let [X0, Y0] ∈ MD. Then (i), (ii), (iii) and (iv) hold:

  1. If 2|X0, then 4|X0, 8|Y0.

  2. If 2|Y0, then 4|X0, 8|Y0.

  3. If 3|X0, then 9|Y0.

  4. If 3|Y0, then 3|X0, 9|Y0.

Proof. The conclusions (i)–(iv) immediately follow from Y02=X0321633D. □

They will be used in Section 4 and Section 5. □

4. Method for determining the set QD

The next lemma will be needed in the proof of Theorem 4.2.

LEMMA 4.1

Letξ0, η0, e ∈ ℤ be such that

(4.1)ξ036e2(mod96)andη09eξ0108e3(mod1728).
Then we have:
  1. ξ0 ≡ 0 (mod 12) andη0 ≡ 0 (mod 216).

  2. There exists exactly one e ∈ {0, 1, 2, 3} satisfying (4.1).

Proof. (i) Since the validity of the congruence ξ0 ≡ 0 (mod 12) is evident, we only prove that η0 ≡ 0 (mod 216). First, observe that η0 ≡ 90 – 108e3 (mod 216). Further, ξ0 ≡ 36e2 (mod 96) is equivalent to 9ξ0 ≡ 324e2 (mod 864). Hence, 9ξ0 ≡ 108e2 (mod 216). This, together with η0 ≡ 90 – 108e3 (mod 216), yields η0 ≡ 0 (mod 216).

(ii) Let ξ0, η0 ∈ ℤ satisfy (4.1) for some e ∈ {0, 1, 2, 3}. Suppose that e is not unique. Then it follows from ξ0 ≡ 36e2 (mod 96) that e ∈ {1, 3} and that ξ0 ≡ 36 (mod 96). On the other hand, using η0 ≡ 90 – 108e3 (mod 1728), we obtain 9ξ0 – 108 ≡ 27ξ0 – 2916 (mod 1728), which yields ξ0 ≡ 60 (mod 96), a contradiction. □

The following Theorem 4.2 provides the necessary and sufficient condition for QD ≠ ∅. In addition, Theorem 4.2 makes it possible to determine a particular polynomial in QD.

THEOREM 4.2.

Let 0 ≠ D ∈ ℤ and let MD ≠ ∅. Then QD ≠ ∅ if and only if there exists an [X0, Y0] ∈ MD such that the elliptic equation

(4.2)η2=ξ3108X0ξ+432Y0
has at least one integer solution [ξ0, η0] satisfying conditions (4.3)–(4.5)
(4.3)36e2ξ00(mod96),
(4.4)108e39eξ0+η00(mod1728),
(4.5)432e4ξ0272e2ξ0+16eη0+144X00(mod110592).
for some e ∈ {0, 1, 2, 3}. In this case.,
g(x)=x4+ex3+36e2ξ096x2+108e39eξ0+η01728x+432e4ξ0272e2ξ0+16eη0+144X0110592QD
and
rg(x)=x4ξ096x2+η01728x+144X0ξ02110592.

Proof. First, assume that QD ≠ ∅. Then there exists an f(x) = x4 + ax3 + bx2 + cx + dQD such that rf(x) = x4 + (R/8)x2 + (S/8)x + T/256 ∈ ℚ[x] where R, S, T are integers satisfying (2.5). Further, from Theorem 3.1 it follows that there exists a [X0, Y0] ∈ MD such that R2 + 3T = X0 and R3 – 9RT + 108S2 = Y0. Substituting 3T = X0R2 into R3 – 9RT + 108S2 = Y0, we obtain

(4.6)4R33X0R+108S2Y0,

and multiplying (4.6) by 432, we get

(4.7)(216S)2=(12R)3108X0(12R)+432Y0.

Put ξ0 = –12R and η0 = 216S. Now, (4.7) implies immediately that [ξ0, η0] is an integer solution of (4.2).

Finally, we have to prove that [ξ0, η0] satisfies (4.3)–(4.5) for some e ∈ {0, 1, 2, 3}. Since a ∈ ℤ, there exist uniquely determined w ∈ ℤ and e ∈ {0, 1, 2, 3} such that a = 4w + e. Substituting a = 4w + e into the first equation of (2.5), we obtain R ≡ –3e2 (mod 8) and –12R ≡ 36e2 (mod 96) follows. This together with ξ0 = –12R yields ξ0 ≡ 36e2 (mod 96). Hence, (4.3).

Further, from the second equation of (2.5), it follows

(4.8)216S=1728c864ab+216a3.

Putting a = 4w + e, 8b = R + 3a2, ξ0 = –12R and η0 = 216S into (4.8), we obtain

(4.9)η0=1728(c4w33ew2)+36w(ξ036e2)+9eξ0108e3.

Reducing (4.9) by modulus 1728 and using ξ0 ≡ 36e2 (mod 96), we get (4.4).

Finally, the third equation of (2.5) implies

(4.10)432T=110592d27648ac+6912a2b1296a4.

For the left-hand side of (4.10), we have 432T=144(X0R2)=144X0ξ02 and the right-hand side of (4.10) can be rewritten, substituting a = 4w+e, 8b = R+3a2, 8c = S+4aba3, ξ0 = –12R and η0 = 216S into

110596(dw4ew3)64w(η09eξ0+108e3)1152w2(36e2ξ0)16eη0+72e2ξ0432e2.

Since η0 ≡ 9eo – 108e3 (mod 1728) and ξ0 ≡ 36e2 (mod 96), we get

144X0ξ0216eη0+72e2ξ0432e4(mod110592).

Hence, (4.5).

Conversely, assume that there exists a [X0, Y0] ∈ MD such that equation (4.2) has an integer solution [ξ0, η0] satisfying (4.3)–(4.5) for some e ∈ {0, 1, 2, 3}. Put

(4.11)R=ξ012,S=η0216,T=144X0ξ02432.

Then, by part (i) of Lemma 4.1, we have R, S ∈ ℤ. We now prove that T ∈ ℤ. From the first and third equation in (4.11) we obtain T = (X0R2)/3. First we show that

(4.12)X00(mod3)R0(mod3).

Let 3|X0. Then, by part (iii) of Lemma 3.3, we have 9|Y0. Further, by (4.11), we have 3|ξ0 and 33|η0. Since η02=ξ02108X0ξ0+432Y0, we also have 0η02ξ03 (mod 35) and ξ0 ≡ 0 (mod 32) follows. This together with ξ0 = – 12R yields 3|R.

Let 3|R. Since ξ0 = –12R, we have 32|ξ0 and 36|ξ03 follows. Next, by (4.11), 36|η03. Since η02ξ03108X0ξ0+432Y0, we have 432Y0 ≡ 0 (mod 35), and Y0 ≡ 0 (mod 32) follows. By part (iv) of Lemma 3.3, we get 3|X0. This proves (4.12).

Further, suppose that X0 ≡ 2 (mod 3). Then from [X0, Y0] ∈ MD it follows that Y022 (mod 3), which is a contradiction. Combining this fact with (4.12), we get

(4.13)X00(mod3)R1(mod3)orR2(mod3)R21(mod3).

Clearly, in both cases (4.12) and (4.13), we have X0R2 ≡ 0 (mod 3). Hence, T ∈ ℤ.

Consider now the polynomial

r(x)=x4+R8x2+S8x+T256[x].

We prove that the discriminant Dr of r(x) is equal to D. First, direct calculation verifies that

Dr=(R2+3T)3(R39RT+108S2)221633.

On the other hand, substituting ξ0 = -12R, η0 = 216S into η02=ξ03108X0ξ0+432Y0 we obtain

(216S)2=(12R)3108X0(12R)+432Y0.

Hence, we get R3 –3R(X0R2)+108S2 = Y0. Since, X0R2 = 3T, we have R3–9RT +108S2 = Y0. This, together with Y02X0321633D yields

D=(R2+3T)3(R39RT+108S2)221633.

Hence, Dr = D.

Finally, let e ∈ {0, 1, 2, 3} satisfy (4.3)–(4.5). Then, by part (ii) of Lemma 4.1, e is uniquely determined. Put g(x) = r(x + e/4). Then we obtain after some calculation that

g(x)=x4+ex3+R+3e28x2+eR+2S+e316x+2e2R+8eS+T+e4256=x4+ex3+36e2ξ096x2+108e39eξ0+η01728x+432e4ξ0272e2ξ0+16eη0+144X0110592

and

rg(x)=r(x)=x4ξ096x2+η01728x+144X0ξ02110592.

Since Dg = Drg = Dr = D, we have g(x) ∈ QD, as desired. The proof is complete. □

Before proceeding, the following notations will be adopted. For any [X0, Y0] ∈ MD, let ED(X0, Y0) denote the set of all [ξ0, η0] where ξ0, η0 ∈ ℤ and Y02X03108X0ξ0+432Y0. Next, let ED denote the set of all [X0, Y0, ξ0, η0, e] where [X0, Y0] ∈ MD, [ξ0, η0] ∈ ED(X0, Y0) and e ∈ {0, 1, 2, 3} satisfy (4.3)–(4.5).

COROLLARY 4.3.

Let 0 ≠ D ∈ ℤ and let f(x) = x4 + ax3 + bx2 + cx + d ∈ ℤ[x]. Then f(x) ∈ QD if and only if there exists [X0, Y0, ξ0, η0, e] ∈ ED and w ∈ ℤ such that

a=4w+e,b=6w2+3ew+36e2ξ096,c=4w3+3ew2+36e2ξ048w+108e39eξ0+η01728,d=w4+ew3+36e2ξ096w2+108e39eξ0+η01728w+432e4ξ0272e2ξ0+16eη0+144X0110592.
PROPOSITION 4.4.

Let 0 ≠ D ∈ ℤ and let MD ≠ ∅. Then (i), (ii) and (iii) hold:

  1. ED (X0, Y0) is a finite set for any [X0, Y0] ∈ MD.

  2. ED is a finite set.

  3. QD/∼ has only finitely many equivalence classes for any QD ≠ ∅.

Proof. (i) Put α = –108X0 and β = 432Y0. Then 4α3+27β2=2839(X03+Y02)=224312D0.

Conclusion (i) now follows from Theorem 1.2.

  1. Conclusion (ii) is a direct consequence of Theorem 1.1 and part (i) of Proposition 4.4.

  2. Let φ : EDQD/∼ be the mapping defined by φ(X0, Y0, ξ0, η0, e) = {fw(x); w∈ℤ}, where

    f0(x)=x4+ex3+36e2ξ096x2+108e39eξ0+η01728x+432e4ξ0272e2ξ0+16eη0+144X0110592.

Then φ is bijective. Injectivity of φ is evident and surjectivity of φ immediately follows from Corollary 4.3. Hence, #QD/∼= #ED. This proves (iii). □

Remark 4.5.

Let [X0, Y0], [X0*,Y0*]MD and let [X0*,Y0*]. By an example we will prove that the set ED(X0,Y0)ED(X0*,Y0*) can be nonempty. For D = –23, we have [64, 6400], [–320, –2816] ∈ M–23 and [96, ±1728] ∈ E–23(64, 6400) ∩ E–23(–320, –2816).

Now we are ready to formulate the method for determining the set QD. It can be formally divided into five steps as follows:

  1. Let 0 ≠ D ∈ ℤ. First we find the set MD of all integer solutions [X0, Y0] of Mordell’s equation Y2 = X3 – 21633D. By Theorem 1.1, MD is a finite set and Theorem 3.1 states that, if MD = ∅, then QD = ∅.

  2. Let MD ≠ ∅. Next we find, for any [X0, Y0] ∈ MD, the set ED(X0, Y0) of all integer solutions [ξ0, η0] of the elliptic equation η2 = ξ3 – 108X0+432Y0. By part (i) of Proposition 4.4, ED (X0, Y0) is a finite set for any [X0, Y0] ∈ MD and Theorem 4.2 says that, if ED(X0, Y0) = ∅ for any [X0, Y0] ∈ MD, then QD = ∅.

  3. In step (iii), we establish the set ED. By part (ii) of Proposition 4.4, ED is a finite set and Corollary 4.3 states that QD ≠ ∅ if and only if ED ≠ ∅.

  4. Let ED ≠ ∅ and let #ED = n. In this step, we assign to each [X0, Y0, ξ0, η0, e] ∈ ED the polynomial

    g(x)=x4+ex3+36e2ξ096x2+108e39eξ0+η01728x+432e4ξ0272e2ξ0+16eη0+144X0110592.

    In this way, we obtain the full system of representatives GD = {g1(x),…, gn(x)} of QD/∼. By part (iii) of Proposition 4.4, QD/∼ is a finite set.

  5. Finally, applying Corollary 2.3 to each gi(x) ∈ GD, i ∈ {1,…,n}, we obtain the n sets {fi,w(x); w ∈ ℤ} where

    fi,w(x)=x4+gi(w)3!x3+gi(w)2!x2+gi(w)1!x+gi(w).

    Hence, we get

    QD=i=1n{fi,w(x);w}.

The below example illustrates our method.

Example 4.6.

Let D = –87. Then we have

M87={[320,±11008],[92,±12376],[448,±15616]}.

Hence,

E87(320,11008)={[80,±1216],[48,±1728],[240,±5184],[384,±8640],[8592,±796608]},E87(320,11008)=,E87(92,12376)={[156,0]},E87(92,12376)={[156,0]},E87(448,15616)={[156,±3240],[96,±1728]},E87(448,15616)=.

Further, we have

E87={[320,11008,240,5184,2],[320,11008,240,5184,2],[448,15616,156,3240,1],[448,15616,156,3240,3]}.

Hence, it follows that #E–87 = #Q–87//∼= 4 and that G–87 = {g1(x), g2(x), g3(x), g4(x)} where

g1(x)=x4+2x3x2+x,g2(x)=x4+2x3x25x3,g3(x)=x4+x3+2x2x,g4(x)=x4+3x3+5x2+6x+3.

Finally,

f1,w(x)=x4+(4w+2)x3+(6w2+6w1)x2+(4w3+6w22w+1)x+w4+2w3w2+w,f2,w(x)=x4+(4w+2)x3+(6w2+6w1)x2+(4w3+6w22w5)x+w4+2w3w25w3,f3,w(x)=x4+(4w+1)x3+(6w2+3w+2)x2+(4w3+3w2+4w1)x+w4+w3+2w2w,f4,w(x)=x4+(4w+3)x3+(6w2+9w+5)x2+(4w3+9w2+10w+6)x+w4+3w3+5w2+6w+3,

and

Q87=i=14{fi,w(x);w}.

Applying the method, the validity of Theorem 4.7 can be verified.

THEOREM 4.7.

Let 0 ≠ D ∈ ℤ and let 1 ≤ |D| ≤ 1000. Then we have:

  1. If 1 ≤ D ≤ 1000, then QD ≠ ∅ if and only if

    D{5,12,20,21,32,40,45,48,49,60,77,81,85,96,104,112,117,125,140,144,148,165,169,189,192,216,221,224,229,252,256,257,260,272,285,288,320,321,333,357,361,392,400,404,432,437,468,469,473,480,488,500,512,525,528,533,544,549,564,572,580,592,605,621,629,656,672,697,725,729,733,761,768,785,788,792,816,832,837,864,892,896,900,916,957,981,985}
  2. If –1 ≥ D ≥ –1000, then QD ≠ ∅ if and only if

    D{3,16,23,27,31,44,59,76,83,87,107,108,112,135,139,140,175,176,199,211,231,236,240,247,255,256,268,275,279,283,288,304,331,332,335,351,367,400,416,428,432,448,464,475,491,500,507,527,556,560,563,575,588,608,643,671,684,688,695,731,751,783,800,816,844,848,863,864,891,931,944,959,972,976,983}.

5. Structure of the set QD

In this section, we establish some results related to the structure of the set QD. For e ∈ {0, 1, 2, 3}, put QD(e) = {x4 + ax3 + bx2 + cx + dQD; ae (mod 4)}. Then QD(0), QD(1), QD(2), QD(3) are pairwise disjoint sets, and

QD=e=04QD(e).

Proposition 5.1 gives the important link between the sets QD(1) and QD(3).

PROPOSITION 5.1

Let 0 ≠ D ∈ ℤ and let e ∈ {1, 3}. Then there exists a one-to-one correspondence between the sets QD(1) and QD (3) given by the relation

[X0,Y0,ξ0,η0,1]ED[X0,Y0,ξ0,η0,3]ED.
Consequently, QD(1) ≠ ∅ if and only if QD(3) ≠ ∅.

Proof. First observe that [ξ0, η0] ∈ ED (X0, Y0) if and only if [ξ0, –η0] ∈ ED (X0, Y0).

If e ∈ {1, 3}, then, by (4.3), ξ0 ≡ 36 (mod 96), which is equivalent to

(5.1)18ξ0648(mod1728).

Next, if e = 1, then, by (4.4), 9ξ0η0 – 108 ≡ 0 (mod 1728). Using (5.1), this congruence can be written in the equivalent form

1719ξ0+η0+10827ξ0+9418ξ0+η0+10827ξ0+η0+5400(mod1728).

Hence, [ξ0, η0] satisfies (4.4) for e = 1 if and only if [0, –η0] satisfies (4.4) for e = 3.

Furthermore, if e ∈ {1, 3}, then, by (4.3), ξ0 ≡ 36 (mod 96), which is equivalent to

(5.2)1152ξ0+414720(mod110592).

Next, if e = 1 then, from (4.4), it follows that η0 ≡ 9ξ0 – 108 (mod 1728) if and only if

(5.3)576ξ0+64η069120(mod110592).

Subtracting (5.3) from (5.2), we now obtain

(5.4)576ξ0+64η0345600(mod110592).

Finally, if e = 1, then (4.5) yields

(5.5)432ξ0272ξ0+64η0+144X00(mod110592).

Subtracting (5.4) from (5.5), we obtain

34993ξ02648ξ0+48η0+144X00(mod110592).

This proves that [X0, Y0, ξ0, η0, 1] ∈ ED if and only if [X0, Y0, ξ0, –η0, 3] ∈ ED, which implies QD(1) ≠ ∅ if and only if QD(3) ≠ ∅. □

For the remaining cases e ∈ {0, 2}, we can prove Proposition 5.2.

PROPOSITION 5.2.

Let 0 ≠ D ∈ ℤ and let e ∈ {0, 2}. Then we have:

  1. [X0, Y0, ξ0, η0, 0] ∈ ED ⇔ [X0, Y0, ξ0, –η0, 0] ∈ ED.

  2. [X0, Y0, ξ0, η0, 2] ∈ ED ⇔ [X0, Y0, ξ0, –η0, 2] ∈ ED.

Proof. Since part (i) of Proposition 5.2 immediately follows from (4.3)–(4.5), we prove (ii). Let [X0, Y0, ξ0, η0, 2] ∈ ED. Then (4.3)–(4.5) yields that ξ0 ≡ 48 (mod 96), η0 ≡ 18ξ0 – 864 (mod 1728), and

(5.6)6912ξ02288ξ0+32η0+144X00(mod110592).

Since ξ0 ≡ 48 (mod 96) is equivalent to 18ξ0 ≡ 864 (mod 1728), we have –η0 ≡ – 18ξ0 + 864 ≡ 1710ξ0 – 864 = 18ξ0 – 864 + 94 · 18ξ0 ≡ 18ξ0 – 864 + 47 · 1728 ≡ 18ξ0 – 864 (mod 1728). Further, ξ0 ≡ 48 (mod 96) is equivalent to 1152ξ0 ≡ 55296 (mod 110592) and η0 ≡ 18ξ0 – 864 (mod 1728) is equivalent to 64η0 ≡ 1152ξ0 – 55296 (mod 110592). Hence, we obtain 64η0 ≡ 0 (mod 110592) and 32η0 ≡ – 32η0 (mod 110592) follows. Now we see that (5.6) is equivalent to 6912 – ξ02 – 288ξ0 – 32η0 + 144X0 ≡ 0 (mod 110592). This proves (ii). □

Remark 5.3.

If QD ≠ ∅, then any of the below seven cases can occur:

  1. QD(0)≠ ∅, QD(1) ∪ QD(3) ≠ ∅, QD(2) ≠ ∅, D = –23,

  2. QD(0) ≠ ∅, QD(1) ∪ QD(3) ≠ ∅, QD(2) ≠ ∅, D = 32,

  3. QD(0) ≠ ∅, QD(1) ∪ QD(3) ≠ ∅, QD(2) ≠ ∅, D = 5,

  4. QD(0) = ∅, QD(1) ∪ QD(3) ≠ ∅, QD(2) ≠ ∅, D = –87,

  5. QD(0) ≠ ∅, QD(1) ∪ QD(3) = ∅, QD(2) ≠ ∅, D = –27,

  6. QD(0) = ∅, QD(1) ∪ QD(3) ≠ ∅, QD(2) = ∅, D = 12,

  7. QD(0) = ∅, QD(1) ∪ QD(3) = ∅, QD(2) ≠ ∅, D = –3.

The above values of D are the least, in absolute value, for which the case occurs.

6. Even and odd solutions of Mordell’s equation Y2 = X3 – 21633D

Some basic arithmetic properties of integer solutions [X0, Y0] of the Mordell’s equation Y2 = X3 – 21633D have already been presented in Lemma 3.3. Combining part (i) and (ii) of Lemma 3.3, we immediately get X0 ≡ 0 (mod 2) ⇔ Y0 ≡ 0 (mod 2). Hence, the following two definitions are possible:

  1. A solution [X0, Y0] ∈ MD is called even, if X0and Y0are even.

  2. A solution [X0, Y0] ∈ MD is called odd, if X0and Y0are odd.

Next, for any 0 ≠ D ∈ ℤ, let

D={[X0,Y0]MD:X0Y00(mod2)},OD={[X0,Y0]MD:X0Y01(mod2)}.

Then DOD=0 and DOD=MD. Finally, for any positive integer n, put

εn=D=1n#D,εn=D=1n#D,on=D=1n#OD,andon=D=1n#OD.

Computer investigation of the values εn, ε–n, on and o–n for n ≤ 1000 reveals a significant difference between the numbers of even and odd solutions in the investigated range. We have found

ε1000=1572,ε1000=1090,o1000=100,ando1000=44.

Hence, it follows that there exist approximately 95% even and only 5% odd integer solutions of Y2 = X3 – 21633D for 0 ≠ |D| ≤ 1000. This surprising fact inspires the study of even solutions in detail. As the main result of this section, we prove that, for any even solution [X0, Y0], equation (4.2) can be replaced by another elliptic equation whose integer coefficients are substantially smaller in the absolute value than in (4.2).

We begin by recalling the well-known proposition concerning the solubility of linear congruences. See, for example, Hardy and Wright [4: p. 62, Theorem 57].

PROPOSITION 6.1.

Let a, b, m ∈ ℤ, m > 1 and let g = gcd(a, m). Then the congruence axb (mod m) is soluble if and only if g|b.

Using Proposition 6.1, we now prove Lemma 6.2.

LEMMA 6.2.

If X0, Y0 ∈ ℤ, then the congruence 3αX0 + Y0 – 4α3 ≡ 0 (mod 27) holds for at most one α ∈ {0,1, 2}.

Proof. The proof consists of three steps. (i) First, suppose that Y0 ≡ 0 (mod 27) and 3X0 + Y0 – 4 ≡ 0 (mod 27). Then Y0 ≡ 4 – 3X0 ≡ 0 (mod 27), which yields 3X0 ≡ 4 (mod 27). Hence, gcd(3, 27) = 3 ∤ 4, which is a contradiction.

(ii) Further, suppose that Y0 ≡ 0 (mod 27) and 6X0 + Y0 – 32 ≡ 0 (mod 27). Then Y0 ≡ 5 – 6X0 ≡ 0 (mod 27), which yields 6X0 ≡ 5 (mod 27). Hence, gcd(6, 27) = 3 ∤ 5, which is a contradiction.

(iii) Finally, suppose that 3X0 + Y0 – 4 ≡ 0 (mod 27) and 6X0 + Y0 – 32 ≡ 0 (mod 27). Then Y0 ≡ 4 – 3X0 ≡ 5 – 6X0 (mod 27), which yields 3X0 ≡ 1 (mod 27). Hence, gcd(3, 27) = 3 ∤ 1, which is a contradiction. Combining (i)–(iii) proves the lemma. □

We are now ready to prove the main result of this section.

THEOREM 6.3.

Let 0 ≠ D ∈ ℤ and let [X0, Y0] ∈ D.

  1. If 3αX0 + Y0 – 4α3 ≡ 0 (mod 27) does not hold for any α ∈ {0, 1, 2}, then the system

    (6.1)R2+3T=X0andR39RT+108S2=Y0
    is not solvable in integers.
  2. If 3αX0 + Y0 – 4α3 ≡ 0 (mod 27) holds for some α ∈ {0, 1, 2}, thenα is uniquely determined, X0 – 4α2 ≡ 0 (mod 12), 3αX0 + Y0 – 4α3 ≡ 0 (mod 108) and the set K of all integer solutions of (6.1) can be obtained from the set L of all integer solutions of the elliptic equation

    (6.2)η2=ξ2αξ2X04α212ξ+3αX0+Y04α3108.

    Moreover, between K and L, there exists a one-to-one correspondence.

Proof. (i) Let [R0, S0, T0] be an arbitrary integer solution of (6.1). Since R0 ∈ ℤ, there exits a uniquely determined r ∈ ℤ and α ∈ {0, 1, 2} such that R0 = 3r + α.

First, we prove

(6.3)X04α2122|X0.

Let 2|X0. Then, by part (i) of Lemma 3.3, 4|X0, and thus 4|X0 – 4α2. Next, the first equation in (6.1) yields X0R02X0α20 (mod 3), and 3|X0 – 4α2 follows. Hence, 12|X0 – 4α2, which means (X0 – 4α2)/12 ∈ ℤ. Because the validity of the converse implication is evident, we get (6.3).

Further, substituting 3T0=X0R02 into R039R0T0+108S02Y0, we obtain

(6.4)4R033X0R0+108S02=Y0.

Since R0 = 3r + α, (6.4) can be written in the equivalent form

(6.4)108S02=108r3108αr2+9(X04α2)r+3αX0+Y04α3.

Reducing (6.5) by the modulus 108, using (6.3), we get 3αX0 + Y0 – 4α3 ≡ 0 (mod 108). This proves (i).

(ii) Assume that there exists an α ∈ {0,1, 2} such that 3αX0 + Y0 – 4α3 ≡ 0 (mod 27). Then Lemma 6.2 states that α is uniquely determined. Since [X0, Y0] ∈ D, by part (i) of Lemma 3.3, 4|X0 and 8|Y0, which yields 4|3αX0 + Y0 – 4α3 for any α ∈ {0, 1, 2}. Hence,

3αX0+Y04α30(mod27)3αX0+Y04α30(mod108).

We now prove that

(6.6)3αX0+Y04α30(mod27)X04α20(mod12).

First, by part (i) of Lemma 3.3, X0 – 4α2 ≡ 0 (mod 4). Next, from the assumption 3αX0 + Y0 – 4α3 ≡ 0 (mod 27), we get Y0α3α (mod 3). Hence, Y02α2 (mod 3). Furthermore, it is clear from [X0, Y0] ∈ D that Y02X03 (mod 3), which, together with X03X0 (mod 3), yields X0α2 (mod 3). Hence, X0 – 4α2 ≡ 0 (mod 3), and (6.6) follows.

Let [R0, S0, T0] ∈ K. Then there exist a uniquely determined r ∈ ℤ and an α ∈ {0,1, 2} such that R0 = 3r + α. In much the same way as in the proof of part (i), we get (6.5). Hence,

(6.7)S02=r3αr2+X04α212r+3αX0+Y04α3108

where (X0 – 4α2)/12 and (3αX0 + Y0 – 4α3)/108 are integers. Put [ξ0, η0] = [–r, S0]. Then it follows from (6.7) that [ξ0, η0] ∈ L.

Conversely, let [ξ0, η0] ∈ L where α ∈ {0,1, 2} satisfies 3αX0 + Y0 – 4α3 ≡ 0 (mod 27). Put R0 = –3ξ0 + α and S0 = η0. Then R0, S0 ∈ ℤ, and substituting ξ0 = (αR0)/3, η0 = S0 into (6.2) some calculation will yield

(6.8)R03+3R0(R02X0)+108S02=Y0.

Next, put T0=(X0R02)/3. Since R0 = –3ξ0 + α, we have T0 = 2α0 – 3ξ02 + (X0α2)/3, and by (6.6), T0 ∈ ℤ. Substituting R02X0=3T0 into (6.8), we obtain R039R0T0+108S02=Y0. This proves that [R0,S0,T0]=[3ξ0+α,η0,2αξ03ξ02+(X0α2)/3]K.

Moreover, it is evident that the mapping ψ: LK defined by

(6.9)ψ(ξ0,η0)=[3ξ0+α,η0,2αξ03ξ02+(X0α2)/3]

is bijective. The proof is complete. □

Example 6.4.

Let D = –87. Then [X0, Y0] = [–320,11008] ∈ 87, and for α = 1, we have 3αX0 + Y0 – 4α3 = 10044 ≡ 0 (mod 27). By Theorem 6.3, the set K of all integer solutions of the system R2 + 3T = –320 and R3 – 9RT + 108S2 = 11008 can be obtained using the set L of all integer solutions of the elliptic equation η2 = ξ3ξ2 + 27ξ + 93. Since L = {[–1, ±8], [7, ±24], [11, ±40], [239, ±3688]}, (6.9) yields

K={[4,±8,112],[20,±24,240],[32,±40,448],[716,±3688,170992]}.

On the other hand, by Section 4, the set K can also be determined by means of the set E–87(–320,11008) = {[–80, ±1216], [–48, ±1728], [240, ±5184], [384, ±8640], [8592, ±796608]} of all integer solutions of the elliptic equation η2 = ξ3 + 34560ξ + 4755456.

Remark 6.5.

If [X0, Y0] ∈ OD, then X04α20 (mod 12) for any α ∈ {0,1, 2}. On the other hand, by examples, it can be proved that both cases (i) and (ii) can occur:

(i)3αX0+Y04α30(mod108),(ii)3αX0+Y04α30(mod108).
  1. If D = –23, then [X0, Y0] = [489,12555] ∈ O23, and for any α ∈ {0, 1, 2}, we have 3αX0+Y04α30 (mod 108).

  2. If D = –107, then [X0, Y0] = [9241, 888445] ∈ O107, and for α = 1, we get 3αX0 + Y0 – 4α3 = 916164 ≡ 0 (mod 108).

7. Some conjectures related to Mordell’s equation

Let 0 ≠ D ∈ ℤ and let

μD={0ifMD=/0,1ifMD/0.

Next, for any positive integer n, put

σn=D=1nμDandσn=D=1nμD.

Computer investigation of the values σn and σ–n for n ≤ 1000 yields

σ10001000=2801000=0.280,σ10001000=4261000=0.426,σ1000σ1000=2804260.657.

Hence, the following conjectures can be made:

(7.1)limnσnn=270.286,limnσnn=370.429,limnσnσn=230.667.

The conjectures (7.1) lead to another interesting question, namely, whether some similar hypotheses can also be stated for the case of general Mordell’s equation. It is clear that, to formulate such hypotheses, much computation will be needed.

Thanks to the computations made by M. A. Bennett and A. Ghadermarzi [1], all integer solutions of Y2 = X3 + k are determined for any 0 ≠ |k| ≤ 107. Based on their results, some new conjectures can be formulated. The following notations will be useful.

For 0 ≠ k ∈ ℤ, let M(k) denote the set of all integer solutions of Y2 = X3 + k, and let

μ(k)={0ifM(k)=0,1ifM(k)0.

Next, for any positive integer n, put

σ(n)=k=1nμ(k)andσ(n)=k=1nμ(k).

By inspecting Table 1 and Table 2 in [1: pp. 642–643], we get

σ(107)107=13329341070.133,σ(107)107=8346041070.083,σ(107)σ(107)=83460413329340.626.

Hence, the following conjectures can be made:

(7.2)limnσ(n)n=215=0.13¯,limnσ(n)n=112=0.083¯,limnσ(n)σ(n)=58=0.625.

Our surmises can be proposed to the reader as Problem 7.1.

Problem 7.1.

Prove or disprove (7.1) and (7.2).

8. Conclusion

The results presented in this paper make it possible to determine the set QD of all monic quartic polynomials x4 + ax3 + bx2 + cx + d ∈ ℤ[x] with a given discriminant 0 ≠ D ∈ ℤ. That provides an opportunity to study the following problem: Establish a general method of deciding in a finite number of steps whether, for a given 0 ≠ D ∈ ℤ, the following statement holds: Let p be an arbitrary prime. Then all polynomials in QD have the same type of factorization over the Galois fieldFp. The validity of an analogous statement for the case of cubic polynomials has been recently examined in [10]–[15] with relatively closed results obtained.


Dedicated to the eminent Czechoslovak mathematician Ladislav Skula

(Communicated by Milan Paštéka)


Acknowledgement

The author thanks the anonymous referee for carefully reading the manuscript.

REFERENCES

[1] [1] BENNETT, M. A.—GHADERMARZI, A.: Mordell’s equation: a classical approach, LMS J. Comput. Math. 18.1 (2015), 633–646.10.1112/S1461157015000182Search in Google Scholar

[2] [2] DICKSON, L. E.: History of the Theory of Numbers, Vol. II, Chelsea, New York, 1952.Search in Google Scholar

[3] [3] EVERTSE, J. H.—GYÖRY, K.: Discriminant Equations in Diophantine Number Theory, Cambridge University Press, Cambridge, 2017.10.1017/CBO9781316160763Search in Google Scholar

[4] [4] HARDY, G. H.—WRIGHT, E. M.: An Introduction to the Theory of Numbers, 6th ed., Oxford University Press, New York, 2008.Search in Google Scholar

[5] [5] HEMER, O.: On the Diophantine Equation y2k = x3, Doctoral Dissertation, Uppsala (1952).10.1007/BF02589282Search in Google Scholar

[6] [6] GAUTHIER, S.—LÉ, F.: On the youthful writings of Louis J. Mordell on the Diophantine equation y2k = x3, Archive for History of Exact Sciences 73 (2019), 427–468.10.1007/s00407-019-00231-1Search in Google Scholar

[7] [7] GEBEL, J.—PETHÖ, A.—ZIMMER, G. H.: On Mordell’s equation, Compos. Math. 110 (1998), 335–367.10.1023/A:1000281602647Search in Google Scholar

[8] [8] GYÖRY, K.: Sur les polynômes à coefficients entiers et de discriminant donné, Acta Arith. 23 (1973), 419–426.10.4064/aa-23-4-419-426Search in Google Scholar

[9] [9] GYÖRY, K.: Polynomials and binary forms with given discriminant, Publ. Math. Debrecen 69 (4) (2006), 473–499.10.5486/PMD.2006.3745Search in Google Scholar

[10] [10] KLAŠKA, J.—SKULA, L.: Mordell’s equation and the Tribonacci family, Fibonacci Quart. 49 (4) (2011), 310–319.Search in Google Scholar

[11] [11] KLAŠKA, J.—SKULA, L.: Law of inertia for the factorization of cubic polynomials – the real case, Util. Math. 102 (2017), 39–50.Search in Google Scholar

[12] [12] KLAŠKA, J.—SKULA, L.: Law of inertia for the factorization of cubic polynomials – the imaginary case, Util. Math. 103 (2017), 99–109.Search in Google Scholar

[13] [13] KLAŠKA, J.—SKULA, L.: Law of inertia for the factorization of cubic polynomials – the case of discriminants divisible by three, Math. Slovaca 66 (4) (2016), 1019–1027.10.1515/ms-2015-0199Search in Google Scholar

[14] [14] KLAŠKA, J.—SKULA, L.: Law of inertia for the factorization of cubic polynomials – the case of primes 2 and 3, Math. Slovaca 67 (1) (2017), 71–82.10.1515/ms-2016-0248Search in Google Scholar

[15] [15] KLAŠKA, J.—SKULA, L.: On the factorizations of cubic polynomials with the same discriminant modulo a prime, Math. Slovaca 68 (5) (2018), 987–1000.10.1515/ms-2017-0161Search in Google Scholar

[16] [16] LONDON, J.—FINKELSTEIN, M.: On Mordell’s Equation y2k = x3, Bowling Green, Ohio Bowling Green State University, 1973.Search in Google Scholar

[17] [17] MORDELL, L. J.: A statement by Fermat, Proc. Lond. Math. Soc. (2) 18 (1920), pp. v–vi.Search in Google Scholar

[18] [18] SIEGEL, C. L.: Uber einige Anwendungen diophantischer Approximationen, Abh. Preuss Akad. Wiss., 1929, pp. 1–41.Search in Google Scholar

[19] [19] SMART, N. P.: The Algorithmic Resolution of Diophantine Equations, Cambridge University Press, Cambridge, 1998.10.1017/CBO9781107359994Search in Google Scholar

Received: 2020-11-16
Accepted: 2021-01-25
Published Online: 2022-02-16
Published in Print: 2022-02-16

© 2022 Mathematical Institute Slovak Academy of Sciences

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Downloaded on 13.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/ms-2022-0003/html
Scroll to top button