## ON THE COMPUTATION OF THE KUMMER RATIO OF THE CLASS NUMBER FOR PRIME CYCLOTOMIC FIELDS glyph[star]

ALESSANDRO LANGUASCO†, PIETER MOREE , SUMAIA SAAD EDDIN‡, AND ALISA SEDUNOVA ∗ ∗

ABSTRACT. Let z q be a primitive q th root of unity with q an arbitrary odd prime. The ratio of Kummer's first factor of the class number of the cyclotomic number field Q ( z q ) and its expected order of magnitude (a simple function of q ) is called the Kummer ratio and denoted by r ( q ) . It is known that typically r ( q ) is close to 1, but nevertheless it is believed that it is unbounded, but only large on a very thin sequence of primes q .

We propose an algorithm to compute r ( q ) requiring the evaluation of O ( q log q ) products and O ( q ) logarithms. Using it we obtain a new record maximum for r ( q ) , namely r ( 6766811 ) = 1 709379 . . . . (the old record being r ( 5231 ) = 1 556562 . . . . ). The program used and the results described here, are collected at the following address http://www.math.unipd.it/~languasc/ rq-comput.html .

glyph[star] : This is a (preliminary) report about the computational part of a joint project with Pieter Moree, Sumaia Saad Eddin and Alisa Sedunova .

## 1. INTRODUCTION

Let z q be a primitive q th root of unity with q an arbitrary odd prime. Put

$$G ( q ) \coloneqq 2 q \left ( \frac { q } { 4 \pi ^ { 2 } } \right ) ^ { \frac { q - 1 } { 4 } }.$$

The function G q ( ) is of super-exponential growth in q ; for example G ( 439 ) ≈ 10 117 , G ( 3331 ) ≈ 10 1607 , G ( 9689 ) ≈ 10 5792 and G ( 2918643191 ) ≈ 10 5741570411 .

Kummer [7] proved that h 1 ( q ) , the first factor of the class number of the cyclotomic number field Q ( z q ) , is a positive integer and conjectured that h 1 ( q ) ∼ G q ( ) as q → + ¥ . We define the Kummer ratio as

$$r ( q ) \coloneqq \frac { h _ { 1 } ( q ) } { G ( q ) }.$$

## 2. FIRST METHOD: USING THE DIGAMMA FUNCTION

Recall that q is an odd prime and let c be a primitive odd Dirichlet character mod q . Using Hasse's theorem [6] we have that

$$h _ { 1 } ( q ) = G ( q ) \prod _ { \chi \text{odd} } L ( 1, \chi )$$

Date : October 22, 2019, 01:45.

2010 Mathematics Subject Classification. Primary 11-04; secondary 11Y60.

Key words and phrases. Kummer ratio, class number, cyclotomic fields.

†: corresponding author ; Università di Padova, Dipartimento di Matematica, 'Tullio Levi-Civita', Via Trieste 63, 35121 Padova, Italy. email : alessandro.languasco@unipd.it.

‡: Institute of Financial Mathematics and Applied Number Theory, Johannes Kepler University, Altenbergerstrasse 69, 4040 Linz, Austria. email : Sumaia.Saad\_Eddin@jku.at.

∗ : Max-Planck-Institut für Mathematik, Vivatsgasse 7, D-53111 Bonn, Germany. email : moree@mpimbonn.mpg.de, alisa.sedunova@phystech.edu.

and, by (1), it follows that

Recalling eq. (3.1) of [2], i.e. ,

$$L ( 1, \chi ) = - \frac { 1 } { q } \sum _ { a = 1 } ^ { q - 1 } \chi ( a ) \ \psi ( \frac { a } { q } ),$$

where y ( x ) = ( G ′ / G )( x ) is the digamma function, inserting (3) into (2) we can also write

$$r ( q ) = \left ( \frac { - 1 } { q } \right ) ^ { \frac { q - 1 } { 2 } } \prod _ { \chi \text{odd} } \sum _ { a = 1 } ^ { q - 1 } \chi ( a ) \ \psi ( \frac { a } { q } ).$$

Computationally it is more convenient to work with log r ( q ) rather than r ( q ) , which leads to

$$\log r ( q ) = \frac { q - 1 } { 2 } \left ( i \pi - \log q \right ) + \sum _ { \chi \text{odd} } \log \left ( \sum _ { a = 1 } ^ { q - 1 } \chi ( a ) \ \psi ( \frac { a } { q } ) \right ),$$

in which the last logarithm is a complex one. It is clear that the sum over odd Dirichlet characters in (5) has an imaginary part equal to -p ( q -1 ) / 2; hence

$$\log r ( q ) = - \frac { ( q - 1 ) } { 2 } \log q + \sum _ { \chi \text{odd} } \log \left | \sum _ { a = 1 } ^ { q - 1 } \chi ( a ) \ \psi ( \frac { a } { q } ) \right |.$$

Since in this formula only odd Dirichlet characters appear, we can embed a decimation in frequency strategy 1 in the Fast Fourier Transform (FFT) algorithm to perform the sum over a , see section 4, thus replacing the digamma function with the cotangent one.

Recalling that FFT on a vector of length N costs O ( N log N ) arithmetic operations, we get that computing r ( q ) via (6) has a computational cost of O ( q log q ) products plus the cost of computing ( q -1 ) / 2 values of the cotangent and logarithm functions. This is good but, at least from a theoretical point of view, we can do slightly better, see the next section, even if the performances in the practical computations of Table 2 below are very similar.

## 3. SECOND METHOD: USING A GENERALIZED BERNOULLI NUMBER

Let c be a primitive odd Dirichlet character mod q with q an odd prime. We define the first c -Bernoulli number B 1 , c (see Proposition 9.5.12 of Cohen [1]) as

$$B _ { 1, \chi } \coloneqq \frac { 1 } { q } \sum _ { a = 1 } ^ { q - 1 } a \chi ( a ).$$

By eq. (2.1) of Shokrollahi [12] we have

$$h _ { 1 } ( q ) = 2 q \prod _ { \chi \text{odd} } \left ( - \frac { B _ { 1, \chi } } { 2 } \right ).$$

Inserting (7) and (8) into (1), we obtain

$$r ( q ) = \left ( \frac { q } { 4 \pi ^ { 2 } } \right ) ^ { - \frac { q - 1 } { 4 } } \prod _ { \chi \text{odd} } \left ( - \frac { B _ { 1, \chi } } { 2 } \right ) = \left ( - \frac { \pi } { q ^ { \frac { 3 } { 2 } } } \right ) ^ { \frac { q - 1 } { 2 } } \prod _ { \chi \text{odd} } \sum _ { a = 1 } ^ { q - 1 } a \chi ( a ),$$

1 We use here this nomenclature since it is standard in the literature on the Fast Fourier Transform, but clearly it can be translated into number theoretic language using suitable properties of Dirichlet characters.

$$r ( q ) = \prod _ { \chi \text{odd} } L ( 1, \chi ).$$

which leads to

$$\log r ( q ) = \frac { q - 1 } { 2 } \left ( \log \pi - \frac { 3 } { 2 } \log q + i \pi \right ) + \sum _ { \chi \text{odd} } \log \left ( \sum _ { a = 1 } ^ { q - 1 } a \chi ( a ) \right ),$$

in which the last logarithm is a complex one. Moreover it is clear that the sum over the odd Dirichlet characters in (10) has an imaginary part equal to -p ( q -1 ) / 2; hence

$$\log r ( q ) = \frac { q - 1 } { 2 } \left ( \log \pi - \frac { 3 } { 2 } \log q \right ) + \sum _ { \chi \text{odd} } \log \left | \sum _ { a = 1 } ^ { q - 1 } a \chi ( a ) \right |.$$

This formula for log r ( q ) , unlike (6), does not require the computation of values of a special function: it is enough to use the sequence a = 1 , . . . , q -1. Since only the sum over odd Dirichlet characters is needed, we can embed a decimation in frequency strategy in the Fast Fourier Transform (FFT) algorithm to perform the sum over a , see section 4. As in the previous section, it is easy to see that computing r ( q ) via (11) has a computational cost of O ( q log q ) arithmetic operations plus the cost of computing ( q -1 ) / 2 values of the logarithm function and products: so far, this is the fastest known algorithm to compute r ( q ) .

This way we were able to get a new record maximal value for r ( q ) , namely

$$r ( 6 7 6 6 8 1 1 ) = 1. 7 0 9 3 7 9 0 4 2 \dots,$$

see Table 3; such a result was also double-checked using the method of section 2. The previously known record r ( 5231 ) = 1 556562 . . . . is due to Shokrollahi [12]. We will see more on these computations in section 4.

## 4. COMPARING METHODS, RESULTS AND RUNNING TIMES

First of all we notice that PARI/Gp, v. 2.11.2, has the ability to generate Dirichlet L -functions (and many other L -functions) and hence, using (2), the value of log r ( q ) = GLYPH&lt;229&gt; c odd log | L ( 1 , c ) | can be obtained with few instructions of the gp scripting language. This computation has a linear cost in the number of calls of the lfun function of PARI/Gp and it is, at least on our Dell Optiplex desktop machine, slower than both the approaches we are about to describe below.

The other approaches we can use to compute log r ( q ) are the following:

- a) use formula (5) and the y -values;
- b) use formula (10) and the first c -Bernoulli number.

This way we can double check the computation we will perform. In both (5) and (10) we remark that, since q is prime, it is enough to determine a primitive root g of q , which leads to the Dirichlet character c 1 mod q uniquely determined by c 1 ( g ) = e 2 p i / ( q -1 ) . The set of non-trivial characters mod q is then { c j 1 : j = 1 2 , , . . . , q - } 2 . Hence, if, for every k ∈ { 0 , . . . , q - } 2 , we denote g k ≡ ak ∈ { 1 , . . . , q - } 1 , every summation in (5) and (10) is of the type GLYPH&lt;229&gt; q -2 k = 0 e 2 p i jk / ( q -1 ) f ( ak / q ) , where j ∈ { 1 , . . . , q - } 2 is odd and f is a suitable function. As a consequence, such quantities are the Discrete Fourier Transform (DFT) of the sequence { f ( ak / q ) : k = 0 , . . . , q - } 2 . This observation is due to Rader [11] and it was used in [2] and in [8] to speed up the computation of similar quantities via the use of Fast Fourier Transform dedicated software libraries.

In our case we can also use the decimation in frequency strategy: following the line in section 4.1 of [8], letting e x ( ) : = exp 2 ( p ix ) , m = ( q -1 ) / 2, for every j = 0 , . . . , q -2, j = 2 t + glyph[lscript] ,

glyph[lscript] ∈ { 0 1 , } and t ∈ Z , we have that

$$\ell \in \{ 0, 1 \} & \text{ and } t \in \mathbb { Z }, \, \text{ we have that} \\ \sum _ { k = 0 } ^ { q - 2 } e \Big ( \frac { j k } { q - 1 } \Big ) f \Big ( \frac { a _ { k } } { q } \Big ) & = \sum _ { k = 0 } ^ { m - 1 } e \Big ( \frac { t k } { m } \Big ) e \Big ( \frac { \ell k } { q - 1 } \Big ) \Big ( f \Big ( \frac { a _ { k } } { q } \Big ) + ( - 1 ) ^ { \ell } f \Big ( \frac { a _ { k + m } } { q } \Big ) \Big ) \\ & = \begin{cases} \sum _ { k = 0 } ^ { m - 1 } e \Big ( \frac { t k } { m } \Big ) b _ { k } & \text{ if } \ell = 0 ; \\ \sum _ { k = 0 } ^ { m - 1 } e \Big ( \frac { t k } { m } \Big ) c _ { k } & \text{ if } \ell = 1, \end{cases}$$

where t = 0 , . . . , m -1,

$$b _ { k } \coloneqq f \left ( \frac { a _ { k } } { q } \right ) + f \left ( \frac { a _ { k + m } } { q } \right ) \quad \text{and} \quad c _ { k } \coloneqq e \left ( \frac { k } { q - 1 } \right ) \left ( f \left ( \frac { a _ { k } } { q } \right ) - f \left ( \frac { a _ { k + m } } { q } \right ) \right ).$$

Since we just need the sum over the odd Dirichlet characters for f ( x ) = x of f ( x ) = y ( x ) , instead of computing an FFT transform of length q -1 we can evaluate an FFT of half a length, applied on a suitably modified sequence according to (12). Clearly this leads to a gain in speed and a reduction in memory occupation in running the actual computer program. In case f ( x ) = y ( x ) we can simplify the expression e k ( / ( q -1 )) ( y ( ak / q ) -y ( ak + m / q ) ) for ck , where m = ( q -1 ) / 2 and k = 0 , . . . , m -1, in the following way. Recalling that 〈 g 〉 = Z ∗ q , ak ≡ g k mod q and g m ≡ q -1 mod q , we can write

$$\psi \Big ( \frac { a _ { k + m } } { q } \Big ) = \psi \Big ( \frac { q - a _ { k } } { q } \Big ) = \psi \Big ( 1 - \frac { a _ { k } } { q } \Big )$$

and hence, using the well-known reflection formula y ( 1 -x ) -y ( x ) = p cot ( p x ) , we obtain

$$\psi \Big ( \frac { a _ { k } } { q } \Big ) - \psi \Big ( 1 - \frac { a _ { k } } { q } \Big ) = - \pi \cot \Big ( \frac { \pi a _ { k } } { q } \Big ),$$

for every k = 0 , . . . , m -1. Inserting the last relation in the definition of ck in (12) we can replace in the actual computation the digamma function with the cotangent one. The case f ( x ) = x is easier; using again 〈 g 〉 = Z ∗ q , ak ≡ g k mod q and g m ≡ q -1 mod q , we can write that ak + m ≡ q -ak mod q ; hence

$$a _ { k } - a _ { k + m } = a _ { k } - ( q - a _ { k } ) = 2 a _ { k } - q$$

so that in this case we obtain ck = ( e k / ( q -1 ))( 2 ak / q -1 ) for every k = 0 , . . . m -1, m = ( q -1 ) / 2.

- 4.1. Computations with trivial summing over a (slower, more decimal digits available). Unfortunately in libpari the FFT-functions work only if q = 2 glyph[lscript] + 1, for some glyph[lscript] ∈ N . So we had to trivially perform these summations and hence, in practice, this part is the most time consuming one in both approach a) and b), as the cost is quadratic in q .

Being aware of such limitations, we performed the computation of r ( q ) with these three approaches for every q prime, 3 ≤ q ≤ 1000, on a Dell OptiPlex-3050, equipped with an Intel i5-7500 processor, 3.40GHz, 16 GB of RAM and running Ubuntu 18.04.2, using a precision of 30 decimal digits, see Table 1. The results coincide up the desired precision. We also computed the values of r ( q ) for q = 1451, 2741, 3331, 4349, 4391, 5231, 6101, 6379, 7219, 8209, 9049, 9689 and these are given in Table 2. These numbers were chosen to extend the number of available decimals for the known data (see Fung-Granville-Williams [4] and Shokrollahi [12]). In this case, in the fifth column of Table 2 we also reported the running time of the direct approach, i.e. using (2), the third and fourth columns are respectively the running times of the approaches a) and b). For these values of q it became clear that the computation time spent in

performing the sums over a was the longest one. This means that inserting an FFT-algorithm in the approaches a) and b) is fundamental to further improve their performances. We will say more on this in the next section.

## 4.2. Computations summing over a via FFT (much faster, less decimal digits available).

- As we saw before, as q becomes large, the time spent in summing over a dominates the overall computational cost. So we implemented the use of the FFT in both the approaches a) and b), by using the fftw [3] library in our C programs. The performance of this part was extremely good in the sense that it was a factor 1000 faster than the same one trivially performed.
- 4.3. Data for the scatter plots. We were able to compute the long double precision r ( q ) -values for every prime 3 ≤ q ≤ 2 · 10 6 and we provide here the scatter plot, see Figure 1, of such values. The minimal value is r ( 3 ) = 0 6045997880 . . . . and the maximal one is r ( 305741 ) = 1 661436 . . . . The data were obtained in about three days of computation time on the Dell OptiPlex machine mentioned before.
- 4.4. Computations for larger q . If bq + 1 is prime for many small b , there is a good chance that r ( q ) will be large. Promising, using this criterion, seemed q = 4178771, 6766811, 28227761, 193894451, 75743411, 212634221, 251160191, 405386081, 538906601, 964477901 and also 1139803271, 1217434451, 1806830951, 2488788101, 2830676081, 2918643191 . For these q the r ( q ) 's were evaluated using the quadruple precision, see Table 3. We remark that the quadruple precision computation performances are affected by a lack of hardware support of the FLOAT128 type of the C programming language.

In Table 4 we evaluated some further cases with potential large r ( q ) value: q = 4151292581, 6406387241, 7079770931, 9854964401 with the long double precision. These computations were performed on an Intel(R) Xeon(R) CPU E5-2650 v3 @ 2.30GHz, with 160 GB of RAM and running Ubuntu 16.04. The computation for q = 9109334831 was performed on the CAPRI ('Calcolo ad Alte Prestazioni per la Ricerca e l'Innovazione') infrastructure of the University of Padova.

The PARI/Gp scripts and the C programs used and the computational results obtained are available at the following web address: http://www.math.unipd.it/~languasc/rq-comput. html .

## 5. FURTHER COMPUTATION ON THE EULER-KRONECKER CONSTANTS

The Euler-Kronecker constant for the prime cyclotomic field Q ( z q ) and for Q ( z q + z -1 q ) , the maximal real subfield of Q ( z q ) , see, e.g. , [8], are defined as glyph[negationslash]

$$\mathfrak { G } _ { q } \coloneqq \gamma + \sum _ { \chi \neq \chi _ { 0 } } \frac { L ^ { \prime } ( 1, \chi ) } { L ( 1, \chi ) }, \quad \text{ respectively} \quad \mathfrak { G } _ { q } ^ { + } \coloneqq \gamma + \sum _ { \overline { \chi } \neq \chi _ { 0 } } \frac { L ^ { \prime } ( 1, \chi ) } { L ( 1, \chi ) }. \quad \text{ (13)}$$

glyph[negationslash]

Using formula (22) of [8] we have that

$$\mathfrak { G } _ { q } - \mathfrak { G } _ { q } ^ { + } = \sum _ { \chi \text{odd} } \frac { L ^ { \prime } ( 1, \chi ) } { L ( 1, \chi ) } = \frac { q - 1 } { 2 } \left ( \gamma + \log ( 2 \pi ) \right ) + \sum _ { \chi \text{odd} } \frac { 1 } { B _ { 1, \overline { \chi } } } \sum _ { a = 1 } ^ { q - 1 } \overline { \chi } ( a ) \log \left ( \Gamma ( \frac { a } { q } ) \right ),$$

which can be easily implemented since the function log G is available in the C programming language. The scatter plot of Figure 2 represents the normalized values of G q -G + q , q prime, 3 ≤ q ≤ 2 · 10 6 and creating it took about three days on the Dell Optiplex machine mentioned before.

Since in (14) we can use a decimation in frequency strategy, we also remark that letting f = log G into (12) leads to simplify the form of ck = ( e -k / ( q -1 )) ( log G ( ak / q ) -log G ( ak + m / q ) ) , where m =( q -1 ) / 2 and k = 0 , . . . , m -1, in the following way 2 . Recalling 〈 g 〉 = Z ∗ q , ak ≡ g k mod q and g m ≡ q -1 mod q , we can write that

$$\log \left ( \Gamma ( \frac { a _ { k + m } } { q } ) \right ) = \log \left ( \Gamma ( \frac { q - a _ { k } } { q } ) \right ) = \log \left ( \Gamma ( 1 - \frac { a _ { k } } { q } ) \right )$$

and hence, using the well-known reflection formula G ( x ) G ( 1 -x ) = p / sin ( p x ) , we obtain

$$\log \left ( \Gamma ( \frac { \mu } { q } ) \right ) - \log \left ( \Gamma ( 1 - \frac { a _ { k } } { q } ) \right ) = 2 \log \left ( \Gamma ( \frac { \mu _ { k } } { q } ) \right ) + \log \left ( \sin ( \frac { \pi a _ { k } } { q } ) \right ) - \log \pi,$$

for every k = 0 , . . . , m -1, m =( q -1 ) / 2. Inserting the last relation in the definition of ck in (12) we obtain the actual sequence we used for performing the computation previously mentioned.

Acknowledgements . Part of the calculations described here were performed using the University of Padova Strategic Research Infrastructure Grant 2017: 'CAPRI: Calcolo ad Alte Prestazioni per la Ricerca e l'Innovazione'. The first author, A. Languasco, thanks Luca Righi (University of Padova) for his help in using CAPRI. The second and fourth authors, P. Moree and A. Sedunova, thank the Max Planck Institute for Mathematics for providing excellent conditions while they were working on the project this paper is part of. The third author, S. Saad Eddin, is supported by the Austrian Science Fund (FWF): Project F5505-N26 and Project F5507-N26, which are part of the special Research Program 'Quasi Monte Carlo Methods: Theory and Application'.

## REFERENCES

- [1] H. Cohen, Number Theory. Volume II: Analytic and Modern Tools , Graduate Texts in Mathematics, vol. 240, Springer, 2007.
- [2] K. Ford, F. Luca, P. Moree, Values of the Euler f -function not divisible by a given odd prime, and the distribution of Euler-Kronecker constants for cyclotomic fields , Math. Comp. 83 (2014), 1447-1476.
- [3] M. Frigo, S. G. Johnson, The Design and Implementation of FFTW3 , Proceedings of the IEEE 93 (2), 216-231 (2005). The C library is available at http://www.fftw.org .
- [4] G. Fung, A. Granville, H. C. Williams, Computation of the first factor of the class number of cyclotomic fields , J. Number Theory 42 (1992), 297-312.
- [5] Gnu Scientific Library, version 2.5 , 2018. Available from http://www.gnu.org/software/gsl/ .
- [6] H. Hasse, Über die Klassenzahl abelscher Zahlkörper , Akademie-Verlag, Berlin, 1952, reprinted with an introduction by J. Martinet, Springer-Verlag, 1985.
- [7] E. E. Kummer, Memoire sur la théorie des nombres complexes composes de racines de l'unité et des nombres entiers , J. Math. Pures Appl. (1851), 377-498, reprinted in his Collected papers, Vol. I, Springer-Verlag, 1975, pp. 363-484.
- [8] A. Languasco, A note on the computation of the Euler-Kronecker constants for prime cyclotomic fields , Arxiv, 2019, http://arxiv.org/abs/1903.05487 .
- [9] P. Moree, Irregular Behaviour of Class Numbers and Euler-Kronecker Constants of Cyclotomic Fields: The Log Log Log Devil at Play , Irregularities in the Distribution of Prime Numbers. From the Era of Helmut Maier's Matrix Method and Beyond (J. Pintz and M.Th. Rassias, eds.), Springer, 2018, pp. 143-163.
- [10] The PARI Group, PARI/GP version 2.11.2 , Bordeaux, 2019. Available from http://pari.math. u-bordeaux.fr/ .
- [11] C. M. Rader, Discrete Fourier transforms when the number of data samples is prime , Proc. IEEE 56 (1968), 1107-1108.
- [12] M. A. Shokrollahi, Relative class number of imaginary abelian fields of prime conductor below 10000, Math. Comp. 68 (1999), 1717-1728.

2 The minus sign here present in the twiddle factor e ( -k / ( q -1 )) comparing with the one in (12) depends on the fact that we are now summing over the conjugate Dirichlet character c instead over c .

## KUMMER RATIO FOR Q ( z q )

TABLE 1. Values of r ( q ) for every odd prime up to 1000 with a precision of 30 digits; computed with PARI/Gp, v. 2.11.2, with a trivial way of executing the sum over a . Total computation time: 8s., 947ms. using the first c -Bernoulli number. [s=seconds; msec=milliseconds]

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01152v4-figure-5f88021a-8032-4158-9a50-b968e01426cb.jpg)

| q       | r ( q )                                                                         | q       | r ( q )                                                                         | q       | r ( q )                                                                           |
|---------|---------------------------------------------------------------------------------|---------|---------------------------------------------------------------------------------|---------|-----------------------------------------------------------------------------------|
| 3       | 0 . 60459978807807261686469275254 . . .                                         | 271     | 0 84120880901441103034587178906 . . .                                           | 619     | . 80463918636548231818097049238 . .                                               |
| 5       | 0 . 78956835208714868950675927999 . . .                                         | 277     | . 1 . 22287167700803659996325347046 . . .                                       | 631     | 0 . 1 . 13964698072442766479584447731 . . .                                       |
| 7       | 0 . 95667518575084187547950733813 . . .                                         | 281     | 1 . 09072312671446411507457756821 . . .                                         | 641     | 1 . 34299156432328475445263673245 . . .                                           |
| 11      | 1 . 10916191287000575896982175317 . . .                                         | 283     | 0 . 98730045924989351176735192970 . . .                                         | 643     | 1 . 01836205611360685304417553494 . . .                                           |
| 13      | 1 . 07714905620985756748597815892 . . .                                         | 293     | 1 . 28843023595237283191056798455 . . .                                         | 647     | 0 . 90233667317118875595490772209 . . .                                           |
| 17      | 0 . 85539034568765268115905873936 . . .                                         | 307     | 0 . 91358725220199482220514916899 . . .                                         | 653     | 1 . 27087727805772466468796098338 . . .                                           |
| 19      | 0 . 70704004900384729070674621978 . . .                                         | 311     | 1 . 14589374542647302213444142687 . . .                                         | 659     | 1 . 39106317898226550143998268526 . . .                                           |
| 23      | 1 . 27303069939685502234405162961 . . .                                         | 313     | 0 . 93893317675819166180673984422 . . .                                         | 661     | 0 . 83544430975232146569368385972 . . .                                           |
| 29      | 1 . 19507225854723141702138692301 . . .                                         | 317     | 0 . 80671823188984812849457198577 . . .                                         | 673     | 1 . 03660206982398637181187353211 . . .                                           |
| 31      | 0 . 88988962107854407891985181571 . . .                                         | 331     | 0 . 81356274956051845902331649336 . . .                                         | 677     | 0 . 92424013312497364401792044662 . . .                                           |
| 37      | 0 . 89617354245182624263930105684 . . .                                         | 337     | 0 . 86111511521922591268832255791 . . .                                         | 683     | 1 . 13528281402409476998254691134 . . .                                           |
| 41      | 1 . 01095149281551337376703651618 . .                                           | 347     | 1 . 08517941758105267446483313058 . . .                                         | 691     | 0 . 76921427957454050696406411036 . . .                                           |
| 43      | . 1 . 00032807083987921579084335194 . . .                                       | 349     | 0 . 98395731344877010445591239132 . . .                                         | 701     | 0 . 92089882867969861041626254382 . . .                                           |
| 47      | 0 . 99510419475843763320461794597 . . .                                         | 353     | 0 . 88603505661744604503087815775 . . .                                         | 709     | 1 . 05648934917801861606174800341 . . .                                           |
| 53      | 1 . 00231549556080469808835403497 . . .                                         | 359     | 1 . 16002644446708254566916432735 . . .                                         | 719     | 1 . 20306325855333927681117243729 . . .                                           |
| 59      | 1 . 03111995957758588341749868917 . .                                           | 367     | 0 . 90864101877936912063265319825 . . .                                         | 727     | 0 . 99856921422780328631340639607 . . .                                           |
| 61      | . 0 . 91541689757636152038607840584 . . .                                       | 373     | 1 . 07507614420133257646267035535 . . .                                         | 733     | 0 . 98014910177261986736078022621 . . .                                           |
| 67      | 1 . 03230196304201968151553976333 . .                                           | 379     | 0 . 72144618647138444694426995688 . . .                                         | 739     | 1 . 10263546824053086630671245464 . . .                                           |
| 71      | . 0 . 94652474710362368092900546271 . . .                                       | 383     | 0 . 83243809267428713960470760380 . . .                                         | 743     | 1 . 03495494096205775904091176831 . . .                                           |
| 73      | 1 . 28217793230760538382246761185 . .                                           | 389     | 0 . 84997782896854503971627563496 . . .                                         | 751     | 1 . 01856200583585073878095848975 . . .                                           |
| 79      | . 0 . 84579459612002975504552940763 . . .                                       | 397     | 0 . 99757781120158579094253246616 . . .                                         | 757     | 0 . 96706876118708598545541455447 . . .                                           |
| 83      | 1 . 22326926548441461619507621390 . . .                                         | 401     | 1 . 13998328316447070631384278931 . . .                                         | 761     | 1 . 46958285813141552491322656984 . . .                                           |
| 89      | 1 . 28632147461922346234453694590 . . .                                         | 409     | 1 . 19919809743909540748744244798 . . .                                         | 769     | 0 . 89892230367392111314972716474 . . .                                           |
| 97      | 0 . 90467614287023765066781857933 . . .                                         | 419     | 1 . 18974458882376935926766971002 . . .                                         | 773     | 1 . 06810947197031447130333305035 . . .                                           |
| 101     | 1 . 11049958753586448051923888082 . . .                                         | 421     | 0 . 86457966530711741177342869535 . . .                                         | 787     | 0 . 97178232843986336686451556474 . . .                                           |
| 103     | 1 . 05565198833718743186163713482 . . .                                         | 431     | 1 . 13754261103593462461717085623 . . .                                         | 797     | 1 . 03075130387360942943641986527 . . .                                           |
| 107     | 0 . 99260767792672501309519612375 . . .                                         | 433     | 1 . 07176135182041771385450595205 . . .                                         | 809     | 1 . 31970441406018712251949567645 . . .                                           |
| 109     | 0 . 91554283885230186850667500246 . . .                                         | 439     | 0 . 68484134061729762055005895626 . . .                                         | 811     | 0 . 80283817264815420707856818905 . . .                                           |
| 113     | 1 . 16185573635061808057761114590                                               | 443     | 1 . 41089988430397986980906568345 . . .                                         | 821     | 1 . 06528437036549643319352814656 . . .                                           |
| 127     | . . . 1 . 06269835499717635407980190888 . . .                                   | 449     | 0 . 90539643658614424895891547460 . . .                                         | 823     | 0 . 96769318476182048656465705918 . . .                                           |
| 131     | 1 . 27897699389762867270592988247 . . .                                         | 457     | 0 . 83734634190585621778636791343 . . .                                         | 827     | 0 . 86555993675758442057691969953 . . .                                           |
| 137     | 1 . 00188853650420792851571142833 . . .                                         | 461     | 1 . 03119557377397403645284724907 . . .                                         | 829     | 0 . 82250033541615549748400919644 . . .                                           |
| 139     | 0 . 87166115187392327886708542130 . . .                                         | 463     | 0 . 96134625111959841778686635231 . . .                                         | 839     | 0 . 91871090540765761610044317666 . . .                                           |
| 149     | 1 . 04886527642691194564791006449 . . .                                         | 467     | 0 . 89740454859192836870657083737 . . .                                         | 853     | 1 . 08223582880253347548004283614 . . .                                           |
| 151     | 1 . 09613526050530812035603232922 . . .                                         | 479     | 1 . 10506715780642069705910978939 . . .                                         | 857     | 1 . 05075311490694694576396382027 . . .                                           |
| 157     | 0 . 74304505329108896600523002862 . . .                                         | 487     | 1 . 13041022782656063139453697156 . . .                                         | 859     | 0 . 88080094180568178476397675721 . . .                                           |
| 163     | 0 . 95167392369442992883081838307 . . .                                         | 491     | 1 . 27221465691304968352754354985 . . .                                         | 863     | 1 . 05694231206444764180240401289 . . .                                           |
| 167     | 0 . 85404891714098835186838607451 . . .                                         | 499     | 0 . 82979024959465063669881382680 . . .                                         | 877     | 0 . 72289398522705741218284637854 . . .                                           |
| 173     | 1 . 25750311100604863256476652232 . . .                                         | 503     | 1 . 09956174719578329093362210466 . . .                                         | 881     | 1 . 09738994199075350184435336371 . . .                                           |
| 179     | 1 . 31898955218699008540672120548 . . .                                         |         | 1 . 39692082719612661320417410651 . . .                                         | 883     | . . .                                                                             |
| 181     | 1 . 01646725307901783240856438797 . . .                                         | 509     | 0 . 74488579181918272860910159248 . . .                                         | 887     | 1 13318227639393214982039012685 . 0 . 96917974196790823108417719936 . . .         |
|         |                                                                                 | 521     |                                                                                 | 907     | .                                                                                 |
| 191     | 1 . 29850955347246763676155271715 . . .                                         | 523     | 0 . 99514847873992894203802693220 . . .                                         |         | 0 . 90262558866311480478051623607 . . .                                           |
| 193     | 1 . 17384956614280523683625176108 . . .                                         | 541     | 0 . 94472655782952981521345779529 . . .                                         | 911     | 1 . 07798557536304873099351043701 . . .                                           |
| 197 199 | 0 . 87142685805870225854275086741 . . . . . .                                   | 547     | 0 . 73868505476195458996166613201 . . . 1 . 01800618130970440243478675140 . . . | 919     | 1 . 04003346554199950901317303457 . .                                             |
| 211     | 0 79775765981803261703336410970 . 0 . 70965810384577007739153826881 . . .       | 557 563 | 0 . 92322125091337523644162001846 . . .                                         | 929 937 | 1 . 04414904452989167744813201725 . . . 0 . 90017934857750019784132262523 . . . . |
| 223     | 0 . 90016736774009107389420074861 . . .                                         | 569     | 0 . 86644384514357385272705168284 . . .                                         | 941     | 1 . 09400867179752235523397214847 . .                                             |
|         | 0 . 76298839763137122603762871170 . . .                                         | 571     | 0 . 99662480636851972762309151349 . . .                                         | 947     | 1 . 22587448270510743026091490433 . .                                             |
| 227     |                                                                                 |         |                                                                                 |         | 1 . 16083173031283885682226845601 . .                                             |
| 229     | 0 . 72414574142010494620086404196 . . .                                         | 577     | 0 . 91370293804018510239277389204 . . .                                         | 953     | . .                                                                               |
| 233     | 1 . 43102216731058063469583770264 . . .                                         | 587     | 0 . 81252459850672121660374173954 . . .                                         | 967     | 0 . 72860004404668861481436825047 . . .                                           |
| 239     | 1 . 18520259221018381028526578871 . . .                                         | 593     | 1 . 07734617489664930780759188442 . . .                                         | 971 977 | 1 . 07939115916440046258710381604 . . .                                           |
| 241     | 1 . 11908192699651325481120769078 . . . 1 . 18041694425392859170387583509 . . . | 599 601 | 0 . 96408773834723069779571268471 . . .                                         |         | 0 . 83890885880371282354125472475 . . .                                           |
| 251     |                                                                                 |         | 0 . 92827339751824097250854300550 . . .                                         | 983     | 0 . 78867677202973854047246566763 . . .                                           |
| 257     | 0 . 90559625735496576640913464538 . . .                                         | 607     | 0 . 83637312705251443247667799101 . . .                                         | 991     | 0 . 90943936153505129760069639750 . . .                                           |
| 263     | 0 . 93717078166852960654064932319 . . . .                                       | 613     | 0 . 87703659303472148910355020294 . . .                                         | 997     | 0 . 85575754491350654466545217865 . . .                                           |
| 269     | 1 . 01052429941342866041104883014 . .                                           | 617     | 0 . 84246084541946716141445378848 . . .                                         |         |                                                                                   |

|    q | r ( q )             | time y -version   | time Bernoulli-version   | time direct version   |
|------|---------------------|-------------------|--------------------------|-----------------------|
| 1451 | 1 . 489316072 . . . | 292ms.            | 299ms.                   | 2s. 083ms.            |
| 2741 | 1 . 498121015 . . . | 1s. 019ms.        | 1s. 061ms.               | 5s. 231ms.            |
| 3331 | 0 . 642429297 . . . | 1s. 496ms.        | 1s. 566ms.               | 7s. 236ms.            |
| 4349 | 1 . 518570512 . . . | 2s. 536ms.        | 2s. 657ms.               | 10s. 572ms.           |
| 4391 | 1 . 507776410 . . . | 2s. 586ms.        | 2s. 710ms.               | 11s. 276ms.           |
| 5231 | 1 . 556562248 . . . | 3s. 361ms.        | 3s. 385ms.               | 14s. 730ms.           |
| 6101 | 1 . 511405291 . . . | 4s. 964ms.        | 5s. 086ms.               | 17s. 490ms.           |
| 6379 | 0 . 673523026 . . . | 5s. 436ms.        | 5s. 699ms.               | 19s. 859ms.           |
| 7219 | 0 . 658084090 . . . | 6s. 978ms.        | 7s. 293ms.               | 23s. 966ms.           |
| 8209 | 0 . 672045039 . . . | 8s. 950ms.        | 9s. 416ms.               | 27s. 322ms.           |
| 9049 | 0 . 667614244 . . . | 10s. 870ms.       | 11s. 461ms.              | 32s. 610ms.           |
| 9689 | 1 . 524371504 . . . | 12s. 487ms.       | 13s. 113ms.              | 36s. 855ms.           |

TABLE 2. A few other values of r ( q ) with a precision of 10 digits; computed with PARI/Gp, v. 2.11.2, with a trivial way of executing the sum over a . [s=seconds; msec=milliseconds]

| q                                                                                                                                                               | r ( q )                                                                                                                                                                                                                                                                                                                         | total time (long double)                                                                                                                                                                                                                                 | total time (quadruple)                                                                                                                                                                                                                                                           |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| 4178771 6766811 28227761 75743411 193894451 212634221 251160191 405386081 538906601 964477901 1139803271 1217434451 1806830951 2488788101 2830676081 2918643191 | 1 . 611588128 . . . 1 . 709379041 . . . 1 . 528720351 . . . 1 . 645759517 . . . 1 . 548501406 . . . 1 . 652149469 . . . 1 . 611898472 . . . 1 . 545118923 . . . 1 . 693680145 . . . 1 . 612596619 . . . 1 . 398836497 . . . 1 . 707310115 . . . 1 . 621464926 . . . 1 . 662760638 . . . 1 . 616923086 . . . 1 . 693092857 . . . | 2s. 471ms. 4s. 200ms. 16s. 267ms. 1m. 05s. 078ms. 2m. 18s. 739ms. 2m. 44s. 018ms. 3m. 02s. 783ms. 4m. 31s. 634ms. 6m. 44s. 584ms. 12m. 21s. 793ms. 21m. 56s. 349ms. 15m. 00s. 374ms. 23m. 00s. 240ms. 29m. 54s. 639ms. 32m. 46s. 208ms. 38m. 40s. 423ms. | 42s. 007ms. 1m. 01s. 314ms. 4m. 30s. 804ms. 19m. 04s. 250ms. 38m. 27s. 838ms. 48m. 24s. 277ms. 54m. 20s. 875ms. 58m. 49s. 163ms. 114m. 59s. 330ms. 217m. 40s. 768ms. 298m. 29s. 058ms. 277m. 58s. 390ms. 417m. 46s. 183ms. 512m. 29s. 915ms. 552m. 42s. 941ms. 697m. 10s. 305ms. |

TABLE 3. Few other values of r ( q ) ; computed using the first c -Bernoulli number with PARI/Gp, v. 2.11.2. and fftw, v. 3.3.8., with long double and quadruple precisions. The sum over a was performed using the FFT algorithm on the Xeon machine mentioned before. [m=minutes; s=seconds; msec=milliseconds, n.a.=not available]

TABLE 4. Few other values of r ( q ) ; computed using the first c -Bernoulli number with PARI/Gp, v. 2.11.2. and fftw, v. 3.3.8., with long double precision. The sum over a was performed using the FFT algorithm on the Xeon machine mentioned before with the only exception of q = 9109334831 for which we used the CAPRI infrastructure of the University of Padova. [m=minutes; s=seconds]

|          q | r ( q )          | total time (long double)   |
|------------|------------------|----------------------------|
| 4151292581 | 1 . 669735 . . . | 63m. 00s.                  |
| 6406387241 | 1 . 625741 . . . | 75m. 38s.                  |
| 7079770931 | 1 . 688607 . . . | 187m. 20s.                 |
| 9109334831 | 1 . 657855 . . . | 120m. 28s.                 |
| 9854964401 | 1 . 688033 . . . | 177m. 14s.                 |

FIGURE 1. The values of r ( q ) , q prime, 3 ≤ q ≤ 2 10 6 · , plotted using GNUPLOT, v.5.2, patchlevel 7. [max = 1 661436 . . . . attained at q = 305741; min = 0 604599 . . . . attained at q = 3; number of r ( q ) &gt; 1: 74795 (50 22%); number of . r ( q ) &lt; 1: 74137; (49 78%) . total number of data = 148932].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01152v4-figure-a29cf01a-0ecf-42e8-80f2-ace9c5be7524.jpg)

FIGURE 2. The values of ( G q -G + ) q / log q , q prime, 3 ≤ q ≤ 2 · 10 6 , plotted using GNUPLOT, v.5.2, patchlevel 7. [ max = 0 546473 . . . . attained at q = 37189; min = -0 569200 . . . . attained at q = 305741; number of positive data = 74190 (49 81%); . number of negative data = 74742 (50 19%); total number of data . = 148932].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01152v4-figure-1ae90ecc-7946-4feb-b2e2-8cc0d200b289.jpg)
Alessandro Languasco Pieter Moree Sumaia Saad Eddin Alisa Sedunova
Computation of the Kummer ratio of the class number for prime cyclotomic fields
Let $\zeta_q$ be a primitive $q^{\text{th}}$ root of unity with $q$ an
arbitrary odd prime. The ratio of Kummer's first factor of the class number of
the cyclotomic number field $\mathbb{Q}(\zeta_q)$ and its expected order of
magnitude (a simple function of $q$) is called the Kummer ratio and denoted by
$r(q)$. It is known that typically $r(q)$ is close to 1, but nevertheless it is
believed that it is unbounded, but only large on a very thin sequence of primes
$q$. We propose an algorithm to compute $r(q)$ requiring the evaluation of
$O(q\log q)$ products and $O(q)$ logarithms. Using it we obtain a new record
maximum for $r(q)$, namely $r(6766811) =1.709379\dotsc$ (the old record being
$r(5231)=1.556562\dotsc$). The program used and the results described here, are
collected at the following address
\url{http://www.math.unipd.it/~languasc/rq-comput.html}.
  This is a (preliminary) report about the computational part of a joint
project with Pieter Moree, Sumaia Saad Eddin, and Alisa Sedunova.