Next Article in Journal
Exponential Multistep Methods for Stiff Delay Differential Equations
Previous Article in Journal
Some Implicational Semilinear Gaggle Logics: (Dual) Residuated-Connected Logics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Solution to a Cubic Diophantine Equation

Department of Mathematics, University of Auckland, Auckland 1010, New Zealand
Axioms 2022, 11(5), 184; https://doi.org/10.3390/axioms11050184
Submission received: 24 March 2022 / Revised: 15 April 2022 / Accepted: 18 April 2022 / Published: 19 April 2022

Abstract

:
A positive integer, which can be written as the sum of two positive cubes in two different ways, is known as a “Ramanujan number”. The most famous example is 1729 = 10 3 + 9 3 = 12 3 + 1 3 , which was identified by Ramanujan as the lowest such number. In this paper, we consider the homogeneous cubic Diophantine equation x 3 + y 3 = u 3 + v 3 , where there is no restriction on the signs of the integers x , y , u , v . We show that every solution can be written in terms of two parameters in the ring Z 3 . It is also shown that solutions with arbitrarily high values of max ( | x | , | y | , | u | , | v | ) arise amongst the primitive solutions.

1. Introduction

The aim of this paper is to give a new solution to the Diophantine equation
x 3 + y 3 = u 3 + v 3 , x , y , u , v Z .
The result found in this paper gives the general solution in terms of two parameters in the ring R : = Z 3 .
Depending on the signs of the integers in any particular solution of (1), we obtain a solution to one or the other of
x 3 + y 3 = u 3 + v 3 , x , y , u , v Z + ,
x 3 + y 3 + z 3 = w 3 , x , y , z , w Z + .
For example, the solution [ x , y , u , v ] = [ 10 , 1 , 9 , 12 ] to (1) implies the solution [ x , y , u , v ] = [ 10 , 9 , 1 , 12 ] to (2), and the solution [ x , y , u , v ] = [ 6 , 4 , 5 , 3 ] to (1) implies the solution [ x , y , z , w ] = [ 3 , 4 , 5 , 6 ] to (3). The values of x 3 + y 3 , in solutions to (2), are known as “Ramanujan numbers”. The historical event that made (1) into a famous problem was recounted by C. P. Snow in his foreword to [1]. In [2], a method was given for generating an unlimited set of solutions by iteration. This differs markedly from the present paper, which generates all solutions parametrically. The problem was also discussed in [3].
In Section 2, we consider some properties of R, and in Section 3, we present a solution to (1). This is followed by Section 3, where the cubic Diophantine Equation (Section 3) is considered, and this leads to an algorithm for generating all possible solutions in Section 4. Finally, in Section 5, we show how to construct primitive solutions with arbitrarily large values.

2. The Ring Z 3

For ξ = a + b 3 R , the norm is defined as ξ = ξ ξ ¯ = a 2 + 3 b 2 . If ξ = a + b 3 , η = c + d 3 , then ξ η = ( a c 3 b d ) + ( a d + b c ) 3 .
We consider the question: “which odd primes greater than 3 are equal to ξ for some ξ ”?
Lemma 1.
Let p be a prime greater than 3; then, there exists a + b 3 R with p a , p b , such that
a 2 + 3 b 2 0 mod p ,
if and only if p 1 mod 6 .
Proof. 
Choose the integer c such that b c = 1 mod p and multiply (4) by c and write d = a c . It follows that 3 d 2 mod p , implying that, using Legendre symbols,
3 p = 1 p 3 p = ( 1 ) p 1 2 ( 1 ) p 1 2 p 3 = p 3 = 1 , p 1 mod 6 , 1 , p 5 mod 6 .
Denote P as the set of primes that are congruent to 1 mod 6 . For a given p P , denote N = { 1 , 2 , , p 1 2 } .
Lemma 2.
Let p P ; then, there exist integers a , b N such that
a 2 + 3 b 2 0 mod p .
Proof. 
Replace a by a in (5), where a = n p ± a with a N and similarly for b. □
Theorem 1.
Let p P ; then, a unique | a | and | b | exist such that
a 2 + 3 b 2 = p ,
Proof. 
The proof is by induction for p = 7 , 13 , 19 , P . For p = 7 , it is verified that 2 2 + 3 ( 1 ) 2 = p . It remains for us to prove the result for p P , p > 7 , where it is assumed to hold for all q P , such that q < p , From Lemma 2, a , b exist such that a 2 + 3 b 2 = p Q , with Q a product of prime powers. We show that this statement is also true with Q = 1 . It is not possible that p Q because
a 2 + 3 b 2 ( p 1 2 ) 2 + 3 ( p 1 2 ) 2 < p 2 .
The proof continues by considering a series of propositions that make it possible to exclude all possible prime divisors of Q. In each case, for a possible divisor g, we show that there exist a , b such that a 2 + 3 b 2 = p Q with Q < Q . In this proof only, the symbol ⇒ is used to introduce values of a , b , Q .
  • gcd ( a , b ) = g > 1 a = a / g , b = b / q . Q = Q / g 2 .
  • 2 Q .
    (a)
    a , b even. This reduces to 1.
    (b)
    a , b odd.
    • 4 ( a + b ) a = ( a 3 b ) / 4 , b = ( a + b ) / 4 . Q = ( Q / 4 )
    • 4 ( a b ) a = ( a + 3 b ) / 4 , b = ( a b ) / 4 . Q = ( Q / 4 )
  • 3 Q . This implies 3 a a = b , b = a / 3 . Q = Q / 3 .
  • q Q , for q prime and q 5 mod 6 . It follows that case 1 applies with g = q .
  • q Q , for q P .
    (a)
    q a , q b . Reduces to case 1.
    (b)
    q a , q b . ( c + d 3 ) ( c d 3 ) ( a + b 3 ) = p q Q .
    • q ( c + d 3 ) ( a + b 3 ) a = ( a c 3 b d ) / q , b = a d + b c , Q = Q / q
    • q ( c d 3 ) ( a + b 3 ) a = ( a c + 3 b d ) / q , b = a d b c , Q = Q / q
We now consider the primes of R and factorization in R.
Remark 1.
The primes of R are
1 ± 3
3
Prime numbers p 5 mod 6
a ± b 3 , with a , b Z + , and a 2 + 3 b 2 P .
Remark 2.
The product of two primes of R has a unique factorization with the exception of
4 = 2 2 = ( 1 + 3 ) ( 1 3 ) .

3. Solution to Cubic Equation

We first state the solution to an introductory equation.
Lemma 3.
The general solution to
A B = C D , A , B , C , D Z +
is
A = i j , B = k , C = i k , D = j , i , j , k , Z + .
Proof. 
Let i = gcd ( A , C ) , j = A / i , = gcd ( B , D ) , k = B / . □
Consider two Diophantine equations
x 3 + y 3 = u 3 + v 3 ,
a ( a 2 + 3 b 2 ) = c ( c 2 + 3 d 2 )
In each case, we consider primitive solutions.
The next result is easily verified and is given without proof.
Lemma 4.
( x , y , u , v ) is a primitive solution to (9) if and only if ( a , b , c , d ) is a primitive solution to (10), where
a b = T x y , c d = T u v ,
with
T = 1 2 1 1 1 1 , x + y odd , 1 2 1 1 1 1 , x + y even , T 1 = 1 2 1 1 1 1 , x + y odd , 1 2 1 1 1 1 , x + y even .
We now focus on (10), but written in a slightly different form.
a ˜ ( a 2 + 3 b 2 ) = c ˜ ( c 2 + 3 d 2 )
a ˜ = a , c ˜ = c .
We use Lemma 3 to solve (11) and then identify the free parameters by requiring them to satisfy (12).
Lemma 5.
The general solution to (11) is
a ˜ = r Y ,
a 2 + 3 b 2 = X Z ,
c ˜ = r X ,
c 2 + 3 d 2 = Y Z ,
where
X = ξ ξ ¯ = α 2 + 3 β 2 , ξ = α + β 3 ,
Y = η η ¯ = γ 2 + 3 δ 2 , η = γ + δ 3 ,
Z = ζ ζ ¯ = s 2 + 3 t 2 , ζ = s + t 3 ,
and r , s , t , α , β , γ , δ are parameters.
Proof. 
Use Lemma 3 and note that Z = gcd ( α 2 + 3 β 2 , γ 2 + 3 δ 2 ) , which is necessarily of the form s 2 + 3 t 2 . □
To eliminate r , s , t , use the conditions a ˜ = a , c ˜ = c to obtain
a ˜ = r ζ ζ ¯ = r ( γ 2 + 3 δ 2 ) ,
c ˜ = r ξ ξ ¯ = r ( α 2 + 3 β 2 ) ,
a = Re ξ ζ = s α 3 t β ,
c = Re η ζ = s γ 3 t δ ,
so that r , s , t satisfy the homogeneous linear system
γ 2 + 3 δ 2 α 3 β α 2 + 3 β 2 γ 3 δ r s t = 0 ,
with solution
r s t = G 1 3 ( β γ α δ ) 9 ( β 3 δ 3 ) + 3 ( α 2 β γ 2 δ ) α 3 γ 3 + 3 ( α β 2 γ δ 2 ) ,
where G 1 is introduced to reduce gcd ( r , s , t ) to 1, if necessary.

4. Algorithm and Sample Solutions

Using the results of Lemmas 4, 5 and Discussion Section 3, an algorithm can be constructed for finding solutions to
x 3 + y 3 = u 3 + v 3 .
  • Choose α + β 3 , γ + δ 3 R ;
  • Evaluate r , s , t ;
  • Evaluate a , b , c , d ;
  • Evaluate x , y , u , v .
It is found computationally that there are exactly 25 solutions to (2) with x 3 + y 3 1000 , 000 and 31 solutions to (3) with w 100 . These are shown in Table 1.

5. Extreme Values

Theorem 2.
Primitive solutions to (2), with x , y , u , v > 0 , exist with arbitrarily high values of x 3 + y 3 = u 3 + v 3 .
Proof. 
For [ α , β , γ , δ ] = [ n , 1 , n , 0 ] with 3 n , details of the solution are found to be
[ r , s , t ] = [ n , n 2 + 3 , n ] , x = n 3 + 2 n 2 + 3 , y = n 3 2 n 2 3 , u = n 3 + n 2 + 3 n , v = n 3 n 2 + 3 n , x 3 + y 3 = u 3 + v 3 = 2 n 3 ( n 2 + 3 ) ( n 4 + 9 n 2 + 9 ) .
Theorem 3.
Primitive solutions to (3), with x , y , z > 0 , exist with arbitrarily high values of w.
Proof. 
For [ α , β , γ , δ ] = [ n , 1 , 0 , n ] , with 3 n , details of the solution are found to be
[ r , s , t ] = [ 3 n 2 , 3 n 2 + 9 , n 2 + 3 n ] , x = 6 n 4 3 n 3 9 n 2 9 n , y = 8 n 4 + 9 n 3 6 n 2 9 , z = 10 n 4 9 n 3 + 6 n 2 + 9 , w = 12 n 4 3 n 3 + 9 n 2 9 n ,

6. Conclusions

A new parametric solution to the homogeneous cubic equation
x 3 + y 3 = u 3 + v 3 , x , y , u , v Z ,
is derived using a pair of parameters in the ring Z ( 3 ) . As a consequence, it is shown that, amongst solutions of (26), there exist arbitrarily large values of N = x 3 + y 3 , x , y , u , v > 0 . Furthermore, it is shown that, amongst solutions of x 3 + y 3 + z 3 = w 3 , x , y , z , w > 0 , there exist arbitrarily large values of w.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hardy, G.H. A Mathematician’s Apology; Cambridge University Press: Cambridge, UK, 1940. [Google Scholar]
  2. Goormaghtigh, R. On a Cubic Diophantine Equation. Math. Gazette 1937, 427, 21. [Google Scholar]
  3. Silverman, J.H. Taxicabs and sums of two cubes. Am. Math. Mon. 1993, 100, 331–340. [Google Scholar] [CrossRef]
Table 1. Solutions to (2) (left-hand table) and (3) (right-hand table).
Table 1. Solutions to (2) (left-hand table) and (3) (right-hand table).
α β γ δ x y u v x 3 + y 3 α β γ δ x y z w
021−1121910172901103456
0111915216410402111689
11−2−11924271020,683012−13101819
0120331523439,31210−117141720
011−1331693440,03310214172225
21−3−22636173964,23224−3118192128
12−223133401265,728012111152729
111238435112134,37911212174041
11−1−32950853149,38910−1216234144
12−3117555424171,28813−2−13363746
11−2−42257958195,841132029344453
13332259360216,02710−2312195354
021−351583067327,76315−3−115424958
131456616942402,59713−3122515467
21226948765439,10111−5236386169
23−5238737617443,88912−247545770
011271541580515,375144−214237071
13−2−475648251684,019215−434396572
20−214186892704,97712−2038436675
033−293113092805,688133−131337276
417363849423842,75112−6−125487481
22−4−133969720920,67317−5−219606982
011−263892498955,01611−3328537584
152−498355992984,06711−4150616485
12−3−269929929994,688113−326557887
10−2221438488
112−117408689
124225388790
114−332548593
245−119539096
212−345697997
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Butcher, J.C. A New Solution to a Cubic Diophantine Equation. Axioms 2022, 11, 184. https://doi.org/10.3390/axioms11050184

AMA Style

Butcher JC. A New Solution to a Cubic Diophantine Equation. Axioms. 2022; 11(5):184. https://doi.org/10.3390/axioms11050184

Chicago/Turabian Style

Butcher, John C. 2022. "A New Solution to a Cubic Diophantine Equation" Axioms 11, no. 5: 184. https://doi.org/10.3390/axioms11050184

APA Style

Butcher, J. C. (2022). A New Solution to a Cubic Diophantine Equation. Axioms, 11(5), 184. https://doi.org/10.3390/axioms11050184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop