arxiv: v1 [math.nt] 29 Aug 2018

Transcription

1 A CUBIC GENERALIZATION OF BRAHMAGUPTA S IDENTITY arxiv: v1 [mathnt] 9 Aug 018 SAMUEL A HAMBLETON Abstract We give an algebraic identity for cubic polynomials which generalizes Brahmagupta s identity and facilitates arithmetic in cubic fields We also pose a question about a relationship between the elements of a cubic field of fixed trace and fixed norm and rational points of an elliptic curve 1 Introduction Brahmagupta s identity is an ancient Indian algebraic identity with several applications The identity is expressed as (1 (x 1 x +Dy 1 y D(x 1 y +x y 1 = ( x 1 Dy 1( x Dy When D = 1, we obtain a well known result on Pythagorean triples, triples of positive integers (x,y,z correspondingto the lengths of the sides of a right triangle so that x +y = z By (1, if (x 1,y 1,z 1 and (x,y,z are Pythagorean triples, then so is (x 3,y 3,z 3, where x 3 = x 1 x y 1 y, y 3 = x 1 y +x y 1, z 3 = z 1 z and denotes the absolute value Similarly, if (x 1,y 1 and (x,y satisfy the Pell equation ( x Dy = 1, where D is an integer, then by (1 so does (x 3,y 3, where (3 x 3 = x 1 x +Dy 1 y, y 3 = x 1 y +x y 1 This shows that if the Diophantine equation ( has a solution other than (±1,0, then more can be obtained When D > 0 is not a square, it is known that ( has infinitely many solutions Possibly the first known method for solving (, the chakravala method due to Bhaskara II, uses (1 See [7, 8, 9] for a description and history of the chakravala method See [14] for other methods of solving ( Consider the matrix given by ( xj Dy (4 M j = j y j x j The transpose of the matrix is called the Brahmagupta matrix; see [17] Brahmagupta s identity (1 is equivalent to the matrix identity M 1 M = M 3 Taking the determinants of the M j, for j = 1,,3, gives (1, and the left column of M 3 is the vector with entries given by (3 Date: Submitted to J R M S on 7 May, Mathematics Subject Classification Primary 11R16, 11D57; Secondary 11D5, 11G05 Key words and phrases Brahmagupta s identity, cubic fields, elliptic curves 1

2 SAMUEL A HAMBLETON The purpose of this article is to introduce an analogous 3 3 matrix such that matrix multiplication gives a similar identity to (1 for a ternary cubic form and gives (1 when we set specific integer values to the letters in this 3 3 matrix Before we introduce this matrix we will consider Brahmagupta s identity in the context of the quadratic field K = Q( D, where D is the discriminant of K Let O K = Z[ω], where ω = s+ D is the ring of integers of K, and s, equal to 0 or 1, is the remainder of D modulo 4 Then for every α O K, there are integers u,y such that α = u+yω We can express α as α = x+y D, where x = u+sy is the trace of α Taking the norm of α gives (5 N(α = ( x+y ( D x y D = 1 ( x Dy 4 We know that the norms of algebraic numbers are multiplicative and consequently if α,β K, then N(αN(β = N(αβ When we write α,β in the form of (5, and multiply by 4, we obtain (1 In fact, there is a bijection φ between the ring O K and the integer solutions (x,y,n to the equation x Dy = 4n, where D is the discriminant of a quadratic field, given by φ : O K { (x,y,n Z 3 : x Dy = 4n }, φ : α = x+y D (x,y,n(α; see [11] As a consequence of this bijection we are able to carry out basic arithmetic in O K and K using matrices with integer and rational entries respectively If instead of (4 we use ( uj my (6 N j = j u j +sy y j where m = D s 4, we have the convenient correspondence α j = u j + y j ω N j with several useful properties The trace of α j is the trace of N j, the norm of α j is the determinant of N j, α j O K if and only if N j has integer entries, α j K if any only if N j has rational entries, α 1 j can be calculated by finding the inverse of N j and considering the entries in the left column, and the eigenvalues of N j are x j±y j D Gauss [10, Art 34] generalized Brahmagupta s identity for the composition of binary quadratic forms, given by Q(x,y = Ax +Bxy +Cy, where A,B,C Z satisfying gcd(a,b,c = 1, B 4AC = D is the discriminant of a quadratic field, and x, y are indeterminants The binary quadratic form Q(x, y is abbreviated Q = (A,B,C Gauss composition law on the GL (Z classes of binary quadratic forms is known to be equivalent to the ideal class group of the quadratic field K The way that Gauss introduced composition is known as a bilinear transformation,, (7 Q 3 (u 3,y 3 = Q 1 (u 1,y 1 Q (u,y,

3 A CUBIC GENERALIZATION OF BRAHMAGUPTA S IDENTITY 3 where Q j (u j,y j = (A j,b j,c j, for j = 1,,3 Since we can always find β j Z such that B j = β j +s, and C j = β j +sβj m A j, we have u 3 = eu 1 u + e (β β 3 u 1 y + e (β 1 β 3 u y 1 + e ( β β 3 β + y 1 y, A A 1 A 1 A y 3 = A 1 e u 1y + A e u y 1 + β+ e y 1y, β + = β 1 +β +s, β = β 1 β +m, e = gcd ( A 1,A,β +, A 3 = A 1A e, and β 3 is the least non-negative integer satisfying β 3 b 3 (mod A 3, where p,q,r are integers satisfying e = A 1 p+a q +β + r, and b 3 = 1 e (A 1β p+a β 1 q +β r It is relatively easy to show that the coefficients of u 1 u, u 1 y, u y 1, and y 1 y in u 3 and y 3 are integers When A 1 = A = 1, B 1 = B = s, and C 1 = C = m, we have A 3 = 1, B 3 = s, and C 3 = m Letting u = x sy and multiplying (7 by 4 again gives (1 Thus Gauss composition may be considered as a quadratic generalization of Brahmagupta s identity A cubic identity Now we develop the analogous situation, for cubic fields and 3 3 matrices, of the relationship between elements of quadratic fields and (6 A binary cubic form C is a form given by (8 C(x,y = ax 3 +bx y +cxy +dy 3, where a,b,c,d Z, and C is irreducible over Q[x,y] Belabas [1] showed that there is a fast algorithm for compiling tables of binary cubic forms C whose discriminant (9 = b c +18abcd 4ac 3 4b 3 d 7a d is equal to that of the field K = Q(ζ, where C(ζ,1 = 0, and C belongs to the image of the Davenport-Heilbronn map [4] This map gives a bijection between the conjugacy class of the cubic field K and the GL (Z class of the binary cubic forms ofthe same discriminant as the field K belonging to the image of such a map, where the conjugacy class of the cubic field L = Q(α is mapped to the GL (Z class of the binary cubic form C L = 1 disc(l N ((ω ω x (ω 3 ω 3 y, {1,ω,ω 3 } is an integral basis of L, disc(l is the discriminant of L, and N is the norm It was also shown [1] that we may easily give the integral basis of the ring or integers O K of K in terms of the coefficients a,b of C of (8 as B = {1,aζ,aζ +bζ}, so that all α O K be be uniquely represented as (10 α j = u j +x j (aζ+y j ( aζ +bζ,

4 4 SAMUEL A HAMBLETON where u j,x j,y j Z We we call the binary cubic form C belonging to the image of the Davenport-Heilbronn map of K the canonical binary cubic form of the cubic field K Consider the matrix (11 N j = ( uj ady j adx j bdy j x j u j bx j cy j cx j dy j y j ax j u j cy j Thematrixgivenby(11correspondstoα j K satisfying(10, whenc = (a,b,c,d is the canonical binary cubic form of the cubic field K Under the correspondence between α j and N j, the trace of α j is equal to the trace of N j, the norm of α j is equal to the determinant of N j, the product α 1 α = α α 1 of elements of K corresponds to the matrix product N 1 N = N N 1, the sum α 1 + α corresponds to the sum N 1 +N, and the multiplicative inverse α 1 j corresponds to the matrix inverse N 1 j Thus the correspondence between α j and N j facilitates performing arithmetic in the cubic field K since we can do this using matrices It is easy to show that these claims are true once we consider the following result, which shows that multiplication of matrices of the form (11, whether or not (a,b,c,d is a binary cubic form, gives a matrix of the same form We also note that we can choose a, b, c, d so that this generalizes Brahmagupta s identity, Gauss bilinear transformation when the binary quadratic forms satisfy Q 1 = Q = (1,B,C, and some known identities on cubic polynomials Proposition 1 Let N j be given by (11, where a,b,c,d are fixed indeterminants and u j,x j,y j are indeterminants for j = 1,,3 Then the matrix product N 3 = N 1 N is commutative and of the form (11 We obtain Brahmagupta s identity by choosing a = 0, b = 1, c = 0, d = D, and taking the determinants of the N j Choosing a = 0, b = C, c = B, and d = 1 gives the special case of Gauss bilinear transformation for equal principal forms Q 1 = Q = (1,B,C Proof We begin by factorizing the matrix N j as the product N j = SUj T, where T denotes the matrix transpose and ad bd S = b c d, U j = u j x j y j u j 0 x j y j a 0 c 0 0 u j 0 x j y j Let u 3 = u 1 u adx y 1 adx 1 y bdy 1 y, x 3 = u 1 x +u x 1 bx x 1 cx 1 y cx y 1 dy 1 y, y 3 = u 1 y +u y 1 +ax 1 x cy 1 y We must show that S ( U3 T UT 1 SUT = [0]3 3, the 3 3 matrix with zero entries Expanding U3 T U1 T SU T gives adx y 1 bdy 1 y adx 1 y adu 1 y adu 1 x + bdu 1 y bx 1 x cx y 1 cx 1 y dy 1 y bu 1 x adx y 1 + cu 1 y bdy 1 y cu 1 x + du 1 y + adx 1 x + bdx 1 y ax 1 x cy 1 y au 1 x + ady 1 y cu 1 y adx 1 y x 1 x u 1 x cx y 1 dy 1 y cx 1 x + dx 1 y x 1 y x y 1 u 1 y + bx y 1 u 1 x bx 1 x y 1 y ax y 1 u 1 y + ax 1 x Multiplying on the left by S produces [0] 3 3 Since swapping u 1 and u, x 1 and x, y 1 and y does not change u 3, x 3, y 3, it follows that matrix multiplication of the N j is commutative

5 A CUBIC GENERALIZATION OF BRAHMAGUPTA S IDENTITY 5 To show that matrix multiplication generalizes Brahmagupta s identity, let a = 0, b = 1, c = 0, d = D Then for j = 1,,3, ( uj 0 Dy j (1 N j = x j u j x j Dy j, y j 0 u j and expanding the matrix N 3 = N 1 N gives u 1u +Dy 1y 0 D(u 1y +u y 1 N 3 = u x 1 x x 1 +u 1x +Dy 1y u 1u u x 1 u 1x +x 1x D(u 1y +u y 1 u y 1 +u 1y 0 u 1u +Dy 1y The determinants of N 1 and N are (u 1 x 1 ( u 1 Dy 1 and (u x ( u Dy The determinant of N 3 is the product of the determinants of N 1 and N Assuming u 1 x 1 and u x, we obtain Brahmagupta s identity To prove the claim about the special case of Gauss bilinear transformation, let a = 0, b = C, c = B, and d = 1 The determinant of N j is equal to Letting We have (u j Cx j +By j ( u j +Bu jy j +Cy j u 3 = (u 1 u Cy 1 y, x 3 = (u 1 x +u x 1 Cx x 1 +Bx 1 y +Bx y 1 y 1 y, y 3 = (u 1 y +u y 1 +By 1 y u 3 Cx 3 +By 3 = (u 1 Cx 1 +By 1 (u Cx +By Canceling this from the determinant of N 3 = N 1 N gives Q(u 3,y 3 = Q(u 1,y 1 Q(u,y Other choices of a, b, c, and d give identities on well known forms For example, choosing a = 1, b = 0, c = 0, and d = n, and taking the determinant of N j gives the form u 3 +nx 3 +n y 3 3nuxy studied by Lagrange and Carmichael [3] Consequently, it it easy to show that if there is one solution (u,x,y (1,0,0 to the Diophantine equation u 3 +nx 3 +n y 3 3nuxy = 1, then more may be generated via matrix multiplication of the N j Proposition Let K = Q(ζ be a cubic field, where ζ is a root of the polynomial C(x,1, C = (a,b,c,d is the canonical binary cubic form of the cubic field K, and let O K be the ring of integers of K generated by {1,ρ,ω}, where ρ = aζ, ω = aζ +bζ Let M Z = {N j = [a ij ] 3 3 : a ij Z} M Q = {N j = [a ij ] 3 3 : a ij Q} There is a ring isomorphism φ and a field isomorphism ψ given by φ : O K M Z, φ : α j N j, ψ : K M Q, ψ : α j N j, where α j = u j +x j ρ+y j ω and N j is defined by (11 The trace of α j is equal to the trace of N j and the norm of α j is equal to the determinant of N j

6 6 SAMUEL A HAMBLETON Proof First we show that φ is a bijection Let α 1,α O K and assume φ(α 1 = φ(α Then the left columns of the matrices φ(α j, for j = 1,, must be equal It follows that u 1 = u, x 1 = x, y 1 = y and we must have α 1 = α Therefore φ is injective To show that φ is surjective, Let N j M We have φ(u j +x j ρ+y j ω = N j and α j = u j + x j ρ + y j ω O K It follows that φ is bijective By transport of the ring structure of O K onto M, we obtain a ring isomorphism Similarly ψ is a field isomorphism Under these maps 0 and 1 map to the zero matrix and the identity matrix respectively To prove the claims about the trace and norm of α j, let τ be an embedding of K The trace and norm are defined t j = α j +τ (α j +τ (τ (α j, n j = α j τ (α j τ (τ (α j, where without loss of generality we assume τ : (ζ 1,ζ,ζ 3 (ζ,ζ 3,ζ 1 and we define ζ 1 = ζ, and let ζ, ζ 3 be the other two roots of C(x,1 By considering the identities ζ 1 +ζ +ζ 3 = b a, ζ1ζ +ζ1ζ3 +ζζ3 = c a, and ζ 1 +ζ +ζ 3 = b ac a, ζ 3 1 +ζ 3 +ζ 3 3 = b3 +3abc 3a d a 3, ζ 1ζ +ζ 1ζ 3 +ζ ζ 3 = c bd a, ζ 1ζ ζ 3 = d a, t j = 3u j bx j cy j, (ζ +ζ 3ζ1 +(ζ 1 +ζ 3ζ +(ζ 1 +ζ ζ3 = bd+3ad, a it is easy although somewhat tedious to show that the traces of α j and N j coincide and the norm of α j coincides with the determinant of N j We note that there is a concise way to represent the determinant of N j, (11 Delone and Faddeev [5, pp 130] considered various triangular forms, ternary cubic forms The determinant of N j is a ternary cubic form, not a Cayley nor a Dirichlet form but if we replace u u+cy, x x, and then take the determinant, we get the Dirichlet form When we change to variables t,x,y, where t is trace of N j and n is the determinant, instead of u,x,y, we obtain the equation (13 t 3 3tQ(x,y+F(x,y = 7n, where Q and F are respectively the Hessian and Jacobian covariant binary forms, of the binary cubic C = (a,b,c,d, defined using partial derivatives H(C = 1 ( ( Cxx C xy Qx Q, J C yx C (Q,C = y, yy C x C y Q(x,y = det(h(c, F(x,y = det(j (Q,C, and satisfy the syzygy (14 F +7 C = 4Q 3 ; see [13], used by Mordell [16] to show that the elliptic curve y = x 3 +k has finitely many integer solutions Note that the Hessian binary quadratic form is simply (15 Q = ( b 3ac,bc 9ad,c 3bd

7 A CUBIC GENERALIZATION OF BRAHMAGUPTA S IDENTITY 7 and while F is not defined by (14, it can be obtained from (9, (14 and (15 See [] for the relationship between binary cubic forms and the ideal class group of a cubic field It can be shown that when C is the canonical binary cubic form of a cubic field K and n = 1, then as a Diophantine equation (13 has a group law such that the group of integer points (t,x,y is isomorphic to the group of units of K of norm 1 and this group law may be evaluated by matrix multiplication of the N j Similarly we may restrict attention to the set of all algebraic numbers of a cubic field having trace 0 This set will form a group of trace 0 elements under addition It is natural to ask whether there is a group of elements of a cubic field having both trace 0 and norm 1 The answer is affirmative There is a relationship between certain rational points (x,y of (13 and certain elliptic curves but the group law cannot be viewed in terms of matrix multiplication or addition of the N j It shows that rational points of some elliptic curves may be thought of as corresponding to elements of a cubic field with a specific trace and specific norm See [15] for a discussion of a relationship between elliptic curves and some quadratic fields Remark 3 Let C = (a,b,c,d be the canonical binary cubic form of the cubic field K and let Q and F the Hessian and Jacobian covariant binary forms of C respectively There is a bijection between those elements of K of fixed trace t and fixed norm n and the rational points (x,y of the elliptic curve (13 provided that the choice of fixed t and n make (13 a non-singular curve containing a rational point When t 3 = 7n, the curve is singular We now discuss the transformation of (13 to Weierstrass form following Duif [6] In homogeneous coordinates we have (16 Γ(X,Y,Z = F(X,Y 3tZQ(X,Y+ ( t 3 7n Z 3 = 0, where x = X Z, y = Y Z Let P = (P X : P Y : P Z be a point of (16 in homogeneous coordinates Since Γ is smooth, the tangent line at P is (17 Γ X (P(X P X +Γ Y (P(Y P Y +Γ Z (P(Z P Z = 0 The intersections of this line with (16 may either be of multiplicity 3 at P, Case 1, or of multiplicity at P and at another point Q of multiplicity 1, Case In the first case, we map P to (0 : 1 : 0 and the line (17 to the line Z = 0 We fix a point Q P on (17, not on (16 and satisfying det(m 0, where M = ( QX P Z 1 Q Y P Y 0 Q Z P Z 0 ( We send Q to (1 : 0 : 0 The matrix M is invertible since det(m 0 Letting UV ( = M 1 XYZ, we have T ku 3 +pu T +quvt +rv T +sut +wvt +ut 3 = 0, where k,p,q,r,s,w,u Q Dividing by k, which must be non-zero, and letting T = kw r gives the projective Weierstrass form We see the affine Weierstrass form when we then replace W by 1 Now assume that we have a point P of (16 and (17 intersecting (16 with multiplicity Let Q be the other intersection of (16 and (17 We take a new projective tangent line at Q If this new line intersects (16 with multiplicity 3, we find the Weierstrass form using Case 1 renaming Q as P Otherwise we call the

8 8 SAMUEL A HAMBLETON other intersection of (16 and our new tangent line the projective point R Since P, Q, R cannot be collinear, the matrix is invertible Letting ( U VT UT M = ( PX Q X R X P Y Q Y R Y P Z Q Z R Z ( = M 1 XYZ and dividing the resulting equation satisfied by T,U,V by U T, we obtain an equation for which we may easily obtain the Weierstrass form by then dividing this equation by the coefficient of U 3 and then replacing T by 1 Example 31 Let C = (1,1,,1, the canonical binary cubic form of a cubic field K of discriminant 3 We fix t = 0 and n = 1 The covariant forms of C are Q = ( 5, 7,1 and F = ( 11,39,48,5 From a unit of K of norm 1 and trace 0 we have the rational point ( 1,1 of the affine elliptic curve E 1 : F(x,y = 7 In homogeneous coordinates we have P = ( 1 : 1 : 1 To find the Weierstrass form of this elliptic curve, we consider the tangent line at P, given by Y = 1 (7X +9Z Substituting into F(X,Y 7Z 3 shows that there is one other intersection at Q = ( 9 : : 3ofthe tangentlinewith theprojectiveversionofe 1 Thetangent line at Q is 155X 133Y +07Z = 0, which gives R = ( 6968 : 7569 : 931 and now we have ( M = Using this together with F(X,Y = 7Z 3, gives the affine elliptic curve in Weierstrass form, v v 4056uv = 55u u It is clear that the GL (Z class of binary cubic forms of the canonical binary cubic form C corresponding to K will give, once we fix rational numbers t,n such that (13 is non-singular, a class of isomorphic elliptic curves (13 Remark 3 suggests the following question For which elliptic curves E can we find a cubic field K and rational numbers t and n such that the canonical binary cubic form C has covariant forms Q and F which give an elliptic curve (13 defined over Q which is isomorphic to E(Q? The matrices presented in this article are also useful in understanding Voronoi s algorithm; see [1] References [1] K Belabas, A fast algorithm to compute cubic fields, Math Comp 66, (1997, no 19, MR (97m:11159 [] M Bhargava, Higher composition laws (II: On cubic analogues of Gauss composition, Ann of Math ( 159 (004, no, MR08144 (005f:1106b [3] R D Carmichael, The theory of numbers and Diophantine analysis Dover Publications, Inc, New York 1959 MR (1 #413 [4] H Davenport, H Heilbronn, On the density of discriminants of cubic fields (ii, Proc Roy Soc Lond A 3 ( MR58:10816 [5] B N Delone, D K Faddeev, The theory of Irrationalities of the third degree, Translations of Mathematical Monographs, Vol 10 American Mathematical Society, Providence, RI 1964 MR (8 #3955 [6] N Duif, Transforming a general cubic elliptic curve equation to Weierstrass form, A sage implementation 011 [7] A K Dutta, Mathematics in ancient India, I An overview, Resonance, 7, 4, (

9 A CUBIC GENERALIZATION OF BRAHMAGUPTA S IDENTITY 9 [8] A K Dutta, Mathematics in ancient India, II Diophantine equations: The Kuttaka, Resonance, 7, 10, (00 6 [9] A K Dutta, Mathematics in ancient India, III Brahmagupta s Lemma: The Samasabhavana, Resonance, 8, 11, ( [10] C Gauss, Disquisitiones Arithmeticae, Translated and with a preface by Arthur A Clarke Revised by William C Waterhouse, Cornelius Greither and A W Grootendorst and with a preface by Waterhouse Springer-Verlag, New York, 1986 MR (87f:01105 [11] S Hambleton, F Lemmermeyer, Arithmetic of Pell surfaces, Acta Arith 146 (011, no 1, 1 1 MR (01b:11097 [1] S A Hambleton, H C Williams, Cubic fields with geometry, CMS Books in Mathematics, Springer, (to appear in print [13] D Hilbert, Theory of algebraic invariants Translated from the German and with a preface by Reinhard C Laubenbacher Edited and with an introduction by Bernd Sturmfels Cambridge University Press, Cambridge, 1993 MR [14] M J Jacobson Jr, H C Williams, Solving the Pell equation CMS Books in Mathematics/Ouvrages de Mathématiques de la SMC Springer, New York, 009 [15] K J Manasa, B R Shankar, Pell surfaces and elliptic curves, J Ramanujan Math Soc 31, 1 ( [16] L J Mordell, Diophantine equations Pure and Applied Mathematics, Vol 30 Academic Press, London-New York, 1969 MR [17] E R Suryanarayan, The Brahmagupta polynomials, The Fibonacci Quarterly, 34, (1996 no 1, MR School of Mathematics and Physics, The University of Queensland, St Lucia, Queensland, Australia 407 address: sah@mathsuqeduau