Elliptic Diophantine Equations
Elliptic Diophantine Equations
and Applications 2
Editor
Colva M. Roney-Dougal, St Andrews, United Kingdom
Nikos Tzanakis
De Gruyter
Mathematics Subject Classification 2010: 11D25, 11D41, 11D88, 11G05, 11G07, 11G50,
11H06, 11J86, 11Y16, 11Y50, 11-04, 14E05, 14H52, 14Q05, 33E05, 52C07, 68W30
Author
Nikos Tzanakis
University of Crete
Department of Mathematics
Voutes Campus
70013 Heraklion, Crete
Greece
tzanakis@[Link]
ISBN 978-3-11-028091-3
e-ISBN 978-3-11-028114-9
Set-ISBN 978-3-11-028115-6
ISSN 2195-5557
Library of Congress Cataloging-in-Publication Data
A CIP catalog record for this book has been applied for at the Library of Congress.
Bibliographic information published by the Deutsche Nationalbibliothek
The Deutsche Nationalbibliothek lists this publication in the Deutsche Nationalbibliografie;
detailed bibliographic data are available in the Internet at [Link]
2013 Walter de Gruyter GmbH, Berlin/Boston
Typesetting: P T P-Berlin Protago-TEX-Production GmbH, [Link]
Printing and binding: Hubert & Co. GmbH & Co. KG, Gttingen
Printed on acid-free paper
Printed in Germany
[Link]
Preface
This book 17 0 is about elliptic Diophantine equations, the most standard instance
of such an equation being the equation y 2 D x 3 C Ax C B in integers x, y, where
A, B 2 Q and the right-hand side has no multiple roots. Many more Diophantine
equations are elliptic Diophantine equations; in Chapter 1 the term is explained in its
generality.
More specifically, the main theme of this book is the explicit resolution of such
equations and the resolution method that is exposed in detail is called elliptic logarithm method or, briefly, Ellog, in accordance with the terminology and notation
introduced in [55]. The method has two main characteristics which are, first, the exploitation of the group structure with which the points of a non-singular cubic curve
are endowed and, second, the use of linear forms in elliptic logarithms. The method
owes its name to its second characteristic, the most modern. Immediately below we
give more explanations.
The first main characteristic (or ingredient) of the Ellog, which makes possible the
transition from the elliptic Diophantine equation to linear forms in elliptic logarithms
the second main characteristic is the fact that, on the one hand, an elliptic Diophantine equation can be transformed, by an appropriate change of variables, into a
non-singular cubic equation of a special shape1 and, on the other hand, that the set of
points2 of the curve defined by this cubic equation, is endowed with the structure of
a finitely generated abelian group. The operation of addition in this group is a consequence of the simple observation that the third point of intersection with the curve
of a line joining two points of the curve with coordinates in a number field also has
coordinates in this same number field; if the two points coincide, as the line joining
them we understand the tangent at the point. Thus, if we know two distinct rational, say, solutions (points), we can obtain a third one, also with rational coordinates,
which is the third intersection with the curve of the line joining the two known points;
and if we know only one solution (point), as a line joining the two known points
we consider the tangent at the known point. This is the chord and tangent method
for generating new solutions from known ones. At this point I refer the reader to the
beautiful booklet of I. G. Bashmakova [3]. Bashmakova seems to believe that, in the
two problems below, found in Arithmetica, Diophantus applied the chord and tangent
method consciously. I am very cautious about this view; nevertheless, one might ex1
2
viii
Preface
Preface
ix
Preface
friendly) that, in practice, the checking of lattice points in the hyper-cube, mentioned
a few lines above, is impossible.
Is it possible to reduce the upper bound of M to a manageable size (say of a few
decades)? This would lead to an explicit resolution of the Diophantine problem. The
answer is, in principle, positive due to a reduction method developed by B. M. M. de
Weger [72], which is based on the LLL-basis reduction algorithm of LenstraLenstra
Lovsz [27]. Chapter 10 presents everything related to the reduction of the upper
bound of M . The reduction method is then applied to the examples of Chapter 9,
leading to upper bounds for M at most 17. This small upper bound leads then, very
easily, to the complete explicit solution.
The theoretical idea of the method8 goes back to Lang [25], who also explains it in
[26]; its brief explanation is found in [45, Chapter IX.5, Linear forms in elliptic logarithm]. The discussion so far and the contents of this book confirm that from theory
to practice a long way had to be covered; the method became practical only in 1994,
after the work of R. J. Stroeker and the author [54], and, independently, the work of
J. Gebel, A. Petho and H. G. Zimmer [15]. These two 1994 papers would not have
appeared if the work of N. Hirata-Kohno [17] and S. David [12], on lower bounds of
linear forms in elliptic logarithms, had not been previously published. In a correspondence9 during 19911992, I asked S. David if he could make explicit the constants
involved in the results of [17]; I am extremely grateful to him for his accepting my
far from non-trivial challenge, which turned out to be a really heavy work [12] of
more than 130 pages! For a clear and relatively short description of the practical 1994
method I refer to [50, Chapter XIII].
Chapter 11 is special. In it, the resolution of the Weierstrass elliptic equation in S integers is discussed. What do we mean by S -integers? If S is a finite set of primes,
an S -integer is, by definition, a rational number, with the property that the prime decomposition of its denominator allows only primes belonging to S ; in particular, every
usual (rational) integer is an S -integer.
One has to develop a theory of p-adic elliptic logarithms (with p a prime) for the
points of an elliptic curve, in analogy with the theory of elliptic logarithms developed
in Chapter 3; this is done in Section 11.1. In Section 11.2 linear forms in p-adic elliptic
logarithms are introduced and the p-adic version of Ellog, is developed. This section is
inspired by the papers of N. Smart [47] and PethoZimmerGebelHerrmann [36]; the
second paper completes the project set up in the first paper. When [36] was published,
no explicit lower bound for linear forms in p-adic elliptic logarithms the analogue in
the p-adic case of S. Davids theorem existed, except for that in [39] which, however,
is applicable only to elliptic curves of rank at most two.10 Therefore the authors of [36]
Without its details : Lower bounds for linear forms in elliptic logarithms, reduction process and,
in general, all computational aspects.
9 Handwritten letters of the good old days!
10 This result is improved in [18], but still treats the case of two p-adic elliptic logarithms.
8
Preface
xi
turned to the recent, for those days, paper [16]. Very recently,11 N. H. Kohno released
a valuable paper [19], in which the p-adic analogue of S. Davids theorem is proved.
This is what is used in Chapter 11 instead of [16]. I am grateful to Noriko, who, meeting
my desire, worked hard in order to provide me with her theorem before the date that I
had to send my manuscript to the editors.
The chapter includes a specific example with the primes p D 2, 3, 5, 7 involved.
About the style of the book. To what extent should my exposition take for granted
standard (more or less) material found in the literature? I had this speculation mainly
concerning Section 1.2 of Chapter 1, Chapters 2, 3 and 9, Section 10.1 of Chapter 10,
and Section 11.1 of Chapter 11. My decisions are detailed below.
The basic theory of elliptic curves is so beautifully written in various text-books
few of them (only) are included in my bibliography , that my hypothetical contribution could be described by C 1! Therefore, Section 1.2 of Chapter 1 includes only
the very basic facts that will be needed and gives references.
For the theory of heights, in Chapter 2, only a moderate use of p-adic theory is
required. I found that standard texts either include so much material that reference to
them would disorientate (with respect to this books aim) my reader, or they adopt
a point of view not very appropriate for the present book, as they build the relevant
theory by considering extensions of Qp .12
Concerning Weierstrass equations over C and R and the Weierstrass }-function,
treated in Chapter 3, I do the following. Since the basic general theory is so neatly
written, for example, in [1, Chapter 1], I decided that I need not provide any explanation. However, specialisation of the general theory to Weierstrass equations with
real coefficients is absolutely necessary in order to build a practical theory of elliptic logarithms. I decided to discuss this issue, rather in detail, guided by my personal
taste. I adopted a very classical point of view, mainly based on the old (still in print)
nice book [73] and the personal notes of N. Kritikos from A. Hurwitzs 19161917
E. T. H. lectures on elliptic functions.13
The hard core of Chapter 9 is a special very important though case of Sinnou
Davids Theorem. As I already mentioned, this is the result of his memoir [12] of
more than 130 pages. The theorem, in the form appropriate for the needs of this book
(Theorem 9.1.2), is stated only and aspects of its application in practice are discussed.
For the applications of Chapter 10 the main tool is the reduction technique of
B. M. M. de Weger [72], which is based on the LLL-algorithm [27]. The style of [72]
is very appropriate for this book,14 but discusses many more applications than those
11 Actually,
xii
Preface
needed here. Therefore, in Chapter 10, among other issues, I expose the reduction
process focusing on the particular applications of the book.
In Section 11.1 of Chapter 11, p-adic elliptic logarithms and their linear forms play
the fundamental role. The theory on which the construction of such logarithms is based
is described in Chapter IV of J. Silvermans valuable book [45], though from a point
of view somewhat more general than necessary for this book. What I decided to do
was state only the absolutely necessary facts from Silvermans exposition translated
into a language appropriate for practical applications.
As in the case with S. Davids Theorem in Chapter 11, the very recent and extremely
important theorem of N. Hirata-Kohno, mentioned before, is only stated in the form
appropriate for the needs of this book (Theorem 11.2.5).
Although a main characteristic of the book is its use of computational methods, it
is not a book on Computational Number Theory; issues such as to mention only a
few examples the actual computation of MordellWeil bases, the search for rational
points on elliptic curves up to a certain bound, the computation of canonical heights,
various aspects of the implementation of the LLL-algorithm, and/or improvements of
existing methods and algorithms, are beyond the scope of the book. To this rule I
allowed three exceptions: In Chapter 3, first, I did not refrain from discussing in detail
the actual computation of a fundamental pair of periods for a Weierstrass equation
with real coefficients, an issue that fits very well in the framework of the chapter.15
Second, again in Chapter 3, I did not resist the temptation to describe the very clever
algorithm of D. Zagier [74] for the computation of elliptic logarithms. Third, in Chapter 8, I present an algorithm of J. Coates related to the computation of the coefficients
of Puiseux series.
Suggestions for reading this book. The Diophantine problems treated in this book
are classified to five classes: Weierstrass, quartic elliptic, simultaneous Pell, general
elliptic, and Weierstrass in S -integers; let us use for these problems the symbols Pi ,
where i D 5, 6, 7, 8, 11, respectively, with this numbering justified by the chapter
where the corresponding problem is mainly (but not exclusively) discussed. Since
Ellog is applied in the most direct manner to P4 , I would suggest that the reader starts
by understanding the resolution of this problem. In general, in order to understand the
complete and explicit resolution of problem Pi , I suggest the following scheme:
15 Besides,
Preface
xiii
xiv
Preface
I am indebted to my wife Maro for her lifelong support, and for her warm encouragement and patience when I was writing this book; this has been a main factor for its
completion!
Heraklion, Crete, May 12, 2013
Nikos Tzanakis
Contents
Preface
1
vii
Heights
29
: E.C/ 7! C= . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
36
38
40
41
44
54
57
60
69
xvi
8
Contents
78
98
121
125
125
127
131
134
137
165
Bibliography
173
Index
177
Chapter 1
1.1
A general overview
.u, v/ 2 R R.
(1.1)
The characterisation of the above equation as Diophantine comes from the requirement that the unknowns u, v belong to the prescribed ring R and not to the whole C or
R. Solving the Diophantine equation is far different from solving the algebraic equation g.u, v/ D 0, in which the unknowns belong to C. The solutions of the algebraic
equation define a curve or, more precisely, a model C of a curve; we state this by
writing
C : g D 0; g.X , Y / D a specific polynomial in X , Y
and we say that C or, more precisely, the model C is defined by the polynomial
g.X , Y /, or by the equation g D 0. Thus, we view C as the set C.C/ D .u, v/ 2
C C : g.u, v/ D 0 and the elements of C.C/ are called points of (the model)
C . Sometimes we wish to focus our interest to the real part of C , which is the set
C.R/ D .u, v/ 2 R R : g.u, v/ D 0 of real points of (the model) C . In general, if A is a subring of C, we set C.A/ D .u, v/ 2 A A : g.u, v/ D 0 and if
.u, v/ 2 C.A/, we say that .u, v/ is an A-point of (on) C . The fact that the model C
is defined by means of the polynomial g, whose coefficients belong to R, is expressed
by saying that C is defined over R.
Sometimes (actually very rarely) we will need to refer to the projective equation or,
equivalently, to the projective model corresponding to equation (1.1). This results from
the so-called homogenisation of the variables u and v, which consists in considering
the equation
g.U : V : W / D 0,
def
(1.2)
n
g.U : V : W / D W g.U=W , V =W /,
n D maxdegu g, degv g.
Let us come back to Definition 1.1.1 and impose the further conditions that gi 2
RX , Y for i D 1, 2 and the four rational functions Uij , Vij have coefficients in
K. Then we say that C1 and C2 are birationally equivalent over K. If K D Q.R/,
the quotient field of R, then, obviously, any R-point .ui , vi / of Ci is mapped by the
birational transformation of Definition 1.1.1 to a Q.R/-point .uj , vj / of Cj .
Clearly, birational equivalence of models of curves is an equivalence relation. The
equivalence class of a model C will be called a curve, denoted by C. Thus, for this
book, a curve is an equivalence class of a model and will be denoted by a capital calligraphic letter. If C1 , C2 are models of the same curve C, then C1 and C2 share a number
of important properties, a very important one being that they have the same genus and
we thus speak about the genus of the curve C rather than the genus of this or that model.
The notion of genus is not easy to define but, fortunately, an in-depth knowledge of
this notion is not absolutely necessary for the main purpose of this book, which is the
practical resolution of elliptic Diophantine equations. Anyway, the interested reader
can refer, for example, to [45, Chapter II.5], especially Theorem 5.4. A more classical
approach is found in [66, Chapter VI, 5.3]. Another classical analytic point of view
for genus, very much in the spirit of the present book, is found in [4, Chapter III, 21];
see especially Corollary 2. A picturesque geometrical notion of genus, easy for everybody, is presented in [38, Appendix to Chapter 1]. Although this approach to genus
is not practical for the purposes of this book, we recommend the reader unaccustomed
with this notion to have a look at it.
In this book, we will deal with Elliptic Diophantine Equations. These are equations (1.1) defining a model C : g D 0 of an elliptic curve over K, which means
that the curve has genus one and the corresponding projective equation (equivalently
stated, the corresponding projective model) (1.2) has a solution (has a point) in P 2 .K/.
In our applications, K D Q but, as a tool, we will sometimes need to consider
elliptic curves defined over a number field K as well. Also, in our applications R D Z,
except for Chapter 11, in which R D ZS , the ring of S -integers, where S is a finite
set of primes; see the beginning of that chapter.
When we solve an elliptic Diophantine equation, defined by an equation (1.1), we
employ properties of the elliptic curve C, a model of which is C : g D 0. Let us put
g D g1 and C D C1 ; the method of solution that we will apply requires working with
one or two further models Ci : gi D 0 (i D 2, 3) of C, also defined over R, with any
two of C1 , C2 , C3 birationally equivalent by means of transformations Uij , Vij (cf. Definition 1.1.1) having coefficients in R (in other words, the said models are birationally
equivalent over Q.R/). These birational transformations will be fixed during the process of the resolution of the Diophantine equation g1 D 0. We will exclude all the
finitely many (and easily computed) exceptional points of C1 .C/ at which at least one
of the rational functions U1j , V1j is not defined, focusing on non-exceptional points of
C1 .C/, i.e. on points which map to points on Cj for j 1 by means of the birational
map
.u1 , v1 / ! .uj , vj / D .U1j .u1 , v1 /, V1j .u1 , v1 //
.u1 , v1 / D .Uj 1 .u1 , v1 /, Vj 1 .u1 , v1 //
.uj , vj /
For the convenience of our notation we will avoid the use of subscripts in curve notation and denote the various (two or three) birationally equivalent models that will
be used during the resolution of an elliptic Diophantine equation by, for example,
C : g.u, v/ D 0, E : f .x, y/ D 0, D : f1 .x1 , y1 / D 0 instead of Ci : gi D 0
(i D 1, 2, 3); note the use of different letters for the variables of different models. In
accordance with the chosen letters, the birational transformation e.g. from C to E will
then be denoted by .U, V/ and its inverse from E to C by .X , Y/, so that
.u, v/ ! .x, y/ D .X .u, v/, Y.u, v//
.u, v/ D .U.x, y/, V.x, y//
.x, y/.
(1.3)
(1.4)
If we denote by C the curve, models of which are C and E, then, in accordance with
what was previously said, a point P of (on) C is visualised by both a pair .u0 , y0 /
satisfying g.u0 , y0 / D 0 and a pair .x0 , y0 / satisfying f .x0 , y0 / D 0 (and analogously
if a third model of C is used), where .u0 , y0 / and .x0 , y0 / are related by (1.3) and (1.4),
with .u0 , v0 / in place of .u, v/ and .x0 , y0 / in place of .x, y/. It will be convenient to
express this by writing P C D .u.P /, v.P // D .u0 , v0 / and P E D .x.P /, y.P // D
.x0 , y0 /.
To sum up: In general, if we have to solve a certain elliptic Diophantine equation
g.u, v/ D 0, we consider the model C : g.u, v/ D 0 and the elliptic curve C corresponding to the model C , along with one or two further models E : f .x, y/ D 0,
D : f1 .x1 , y1 / D 0 which will remain fixed in the process of the resolution of the Diophantine equation. We will consider fixed birational transformations between any pair
of models C , E, D as in (1.3) and (1.4). We will say that P is a point of (on) C and
will state this by writing P 2 C, if there exist points P C D .u.P /, v.P // 2 C.C/,
P E D .x.P /, y.P // 2 E.C/ such that the rational functions U.u, v/ and V.u, v/
are defined at .u.P /, v.P // (in other words, the point .u.P /, v.P // of the model C is
non-exceptional) and the pairs .u.P /, v.P // and .x.P /, y.P // are related by means
of (1.3) and (1.4), with .u.P /, v.P // in place of .u, v/ and .x.P /, y.P // in place of
.x, y/.1
Consequently, in accordance with these conventions/notations, we have
P 2 E,
(1.5)
Finally, if, as in Chapter 3, we use only one model, say E, then we will omit the
superscript E . Moreover, and in accordance with this rule, the term point P will
mean point of the model E and we will simply write P instead of P E .
We emphasise, once again, that a deep comprehension of the notion of genus is not
absolutely necessary for the purposes of this book. From a practical point of view, the
following fact is very important.
Fact 1.1.2. If C : g.u, v/ D 0 is a model of genus one, where g has rational coefficients, then C is birationally equivalent over Q to a short Weierstrass model
E : y 2 D x 3 C Ax C B with A, B 2 Q by means of a birational transformation (1.5), the coefficients of which are real algebraic numbers of degree at most
mindegu g, degv g. These coefficients, as well as A and B, can be explicitly computed.
A practical general algorithm of M. van Hoeij for computing the short Weierstrass
model and the birational transformation, mentioned in Fact 1.1.2, is described in [22];
an implementation of this algorithm is included in the package algcurves of the computer algebra system MAPLE. Given a non-singular point .u0 , v0 / 2 C.Q/ and the coefficients of g as input, this algorithm returns the coefficients A, B of the equation
of the short Weierstrass model and the coefficients of the birational transformation
(1.3)(1.4), which belong to Q.u0 , v0 /. Although this is not very explicit in [22], it
can be verified by careful scrutiny of [20, 1 and 2.1] and [21, 13.1] on which
the algorithm of [22] is based.
The inclusion A, B 2 Q.u0 , v0 / can be further improved to A, B 2 Q: by an argument found on pages 9395 of [6], a simple transformation .x, y/ 7! .t 2 x, t 3 y/, for a
conveniently chosen t 2 Q, maps the equation y 2 D x 3 C Ax C B to an equivalent
short Weierstrass equation with coefficients in Q.
The above-mentioned algorithm, can be applied to quite exotic examples of models C : g.u, v/ D 0, like those appearing in [56]. If g has a more ordinary form,
one can see more directly the necessary birational transformation between C and
the short Weierstrass model, or even do the necessary computations by hand; see
[10, Section 7.2].
1.2
Indicated titles; only a few with which the author happens to be more acquainted.
theory textbooks devoted to elliptic curves, e.g. [7, Section 5.9], [32, Section 5.5]3 . In
this chapter we will mostly refer to [45, Chapter III].
Firstly, any elliptic curve E over K has a model
E : y 2 C a1 xy C a3 y D x 3 C a2 x 2 C a4 x C a6 ,
(1.6)
(1.7)
By the previous discussion we can form the sum Q1E C Q2E , which is a point, say
P E , of the model E. If the birational transformation (1.4) is also defined on P E , then,
in accordance with the conventions adopted in Section 1.1, P E is the representative
of a point P of E. Then, this will allow us to write Q1 C Q2 D P , i.e. by definition,
def
Q1 C Q2 D P Q1E C Q2E D P E .
If, moreover, the birational transformation (1.5) is defined on Q1E , Q2E and P E , then
the following makes sense:
def
m1 Q1 C C ms Qs D P m1 Q1E C C ms QsE D P E .
(1.8)
and
.m1 Q1 C C ms Qs /C D .U..m1 Q1E C C ms QsE /, V..m1 Q1E C C ms QsE //.
(1.9)
Finally, we give a simple useful lemma, that we will apply in subsequent chapters.
Lemma 1.2.2. If the coefficients a1 , a2 , a3 , a4 , a6 in (1.6) are in Z and .x, y/ 2 E.Q/,
then .x, y/ D .x1 =z 2 , y1 =z 3 /, where x1 , y1 , z 2 Z, z > 0 and gcd.x1 , z/ D 1 D
gcd.y1 , z/.
Proof. Let us write .x, y/ D .x1 =w1 , y1 =w2 /, where x1 , w1 , y1 , w2 are integers with
w1 , w2 positive and gcd.x1 , w1 / D 1 D gcd.y1 , w2 /. Replacing for x and y in (1.6)
we get
w12 .y12 w1 C a1 x1 y1 w2 C a3 y1 w1 w2 / D w22 .x13 C a2 x12 w1 C a4 x1 w12 C a6 w13 /.
We note that w12 is relatively prime to the number in the parenthesis in the right-hand
side, therefore w1 jw2 and we write w2 D zw1 , where z is a positive integer relatively
prime to x1 y1 . Replacing for w2 in the above equation we obtain
w1 .y12 C a1 x1 y1 z C a3 y1 w1 z/ D z 2 .x13 C a2 x12 w1 C a4 x1 w12 C a6 w13 /.
Now we observe that w1 is relatively prime to the number in the parenthesis in the
right-hand side, hence w1 jz 2 , and, on the other hand, z 2 is relatively prime to the
number in the parenthesis in the left-hand side, hence z 2 jw1 . It follows that w1 D z 2
and w2 D zw1 D z 3 , as claimed.
Chapter 2
Heights
In all Diophantine applications of the theory of elliptic curves we need a measure for
the points of the elliptic curves involved in our investigations. This is accomplished
with the notion of height of a point.1 In this chapter we develop the relevant theory
from a point of view appropriate for the purposes of this book.
2.1
We assume familiarity with standard algebraic number theory, as exposed, for example, in [35, Chapters 4 3, 5, 6] and (partially) in [40, Chapter V]. For the convenience
of the reader, we collect the basic facts that we will need without giving special reference to each of them.
We will need to work with various relative finite extensions L=K=Q, where K and
L may vary. The ring of integers of K will be denoted by ZK and analogously for L.
For the elements of the number fields we will use small Greek letters.
For the ideals of ZK we will use small fracture letters (like p, for example) and for
the ideals of ZL we will use capital fracture letters (like P, for example). The symbol
j between ideals means a divisibility relation, and is equivalent to the relation .
If a is a non-zero ideal of ZK , then ZK =a is finite and we write NK=Q .a/ :D
[Link] =a/; this is the absolute norm of the ideal a. If 2 K , then NK=Q .hi/ D
jNK=Q ./j, where NK=Q ./ is the usual element norm of .
By default, the term ideal refers to integral ideal, i.e. ideal of the ring of integers.
But we will also use fractional ideals, for which we will use fracture letters as well; the
fact that an ideal is not necessarily integral will be emphasised by using the adjective
fractional. If a is an ideal (integral or fractional) in K and we want to view it as an
ideal in L, we will write aZL . If 2 K, we will write hi for the principal ideal ZK .
If we view as an element of L rather than as an element of K, then we will write
ZL for the principal ideal that generates in L.
We will write P for the set of rational primes and PK , PL for the set of prime ideals
of ZK and the set of prime ideals of ZL , respectively. If a is a non-zero fractional
ideal in K and p 2 PK , we define the exponent of p in a and denoted by p .a/ as the
exponent of p in the factorisation of a into prime ideals. If 2 K , we will write p ./
instead of p .hi/. Analogous notations will be used for the non-zero ideals of L.
1
10
Chapter 2 Heights
For P 2 PL there exists precisely one p 2 PK , such that PjpZK (the ideal P is
over p or, in equivalent wording, p is below P). This p, in turn, is over a prime p 2 P,
i.e. p is the unique rational prime divisible by p.
With p, p and P related as just above, let us view p 2 PK as an ideal (not necessarily
prime) of ZL , i.e. let us consider the ideal pZL of ZL and decompose it into prime
ideals, as follows:
Y
pZL D
PeL=K .P/ ;
(2.1)
Pjp
where eL=K .P/ is the ramification index of P relatively to L=K (hence eL=K .P/ D
P .pZL ). The following relation is valid if P 2 PL is over p 2 PK :
eL=Q .P/ D eL=K .P/ eK=Q .p/.
(2.2)
Further, the finite field ZL =P is a finite extension of the finite field ZK =p. The degree
of this extension, which is called the degree of P relative to the extension L=K, is
denoted by fL=K .P/. Similarly to (2.2) we have the relation
fL=Q .P/ D fL=K .P/ fK=Q .p/.
(2.3)
Then the ideal norm of P relative to the extension L=K is, by definition,
NL=K .P/ D pfL=K .P/ ;
(2.4)
note that, if in place of L and K we respectively have K and Q, this definition agrees
with the one given above for the absolute norm. The definition of the ideal norm is
extended by multiplicativity to all non-zero fractional ideals of L. Namely, if A is a
non-zero fractional ideal in L and
Y
PP .A/
AD
P2PL
is its prime decomposition, then, by definition, the ideal norm of A relative to the
extension L=K is
Y Y
def
pP .A/fL=K .P/
(2.5)
NL=K .A/ D
p2PK Pjp
and the norm of the zero ideal is defined to be the zero ideal. Two important properties
of the norm are the following: The function NL=K from the group of fractional ideals
of L to the group of fractional ideals of K is a homomorphism; moreover, for any
fractional ideal A of L we have
NL=Q .A/ D NK=Q .NL=K .A//.
(2.6)
If we have the decomposition (2.1), then, for every P dividing p we have the important relation
X
def
dL=K .P/ D L : K, dL=K .P/ D eL=K .P/ fL=K .P/;
(2.7)
Pjp
11
dL=K .P/ is called local degree of L=K at P. Using (2.2), (2.3) and (2.7) we obtain
another useful relation:
X
X
dL=Q .P/ D
eL=Q .P/fL=Q .P/
Pjp
Pjp
Pjp
(2.8)
Pjp
2.2
Easily then, every absolute value j j equivalent to j jp has the same property that characterises nonarchimedean absolute values: j C j maxjj, jj and, if jj jj, then strict inequality holds.
Therefore it makes sense to speak about a non-archimedean place.
12
Chapter 2 Heights
complex
p2PK
1 if is real
d D
2 if is complex.
The proof of the product formula follows immediately from
Y
Y
Y
d
.p/
jjp K=Q D
NK=Q .p/p ./ D NK=Q
pp ./
where
p2PK
p2PK
p2PK
jjd D
o
2EK
(2.11)
2EK
The following lemma is very useful, as it relates the absolute values in relative extensions.
Lemma 2.2.1. Let L=K be an extension of number fields and p 2 PK . Then,
(1) for any 2 K and any P dividing p,
jjP D jjp ,
(2) for any 2 K,
L:KdK=Q .p/
jjp
Y
Pjp
jjPL=Q
.P/
13
Proof. Obviously, it suffices to consider a non-zero 2 K. With the aid of (2.2), (2.3)
and the definitions of dL=Q .P/ and dK=Q .p/ we observe that
dL=Q .P/ D eL=Q .P/fL=Q .P/ D eK=Q .p/eL=K .P/ fK=Q .p/fL=K .P/
D eK=Q .p/fK=Q .p/ eL=K .P/fL=K .P/ D dK=Q .p/ eL=K .P/fL=K .P/
which shows that
eL=K .P/fL=K .P/ D
dL=Q .P/
.
dK=Q .p/
jjPL=Q
.P/
D jjp L=Q
It follows then that
Y
.P/
jjPL=Q
.P/
D jjp
Pjp
dL=Q .P/
L:KdK=Q .p/
D jjp
Pjp
where the right-most equality holds because of (2.8). The proof is now complete.
2.3
14
Chapter 2 Heights
(2.13)
maxjxi jP dP D
P2PL
Y Y
maxjxi jP dP .
p2PK Pjp
By (1) of Lemma 2.2.1, if Pjp, then jxi jP D jxi jp , therefore, in the right-most side of
the above displayed equation, maxi jxi jP D maxi jxi jp D (say) jxi0 jp . Then
HL,fin .x0 : x1 : : xn / D
Y Y
dP
jxi0 jp
D jxi0 jp
p2PK Pjp
Y
Pjp dP
L:K
L:Kdp
D jxi0 jp
(by (2.8))
maxjxi jp dp
p2PK
(2.14)
15
Y Y
2EL
maxj .xi /j
2EK j
Y
L:K
maxj .xi /j
2EK
Y
L:K
d
maxj .xi /j
o
2EK
(2.15)
(2.16)
h. n / D jnjh./.
(2.17)
16
Chapter 2 Heights
Proof. Let K be any number field containing and . Denote by P the set of prime
ideals of ZK and by EK the set of all embeddings K ,! C; note that [Link] / D K :
Q. For every p 2 P, let dp D dK=Q .p/.
We have
Y
Y
H./K:Q D HK ./ D
max1, jjp d p
max1, j ./j.
p2P
2EK
Using the inequality max1, xy max1, x max1, y, valid for every pair
of non-negative real numbers x, y, we obtain from the above displayed inequality
H./K:Q H./K:Q H./K:Q , from which the inequality (2.15) follows immediately on taking logarithms.
In order to bound H. C /K:Q we first observe the following. If p 2 P, then
max1, j C jp max1, maxjjp , jjp max1, jjp max1, jjp .
For 2 EK ,
max1, j . C /j max1, j ./j C j ./j
2 max1, j ./j max1, j ./j.
Note that, if in the last inequality we let run through EK and we multiply the resulting
relations, then the factor 2K:Q will appear in the right-hand side.
Using the above two displayed inequalities in HK,fin . C / and HK,1 . C /,
respectively, we get
H. C /K:Q HK,fin ./ HK,fin ./ 2K:Q HK,1 ./ HK,1 ./
D 2K:Q HK ./HK ./ D .2H./H.//K:Q ,
which completes the proof of (2.16).
Concerning (2.17), this is clear if n 0, so it suffices to prove that h. 1 / D h./.
We use the identity max1, x 1 D x 1 max1, x, valid for every real x > 0. We
have
Y
Y
max1, j 1 jp d p
max1, j . 1 /j
H. 1 /K:Q D HK . 1 / D
D
p2P
1
jp max1, jjp
p2P
Y
p2P
dp
j 1 jp
j 1 j
2EK
j . 1 /j max1, j ./j
2EK
HK ./ D H./K:Q ,
2EK
where for the right-most equality we used the product formula (2.10). Thus H./ D
H. 1 /, as required.
17
Heights over Q. Let us specialise our previous discussion to the important case of
rational numbers.
First, let us note that, for a point .x0 : x1 : : xn / 2 P n .Q/, its absolute height is
Y
H.x0 : x1 : : xn / D maxjx0 j, jx1 j, : : : , jxn j
maxjx0 jp , jx1 jp , : : : , jxn jp ,
p2P
(2.18)
bja1 j
bjan j
,:::,
D bmax1, jx1 j, : : : , jxn j. (2.19)
b1
bn
maxjaj, jbj
.
gcd.a, b/
(2.21)
A third expression for h.q/ results as follows. Let q D a=b, where a, b 2 Z and write
q as a fraction a1 =b1 in lowest terms. Then, specialising to n D 1 the right-most side
of (2.19), we have H.1 : q/ D jb1 j max1, jqj and, since b1 D b= gcd.a, b/, we
conclude that
jbj
h.q/ D log max1, jqj C log
.
(2.22)
gcd.a, b/
18
Chapter 2 Heights
ci D
aj bk
.i 0/.
(2.23)
j CkDi
(2.24)
p2PK
where the right-most equality holds due to the product formula (2.10) and (2.11).
The following proposition is, actually, a generalised and more sophisticated version
of the so-called Gauss Lemma for Polynomials.
Proposition 2.4.1. If P , Q 2 KX , then HK,fin .PQ/ D HK,fin .P /HK,fin .Q/.
Proof. Obviously it suffices to prove our claim when both P and Q are non-zero
polynomials. Let P D .ai /, Q D .bi / and PQ D .ci / (cf. (2.23)). Fix p 2 PK . Let
P D .aj0 /, bk1
Q D .bk0 /.
jaj0 jp D maxj jaj jp , jbk0 jp D maxk jbk jp and set aj1
0
0
Obviously, jaj0 jp 1 for every j and jbk0 jp 1 for every k, with equality holding for
at least one subscript in both cases.
P
Now aj1
bk1
PQ D .ci0 /, where ci0 D j CkDi aj0 bk0 for every i , so that jci0 jp 1.
0
0
Let j1 and k1 be respectively the least j for which jaj0 jp D 1 and the least k for
P
which jbk0 jp D 1. We have cj01 Ck1 D j CkDj1 Ck1 aj0 bk0 . Consider .j , k/ such that
j C k D j1 C k1 . If j < j1 , then jaj0 jp < 1 and, since jbk0 jp 1, this implies
19
jaj0 bk0 jp < 1. If j > j1 , then k < k1 , so that jbk0 jp < 1 and jaj0 bk0 jp < 1. It follows
P
that, in ci0 D j CkDi ai0 bj0 , all summands but aj0 1 bk0 1 have p-absolute value strictly
less than 1, while jaj0 1 bk0 1 jp D 1, leading to the conclusion that jcj01 Ck1 jp D 1. But,
since jci0 jp 1 for every i , we conclude that maxi jci0 jp D 1 so that,
d
Since the above relation holds for every p 2 PK , the proof is now complete.
Proposition 2.4.2. Let be an algebraic number and f .X / D a0 X d C C ad
its minimal polynomial over Z. By this we mean that f has integer coefficients, with
a0 > 0 and gcd.a0 , : : : , ad / D 1; f is irreducible over Q, and f ./ D 0. Denote
by .i/ , i D 1, : : : , d the roots of f in C. Then,
Y
d
1
.i/
max1,
j
j
.
h./ D log a0
d
iD1
1
1
log HK .1 : / D [Link],fin ./ HK,1 .//
d
d
1
log.a0 HK,1 .//.
d
(2.25)
20
Chapter 2 Heights
A, B 2 Q.
(2.26)
g2 D 4A, g3 D 4B,
.x 0 , y 0 / D .x, 2y/.
For points P of E, such that P E 2 E.Q/, the naive or Weil height of P is, by
definition,
def
h.P / D h.x.P //.
In the special important case that P E 2 E.Q/, we put P E D .x, y/, so that P D D
.x 0 , y 0 / D .x, 2y/, with x D a=b, a, b 2 Z, gcd.a, b/ D 1. Then,
X
def
log max1, jxjp D log maxjaj, jbj. (2.27)
h.P / D h.x/ D log max1, jxj C
p
For S. David, whose main result in [12] is extremely important for the application of
Ellog, it is more convenient to use a somewhat different height of P , namely,
hD .P / D h.1 : x 0 : y 0 / D h.1 : x : 2y/.
def
(2.28)
From the point of view of Diophantine equations, it is more convenient to work with
the model E than with the model D, this last model being more appropriate for [12].
Therefore we need to see how the two heights are related. This we do in the following
lemma.
Lemma 2.5.1. Let E as in (2.26) and let P E D .x, y/ 2 E.K/, where K is a number
field. Then
1
5
hD .P / h.P / C 2 log 2 C .h.A/ C h.B//.
2
2
21
Proof. We will make extensive use of Proposition 2.3.2, without special reference
each time we use any of the relations (2.15), (2.16) or (2.17). Moreover, we will make
use of the fact that, if x1 , x2 are algebraic numbers, then
h.1 : x1 : x2 / h.x1 / C h.x2 /
(2.29)
22
Chapter 2 Heights
and
1
max1, jxj, jyj < max1, jxj, jyj max1, jxj, j2yj 2 max1, jxj, jyj,
2
from which it follows that
log 2 C hE .P / hD .P / log 2 C hE .P /.
(2.30)
(2.31)
(2.32)
23
By (2.18),
H..dt /3 : [Link] / : s/ D
(2.33)
where
:D gcd..dt /3 , dt r, s/ gcd.d 3 t 3 , d 3 t r, d 3 s/ d 3 gcd.t 3 , t r, s/ D d 3
and we have used the fact that gcd.r, t / D 1 D gcd.s, t /.
By (2.33) and (2.31),
C1
.dt /3
.dt /3
max1, jxj3=2 H..dt /3 : [Link] / : s/ C2
max1, jxj3=2 . (2.34)
d 2t 2
d 2t 2
D
log
max1,
jxj
C
log
,
gcd.r, d 2 t 2 /
gcd.r, d 2 /
hence
max1, jxj3=2 D
gcd.r, d 2 /3=2
exp. 32 h.x//.
.dt /3
where, obviously,
1
gcd.r, d 2 /3=2
d 3.
d3
2.6
Let now E be an elliptic curve over a number field K and P 2 E. By Fact 1.1.2 of
Section 1.1 there exists a Weierstrass model
C :
y 2 C a1 xy C a3 y D x 3 C a2 x 2 C a4 x C a6
with coefficients in K. Working with the model C of E and its function field K.C /, we
24
Chapter 2 Heights
(2.35)
where h./ denotes the absolute logarithmic height. This limit is independent from the
Weierstrass model. Indeed, let C 0 : W .x 0 , y 0 / D 0 be another Weierstrass model over
K of E. By [45, Proposition 3.1(b), Chapter III.3], x 0 .P / D
2 x.P /C for convenient
, 2 K. Working now with the model C 0 of E and its function field K.C 0 /, we have
0
the degree-two function x 0 2 K.C 0 /, defined by C.K/ 3 P C D .x 0 .P /, y 0 .P // 7!
x 0 .P / 2 K and this function is even (x 0 .P / D x 0 .P /). Therefore, as before, the
0
canonical height of P C can be defined by
0
h.x 0 .2N P C //
1
.
lim
2 N !1
4N
But x 0 is an even, degree-two function of K.C / also, in view of the simple relation
between x 0 .P / and x.P /, therefore, by [45, Proposition 9.1, Chapter VIII.9], the two
limits displayed above are equal and their common value (2.35) is defined as the canonO /, without any model
ical or NronTate height of the point P 2 E, denoted by h.P
indication on P . Thus, taking also into account (2.27), we have
N C
O / D 1 lim h.2 P / .
(2.36)
h.P
2 N !1
4N
It is worthwhile to note, although we will not need it in this book, that for the definition
of the canonical height it is not necessary to use the function x 2 E.K/, but any even
function in E.Q/; see [45, VIII.9].
As already noted, for S. David it is more convenient to work in [12] with models D
like the one on page 20 and adopt hD as the absolute logarithmic height of a point, as
defined in (2.28). As a consequence, in the notation just before the relation (2.28), the
canonical height used by S. David is defined by
hD .2N P /
.
(2.37)
N !C1
4N
As an immediate consequence of the above relation and Lemma 2.5.2 we have proved
the following proposition:
hOD .P / D
lim
g2 D 4A, g3 D 4B,
25
O /
so that D is also a model of E. Then, for any point P 2 E, the canonical heights h.P
O
and hD .P /, defined by (2.36) and (2.37), respectively, are related by
O /.
hOD .P / D 3h.P
A very important property of the canonical height is that, by means of it, a positivedefinite quadratic form is defined as follows: First, one defines the so-called Nron
Tate (or Weil) pairing by
O C Q/ h.P
O / h.Q/.
O
hP , Qi D h.P
The following important properties for the canonical height and the NronTate pairing hold (see [45, Theorem 9.3]):
Using these properties, it is not difficult to see that, if P E is expressed as in (1.7), then
X
O /D 1
h.P
hPi , Pj imi mj .
(2.38)
2
1i,j r
The function hO can be extended to the r-dimensional real vector space E.Q/ R,
where E.Q/=Etors .Q/ sits as an r-dimensional lattice, and defines a positive definite
quadratic form on this vector space; see [45, items 9.4-9.7]. Then, the height pairing
matrix of .P1 , : : : , Pr / is, by definition, the positive definite matrix
H D H.P1 , : : : , Pr / D . 12 hPi , Pj i/rr
(2.39)
(cf. (2.38)).
The following result is crucial for the applications of this book:
Proposition 2.6.2. Let P E be expressed as in (1.7). Then
O / max m2 ,
h.P
i
1ir
where is the least eigenvalue of the height pairing matrix H defined in (2.39).
Proof. According to (2.38) we have
O / D mT Hm,
h.P
where m is the column vector with components m1 , : : : , mr . As H is symmetric, a
def
diagonal matrix of eigenvalues 0 < D 1 < 2 < < r of H and an
26
Chapter 2 Heights
r
X
iD1
r
X
iD1
r
X
i n2i
iD1
n2i D nT n D mT QT Qm D mT m
m2i c1 max m2i ,
1ir
as claimed.
Proposition 2.6.3. Let E : y 2 D x 3 C Ax C B with A, B 2 Q be an elliptic curve
model. Let
D : y12 C a1 x1 y1 C a3 y1 D x13 C a2 x12 C a4 x1 C a6 ,
(2.40)
1
1
1
1
log jj C
logC jj j C logC jb2 =12j C log 2
12
12
2
2
1
C .log 2 C h.// C h.
/ C log j
j C 1.07,
2
(2.41)
where and j are, respectively, the discriminant and j -invariant of the model E,
b2 D a12 C 4a2 , 2 D 1 or 2 according to whether b2 vanishes or not, respectively,
and logC is defined for any real > 0 by logC D log max1, .
Then, for every P E D .x.P /, y.P // 2 E.Q/, we have
O / 1 h.x.P // .
h.P
2
(2.42)
Proof. Note that the two models D and E have equal j -invariants, while their corresponding discriminants 1 and are related by
1 D
12 .
(2.43)
Now we apply Silvermans Theorem 1.1 in [44] to the model (2.40). That theorem
requires that the ai s be algebraic integers. In our proposition, these coefficients are
3
One can easily see that there always exists such a model D.
27
rational integers, which permits to replace h./ and h1 .j /, h1 .b2 =12/, which appear in Silvermans theorem, by log jj and logC jj j, logC jb2 =12j, respectively. We
thus obtain the following inequality
O / 1 h.x1 .P // 0
h.P
2
(2.44)
with
0 D
1
12
log j1 j C
1
12
where 1 and j1 are, respectively, the discriminant and the j -invariant of the model
D. As already noted, (2.43) holds and j1 D j , the j -invariant of E. Therefore,
0 D
1
12
1
12
Now,
h.x1 .P // D h.
2 x.P / C /
log 2 C 2h.
/ C h./ C h.x.P //
therefore,
O / 1 h.x.P // h.P
O / 1 h.x1 .P // C 1 .log 2 C h.// C h.
/
h.P
2
2
2
1
0
(by (2.44)).
C .log 2 C h.// C h.
/
2
On combining the last inequality with (2.45) we obtain the inequality (2.42).
Remark. According to the announcement of Proposition 2.6.3, the Weierstrass model
D must fulfil certain conditions, but otherwise it is arbitrary. What we do in practice
when we have to solve a specific Diophantine equation, is to choose (2.40) in such
a way that the value of is as small as possible. Usually we choose for (2.40) the
(global) minimal Weierstrass model (see [45, Section VIII.8]), which can be computed
by various packages, for example, PARI, MAGMA, SAGE; not rarely, it happens that the
minimal Weierstrass model is (4.1).
The following proposition is useful in order to compute an upper bound for the
canonical height for the point of an elliptic curve, if the coordinates of the point belong
to a number field. Computing the canonical height of such a point is more difficult
compared with the analogous task for a rational point (see [43]); anyway, standard
packages such as PARI and MAGMA refuse to work with points over a number field.
Fortunately, for the needs of this book it suffices that we know a reasonably good
upper bound for the canonical height.
28
Chapter 2 Heights
1
12
log jj C
1
12
where, 1 and j1 are, respectively, the discriminant and j -invariant of the model D,
b2 D a12 C 4a2 , 2 D 1 or 2 according to whether b2 vanishes or not, respectively,
and logC is defined for any real > 0 by logC D log max1, .
Proof. Immediate application of Silvermans theorem [44, Theorem 1.1]. Note that,
in our proposition, the coefficients of the equation defining D are rational integers.4 Therefore, 1 2 Z and j1 , b2 2 Q, so that we can replace h.1 / and
h1 .j1 /, h1 .b2 =12/, appearing in Silvermans theorem, by log j1 j and logC jj1 j,
logC jb2 =12j, respectively.
Chapter 3
3.1
For the basic background of this section we refer to [1]. Let !1 , !2 be two complex
numbers, such that !2 =!1 62 R, which we will call periods and let D Z!1 CZ!2
C be the lattice generated by these numbers, which we call a period lattice. The paral
!
lelogram D x1 !1 C x2 !2 : 0 x1 , x2 < 1 and all its translations by vectors 0!
1
with ! 2 are called period parallelograms of . The lattice has infinitely many
period parallelograms, since there are infinitely many bases .!1 , !2 / of ; actually,
any pair .!10 , !20 / obtained from .!1 , !2 / by a unimodular linear transformation with
integer coefficients is another basis for .
Definition 3.1.1. The Weierstrass }-function corresponding to the lattice is defined
by the series
X
}.z/ D z 2 C
..z !/2 ! 2 /.
(3.1)
!2X0
The Weierstrass }-function defined above is an even, doubly periodic function, with
set of periods . It has a double pole at each ! 2 and is analytic at every z 2 C X;
see [1, Theorem 1.10 ]. Since } is a meromorphic and doubly periodic function, it is,
by definition, an elliptic function. Any basis .!1 , !2 / of is also called a fundamental pair of periods of }. If we translate the closure of a period parallelogram of by
!
a vector 0z, (the end point z need not be a lattice point) and from each pair of parallel sides of the resulting (closed) parallelogram we remove exactly one side (and its
vertices), we obtain a fundamental parallelogram of }. Thus, from the four vertices
which bound a fundamental parallelogram, exactly one is contained in it. Obviously,
any period parallelogram of is a fundamental parallelogram of the corresponding
}-function but, clearly, the converse is not true.
As noted just before Definition 3.1.1, for the function }, there are infinitely many
fundamental pairs of periods, hence there are infinitely many period parallelograms
and fundamental parallelograms. In view of periodicity it is clear that } may be viewed
as an 1-1 function on any of its fundamental parallelograms rather than as a function
on C.
Note that, in a period parallelogram, only the left-lower vertex is included; equally well, however,
instead of this vertex, one could include a different one.
30
Gn D
!2X0
1
;
!n
(3.3)
see [1, Theorem 1.11]. The function } satisfies the differential equation
} 0 .z/2 D 4}.z/3 g2 }.z/ g3 ,
(3.4)
where, by definition,
g2 D 60G4 ,
g3 D 140G6 ;
(3.5)
X
!2
1
,
.z !/3
Qi D .ei , 0/ 2 E.C/,
.i D 1, 2, 3/
(3.7)
The roots (3.6) are distinct, hence the discriminant of the cubic polynomial 4X 3
g2 X g3 is non-zero: D g23 27g32 0; see [1, 1.10]. By [1, Theorem 1.13],
the coefficients G2nC2 in (3.2) are expressible as polynomials in g2 , g3 with positive
rational coefficients, therefore the Weierstrass } function is uniquely defined by the
pair .g2 , g3 / satisfying the condition g23 27g32 0. The converse is also true. More
precisely, if two given complex numbers a2 , a3 given satisfy a23 27a32 0, then
there exists a lattice D Z!1 C Z!2 , with !2 =!1 62 R, such that g2 ./ D a2 and
g3 ./ D a3 ; see [1, Theorem 2.9].
31
We summarise:
Fact 3.1.2. Given a lattice D Z!1 C Z!2 , with !1 =!2 62 R, we define the Weierstrass function }.z; / D }.z/ by means of (3.1), and the parameters g2 ./ D g2
def
and g3 ./ D g3 by means of (3.5) and (3.3). Then, the discriminant ./ D D
g23 27g32 is non-zero and the differential equation (3.4) is satisfied.
Conversely, given two complex numbers a2 , a3 such that a23 27a32 0, there exists
a lattice D Z!1 C Z!2 , with !1 =!2 62 R, such that g2 ./ D a2 and g3 ./ D a3 .
The Weierstrass function }.z; / D }.z/ defined by (3.1) will be called Weierstrass
} function with parameters g2 D a2 and g3 D a3 .
It is sometimes convenient to replace the lattice by its homothetic lattice D
Z C Z , where D !2 =!1 , so that D !1 . The simple equations below relate the
values of g2 , g3 and that correspond to the lattices and . Note, however, that,
instead of writing g2 . /, we simply write g2 . / and similarly for g3 and . We have
g2 D g2 ./ D !14 g2 . /
g3 D g3 ./ D !16 g3 . /
D ./ D !112 . /
(3.8)
(see, for example, [1, 1.11]). Moreover, if is in the upper half-plane, then there is
a Fourier expansion for . /
. / D .2/
12
1
X
t .n/e 2 i n ,
nD1
1
Y
.1 e 2 i n /
(3.9)
(3.10)
nD1
3.2
A, B 2 C, 4A3 C 27B 2 0
(3.11)
32
!1 !1 C!2 !2
,
,
2
2
2
Section 3.3
: E.C/ 7! C=
33
period parallelogram) and let ! be the only pole of } in F. Then the map
.}.z/, 12 } 0 .z// if z !
2 E.C/
F 3 z 7!
O
if z D !
(3.12)
3.3
: E.C/ 7! C=
In this section E will be the elliptic curve model considered in Fact 3.2.1. According
to Section 1.1, E is just one out of many models of an elliptic curve E, therefore, any
point on E, being a representative of a certain abstract point P (say) of E, should be
denoted by P E . However, since in the present chapter we will not make use of any
other model of E except for E, and for the sake of simplifying the notation, we will
omit the superscript E on P , simply writing P instead of P E .
Our main purpose is to define a group isomorphism
the inverse of that we defined in Fact 3.2.1.
34
(a) }.z1 / }.z2 /. Then,
y.P2 / y.P1 / 2
x.P1 / x.P2 /
x.P1 C P2 / D
x.P2 / x.P1 /
1 } 0 .r/ } 0 .z2 / 2
}.z1 / }.z2 /
D
4 }.z1 / }.z2 /
D }.z1 C z2 /,
the last equality being true because of the addition-theorem [73, 20.31]. On the
other hand, by the first determinant formula in [73, 20.3], we have
}.z1 /
} 0 .z1 /
1
}.z2 /
} 0 .z2 /
1 D 0
}.z1 C z2 / } 0 .z1 C z2 / 1
and after some standard calculations we find that
1 } 0 .z2 / } 0 .z1 /
1 } 0 .z1 /}.z2 / } 0 .z2 /}.z1 /
1 0
} .z1 C z2 / D
}.z1 C z2 /
2
2 }.z2 / }.z1 /
2
}.z2 / }.z1 /
y.P2 / y.P1 /
y.P1 /x.P2 / y.P2 /x.P1 /
D
x.P1 C P2 /
x.P2 / x.P1 /
x.P2 / x.P1 /
D y.P1 C P2 /.
This proves that .x.P1 C P2 /, y.P1 C P2 // D .}.z1 C z2 /, 12 } 0 .z1 C z2 //, hence
.P1 C P2 / D .z1 C z2 / C D .P1 / C .P2 /.
(b) }.z1 / D }.z2 /. Then, } 0 .z2 / D } 0 .z1 /. If the minus sign holds, then
.}.z2 /, 12 } 0 .z2 // D .}.z1 /, 12 } 0 .z1 //, hence P2 D P1 and .P1 CP2 / D .O/ D
. On the other hand, .P1 / C .P2 / D .z1 C z2 / C D and this last equality
holds because .}.z2 /, 12 } 0 .z2 // D .}.z2 /, 12 } 0 .z2 // D .}.z1 /, 12 } 0 .z1 // which
implies that z1 D z2 .mod /.
If the plus sign holds, then z1 C D z2 C , hence, P1 D P2 D (say) P and
we put .x.P /, y.P // D .}.r/, 12 } 0 .r// with r 2 P . We have then to show that
.2P / D 2 .P /. If P D .ei , 0/ for some i 2 1, 2, 3 (cf. (3.6)), then, on the
2 !2
one hand, z D z0 2 !21 , !1 C!
, 2 and, on the other hand, 2P D O. Consequently,
2
.2P / D D 2.z0 C/ D 2 .P /. If for every i D 1, 2, 3 we have P .ei , 0/, then
y.P / 0 and } 0 .r/ 0. Then, differentiation of } 0 .r/2 D 4}.r/3 g2 }.r/ g3 D
4}.r/3 C 4A}.r/ C B gives } 00 .r/ D 2.3}.r/2 C A/. By the duplication formula
[73, 20.311],
1 } 00 .r/ 2
3}.r/2 C A 2
}.2r/ D
2}.r/
/ 2}.r/ D
4 } 0 .r/
} 0 .r/
.3x.P /2 C A/2 8xy 2
D x.2P /.
D
4y 2
Section 3.3
: E.C/ 7! C=
35
CA 2
We differentiate the relation }.2r/ D . 3}.r/
/ 2}.r/ above and take into
} 0 .r/
00
2
account the relations } .r/ D 2.3}.r/ C A/, }.r/ D x.P / and } 0 .r/ D
2y.P /. After some elementary calculations we find that } 0 .2r/ D 2y.2P /. Thus,
.x.2P /, y.2P // D .}.2r/, 12 } 0 .2r//, which shows that .2P / D 2r C D 2 .P /,
as required.
Combining our discussion above with Fact 3.2.1 we obtain the following theorem:
2
Theorem 3.3.1. We refer to Fact 3.2.1, in which we put F D , where is the period
parallelogram at the beginning of this section. We identify F D with the group
C=. Then : C= ! E.C/ is a group isomorphism, whose inverse isomorphism
: E.C/ ! C= is as follows: .O/ D and, for a non-zero point P D .x, y/ 2
E.C/, we define .P / D r C , where r is the unique non-zero r 2 with the
property .}.r/, 12 } 0 .r// D .x, y/.
Definition 3.3.2. Fix any fundamental parallelogram F of }. Theorem 3.3.1 implies
that, for any P D .x, y/ 2 E.C/ there exists a unique r 2 F such that P D
.}.r/, 12 } 0 .r//. This r is called elliptic logarithm of the point P (belonging to F).
In analogy with complex logarithms, any point has infinitely many elliptic logarithms, depending on the fundamental parallelogram which we choose. Any two of
them, however, differ by an element of . This situation is similar to that of complex
logarithms, in which any non-zero complex number has infinitely many logarithms
and any two of them differ by an integral multiple of 2 i .
Now a natural problem arises:
Problem 3.3.3. Given a point P D .x.P /, y.P // 2 E.C/ and a fundamental parallelogram F, compute the elliptic logarithm of P belonging to F. In other words,
compute r 2 F such that .}.r/, 12 } 0 .r// D .x.P /, y.P //.
Once again we turn to the classical book [73]. First, we cut the complex plane along
a line segment joining any two of e1 , e2 , e3 and along any half-line with origin the
remaining ei which (half-line) has no common points with the (closed) line segment.
What remains then is a cut plane which we denote by C ; we understand that e1 , e2 , e3
are not points of C . If we fix some z0 2 C, such that 4z03 g2 z0 g3 0, along with
one of the two square roots of 4z03 g2 z0 g3 , let us denote it by w0 , then there exists
an analytic function h on C such that h.z/2 D 4z 3 g2 z g3 for every z 2 C and
h.z0 / D w0 . This
pis a very classical result. For obvious reasons, this analytic function
h is denoted by 4z 3 g2 z g3 , despite the ambiguity of the notation coming from
the fact that, for the chosen z0 there are two choices for w0 . However, in using the
square root notation, we understand that a certain pair .z0 , w0 /, as above, has been
fixed.
R1
, where the integration
Suppose now that w 2 C . Let z.w/ D w p 3 dt
4t g2 tg3
path, except for 1 and, possibly, its end point w, is contained in C and does not pass
36
through any ei , i.e. only its end point w might coincide with a root ei . Then, according
to [73, 20.221], }.z/ D w. This is the so-called integral formula for }.z/, which is
expressed by the following implication:
Z 1
dt
p
) }.z/ D w.
(3.13)
zD
3
w
4t g2 t g3
Using (3.13) we can partially solve Problem 3.3.3. Indeed, setting in (3.13) w D x.P /
we have }.z/ D x.P /. Consequently, } 0 .z/2 D 4x.P /3 g2 x.P /g3 D .2y.P //2 ,
implying "} 0 .z/ D 2y.P /, where " 2 1, 1. Then, since } 0 is an odd function, } 0 ."z/ D 2y.P /. If we choose r 2 F such that r D "z .mod /, then
.x.P /, y.P // D .}.r/, 12 } 0 .r//, as required.
This answer is partial in the sense that, given the point .x, y/ we do not know a
priori how to choose " 2 1, 1. However, in the case that A and B are real numbers
and P 2 E.R/ we will give a very satisfactory answer in Section 3.5.
A, B 2 Q,
4A3 C 27B 2 0.
(3.14)
In terms of the Weierstrass function, we have to find all z 2 C, such that .}.z/, 12 } 0 .z//
is an integral (or rational) point, where } is the Weierstrass function with parameters
g2 D 4A,
g3 D 4B.
(3.15)
37
38
3.4.1 > 0
In Figure 3.1 we translate the closure of the period parallelogram marked out by vertices
, by the vector with initial point 12 .!1 C!2 / and end point 0, and then we remove
the lower and the left side (with their vertices) of the resulting parallelogram. We thus
obtain the fundamental parallelogram, the vertices of which are marked by in Figure 3.1; by its definition, this fundamental parallelogram contains only the upper-right
vertex. We examine for which z 2 P we have }.z/ 2 R:
}.z/ 2 R , }.z/ D }.z/ , }.z/ D }.z/ ) } 0 .z/ D } 0 .z/.
Therefore, if z 2 P and }.z/ 2 R, then the point .}.z/, } 0 .z// is equal to either .}.z/, } 0 .z// or .}.z/, } 0 .z//, hence, either z D z .mod / or z D z
.mod /; equivalently,
1
1
.z z/ 2 .
(3.16)
2
2
Remember now that P includes only the right vertical and the upper horizontal side.
If z is on one of these sides of P, then (3.16) is true with the plus or the minus sign,
respectively. Next, consider points z in the interior of P. Since the only point of 12
lying in the interior of P is 0, it is easy to see that (3.16) is true with the plus or minus
sign if and only if z is, respectively, on the line segment joining !22 with !22 or the
line segment joining !21 with !21 .
Now we consider the small grey closed parallelogram P0 with vertices 0, !21 ,
!1 C!2 !2
, 2 . Our discussion above, along with the fact that the function } is even, imply
2
1 + 2
1 1
2
P
P0
1
2 1
fundamental
parallelogram
1 (1 + 2 )
2
period
parallelogram
1 ( + )
2
1
2
1 ( )
2
1
2
1 ( )
1
2
2
Figure 3.1. The fundamental parallelogram and the period parallelogram: Case > 0.
39
that, as z runs along the contour of P0 , all possible real values of }.z/ are obtained.
Therefore, we study the behaviour of the function } on every side of P0 .
The function .0, !21 3 z 7! }.z/ 2 R is 1-1. Indeed, suppose that z1 z2 are on
the domain of this function and }.z1 / D }.z2 /. Then } 0 .z2 / D } 0 .z1 /. If the plus
sign holds, then z1 D z2 .mod /, clearly impossible. If the minus sign holds, then
.}.z2 /, } 0 .z2 // D .}.z1 /, } 0 .z1 // and, consequently, z2 C z1 2 , again impossible, since 0 < z1 Cz2 < !1 . Observe now that, in view of (3.2), limz!0C }.z/ D C1
and, in view of (3.6), }. !21 / D e1 . Therefore, the continuous function above is 1-1 and
takes on all values in the interval .C1, e1 . It follows that e1 is the largest among the
roots (3.6) of the polynomial 4X 3 g2 X g3 . For, otherwise, some z 2 the interval
.0, !21 / would exist such that }.z/ is one of the numbers (3.6), other than e1 . This
2
would imply that either z !22 or z !1 C!
belongs to , which is impossible.
2
2
Next consider the function } restricted to the side with end points !21 and !1 C!
.
2
With similar arguments we conclude the following: Along this side, the function } is
1-1; at the vertex !21 it takes the value e1 which, as already seen, is the least among
2
/ (D e2 , by (3.6)) which,
the roots (3.6); at the other vertex it takes the value }. !1 C!
2
therefore, must be the second larger root; consequently, e3 D }. !22 / is the least root
of this polynomial.
2
Similarly, as z runs along the side with end points !1 C!
and !22 the values of }.z/
2
vary in 1-1 way from the middle root e2 to the smallest root e3 .
Finally, as z runs along the side with end points !22 and 0, the values of }.z/ vary
in a 1-1 way from the smallest root e3 to 1.
Now we study the behaviour of the function } 0 .
We have seen that }.z/ is strictly decreasing as z runs along the side with end points
0 and !21 , hence, for these zs we have }.z/ e1 . Then, } 0 .z/2 is a non-negative real
number, hence } 0 .z/ is a real number; more precisely, since }.z/ is strictly decreasing,
we have } 0 .z/ 0 with equality holding if z D !21 .
Next, we see that the values of }.z/ as z runs along the points of the upper horizontal
2
side of P0 following the direction from !1 C!
to !22 , coincide with the values of the
2
real function 0, !21 3 x 7! }.x C !12C!2 / 2 e3 , e2 , which is continuous, mapping
0 to e2 and !21 to e3 . Therefore, this is a strictly decreasing function, hence its derivative
takes negative values in the open interval .0, !21 / and becomes zero at 0 and !21 . But
the derivative values of the above function coincide with those of } 0 .z/ as z runs along
the upper horizontal side of P0 , hence } 0 .z/ < 0 for z in the interior of this side and
} 0 .z/ D 0 if z is a vertex.
When z runs along the vertical sides of the parallelogram P0 , the values of }.z/
belong to either the interval e2 , e1 (if z is on the right vertical side) or to the interval
.1, e3 / (if z is on the left vertical side), hence } 0 .z/2 0 and, except if z is an end
point of a side, .}.z/, } 0 .z// 62 E.R/.
Summing up: The only points z in P0 for which .}.z/, 12 } 0 .z// 2 E.R/ are those
of the two horizontal sides. For those zs, we have } 0 .z/ 0, with equality only on
40
...................
!2 C!1
2
. !1 . . . . . . . . . . . . . . . . . . . . . . . !1
2
!2 !1
2
!2
2
!2 C!1
2
!1
2
!1
2
}.z/
e2
&
e3
e2
e1
C1
&
e1
} 0 .z/
>0
<0
C1
# jump
1
infinite
3.4.2 < 0
2
In this case, .1 , 2 / D .!1 , !1 C!
/ is a fundamental pair of periods (see page 37);
2
1
C
note that 2 2 1 2 i R . In Figure 3.2, a fundamental parallelogram P is obtained as
follows: the closure of the period parallelogram marked out by vertices
is translated
!
by the vector 12 01 and then, from the resulting parallelogram, the sides left of the
vertical diagonal are removed, together with their vertices. The parallelogram P is
12 1
fundamental
parallelogram
2 12 1
1
2 1
2 + 12 1
period
parallelogram
Figure 3.2. The fundamental parallelogram and the period parallelogram: Case < 0.
41
1
2
1
2
}.z/
e1
C1
&
e1
} 0 .z/
C1
# jump
1
infinite
marked out by the vertices , but only the vertex 21 is included. We examine for
which z 2 P (Figure 3.2) we have }.z/ 2 R. As in the case of > 0, we see that
}.z/ 2 R is equivalent to 12 .z z/ 2 P \ 12 . It can be easily checked that the last
relation is true with the plus or the minus sign if and only if z is on the vertical or the
horizontal diagonal of P, respectively. By (3.6), }. 21 / D }. !21 / D e1 , hence this is
the only real root of 4X 3 g2 X g3 .
Now, assume that z runs along the horizontal diagonal of P. Then }.z/ 2 R and, by
(3.2), limz!0 }.z/ D C1. Moreover, } is 1-1 on the intervals . 21 , 0/ and .0, 21 ,
respectively. Therefore, as z moves from 21 to 0 along the real axis, }.z/ increases
from e1 to C1. In particular, this implies that } 0 .z/ > 0. On the other hand, as z
runs along the interval .0, 21 /, }.z/ decreases from C1 to e1 which, in particular,
implies that } 0 .z/ < 0.
Next, assume that z runs along the vertical diagonal of P. If z tends to 2 12 1 ,
then }.z/ tends to }. 12 1 / D e1 . For the interior points z of the upper half of this
diagonal we have, by (3.2), }.z/ < 0 and limz!0i }.z/ D 1, hence the values
of }.z/ belong to the interval .1, e1 /. It follows that 4}.z/3 g2 }.z/ g3 < 0.
Consequently, } 0 .z/ 2 i R and no real point on E arises when z is on the open upper
half of the vertical diagonal of P. Similarly the same is true if z is in the open lower
half of this diagonal. In Table 3.2 we schematically summarise our results.
3.4.3
We refer to our discussion, notation etc. following Problem 3.3.3 up to the end of
Section 3.3.
way, so that the real half-line from
Computation of !1 . We cut C in the appropriate
p
e1 to C1 is contained in C and consider 4x 3 g2 x g3 as the
panalytic continuation on C of the positive real-valued function .e1 , C1/ 3 x 7! 4x 3 g2 x g3 .
With w D e1 in (3.13) and integration path the half-line from e1 to C1 we have
42
}.z/ D e1 D }. 12 !1 /, hence
Z 1
dt
!1
p
3
2
e1
4t g2 t g3
.mod /
(observe that !21 !21 .mod /). Denote, temporarily, by I the integral in the
right-hand side of the above relation. Since I 2 RC , this relation can be written as
I D 12 !1 C n!1 for some non-negative integer n. We claim that n D 0. Suppose that
n > 0. Consider the real function
Z C1
dt
p
t 7! I.x/ D
3
x
4t g2 t g3
defined on the interval .e1 , C1/, which is continuous and strictly decreasing. Since
limx!e1 C I.x/ D I , this function takes values in the interval .0, 12 !1 C n!1 /, therefore, there exists x0 > e1 such that I.x0 / D 12 !1 . This says that the relation (3.13)
holds with w D x0 and z D 12 !1 , which implies that x0 D }.!1 =2/ D e1 , a contradiction. Therefore, n D 0 and
Z C1
!1
dt
p
,
(3.17)
D
2
e1
4t 3 g2 t g3
hence, in terms of the coefficients of our Diophantine equation (3.14), we obtain, in
view of (3.15),
Z C1
dt
p
.
(3.18)
!1 D
e1
t 3 C At C B
Computation of !2 . We distinguish the two cases > 0 and < 0. We put
q.t / D 4t 3 g2 t g3 .
Let > 0. In this case we choose C to be the complex plane from which we
have removed the open line segment with end points e2 , e3 and the open half-line on
axis are the intervals
the real axis from e1 to C1. What remains in C from the real p
.1, e3 and e2 , e1 on which the real non-negative function
q.x/ is defined.
p
We consider therefore the analytic function C 3 z 7! i q.z/ 2 C which, for
z 2 .1, e3 [ e2 , e1 , takes values in i RC . Clearly, this is an
p analytic branch in
of
the
square
root
of
q.x/.
Therefore,
if
in
(3.13)
we
take
i
q.t / in place of
C
p
3
4t g2 t g3 , e3 in place of w and, as an integration path, the half-line on the real
axis from e3 to 1, then we have
Z 1
dt
p
, }.z/ D e3 D }.!2 =2/.
zD
e3
i q.t /
It follows that
Z 1
!2
dt
i
.mod /.
p
2
e3
q.t /
43
fined on the interval .1, e3 /. It is an increasing function with values in the interval .0, . 12 C n/s/. Therefore, there exists x0 < e3 such that I.x0 / D s=2, hence
R 1
pdt
D !2 =2 and, consequently, x0 D
iI.x0 / D !2 =2. This means that x0
i
q.t/
e3
1
Case > 0
Z e3
dt
dt
p
p
Di
.
q.t /
.t 3 C At C B/
1
(3.19)
Let < 0. Now we cut the complex plane along the half-line on the real axis
from e1 to C1 and along the segment joining e2 with e3 D e2 . In analogy to the
previous
p case we observe that q.x/ > 0 for x < e1 , the real function .1, e1 3
x 7!
q.x/
2 0, C1/ is well defined and is extended analytically on C . As
R 1
pdt , then }.z/ D e1 D }.!1 =2/. Now, a fundamental pair of
before, if z D e1
i q.t/
R e1
R 1
p dt
pdt
periods is .1 , 2 / D .!1 , .!1 C!1 /=2/, therefore i 1
D e1
D
q.t/
q.t/
2
C m!1 C n !1 C!
D . 12 C m C n2 /!1 C n2 !1 , for some m, n 2 Z. The left-hand
2
C
side, as well
R e1 as !dt2 belongn to i R , therefore, the coefficient of !1 must be zero, n is
D s, where !2 D i s. In the last but one equality, the integral
odd and 1 p
!1
2
q.t/
2
x0 D }. !22 / D }. !1 C!
2
and we conclude that
q.t/
Case < 0
Z e1
dt
Di
.
!2 D 2i
p
p
3
q.t /
.t C At C B/
1
1
Summing up we have the following theorem.
Z
e1
dt
(3.20)
44
anC1
a0 D a,
an C bn
,
D
2
bnC1
b0 D b,
p
D an bn
.n 0/.
2
ds
a2 cos2 s C b 2 sin2 s
2
.
M.a, b/
(3.21)
e2 t e1 e2 C e1 e3 e2 e3
t e2
C1
e1
dt
p
D
q.t /
e3
e2
dt
p
q.t /
and then, the change of variable t D e3 C .e2 e3 / sin2 s transforms the last integral
into
Z
1 2
ds
p
,
2
4 0
.e1 e3 / cos s C .e1 e2 / sin2 s
p
p
which, by (3.21) is equal to 2=M. e1 e3 , e1 e2 /.
45
R e1
e2
p dt
q.t/
which, in
!1 D
(3.22)
(3.23)
(3.24)
46
p
2
(3.25)
p in place of e3 < e2 < e1 , so that, in place
p of e1 e3 we have 2p 12e1 g2 D
2
2
2
4 3e1 C A, in place of e1 e2 we have
3e1 C 2 3e1 C A, and
p 3e1 C 12e1 g2 D p
in place of e2 e3 we have 3e1 C 12e12 g2 D 3e1 C 2 3e12 C A. Therefore,
q
w1 D
2M
w2 D
3e12 C A ,
q
2M
1
2
q
p
3e1 C 2 3e12 C A
i
r
3e12 C A ,
1
2
q
3e1 C 2 3e12 C A
(3.26)
!.
.2w2 , w1 /,
.2w1 , w1 C w2 /.
The first two cases are easily excluded in view of the conditions 1 2 RC and 1
1
2 i RC , so that the third case remains. Instead of .1 , 2 / D .2w1 , w1 C w2 /
2 1
we prefer to take as a fundamental pair of periods .!10 , !20 / D .2 , 1 2 / D
.w1 C w2 , w1 w2 /.
Now we relate !1 , !2 with w1 , w2 . Since the pair .!10 , !20 / generates all periods, we
have !1 D m!10 C n!20 D .m C n/w1 C .m n/w2 for some integers m, n. Since
!1 is a real period, we must have m D n and !1 D m.2w1 /. But 2w1 D !10 C !20 is
a positive real period, while !1 is the least real period. Therefore, m D 1, !1 D 2w1
and, consequently, from (3.26) we obtain
!.
!1 D
(3.27)
r
q
q
1
4
M
3e12 C A , 2 3e1 C 2 3e12 C A
Section 3.5
47
i
.
q
p
p
M 4 3e12 C A , 12 3e1 C 2 3e12 C A
(3.28)
Since .!10 , !20 / D .w1 C w2 , w1 w2 / D ..!1 C !2 /=2, .!1 !2 /=2, we can take the
last pair as a fundamental pair of periods. Note that .!10 C !20 , !10 / D .!1 , .!1 C !2 /=2
is also a fundamental pair of periods.
Case < 0. As a fundamental pair of periods we can take any of
.
!1 C !2 !1 !2
,
/,
2
2
.!1 ,
!1 C !2
/,
2
3.5
We keep the notation of Theorem 3.4.1 and go back to the end of Section 3.3. We
will revisit Theorem 3.3.1 now focusing our interest exclusively on points of E.R/.
Our first task will be, given a point P 2 E.R/, to compute explicitly .P /; see
Theorem 3.5.1 below. Next, in Theorem 3.5.2 we will focus on the set E0 .R/ of real
points .x, y/, such that x e1 , showing that, actually, the restriction of to E0 .R/
establishes a group isomorphism between E0 .R/ and R=Z!1 . Finally, we will show
that, in the case of a positive discriminant, this last isomorphism can be extended to a
two-to-one group epimorphism E.R/ ! R=Z!1 .
Note first that E0 .R/ is indeed a subgroup of E.R/. For, if < 0, then E0 .R/ D
E.R/. If > 0, then E.R/ D E0 .R/ [ E1 .R/, where E1 .R/ is the bounded set of
real points .x, y/, such that e3 x e2 , i.e. the set of points on the egg, in pop
terminology. It is geometrically obvious that any line of the plane having a common
point with the egg must have two common points with it, where in two we also
include double points, in which case the line is tangent to the egg. Consequently, if
a line joins two points belonging to E0 .R/ (with the meaning of two as above),
then the third point of intersection must also lie on E0 .R/ and not on the egg; for
otherwise, the line would have four common points multiplicities taken into account
with the cubic model, which is absurd. From this we conclude that E0 .R/ is closed
under addition of points and, consequently, E0 .R/ is indeed a subgroup of E.R/,
hence a subgroup of E.C/ as well. Moreover, the sum of two points on the egg is a
point on E0 .R/ and the sum of a point on the egg and a point on E0 .R/ is a point
on the egg; in other words the group E.R/=E0 .R/ is isomorphic to Z2 .
First, let P D .x.P /, y.P // 2 E0 .R/. From Tables 3.1 and 3.2 and the conclusions of the relevant subsections we see that there exists a unique real number r in
48
the interval . 12 !1 , 12 !1 such that .x.P /, y.P // D .}.r/, 12 } 0 .r//. We find now an
explicit formula for r. Consider the real integral
Z
Z C1
dt
dt
1 C1
.
p
D
p
zD
3
2 x.P / t 3 C At C B
x.P /
4t g2 t g3
According to the discussion following the integral formula (3.13) for }.z/, with a
convenient " 2 1, 1, we have .x.P /, y.P // D .}."z/, 12 } 0 ."z//. On the other
hand, since P 2 E0 .R/, we have x.P / e1 , therefore 0 < z 12 !1 , in view of
(3.17). From Tables 3.1 and 3.2 we see that } 0 .z/ 0, concluding thus that " D 1
or 1 according to whether y.P / 0 or y.P / 0, respectively. Therefore, in this
case,
Z C1
1
dt
r D "P
,
(3.29)
p
3
2
x.P /
t C At C B
where "P D 1 or 1, according to whether y.P / > 0 or y.P / 0, respectively.
Next, let P D .x.P /, y.P // 2 E1 .R/. Note that this can occur only in the case
of positive discriminant. In (3.7) we defined the points Qi D .ei , 0/ (i D 1, 2, 3).
2
Theorem 3.3.1 and (3.6) imply that .Q1 / D !21 , .Q2 / D !1 C!
and .Q3 / D !22 ,
2
therefore we may assume that P Q2 , Q3 . We refer to Section 3.4.1. From Table
3.1 and the conclusions of that section we see that there exists an r on the line segment
1
1
with end points !2 !
(excluded) and !2 C!
(included), such that .x.P /, y.P // D
2
2
1 0
}
.r//.
We
will
find
this
r.
.}.r/, 2
Observe first that P C Q2 2 E0 .R/ and y.P C Q2 / have the same sign as y.P /.
We have, in view also of our conclusion about points of E0 .R/,
.P / D
..P C Q2 / C Q2 / D .P C Q2 / C .Q2 /
Z C1
!1 C !2
dt
1
C
C ,
p
D "P
2
2
x.P CQ2 /
t 3 C At C B
..P C Q2 / C Q2 / D .P C Q2 / C .Q2 /
Z C1
!1 C !2
dt
1
C m!1 C n!2
C
p
D "P
3
2
2
x.P CQ2 /
t C At C B
1
1
belongs to the line segment with end points !2 !
and !2 C!
as mentioned above.
2
2
!1
!1
!2
This condition is equivalent to 2 <.r/ 2 and i =.r/ D 2 . The second condition
immediately implies n D 0 and the first condition implies that m D 0 if "P D 1 and
m D 1 if "P D 1. Therefore,
Z C1
dt
!2 C "P !1
1
.
(3.30)
C
p
r D "P
2
2
x.P CQ2 /
t 3 C At C B
Section 3.5
49
Note that, if we adopt the convention x.O/ D C1, then formula (3.30) with P D
1
1
D !2 C!
.mod /. Also, with P D Q3 we similarly obtain
Q2 , gives r D !2 !
2
2
x.Q3 C Q2 / D x.Q1 / D e1 and then formula (3.30), in combination with (3.18),
gives r D !22 . Thus, formula (3.30) is valid for all points on P 2 E1 .R/.
Our previous discussion leads to the following theorem:
Theorem 3.5.1. Let E be the short Weierstrass model
E : y 2 D x 3 C Ax C B,
A, B 2 R,
D 16.4A3 C 27B 2 / 0,
50
l0 .P /
l.P / D
l0 .P C Q2 /
if P 2 E0 .R/
if P 2 E1 .R/.
(3.31)
Section 3.5
51
Finally, is also onto, for, if r is a real number and we set P D .}.r/, 12 } 0 .r// (if
r 2 we mean, of course, that P D O) then, by Theorem 3.5.1, .P / D r C and,
consequently, by the definition of l0 , we have l0 .P / D r C Z!1 .
Proof of (b). The map l is a group epimorphism because, on the one hand, its restriction to E0 .R/ is, in view of (a), a group isomorphism and, on the other hand,
the group E.R/=E0 .R/ is isomorphic to Z2 . For example, if P1 , P2 2 E1 .R/ then
P1 C P2 2 E0 .R/ and, therefore, l.P1 C P2 / D l0 .P1 C P2 /. On the other hand,
l.P1 /Cl.P2 / D l0 .P1 CQ2 /Cl0 .P2 CQ2 / D l0 .P1 CP2 C2Q2 / D l0 .P1 CP2 CO/ D
l.P1 C P2 /, as required.
Finally, let P1 , P2 2 E.R/ with l.P1 / D l.P2 /. If P1 , P2 2 E0 .R/, then l0 .P1 / D
l.P1 / D l.P2 / D l0 .P2 /, therefore, by the second preliminary observation, P1 D P2 .
Next, suppose that > 0. If P1 , P2 2 E1 .R/, then, Pi C Q2 2 E0 .R/ for i D 1, 2, so
that l0 .P1 CQ2 / D l.P1 / D l.P2 / D l0 .P2 CQ2 / and, as before, P1 CQ2 D P2 CQ2 ,
from which P1 D P2 . If exactly one among P1 , P2 , say the second, belongs to E1 .R/,
then l0 .P1 / D l.P1 / D l.P2 / D l0 .P2 C Q2 /, hence P1 D P2 C Q2 , which is
equivalent to P2 D P1 C Q2 .
Proof of (c). Let n1 l.P1 / C C nk [Link] / C l.T / C n0 !1 D 0. Since l is a homomorphism, we have
l.n1 P1 C C nk Pk C T / D n1 l.P1 / C C nk [Link] / C l.T / C `!1 ,
for an appropriate integer ` which makes the left-hand side falling in the interval I D
. 12 !1 , 12 !1 . In view of the hypothesis, the right-hand side of the last equality is
equal to .` n0 /!1 ; since the left-hand side belongs to I , this forces ` D n0 and
l.n1 P1 C C nk Pk C T / D 0. Then, in view of (b), n1 P1 C C nk Pk C T is
equal either to O or to Q2 . The second alternative is impossible if Q2 62 E.K/; and if
Q2 2 E.K/, both alternatives imply that n1 P1 C Cnk Pk is a torsion point of E.K/,
which contradicts the hypothesis that P1 , : : : , Pk are Z-linearly independent.
Z
"P C1 dt
x.P /
if P 2 E0 .R/
, xP D
(3.32)
l.P / D
p
2 xP
x.P C Q2 / if P 2 E1 .R/
f .t /
where, as before, "P D 1 if y.P / > 0 and "P D 1 if y.P / 0 (note that y.P /
and y.P C Q2 / have equal signs).
52
0
if P D O
<
r.P /
if P 0 and y.P / 0 .
s.P / D
!1
1
/
2
1
2
Section 3.5
0 < ai1 C
di C1
22
C <
1
,
2
53
y.2i1
1
2
while i N do
if Q D O then goto ()
else
goto ()
if y.Q/ D 0 then s
s C 1i :
2
else
if y.Q/ > 0 then s
s C 1i end if
2
end if
end if
if i < N then Q D 2 Q endif
i
i C1
end while
() u D !1 s.
if P O and y.P / > 0 then u
u !1 end if
l.P / D u.
END
E0 .R/.
Chapter 4
In this chapter we make a general description of the elliptic logarithm method, briefly
Ellog, which will be used in later chapters for solving various types of elliptic equations in integers. In this chapter we will refer to the resolution of any such equation as
the Diophantine problem.
To each Diophantine problem, we will attach a convenient short Weierstrass model
E : y 2 D x 3 C Ax C B,
A, B 2 Q;
(4.1)
this is done separately for each problem in Chapters 5, 6, 7, 8 and 11, respectively.
We denote by r the rank of the elliptic curve E.Q/. We also denote by r0 the least
common multiple of the orders of the non-zero points of Etors .Q/; by Mazurs theorem
(see [29, 30], or [45, Theorem 7.5]), 1 r0 12, r0 11.
We refer to Theorem 1.2.1, which we apply for K D Q. According to that theorem,
if r > 0, then there exist points of infinite order P1E , : : : , PrE in E.Q/, with the
property that every point P E D .x.P /, y.P // 2 E.Q/ is written uniquely in the
form
(4.2)
P E D m1 P1E C C mr PrE C T E
for some integers m1 , : : : , mr and T E D .x.T /, y.T // 2 Etors .Q/.
11
Thus, (4.2) establishes a map P E ! .m, T E /, where m D .m1 , : : : , mr / is a
lattice point in Rr and T E 2 Etors .Q/.
A first important feature of Ellog is that the solutions of the Diophantine problems
in Chapters 5, 6, 7, 8 and 11 are associated in some way, particular to each Diophantine
problem, to points P E $ .m, T E / with m belonging to a (finite) cube of Rr ; that is,
if P E is the point that covers (comes from) a solution of the Diophantine problem
def
and P E is expressed as in (4.2), then M D max1ir jmi j is bounded. Our task will
be to compute an explicit upper bound for M ; a huge one first, and then with the use
of it, a bound of manageable size.
A second important feature of Ellog is that, given any point P E 2 E.Q/, one can
easily (in principle) check whether it covers an actual solution of the Diophantine
problem. Therefore, if an explicit upper bound, say K, of M is known, then, for each
choice K mi K, i D 1, : : : , r and T E 2 Etors .Q/, a point P E 2 E.Q/ is
calculated by means of (4.2) and then, this P E is checked whether it comes from a
solution of the Diophantine problem. In particular, one can explicitly find the set T
consisting of all points T of E, such that T E 2 Etors .Q/ and T covers a solution of
55
the Diophantine problem. As the calculation of the set T is a rather trivial task, we
will henceforth assume that, once we are given a certain Diophantine problem, T has
been computed already from the beginning and we are left with the difficult task of
explicitly computing a finite set P of points of E with the following properties: Each
solution of the Diophantine problem not covered by points of T , is covered by some
point of P, and P E 62 Etors .Q/ for every P 2 P. The last condition is, of course,
equivalent to
r 1 & .m1 , : : : , mr / .0, : : : , 0/.
(4.3)
In this way we will effectively solve the Diophantine problem under consideration.
Here, effectively means after finitely many steps, no matter how many they are.
Take a very modest typical example, in which r D 2, the torsion points are two and,
by some means we calculate K D 1030 , which is rather very small. Then the number of checks will be, approximately, 23 1060 . Of course, there is no hope to solve
explicitly our Diophantine problem this way, but, at least, we have an effective or, if
you prefer, a constructive proof of the fact that there exist finitely many solutions.
Certainly, this is far more satisfactory for our mind than an existential (in other words,
non-constructive) confirmation which would merely assure the existence of finitely
many integer solutions. At the same time, an effective proof of the finiteness of solutions suggests, at least in theory, a method of resolution of the problem. We will
show that one can exploit this effective method to finally obtain an explicit method of
resolution of the Diophantine problem, i.e. a method whose output will be a finite set
S consisting precisely of all solutions of the Diophantine problem.
Below we make a very general description of how an upper bound for M is obtained. We recommend the reader to revisit very quickly Chapter 1 in order to recall
the terminology, conventions and facts that will be used. In accordance with them:
We will view E as a model of an elliptic curve. Referring to (4.1), we put g2 D
4A, g3 D 4B and we adopt the notations and assumptions of Theorem 3.4.1. In
particular, q.X / D 4X 3 g2 X g3 and e1 , e2 , e3 denote the roots of q.X /, where
e1 > e2 > e3 if all three roots are real; otherwise e1 is the only real root. Further, we
put f .X / D q.X /=4 D X 3 C AX C B.
In order that the group relation (4.2) provide us with a numerical relation, we apply
the group epimorphism l : E.R/ ! R=Z!1 of Theorem 3.5.2. We remind here that,
according to this theorem, we view the group R=Z!1 as the interval . !21 , !21 in
which the group operation is the usual addition adjusted by an integral multiple of
!1 , so that the final sum falls in the above interval.
The group relation (4.2) implies the following numerical relation1
0 l.P / D m1 l.P1 / C C mr [Link] / C l.T / C m0 !1 ,
1
(4.4)
Below we omit the indication E from the points when they are an argument of the elliptic logarithm
function l; see Conclusions and remarks 3, page 52.
56
where m0 !1 is the adjusting summand mentioned a few lines above. This is a linear
form in elliptic logarithms of points with algebraic coordinates of degree at most three;
see (2) in Conclusions and remarks, page 52. The non-vanishing of l.P / is justified
by (4.3) and Theorem 3.5.2 (c).
If t 2 is the (group) order of a torsion point T E O, then t l.T / D l.t T / D
l.O/ D 0 2 R=Z!1 , which shows that l.T / D st !1 with 12 < st 12 . Thus, in
general,
s
1
s
1
r0
l.T / D , t jr0 , t > 0, jsj , < .
(4.5)
t
2
2
t
2
Therefore, we can rewrite (4.4) as
s
(4.6)
0 l.P / D m0 C !1 C m1 l.P1 / C C mr [Link] /.
t
Since every l-value is absolutely less than !21 , it follows from (4.6) and (4.5) that
jm0 j 1 C .jm1 j C C jmr j/=2. Therefore, if we set
M D max jmi j,
1ir
then
M0 D max jmi j
0ir
M0 max M , 12 rM C 1 .
(4.7)
(4.8)
Closely related to l.P / is a linear form L.P /, which we will define for each particular case in Chapters 5, 6, 7, 8 and 11. In each case we will compute an explicit upper
bound for jL.P /j in terms only of M ; denote it (temporarily) by Ub.M /. On the other
hand, the application of Theorem 9.1.2 to L.P / will furnish either an explicit upper
bound for M itself, and we are done, or an explicit lower bound for jL.P /j in terms
only of M , which we will temporarily denote by Lb.M /. In the second case we will
thus arrive to
Lb.M / jL.P /j Ub.M /.
Both functions Lb.x/ and Ub.x/ are decreasing for x positive and sufficiently large.
What is very important is that there exists an explicit positive constant K such that
Lb.x/ > Ub.x/ for every x > K. This, combined with the last displayed relation,
clearly implies that M K, and we obtain the required upper bound for M .
Remark on terminology. In the chapters that follow, the term constant always
means number depending on the particular Diophantine problem. Once the Diophantine problem is expressed in concrete numerical terms, the constant has a specific numerical value.
In the constants of the form ci , which are introduced in this book, we use subscripts
i 7 in order to avoid confusion with the standard symbols c4 and c6 , that are used
in the general theory of elliptic curves (see [45, p. 42]).
Chapter 5
v D
3 y C
2 x C ,
(5.2)
58
as claimed.
The main result of this section is the following.
Theorem 5.2. Let E be the elliptic curve represented by the model C in (5.1). Let
E : y 2 D f .x/ be the short Weierstrass model (4.1) which is obtained from C by
means of the transformation (5.2), so that everything in Chapter 4, especially relations
(4.4) through (4.8), refers to this particular model E.
Assume that P 2 E is such that P C D .u.P /, v.P // has integer coordinates with
u.P / > 2
2 max 12 , je1 j, je2 j, je3 j C ,
where, as usual, e1 , e2 , e3 are the roots of f .X /.
Write P E as in (4.2). With l.P / as in (4.6), set L.P / D l.P /. Then, with M , ,
and defined, respectively, by the first relation (4.7) and by Propositions 2.6.2, 2.6.3
and 5.1, we have
p
(5.3)
jL.P /j 2 2 exp. C 12 log M 2 /.
Proof. First we observe that all x 2 maxje1 j, je2 j, je3 j satisfy the inequality
Z 1
p
dt
p
(5.4)
0<
4 2jxj1=2 .
x
f .t /
Indeed, for such an x and t x we have 0 < f .t / D jt e1 jjt e2 jjt e3 j. As t is
larger than the absolutely largest root of the polynomial f .X /,p
it follows that jt e1 j
t je1 j t =2, and likewise for e2 and e3 . Consequently, 1= f .t / 23=2 t 3=2 and
hence, for all N > x,
Z N
Z N
p
dt
0<
p
23=2 t 3=2 dt D 4 2.x 1=2 N 1=2 /.
x
x
f .t /
Letting N tend to infinity we obtain (5.4).
59
Now, in view of (5.2) and the hypothesis about the lower bound of u.P / we have
x.P / > 2 maxje1 j, je2 j, je3 j, therefore, by (5.4),
Z 1
p
dt
0<
p
4 2 x.P /1=2 .
x.P /
f .t /
But, x.P / e1 , therefore P E 2 E0 .R/ and consequently, by the definition of l (see
(3.32)), thepleft-hand side in the above displayed relation is equal to 2jl.P /j, so that
jl.P /j 2 2 x.P /1=2 . By Propositions 5.1 and 2.6.2,
logx.P /
1
2
O / C 1 log M 2 ,
C 12 log h.P
2
Chapter 6
In this chapter we will deal with the integer solutions of Diophantine equations
F .U / D V 2 , where F denotes a quartic polynomial with rational coefficients and
the equation F .x/ D y 2 defines an elliptic curve over Q. This implies, first, that the
genus of the corresponding curve is one, which is equivalent to the non-singularity
of the curve, hence to the non-vanishing of the discriminant of the polynomial F .X /,
and second, that we know a rational solution .U , V / D .u0 , v0 /. In practice, the search
for such a solution can be done by some powerful computational tool like ratpoints
[53] due to M. Stoll [52], or the routine Points of MAGMA; of course, although it is
highly probable that the search will be successful, there is no guarantee for this. Anyway, assuming this, we then reduce our equation to solving in integers the equation
Q.u/ D v 2 , where Q.u/ D F .u C u0 /, .u, v/ D .U u0 , V / and the constant term
of the polynomial Q is a square of an integer. Therefore, we will consider the elliptic
curve E, a model of which is on the quartic elliptic model
C : g.u, v/ D 0,
g.u, v/ D v 2 Q.u/,
(6.1)
a, b, c, d , e 2 Q , a, e > 0 , discr.Q/ 0.
(6.2)
Before proceeding we note that we may assume that a is not a perfect square, for, in
this case, the equation Q.u/ D v 2 in integers can be solved very efficiently by totally
elementary means; see [37].
According to our discussion in Chapter 4, to P C D .u.P /, v.P //, the generic point
on C with integer coordinates, we will attach a linear form in elliptic logarithms L.P /;
see page 56 immediately after (4.8). We intend to prove Theorem 6.8, which gives an
explicit upper bound for jL.P /j, as a first main step towards the explicit determination
of all points P C as above. Our presentation is based on the authors paper [62].
Throughout the chapterpwe keep the notations, assumptions and results of Chapter
4. Wherever the symbol
is used, always, care is taken that the radicand be nonp
negative and then denotes the non-negative square root of the radicand.
The short Weierstrass model (4.1) with
1
A D c 2 C bd 4ae 2 ,
3
BD
8
2 3 1
c bcd ace 2 C b 2 e 2 C ad 2
27
3
3
(6.3)
61
is birationally equivalent to (6.2) and the birational transformations (1.3) and (1.4)
between the two models (cf. Fact 1.1.2) are the following:
6ev C cu2 C 3du C 6e 2
,
3u2
ebu3 C 2ecu2 C 3edu C duv C 4e 2 v C 4e 3
y D Y.u, v/ D
u3
(6.4)
(6.5)
x D X .u, v/ D
1
v D V.x, y/ D .32c 3 e 3 108b 2 e 5 C 108bcde 3 C 27 d 3 y C 216 be 4 y
3
108 de 2 cy C 216e 3 x 2 c 81ex 2 d 2 C 9ec 2 d 2 C 216e 3 x 3
27bed 3 108e 3 y 2 /=.6ye 3dx C cd 6be 2 /2 .
(6.6)
p
sgn.d a C eb/
D
p
sgn.8e 3 a C 4e 2 c d 2 /
p
a C eb 0
.
p
if d a C eb D 0
if d
p
R.u/ D beu3 C 2ceu2 C 3deu C 4e 3 C .4e 2 C du/ Q.u/.
(6.7)
Then, for sufficiently large u > 0 we have Q.u/ > 0 and sgn.R.u// D .
p
p
p
Proof. If p
d a C eb D 0 and 8e 3 a C 4e 2 c d 2pD 0, then b D d a=e and
c D 2e a C d 2 =.4e 2 /. Consequently, Q.u/ D . au2 .d=2e/up e/2 , which
contradicts
p the hypothesis discr.Q/ 0. Therefore, at least one of d a C eb D 0
3 a C 4e 2 c d 2 D 0 is non-zero.
and 8e
p
p
3
Suppose first that d a C be 0. Then
p limu!1 R.u/=u D d a C be, which
proves our claim in this case. Next, let d a C be D 0. Then
r
r
b
b
3de
4e 3
R.u/
2
D
.be
C
d
a
C
a
C
.
C
/u
C
.2ce
C
4e
C
/
C
C
u2
u
u
u
u2
p
p
As u tends to infinitypthe right-hand side tends to bd=.2
a/C.2ce C4e 2 a/C0C0,
p
which, in view of d a C be D 0, is equal to .8e 3 a C 4e 2 c d 2 /=.2e/, hence of
the same sign as .
62
Now, let
p
c
(6.8)
x0 D 2e a C .
3
Following the notation in Chapters 3 and 4, we set f .X / D X 3 C AX C B and
denote by e1 , e2 , e3 the roots of f .X /, enumerating them as explained on page 55. The
following proposition describes the position of x0 relatively to e1 , e2 , e3 .
p
p
Lemma 6.2. If d a C eb 0, then either x0 > e1 or e3 < x0 < e2 . If d a C
eb D 0, then x0 D e1 and D C1.
p
Proof. p
It is easy to check that f .x0 / D x03 C Ax0 C B D .d a C be/2 .
If d a C be 0, then f .x0 / > 0 and, consequently, either x0 > e1 , or (in the
case that
p e2 , e3 2 R) e3 < x0 < e2 .
If d a C eb D 0, then x0 is a root of f .X / and we set f .X / D .X x0 /h.X /,
X CAp
C x02 . In particular, e2 , e3 are the roots of h.X /. But
where h.X / D X 2 C x0p
2
3
h.x0 / D p3x0 C A D a .8e a p
C 4ce 2 d 2 /=e; for the last equality we used
b D d a=e. By Lemma 6.1, 8e 3 a C 4ce 2 d 2 0, hence x0 e2 , e3 and,
if e2 , e3 2 R,
consequently, x0 D e1 . Moreover, if e2 , e3 62 R, then h.x0 / > 0; and p
then x0 > e2 , so that again h.x0 / > 0. Hence, D sgn.8e2 a C 4ce a C bd / D
sgn.h.x0 // D C1.
Lemma 6.3. Let u0 1 be such that for u u0 the conclusion of Lemma 6.1 holds.
For u u0 define
p
cu2 C 3du C 6e 2 C 6e Q.u/
.
(6.9)
x.u/ D
3u2
Then, in the interval .u0 , C1/, x is strictly increasing or strictly decreasing, according to whether D 1 or C1, respectively. Moreover, limu!C1 x.u/ D x0 .
R.u/
Proof. For the first claim, just observe that dx.u/
D 3p
and combine this with
du
u
Q.u/
Lemma 6.1. The claim about the limit is obvious.
63
p
d a C be 0 and x0 > e1 . Then, x is strictly increasing when D 1 and
strictly decreasing when D C1, with limiting value (in both cases) x0 , as u !
C1. Hence, for an appropriate u u0 , it is indeed true that u u ) x.u/ 2
.x0 , C1/ .e1 , C1/.
p
d a C be D 0 and x0 D e1 . Then, necessarily, D C1, hence x is strictly
decreasing with limiting value x0 , as u ! C1. Same conclusion as in the previous
bullet.
p
d a C be 0 and x0 2 .e3 , e2 /. Then x is strictly decreasing when D C1 and
strictly decreasing when D 1 with limiting value (in both cases) x0 2 .e3 , e2 /, as
u ! C1. Therefore, for an appropriate u u0 , it is true that u u ) x.u/ 2
.x0 , e2 / .e3 , e2 / if D C1 and x.u/ 2 .e3 , x0 / .e3 , e2 / if D 1; hence, in
both cases, x.u/ 2 .e3 , e2 /.
Let u be as in the above bullets. Then, for u u0 , the interval .u, C1/ is
mapped by x onto the open interval with end points x.u/ and x0 , which is a subset
of .e1 , C1/ in the case of the first two bullets and a subset of .e3 , e2 / in the case
of the last one. Therefore, in the right-hand side of (6.10)p
we can make the change of
variable t D x.w/ D .6ez C 3t d C ct 2 6e 2 /=.3t 2 /, z D Q.w/,
so that
dt D
dt
f .t /
D
dw
dw
,
D p
z
Q.w/
64
and y D R.u/=u3 , where R.u/ is given in (6.7). In view of Lemma 6.1, R.u/ > 0,
as claimed.
From now on, P will denote a point on E such that v.P / > 0 and u.P / is an integer
u , where u is as in Proposition 6.4.
Also, we define P0 2 E by
p
(6.11)
P0E D .x.P0 /, y.P0 // D .x0 , .be C d a//
(cf. beginning of proof of Lemma 6.2).
We will use, implicitly, Proposition 6.5 as well as the conclusions of the bullets
in the proof of Proposition 6.4. We will also use the following lemma.
Lemma 6.6. If e2 , e3 2 R and e3 x e2 , then
Z x
Z C1
dt
dt
p
p
D
,
e2
x0
f .t /
f .t /
where
.e1 e2 /.e2 e3 /
2 e1 , C1/
x 0 D e2 C
e2 x
In particular, in the case that x D x.P / 2 e3 , e2 , we have x 0 D x.P C Q2 /, where
Q2 is the point of E with Q2E D .e2 , 0/ (see (3.7)).
Proof. For t 2 x, e2 / we make the change of variable
.e1 e2 /.e2 e3 /
.e1 e2 /.e2 e3 /
t D e2 C
.
so that t 0 D e2 C
0
e2 t
e2 t
The function t 7! t 0 is strictly increasing in the interval x, e2 /, with values in
x 0 , C1/. A simple calculation now proves the equality of the integrals, as stated
in the announcement of the lemma.
Concerning the identity x 0 D x.P C Q2 /, its proof is a matter of an easy (even by
hand) symbolic calculation.
p
Now we proceed further, supposing first that eb C d a 0. According to
Lemma 6.2, either x0 > e1 , or e2 , e3 2 R and x0 2 .e3 , e2 /.
If x0 > e1 , then (first bullet in the proof of Proposition 6.4), x.P / D x.u.P // >
e1 , so that we can write
Z x.P /
Z C1
Z C1
Z C1
Z C1
dt
dt
dt
dt
dt
D
D
.
p
p
p
p
p
f .t /
f .t /
f .t /
f .t /
f .t /
x0
x0
x.P /
x.P0 /
x.P /
(6.12)
If e2 , e3 2 R and e3 < x0 < e2 , then, for t 2 e3 , e2 / we write
t 0 D e2 C
.e1 e2 /.e2 e3 /
e2 t
65
(cf. proof of Lemma 6.6). By the third bullet in the proof of Proposition 6.4, x.P / D
x.u.P // 2 .e3 , e2 / and, using Lemma 6.6, we compute
Z x.P /
Z e2
Z e2
dt
dt
dt
p
p
p
D
x0
x0
x.P /
f .t /
f .t /
f .t /
Z C1
Z C1
dt
dt
D
p
p
0
x0
x.P /0
f .t /
f .t /
Z C1
Z C1
dt
dt
D
.
(6.13)
p
p
x.P0 CQ2 /
x.P CQ2 /
f .t /
f .t /
p
Next, suppose that be C d a D 0. Then, by Lemma 6.2, x0 D e1 , D C1 and
the second bullet in the proof of Proposition 6.4 x.P / D x.u.P // > e1 . Therefore
we can write
Z x.P /
Z C1
Z C1
dt
dt
dt
p
p
p
D
.
(6.14)
x0
e1
x.P /
f .t /
f .t /
f .t /
All three integrals in the left-hand sides of relations (6.12)(6.14) are, by Proposition 6.4, equal to
Z C1
dw
p
u.P /
Q.w/
and, moreover, in (6.14) we know that D C1.
On the other hand, we look at the right-hand sides of relations (6.12)(6.14). For
both relations (6.12) and (6.13), the right-hand side is equal to 12 .l.P0 /"P0 Cl.P /"P /,
in view of the relation (3.32). Further, in this case, by the definition of P0 , we have
y.P0 / > 0, hence "P0 D 1. Also, since u.P / > u0 and v.P / > 0, Proposition 6.5
asserts that y.P / > 0, hence "P D 1 and, consequently, the right-hand side of both
(6.12) and (6.13) is equal to 12 l.P0 / C 12 l.P /. In the case of the relation (6.14), the
first integral is equal to !1 , in view of (3.18). Then, using (3.32) as before, we conclude
that the right-hand side of (6.14) is !1 C 12 l.P /.
A straightforward combination of our conclusions above proves the following
proposition.
Proposition 6.7. Let P 2 E be such that u.P / u and v.P / > 0.
66
Assume that P 2 E is such that P C D .u.P /, v.P // has integer coordinates with
v.P / > 0 and u.P / u and write P E as in (4.2). With l.P / as in (4.6), define
L.P / by
(6.15)
67
p
Since v.P / > 0, we have v.P / D Q.u.P // and then, in view of (6.4) and (6.9),
x.P / D x.u.P //. Let x.P / D m=n, where m, n are relatively prime integers. But
also, x.P / D .x.u.P // 3u.P /2 /=.3u.P /2 / where the numerator and denominator
are integers; actually, the numerator is equal to c u.P /2 C3d u.P /C6e 2 C6e v.P /.
Therefore, jmj 3jx.u.P //j u.P /2 3c7 u.P /2 , jnj 3u.P /2 3c7 u.P /2
and, consequently, for the logarithmic height of the rational number x.P / we have
(see (2.21))
h.x.P // D log maxjmj, jnj log.3c7 u.P /2 / D log.3c7 / C 2 log u.P /. (6.17)
Then,
1
1
log.3c7 / h.x.P //
(by (6.17)/
2
2
1
O /
log.3c7 / C h.P
(by Proposition 2.6.3)
2
1
(by Proposition 2.6.2).
log.3c7 / C M 2
2
Now, a straightforward combination of the last inequality with (6.16) proves (6.15).
log.u.P //1
Remark. The value of c7 given in the proof was obtained under the assumption that u
is so large that Q.u/ 2au4 . Obviously, this choice is arbitrary; we could have chosen
any constant > 1 and then assumed that u is so large that Q.u/ constant u4 .
Solutions .u, v/ with u < 0. On the first glance, Theorem 6.8 seems to treat integer
solutions .u, v/ to (6.1) with large positive u, only. However, as we will show below,
with a very little extra effort we can treat all integer solutions .u, v/ with large juj
and obtain an exactly analogous upper bound in which the assumption concerning
the large unknown point P C D .u.P /, v.P // is that ju.P /j is sufficiently large
rather than u.P / is sufficiently large.
It is clear what one should do; it suffices to replace the coefficients a, b, c, d , e of
the given equation by a, b, c, d , e, respectively, and consider the equation v 2 D
au4 bu3 C cu2 du C e 2 in integer .u, v/ with sufficiently large u > 0. How does
the replacement of .a, b, c, d , e/ by .a, b, c, d , e/ affect the crucial quantities (=
parameters and auxiliary functions) involved in the discussion so far? This is shown in
Table 6.1. The symbols of the quantities involved in our discussion are on the left,
while on the right the analogous symbols are listed with a bar on them, denoting
the analogous quantities which result when the coefficients b, d are replaced by
b, d , respectively. We avoid the bar indication if the replacement does not affect
the quantity.
68
Table 6.1. Parameters and auxiliary functions for the solution of the quartic elliptic equation.
0 10 Q.u/ D au4 C bu3 C cu2 C du C e 2
A, B (see (6.3))
x0 (see (6.8))
if d a C be 0
p
D
if d a C be D 0
A, B
x0
x.u/ D
x.u/ D
Choose u0 1 so that
Choose u0 1 so that
Choose u u0 so that
In the case that x0 e1 :
u > u ) x.u/ > e1
Choose u u0 so that
In the case that x0 e1 :
u > u ) x.u/ > e1
P0E if d a C be 0
E
p
P0 D
if d a C be D 0
P0E
p
l.P0 / if d a C be 0
p
l.P 0 / D
l.P0 /
if d a C be D 0
a//
l.P0 /
L.P /
if either x0 > e1
or e2 , e3 2 R and x0 2 .e3 , e2 /
then L.P / D l.P / l.P0 /
L.P / D l.P / l.P0 /
if x0 D e1 then
L.P / D l.P / C 2!1
L.P /
if either x0 > e1
or e2 , e3 2 R and x0 2 .e3 , e2 /
then
p
l.P / C l.P0 / if d a C be 0
p
L.P / D
l.P / l.P0 / if d a C be D 0
if x D e then
0
1
L.P / D l.P / C 2!1
To continue with the resolution of the quartic elliptic equation, one can safely avoid
Chapters 7 and 8 and go directly to Chapter 9, Sections 9.1, 9.2 and 9.4.
Chapter 7
W 2 D2 V 2 D A2
(7.1)
of simultaneous Pell equations, where .U , V , W / are positive integers, under the assumption that we know a rational solution .U0 , V0 , W0 / with V0 > 0 and U0 , W0 0,
not both zero. This assumption is not really restrictive in practice. All known examples possess at least one small solution, easily found by a direct search; see e.g. the
table in [60, Appendix]. An alternative method to discover a rational, or even integral,
solution .U0 , V0 , W0 / is to consider the equation .D1 V 2 C A1 /.D2 V 2 C A2 / D Z 2
(Z D U W ) and search for rational solutions. For example, if we take .D1 , A1 / D
.7, 2/ and .D2 , A2 / D .32, 23/ (an example from Z. Y. Chen, mentioned in [60, Appendix]), then the MAGMA routine Points applied to the corresponding curve1 immediately returns the points .V , Z/ D .1, 9/, .271, 1099161/, which furnish the solutions
.U , V , W / D .3, 1, 3/, .717, 271, 1533/.
In [60] an alternative method for the explicit resolution of (7.1) is developed, which
requires the solution of a number of quartic Thue equations; moreover, that paper
contains a rich bibliography on simultaneous Pell equations, to which we refer the
interested reader. Anyway, choosing between Ellog and Thue equations for the resolution of a certain Diophantine problem does not admit a general answer; this issue is
discussed in the beautiful paper [58].
Now, let us come back to (7.1). To each solution .U , V , W / we will associate a
point P on a certain short Weierstrass model, along with a linear form L.P / and then
we will prove Theorem 7.1, which will provide us with an explicit upper bound for
jL.P /j (cf. page 56 immediately after (4.8)). This chapters exposition is based on the
authors paper [63].
Throughout the chapter we keep the notations, assumptions and results of Chapter 4.
Clearly, for a solution .U0 , V0 , W0 / as above,
p
p
(7.2)
sgn.U0 V0 D1 / D sgn.A1 /, sgn.W0 V0 D2 / D sgn.A2 /.
Moreover, we will assume that A1 D2 A2 D1 0, otherwise the system is rather
trivially solved. Indeed, if A1 D2 A2 D1 D 0, then (7.1) implies D2 U 2 D D1 V 2 ,
so that A2 =A1 D D2 =D1 D q 2 for some q 2 Q and, consequently, the solutions
1
70
U 0 V V0 U
,
W0 V C V0 W
(7.4)
(7.5)
so that
V0 .U C uW /
.
(7.6)
U0 uW0
Then we substitute for V in (7.4) to obtain a quadratic equation in U=W . One solution
of that equation is, obviously, U=W D W0 =U0 , hence the other solution is
V D
A2 U0 u2 2A1 W0 u C A1 U0
U
.
D
W
A2 W0 u2 2A2 U0 u C A1 W0
(7.7)
A2
,
80W 2
(7.8)
U02 D1 V02 D A1 ,
W02 D2 V02 D A2 ,
(7.9)
imply that
def
A2 W0 u2 2A2 U0 u C A1 W0
,
W
p
D sgn.A1 / D sgn.U0 V0 D1 /,
vD
(7.10)
(7.11)
(7.12)
.a, b, c, d / D .A22 , 4A2 U0 W0 , 4A1 W02 2A1 A2 C 4A2 U02 , 4A1 U0 W0 , jA1 j/.
(7.13)
Equation (7.10) has been studied in Chapter 6. There we focused our interest on integer
solutions .u, v/ with arbitrarily large u > 0, whilst now the u and v in (7.5) and (7.11),
71
respectively, are certainly rationals but, in general, not integral and, moreover, as we
will see, u belongs to a finite interval instead of approaching infinity. Nevertheless, we
will exploit part of the results of our study in Chapter 6 in order to solve our present
problem.
In (7.11) we replace u by its expression in (7.5) in order to express v as a rational
function of U , V , W (we also take into account (7.9)):
A2 U0 U C .D2 A1 D1 A2 /V0 V C A1 W0 W
.
(7.14)
v D 2V02
.V0 W C W0 V /2
Conversely, we need also express U , V , W as rational functions of u, v. By the definition of v,
A2 W0 u2 2A2 U0 u C A1 W0
W D
(7.15)
v
and this, combined with (7.8) and (7.7) gives
A2 U0 u2 2A1 W0 u C A1 U0
.A2 u2 C A1 /V0
, U D
.
(7.16)
V D
v
v
Now, using (6.4), (6.5) and (6.6), we establish a birational transformation between
(7.10) and
def
(7.17)
E : y 2 D x 3 C Ax C B D f .x/,
with A and B expressed by the formulas (6.3) as functions of a, b, c, d , e taken from
(7.13).
R
As in Chapter 6, we will express an elliptic integral pdt , with appropriate upf .t/
per and lower limits, in terms of another integral in which the variable V is directly
involved. This will be accomplished with the aid of a function x.V /, defined on an interval V , C1/, for a sufficiently large positive V . Below we proceed to construct
this function.
V0 , W0 / be any positive real solution to (7.1), so that U D
p Let .U , V , W / .U0 ,p
A1 C D1 V 2 and W D A2 C D2 V 2 . Then (7.5) and (7.11) suggest that we define
the functions
p
U0 V V0 A1 C D1 V 2
(7.18)
V 7! u.V / D
p
W0 V C V0 A2 C D2 V 2
and
A2 W0 u.V /2 2A2 U0 u.V / C A1 W0
V 7! v.V / D
p
,
(7.19)
A2 C D2 V 2
p
A1 C D1 V 2
which satisfypthe relation v.V /2 D Q.u.V //. The substitutions U
2
and W
A2 C D2 V in (7.5) makes u D u.V /, v D v.V / (by (7.15)), and
consequently, in view of (7.14),
p
p
2
2
2 A2 U0 A1 C D1 V C .A1 D2 A2 D1 /V0 V C A1 W0 A2 C D2 V
v.V / D 2V0
p
.
.W0 V C V0 A2 C D2 V 2 /2
(7.20)
72
After some standard computations, in which use is made of the definition of (see
(7.12)), as well as of the relations (7.2) and (7.3), we verify the following (dashes
denote derivatives):
p
U0 V0 D1
def
p .
(a) u1 D lim u.V / D
V !C1
W0 C V0 D2
p
p
p
. D2 U0 C D1 W0 /.U0 V0 D1 /V0
3
0
.
p
p
(b)
lim V u .V / D
V !C1
D1 D2 .W0 C V0 D2 /
(c)
lim v.V / D 0.
V !C1
p
p
p
2 jU V
D
j.U
D
C
W
D1 /
2V
0
0
1
0
2
0
0
(d)
lim V 2 v0 .V / D
p
.
V !C1
W0 C V0 D2
(e) Q.u1 / D 0.
p
p
p
p
8V03 D1 D2 .U0 V0 D1 /.U0 D2 C W0 D1 /
p
.
(f) Q .u1 / D
W0 C V0 D2
We note first that, by (e), the polynomial Q has four real roots which,
p in view
p of
(a), are actually the algebraic conjugates of u1 in the quartic field Q. D1 , D2 /;
temporarily, let us denote these roots by r4 < r3 < r2 < r1 .
If D 1, then, by (f), Q is strictly increasing in an interval with u1 as its centre,
therefore, u1 D r1 or r3 and Q.u/ > 0 in an interval I with u1 as left end point
and right end point appropriately close to u1 (or even C1 if u1 D r1 ). If D C1,
then an analogous argument shows that u1 D r2 or r4 and Q.u/ > 0 in an interval I
with u1 as right end point and left end point appropriately close to u1 (or even 1
if u1 D r4 ).
Then, with the interval I as above, it is clear that
p
Q.u/ C cu2 C 3du C 6e 2
def 6e
,
(7.21)
I 3 u 7! x.u/ D
3u2
p
ebu3 C 2ecu2 C 3edu C 4e 3 C .du C 4e 2 / Q.u/
def
I 3 u 7! y.u/ D
(7.22)
u3
are meaningful as real-valued functions; note that the formula for x.u/ is the same as
that in (6.9), but the domains of definition differ: There, the domain is a neighbourhood
of C1, here the domain is thepinterval I . We also remark that, if in the right-hand
sides of (6.4) we replace v by Q.u/, then we obtain the right-hand sides of (7.21)
and (7.22), respectively, except that is defined differently. Note, however, that a
different choice of in (6.4), (6.5) and (6.6) does not affect the fact that these relations
establish a birational transformation between v 2 D Q.u/ and y 2 D x 3 C Ax C B;
cf. the comment just before the announcement of Lemma 6.1.
0
73
Next we note that, in view of (c) and (d) on page 72, if V is sufficiently large, then
v.V / > 0.
Having in mind the above remarks, let .u, v/ D .u.V /, v.V //. In view of (a) and
(b) on page 72, if V is sufficiently large, then u D u.V / 2 I and v Dpv.V / > 0.
Since v 2 D v.V /2 D Q.u.V // D Q.u/, we conclude then that v D Q.u/
p and,
consequently, from (7.21) and (6.4) we see that x.u/ D X .u, v/
p D X .u, Q.u//.
Similarly, (7.22) and (6.4) imply that y.u/ D Y.u, v/ D Y.u, Q.u//. Therefore,
for sufficiently large V , the functions
def
V 7! x.V / D .x u/.V /,
def
V 7! y.V / D .y u/.V /
v.V / D
(7.23)
p
Q.u.V // ;
x0 D
D
and
V !C1
4V02
p
Q.u1 //
(7.24)
p
D2 U02 C D1 D2 V02 C D1 W02 C 3 D1 D2 U0 V0
p
p
C 3D1 D2 V0 W0 C 3 D1 D2 U0 W0 > 0
p
lim y.V / D Y.u1 , Q.u1 //
(7.25)
V !C1
p
p
p
p
p
D 8V03 D1 D2 .W0 D1 C U0 D2 /.U0 C V0 D1 /.W0 C V0 D2 / > 0.
def
y0 D
74
p
Note that .x0 , y0 / is a point of E.R/; more precisely, a point of E.Q. D1 D2 //.
Indeed, look at (7.24) and (7.25) and, once again, remember that any .x, y/ D
.X .u, v/, Y.u, v// with real u, v and v 2 D Q.u/ is a point of E.R/. The point
p
def
P0E D .x0 , y0 / 2 E.Q. D1 D2 //
(7.26)
plays a role in our method; see Theorem 7.1 below. In the applications we will find
important the property that 2P0E 2 E.Q/; more specifically,
2P0E D . 43 .D1 D2 V02 C D1 W02 C D2 U02 /V02 , 8D1 D2 U0 W0 V04 /.
(7.27)
Coming back to x.V / and y.V /, we conclude that, according to our above discussion, if V 2 .V , C1/ for a sufficiently large V , then x.V / 2 I.x0 /, where
I.x0 / is an interval .x1 , x0 / or .x0 , x1 /, according to whether D C1 or D 1,
respectively, for an
p appropriate x1 . Moreover, for V 2 .V , C1/ we have y.V / > 0,
whence y.V / D f .x.V //.
Rx
Now we focus our interest on the integral x 0 pdt when x belongs to the interval
f .t/
I.x0 / as above. According to our previous discussion, the interval with end points x
(included) and x0 (excluded) is the image under the function x D x u of an interval
V , C1/, where x D x.V /. Therefore, for t 2 x, x0 / we can make p
the change of
variable
t D x.T / D x.u.T //, where T 2 V , C1/. Then y.T / D f .x.T // D
p
f .t / and
dx du
dt
D
.
dT
du dT
From (7.21) we calculate
p
eb u3 C 2ec u2 C 3ed u C 4e 3 C .d u C 4e 2 / Q.u/
dx
p
D
du
u3 Q.u/
p
f .t /
y
y.u/
.
D D
D (by (7.22)) p
v
v
Q.u/
From (7.18) we compute
p
p
du
V0 .A2 U0 A1 C D1 T 2 C .A1 D2 A2 D1 /V0 T C A1 W0 A2 C D2 T 2
p
p
p
D
dT
.V0 A2 C D2 T 2 C W0 T /2 A1 C D1 T 2 A2 C D2 T 2
v
D (by (7.20))
p
.
2V0
.A1 C D1 T 2 /.A2 C D2 T 2 /
On combining the above relations we obtain
p
f .t /
dt
1
p
D
dT
2V0
.A1 C D1 T 2 /.A2 C D2 T 2 /
and, consequently,
Z C1
Z x0
1
dt
dT
D
.
p
p
2V0 V
x.V /
f .t /
.A1 C D1 T 2 /.A2 C D2 T 2 /
(7.28)
75
Now we examine more closely the short Weierstrass model E defined by (7.17).
First we note that the polynomial f .X / has the following three rational roots:
4
4
1 D .D2 X02 2D1 D2 Y02 C D1 Z02 /Y02 D .A1 D2 C A2 D1 /Y02
3
3
4
4
2 D .2D2 X02 C D1 D2 Y02 C D1 Z02 /Y02 D .A2 D1 2A1 D2 /Y02
3
3
4
4
3 D .D2 X02 C D1 D2 Y02 2D1 Z02 /Y02 D .A1 D2 2A2 D1 /Y02 .
3
3
Note that, with our ordinary notation, the roots of f .X / are denoted by e1 , e2 , e3 , where
e1 > e2 > e3 ; therefore, .e1 e2 e3 / is a permutation of .1 2 3 /. Now we consider
the point P0 with P0E D .x0 , y0 / 2 E.R/ defined in (7.26). We easily check that
x0 i > 0 for i D 1, 2, 3, therefore, P0E 2 E0 .R/ and, consequently, according to
our comment immediately after (7.25), we may assume that, if V is sufficiently large,
then the point PV , which is defined by P E D .x.V /, y.V //, belongs to E0 .R/ and
y.V / > 0. Then, we can write the left-hand side of (7.28) as follows:
Z C1
Z C1
Z x0
dt
dt
dt
p
p
p
D
D 2.l.P / l.P0 // (by (3.32)).
x.V /
x.V /
x0
f .t /
f .t /
f .t /
Easily, the right-hand side of (7.28) is absolutely less than apconstant times V , provided
that V is not very small. For example, if V maxp
iD1,2 2jAi j=Di , then the right1
hand side of (7.28) is absolutely less than V =.V0 D1 D2 / and, consequently,
jl.P / l.P0 /j <
V 1
.
p
2V0 D1 D2
(7.29)
u.V / (cf. (7.18)) and v.V / (cf. (7.19)) satisfy Q.u.V // > 0 and v.V / > 0, respectively.
V belongs to an infinite interval .V , C1/ in which the functions x and y, defined
in (7.23), have the following properties: The function x is strictly monotonous, x.V /
is larger than the maximum of the i s displayed at page 75 and y.V / > 0.
p
p
For example, this constant can be chosen to be 2 Di if Ai > 0 and Di if Ai < 0.
76
Let E be the Weierstrass model defined in (7.17) and consider that, everything in
Chapter 4, in particular the relations (4.1) and (4.4)(4.8), refers to this particular
model E. Consider also the points P E D .x.V /, y.V // 2 E.Q/ and P0E defined in
(7.26). Express P E as in (4.2) and l.P / by (4.6), and set L.P / D l.P / l.P0 /. Then,
jL.P /j
1
p
exp. 12 c8 C M 2 /,
2V0 D1 D2
(7.30)
where 12 0 and are as in Propositions 2.6.2 and 2.6.3, respectively, and c8 is a positive constant which can be explicitly calculated as explained below, in the paragraph
immediately before relation (7.32).
Proof. Firstly, by what we have already seen on pages 73 and 74, .x.V /, y.V // and
.x0 , y0 / are, indeed, points of E.R/. Further, we saw on page 75 that .x0 , y0 / 2 E0 .R/
and, if V is sufficiently large,
then also .x.V /, y.V // belongs to E0 .R/. More prep
cisely, .x0 , y0 / 2 E0 .Q. D1 D2 // and, since U , V , W are integers, .x.V /, y.V // 2
E0 .Q/.
Next, a quick look at our previous discussion after page 71 and before the relation
(7.29) easily reminds that, indeed, for sufficiently large V all conditions marked by
bullets are fulfilled, so we proceed to the main part of the proof, noting the following:
If u and v are given by (7.5) and (7.14), respectively, then
x.V / D
(7.31)
hence
log V 1
1
1
c8 h.x.V //
2
2
Here we need the requirement set by the first bullet in the announcement of the theorem.
(7.32)
77
for some explicitly computable positive constant c8 ; in the notation of the previous
lines,
c8 D log maxconst1 , const2 .
O / 1 h.x.V // ; by Proposition 2.6.2, h.P
O / M 2 and
By Proposition 2.6.3, h.P
2
combining these with (7.32) gives
O / 1 c8 C M 2 .
log V 1 12 c8 12 h.x.V // 12 c8 C h.P
2
p
By the definition of L.P / and (7.29) we have L.P / < V 1 =.2V0 D1 D2 / and, combining this with the last displayed inequality, gives the inequality (7.30).
To continue with the resolution of the simultaneous Pell equations, one can safely
avoid Chapter 8 and go directly to Chapter 9, Sections 9.1, 9.2 and 9.5.
Chapter 8
(8.1)
where g is a polynomial with coefficients in Z, irreducible over Q and the curve represented by the model (8.1) is of genus one and possesses a non-singular point with
rational coordinates.
Throughout the chapter we keep the notations, assumptions and results of Chapter 4.
This chapter is based mainly on the joint work [56] of R. Stroeker and the author.
It is well known that the equation (8.1) has finitely many solutions (see, for example, [2] or [41]), but here, in analogy to Chapters 5, 6 and 7, we intend to prove
Theorem 8.7.2, which gives an explicit upper bound for jL.P /j, where L.P / is an
appropriate linear form in elliptic logarithms; cf. page 56 immediately after (4.8).
Without loss of generality we assume that degu g degv g D n. We will make the
further assumption that C.R/ possesses infinite (= unbounded) branches along both
the u-direction and the v-direction. This is very natural to assume, for, otherwise, the
problem of solving g D 0 in integers is, obviously, trivial.
Before proceeding to the study of (8.1), we collect all the assumptions that we set
so far, which will be implicitly understood throughout the whole chapter.
The curve represented by the model (8.1) is of genus one and a non-singular point
with rational coordinates is known.
degu g degv g D n.
C.R/ possesses infinite branches along both the u-direction and the v- direction.
A, B 2 Q.
(8.2)
79
Proof. We refer to van Hoeijs [22], where an algorithm is presented for calculating
the birational transformation between C and a short Weierstrass model E, provided
that a non-singular point .u0 , v0 / of C with algebraic coordinates is known. According
to Section 3.2 of that paper, the coefficients of the transformation which is constructed,
as well as the coefficients of the Weierstrass model E, belong to the field Q.u0 , v0 /;
this is even more clear if one looks at van Hoeijs papers [20] and [21] on which the
algorithm of [22] is based. Since we have assumed that a non-singular point of C with
rational coordinates is known, it follows that, if we insert this particular point in van
Hoeijs algorithm [22], we will get a birational transformation with rational coefficients between C and a short Weierstrass model E which will be defined over Q.
The birational transformation of Proposition 8.1.1 will be denoted by
C 3 .u, v/
.X ,Y/
!
.U,V/
.x, y/ 2 E ,
(8.3)
meaning that
u.P / D U.x.P /, y.P //,
(8.4)
With the aid of MAPLE we find out that the genus of the curve C : g D 0 is one. Obviously, u D .0, 0/ is a point, so that C is a model of an elliptic curve over Q. The MAPLE
implementation of van Hoeijs algorithm [22] gives the birational transformation (8.3)
between C and the model
120408061761
E : y 2 D f .x/ D x 3 17846163x C
.
(8.5)
4
The functions X and Y are the following:
X .u, v/ D 2710u 9225uv 3 43821uv 5904u2 C 21867uv 2 C 726v 4 C 2169u3
C 8883u2 v 2 48429vu2 C 2349v 4 u/=u.u 1/2
Y.u, v/ D 79947u4 C 14069376vu3 829440u3 v 2 C 1802709u3 C 137565u2
6952716v 3 u2 C 140667858vu2 39726990u2 v 2 1224558v 4 u2
22481070uv 2 4683906v 4 u C 11958534uv 3 C 34375266uv
C 81029u 395286v 4 /=.u.u 1/3 ,
(8.6)
80
(8.7)
1
X
(8.8)
kD i
We will need them only at the final stage of the resolution, in Section 10.2.4.
We note that the terminology Puiseux series is absent from these two books.
81
n
Y
.v vi .u//
iD1
n
v in
g.u, v/.
(d) Each series (8.8) converges for juj > B0 , where B0 is the maximum modulus of
the roots of the polynomial Resv .g, @g
/ 2 Zu.
@v
P
k
(e) For each i the function t 7! vi .t i / D 1
kD i k,i t is analytic and one-to1=i
one into the punctured disk with centre at the origin and radius B0
Proof. We do not present here a proof in the strict sense of the word. The proofs of
statements (a) through (d) are scattered in the bibliography.
Statements (a), (b), and (c) are classical results about Puiseux series. They can be
found in classical books such as [4, Chapter II] and [66, Chapter IV], though in slightly
different form. In these books the notion of parameterisation is used in order to express
the solutions .u, v/. More specifically, instead of (8.8), the parameterisation
u D t i
v D vi .t i / D
1
X
k,i t k
kD i
1
X
kD i
.k,i / mk uk=i ;
82
in particular,
vi .u, 0, id/ D vi .u/
.i D 1, : : : , n/.
The series vi .u, m, / also satisfies g.u, vi .u, m, // D 0, therefore, by Fact 8.2.1(b),
it must coincide with another series (8.8), say with vj .u/ D vj .u, 0, id/. We then say
that vj .u/ and vi .u/ belong to the same conjugacy class. The n series (8.8) are thus
partitioned into disjoint conjugacy classes (see after relation (3.5) in [67]). Therefore,
the computation of a smaller set of series (8.8) composed of representatives of each
conjugacy class is sufficient for the computation of all Puiseux series (8.8).
ExampleStep 2 . (continued from the end of Step 1, page 80)
With the aid of MAPLE we compute two conjugacy classes of Puiseux series solving
g.u, v/ D 0. The first class consists of the series
1 1105 1=2
6529121 3=2 390266 2
v1 .u/ D i u1=2 C
iu
iu
u C ,
138u1 C
C
2
24
1152
9
p
where i D 1, and
1 1105 1=2
6529121 3=2 390266 2
138u1
iu
iu
u C ,
v2 .u/ D i u1=2 C C
2
24
1152
9
The second class consists of the series10 5
1
1271288 4=3
275 1=3 3317 2=3
v3 .u/ D u1=3 C
C 92u1
u
u
u
3
9
81
729
9364963 5=3
C ,
u
6561
1 275 1=3 3317 2 2=3
1271288 4=3
C 92u1
C
!u
! u
!u
3
9
81
729
9364963 2 5=3
! u
C ,
6561
v4 .u/ D ! 2 u1=3
Computing the number field K . The extension K=Q mentioned in Fact 8.2.1(a),
can be computed by an algorithm of J. Coates which is implicitly contained in the
proof of [8, Lemma 3].
83
n
X
pj .u/v nj .
j D0
For any i D 1, : : : , n it suffices to show that there exists an index k0,i 0 such that,
kC 2 Q. , : : : , Ck0,i / for all k k0,i . Therefore, let us fix i 2 1, : : : , n and
subsequently, in order to simplify notation, omit the subscript i from
i , k,i , i and
vi .u/ in (8.8). Compute the integers D i and
D
i by the first few steps in the
construction of the Newton polygon.
Choose a non-negative integer N such that
degpj C
.n j / C N 0 .j D 0, : : : , n/
with equality for at least one subscript j , and put
Pj .X / D pj .X / X .nj /CN
and
G.X , Y / D
n
X
.j D 0, : : : , n/
Pj .X / Y nj .
j D0
1
X
kC X k .
(8.9)
kD0
Before giving Coates algorithm, let us make a notation remark: For any polynomial h
over any field, in one or more variables X , : : :, when we write r D ordX .h/, we mean
that X r appears in some monomial of h with non-zero coefficient, and r is minimal
with respect to this property.
10 6 Computing the finite extension Q. , C1 , : : :/
DESCRIPTION: Compute k0 0 such that Ck 2 Q. , : : : , Ck0 /
INPUT: G.X , Y / as above, m D degX G,
D .2n 2/m C 1.
OUTPUT: k0 .
INITIAL VALUES: G0 .x, y/
G.x, y/, k
0.
Hk .Y /
Gk .0, Y /.
if degHk D 1 then go to ()
else
Define uk by Hk .uk / D 0
r
ordX .Gk .X , uk C X Y //
X r Gk .X , uk C X Y /
Gk .X , Y /
k
kC1
goto ()
endif
() k0
k
()
END
8k k0 .
84
By the proof of [8, Lemma 3], the polynomials Hk .k D 0, 1, 2, : : : / are not identically zero, their degrees are in non-increasing order and degHk0 D 1 for some k0
,
so that for all k k0 the polynomial Hk is of degree 1. Now for each k 0, the
algebraic number kC in (8.9) is a root of Hk , and therefore assumes one of the
possible values of uk . By the linear character of Hk for all k k0 , it follows that
for k > k0 , the coefficient kC is uniquely determined by the previous s and
kC 2 Q. , C1 , : : : , Ck0 /.
@g
/ D .1/n.n1/=2 p0 .u/ discv .g.u, v//;
@v
(8.10)
85
Notations and assumptions. From now on and until the end of the chapter, E will
denote the elliptic curve, a representative of which is the model C defined in (8.1),
another one being the model E referred to in Proposition 8.1.1. We will denote by P
the generic point of E such that the coordinates of P C D .u.P /, v.P // are integers
and u.P / > B0 ; later we will require that u.P / be larger than an explicit constant
larger than B0 . Then, by Lemma 8.3.1, v.P / D vs .u.P // for some s 2 1, : : : , n,
where vs is as in Fact 8.2.1(a), with all its coefficients real algebraic numbers.
For a positive real number u, when we write uk=s we mean the unique s -th real
root of uk .
Referring to the birational transformation (8.3), we have the following proposition:
Proposition 8.3.2. The limit limu!C1 X .u, vs .u// exists in R [ 1. If this limit
is 1, then it is a real algebraic number that can be explicitly computed.
Next, define
x.u/ D X .u, vs .u//, y.u/ D Y.u, vs .u//.
Then, there exists a B1 B0 such that, in the interval .B1 , C1/ the functions x and
y are continuous, x is strictly monotonous and y does not change sign.
Proof. First note that g.u, v/ cannot be a factor of either the numerator or the denominator of the rational function X .u, v/. For, otherwise, the model C could be injectively
mapped into a straight line, which is impossible for a curve of genus 1.
1=
Next, put u D t s with t 2 R and 0 < t < B0 s . Then, by Lemma 8.3.1 and
Fact 8.2.1(a),
x.u/ D X .t
s
1
X
kD s
k,s u
k=s
00
t
C 0 t
C 00 t
C : : :
D x.t /, (8.11)
/D
t C 0 t 0 C 00 t 00 C : : :
if D ,
<=
lim x.u/ D 0
if > ,
u!1
:
sgn.= /1 if < .
By Proposition 8.1.1, the rational function X is explicitly computable and the same is
true for the k,s due to the algorithm on page 83. It follows that , and x0s can be
explicitly computed.
Next, since the polynomial g does not divide neither denominator of X or Y, there
exists a positive constant M1 B0 such that, if g.u, v/ D 0 and u > M1 , then
.u, v/ is not a zero of neither denominator of X or Y. This implies that in the interval .M1 , C1/, the functions x and y are continuous, because X and Y are rational
functions and vs is continuous (cf. Fact 8.2.1(e)).
86
Concerning the strict monotonicity of the function x near C1, it suffices to show
the strict monotonicity of the right-hand side of (8.11) near 0C . But this function in t is
0
00
of the form t C 0 t C 00 t C , where , 0 , 00 , : : : are non-zero real numbers
and D
<
0 <
00 < : : :. As t ! 0C the derivative of this function is
non-zero with constant sign, and this proves our claim.
0
00
In complete analogy, we can express y as a series t C 0 t C 00 t C , where
0
00
0
00
, , , : : : are non-zero real numbers and
<
<
, : : :. Sufficiently close to 0C ,
such a series, clearly, does not change sign.
Remark. In specific numerical examples the constant B1 can be explicitly calculated,
as described in Section 8.5.
Now we define a point P0 that will be involved in the linear form L.P / of the main
result of this chapter, namely, Theorem 8.7.2.
def
Definition 8.3.3. Let x0s D limu!1 X .u, vs .u// and denote by P0 D P0s the point
of E with
87
C1
u.P /
g.u/ du
D"
gv .u, vs .u//
x0s
x.P /
dx
p
.
f .x/
(8.12)
and
(8.13)
3 x.u/2 C A 0
x .u/.
y0 .u/ D " p
2 f .x.u//
(8.14)
The relation g.u, vs .u// D 0 implies gu .u, vs .u// C gv .u, vs .u//vs0 .u/ D 0. Solving
for vs0 .u/ and substituting in the relation which results on equating the right-hand sides
of (8.13) and (8.14), gives
p
f .x.u//
0
.
(8.15)
x .u/ D " g.u/
gv .u, vs .u//
In view of the strict monotonicity of the function x we can make the change of variable
in the integral below
Z
Z C1
Z x0s
dx xDx.u/ C1 x0 .u/ du .8.15/
g.u/ du
p
D
p
D "
,
u.P /
u.P / gv .u, vs .u//
x.P /
f .x/
f .x.u//
as claimed.
For Ellog it is important that the integrand in the left-hand side of (8.12) tends to
zero as u ! C1. This is made precise in the proposition below.
Proposition 8.4.2. There exists a sufficiently large constant B2 B1 such that, if
u B2 , then
g.u/
1
,
(8.16)
g .u, v .u// c9 juj
v
s
where s1 and c9 is a positive constant independent from u.
88
hence,
x0 .t /
dt D 2 C At 4 C Bt 6 C O.t 8 /
y.t /
R
and we see that there are no singularities in the expansion of dx=y. This, by definition, means that the elliptic integral
is of the first kind (see [4, 24]). Therefore, the
R
same is true for the integral .g.u/=gv .u, v// du, implying that, for any parameterisation .u, v/ D .u.t /, v.t // of g.u, v/ D 0, the t -expansion of
dx
D
y
g.u/
du
gv .u.t /, v.t // dt
has no negative powers of t . Applying this conclusion to the parameterisation
u.t / D t s ,
v.t / D vs .t s / D
1
X
k,s t k ,
kD s
and taking into account that du=dt D s t s 1 , we are led to the inequality
ord t
g.t s /
s C 1.4
gv .t s , vs .t s //
D t 1Cs C .0 C 1 t C 2 t 2 C /,
0, 0 0.
D .1 C /=s 1=s .
For sufficiently large u, the absolute value of the right-hand side is bounded above by
c9 juj1
, as claimed.
Remark. In specific numerical examples the constants B2 , c9 and of Proposition 8.4.2
can be made explicit according to Section 8.6.
4
The function ord has been defined just before the algorithm on page 83.
8.5
89
In Proposition 8.3.2 we introduced the constant B1 , the existence of which is guaranteed by a non-constructive proof. In any specific equation, however, we must explicitly
compute such a B1 , and in this section we show how we can do this. First we state a
simple lemma.
Lemma 8.5.1. Let F be a polynomial in two variables with real coefficients and let
V : R ! R be a continuous function, such that F .u, V .u// D 0 for juj > U0 ,
where U0 is a positive constant. Let R be the set of all real roots of the polynomial
F .0, Y /, and S D u : juj > U0 & V .u/ 2 R. Put .Umin , Umax / D .U0 , U0 /
or .Umin , Umax / D .min S, max S/ according to whether S is or is not empty. Then
the function V keeps a constant sign in the interval .Umax , C1/ and does so in the
interval .1 , Umin /.
Proof. Suppose that u2 > u1 > Umax and V .u2 /V .u1 / < 0. Then, by the continuity
of V , there exists a u0 2 .u1 , u2 /, such that V .u0 / D 0. This means that u0 2 S , hence
u0 Umax . But u0 > u1 > Umax , a contradiction.
Similarly we arrive at a contradiction if we assume that u2 < u1 < Umin and
V .u2 /V .u1 / < 0.
We will use this lemma in order to show how to compute M1,min and M1,max absolutely larger than B0 , such that y keeps a constant sign in the interval .M1,max , C1/
and does so in the interval .1 , M1,min /.
Let D.u, v/ be the square-free part of the product of the denominators of X and Y.
Since g is an absolutely irreducible polynomial and does not divide neither denominator of X or Y, the resultant Resv .D, g/ is a non-zero polynomial in u; we denote
by R the maximum absolute value of its real roots. Let any u with juj > maxB0 , R.
After clearing out the denominator in y.u/ D Y.u, vs .u//, we obtain the polynomial
relation H1 .u, vs .u/, y.u// D 0. On the other hand, we also have g.u, vs .u// D 0 and
we can eliminate vs .u/ from the last two relations, using resultants. We thus obtain
a polynomial relation R1 .u, y.u// D 0. Now, applying Lemma 8.5.1 with F D R1
and V D y we compute values Umin , Umax (in the notation of the lemma) and we take
M1,max D dmaxUmax , B0 , Re and M1,min D bminUmin , B0 , Rc.
Next we show how to compute M2,mi n and M2,max such that in both intervals,
.M2,max , C1/ and .1 , M2,mi n /, x is strictly monotonous.
Again, let any u with juj > B0 . Clearing out the denominator in x.u/ D
X .u, vs .u// gives a polynomial relation H2 .u, vs .u/, y.u// D 0. From this relation
and g.u, vs .u// D 0 we eliminate vs .u/, using resultants. Thus we obtain a polynomial relation R2 .u, x.u// D 0. Differentiating with respect to u gives
@R2
@R2
.R2 .u, x.u// C
.R2 .u, x.u// x0 .u/ D 0,
@u
@x
(8.17)
90
91
cients with at most 33 decimal digits; also, degu .R30 / D 10 and degx0 .R30 / D 5.
If we assume that juj 540, then, necessarily, R30 D 0 and we apply Lemma
8.5.1 with F D R30 , V D x0 , U0 D B0 D 484. Now F .0, x0 / D .271x0 299/
times a large constant, therefore R D 299=271. The approximate solutions of
R30 .x, 299=271/ D 0 are 932.85436, 738.24160, 173.11938 and 0, therefore,
S D 932.85436 : : : , 738.24160 : : :. Consequently, M2,mi n D 933, M2,max D
484 and Lemma 8.5.1 implies that x0 is non-zero in the interval .1 , 933/ and has a
constant sign, hence x is strictly monotonous. Similarly, x is strictly monotonous in the
interval .484 , C1/. Finally, by the displayed formula for B1 we obtain B1 D 2340.
Therefore, from now on, we will assume that juj > B1 D 2340.
Example continued on page 91.
8.6
At first we work as in Section 8.5. We denote by I D I.u/ the left-hand side of (8.16)
and clear out the denominator to obtain a polynomial relation H4 .u, vs .u/, I/ D 0.
Using resultants we eliminate vs .u/ from the last relation and g.u, vs .u// D 0. Thus
we obtain a polynomial relation
I m C q1 .u/ I m1 C C qm1 .u/ I C qm .u/ D 0,
(8.18)
92
where H40 .u, vs .u/, I/ is the sum of 84 monomials, degu .H40 / D 7, degvs .H40 / D 12
and, of course, H40 is linear in I. Then,
Resvs .H40 .u, vs .u/, I/, g.u, vs .u// D
constant u12 .u 1/16 .81000u C 19902511/.784u 17846163/h.u/R4 .u, I/,
where h.u/ is a second degree polynomial. In the right-hand side, the roots of the
non-constant factors, except for R4 .u, I/ belong to the interval .696 , 22763/. Also,
R4 .u, I/ is the sum of 12 monomials, degu .R4 / D 7 and degI .R4 / D 5. If we assume
that juj 22763, then R4 .u, I/ D 0 and we write the last relation as follows:
I 5 C q2 .u/ I 3 C q3 .u/ I 2 C q4 .u/ I C q5 .u/ D 0,
(8.19)
where
q.u/ D26244u4 C 20515275u3 C 1099852416u2 2071134904704u
374185039449856.
Now we work as follows. Consider q2 .u/. Its numerator has no real roots and those
of the denominator belong to the interval .484 , 318/. But we have already assumed
that juj 22763 and we compute that q2 .22763/ < 0 and q2 .22763/ > 0. Hence, if
u 22763, then q2 .u/ < 0 and if u 22763, then q2 .u/ > 0. In an analogous way
we check that, if juj 22763, then q3 .u/ > 0, and if u 22763, then both q4 .u/
and q5 .u/ are negative, and if u 22763, then both q4 .u/ and q5 .u/ are positive.
In view of the above discussion, if u 22763, then qi .u/ > 0 for i D 2, 3, 4, 5 and,
consequently, in (8.19), I < 0. Setting I D J < 0 we obtain the equation
J 5 C q2 .u/ J 3 q3 .u/ J 2 C q4 .u/ J q5 .u/ D 0,
(8.20)
where now the strictly negative coefficients are q3 .u/ and q5 .u/. By Cauchys rule,
0 < J < max.2q3 .u//1=3 , .2q5 .u//1=5 .
We have
q3 .u/ D
32 3u4 C 1355u3 4
u < .2=27/3 u4
9
q.u/
and
q5 .u/ D
93
32
u4
u7 < .23 =313 / u7 .
243 q.u/
jq4 .u/j D
1175056
u4
juj7
243
jq.u/j
jq5 .u/j D
32u4
juj7 .
243 jq.u/j
From these relations we easily obtain the following upper bound for I:
0 < I < 0.11 juj7=5 .
Next, let I < 0. We put I D J with J > 0, so that the equation (8.20) holds. Since
u 22763, the strictly negative coefficients in (8.20) are now qi .u/ with i D 2, 3, 4.
By Cauchys rule,
0 < J < max.3jq2 .u/j/1=2 , .3jq3 .u/j/1=3 , .3jq4 .u/j/1=4 .
Working as we did a few lines above, we obtain 0 < I D J < 0.11juj4=3 .
Combining all our partial results above we come to the following conclusion:
If juj 22763, then the left-hand side of (8.16) is bounded above by 0.11 juj4=3 ,
hence, in the notation of Proposition 8.4.2, B2 D 22763, D 1=3 D 1=s and
c9 D 0.11.
Remark. An easy search for integer solutions of g.u, v/ D 0 in the range juj < 22763
reveals no further solutions than those already found, namely, .u, v/ D .0, 0/, .1, 3/
and .243, 3/.
Example continued on page 96.
94
95
Section 8.7 The linear form L.P / and its upper bound
Z
D
C1
x.P CQ2 /
dx
p
f .x/
C1
x.Q2 CP0 /
dx
f .x/
1 C1 g.u/ du
jl.P / C "P l.P0 /j D
.
2 u.P / gv .u, vs .u//
(8.22)
Now we will prove a proposition which will permit us to compute an upper bound for
the right-hand side of (8.22) in terms of u.P /.
Proposition 8.7.1. Let us write the relation g.u, v/ D 0 in the form
v n C a1 .u/ v n1 C C an1 .u/ v C an .u/ D 0
(8.23)
where the ai s are polynomials in u. Let B3 be a constant larger than every root of
every non-zero polynomial ai .
If P C 2 C.Z/ and u.P / > maxB2 , B3 , then
h.x.P // c10 C c11 log ju.P /j,
(8.24)
where, h./ denotes absolute logarithmic height and c10 , c11 are explicitly computable
positive constants, independent from P .
Proof. We have v.P / D vs .u.P //. We write X D F1 =F2 for some relatively prime
polynomials with rational integer coefficients. For simplicity in the notation of this
proof let us put u.P / D u and v.P / D v, so that
h.x.P // D h.X .u.P /, v.P // log maxjF1 .u, v/j, jF2 .u, v/j D log jFj .u, v/j
for the proper choice of j D 1, 2.
Let 1 i n be such that ai is a non-zero polynomial. Then ai .u/ does not
change sign as u runs through the values > B2 . Therefore, for u > B2 , it makes sense
to distinguish all the (strictly) negative coefficients, say am1 .u/, am2 .u/, : : : , amk .u/,
in the left-hand side of (8.23). If v 0, Cauchys rule8 implies that
0 v max .kjami j/1=mi
1ik
and, clearly, the right-hand side is bounded by 0 u , where 0 and 0 are explicit
positive constants ( 0 2 Q). If v < 0, we put v D w with w > 0 and we rewrite
(8.23) as w n C b1 .u/ w n1 C C bn1 .u/ w C bn .u/ D 0, where bi D ai for
00
i D 1, : : : , n. As before, Cauchys rule implies a relation of the form 0 < w 00 u
and, combining the two upper bounds we obtain an explicit bound jvj u .
8
See Footnote 5.
96
P
Next, we write Fj .u, v/ D .k,l/ ak,l uk v l and let c11 D max.k,l/ k C l, where
the maximum P
runs over all pairs .k, l/ for which ak,l 0. Then, jFj .u, v/j c10 uc11 ,
where C10 D .k,l/ jak,l j l and (8.24) holds with c10 D log C10 .
ExampleStep 6 . (continued from the end of Step 5, page 93)
We have a1 D 0, a2 D u, a3 D 0, a4 D 271
u and a5 D u2 , hence we can take
3
B3 D 1.
We put .u.P /, v.P // D .u, v/ and, according to Proposition 8.7.1, we assume that
juj maxB2 , B3 D B2 D 22730.
Assume u > 22730.
If v 0, then Cauchys
rule
implies 0 v max.2ja4 j/1=4 , .2ja5 j/1=5 D
p
542 1=4 p
5
5
max. 3 u/ , 2 u2=5 D 2 u2=5 .
and, clearly, 30000juj3 is an upper bound for the absolute value of the denominator
u.u 1/2 of X .u, v/. Thus,
h.x.P // D log jx.P /j
D log [Link] .u.P /, v.P //j/ , [Link] .u.P /, v.P //j/
log.30000juj3 / D log.30000/ C 3 log ju.P /j
and, consequently, c10 D log.30000/ and c11 D 3.
Example continued on page 118.
Theorem 8.7.2. Let E be the elliptic curve represented by the model C in (8.1), subject to the conditions at the beginning of this chapter.
def
Let E : y 2 D x 3 C Ax C B D f .x/ be the short Weierstrass model with A, B
defined in Proposition 8.1.1 and consider that everything in Chapter 4, in particular,
equation (4.1) and relations (4.4)(4.8), refer to this particular model E.
Let X be the rational function defined in Section 8.1 and let P0 be the point of E in
Definition 8.3.3.
Section 8.7 The linear form L.P / and its upper bound
97
(8.26)
For B2 see Proposition 8.4.2 and Section 8.6; for B3 see Proposition 8.7.1.
also Section 8.6.
10 See
Chapter 9
<0
" D 1
1
:
1
where
and
m00
99
(9.2)
where the meaning of the various quantities involved in it will be made explicit by
distinguishing three cases.
Let " D 0. Remember the definition of r0 at the beginning of Chapter 4, a few lines
after (4.1), and the relation (4.5), according to which, t jr0 . Therefore, we can write
m00 C
s
n0
D ,
t
r0
n0 2 Z.
On putting
k D r,
d D 1,
ni D mi .i D 1, : : : , k/,
where
(9.3)
Therefore, by (9.3),
def
N D maxjni j, i D 0, 1, : : : , k N0 .
Let " D 1 and the points P1E , : : : , PrE , P0E are Z-linearly independent. Clearly, a
necessary condition for this case is that P0E 62 E.Q/ but, certainly, this condition is
by no means sufficient. In this case, m00 D m0 and, in analogy with the case " D 0,
we put
n0
s
m 0 C D , n0 2 Z
t
r0
100
and
k D r C 1,
d D 1,
ni D mi .i D 1, : : : , k 1/, `k D `0 , nk D ",
N D maxjni j, i D 0, 1, : : : , k N0 ,
where N0 is again given by the right-most side of (9.3) with D 1.
Let " D 1 and the points P1E , : : : , PrE , P0E are Z-linearly dependent. In this case
there exist integers , d1 , : : : , dr with 1 and gcd., d1 , : : : , dr / D 1, such that
P0E D d1 P1E C C dr PrE C T0 .
Since the points PiE (i D 0, 1, : : : , r) are known, we will consider the di s as integers that are explicitly known. Applying the homomorphism l to the last displayed
relation, we get `0 d1 `1 C dr `r C st 0 !1 .mod Z!1 /, where2 l.T0 / D s0 =t0 .
0
From this it follows that
s0
`0 D d1 `1 C dr `r C !1 C 0 !1 ,
t0
where 0 is an appropriate explicitly known integer. We note that the integers s0 , t0
satisfy conditions similar to those in (4.5) with s0 and t0 in place of s and t , respectively; in particular, t0 jr0 and js0 =t0 j 1=2.
Now our linear form becomes (note that m00 D m0 in this case)
L.P / D .m0 C
0
s0
s
d1
dr
C
C
/!1 C .m1 C /`1 C C .mr C /`r . (9.4)
t
t0
s
d0
s0
n0
C
C
D
t
t0
dr0
ni D d mi C di .i D 1, : : : , k/,
so that L.P / is once again expressed by (9.2), and we compute now an upper bound
for max0ik jni j. Obviously, max1ik jni j dM C max1ik jdi j. For jn0 j we
take into account that js=t j, js0 =t0 j 1=2 (cf. (4.5)) and we compute
s
s0
jn0 j D jd m0 r0 C dr0 C r0 d0 C r0 j dr0 jm0 j C 12 dr0 C r0 d0 C 12 r0
t
t0
dr0 maxM , 12 rM C 1 C 12 .d C 1/r0 C jd0 jr0 ,
2
3
101
N D maxjni j, i D 0, 1, : : : , k N0 ,
where
N0 D maxdM C max jdi j, dr0 maxM , 12 rM C1C 12 .d C1/r0 Cjd0 jr0 . (9.5)
1ik
(9.6)
Remark. A remark concerning the points to which the `i s refer is necessary. According to Conclusions and remarks (2) on page 52, if `i D [Link] / and it happens that
Pi 2 E1 .Q/,4 , then `i is the elliptic logarithm of the point PiE C Q2E , where, as
usually, Q2E D .e2 , 0/. Therefore, if e2 62 Q, then e2 belongs to either a quadratic or
a cubic number field, hence this will be the case with the coordinates of PiE C Q2E ,
E
as well; anyway,
p0 /, the coordinates of P0
p D 3. In the case that `i D `rC1pD l.P
belong to Q. a/ if we apply Theorem 6.8 or to Q. D1 , D2 / if we apply Theorem
7.1. In the case of Theorem 6.8, if P0E 2 E1 .R/, then l.P0 / is the elliptic logarithm
of the point P0E C Q2E , the coordinates of which belong to a number field of degree
D 6 (D 5). In the case of Theorem 7.1, P0E 2 E0 .R/ and Q2E 2 E.Q/, as
E
E
already
p noted
pon page 75. Therefore, the coordinates of P0 C Q2 belong, anyway,
to Q. D1 , D2 /, so that D is a divisor of 4 in this case. If we apply Theorem 8.7.2
we cannot be so specific concerning the number field generated by the coordinates of
P0E C Q2E in the case that P0E 2 E1 .R/.
Our next task is to compute an explicit upper bound for M by using a very important
result of S. David, namely, Theorem 9.1.2. The following strategy will be followed.
Theorem 9.1.2 will provide us with a lower bound for L.P / as a function of N
hence, as a function of M as well which we will combine with the upper bounds
of jL.P /j obtained in Theorems 5.2, 6.8, 7.1 and 8.7.2. The upper bounds given in
those theorems are again functions of M . However, as it will turn out, if M is very
4
102
large, then the upper bound function is smaller than the lower bound function and
this contradiction will lead us to an upper bound for M .
We proceed now to our task. First, let us note that, for the application of Theorem 9.1.2, we will need to compute a fundamental pair of periods .$1 , $2 / as in
Lemma 9.1.1 below. By Theorem 3.4.1 and Section 3.4.4 we can explicitly calculate a fundamental pair of periods .!10 , !20 / with !10 D !1 . From this pair it is possible
to compute the pair .$1 , $2 /.
Lemma 9.1.1. For the Weierstrass function } associated to (4.1) there exists a fundamental pair .$1 , $2 / of periods with the property that, if we set Q D $2 =$1 , then
the following inequalities are satisfied:
= Q > 0,
j j
Q 1,
1
j< Q j .
2
M D
a
c
b
d
/ 2 GL2 .Z/
a C b
.
c C d
Now we proceed to the proof. Let .!1 , !20 / be a fundamental pair of periods for }
obtained by Theorem 3.4.1, where, in the notation of that theorem, !20 D !2 in the case
of positive discriminant, and !20 D .!1 C !2 /=2 in the opposite case. Let D !20 =!1
(note that = > 0). By [45, Proposition 12.1], there exists
a0 b0
2 SL2 .Z/
M0 D
c0 d0
such that M0 satisfies both j<.M0 /j 12 and jM0 j 1. Actually, M0 is a
word in SL2 .Z/ generated by
0 1
1 1
SD
and T D
.
1 0
0 1
Put $1 D c0 !20 C d0 !1 and $2 D a0 !20 C b0 !1 . We have det M0 D 1, therefore
.$1 , $2 / is a fundamental pair of periods. Moreover,
def
Q D
$2
a0 C b0
D
D M0 ,
$1
c0 C d0
103
Set
D D K : Q,
where K is the least number field containing the coordinates of all points of E, the
elliptic logarithms of which are `1 , : : : , `k , respectively; cf. the remark on page 101.
Set
hE D max1, h.1 : 4A : 4B/, [Link] /,
where jE D 28 33 A3 =.4A3 C 27B 2 / is the j -invariant of E and h. / means absolute
logarithmic height as defined in Section 2.3, especially page 17 ff.
Compute first a fundamental pair of periods .!1 , !20 / and then a fundamental pair of
periods .$1 , $2 /, as explained in Lemma 9.1.1.
O i / using a software package, e.g. PARI or MAGMA, if
For i D 1, : : : , k compute h.P
O i / using Proposition 2.6.4.
this is possible; if not, compute an upper bound for h.P
Note that the routines of the above mentioned packages for the canonical height5
can do the job always6 if the point has rational coordinates, but refuse7 to compute
the canonical height of a point with irrational coordinates. In such a case, we confine
O 0 / using Proposition 2.6.4.
ourselves to computing an upper bound for h.P
3j!1 j2
H0 max hE ,
,
Dj$1 j2 = Q
3`2i
O
, [Link] /
Hi max hE ,
Dj$1 j2 = Q
7
8
.i D 1, : : : , k/.
104
Theorem 9.1.2 (Sinnou David). All points P E 2 E.Q/ for which L.P / D 0 can
be explicitly calculated. More specifically, such points can result only in either of the
following two cases: (1) If P is a torsion point, or (2) If the situation described in the
third bullet, page 100, occurs, and .d , m1 , : : : , mr / D .1, d1 , : : : , dr /.
If L.P / 0, then, either N c12 , where
c12 D maxexp.ehE /, jdr0 j, [Link] =D/, i D 0, : : : , k,
or
jL.P /j > exp c13 .log N C c14 /.log log N C c15 /kC2 ,
where
c13
c14 D 1 C log D ,
c15 D 1 C hE C log D.
2 C13kC23.3
k
Y
(9.7)
Hi ,
iD0
Proof. In the case of the first bullet, page 99, the relation L.P / D 0 means that
r0 m1 l.P1 / C C r0 mr [Link] / C n0 !1 D 0. By Theorem 3.5.2 (c) it follows then that
m1 D D mr D 0, hence P is a torsion point.
In the case of the second bullet, page 99, the relation L.P / D 0 is impossible.
Indeed, in this case, L.P / D 0 means that r0 m1 l.P1 / C C r0 mr [Link] / r0 l.P0 / C
n0 !1 D 0 and the points P1E , : : : , PrE , P0E are Z-linearly independent. Since the
coefficient of l.P0 / is non-zero, this contradicts Theorem 3.5.2 (c).
In the case of the third bullet, page 100, L.P / D 0 means (cf. (9.4)) that r0 .d m1 C
d1 /`1 C C r0 .d mr C dr /`r C n0 !1 D 0 and Theorem 3.5.2 (c) implies that mi D
di =d for i D 1, : : : , r. In particular, d is a common divisor of d1 , : : : , dr and, since
gcd.d , d1 , : : : , dr / D 1, we conclude that d D 1 and mi D di for i D 1, : : : , r.
If L.P / 0, then the conclusion of the theorem results from an almost straightforward, though careful (see remarks below), application of [12, Thorme 2.1].
Remarks on the application of S. Davids theorem.
(1) In S. Davids theorem [12, Thorme 2.1] the k elliptic logarithms appearing in
the linear form may come from different elliptic curves defined over a number field
K. In our case, all elliptic curves coincide with the elliptic curve (4.1) defined over
Q. We remind, however, that some `i s may be elliptic logarithms of points with
non-rational coordinates, as noted on page 101.
(2) For the parameter E appearing in [12, Thorme 2.1] we chose E D e.
(3) A detail concerning canonical heights. According to our discussion in Chapter 2,
the canonical height in [12, Thorme 2.1], which on page 20 ff. we denoted by hOD ,
is, by Proposition 2.6.1, three times the canonical height defined by J. Silverman
in [45, section VIII.9], i.e. the canonical height that we adopt in this book. This
O i / is multiplied by a
explains why in the definition of Hi the canonical height h.P
factor 3, which does not appear in S. Davids theorem.
105
Now we are ready state the theorem from which an explicit upper bound of M is
obtained.
Theorem 9.1.3. Let L.P / be the linear form appearing in any of Theorems 5.2, 6.8,
7.1 or 8.7.2. Write L.P / in the form (9.2) following the instructions in the appropriate
bullet (page 99 ff.) and consider the constants and (cf. (9.6)), as well as the
constants c12 , c13 , c14 , c15 resulting from Theorem 9.1.2.
If L.P / 0, then, either M c12 , or
M 2 c18 c13 .log.M C/Cc14 /.log log.M C/Cc15 /kC2 C Cc18 log c16 Cc17 ,
(9.8)
where,
8
..2V0 D1 D2 / , 2 c8 , 1/
:
1
.c9 =.1 C / , 2 c10 , .c11 =.2 // in case of Theorem 8.7.2.
Proof. Let the linear form L.P / considered in any of the Theorems 5.2, 6.8, 7.1 or
8.7.2 and define N according to the appropriate bullet, page 99 ff., as the case may
be.
If M > c12 , the left inequality (9.6) implies N > c12 and then Theorem 9.1.2
implies the lower bound (9.7) for jL.P /j. In view of the right inequality (9.6), this
lower bound is also valid after replacing N by M C . Then the resulting lower
bound, combined in a straightforward manner with the upper bound for jL.P /j already
obtained from Theorem 5.2 or 6.8 or 7.1 or 8.7.2, depending on the case, leads to the
inequality (9.8).
Concluding remark. In (9.8), the left-hand side tends to infinity faster than the righthand side. Therefore, if M is sufficiently large, then this relation cannot hold. Consequently, if all constants in (9.8) are explicitly known, then an upper bound for M can
be explicitly calculated; cf. Chapter 4, immediately after (4.8).
9.2
Computational remarks
In sections 9.3, 9.4, 9.5 and 9.6 we show by concrete examples how to apply Theorem 9.1.3 and the Concluding remark following it. In the present section we give some
practical information concerning our concrete numerical computations.
We make a combined use of MAGMA, PARI and a few rather simple codes based on
MAPLE written by the author. In all examples we started doing all calculations with a
precision of 100 decimal digits. With the exception of Section 9.4, in all other sections,
the numbers resulting from the calculations are of a size much less than 10100 , so that
106
we consider them reliable. In Section 9.4, however, the value of the constant c13 turned
out to be considerably larger that 10100 and we repeated all the computations using a
200 decimal digits precision. All packages nowadays work with this and much larger
precision without the slightest problem. Anyway, in the text we write the results of
our computations correct to 16 decimal digits.
Computations related to elliptic curves are done using MAGMA and/or PARI. Given
the vector a1 , a2 , a3 , a4 , a6 of the coefficients of an elliptic curve (see [45, Chapter III.1]), the routines EllipticCurve of MAGMA and ellinit of PARI create an object
e containing all necessary information to work with. Minimal models are calculated
by any of the routines ellminimalmodel of PARI or MinimalModel of MAGMA.
For the computation, even with a very high precision, of a fundamental pair of
periods for our elliptic curves (or, rather, for the }-functions associated to them) we
can easily write a simple code according to the instructions in Section 3.4.4. If one
does not feel comfortably with programming, any of the commands Periods.e/ of
MAGMA or [Link] of PARI does the job.
Elliptic logarithms can be calculated by the algorithm on page 53, for which it
is not difficult to write a code. Alternatively, the routines ellpointtoz of PARI or
EllipticLogarithm of MAGMA do the same thing. We remind the reader that, if for
some point P we have P E 2 E1 .Q/, then l.P / is the elliptic logarithm of P C Q2
(where Q2E D .e2 , 0/) rather than the elliptic logarithm of P , which is not a real number; cf. Conclusions and remarks (1) on page 51. Probably, the point P E C Q2E has
non-rational coordinates; in this case, its coordinates belong to a quadratic or a cubic
number field. But computing the elliptic logarithm of such a point is not really a problem, because the algorithm on page 53 works perfectly with real numbers. Also, the
routine ellpointtoz of PARI accepts input points on the elliptic curve with coordinates
of type real, but the routine EllipticLogarithm of MAGMA would complain. In a
symbolic calculation package like, for example, MAPLE, we can write a program for
calculating with points on an elliptic curve over a number field, by implementing the
Group law algorithm 2.3 of [45, Chapter III.2].9
The computation of the canonical height is quite a sophisticated task; see, for example, [42, 43, 75, 76].10 For the needs of our examples we used either of the routines
CanonicalHeight of MAGMA or ellheight of PARI.
Warning! The definition of the canonical height given in Section 2.6 agrees with
the one given by Silverman [45, Chapter VIII]. The canonical height computed by the
above routines of PARI and MAGMA is twice as large, because the factor 1=2 in (2.36) is
missing from the definition of the canonical height adopted by PARI and MAGMA. Analogously, for the height pairing matrix, either of the routines HeightPairingMatrix or
ellheightmatrix of MAGMA or PARI, respectively, is used. Note that the eigenvalues
9
10 The
107
of the height pairing matrix are not affected by the normalisation for the canonical
height that one may adopt.
Warning! The routines ellheight and ellheightmatrix of PARI demand that their
input data come from a minimal Weierstrass model.
The extra point P0 that makes its appearance possibly in Theorem 6.8 and, anyway, in
Theorems 7.1 and 8.7.2, has probably irrational coordinates and, in this case, neither
PARI nor MAGMA know how to compute its canonical height. Should we need to know
O 0 /, we have to implement for ourselves the algorithms of [43].
the precise value of h.P
But for the application of Theorem 9.1.2 we only need a reasonably good upper bound
O 0 /, which we obtain from Proposition 2.6.4.
of h.P
Computing the constants c12 , : : : , c18 is a straightforward but quite boring task.
Therefore, we preferred to write a rather simple MAPLE code for their computation.
A considerably more sophisticated code running, for example, under MAGMA, which
would receive much less data only, say, the coefficients of C and E and return the
numerical values of above ci s automatically, is possible. We did not attempt to write
such a code as we prefer a more transparent process which is more easily under our
control.
9.3
Taking in (5.2)
.
, ,
, / D .1, 5=12, 1=2, 7=24/,
we obtain the model
143623
3409
x
D f .x/.
48
864
The roots of f .X / are e3 D 41=6 < e2 D 31=12 < e1 D 113=12. Note also that,
in Proposition 5.1,
D 12.
E : y2 D x3
We will work with the elliptic curve E, two models of which are C and E, searching
for all points P on E such that P C 2 C.Z/. In order to apply Theorem 5.2, we need
that u.P / > 18. Also, we note that there are no points .u, v/ 2 C.R/ with u < 7.11,
therefore we first compute all P with integer u.P / between 7 and 18; these are
exactly the following:
.u.P /, v.P // D.7, 1/, .7, 7/, .6, 5/, .6, 10/, .4, 5/, .4, 8/, .3, 1/
.9, 5/, .13, 43/, .13, 29/, .14, 50/, .14, 35/.
From now it will be understood that P
(9.9)
108
We compute E.Q/ Z2 Z2 ZZ, so that the rank r D 2. The torsion subgroup
consists of the points O and QiE D .ei , 0/ (i D 1, 2, 3). Thus,
Etors .Q/ D hQ1E , Q2E i D h.113=12, 0/, .31=12, 0/i
&
r0 D 2.
Let now g.u, v/ D 0, where u, v 2 Z and u > 18, and consider the point P C D
.u, v/ 2 C.Z/. Then P E 2 E.Q/ and, specialising (4.2) to our case, we have
P E D m1 P1E C m2 P2E C T E ,
implying
s
L.P / D l.P / D .m0 C / !1 C m1 l.P1 / C m2 l.P2 /,
2
which is the specialisation of (4.6) to our case. By Theorem 9.1.2, L.P / D 0 implies
m1 D m2 D 0, i.e. P E 2 Etors .Q/ D Q1E , Q2E , Q1E C Q2E , O. Only for P D
Q1 , Q2 we obtain a point P C with integer coordinates, namely Q1C D .9, 5/ and
Q2C D .3, 1/, already listed in (9.9). Therefore, in what follows we assume that
L.P / 0. The linear form L.P / falls under the scope of the first bullet, page 99,
so that
k D 2,
d D 1,
n0 D 2m0 C s, jsj 1,
.n1 , n2 / D .m1 , m2 /,
!2 1.0599189920498769 i ,
noting that
D !2 =!1 1.2625682091829191 i
satisfies = > 0, j j > 1 and j< j < 12 ; hence, in Lemma 9.1.1 we can take
.$1 , $2 , Q / D .!1 , !2 , /.
We calculate the canonical heights
O 1 / 1.0547377027055837,
h.P
O 2 / 0.3611721523399011
h.P
109
`2 D l.P2 / 0.2755755724720115.
Next we compute
c12 6.4272819599100395 1028 ,
c14 D 1,
c15 25.402522816413709.
1
2
Assuming M > c12 , we insert into (9.8) the values of , , , and c13 , : : : , c18 that
we computed, and conclude that B.M / > 0, where
B.M / D 9.11632 1073 .log.2M C 3/ C 1/[Link].2M C 3// C 25.4025234
C 7.32132 0.61064M 2 .
But, for all M > 1.1 1041 , we check that B.M / < 0, which implies that
M maxc12 , 1.1 1041 D 1.1 1041 .
(9.10)
We cannot obtain an upper bound for M essentially better than the above using Theorem 9.1.3; indeed, we check that B.1.07 1041 / > 0 which shows that with a little
smaller bound for M we do not arrive at a contradiction. We will succeed to make
an impressive reduction of the huge upper bound (9.10) in Section 10.2.1 of the next
chapter.
110
31
,
3
5 4 1 3
u C u C u2 C 1 v 2 .
2
2
BD
685
,
108
so that
E : y 2 D f .x/ D x 3 C Ax C B.
The approximate values of the roots of f .X / are
e1 3.4860779869512022 > e2 0.6390559614483767
> e3 2.8470220255028254.
We will work with the elliptic curve E, two models of which are C and E. The rank
of E is 3 and the torsion subgroup Etors .Q/ is trivial; thus
r D 3, r0 D 1.
The following points form a basis for E.Q/:
P1E D .2=3, 1=2/,
Let now .u, v/ be an integer solution of the equation Q.u/ D v 2 and P C D .u, v/ 2
C.Z/. Then P E 2 E.Q/ and, specialising (4.2) to our case, we have
P E D m1 P1E C m2 P2E C m3 P3E
implying
l.P / D m0 !1 C m1 l.P1 / C m2 l.P2 / C m3 l.P3 /,
which is the specialisation of (4.6) to our case. We also compute a fundamental pair
of periods
!1 1.3856089566101190,
!2 1.5968599642739782 i ,
with
D !2 =!1 1.1524607694372033 i
satisfying = > 0, j j > 1 and j< j < 12 , hence, in Lemma 9.1.1 we can take
Q D .!1 , !2 , /.
.$1 , $2 , /
The approximate values of the canonical heights are
O 1 / 0.55387253844458413, h.P
O 2 / 0.6064693463824573,
h.P
O 3 / 0.6518823598427176
h.P
111
112
Table 9.1. Parameters and auxiliary functions for the solution of the quartic elliptic equation
example (cf. Table 6.1).
Q.u/ D 5 u4 C 1 u3 C u2 C 1
Q.u/ D 5 u4 1 u3 C u2 C 1
D1
D 1
u2 C 6 C 6 .Q.u//1=2
3u2
u0 D 1
u D 1
u D 200, c7 D 3.498
p
P0E D . 1 C 10, 1 /
x.u/ D
x.u/ D
u
E
P0
u2 C 6 C 6 .Q.u//1=2
3u2
u0 D 17
u D 17
D 17, c 7 D 3.496
p
D 1 C 10, 1
3
l.P 0 / D l.P0 /
Using the routine IsLinearlyIndependent of MAGMA, we see that the points P0E ,
P1E , P2E , P3E are Z-linearly independent,11 so that we are in the situation described
in the second bullet, page 99. Therefore, in (9.2),
k D 4, d D 1, r0 D 1, .n0 , n1 , n2 , n3 / D .m0 , m1 , m2 , m3 /, n4 D 1, `4 D `0 .
In (9.3) N0 D 32 M C
3
2
hence, in (9.6),
., / D .3=2, 3=2/.
We compute
`0 D l.P0 / 0.6737142222951107.
O 0 /. Since P E is not a
On applying Theorem 9.1.2, we will also need to compute h.P
0
rational point we confine ourselves to a reasonably good upper bound of its canonical
height (see Section 9.2) which we obtain from Proposition 2.6.4. In the notation of
that proposition we take D : y 2 C y D x 3 C x 2 p
10x 10, which is, actually,
D
thepminimal model of our elliptic curve. Then, P0 D . 10, 0/. By Proposition 2.4.2,
h. 10/ D .log 10/=2 and then, by Proposition 2.6.4,
O 0 / 4.12.
h.P
Finally, following the section Preparatory to Theorem 9.1.2 we see that
D D 6.
11 For
a more transparent way to prove this fact see the end of this section.
113
Indeed, in L.P / we have the elliptic logarithms of the points PiE C Q2E , the coordinates of which belong
to Q.e2 /, and the elliptic logarithm of the point P0E , with
p
coordinates in Q. 10/. Therefore,
the least number field containing the coordinates
p
of all these points is Q.e2 , 10/, the degree of which is 6.
We also compute
c12 9.5771564156066246 1021 ,
c14 2.7917594692280550,
c17 1.1754018326322600,
c18 D 1.
Assuming M > c12 , we insert into (9.8) the values of , , , and c13 , : : : , c18 that
we computed, and conclude that B.M / > 0, where
B.M / D 5.607 10165 log.1.5M C 1.5/ C 2.79176
[Link].1.5M C 1.5/ C 21.41156 C 5.96764 0.8549 M 2 .
But, for all M 2.3 1088 , we check that B.M / < 0, which implies that
jM j maxc12 , 2.3 1088 D 2.3 1088 .
(9.11)
We cannot obtain an upper bound for M essentially better than the above using Theorem 9.1.3; indeed, we check that B.2.2 1088 / > 0 which shows that a little smaller
bound for M does not lead to a contradiction. A very big reduction of the huge upper
bound (9.11) will be accomplished in Section 10.2.2 of the next chapter.
Z-linear independence of P0E , P1E , P2E , P3E . If these points were Z-linearly deE 2 E.Q/. We show that
pendent, then a positive integer d would exist, such that dP0p
E D .1=3 C , 1=2/, where D
10, our claim is obvious
this is impossible. Since P0
for d D 1. Now we proceed by induction. Actually, we show something stronger:
If RE D .m C n, p C q/ 2 E.Q.//, m, n, p, q 2 Q, then RE C P0E 62 E.Q/.
We prove this as follows. From f .m C n, p C q/ D 0 we obtain
108m3 C 3240mn2 1116m 685 108p 2 1080q 2 D 0
1080n3 C 324m2 n 1116n 216pq D 0.
(9.12)
114
the contrary and consider the two equations (9.12) along with the two equations
p11 .m, n, p, q/ D 0 and p21 .m, n, p, q/ D 0. Elimination of q from these four equations leads to three polynomial equations in m, n, p with integer coefficients. A new
elimination of p from the last three equations leads to two polynomial equations, say
gi .m, n/ D 0, i D 1, 2 with integer coefficients. Specifically,
g1 .m, n/ D n8 .9m2 C 30n2 31/6 h1 ,
10t 2 60t 3
,
3.10t 2 C 3/
nD
10t 2 C 2t 3
.
10t 2 C 3
(9.13)
A1 D 14,
D2 D 7, A2 D 6.
p
p
We check for which V D 1, : : : , 1000 both A1 C D1 V 2 and A2 C D2 V 2 are
integers; only V D 1, 5 have this property, giving us the solutions .U , V , W / D
.5, 1, 1/, .17, 5, 13/. As a particular solution .U0 , V0 , W0 / we choose
.U0 , V0 , W0 / D .5, 1, 1/
115
(9.14)
123412480
326848
xC
.
3
27
As noted a few lines after (7.28), f .X / has three rational roots, which are
1048
128
920
< e2 D
< e1 D
.
3
3
3
The functions u and v, defined in (7.18) and (7.19), are
p
5V 14 C 11V 2
,
p
u.V / D
V C 6 C 7V 2
p
p
4.7 6 C 7V 2 C 82V 15 14 C 11V 2 /
.
v.V / D
p
.V C 6 C 7V 2 /2
e3 D
Plotting v.V / and Q.u.V // for V > 1000 we see that v.V / and Q.u.V // are positive,
as the second bullet in the announcement of Theorem 7.1 requires. Also, according
to (7.21) and (7.22),
p
376
28 Q.u/ 280u C 392
x.u/ D
C
3
u2
p
1680u3 10528u2 11760u C 10976 C .784 280u/ Q.u/
y.u/ D
.
u3
The functions V 7! x.V / and V 7! y.V / are defined in (7.23) by
x.V / D .x u/.V /,
y.V / D .y u/.V /.
An explicit expression of x.V / and y.V / is an easy matter for a symbolic computation package, like MAPLE, for example, but this is a too long expression to print here.
Besides, what we only need to check is (see third bullet in the announcement of Theorem 7.1) that, in the interval .1000, C1/ the function V 7! y.V / is positive and
the function V 7! x.V / is strictly monotonous, with values exceeding e1 . Plotting
116
these two functions for V > 1000 we see that this is indeed true. More precisely, the
function V 7! x.V / is strictly decreasing and its limiting value is
p
p p
p
920
4164 2300 11 C 924 7 420 11 7
p
>
D e1 .
3
54 C 15 11
Finally, by (7.24) and (7.25) we compute P0E D .x0 , y0 /, where
p
p
p p
1052
C 44 7 C 140 11 C 20 11 7
x0 D
3
p p
p
p
p
p
y0 D 8 11 7 .5 7 C 11/.1 C 7/.5 C 11/.
Now we calculate the constant c8 following the guidelines of the paragraph before the
relation (7.32). From (7.31) we have
x.V / D P1 .U , V , W /=P2 .U , V , W /
D 2.840U C 4592V C 392W 479V 2 C 130U V
C U 2 308V W C 140U W C 196W 2 /=.U 5V /2 .
We obtain an upper bound for P1 .U , V , W / of the form constant times V 2 , where the
constant can be obtained using (9.14) and V 0.001V 2 . We obtain a much better
bound of this shape if we consider the function
p
p
jP1 . 14 C 11V 2 , V , 6 C 7V 2 /j
V2
and observe that this is strictly decreasing in the interval .1000, C1/ taking a value
< 3503.22
at V Dp1000. Working similarly for the denominator we see that
p
jP2 . 14 C 11V 2 , V , 6 C 7V 2 /j=V 2 is a strictly increasing function with limiting value < 2.834 as V ! C1. Hence, we can take
c8 D log.3503.22/.
Now it is time to look for the arithmetic-geometric details of the model E. Our
calculations have been made with a precision of 100 decimal digits; below we exhibit
the numerical values in 16 digits.
The free part of the group E.Q/ is generated by
P1E D .1976=3, 14784/,
Also, Etors .Q/ Z2 Z2 is generated by Q1E D .e1 , 0/ and Q2E D .e2 , 0/, so that
Etors .Q/ D h.920=3, 0/, .128, 0/i
&
r D 2, r0 D 2.
!2 D 0.1389333214997840 i
117
$2 D !1 .
For the point P E D .x.V /, y.V // coming from the unknown integer solution
.U , V , W / of (7.1), with V > 1000 and U , W 0, we write
P E D m1 P1E C m2 P2E C T E ,
T E 2 Etors .Q/,
n0 D 2m0 C s, s 2 1, 0, 1,
`2 D l.P2 / 0.0606534846504600.
O 2 / 1.0456910246342577,
h.P
118
c15 30.7739853633026035.
c17 4.0807189122684443,
c18 D 1.
Assuming M > c12 , we insert into (9.8) the values of , , , and c13 , : : : , c18 that
we computed, and conclude that B.M / > 0, where
B.M / D 1.6559 1074 .log.2M C 3/ C 1/.[Link].2M C 3// C 30.774/4
C 6.702435 0.9301430 M 2
But, for all M 1.7 1041 , we check that B.M / < 0, which implies that
jM j maxc12 , 1.7 1041 D 1.7 1041 .
(9.15)
Since B.1.6 1041 / > 0, the bound (9.15) is essentially the best upper bound that can
be obtained from Theorem 9.1.3; its reduction to an upper bound of manageable size
will be accomplished in Section 10.2.3 of the next chapter.
c17 D
1
2
c18 D 92 .
119
$2 D !1 C !2 ,
O 2 / 1.9162842498242523,
h.P
`2 D l.P2 / 0.0038280584521748.
120
c14 D 1,
c15 D 26.51416.
Now we have all data to apply Theorem 9.1.3. Inserting into (9.8) the values of the
constants , , , and c13 , : : : , c18 which we have already computed, we see that
B.M / > 0, where
B.M / D 4.689 1074 .log.M C 1/ C 1/ .[Link].M C 1// C 26.5142/4
0.6127 0.46459 M 2 .
But, for all M 3.1 1041 , we check that B.M / < 0, which implies that
jM j maxc12 , 3.1 1041 D 3.1 1041 .
(9.16)
Since B.3 1041 / > 0, the bound (9.16) is essentially the best upper bound that can
be obtained from Theorem 9.1.3; its reduction to an upper bound of manageable size
will be accomplished in Section 10.2.4 of the next chapter.
Chapter 10
In this chapter we will show how one can reduce the upper bound for M , which was
obtained from the application of Theorem 9.1.3 and the Concluding remark following
it. In our examples, more specifically, in Sections 9.39.6 we obtained an upper bound
for M so large that there is no hope to check which points P cover solutions of the
Diophantine problem under consideration; cf. relevant discussion on page 55, below
(4.3). In the terminology therein, the Diophantine problems in our examples have been
solved effectively but not explicitly. In the present chapter we aim to show how one
can use the superficially non-practical huge upper bound for M in order to make an
enormous jump down to a new really small upper bound. Sounds strange but, still,
it is true! This new bound is so small that it allows one to check directly which points
P cover actual solutions and leads provably to the complete set of them.
We keep, of course, all notations etc. of Chapter 9; there we wrote L.P / for the
linear form (9.2). Let us consider the following slightly different linear form, which is
more appropriate for the purposes of the present chapter,
dr0
L.P / D n0 C n1 1 C C nk k ,
!1
r0 `i
.i D 1, : : : , k/
i D
!1
D .P / D
(10.1)
The upper bound for jL.P /j given by (5.3), (6.15), (7.30) and (8.25) can be written
1 . C c M 2 // (for c , c and c see Theouniformly as jL.P /j c16 exp.c18
17
16 17
18
rem 9.1.3). Moreover, by (9.6)),
N D max jni j M C ,
0ik
so that M 2
3 N 2 , for some explicit positive constant
3 . It follows then that
jj
1 exp.
2
4 N 2 /,
(10.2)
where
1 D
dr0 c16
,
!1
1
2 D c18
. C c17 /,
1
4 D c18
3 .
(10.3)
122
i1
X
ij bj ,
j D1
ij D
hbi , bj i
jbj j2
where hi denotes the usual inner product. The fact that .b0 , b1 , : : : , bk / is an LLLreduced basis means, by definition, that
j
ij j 1=2
1
.0 j < i k/
b0 , b1 , : : : , bk is a basis of RkC1 but, in general, one cannot expect that it is a basis for .
123
and
.0 < i k/.
Starting from any basis of a given lattice in our case, from the basis consisting of the
column vectors of M , the LLL-algorithm of [27] computes an LLL-reduced basis
very effectively; in our case, this means that the number of arithmetic operations that
are needed are O..k C 1/4 log C /, in view of [27, Proposition 1.26].
LLL-reduced bases, in general, have a number of important properties, extremely
useful for various applications, among them simultaneous approximation [27, Proposition 1.3] and factorisation of polynomials with rational coefficients [27, Sections 2
and 3]. The close connection of LLL-bases to the explicit solution of Diophantine
equations was described for the first time2 in [72]. Applications to the explicit resolution of various classes of Diophantine equations are found, for example, in3 [64] (Thue
equations), [65] (ThueMahler equations), [13] (norm form equations), [14] (relative Thue equations), [54, 15] (Weierstrass equations), [49] (Weierstrass equations
over number fields), [62] (quartic elliptic equations), [59] (general cubic equations),
[57] (elliptic binomial equations), [63] (simultaneous Pell equations), [56] (general elliptic equations), [48] (triangularly connected decomposable form equations), [47, 36]
(S -integral solutions of Weierstrass equations); see also the book [50] for a concise
overview of applications to Diophantine equations.
An important property of the LLL-bases that we will use in this section is the following: Once we know b0 , we can have an explicit lower bound for the length of any
non-zero vector of . More specifically,
jxj 2k=2 jb0 j
(10.4)
Suppose now that we want to solve the inequality (10.2), with as in (10.1), where
.n0 , n1 , : : : , nk / 2 ZkC1 and jni j < B1 .N / for i D 0, 1, : : : , k. Here we imagine
B1 .N / as a very large integer.
We consider the following x 2 :
1
1 0
0
1 0
n1
n1
n1
C
B . C B . C B
..
C
B . C B . C B
def
C,
C
B
B
C
B
.
x D M B . C D B . C D B
C
A
@ nk A @ nk A @
nk
n0
n1 C 1 C C nk C k C n0 C
Q
where Q D n1 C 1 C C nk C k C n0 C . For i D 1, : : : , k we have C i D
C i C i , for some i absolutely less than 1, hence, Q D C .n1 1 C C nk k /
and, consequently,
Q C jj C kN < C
1 e 2 4 N C kB1 .N /.
jj
2
2
3
124
jb0 j
k
X
2
jni j2 C Q 2 kB1 .N /2 C .C
1 e 2 4 N C kB1 .N //2 .
iD1
4 N 2
2 C log.
1 C / log 2k jb0 j2 kB1 .N / kB1 .N /.
(10.5)
In order that the left-hand side makes sense, the quantity inside the brackets must be
a positive real number, which is equivalent to 2k jb0 j2 kB1 .N /2 > k 2 B1 .N /2 .
This, in turn, is equivalent to 2k jb0 j2 > .kB1 .N / C 12 /2 14 and, clearly, a sufficient
condition for the last inequality is
jb0 j > 2k=2 .k C 12 /B1 .N /.
(10.6)
The question is: How can we guarantee the validity of the last condition? Now is the
moment for choosing C . Remember that the LLL-reduced basis of is almost orthogonal. Heuristically, this implies that the volume of the parallelepiped formed by
, : : : , bk the fundamental volume of the lattice is approxithe vectors b0 , b1Q
mately equal to kiD0 jbi j. But the volume of is equal to the absolute value of the
Q
determinant of M , which is C , hence kiD0 jbi j is of the size of C . On the other
hand, all vectors of an LLL-reduced basis have lengths of same size, hence we expect that jb0 j is of the size of C 1=.kC1/ . Therefore, if we had chosen the integer C
somewhat larger than
1 kC1
B1 .N /kC1 , 013
(10.7)
2k.kC1/=2 k C
2
it would be reasonable to expect that condition (10.6) is fulfilled; if not, we would try
a new choice for C , somewhat larger than the previous one; in practice this always
works. Summing up, we have proved the following proposition:
Proposition 10.1.1. Let D n0 C n1 1 C nk k , where 1 , : : : , k are specific
real numbers and n0 , n1 , : : : , nk are unknown integers satisfying the two conditions
N D max0ik jni j B1 .N / and jj
1 exp.
2
4 N 2 /, where we understand
that B1 .N / is a specific very large integer and
1 ,
2 ,
4 specific positive real numbers.
Choose a positive integer C of size (10.7) and consider the lattice defined at
the beginning of this section. Compute an LLL-reduced basis .b0 , b1 , : : : , bk /. If b0
satisfies (10.6), then N satisfies (10.5).
The very important result of this proposition is that the new upper bound for N is
of the size of
1=2
k
1
1=2
,
4
log 2 C log k C
C log B1 .N /
2 C log
1 C .k C 1/
2
2
p
i.e. the size of the new upper bound for N is somewhat larger than log.B1 .N // !
125
10.2
Examples
In this section we apply the reduction technique described in Section 10.1 to the linear form L.P / already obtained in each of Sections 9.3, 9.4, 9.5 and 9.6. In order to
compute an LLL-reduced basis for the lattice in each of our examples, we can use
e.g. either of the following:
The routine LLL of MAPLE with input data the column vectors of M .
The routine qflll of PARI with input data M and option 1, to declare that the
routine must view its data as integers. Actually, qflll returns the transformation
matrix from the reduced to the initial basis, so that the LLL-reduced basis is formed
by the column vectors of the matrix product M qflll.M /.
MAGMA
As b0 we can take either the first row of the matrix [Link] / (MAGMA), or the first
column of the matrix M qflll.M / (PARI), or the first vector of the vector array
returned by LLL.c1 , c2 , c3 / (MAPLE), where c1 , c2 , c3 are the columns of M regarded
as vectors.
10.2.1
Weierstrass equation
We will make use of some of the numerical results that we obtained in Section 9.3. We
have k D 2, d D 1, r0 D 2, D 2, D 3 and we have computed all the parameters
involved in the definition of
1 and
2 (see (10.3)). Concerning
3 , this must be chosen
so that M 2
3 N 2 (cf. just above (10.2)). In our case, N M C D 2M C 3,
49
N 2 . For the right-most inequality it suffices to assume
hence M 2 .N 3/2 =4 400
49
that N 10, which certainly holds if M 10. Thus,
3 D 400
and now, in view of
(10.3), we can choose
1 D 6.73841,
2 D 6.2815,
4 D 0.0748.
2`1
2`2
2
C n2
D n0 C n1 1 C n2 2 ,
L.P / D n0 C n1
!1
!1
!1
D n0 C n1 .0.85484 : : :/ C n2 .0.65652 : : :/,
126
we do them with a precision of 150 decimal digits. Then C is the integer resulting
from C if we truncate its decimal digits. In this way we calculate
8548414188
0929152369
C 1 D
,
130 digits
C 2 D
6565274509
9728533696
130 digits
Z3
T 2 Etors .Q/
have the property that P C has integer coordinates. In other words, given a point P E D
.x, y/ as above, we check whether or not the point P C D .u, v/ with4
1
7
5
uDx , vDy x
12
2
24
has integer coordinates. A simple computer program can do the job very easily. Our
computational results are collected in Table 10.1.
4
127
P E D .x, y/
12 5 P C D .u, v/
1, 1, O
1, 0, O
1, 1, O
0, 2, O
0, 1, O
0, 1, O
0, 2, O
1, 1, O
1, 0, O
1, 1, O
1, 0, .41=6, 0/
0, 1, .41=6, 0/
0, 1, .41=6, 0/
1, 0, .41=6, 0/
1, 2, .31=12, 0/
0, 3, .31=12, 0/
0, 1, .31=12, 0/
0, 0, .31=12, 0/
0, 1, .31=12, 0/
0, 3, .31=12, 0/
1, 2, .31=12, 0/
1, 1, .113=12, 0/
0, 1, .113=12, 0/
0, 0, .113=12, 0/
0, 1, .113=12, 0/
1, 1, .113=12, 0/
.293=12, 225=2/
.79=12, 4/
.1733=12, 3465=2/
.161=12, 36/
.67=12, 15=2/
.67=12, 15=2/
.161=12, 36/
.1733=12, 3465=2/
.79=12, 4/
.293=12, 225=2/
.3233=12, 4420/
.581=12, 663=2/
.581=12, 663=2/
.3233=12, 4420/
.2417=12, 2856/
.73613=12, 960925=2/
.173=12, 85=2/
.31=12, 0/
.173=12, 85=2/
.73613=12, 960925=2/
.2417=12, 2856/
.269=12, 195=2/
.43=12, 13=2/
.113=12, 0/
.43=12, 13=2/
.269=12, 195=2/
.24, 100/
.7, 1/
.144, 1805/
.13, 29/
.6, 5/
.6, 10/
.13, 43/
.144, 1660/
.7, 7/
.24, 125/
.269, 4285/
.48, 307/
.48, 356/
.269, 4555/
.201, 2957/
.6134, 477395/
.14, 50/
.3, 1/
.14, 35/
.6134, 483530/
.201, 2755/
.22, 109/
.4, 8/
.9, 5/
.4, 5/
.22, 86/
10.2.2
Quartic equation
1 D 1.8256,
2 D 5.0395,
4 D 0.274518.
128
Choice of C : According to (9.11), B1 .N / D 2.3 1088 and (10.7) suggest that the
integer C be of the size of 1.2162 10448 . We choose C D 10450 and ask our computer
to work with a precision of 460 decimal digits.
The linear form to which we will apply the reduction process is
1
`1
`2
`3
`0
D
C n2
C n3
C n4
L.P / D n0 C n1
!1
!1
!1
!1
!1
D n0 C n1 1 C n2 2 C n3 3 C n4 4
D n0 C n1 .0.03977 : : :/ C n2 .0.25601 : : :/
C n3 .0.41647 : : :/ C n4 .0.48622 : : :/.
In the present case is a sublattice of Z5 and
0
1
0
0
0
B 0
1
0
0
B
M D B
0
0
1
0
B
@ 0
0
0
1
C 1 C 2 C 3 C 4
0
0
0
0
C
1
C
C
C,
C
A
where
3977480557
4553315116
C 1 D
,
C 2 D
2560106453
4697144898
C 3 D
4164734877
6333939317
,
C 4 D
4862224793
9803530989
449 digits
450 digits
450 digits
450 digits
b0 D
C
A,
82327895474728509176238354944130475886404355810644437312762526348821779131423108393917712
15471125794154247726983556666031672982938641013672184463414407364168745466146688236207559
B 111795263890064364945667733492720885363502478756971653430085800804576405203867180516532531
@
244350433420305070938288260659067994439387630719216219811147801739708020305288497923338381
491414271890840780959306036704856290884457656764343679435212935777287084290581917008213944
with approximate length 5.664 1089 , hence of the size of C 1=.kC1/ D 1090 , as expected. We check that b0 satisfies the condition (10.6), therefore, by Proposition 10.1.1,
N satisfies (10.5), which implies that N 55. We set now B1 .N / D 55 and repeat
the process, by choosing C D 1016 . We obtain
0
1
151
B 77 C
B
C
C
new b0 D B
B 378 C
@
62 A
984
which satisfies (10.6) and Proposition 10.1.1 asserts that we can apply relation (10.5).
Doing so we obtain the new upper bound N 11, hence also M 11, and we check
129
which points
P E D m1 P1 C m2 P2 C m3 P3 ,
have the property that P C has integer coordinates. In other words, we check for which
points P E D .x, y/ as above, the point P C D .u, v/ with5
4.3x C 2/ 216x 3 C 216x 2 59 108y 2 C 108y
.u, v/ D
,
(10.8)
3.2y 1/
27.2y 1/2
has integer coordinates. Our computational results, obtained under the restriction y
1=2, are collected in Table 10.2.
Table 10.2. All points P E D n1 P1E C n2 P2E C n3 P3 with P C D .u, v/ 2 Z Z.
n1 , n2 , n3
1, 1, O
1, 0, 1
1, 0, 0
1, 0, 1
1, 2, 0
0, 2, 0
0, 0, 1
0, 0, 1
0, 1, 0
12 5
P E D .x, y/
.19=3, 27=2/
.43=12, 13=8/
.2=3, 1=2/
.13=3, 11=2/
.1433=507, 2343=4394/
.262=75, 109=250/
.8=3, 3=2/
.8=3, 3=2/
.5=3, 5=2/
P C D .u, v/
.1, 2/
.4, 25/
.0, 1/
.2, 7/
.130, 26701/
.130, 26701/
.2, 7/
.4, 25/
.1, 2/
130
.u, v/ D .0, 1/ does not appear is because it corresponds to the zero point O 2 E, as
a computation in the function field Q.C / shows.8
We have thus proved the following proposition:
Proposition 10.2.2. All integer solutions of the equation 52 u4 C 12 u3 C u2 C 1 D v 2
are .u, v/ D .0, 1/ and those listed in the right-most column of Table 10.2.
Singular values of the birational transformation. We must check the finitely many
points .x, y/ 2 E, such that y is a zero of the denominator of the rational transformation E ! C separately, as to whether they correspond to a point on C which otherwise
we might miss. In our example, we have to check only the point .x, y/ D .2=3, 1=2/.
Below we explain why this point is mapped to the point .2=9, 81=83/ 2 C , hence to
no integer solution of our Diophantine equation.
The finitely many points .u, v/ 2 C , such that u is an integer and, at the same time,
a zero of a denominator of the rational transformation C ! E can be trivially found.
In this sense, it is not really necessary to find out to which points .x, y/ 2 E such
points .u, v/ are mapped; nevertheless, below we discuss this issue as well.
First, we work in the function field Q.E/ viewing u and v in (10.8) as functions. In
view of y 2 D x 3 .31=3/x685=108 we have y 2 1=4 D .3xC2/.9x 2 6x89/=27,
hence
27
2y C 1
3x C 2
D
2
2y 1
4 9x 6x 89
and, therefore, by (10.8), the value of the function u at the point .2=3, 1=2/ is
u.2=3, 1=2/ D 2=9. Next, we consider the function v. By (10.8) we have
8
8.3x 1/.3x C 2/2 27.2y 1/2
3x C 2 2
D .3x 1/
1,
vD
27.2y 1/2
27
2y 1
from which we see that v.2=3, 1=2/ D 83=81. This proves that, in the birational
transformation from E to C , the point .x, y/ D .2=3, 1=2/ 2 E is mapped to the
point .u, v/ D .2=9, 83=81/ 2 C .
Next, we will show that, in the inverse transformation from C to E, the point
.u, v/ D .0, 1/ 2 C is mapped to O 2 E.
We work in the function field Q.C /. By v 2 1 D .5=2/u4 .1=2/u3 C u2 we have
5u2 C u C 2
vC1
D
.
u2
2.v 1/
Then, by this relation and (10.9),
15u2 C 3u C 5 C v
vC1 1
D
,
xD2
C
u2
3
3.v 1/
yD
8
2 5u2 C u C 2 2 1
20u2 C 3u C 4 C uv C 4v
4 vC1 2 1
2 C C D
C C D
.
u u
u 2
u
v1
u 2
2u.v 1/
131
Taking into account these expressions of x and y and writing projectively the functions
u, v and x, y, using u D U=W , v D V =W and x D X=Z, y D Y =Z, we can express
the rational map from the projective curve C to the projective curve E as follows:
.U : V : W / 7! 2.15U 2 C 3U W C 5W 2 C V W /U :
3.20U 2 C 3U W C 4W 2 C U V C 4V W /W :
6U W .V W / ,
from which we see that the point .0 : 1 : 1/ 2 C is mapped to the point .0 : 1 : 0/ D
O 2 E.
10.2.3
1 D 0.7496,
2 D 9.5675,
4 D 0.1139.
Choice of C : According to (9.5), B1 .N / D 1.71041 and (10.7) suggest that the integer
C be of the size of 6.14125 10125 . We choose C D 10126 and ask our computer to
work with a precision of 140 decimal digits.
The linear form to which we will apply the reduction process is
1
`1
`2
L.P / D n0 C n1
D
C n2
D n0 C n1 1 C n2 2
!1
!1
!1
D n0 C n1 .0.52627 : : :/ C n2 .0.79791 : : :/.
The lattice is a sublattice of Z3 and
0
1
1
0
0
1
0 A,
M D @ 0
C 1 C 2 C
where
C 1 D
5262784452
0127775724
,
126 digits
C 2 D 7979116877
4671848557
126 digits
132
The vector b0 satisfies the condition (10.6), therefore, by Proposition 10.1.1, N satisfies (10.5), which implies that N 42. We set now B1 .N / D 42 and repeat the
process, by choosing C D 108 . We obtain
0
1
96
37 A .
new b0 D @
192
The new b0 satisfies (10.6), hence we obtain the new upper bound N 14.
Given a point .x, y/ 2 E we obtain a solution .U , V , W / to the system of simultaneous Pell equations as follows:
From .x, y/ we obtain .u, v/ using (6.5) and (6.6), where is as in (7.12), not as in
Lemma 6.1.9
We use (7.15) to obtain the value of W and (7.16) to obtain V and U .
VD
W D
T 2 Etors .Q/,
have the property that P E D .x, y/ is not a zero of the denominator of the above
expressions of U , V , W and maps to an integer solution .U , V , W /. Our results are
summarised in Table 10.3.
The exceptional values of .x, y/. Now we must check the points .x, y/ 2
def
E.Q/ which are zeros of the denominator q.x, y/ D 908119808 166320y
28404x 2 27y 2 C 54x 3 in the expressions of U , V , W above Table 10.3. Computing Resy .q.x, y/, f .x/ y 2 / we find out that it is the product of a quartic irreducible
polynomial in x times .3x 1052/2 , hence the only points .x, y/ 2 E.Q/, for which
q.x, y/ D 0, are .1052=3, 3080/.
133
P E D .x, y/
12 3
1, 0, O
0, 1, O
1, 1, O
1, 0, .1048=3, 0/
0, 1, .1048=3, 0/
1, 1, .1048=3, 0/
1, 0, .128=3, 0/
0, 1, .128=3, 0/
1, 1, .128=3, 0/
1, 0, .920=3, 0/
0, 1, .920=3, 0/
1, 1, .920=3, 0/
P C D .u, v/
.1976=3, 14784/
.1052=3, 3080/
.1304=3, 6272/
.848=9, 101024=27/
.1352=75, 202048=125/
.64=3, 2624/
.376=3, 4032/
.880=3, 3360/
.664=3, 4224/
.2396=3, 20664/
.12728=3, 275520/
.4979=3, 66297/
.17, 5, 13/
.5, 1, 1/
.17, 5, 13/
.17, 5, 13/
.5, 1, 1/
.17, 5, 13/
.17, 5, 13/
.5, 1, 1/
.17, 5, 13/
.17, 5, 13/
.5, 1, 1/
.17, 5, 13/
We compute10 the Puiseux series of f .x/ y 2 around x D 1052=3. This is equiv/ y 2 around t D 0. We find
alent to computing the Puiseux series of f .t C 1052
3
y1 .t / D 3080
211
2603 2
tC
t C O.t 3 /
5
22000
211
2603 2
t
t C O.t 3 /.
5
22000
Inserting x D t C 1052=3, y D yi .t / (i D 1, 2) in the expressions of U , V and W and
writing the resulting expressions as series in t , which we denote by Ui .t /, Vi .t /, Wi .t /,
we obtain
55389191
25665908007097 2
18955
t
t C O.t 3 /
U1 .t / D
2603
948585260
76050356898824000
y2 .t / D 3080 C
V1 .t / D
W1 .t / D
4903
3893357
14906658779053
t
t 2 C O.t 3 /
2603 189717052
167310785177412800
11297
18587273
3883710777647
C
tC
t 2 C O.t 3 /
2603
298126796
23901540739630400
and
10 All
U2 .t / D 5
1
27 2
t
t C O.t 3 /
140
616000
V2 .t / D 1
1
23
t
t 2 C O.t 3 /
308
1355200
134
W2 .t / D 1
1
139 2
tC
t C O.t 3 /.
44
1355200
1 D 0.7032,
2 D 2.3588,
4 D 0.0836.
Choice of C : According to (9.6), B1 .N / D 3.11041 and (10.7) suggest that the integer
C be of the size of 3.724 10126 . We choose C D 10127 and work with a precision of
140 decimal digits.
The linear form to which we will apply the reduction process is
`1
`2
1
C n2
D n0 C n1 1 C n2 2
D
L.P / D n0 C n1
!1
!1
!1
D n0 C n1 .0.25260 : : :/ C n2 .0.03262 : : :/.
The lattice is a sublattice of Z3 and
0
1
1
0
0
1
0 A,
M D @ 0
C 1 C 2 C
where
2526017522
97484775902
C 1 D
,
128 digits
C 2 D
3262684039
70180087887
.
127 digits
135
The vector b0 satisfies the condition (10.6), therefore, by Proposition 10.1.1, N satisfies (10.5), which implies that N 48. We set now B1 .N / D 48 and repeat the
process, by choosing C D 108 . We obtain
0
1
329
new b0 D @ 150 A .
175
The new b0 satisfies (10.6), hence we obtain the new upper bound N 13.
We check, therefore, which points P , where
P E D m1 P1 C m2 P2 ,
have the property that P E D .x, y/ maps via the transformation (8.7) to a point
.u, v/ 2 C with integer coordinates. The computer search does not include singular
points .x, y/ that are zeros of the denominators in (8.7). For non-singular points our
search is summarised in Table 10.4.
Table 10.4. All points P E D n1 P1E C n2 P2E with P C D .u, v/ 2 Z Z.
12 3 n1 , n2
P E D .x, y/
P C D .u, v/
0, 1
1, 2
1, 1
.68292, 35626473=2/
.182145033=58564, 969834918357=14172488/
.2710, 81029=2/
.0, 0/
.243, 3/
.0, 0/
The exceptional values of .x, y/. The irreducible factors of the denominators in the
transformation (8.7) are x C 4365 and x 2169, hence we must check separately
the points .x, y/ D .4365, 315171=2/, .2169, 79947=2/ 2 E.Q/. Instead of
computing the Puiseux expansions of f .x/ y 2 around x D 4365, we compute the
Puiseux series of f .t 4365/ y 2 around t D 0. We find
y1 .t / D
y2 .t / D
315171 161784
1785281093 2
tC
t C O.t 3 /
2
1297
19636425657
161784
1785281093 2
315171
C
t
t C O.t 3 /.
2
1297
19636425657
Then, by (8.7),
U.t 4365, y1 .t // D
13740972520149225
C O.t /,
23841086865992
V.t 4365, y1 .t // D
6346275
C O.t /,
313874
136
hence the point .x, y/ D .4365, 315171=2/ 2 E is mapped to the point .u, v/ D
.13740972520149225=23841086865992, 6346275=313874/ 2 C , with non-integer coordinates. We compute also
492331523
C O.t /
U.t 4365, y2 .t // D 1226330361 t 2 865647 t 1 C
1682209
16679
V.t 4365, y2 .t // D 35019 t 1 C
C O.t /,
1297
which shows, if we write the rational map (8.7) projectively, that the projective point
.4365 : 315171=2 : 1/ 2 E is mapped to the projective point .1 : 0 : 0/, the
point at infinity of C , which does not furnish a solution to our Diophantine equation.
Similarly, working with the points .2169, 79947=2/, we find that .2169, 79947=2/
is mapped to the non-integer point
.19460428805563741=56309567282892, 1324909=238854/ 2 C ,
and .2169, 79947=2/ is mapped to the point at infinity.
We have, thus, proved the following proposition:
Proposition 10.2.4. The equation 3v 5 C 3uv 3 271uv 3u2 D 0 has .u, v/ D
.0, 0/, .243, 3/ as its only integer solutions.
Chapter 11
a1 , a2 , a3 , a4 , a6 2 Z
(11.1)
11.1
In this section we will study the curve C defined by (11.1) from the point of view of
its formal group. Our exposition is based mainly on [45, Chapter IV].
We fix a (rational) prime p and, as usual, we denote by Qp the field of p-adic
numbers (the completion of Q under the metric induced by the p-adic absolute value
j jp ) and by Zp the ring of p-adic integers, i.e. those x 2 Qp with jxjp 1. The
138
(only) maximal ideal of Zp is pZp which, for simplicity, we will denote by M. The
natural homomorphism Zp ! Zp =M is denoted by t ! tQ. The residue field Zp =M
is isomorphic to Fp D Z=pZ.
The curve CQ =Fp defined by v 2 C aQ 1 uv C aQ 3 v D u3 C aQ 2 u2 C aQ 4 u C aQ 6 is the
reduction of C mod p. The reduction mod p homomorphism
[Link] / 3 Q 7! QQ 2 CQ .Fp /
is defined as follows. Any point Q 2 [Link] / can be written in projective coordinates
.U0 : V0 : W0 / with all three coordinates in Zp and not all zero. Then QQ D .UQ0 : VQ0 :
WQ 0 / is a point belonging to CQ .Fp /.
The curve CQ =Fp may have singular points but, anyway, the subset of the nonsingular points of CQ .Fp /, denoted by Cns .Fp /, is a group. We define
C0 .Qp / D Q 2 [Link] / : QQ 2 Cns .Fp /,
which is a subgroup of [Link] /, and the kernel of the reduction
Q
C1 .Qp / D Q 2 [Link] / : QQ D O;
for the above we refer the reader to [45, beginning of Section VII.2 and Proposition
2.1]. Note that OQ 2 Cns .Fp /, hence C1 .Fp /
C0 .Fp /.
In order to define the formal group associated to C =Qp we need the model of (11.1)
defined by
def
w D z 3 C a1 zw C a2 z 2 w C a3 w 2 C a4 zw 2 C a6 w 3 D h.z, w/.
(11.2)
139
is expressed as a convergent in Qp power series in z and w.z/ D h.z, w.z//; see [45,
Chapter IV, Proposition 1.1]. This makes it possible to express u and v as convergent
power series of Qp as follows:
z
(11.3)
D z 2 a1 z 1 a2 a3 z .a4 C a1 a3 /z 2 2 Qp
w.z/
1
D z 3 C a1 z 2 C a2 z 1 C a3 C .a4 C a1 a3 /z 2 Qp .
v.z/ D
w.z/
u.z/ D
z3 D z3 .z1 , z2 / D z1 z2
a1 C a3 2 C a2 C 2a4 C 3a6 2
2M
1 C a2 C a4 2 C a6 3
(the fact that , and z3 belong to M is not obvious at first glance; cf. [45, page 119])
and, finally,
F.z1 , z2 / D i.z3 .z1 , z2 //
D z1 C z2 a1 z1 z2 a2 .z12 z2 C z1 z22 /
C .2a3 z13 z2 C .a1 a2 3a3 /z12 z22 C 2a3 z1 z23 / C 2 M.
1
2
140
Let
1 if p > 2
D
.
(11.7)
2 if p D 2
Since D 1 for p > 2, the above claim is already proven a few lines above; it suffices
therefore to check for p D 2, when D 2. In this case, if z 2 M2 , then it is easily seen
that 2 .z k =k/ k C 1 for every k 2, which is sufficient for the proof of the claim.
The (p-adic) exponential function on C D CO .M/ is formally defined6 as the unique
power series expC z satisfying
logC .expC z/ D expC .logC z/ D z.
P
By [45, Chapter IV, Proposition 5.5], expC z D z C 1
kD2
every k. In analogy with the logarithmic function:
kD2
zk ,
where k 2 Z for
1
X
k k
expC z D z C
z ,
k
(11.8)
k
k
(11.9)
141
To prove this, it suffices to show that the value of p at each summand of the infinite
sum is at least . We have p .k t k =k/ kp .t / p .k/ k p .k/. On the
other hand,
k
k
k
k
k
p .k/ D
C
,
C < C 2 C D
2
p
p
p
p
p1
k
1
hence, p .z k =k/ > k p1
k. p1
/. The observation that the right-most
side is k=2 if p 3 and k if p D 2 completes the proof of our claim.
In view of the above discussion, we also conclude that:
is a group isomorphism, the inverse isomorphism being the function expC .7 This means
in practice that, if z1 , z2 2 M , then
logC F.z1 , z2 / D logC z1 C logC z2
Lemma 11.1.2.
(i) Let Q D .u.Q/, v.Q// 2 C.Q/ and a prime p such that p .u.Q// < 0. Then,
p .u.Q// is an even number and p .z.Q// D 12 p .u.Q//, hence jz.Q/jp D
1=2
ju.Q/jp .
(ii) If Q D .u.Q/, v.Q// 2 C.Q/ and p .u.Q// < N for some positive N 2 Z,
then, also, p .u.Q// < N .
(iii) For i D 1, 2, let Qi D .[Link] /, [Link] // 2 C.Q/ be such that p .[Link] // < N
for some positive integer N . Then, p .u.Q1 C Q2 // < N .
(iv) For i D 1, : : : , k, let Qi D .[Link] /, [Link] // 2 C.Q/ satisfy p .[Link] // < N
for some positive integer N . Then, p .u.m1 Q1 C C mk Qk // < N for any
integers m1 , : : : , mk .
Proof. (i) In view of Lemma 1.2.2 we can write u.Q/ D U=W 2 and v.Q/ D V =W 3 ,
where U , V , W 2 Z and .U , W / D 1 D .V , W /. If p .u.Q// < 0, then, necessarily, p .W / 1 and, consequently, p .U / D 0 D p .V /. But then, p .u.Q// D
7
142
(11.11)
or, equivalently,
lp .n1 Q1 C C nk Qk / D n1 lp .Q1 / C nk lp .Qk /.
(11.12)
143
It follows that p .um / is strictly less than the value of p at any other term of m .u/;
therefore p .m .u// D m2 p .u/.
144
Now we estimate p . m .u/2 /. The value of p at any term different from the leading
is at least .m2 2/p .u/, while
2 1
p .m2 um
/ D .m2 1/p .u/C2p .m/ .m2 1/p .u/p .u/ D .m2 2/p .u/.
Therefore, p .
2
m .u/ /
2
m .u/ /
145
2
u C /, where
, ,
, are appropriate rational numbers. If S0 is the set of finite
primes dividing the product of the denominators of
, ,
, , then, clearly, in order
to compute all S -integer points .u, v/ on C , it suffices to compute all S [ S0 -integer
points .u0 , v0 / on the minimal model C0 .
Along with the model C we will need to consider a short Weierstrass model
def
E : y 2 D f .x/ D x 3 C Ax C B,
A, B 2 Z
(11.13)
of the same elliptic curve, which we will denote by E. Obviously, everything in Section
11.1 applies if in place of the model (11.1) we have the model (11.13).
A change of variable from C to E is of the form
u D
2 x C ,
v D
3 y C
2
x C ,
(11.14)
with
, ,
, appropriate rational numbers. The numbers
, ,
and can be chosen
subject to the additional property that, in the inverse transformation
.x, y/ D .
2 u
2 ,
3 v
3 u C
3 .
//
D .
0 2 u C 0 ,
0 3 v C
0 u C 0 /,
(11.15)
the coefficients
0 , 0 ,
0 , 0 are integers; in particular
1 2 Z. For example, with
1
a2 C 4a2
a1
a3 C 4a1 a2 12a3
, D 1
,
D , D 1
6
12
2
24
we transform equation (11.1) into the equation (11.13), with
146
2 maxje1 j, je3 j
D
2je1 j
2
.
3 q
/
if q 1
if e2 , e3 2 R
.
if e2 , e3 62 R
jv.P /jq D j
j3q jy.P /jq ,
and
147
We have q .x.P // D q .
0 2 u.P / C 0 / and q .
0 2 u.P // < q .0 /. Indeed, the last
inequality is equivalent to 2q .
/ C q .u.P // < q ./ 2q .
/ which is true since,
by hypothesis, q .u.P // < bq . Therefore,
q .x.P // D q .
0 2 u.P // D 2q .
/ C q .u.P // D 2q .
/ 2
and
Analogously, q .y.P // D q .
0 3 v.P / C
0 u.P / C 0 / and q .
0 3 v.P // is strictly
less than both q .
0 u.P // and q .0 /, because of the relation q .u.P // < bq . For example, q .
0 3 v.P // < q .
0 u.P // is equivalent to 3q .
/ C q .v.P // < q .
/
3q .
/ C q .u.P //, hence equivalent to > q .
/, since q .u.P // D 2
and q .v.P // D 3. We see that the last inequality is true on replacing by
q .u.P //=2.
We thus conclude that
q .y.P // D q .
0 3 v.P // D 3q .
/ C q .v.P // D 3q .
/ 3
and
so that
jz.P C /jq D
j
j2q jx.P /jq
ju.P /jq
E
D
D j
j1
q jz.P /jq .
jv.P /jq
j
j3q jy.P /jq
148
0<
x.P /
dt
f .t /
Z
x.P /
p
23=2 t 3=2 dt D 4 2 .x.P /1=2 X 1=2 /.
149
i D 0. In any case, the relation jx.P /jpi jx.P /jq implies that
log q
(i D 1, : : : , s 1).
i max 0, .2q .
/ bq / log
p
i
M D max jmi j.
1ir
(11.18)
Our main purpose is to obtain an upper bound for M . We will need to investigate
Weil heights of points with S -integer coordinates. First, a few general remarks have
their place.
If x is a non-zero rational number and p is any (finite) prime, p divides the numerator of x if jxjp < 1, while p divides the denominator of x if jxjp > 1. Therefore, for
any S -integer x, the formula (2.20) of the absolute logarithmic height of x becomes
X
log max1, jxjp for x 2 ZS .
h.x/ D
p2S
11 Cf.
150
Also,
Q the product formula (2.10) specialised to a non-zero rational x, takes the form
p jxjp D 1, where we understand that p runs through all (finite) primes and 1. In
the particular case that x is an S -integer we have
Y
Y
Y
1D
jxjp
jxjp
jxjp ,
p2S
p62S
p2S
where the right-most inequality holds because, if p 62 S , then jxjp 1, with strict
inequality if p divides the numerator of x. The inequality above, in particular, proves
the following simple but very useful lemma:
Lemma 11.2.3. Let x be any S -integer and suppose that for some q 2 S we have
jxjq D maxp2S jxjp . Then jxjq 1 and h.x/ s log jxjq .
Remarks.
(1) The isogeny of C and E and relation (11.18) imply that
P C D m1 P1C C mr PrC C T C .
(11.19)
(2) Combining Proposition 11.1.7 and Lemma 11.1.6 we see that for every q 2 S X
1 we can choose a positive integer tq satisfying q .[Link] Pi // < bq for every
i D 1, : : : , r. Then, q .x.P // D q .u.P // 2q .
/ (cf. displayed relation at the
beginning of the proof of Lemma 11.2.1(i)). Therefore, q .x.P // < bq 2q .
/,
which is equivalent to jx.P /jq > q 2q ./bq . Consequently, by Lemma 11.2.1(i),
we conclude that, for every i D 1, : : : , r, the p-adic elliptic logarithms lq .PiC /
and lq .PiE / are meaningful and satisfy jlq .P E /jq D j
jq jlq .P C /jq .
We will need to be more selective about our choice of tq because of our application
of Theorem 11.2.5, which gives an explicit lower bound for a linear form in p-adic
elliptic logarithms. In that theorem the number of elliptic logarithms involved in the
linear form has an essential effect on the size of the bound; the smaller their number,
the better the lower bound. Therefore, we will take tq with the extra property of being
a multiple of the order of the Q-torsion subgroup, so that tq T is the zero point for
every torsion point T . We have to impose one further condition on tq : The original
paper [19] of N. Hirata-Kohno, of which Theorem 11.2.5 is an immediate corollary,
requires, for convergence reasons, that jlq .tq PiE /jq < q
q , where q is defined in
Theorem 11.2.5. Using Lemma 11.2.1(i), Proposition 11.1.4, and Lemma 11.1.2(i),
we see that the last condition can be written as
q
q > j
jq jlq .tq PiC /jq D j
jq [Link] PiC /jq D j
jq [Link] PiC /j1=2
,
q
and this is equivalent to q .tq PiC / < 2q .
/ 2q . Summing up we have the following instructions for choosing tq :
Choice of tq . If q 1 we choose the positive integer tq such that (i) q .tq PiC / <
minbq , 2q .
/ 2q for every i D 1, : : : , r, and (ii) tq T C D O for every torsion
point T C 2 C.Q/.
151
tq lq .P E / if q 1
def
def
Lq .P E / D
and Lq .P C / D tq lq .P C / if q 1,
l.P /
if q D 1
where, in case q D 1, the map l is that of Theorem 3.5.2, and the integer m0 ,
satisfying jm0 j 12 rM , is chosen according to the discussion in Chapter 4 between
relations (4.4) and (4.8).
Then, the above definitions are meaningful, i.e. the q-adic elliptic logarithms
lq .P C / and lq .P E / are meaningful, and
Lq .P C / D m1 lq .tq P1C / C C mr lq .tq PrC /
if q 1
if q 1
if q D 1.
Moreover,
jLq .P E /jq c19 exp.c20 M 2 /,
where
c19
jtq jq e =s
D p =s
2 2e
if q 1
,
if q D 1
(11.20)
c20 D =s
(11.21)
152
Then, jLq .P E /jq D jtq jq jlq .P E /jq D jtq jq jz.P E /jq , by Proposition 11.1.4 and,
analogously, jLq .P C /jq D jtq jq jlq .P C /jq D jtq jq jz.P C /jq . By Lemma 11.2.1,
jz.P C /jq D j
1 jq jz.P E /jq , which proves the relation (11.21).
Now we estimate jLq .P E /jq . We saw above that q .u.P // < bq , therefore, by
. Thus, finally,
Lemma 11.2.1(i), jz.P E /jq D jx.P /j1=2
q
jLq .P E /jq D jtq jq jx.P /j1=2
jtq jq e h.x.P //=2s ,
q
(11.22)
where the right-most inequality is implied by Lemma 11.2.3. Using Propositions 2.6.3
and 2.6.2 we obtain successively,
1
1
O // 1 . M 2 /,
h.x.P // . h.P
2s
s
s
therefore by (11.22),
jLq .P E /jq jtq jq e =s e .=s/M ,
2
as claimed.
(ii) Let q D 1. The hypothesis P C 62 P1 implies jx.P /j > b1 so that, by relation
(3.32) of Chapter 3 and Lemma 11.2.1(ii), we have
Z
p
1 C1 dt
E
2 2 x.P /1=2 .
jL1 .P /j1 D jl.P /j D
p
2 x.P / f .t /
By Lemma 11.2.3 we have h.x.P // s log jx.P /j D s log x.P /. Then,
log x.P /1=2
1
h.x.P //
2s
1
O //
. h.P
s
1
. M 2 /
s
Ifp
we combine this with the above displayed upper bound of jl.P /j we obtain jl.P /j <
2
2 2e =s e .=s/M , as claimed.
Finally, by definition, L1 .P E / D l.P / and, according to Chapter 4, relations (4.4)
through (4.8),
l.P / D m1 l.P1 / C C mr [Link] / C l.T / C m0 !1 ,
where m0 is an appropriate integer satisfying jm0 j 12 rM C 1.
When we solve an elliptic Diophantine equation over Z, Theorem 9.1.3 was the basic tool for obtaining a first, very large, upper bound for M . Now that we are interested
in solutions over ZS , the analogous tool is Theorem 11.2.6, below, which is heavily
based on the following very important theorem.
153
3
if q D 2
q D
.
1=.q 1/ if q > 2
O i /
ai D max1, h, tq2 h.P
H D .q
q max jlq .tq PiE /jq /1 ,
1ir
.i D 1, : : : , r/,
d D max1, 1= log H
and
g D max1, h, log a1 , : : : , log ar , log d .
If M
eg ,
then
jLq .P E /jq exp.c23 log M /,
where
c23 D 24r
2 C3r
.r C 1/2r
2 C9rC4
(11.23)
r
Y
ai .
iD1
Proof. This theorem is an almost straightforward specialisation to our case of N. Hirata-Kohnos result [19] which gives an effective lower bound for the p-adic absolute
value of non-vanishing linear forms in p-adic elliptic logarithms. In the notation of that
result, ui D lq .tq PiE / and expp .ui / D .[Link] / : [Link] / : 1/ (projectively). The nonvanishing of our linear form follows from its definition: Lq .P E / D tq lq .P E / 0
by Proposition 11.1.4.
Theorem 11.2.6. With the assumptions and notations of Theorem 11.2.4, the following holds:
Let q D 1. Referring to the model E in (11.13), follow the instructions in Preparatory to Theorem 9.1.2 , page 103 to the linear form L1 .P E / in order to compute
hE and H0 , H1 , : : : Hr , having in mind that, in the notation of the instructions,
L.P / D L1 .P E / and in (9.2), k D r, d D 1, ni D mi for i D 0, 1, : : : , r,
`i D [Link] / for i D 1, : : : , r. Also, let and be as in Propositions 2.6.2, and 2.6.3,
respectively. Then, either
def
D 32 r0 ,
remind that r0 is the lcm of the orders of non-zero points of Etors .Q/.
154
or
M 2 < 32 s log 2 C C c13 [Link].M C / C 1/.log log.M C / C 1 C hE /rC2 ,
(11.24)
where
r
Y
2
2
Hi .
c13 D 2.9 106rC12 42.rC1/ .r C 2/2r C13rC23.3
iD0
Let q 1. If M e g , then
M 2 s [Link] jq / C C s c23 log M ,
(11.25)
155
For reasons explained at the very beginning of Section 11.2, we choose to work
with the linear form Lq .P C /, rather than with Lq .P E /. On the other hand, the elliptic
curve in N. Hirata-Kohnos Theorem [19] is defined by a short Weierstrass model E
with integer coefficients. Therefore we will have to pass from the model C to the
model E and vice versa; in our transition from one model to the other, the basic tool
is Lemma 11.2.1.
Simplifying the notation, we rewrite the linear form Lq .P C / as follows:
L D Lq .P C / D m1 1 C C mk r ,
def
c21 D q .
/
log c19
,
log q
c22 D
c20
.
log q
(11.27)
Forgetting now the specific meaning of the i s, except that they all belong to M ,
and the specific values of c21 and c22 , except that c22 > 0, we will solve the following
computational problem:
Problem. Let 0 , 1 , : : : , r 2 M and consider the linear form
L D m1 1 C C mr r ,
where m1 , : : : , mr 2 Z. Let M D max1ir jmi j and assume the following:
(i) M Mq for some specific large bound Mq .
(ii) q .L/ c21 C c22 M 2 , for some explicit parameters c21 and c22 > 0.
(iii) (Without loss of generality) q .r / D min1ir q .i /.
Under the assumptions (i)(iii), find an upper bound of M considerably smaller
than Mq .
For any 2 Zq and any integer n > 1 we denote by .n/ the n-th rational
P approxii
mation of . Thus, .n/ 2 Z and q . .n/ / > n. In other words, if PD 1
iD0 ai q ,
n
where ai is an integer with 0 ai < p for every i 0, then .n/ D iD0 ai q i .
First we consider the case r 2. We put
i D
r
LDD
r1
X
iD1
i
r
.i D 1, : : : , r 1/,
mi i C mr ,
.n/
D
r1
X
mi i.n/ C mr .
iD1
156
1
C
C
C.
C
A
Proposition 11.3.1.
0
q ./ n
1
m1
B . C
B . C 2 .n/ .
@ . A
mr
.n/
Proof. Suppose first that q ./ n. From q .i i / > n for every i D 1, : : : , r,
we see that q . .n/ / n, hence q ..n/ / n. Since .n/ 2 Z, it follows that
P
.n/
C mr D xr q n for some xr 2 Z. Consequently,
.n/ D r1
iD1 mi i
0
1 0
1
m1
m1
B . C B . C
B . C B . C
C B
C
Z .n/ B
B . C D B . C,
@ mr1 A @ mr1 A
xr
mr
1
C
C
C.
C
A
P
.n/
It follows that mi D xi for i D 1, : : : , r 1 and mr D r1
C xr q n . ConseiD1 xi i
quently,
r1
r1
X
X
.n/
.n/
xi i C .
xi i C xr q n / D xr q n ,
.n/ D
iD1
iD1
157
n 1 C q .r / c21
.
c22
158
the volume formed by the initial basis of the lattice which, in turn is equal to the
determinant of Z .n/ . But this is q n , thereforepwe expect that jb1 j is of the size of q n=r .
Since we want that jb1 j satisfies jb1 j > 2r=2 r Mq , we must choose n so that
&
'
p
r log.2r=2 rMq
n
.
(11.28)
log q
11.4 Example
We will compute all S -integral points on
C : v 2 C v D u3 7u C 6,
(11.29)
for S D 2, 3, 5, 7, 1. Note that C is a minimal model. The rank over Q is 3 and the
torsion subgroup is trivial. A MordellWeil basis is
P1C D .1, 0/,
D , D 0,
D 0, D
2
2
and
P1E D .4, 4/, P2E D .8, 4/, P3E D .0, 20/.
In Proposition 2.6.3, in place of the curve D we take our present curve C and
we compute 2.87085138. Also, applying Proposition 2.6.2, we obtain
0.48599751, therefore c20 D 3=5. We need and for calculating the parameters
c19 (depending on the prime q 2 S ) of Theorem 11.2.6.
First we compute all points P C 2 P; cf. Lemma 11.2.2. We compute b1
24.09944 < 24.1 and .b2 , b3 , b5 , b7 / D .5, 1, 1, 1/. Following the steps of the
constructive proof of Lemma 11.2.2, it is easy to see that, if P C 2 P, then the height
of x.P / (maximum of the numerator and denominator of jx.P /j) is less than 121462.
Using the routine Points of MAGMA, we explicitly compute 274 points only 119 of
which have S -integral coordinates; thus P is a subset of an explicit set P consisting
of 119 S -integral points.
159
Now we assume that P C 62 P so that we can apply Theorem 11.2.6. We denote by q that prime in S for which jx.P /jq D maxp2S jx.P /jp and we distinguish cases according to q D 2, 3, 5, 7, 1. In the case that q 1, the routine
pAdicEllipticLogarithm of MAGMA is used for computing the lq -values.
q D 2. The instructions for the choice of tq on page 150 lead us to t2 D 80, so that
L2 .P C / D m1 l2 .80P1C / C m2 l2 .80P2C / C m3 l2 .80P3C /.
Theorem 11.2.6, in particular relation (11.25), provides us with the upper bound
M M2 D 1029 .
Now we reduce this large upper bound M2 . We need first to calculate the parameters
c21 and c22 according to (11.27); we actually need an approximation (not a high one)
from below:
c21 3.17164739, c22 0.865617.
According to (11.28) we must work with a 2-adic precision at least O.2316 /. We prefer to work with a somewhat larger precision, choosing n D 320 in order to construct
the matrix Z .n/ . A first low precision computation shows that 2 .l2 .80PiC // D 6 for
i D 1, 2 and 2 .l2 .80P3C // D 5. Therefore we put i D l2 .80PiC / for i D 1, 2, 3,
so that i D i =3 for i D 1, 2. Using MAGMA we compute
.320/
1
D 224782
369609
2,
96 digits
.320/
2
D 163742
694465
2.
96 digits
q D 3. The P
instructions for the choice of tq on page 150 lead us to t3 D 7, so that
L3 .P C / D 3iD1 mi l3 .7PiC /.
13 See
160
161
that, in the first reduction step, the upper bound falls to M1 D 51 and the next step
reduces it further to M1 D 13. As expected, all 119 points in P (see the beginning
of this section) are rediscovered by our search, while 9 points from those found do
not belong to P.
Summing up, we conclude that M maxM2 , M3 , M5 , M7 , M1 D 13 and we
check all points P C D m1 P1C C m2 P2C C m3 P3C with M D max1i3 jmi j for
S -integrality. The results of our search is summarised in Table 11.1.
Table 11.1. All S -integral points P C D m1 P1C C m2 P2C C m3 P3C on the curve (11.29),
S D 2, 3, 5, 7, 1.
m1 , m2 , m3
4, 2, 0
4, 0, 1
3, 2, 0
3, 2, 2
3, 1, 1
3, 1, 3
3, 0, 1
3, 1, 0
3, 1, 1
2, 3, 2
2, 3, 2
2, 2, 0
PC
1541761
, 1882859327
153664
60236288
11948 176534
3969 , 250047
6169 641312
, 531441
6561
2759 60819
1024 , 32768
391 7564
,
25 125
19849 1787743
, 729000
8100
4537 305425
, 216
36
47 9191
256 , 4096
.816, 23309/
2985362173625
,
4096
262144
6142 480700
,
81
729
207331217
.406, 8181/
P C
1541761
1822623039
,
153664
60236288
11948
73513
3969 , 250047
6169 109871
,
6561 531441
2759 93587
1024 , 32768
391 7689
, 125
25
19849 1058743
,
8100 729000
4537 305641
, 216
36
47
256 , 132871
4096
.816, 23310/
207331217
, 2985362435769
4096
262144
6142 481429
,
81
729
2, 1, 1
13 804
,
49 343
848961 782275425
, 175616
3136
33 81
, 64
16
4 62
, 27
9
2, 1, 2
.93, 897/
.93, 896/
2, 2, 1
2, 2, 4
2, 1, 0
.406, 8180/
13
, 1147
49
343
848961 782099809
, 175616
3136
33 17
,
16 64
4 35
,
9 27
162
PC
m1 , m 2 , m3
43
, 132
49 343
2, 0, 1
2, 0, 0
2, 0, 1
2, 0, 2
2, 1, 1
.14, 52/
2, 1, 0
2, 1, 1
P C
43 475
, 343
49
33
74 , 8
, 209 /
. 106
49 343
, 329
69
25
125
25 64
,
9 27
24 18
,
25 125
1, 1, 0
.3, 0/
2, 2, 0
2, 3, 1
1, 3, 0
1, 2, 0
1, 2, 1
1, 2, 2
1, 1, 1
1, 1, 2
1, 0, 2
, 552
343
33304 5642994
,
15625 1953125
221 2624
,
49 343
8159 390925
2916 , 157464
1343 36575
576 , 13824
26 101
, 125
25
25 111
, 8
4
66 252
25 , 125
31 116
,
9 27
.3, 1/
.4, 7/
17
, 39
16
64
151
1845
,
64
512
1, 0, 1
.8, 21/
.8, 22/
1, 0, 0
.1, 1/
.1, 0/
1, 0, 1
1, 0, 2
.3, 4/
26 28
,
9 27
.4, 6/
, 25
16
64
1333
,
151
64 512
204
69
,
25 125
25 91
, 27
9
24 143
, 125
25
625 14839
, 512
64
17
74 , 25
8
49
625
1, 1, 1
2, 1, 3
.14, 51/
106
, 15351
64
512
33304 7596119
15625 , 1953125
221 2967
, 343
49
8159 233461
2916 , 157464
1343 50399
576 , 13824
26
24
, 125
25
25 119
, 8
4
66 377
25 , 125
31 143
, 27
9
2, 1, 2
.3, 3/
26
55
,
9
27
163
PC
m1 , m2 , m3
P C
58 559
,
81 729
9 7
,
4 8
58
, 1288
81
729
9 15
, 8
4
1, 1, 2
1, 1, 1
1, 1, 0
.2, 3/
.2, 4/
1, 1, 1
.11, 35/
.11, 36/
.52, 374/
.52, 375/
1, 2, 1
1, 2, 0
1, 2, 1
1, 3, 3
1, 3, 1
1, 3, 1
0, 3, 2
0, 2, 1
9
, 133
16
64
172 1079
, 729
81
338776 197184149
, 3375
225
1219
9828
,
625
15625
13961 1648791
,
100
1000
1793 68991
,
256 4096
68 1079
49 , 343
, 69
16 64
172 350
,
81 729
338776 197180774
, 3375
225
1219
5797
,
625
15625
13961 1649791
,
100
1000
1793 73087
, 4096
256
68 1422
49 , 343
.21, 95/
.21, 96/
0, 2, 0
7
9
0, 2, 2
9
27
114 720
, 343
49
, 17
9 27
114 377
,
49 343
0, 1, 2
.342, 6324/
.342, 6325/
0, 2, 1
,
44
7
, 21
4
8
1
1
0, 1, 1
, 13
8
0, 1, 0
.2, 1/
.2, 0/
0, 1, 1
.1, 3/
.1, 4/
0, 1, 2
0, 0, 3
0, 0, 2
0, 0, 1
.37, 224/
1525
, 438957
2401
117649
49
32
, 125
25
.0, 2/
.37, 225/
, 321308
1525
2401 117649
49
93
, 125
25
.0, 3/
List of symbols
Reference
Chapter 1
gD0
Uij , Vij
Eq. (1.2)
Def. 1.1.1
C,E
p. 4
X , Y, U, V
PC,PE
.x.P /, y.P //
.u.P /, v.P //
a1 , a2 , a3 , a4 , a6
Etors
r
P1 , : : : , Pr
P
T
p. 4 ff.
homogenisation of g D 0
birational transformation between curves
Ci , Cj
birationally equivalent models usually representing elliptic curve with E short Weierstrass
model
birational transformation between curves
E, C
representatives of point P on models C , E
coordinates of point P on model E
coordinates of point P on model C
coefficients of Weierstrass equation
torsion subgroup of E
rank of elliptic curve
MordellWeil basis
generic point of elliptic curve
generic torsion point of elliptic curve
Eq. (1.6)
Thm. 1.2.1
Eq. (1.7)
166
Symbol
p
p
eL=K .P/
fL=K .P/
dL=K .P/
j jp
NL=K ./
HK ./
HK,fin ./
HK,1 ./
HK,fin .P /
H./
h./
h.P /
hD .P /
O /
h.P
O
hD .P /
H
, , , ,
List of symbols
Reference
Chapter 2
Heights
Sec. 2.1
p. 17
Eq. (2.1)
above Eq. (2.3)
Eq. (2.7)
Eq. (2.9)
Eqs. (2.4), (2.5)
Sec. 2.3
Sec. 2.6
Prop. 2.6.1
Eq. (2.39)
Prop. 2.6.2
Prop. 2.6.3
Sec. 2.4
Def. 2.3.1
Sec. 2.5
Chapter 3
.!1 , !2 /
Sec. 3.1
g2 , g3
e1 , e2 , e3
A, B
g2 , g3
Eq. (3.6)
Fact 3.1.2
below Fact 3.1.2
Sec. 3.2
Fact 3.2.1
Sec. 3.3
167
List of symbols
Symbol
Reference
Def. 3.3.2
"P
Sec. 3.4.4
Sec. 3.5
Thm. 3.5.1
l0
l
Thm. 3.5.2
Eq. (3.31)
Chapter 4
A, B; E
m1 , : : : , mr
Eq. (4.1)
Eq. (4.2)
M
m0
M0
r0
e1 , e2 , e3
Qi
s=t
L.P /
!1
!2
Sec. 3.4
.!1 , !2 /
p. 37
.1 , 2 /
M.a, b/
Chapter 5
Eq. (5.1)
, , ,
E
Prop. 5.1
Thm. 5.2
L.P /
168
Symbol
List of symbols
Reference
Chapter 6
Eq. (6.1)
Q, a, b, c, d , e
A, B
Eq. (6.2)
Eq. (6.3)
E
X,Y
Eq. (4.1)
Eq. (6.4)
U, V
Lemma 6.1
x0
u0
u
x.u/
Eq. (6.8)
Lemma 6.3
Prop. 6.4
Eq. (6.9)
P0
u , c7
L.P /
Eq. (6.11)
Thm. 6.8
Q, , x0 , u0
x.u/, u , P 0 , u
c 7 , L.P /
Table 6.1
Chapter 7
A1 , D1 ; A2 , D2
.U0 , V0 , W0 /
Q
Eq. (7.1)
below Eq. (7.1)
Eq. (7.10)
a, b, c, d , e
A, B
E
u1
u.V /, v.V /
x.u/, y.u/
Eq. (7.12)
Eq. (7.13)
Eq. (6.3)
Eq. (7.17)
p. 72, item (a)
Eqs. (7.18), (7.19)
Eqs. (7.21), (7.22)
169
List of symbols
Symbol
Reference
x.V /, y.V /
x0 , y0
P0E
1 , 2 , 3
V
L.P /
Eq. (7.23)
Eqs. (7.24), (7.25)
Eq. (7.26)
p. 75
Thm. 7.1
c8
Eq. (7.30)
Chapter 8
Eq. (8.1)
Prop. 8.1.1
X , Y; U, V
vi
k,i ,
i , i
B0
x, y
Fact 8.2.1(a)
Fact 8.2.1(d)
Prop. 8.3.2
B1
P0E
B2
Def. 8.3.3
Prop. 8.4.2
c9 ,
B3 , c10 , c11
Prop. 8.7.1
L.P /
Thm. 8.7.2
Chapter 9
L.P /
Eq. (9.2)
`0 , `1 , : : :
"
m00
p. 103
Sec. 9.1
170
List of symbols
Symbol
Reference
k, n0 , n1 , : : : , nk , d
bullets
Eq. (9.2)
p. 99 ff.
pp. 99100
third bullet p. 100
,
$1 , $2 , Q
Eq. (9.6)
Lemma 9.1.1
K, D, jE , hE
H0 , H1 , H2 : : :
c12 , c13 , c14 , c15
c16 , c17 , c18
p. 103
Thm. 9.1.2
Thm. 9.1.3
Chapter 10
Reduction of M
Eq. (10.1)
below Eq. (10.1)
Eq. (10.1)
above Eq. (10.2)
D .P /
1 , : : : , k
n0 , n1 , : : : , nk
N
1 ,
2 ,
4
b0
B1 .N /
1
S
s
C
z.Q/
!.z/
CO .M/ or C
logC
expC
beginning of Ch. 11
Eq. (11.1)
Sec. 11.1
p. 140
Eq. (11.6)
Eq. (11.9)
Eq. (11.7)
171
List of symbols
Symbol
Reference
lp
E
Definition 11.1.3
Eq. (11.13)
, , ,
bq , b1
Pq
p. 146
Lemma 11.2.2
Eq. (11.17)
tq
p. 150
Lq ./
c19 , c20
h, q , H , d , g
a1 , : : : , ar
c23
, , c12 , c13
Thm. 11.2.4
c21 , c22
Mq
Eq. (11.27)
Problem at p. 155
.n/
p. 155
Thm. 11.2.5
Thm. 11.2.6
Z .n/ , .n/
b1
Bibliography
[1] T. M. Apostol, Modular Functions and Dirichlet Series in Number Theory, Graduate
Texts in Mathematics 41, Springer-Verlag, New York Heidelberg Berlin, 1976.
[2] A. Baker, J. Coates, Integer points on curves of genus 1, Proc. Camb. Phil. Soc. 67
(1970), 595602.
[3] I. G. Bashmakova, Diophantus and Diophantine Equations, Dolciani Mathematical Expositions, No. 20, The Mathematical Association of America, New York, 1977.
[4] G. A. Bliss, Algebraic Functions, American Mathematical Society, New York, 1933.
[5] J.-B. Bost, J.-F. Mestre, Moyenne arithmtico-gometrique et priodes des courbes de
genre 1 et 2, Gaz. Math. Soc. Math. France 38 (1988), 3664.
[6] J. W. S. Cassels, Lectures on Elliptic Curves, London Math. Soc. Student Texts 24, Cambridge University Press, Cambridge New York Port Chester Melbourne Sydney, 1991.
[7] J. S. Chahal, Topics in Number Theory, The University Series in Mathematics, Plenum
Press, New York London 1988.
[8] J. Coates, Construction of rational functions on a curve, Proc. Camb. Phil. Soc. 68
(1970), 105123.
[9] H. Cohen, A Course in Computational Algebraic Number Theory, Graduate Texts in
Mathematics 138, Springer-Verlag, New York Heidelberg Berlin, 1995.
[10] H. Cohen, Number Theory-vol. 1: Tools and Diophantine Equations, Graduate Texts in
Mathematics 239, Springer-Verlag, New York Heidelberg Berlin, 2007.
[11] J. E. Cremona, Algorithms for Modular Elliptic Curves, Cambridge University Press,
Cambridge 1997; also available online from the authors home page.
[12] S. David, Minorations de formes linaires de logarithmes elliptiques, Soc. Math. France,
Mmoire 62 (Suppl. Bull. S.M.F.) 123 (1995), fasc. 3, 143 pp.
[13] I. Gal, An efficient algorithm for the explicit resolution of norm form equations, Publ.
Math. Debrecen 56 (2000), 375390.
[14] I. Gal, M. Pohst, On the resolution of relative Thue equations, Math. Comp. 71 (2001),
429440.
[15] J. Gebel, A. Petho, H. G. Zimmer, Computing integral points on elliptic curves, Acta
Arith. 67 (1994), 171192.
[16] L. Hajdu, T. Herendi, Explicit bounds for the solutions of elliptic equations with rational
coefficients, J. Symbolic Comp. 25 (1998), 361366.
174
Bibliography
[17] N. Hirata-Kohno, Formes linaires de logarithmes de points algbriques sur les groupes
algbriques, Invent. math. 104 (1991), 401433.
[18] N. Hirata-Kohno, Linear forms in two elliptic logarithms in the p-adic case, Kyushu
J. Math. 64 (2010), 239260.
[19] N. Hirata-Kohno, Linear forms in p-adic elliptic logarithms, (preprint)
[Link] hirata/[Link].
[20] M. van Hoeij, An algorithm for computing an integral basis in an algebraic number field,
J. Symbolic Comp. 18 (1994), 353363.
[21] M. van Hoeij, Computing parameterizations of rational algebraic curves, ISSAC 94 Proceedings (1994), 187190.
[22] M. van Hoeij, An algorithm for computing the Weierstrass normal form, ISSAC 95
Proceedings (1995), 9095.
[23] A. Hurwitz, A course in the Theory of Functions. Part III-Elliptic Functions Notes (in
Greek) by Nikos Kritikos from a course given at ETH-Zrich in 19161917. Department
of Mathematics, University of Crete, Heraklion 1981.
[24] A. W. Knapp, Elliptic Curves, Mathematical Notes 40, Princeton University Press,
Princeton, New Jersey 1992.
[25] S. Lang, Diophantine approximation on toruses, Amer. J. Math., 86 (1964), 521533.
[26] S. Lang, Elliptic Curves: Diophantine Analysis, Grundlehren Math. Wiss. 231, Springer,
Berlin, 1978.
[27] A. K. Lenstra, H. W. Lenstra, Jr., L. Lovsz, Factoring polynomials with rational coefficients, Math. Ann. 261 (1982), 515534.
[28] D. Masser, Algebraic Numbers, Heights, Transcendental Numbers, Notes for a series of
lectures in the Mathematics Department of the University of Crete, AugustSeptember
1999 (unpublished).
[29] B. Mazur, Modular curves and the Eisenstein ideal, Inst. Hautes tudes Sci. Publ. Math.
47 (1977), 33186.
[30] B. Mazur, Rational isogenies of prime degree (with an appendix by D. Goldfeld), Invent. Math. 44 (1978), 129162.
[31] M. Mignotte, Mathmatiques pour le calcul formel, Presses Universitaires de FranceMathmatiques, Paris, 1989.
[32] R. A. Mollin, Algebraic Number Theory, Chapman & Hall, CRC Press, London, Boca
Raton, 1999.
[33] L. J. Mordell, On the rational solutions of the indeterminate equations of the 3rd and 4th
degree, Proc. Camb. Phil. Soc. 21 (1922), 179192.
[34] L. J. Mordel, Diophantine Equations, Pure and Applied Mathematics 39, Academic
Press, London & New York, 1969.
[35] W. Narkiewicz, Elementary and Analytic Theory of Algebraic Numbers, Second edition, Polish Scient. Publishers Warszawa, Springer-Verlag, Berlin Heidelberg New York
London Paris Tokyo Hong Kong, 1990.
Bibliography
175
176
Bibliography
Index
elliptic
integral
of first kind, 88
Elliptic Logarithm Method, 54
Ellog, vii, ix, x, xii, 54
embedding
complex, 11
real, 11
equation
algebraic, 1
Diophantine
elliptic, vii, 37
elliptic
general, 78
Pell, 69
projective, 1
quartic, 60
Thue
quartic, 69
Weierstrass, 6, 31, 36, 57
short, 6
equivalent
absolute values, 11
models
birationally, 25
exponent
of prime ideal, 9
formula
addition
for }, 33
duplication
for }, 33
integral, 36
product, 12, 14, 16, 18, 150
function
elliptic, 29
178
exponential
p-adic, 140
logarithmic
p-adic, 140
} of Weierstrass, 2949
function-field, 2
genus, 3
group
formal, 137, 140
height
absolute
of projective point, 1517
absolute logarithmic
of algebraic number, 1519
of projective point, 15
canonical, 2328, 106
naive, 20
Nron-Tate, 24
of polynomial, 18
over Q, 17
Weil, 20, 149
homogenisation, 1
infinity
point at, 2
solution at, 2
integer
p-adic, 137
S -, x, 137161
invariant differential, 139
kernel
of reduction, 138
lemma
Gausss, 18
linear form
in elliptic logarithms, viix, 56
linear form L.P /, 58, 66, 76, 9799, 104,
105, 121
logarithm
elliptic, 35, 51, 52, 101, 106
computation of, 5253, 106
p-adic
elliptic, xii, 142
Index
matrix
height-pairing, 25, 106
method
chord-tangent, vii, 6
elliptic logarithm, vii
model
of curve, 1
projective, 1
quartic elliptic, 60
Weierstrass, 6, 57
minimal, 27, 107, 144
short, 56
norm
absolute, 9
ideal, 10
number
p-adic, 137
numerator, 137
pairing
Nron-Tate, 25
Weil, 25
parallelogram
fundamental, 29, 38, 40
period-, 29, 38, 40
period lattice, 29
periods, 29
computation of, 4447
explicit expressions, 4144
fundamental pair of, xii, 29, 37, 44,
102, 106
place, 11
archimedean, 11
non-archimedean, 11
point
at infinity, 6
of a model, 1
zero, 6
polynomial
division-, 143
prime
finite, 137
infinite, 12, 107
Puiseux series, 80
179
Index
ramification index, 10
reduction
mod p, 138
p-adic, 154, 158
reduction method
de Wegers, 154
resolution
effective, ix, 55
explicit, x, 55
routine
CanonicalHeight, 103, 106
ellheight, 103, 106, 107
ellheightmatrix, 106, 107
ellinit, 106
EllipticCurve, 106
EllipticLogarithm, 103, 106
ellminimalmodel, 106
ellpointtoz, 103, 106
HeightPairingMatrix, 106
MinimalModel, 106
omega, 106
pAdicEllipticLogarithm, 144, 159
Periods, 106
Points, 60, 69, 149, 158
puiseux, 81
ratpoints, 60, 149
software
MAGMA,
xiii, 105107
theorem
S. Davids, ixxi, 20, 102, 104
N. Hirata-Kohnos, x, xii, 153
implicit function, 80
Lutz-Nagell, 7
B. Mazurs, 7
Mordell-Weil, 6
z-coordinate, 138









