Let C be a hyperelliptic curve of the genus defined on . The Mordell conjecture proved by Fattings, gives that the set of rational points is finite. The hyperelliptic curves can be subdivided into two types: those with real models and those with imaginary models. Our objective is to compute the rational points in the case of real hyperelliptic curves of genus 3 whose Jacobian have a Mordel-Weil rank equal to 0. This fits into the particular case where has been considered by Chabauty and the techniques developed by Coleman in 1980 , allow us to use the p-adic integration to explicitly computation the set of rational points .
The AMS classification: MSC2020, 14G05, 11G30.
Although there is an algorithm by Jennifer Balakrishnan computing Coleman integrals in the case of real hyperelliptic curves in , this one is unfortunately not available in Sage. We then use the birational transformation of real hyperelliptic curves into imaginary hyperelliptic curves, because in the latter case, Balakrishnan, Bradshaw and Kedlaya described and implemented in Sage the explicit computation of Coleman integrals on imaginary hyperelliptic curves . We start with an overview of hyperelliptic curves and the transition from a real hyperelliptic curve to an imaginary hyperelliptic curve in Section 1. In Section 2, we recall the Chabauty-Coleman method and explicit Coleman integration. In Section 3, we describe our algorithm, which is subdivided into three sub-algorithms: the first one that transforms a real curve to an imaginary curve; the second computes the set of rational points of the imaginary hyperelliptic curve found in Algorithm 3. Algorithm 4 is this of Maria Inés de Frutos Fernaandez and Sachi Hashimoto in  and its implementation in . To the latter, we added a function allowing it to take the hyperelliptic curves of genus 3 and of rank 0 given by non-monic polynomials. Finally, the third algorithm constructs the rational points of the curve C which is real hyperelliptic from the rational points of the imaginary hyperelliptic curve birationally equivalent to C. In Section 4, we present the results obtained on a list of curves taken from the database of hyperelliptic curves of genus 3 .
2. Background on Hyperelliptic Curves
In this section we recall the definition of hyperelliptic curves, the different types of hyperelliptic curves and how to pass from a real hyperelliptic curve to an imaginary hyperelliptic curve. For more detail in this section, we refer the reader to  . A hyperelliptic curve of genus over is defined by an equation:
where f is a polynomial of of degree or . This defines a smooth irreducible algebraic curve in the affine plane. Considering its smooth projective model which is obtained by adding one or two points at infinity, we distinguish two types of hyperelliptic curves which differ by the number of rational points at infinity:
- Real hyperelliptic curves: they are those which have two points at infinity. In this case, the polynomial f is of degree ;
- Imaginary hyprelliptic curves: they are those which have a one point at infinity. In this case, the polynomial f is of degree .
In this article, we are interested in real hyperelliptic curves.
A real hyperelliptic curve C of genus g over is defined by an absolutely irreducible equation of the form:
where . The two points at infinity correspond to the two square roots of the leading coefficient of f. These points are rational if and only if this leading coefficient is a square, so that its square roots are rational numbers.
Let us soppose that is a non-zero square in . We set . The points at infinity of C are given in projective coordinates by and , which are rational .
The set of rational points on C denoted by , is defined by:
Let C be a real hyperelliptic curve of genus g over with equation . If there exists a ramified prime divisor of degree 1 in of C, then the curve C is birationally equivalent to an imaginary hyperelliptic curve of genus g. Indeed, if we denote a rational point of the curve C such that P is equal to its image under the hyperelliptic involution i.e.: that , then P is said to be ramified and in this case the divisor is the ramified prime divisor, then we can make the following change of variable:
In the equation of C, we thus find the equation of .
Therefore P is ramified if . By substituting a and in the equation of C, we get . We conclude that if the polynomial has at least one root in , then the real hyperelliptic curve of equation admits at least one rational ramified point and therefore it is birationally equivalent to an imaginary hyperelliptic curve via the change of variables below.
We are interested in real hyperelliptic curves of genus 3 having at least one rational point of ramification and whose Jacobians have Mordell-Weil rank 0.
3. Background Chabauty-Coleman Method and Coleman Integration
In this section we recall the Chabauty-Coleman method used to compute the rartional points in part 2 of our algorithm. We also give a brief reminder on Coleman integrals. For more details see .
3.1. Chabauty-Coleman Method
Let C be a smooth projective curve over the rational numbers of genus at least 2. From the work of Faltings, we know that is finite, but the proof of Faltings does not explicitly give the set . However, before Faltings’ work, Chabauty considered the following configuration. Let p be a prime number and . We consider the following inclusion:
Let be the p-adic closure of and define . Chabauty proved the following case of Mordell’s conjecture:
Theorem 2.1. Let be a curve of genus such that the Mordell-Weil the rank of the Jacobian J of C over Q is less than g, and let p be a prime number, then is finite.
Chabauty’s result was then reinterpreted and made effective by Coleman, who showed what follows:
Theorem 2.2. Let C be as above and suppose that p is a prime number of good reduction for C. If , then:
The Chabauty-Coleman method uses the theory of p-adic integrations on curves to construct the p-adic integrals of 1-forms on which vanish on and which allow in practice to explain the set of rational points which are contained in the intersection of the set of zeros of the integrals of g-r-independent forms on the base change of C to .
3.2. Computing Coleman Integrations
In order to Compute , we need to evaluate for where is the g-dimensional vector space of regular 1-forms on C. Let be a differential basis of with , then we can write
where and is in the Monsky-Washnitzer (MW)
weak completion of the coordinate ring of C deprived of the Weierstrass points . And let . So:
Suppose that p is a prime of good reduction for C and let be the reduction of C modulo p. Then there is a natural reduction which sends the points of to the points of . If , denote by the reduction of P modulo p. We say that two points are in the same residue disk if they have the same reduction modulo p. To compute , we consider two cases: P and Q are in the same residue disk and P and Q are in the residue disks different. Algorithms in both cases are given below (Algorithm 1, Algorithm 2). For more details see .
4. The Algorithm
In this section, we specialize in the case that interests us, where C is a hyperelliptic curve of genus 3 given by a model of even degree defined by:
where is monic of degree 8 (to avoid the leading coefficient of the polynomial f is not a square in ) and having at least one root in . We assume that the Jacobian J of C has a Mordell-Weil rank 0 over . Our goal is to compute . Since our hyperelliptic curve is given by a polynomial having at least one root in . Then from Section 1, we can find at least one hyperelliptic curve of odd degree model isomorphic to C. We have implemented the algorithm bellow on Sage (Algorithm 3), which determines, the curve isomorphic to C.
Algorithm 1. are in the same residue disk (tiny Integrals).
Algorithm 2. are in different residue disks, we use Frobenius to .
Algorithm 3. Transformation of a genus 3 real hyperelliptic curve in a genus 3 imaginary hyperelliptic curve.
The hyperelliptic curve obtained in this algorithm is given by an equation of the form:
where not necessarily monic of degree 7.
Such a model has a unique point at inifinity which we denote by . The jacobian of also has a rank of Mordell-Weil 0 on .
Let a prime number such that has a good reduction mod p (donte that p must not be divisible by the dominant coefficient of F). We now want, probably to count rational points of so that we can count the set of rational points of C. For this, we use the algorithm described in  which we transcribe here (Algorithm 4).
Note that the output of this algorithm is a finite subset containing , which is returned as three separate subsets:
● The set of rational points of .
● The set of points Q in such that is 2-torsion.
● The set of points Q in such that is an n-torsion point for.
In this paper we are interested that at the set of rational points of .
Note that the implementation of this algorithm only took into account the hyperelliptic curves defined by a polynomial monic . Since, in most of the cases that have processed the hyperelliptic curve it was defined by a polynomial no monic. So we had to add to the code in Sage of  a function allowing this algorithm to take into account polynomial not necessarily monic.
Here we present the code for this function in SageMath Software .
We give below the steps of Algorithm 4, for more details see .
Step 1 (Required precision.) We need to choose the p-adic precision N and the t-adic precision M to guarantee that, in Step 3, we will obtain all the roots of in . Set and is sufficient for p prime .
Algorithm 4. Chabauty-Coleman method for genus 3 rank 0 hypperelliptic curve.
Step 2 (Annihilator) A basis of the space of differentials ; is given by where . For each , define:
where denotes the point at infinity. The functions are zero on all rational points of , but not identically zero.
Step 3 (Searching in residue discs.) For each point , we compute the set of -rational points P reducing to such that for . To perform this computation, we consider two different cases:
1) If there is a point reducing to , then we expand each in terms of a uniformizer t at P and we formally integrate to obtain three power series , that parametrize the integrals of the between P and any other point in the residue disc.
2) Otherwise, we start by finding a -point P reducing to (note that cannot be in this case). If we can take where is the Hensel lift of to a root of . Otherwise, we can take where is any lift of to and is obtained from by applying Hensels Lemma to . Then we set ; where each parametrizes the integral of between P and any other point in the residue disc.
Step 4 (Identifying the rational points.) Now, for each of the points Q found in Step 3, we attempt to reconstruct Q as a Q-rational point, using Sage. If this is not possible, then must be a torsion point in , because J has rank 0.
Once we have the rational points of we can reconstruct the rational points of C. Note that C has the same number of rational points as and being a curve given by a model of even degree, it has two points at infinity that we note and .
We use above algorithms to determine the rational points of C from the rational points of (Algorithm 5).
Note also that, if after reconstruction, we do not have the at inifinity, we must
Algorithm 5. Computation of the rational points of C.
add them because they are known in advance in our case these are the coordinate points and .
We run our algorithms in Sage on a list of 47 hyperelliptic curves of genus 3 and rank 0, obtained from the database  giving the models of even degree polynomials having at least one root in .
Our implementation proves that for each of the studied curves, the set of rational points is equal to the set of rational points of naive height at most 105.
Figure 1 shows how many of the curves in our list have certain number of rational points.
We observe that the maximun of the points is 4 and that a majority of the curves have three rational points.
We conclude this section, showing how the algorithms work on a particular curve and we show an example of a curve which does not have rational points at infinity.
Consider a hyperelliptic curve C of genus 3 given by:
The function RankBound in Magma shows that the jacobian J of C has a rank of Mordell-Weil 0. Then we can apple the algorithms described above to compute the number of rational points of C. We proceed as follows:
Let’s pose . It is easy see that f vanishes into . Therefore, the curve C is -isomorphic a imaginary hyperelliptic curve that denote .
Step 1: Compute of the equation of C'
Let the equation of the curve . By Algorithm 3, we have with , then is given by:
Figure 1. Number of hyperelliptic curves of genus 3 given by a model of even degree in the list and n points rationels.
Step 2 Computing rational points of C'
Using Magma we find that the set of rational points of with a height bounded by 105 is:
Since the curve has good reduction modulo 7 > 6, we run the Algorithm 4 using this prime. The points of are as follows:
After Hensel lifting each of these points to a point of , we write and in local coordinates, we obtained:
Then we use the PARI/GP function polrootspadic to compute the common zeros of the , and . We notice, that they have the common zeros in the discs and which correspond respectively to the points and .
Therefore, we have shown that .
Step 3 Determination of the rational points of C
We now move on to the phase of constructing the rational points of C from using Algorithm 5. We start adding the point used for the variable change i.e.: and the points at infinity and . Let , we compute the point such that
and Hence . Since Q is the only point in
whose the x-coordinate is different from , Algorithm 5 stops here. Therefore:
We check that we have the same result using Magma with a height of 105.
In our list, we found case of a hyperelliptic curve whose polynomial f is not monic and its leading coefficient is not a square in . So the curve has no pints at infinity. We present it in the following example.
Let be the hyperelliptic curve:
This curve does not have rational points at infinity since the leading coefficient of the polynomial defining the curve is not a square in .
By Algorithm 3, the curve C is transformed into the imaginary hyperelliptic curve given by:
Using Magma, we have:
Since the curve has good reduction modulo 7, we run Algorithm 4 and we find that .
Since has only one rational point, then C will only have one rational point too which will be none other than the point used for the change of variable. By Algorithm 5, we get .
In this paper, we have presented the computation of rational points on hyperelliptic curves of genus 3 given by an even degree model and whose Jacobian has a Mordell-Weil rank 0, using the isomorphism between the real hyperelliptic curve and the imaginary hyperelliptic curve. Our calculations show that the rational points on C are the same as when we search on Magma with a height of at most 105.
 Balakrishnan, J.S., Bianchi, F., Cantoral-Farfán, V., Çiperiani, M. and Etropolski, A. (2018) Chabauty-Coleman Experiments for Genus 3 Hyperelliptic Curves. In: Balakrishnan, J., Folsom, A., Lalín, M. and Manes, M., Eds., Research Directions in Number Theory, Springer, Cham, 67-90.
 Balakrishnan, J.S., Bradshaw, R.W. and Kedlaya, K.S. (2010) Explicit Coleman Integration for Hyperelliptic Curves. In: Hanrot, G., Morain, F. and Thomé, E., Eds., Algorithmic Number Theory, ANTS 2010. Lecture Notes in Computer Science. Springer, Berlin, 16-31.