Teoría Analítica de Números
Teoría Analítica de Números
TEORÍA ANALÍTICA DE
NÚMEROS
¿Por qué los números son hermosos? Es como
preguntar por qué la novena sinfonía de Beethoven es
hermosa. Si no ves por qué, nadie puede explicártelo.
Yo sé que los números son hermosos. Si no lo son,
nada lo es.
Paul Erdős
Índice General
Introducción vii
v
vi ÍNDICE GENERAL
Bibliografía 283
Este libro supone que el lector tiene algunos conocimientos básicos de álgebra
y análisis matemático (bastantes menos que los contenidos en mis libros de
Álgebra [Al] y Análisis matemático [An]) y buena parte de él puede verse como
una muestra de las aplicaciones de la teoría expuesta en mi libro de Funciones
de variable compleja [VC]. En algunas ocasiones necesitaremos también algunos
resultados del capítulo II de mi libro de Análisis avanzado [AA] o de mi libro
de Teoría algebraica de Números [TAl].
vii
viii Introducción
resultados algebraicos sobre números naturales suelen tratarlos “de forma indi-
vidualizada”. Por ejemplo, la ley de reciprocidad cuadrática de Gauss [Al 12.14]
es un resultado aplicable a cada par de primos impares individuales. Esto con-
trasta con afirmaciones como la siguiente:
Esta afirmación habla sobre los números libres de cuadrados, pero no dice
nada sobre ninguno de ellos en particular. Los resultados analíticos que vamos a
obtener en este libro comparten mayoritariamente esta naturaleza “estadística”.
Vamos a desarrollar esta idea comentando algunos de los resultados de los que
nos vamos a ocupar en los capítulos siguientes.
El estudio de los números primos es uno de los temas fundamentales de la
teoría de números (tanto algebraica como analítica). La ley de reciprocidad
cuadrática, que acabamos de citar, es un ejemplo de la clase de afirmaciones
que el álgebra puede probar sobre los números primos. El análisis matemático,
por su parte, permite obtener muchos resultados que ponen en evidencia las
regularidades ocultas en la apariencia caótica de la distribución de los números
primos en el seno de los números naturales.
Llamaremos {pn }∞ n=1 a la sucesión de los números primos o, dicho con otras
palabras, llamamos pn al primo n-simo. Sus primeros términos son:
43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, . . .
Se trata de una sucesión muy irregular, en el sentido de que no hay forma de
predecir exactamente cuándo, al ir enumerando los números naturales “aparece”
el siguiente primo. Puede ocurrir que pn+1 = pn +2 y puede ocurrir que haya que
recorrer bastantes números naturales a partir de pn hasta que aparezca pn+1 .
La gráfica muestra la sucesión pn+1 − pn , y en ella vemos que abundan más
las distancias cortas, pero que de vez en cuando se dan también distancias más
largas.
ix
que asigna a cada número real x el número de primos menores o iguales que x.
He aquí su gráfica:
1000
500
π(x)
1000
500 x/ log x
2000 4000 6000 8000 10 000
0.4
0.3
0.2
π(x)/x
0.1
1/ log x
f (x)
lím = 1.
x→+∞ g(x)
f (x) − g(x)
lím = 0,
x→+∞ f (x)
lo que significa a su vez que el error relativo de aproximar f (x) por g(x) tiende
a 0. Discutiremos esta noción con más detalle en la sección 2.2.
Sin embargo, x/ log x no es la función más razonable para aproximar π(x).
Si de verdad 1/ log x aproxima a la densidad de los números primos en cualquier
intervalo, a la hora de aproximar el número de primos entre 1 000 y 1 100 será
más fiable multiplicar la densidad aproximada 1/ log(1 100) por la longitud del
intervalo, o sea, 100, lo cual nos da
100
= 14.28,
log(1 100)
xi
0.10
60 000
0.09
40 000
0.08
20 000 0.07
0.06
200 000 400 000 600 000 800 000 1 × 106 200 000 400 000 600 000 800 000 1 × 106
102 12 3 22 25 29 4 16
107 6.6 44 158 620 421 664 579 664 917 338 0.05
108 5.8 332 773 5 428 682 5 761 455 5 762 208 753 0.013
109 5.1 2 592 591 48 254 943 50 847 534 50 849 234 1 700 0.0033
1010 4.6 20 758 030 434 294 482 455 052 512 455 055 613 3 101 0.00068
El teorema de los números primos está relacionado con muchos otros patrones
regulares que pueden derivarse de la sucesión de los números primos. Ya hemos
señalado uno de ellos, a saber, el parecido de la gráfica de pn /n con la función
log x. Observemos ahora la gráfica de la función1
P
p,
p≤x
es decir la suma de todos los primos menores o iguales que un número real x
dado:
2
x
La que queda ligeramente por debajo es 2 log x . Sucede que para probar que
estas dos funciones son asintóticamente equivalentes a la gráfica de la suma de
los primos no nos bastará con el teorema de los números primos tal y como
lo hemos enunciado, sino que necesitaremos una versión fuerte que estime la
diferencia entre π(x) e Il(x).
Veamos por último lo que sucede si, en vez de 500
sumar los primos, sumamos los logaritmos de los
primos, es decir, si consideramos la función 400
P
ϑ(x) = log p. 300
p≤x
200
Pruebas de trascendencia
1
2 Capítulo 1. Pruebas de trascendencia
con a0 M0 6= 0.
Vamos a encontrar un primo p tal que a0 M0 no sea divisible entre p, mientras
que la suma
Pr ki
P
ai Mki
i=1 k=1
con lo que tendremos una contradicción, pues en (1.1) los dos primeros sumandos
son un entero no divisible entre p, luego no nulo, mientras que el tercero tiene
módulo menor que 1. Para conseguir todo esto definimos primeramente
r ki
r Q
fi (z) = b0 + b1 z + · · · + bN z N = bN
Q Q
f (z) = (z − βki ),
i=1 i=1 k=1
r
P
donde N = ki y b0 6= 0, ya que las raíces son no nulas. Podemos suponer
i=1
que bN > 0.
Sea
+∞ (N −1)p−1 p−1 p
f (z)e−z
Z
bN z
M0 = dz,
0 (p − 1)!
donde p es un número primo y la integración se realiza sobre el semieje real
positivo.
Vamos a probar que la integral es finita y que, si p es suficientemente grande,
se trata de un número entero no divisible entre p. Notemos que
(NP
+1)p
(N −1)p−1 p−1 p (N −1)p−1 p p−1
bN z f (z) = bN b0 z + cs−1 z s−1 ,
s=p+1
1.2. El teorema de Lindemann-Weierstrass 5
(N +1)p
(N −1)p−1 p
X (s − 1)! (N −1)p−1 p
= bN b0 + cs−1 = bN b0 + pC,
s=p+1
(p − 1)!
para un cierto C ∈ Z, donde hemos hecho uso de la identidad de Euler [An 8.60]
Z +∞
n! = z n e−z dz.
0
βki M
- ki
ki
-
M0
βki M
- ki
ki
?
M0 R
6 Capítulo 1. Pruebas de trascendencia
Ahora bien, es fácil ver que la integral sobre el segmento vertical tiende a 0
con R, luego la integral que define Mki es finita, y al sumarle la integral que
define a ki da exactamente M0 . Así pues, Mki + ki = eβki M0 , y tenemos la
descomposición buscada.
Si en la definición de Mki descomponemos en factores el polinomio f (z)
obtenemos
kj
r Q
bN p−1 p−1
(z − βtj )p e−z+βki
Q
Z +∞ N z
j=1 t=1
Mki = dz.
βki (p − 1)!
El camino de integración es z = u + βki , luego dz = du. Al hacer el cambio
la integral se convierte en
kj
r Q
bN p−1
(u + βki )p−1 up e−u ∗
(u + βki − βtj )p
Q
Z +∞ N
j=1 t=1
Mki = du,
0 (p − 1)!
donde el asterisco en el producto indica que falta el factor (t, j) = (k, i), que
hemos extraído como up . Podemos redistribuir los coeficientes bN :
kj
r Q
(bN u + bN βki )p−1 up e−u ∗
(bN u + bN βki − bN βtj )p
Q
Z +∞
j=1 t=1
Mki = du.
0 (p − 1)!
Es fácil ver que, puesto que bN es el coeficiente director de un polinomio
cuyas raíces son los βij , los números αij = bN βij son enteros algebraicos. Con
esta notación:
kj
r Q
(bN u + αki )p−1 up e−u ∗
(bN u + αki − αtj )p
Q
Z +∞
j=1 t=1
Mki = du.
0 (p − 1)!
Sumando obtenemos que
r ki +∞
up Φ(u)e−u
X X Z
ai Mki = du,
i=1 0 (p − 1)!
k=1
donde
r ki kj
r Q
∗
(bN u + αki )p−1 (bN u + αki − αtj )p .
P P Q
Φ(u) = ai
i=1 k=1 j=1 t=1
Si consideramos una cuerpo numérico normal K que contenga a todos los αij ,
resulta que un automorfismo de K permuta los números α1i , . . . , αki i , y se ve
claramente que entonces deja invariante a Φ(u). Esto significa que Φ(u) ∈ Q[u],
y como los αij son enteros algebraicos, en realidad Φ(u) ∈ Z[u]. Digamos que
(NP
+1)p
up Φ(u) = ds−1 us−1 .
s=p+1
1.2. El teorema de Lindemann-Weierstrass 7
Entonces
r ki (N +1)p +∞ (N +1)p
(s − 1)!
Z
X X X ds−1 X
ai Mki = us−1 e−u du = ds−1 = pC,
i=1 s=p+1
(p − 1)! 0 s=p+1
(p − 1)!
k=1
donde los exponentes β1 , . . . , βN son distintos dos a dos (es decir, donde los
coeficientes Bi se obtienen multiplicando un Aki para cada i y después sumando
todos los productos que acompañan a un mismo exponente), se cumple que al-
guno de los coeficientes Bi es no nulo.
donde los coeficientes bi son enteros racionales obtenidos sumando los ch1 ,...,hn
que acompañan a un mismo exponente, es decir, según las hipótesis del teorema
anterior, por lo que alguno de ellos es no nulo (y claramente ha de haber al
menos dos no nulos). Los números bi son números algebraicos distintos, luego
tenemos una expresión como la original pero con coeficientes enteros.
A partir de ahora suponemos (1.3) con ci ∈ Z y donde α1 , . . . , αn son nú-
meros algebraicos distintos.
Sea f (x) ∈ Q[x] el producto de los polinomios mínimos de los números αi
(sin repetir dos veces el mismo factor). Sea m ≥ n el grado de f , sean γ1 , . . . , γm
todas las raíces de f . Llamemos µ = m(m − 1) . . . (m − n + 1), al número de
n-tuplas posibles (i1 , . . . , in ) de números distintos comprendidos entre 1 y m.
Entonces
(c1 eγi1 + · · · + cn eγin ) = 0,
Q
Sumando resulta
r
(c1 eγi1 + · · · + cn eγin ) Bk (eγ1k + · · · + eγtk k )
Q P
=
k=1
q
Ai (eβ1i + · · · + eβti i ) = 0,
P
=
i=1
10 Capítulo 1. Pruebas de trascendencia
donde
r
P
Ai = pik Bk ∈ Z.
k=1
Notemos que todos los números βij son distintos, pues son raíces de polino-
mios irreducibles distintos. Por construcción, los exponentes βij son todas las
sumas distintas de exponentes γij que aparecen al efectuar el producto de la
izquierda de la ecuación. Podemos aplicar el teorema anterior y concluir que
alguno de los coeficientes Ai es no nulo. Eliminando los nulos podemos suponer
que ninguno lo es. En resumen tenemos
q
Ai (eβ1i + · · · + eβti i ) = 0,
P
i=1
donde los coeficientes son enteros racionales no nulos y los exponentes de cada
sumando son familias de números conjugados correspondientes a polinomios
irreducibles distintos si (x) de grado ti . Distinguimos dos casos:
1) Algún βki = 0. Pongamos por ejemplo i = 1. Esto significa que s1 (x) = x,
luego además t1 = 1 y la ecuación se reduce a
q
Ai (eβ1i + · · · + eβti i ) = 0,
P
A1 +
i=2
donde los exponentes son todos no nulos, y esto contradice al teorema 1.3.
2) Todos los βki son distintos de 0. Dividimos la ecuación entre eβk1 para
k = 1, . . . , t1 , con lo que obtenemos las ecuaciones
q ti
eβti −βk1 = 0,
P P
Ai k = 1, . . . , t1 .
i=1 t=1
q ki
eαki = 0,
P P
A0 + Ai
i=1 k=1
2. El número π es trascendente.
Si π fuera algebraico también lo sería iπ, y el número eiπ = −1 sería
trascendente.
(β − 1)eiα + (β + 1)e−iα = 0,
Ejercicio: Probar que las funciones arcsen, arccos y arctan toman valores trascen-
dentes sobre números algebraicos (salvo casos triviales).
12 Capítulo 1. Pruebas de trascendencia
Las primeras pruebas del teorema fueron publicadas por S. Lang y K. Ra-
machandra, si bien un caso particular fue probado por Alaoglu y Erdős en 1944
en un artículo en el que afirmaban que Siegel conocía una prueba.
Vamos a necesitar algunos resultados adicionales sobre cuerpos numéricos.
Si K es un cuerpo numérico normal1 de grado h y σ1 , . . . , σh son sus automor-
fismos, podemos definir la norma y la traza [Al 7.39]:
N : K −→ Q, Tr : K −→ Q
mediante N(α) = σ1 (α) · · · σh (α), Tr(α) = σ1 (α) + · · · + σh (α). La traza define
a su vez una forma bilineal [Al 10.1] K ×K −→ Q dada por (α, β) 7→ Tr(αβ), de
modo que si α1 , . . . , αh es una Q-base de K, podemos considerar su base dual
β1 , . . . , βh , caracterizada por que
1 si i = j,
Tr(αi βj ) =
0 si i 6= j.
Así, si α = a1 α1 + · · · + ah αh es un elemento arbitrario de K, tenemos que
ai = Tr(αβi ). Definimos |α| como el máximo de los módulos de los conjugados
de α, y así
|ai | = | Tr(αβi )| ≤ h|αβi | ≤ h|βi | |α|.
Fijamos la base α1 , . . . , αh con la condición adicional de que sea una entera
[Al 10.15], es decir, una base del anillo de los enteros de K como Z-módulo. Así,
si llamamos cK al producto de h por el máximo de los módulos de los conjugados
de los elementos de su base dual, tenemos que las coordenadas de un α ∈ K en
dicha base están acotadas por:
|ai | ≤ cK |α|.
Por último, si α ∈ K es entero no nulo, todos sus conjugados también lo
son, luego también lo es N(α) y, al ser también un número racional, es un
entero racional. Por consiguiente:
h−1
1 ≤ | N(α)| ≤ |α| |α|. (1.4)
Aquí hemos acotado por |α| todos los conjugados de α excepto él mismo.
Vamos a necesitar una generalización del siguiente resultado elemental:
Teorema 1.8 (Lema de Siegel) Sea (ajk ) una matriz M ×N con coeficientes
enteros racionales tal que M < N y de modo que todos los coeficientes estén
acotados en módulo por A ≥ 1. Entonces el sistema de ecuaciones lineales
aj1 x1 + · · · + ajN xN = 0, 1 ≤ j ≤ M,
tiene una solución entera no trivial tal que
|xk | ≤ (N A)M/(N −M ) ,
para 1 ≤ k ≤ N .
1 En realidad la hipótesis de normalidad no es necesaria.
14 Capítulo 1. Pruebas de trascendencia
yj = aj1 x1 + · · · + ajN xN , 1 ≤ j ≤ M.
Bj + Cj = |aj1 | + · · · + |ajN | ≤ N A,
y claramente −Bj H ≤ yj ≤ Cj H.
Ahora bien, el número de N -tuplas (x1 , . . . , xN ) tales que 0 ≤ xk ≤ H
es (H + 1)N , mientras que sus M -tuplas asociadas (y1 , . . . , yM ) varían en un
conjunto de a lo sumo (Cj H +Bj H +1)M ≤ (N AH +1)M < (H +1)N elementos,
luego tiene que haber dos N -tuplas distintas con la misma imagen. Su diferencia
cumple el teorema.
Ahora generalizamos el resultado para matrices de enteros algebraicos:
Teorema 1.9 Sea (αkl ) una matriz p × q con coeficientes enteros en un cuerpo
numérico K tal que p < q y de modo que |αkl | ≤ A. Entonces el sistema de
ecuaciones lineales
αk1 ξ1 + · · · + αkq ξq = 0, 1 ≤ k ≤ p,
tiene una solución entera (en K) no trivial tal que |ξl | ≤ c(cqA)p/(q−p) , para
1 ≤ l ≤ q, donde c es una constante que depende de K y de una base entera
β1 , . . . , βh de K, pero no de la matriz.
ξl = xl1 β1 + · · · + xlh βh , 1 ≤ l ≤ q,
Entonces
q
P q
P h
P P q
h P h
P
αkl ξl = αkl xlr βr = xlr aklru βu
l=1 l=1 r=1 r=1 l=1 u=1
h P
P q
h P
= aklru xlr βu ,
u=1 r=1 l=1
1.3. El teorema de las seis exponenciales 15
luego (ξ1 , . . . , ξq ) será solución del sistema de ecuaciones si y sólo si las coorde-
nadas (xl1 , . . . , xlh ) son solución del sistema de M = hp ecuaciones con N = hq
incógnitas
Ph Pq
aklru xlr = 0, 1 ≤ u ≤ h, 1 ≤ k ≤ p.
r=1 l=1
Por el teorema anterior este sistema de ecuaciones tiene una solución entera no
trivial tal que
Los (ξ1 , . . . , ξq ) con estas coordenadas son enteros de K no todos nulos que
cumplen el sistema de ecuaciones y además
Teorema 1.10 Sea (αkl ) una matriz p × q con coeficientes en un cuerpo numé-
rico K tal que p < q y de modo que |αkl | ≤ A. Entonces el sistema de ecuaciones
lineales
αk1 ξ1 + · · · + αkq ξq = 0, 1 ≤ k ≤ p,
tiene una solución entera (en K) no trivial tal que |ξl | ≤ c(cqdA)p/(q−p) , para
1 ≤ l ≤ q, donde d es una cota del mínimo común denominador de los coefi-
cientes de cada ecuación y c es una constante que depende de K y de una base
entera β1 , . . . , βh , pero no de la matriz.
donde los aij serán enteros de K que vamos a elegir de modo que F se anule en
todos los puntos de la forma
k1 y1 + k2 y2 + k3 y3 , 1 ≤ ki ≤ n,
para un cierto n elegido adecuadamente, al igual que r.
Las condiciones F (k1 y1 + k2 y2 + k3 y3 ) = 0 forman un sistema de n3 ecuacio-
nes lineales con r2 incógnitas (las aij ). Pretendemos aplicar el teorema anterior,
por lo que necesitamos que r2 > n3 . Concretamente, tomaremos r2 = (4n)3 ,
Notemos que siempre podemos encontrar números que cumplan esto con n arbi-
trariamente grande. Los coeficientes del sistema de ecuaciones son los números
algebraicos
(ex1 y1 )ik1 (ex1 y2 )ik2 (ex1 y3 )ik3 (ex2 y1 )jk1 (ex2 y2 )jk2 (ex2 y3 )jk3
cuyo mínimo común denominador está acotado por d6rn y, si ec0 es una cota2
de los conjugados de las seis exponenciales exu yv , entonces los conjugados de
los coeficientes están acotados por ec0 rn . Así, el teorema anterior nos da que
existen enteros aij ∈ K que hacen que F tenga los ceros requeridos y de modo
que
3
/(r 2 −n3 ) 5/2
|aij | ≤ C(Cr2 d6rn ec0 rn )n = C(4Cn3 d6rn ec0 rn )1/63 ≤ ec1 n .
Detallamos la última desigualdad, si bien en lo sucesivo omitiremos estas
comprobaciones rutinarias: Teniendo en cuenta que r = 8n3/2 , tenemos que
1
log |aij | ≤ log C + log 4C + 3 log n + 48 log d n5/2 + 8c0 n5/2 ,
63
es decir,
5/2 C1 log n
log |aij | ≤ C1 + C2 log n + C3 n = + C2 5/2 + C3 n5/2 ,
n5/2 n
y la expresión entre paréntesis tiende a C3 , luego está acotada por una cons-
tante c1 .
Como x1 , x2 son linealmente independientes sobre Q, los coeficientes ix1 +jx2
son distintos dos a dos, luego el polinomio
r
aij X (ix1 +jx2 )
P
P (X) =
i,j=1
| N(d6r(s+1) F (w))| ≥ 1.
donde la función del miembro derecho es entera, ya que los polos simples del
producto se cancelan con ceros de F . Consideremos el disco de centro 0 y radio
R = s3/2 .
Notemos que |w| ≤ 3(s + 1) máx{|yi |}, por lo que, tomando n es suficien-
temente grande (y, por consiguiente, s también), podemos exigir que |w| < R.
Más aún, podemos exigir que si |z| = R y 1 ≤ ki ≤ s + 1, entonces
|z − (k1 y1 + k2 y2 + k3 y3 )| ≥ R/2.
luego
1
log |F (w)| ≤ log FR + s3 (log c2 − 2 log s).
Nos falta estimar FR , que depende de s:
5/2 5/2
+c3 8n3/2 s3/2 3
FR ≤ r2 ec1 n ec3 rR ≤ 16n3 ec1 n ≤ ec 4 s .
(log α)−k ρki eρi l = (log α)−k (a + bβ) log α)k el(a+bβ) log α = (a + bβ)k αal γ bl ∈ K,
con 1 ≤ l ≤ m, 1 ≤ a, b ≤ q, 0 ≤ k ≤ n − 1.
Sea c1 un número natural no nulo tal que c1 α, c1 β y c1 γ sean enteros en K.
En cada coeficiente, al desarrollar el binomio (a + bβ)k aparecen monomios de
α, β y γ con grado a lo sumo
q mr
|S(x)| ≤ r! cr10 r(r+3)/2 cr11
√r mr −mr/2
≤ rr cr10 r(r+3)/2 cr11 2m r = cr12 r(3r+3−mr)/2 .
Pero combinando esto con (1.6) resulta c−r 5 < cr14 r(3h−r)/2 , o lo que es lo
(r−3h)/2 r r r
mismo, r < c14 c5 = c15 . Tomando logaritmos es fácil llegar a que
1 3h
− log r < log c15 .
2 2r
Hemos probado que esto se cumple para una constante c15 y para valores de r
arbitrariamente grandes (pues r ≥ n), pero esto es claramente contradictorio,
pues el miembro de la izquierda tiende a +∞ cuando r tiende a +∞.
1.4. El teorema de Gelfond-Schneider 21
Por ejemplo, ahora podemos asegurar que un logaritmo logα β con α, β alge-
braicos es racional o trascendente, pues si si fuera algebraico irracional entonces
β = αlogα β sería trascendente. Por ejemplo, es fácil ver que log2 3 es trascen-
dente, ya que no puede ser un número racional a/b. En tal caso 2a/b = 3, luego
2a = 3b , lo cual es imposible, con a, b 6= 0.
En particular vemos que un número algebraico elevado a un número tras-
cendente puede ser algebraico. Por otra parte, por una mera cuestión de car-
dinalidad, es claro que αβ tiene que ser trascendente para muchos números
trascendentes β.
Capítulo II
Funciones aritméticas
23
24 Capítulo 2. Funciones aritméticas
f (d)f −1 (n/d) = 0,
P
d|n
Ejemplo En [Al 6.13] demostramos que todo subgrupo finito del grupo multi-
plicativo de un cuerpo es cíclico (lo que en particular implica que existen raíces
primitivas módulo todo primo p), para lo cual nos apoyamos en los teoremas de
estructura de los grupos abelianos finitos. Vamos a dar una prueba mucho más
elemental basada en el teorema anterior y en la observación siguiente:
Si G es un subgrupo de orden q de un cuerpo F y n | q, el polinomio xn − 1
tiene n raíces distintas en G.
En efecto, considerando congruencias módulo P = xn − 1 en F [x] vemos que
Así pues, para cada n | q, resulta que G tiene φ(n) elementos de orden n y, en
particular, tiene φ(q) ≥ 1 elementos de orden q.
Veamos un segundo ejemplo de aplicación de la fórmula de inversión:
(f + g)0 = f 0 + g 0 , (f ∗ g)0 = f 0 ∗ g + f ∗ g 0 .
Notemos que c01 = log (donde el miembro derecho es la función n 7→ log n).
Introducimos ahora una función aritmética que tendrá gran importancia en
lo sucesivo:
f ◦ (g ◦ F ) = (f ∗ g) ◦ F,
entonces
f −1 (n)G(x/n).
P
F (x) =
n≤x
f −1 ◦ G = f −1 ◦ (f ◦ F ) = (f −1 ∗ f ) ◦ G = 1 ◦ F = F.
En particular:
entonces
P
F (x) = µ(n)f (n)G(x/n).
n≤x
Veamos un ejemplo:
Teorema 2.14 El número de fracciones irreducibles 0 < a/b < 1 con 0 < b ≤ n
viene dado por
1 P
f (n) = µ(m)E[n/m](E[n/m] − 1).
2 m≤n
donde hemos usado que E[n/m] = E[x/m]. En efecto, es evidente que se cumple
E[n/m] ≤ E[x/m] y, si la desigualdad fuera estricta, tendríamos que
n/m < E[n/m] + 1 ≤ E[x/m] ≤ x/m,
luego n < m(E[n/m] + 1) ≤ x, luego E[x] + 1 = n + 1 ≤ x, contradicción. El
teorema anterior nos da que
P
F (x) = µ(m)G(x/m),
m≤x
Así, en el miembro derecho sólo están los términos π(x1/m )/m correspondientes
a índices m impares. A su vez,
∞ ∞
1 1/3 1 ∗ 1/3 1 ∗ 1/6 X π(x1/3m ) X π(x1/6m )
f1 (x ) = π (x ) − π (x ) = − ,
3 3 6 m=1
3m m=1
6m
luego
1 1 1
f2 (x) = π ∗ (x) − π ∗ (x1/2 ) − π ∗ (x1/3 ) + π ∗ (x1/6 )
2 3 6
∞ ∞ ∞ ∞
X π(x1/m ) X π(x1/2m ) X π(x1/3m ) X π(x1/6m )
= − − + .
m=1
m m=1
2m m=1
3m m=1
6m
Vemos que en la primera expresión para f2 (x) aparecen todos los índices m
libres de cuadrados divisibles únicamente entre 2 y 3 y el signo es µ(m), mientras
que en la segunda expresión aparecen todos los índices m no divisibles ni entre
2 ni entre 3 (como los múltiplos de 6 se restan dos veces, luego se suman para
que sólo sean sustraídos una vez).
Una inducción rutinaria muestra que esto es cierto en general, es decir, que
fk (x) admite una expresión con sumandos µ(m) ∗ 1/m
m π (x ), donde m recorre todos
los números divisibles únicamente entre los primos p1 , . . . , pk (libres de cuadra-
dos, si se quiere, porque en caso contrario µ(m) = 0) y también una segunda
expresión con sumandos (1/m)π(x1/m ) para todos los índices no divisibles entre
ninguno de los primos p1 , . . . , pk .
Así, cuando k es suficientemente grande, la primera expresión para fk coin-
cide con el miembro derecho de la fórmula del enunciado (pues todos los suman-
dos que faltarían por aparecer serían ya nulos) y la segunda expresión se reduce
a π(x).
f (x)
lím = 1.
x→a g(x)
o también a que
Las dos últimas expresiones se interpretan como que el error relativo cometido
al aproximar una por la otra (es decir, el error absoluto dividido entre el valor
real que queríamos aproximar) tiende a 0.
Ejemplo 1 Se cumple que x2 + x ∼ x2 (en +∞). El error absoluto cometido al
aproximar una por la otra tiende a infinito, pero dicho error es cada vez menor
en relación al valor que toman ambas funciones.
10 000
100
8000
80
60 6000
40 4000
20 2000
2 4 6 8 10 20 40 60 80 100
Siempre que empleemos esta notación habrá que entender que n recorre los
números naturales no nulos menores o iguales que el número real x. Equivalen-
temente,
E[x]
P P E[x](E[x] + 1)
n= n= .
n≤x n=1 2
Se cumple entonces que
P x2
n∼ ,
n≤x 2
pues
P
n
|E[x]2 + E[x] − x2 |
n≤x (x − E[x])(x + E[x]) + E[x]
− 1 = ≤
2 2
x /2 x x2
1(x + x + 1) + x + 1 3x + 2
≤ = → 0.
x2 x2
50
4
40
2
30
20 2 4 6 8 10
10 -2
-4
2 4 6 8 10
2000 20 40 60 80 100
-20
1000
-40
20 40 60 80 100
Así pues, x2 /2 es una función sencilla que describe bien el crecimiento de las
sumas finitas de números naturales.
En los dos ejemplos anteriores hemos comparado dos funciones que tienden
a +∞ en +∞. Será el caso más habitual, pero también tiene interés estudiar el
modo en que una función tiende a 0 en un punto finito o infinito. Por ejemplo,
podemos afirmar que
sen x ∼ x alrededor de 0.
36 Capítulo 2. Funciones aritméticas
0.5
-0.5
-1.0
Es frecuente usar una notación un tanto abusiva para expresar este hecho,
consistente en escribir O(g(x)) para referirse a una función arbitraria que crezca
como g(x), es decir, una función arbitraria que permanece acotada cuando se
divide entre |g(x)|.
Ejemplo 1 Hemos dicho que esta definición es más débil que la anterior porque,
claramente, si f (x) ∼ g(x) entonces el crecimiento de f es del orden del de g
(pues la definición se cumple con cualquier constante mayor que 1), y el recíproco
no es cierto en general. Sin embargo, con el concepto de orden de crecimiento
podemos expresar relaciones más precisas que con el de equivalencia asintótica.
Por ejemplo, decir que
n = O(x2 /2)
P
n≤x
es más débil que decir que ambas funciones son asintóticamente equivalentes,
pues sólo expresa que el cociente está acotado, no que tienda a 1. De hecho,
esto es equivalente a
n = O(x2 ),
P
n≤x
pues al dividir entre x2 /2 una función que permanece acotada al dividirla entre
x obtenemos una función que permanece acotada al dividirla entre 2/x, luego
también si la dividimos entre 1/x.
Vemos que la expresión (2.1) no sólo dice que las dos primeras funciones son
asintóticamente equivalentes, sino que nos acota el error absoluto de aproximar
una por la otra, y la expresión (2.2) no sólo indica explícitamente que el cociente
de la izquierda tiende a 1, sino que nos informa de la rapidez con la que tiende
a 1 (al menos tan rápidamente como 1/x tiende a 0).
En general, cuando comparamos dos funciones que tienden a +∞, lo que
expresa la definición es que f (x) deja de tender a +∞ cuando se la divide entre
g(x) o, en otras palabras, que g(x) crece con la rapidez suficiente como para
compensar el crecimiento de f (x).
P
Así, la serie n tiende a +∞, pero deja de tender a +∞ y permanece
n≤x
acotada cuando se la divide entre x2 , o el error de aproximar dicha serie por
x2 /2 tiende a +∞, pero deja de tender a +∞ cuando se divide entre x.
Ejemplo 2 Consideremos ahora el caso de una serie convergente. Se cumple
que
X 1
= 1 + O(1/2x ).
2n
n≤x
luego P
(1/2n ) − 1
n≤x
= 2x−E[x] ≤ 2,
1/2x
38 Capítulo 2. Funciones aritméticas
pero con la estimación que hemos hecho no sólo sabemos que la serie converge
a 1, sino que tenemos información sobre la velocidad con la que converge (al
menos tan rápidamente como 1/2x converge a 0).
Ejemplo 3 Considerando la serie de Taylor del seno es fácil ver que, alrededor
de 0,
x3
sen x = x − + O(x5 ), (2.3)
6
lo cual significa, por definición, que el crecimiento de la función sen x − x + x3 /6
es como el de x5 , es decir, que el cociente permanece acotado. De hecho, se
cumple que
sen x − x + x3 /6 1
lím = .
x→0 x5 5!
Si dividiéramos entre x4 el límite sería 0 y si dividiéramos entre x6 el límite
sería infinito, por lo que 5 es “el exponente justo” que da un cociente acotado
que no tiende ni a 0 ni a +∞. De (2.3) se sigue inmediatamente que
sen x x2
=1− + O(x4 ),
x 6
pues si dividimos entre x una función que permanece acotada al dividirla en-
tre x5 , obtenemos una función que permanece acotada al dividirla entre x4 , y
esto implica a su vez que sen x ∼ x, por lo que (2.3) proporciona información
más precisa sobre el modo en que sen x tiende a 0 en a = 0.
En lugar de (2.3), podríamos haber escrito
x3
sen x = x − + O(x4 ),
6
lo cual, aunque es cierto, es menos preciso, porque la función que estamos repre-
sentando por O(x4 ) (que es la misma que en (2.3) representamos por O(x5 ), a
saber, sen x − x + x3 /6), no sólo es una función que permanece acotada cuando
se divide entre x4 , sino que de hecho tiende a 0. Para hacer constar este hecho
usaremos una o minúscula:
x3
sen x = x − + o(x4 ).
6
2.2. Orden de crecimiento 39
eβx x
lím α
= lím e(β log x −α) log x = +∞,
x→+∞ x x→+∞
40 Capítulo 2. Funciones aritméticas
• Por debajo de todos los crecimientos potenciales O(xα ) están los creci-
mientos logarítmicos O(logβ x). En efecto, aplicando la regla de L’Hôpital,
β β
logβ x
log x 1
lím = lím = lím = 0,
x→+∞ xα x→+∞ xα/β x→+∞ (α/β)x1+α/β
e2x ex ex/2
100 x2 x1.1 x
80
x0.9
60
40
√
x
20
log2 x
log x
√
log x
0 20 40 60 80 100
Para la segunda parte observamos que si n ≤ t < n+1 entonces C(t) = C(n),
luego Z n+1
C(n)(f (n) − f (n + 1)) = − C(t)f 0 (t) dt,
n
y C(t) = 0 si t < n1 . Así, la integral desde n1 hasta N = E[x] es el sumatorio
de la expresión que hemos obtenido, mientras que
Z x
− C(t)f 0 (t) dt = −C(x)(f (x) − f (E[x])) = C(x)f (E[x]) − C(x)f (x),
N
X 1 π2
2. 2
= + O(1/x),
n 6
n≤x
X µ(n) 6
3. = 2 + O(1/x),
n2 π
n≤x
42 Capítulo 2. Funciones aritméticas
X 3 2
4. φ(n) = x + O(x log x).
π2
n≤x
lo que en particular implica que γ cumple (2.4), por lo que se trata de la cons-
tante de Euler que ya conocemos.
2) Ahora tomamos f (t) = 1/t2 , con lo que
Z x Z x Z x
X 1 E[x] E[t] E[x] t − E[t] 1
= + 2 dt = − 2 dt + 2 dt
n2 x2 1 t 3 x 2
1 t3
1 t 2
n≤x
x
t − E[t]
Z
E[x] 2
=2+ 2 − −2 dt.
x x 1 t3
2.3. Cálculo de órdenes 43
2 x−1 x+1
− ≤ ≤ = O(1/x).
x x2 x2
3) En 2.8 hemos probado que la serie converge a 6/π 2 . Ahora se trata de
estimar las aproximaciones por sumas parciales:
π2 X 1
6 X µ(n) X µ(n) X 1
π2 − ≤ −
= = = O(1/x).
2 2
n n>x n n>x n
2 6 n2
n≤x n≤x
2
n − x2 , la relación (2.1) equivale a que |f (x)| ≤ Cx.
P
Si llamamos f (x) =
n≤x
Usando los apartados anteriores concluimos que
2
x2 X µ(d)
X X x Xx
φ(n) = µ(d) + f (x/d) ≤ + C
2d2 2 d2 d
n≤x d≤x d≤x d≤x
x2
6 3
= 2
+ O(1/x) + Cx(log x + O(1)) = 2 x2 + O(x) + O(x log x) + O(x)
2 π π
3
= 2 x2 + O(x log x).
π
Hemos usado que si dividimos entre x log x la suma de dos funciones que perma-
necen acotadas cuando se dividen entre x y una que permanece acotada cuando
se divide entre x log x, obtenemos la suma de dos funciones que tienden a 0 y
otra acotada, luego la suma total está acotada.
Veamos un par de aplicaciones:
E[y 2 ] = L((y/m)2 ).
P
m≤y
L(y 2 ) = µ(m)E[(y/m)2 ].
P
m≤y
luego
X µ(m) 6y 2 6y 2
L(y 2 ) = y 2 + O(y) = + y 2
O(1/y) + O(y) = + O(y).
m2 π2 π2
m≤y
Equivalentemente:
6x √
L(x) = 2
+ O( x).
π
Esto implica en particular que la probabilidad es la indicada.
2.3. Cálculo de órdenes 45
o, equivalentemente, el de la función
1 X 6 2 1 6 log x 1
2
2 φ(n) − 1 = 2
+ 2
O(x log x) − 2
= 2
+ O − 2.
x π x x π x x
n≤x
30
20
10
Por otra parte, es obvio que d(n) puede tomar valores arbitrariamente grandes.
La mayor dificultad a la hora de estimar el crecimiento de d(n) es que de-
pende de la descomposición en factores primos de n, y la factorización de un
número no tiene nada que ver con la del número anterior o posterior, sino que
después de un número con muchos divisores puede venir un primo o viceversa.
Por ello el estudio del orden de crecimiento de d(n) es un problema más com-
plicado que el de todas las funciones que hemos abordado hasta ahora. Pese a
ello, estamos en condiciones de aproximarlo con precisión. De hecho, lo difícil
es “no pasarse”, porque vamos a ver, por ejemplo, que
d(n) = o(n ), para todo > 0, (2.5)
de modo que si dividimos d(n) entre n , no sólo hacemos que la sucesión pase
a estar acotada, sino que de hecho la hemos hecho tender a 0. Una estimación
más fina es la siguiente:
Teorema 2.22 Para todo > 0 existe un n0 tal que, si n ≥ n0 , entonces
log n
log d(n) ≤ (1 + ) log 2 .
log log n
Equivalentemente,
log d(n) log log n
lím ≤ log 2,
n log n
o también log n (1+) log n (1+) log 2
d(n) ≤ e(1+) log 2 log log n = 2 log log n =n log log n .
De aquí se sigue (2.5), pues
d(n) (1+) log 2
≤ e( log log n −) log n → 0.
n
Demostración: Dado n = pe11 · · · perr , consideremos un δ > 0 que luego
especificaremos. Entonces
r
d(n) Y ei + 1
= .
nδ i=1
pei i δ
La serie de Taylor de ex muestra que, si x > 0, entonces 1 + x ≤ ex , por lo que
ei δ log 2 ≤ eei δ log 2 = 2ei δ ≤ pei i δ .
ei + 1 ei + pei i δ ei 1
ei δ
≤ ei δ
= 1 + ei δ ≤ 1 + ≤ e1/(δ log 2) .
pi pi pi δ log 2
2.3. Cálculo de órdenes 47
ei + 1 ei + 1
ei δ
≤ .
pi 2 ei
Por lo tanto,
d(n) Y ei + 1 Y ei + 1 Y 1/δ
= ei δ ei δ
≤ e1/(δ log 2) ≤ e2 /(δ log 2) .
nδ pi pi
pi <21/δ p ≥21/δ
i p≤21/δ
Equivalentemente,
21/δ
log d(n) − δ log n ≤ .
δ log 2
Ahora tomamos, concretamente,
(1 + /2) log 2
δ= ,
log log n
con lo que
log log n log log n
21/δ = 2 (1+/2) log 2 = e 1+/2 = (log n)1/(1+/2) .
Por lo tanto,
log2 log n 3
/2
≤ 1+ log 2 (2.6)
(log n) 1+/2 2 2
d(n)
log x
e(1.2) log 2 log log x
log x
elog 2 log log x
lo primero que cabe pensar es que el teorema es falso, pues deberíamos ver que,
a partir de cierto valor, la función d(n) queda por debajo de la mayor de las dos
curvas. Si extendemos la gráfica hasta, digamos, n = 106 , la imagen es similar.
La razón es que hemos probado que la función d(n) debe quedar por debajo
de la curva mayor para todo n suficientemente grande, pero ¿cómo de grande?
El único punto de la prueba que requiere que n sea grande es (2.6) que, para
= 0.2, es, concretamente
log2 log n
≤ 0.037.
(log n)0.09
log2 x
Ahora bien, llamando x = log n, ésta es la gráfica de :
x0.09
70
60
50
40
30
20
10
En particular,
P
d(k)
k≤x
lím − log x = 2γ − 1 = 0.1544313298 . . .
x→+∞ x
Demostración: Observemos que d(k) P es el número de pares de números
naturales (u, v) tales que uv = k, luego d(k) es el número de pares de
k≤x
números naturales no nulos tales que uv ≤ x.
Equivalentemente, es el número de pa- 14
contradicción.
50 Capítulo 2. Funciones aritméticas
√ √ √ 2 √ √
= 2x(log x + γ + O(1/ x)) − 2( x /s + O( x)) + O( x))
√
= x log x + (2γ − 1)x + O( x).
Demostración: Por el teorema anterior existe una constante c > 0 tal que
X
d(k) ≤ cx log x.
k≤x
X X X X X X X
d2 (n) = d(n) 1= d(km) ≤ d(k)d(m)
n≤x n≤x k|n k≤x km≤x k≤x m≤k/x
X X X x X d(k)
= d(k) d(m) ≤ c d(k) log(x/k) ≤ cx log x
k k
k≤x m≤x/k k≤x k≤x
Z x Z x
1 1 log t
≤ cx log x +c t log t dt = c log x + c dt
x 1 t2 1 t
c
= c log x + log2 x = O(log2 x).
2
2.3. Cálculo de órdenes 51
xd m0 d m0 d xd xd 2xd 2x
N (m, n) ≤ E[ − 0 ] − E[− 0 ] = E[ ]< +1≤ = .
mn n n mn mn mn dm0 n0
Por consiguiente,
X X1 X 1
d(r)d(r + k) < 4x
d m0 n0
1≤r≤x−k d|k m0 ,n0 ≤x
X 1 X 1 2 X1
≤ 4x ≤ 4x (log x + 1)2 ,
d m d
d|k 1≤m≤x d|k
Rm P
donde hemos acotado 1/m ≤ m−1 (1/t) dt, con lo que (1/m) ≤ log x.
Esto termina la prueba de (2.7). 1≤m≤x
aplicando el teorema 2.20 igual que antes. Por otra parte, llamando k = n − m
y aplicando (2.7),
E[x−1]
X d(m)d(n) X 1 E[x−k]
X
= d(m)d(m + k)
n−m k m=1
m<n≤x k=1
E[x−1]
X 1X1 X 1
< 4x(log x + 1)2 ≤ 4x(log x + 1)2
k d 0
d2 d0
k=1 d|k dd ≤x
X 1 X 1
≤ 4x(log x + 1)2
d2 0 d0
d≤x d ≤x
∞
X 1 X 1
≤ 4x(log x + 1)2 = O(x log3 x),
d2 0 d0
d=1 d ≤x
P 1
donde nuevamente hemos acotado d0 ≤ 1 + log x = O(log x).
d0 ≤x
121212 = 212121.
n/d (n/d)!
m(d) = PRr1 /d,...,rk /d = .
(r1 /d)! · · · (rk /d)!
1X Xh 1X X h/l
= m̄(l) µ(d/l) = m̄(l) µ(t)
n d n t
l|h l|d|h l|h t|(h/l)
1X 1X
= m(h/l)(N ∗ µ)(h/l) = m(d)φ(d).
n n
l|h d|h
La distribución de los
números primos
3.1 Preliminares
Empezamos probando algunos hechos variados que vamos a necesitar en
las secciones siguientes. Por ejemplo, en la introducción hemos presentado dos
versiones del teorema de los números primos:
Z x
x dt
π(x) ∼ ∼ Il(x) = .
log x 2 log t
luego Z x
Il(x) 2 log x log x dt
=1− + .
x/ log x log 2 x x 2 log2 t
57
58 Capítulo 3. La distribución de los números primos
entonces
X bn x x
=α +o .
log n log x log x
2≤n≤x
P
Demostración: Llamando C(x) = bn , por 2.20 tenemos que
2≤n≤x
Z x
X bn 1 C(t) 1
= C(x) + dt.
log n log x 2 t log2 t
2≤n≤x
Por lo tanto,
Como el cociente C(t)/t converge a α, está acotado, digamos por M , con lo que
y en la prueba del teorema anterior hemos visto que el miembro derecho tiende
a 0.
Un resultado elemental que vamos a necesitar y que tiene interés en sí mismo
es la descomposición en factores primos de los factoriales:
3.1. Preliminares 59
pj(n,p) ,
Q
n! = (3.1)
p
E[n/pm ].
P
donde j(n, p) =
m≥1
Basta tener en cuenta que en n! hay E[n/p] múltiplos de p, de los cuales
E[n/p2 ] son múltiplos de p2 , etc.
Otro hecho interesante es la cota siguiente del producto de los primeros
primos:
p < 4n .
Q
Teorema 3.3
p≤n
Si n = 2m + 1, entonces, como
2m + 1 (2m + 1)!
= ,
m m!(m + 1)!
Por lo tanto,
p ≤ 4m+1 4m = 4n .
Q Q Q
p= p
p≤n p≤m+1 m+1<p≤n
luego
π(x) ϑ(x) 1
= + O( ). (3.2)
x/ log x x log x
Esto ya muestra que el teorema de los números primos es equivalente a ϑ(x) ∼ x.
Nos ocupamos ahora de la segunda función de Chebyshev. Observemos en
primer lugar que, según la definición de la función de Mangoldt, se cumple que
P P
ψ(x) = Λ(n) = log p,
n≤x pm ≤x
donde en la segunda expresión hay que entender que hay un sumando log p por
cada potencia pm ≤ x. Equivalentemente,
E[ log x
P
ψ(x) = log p ] log p.
p≤x
62 Capítulo 3. La distribución de los números primos
Notemos que la serie es finita, pues los sumandos son nulos a partir del momento
en que x1/m < 2, es decir, cuando m > log x/ log 2.
1
El número de primos p ≤ x1/m es a lo sumo x1/m , y log p ≤ m log x ≤ log x,
luego, para m ≥ 2, tenemos que
luego
1 1/2
ϑ(x1/m ) ≤ x log2 x,
P
m≥2 log 2
(puesto que el número de sumandos es a lo sumo log x/ log 2). Más brevemente:
o también
ψ(x) = ϑ(x) + O(x1/2 log2 x). (3.3)
Esto implica a su vez que
3. ϑ(x) ∼ x.
4. ψ(x) ∼ x.
En 3.34 probaremos que todas las afirmaciones del teorema anterior son
ciertas, pero podemos probar otras propiedades más débiles de las funciones de
Chebyshev que, no obstante, tienen repercusiones notables sobre la distribución
de los números primos. Por ejemplo, tomando logaritmos en la desigualdad del
teorema 3.3 (aplicado a n = E[x]) obtenemos el teorema siguiente:
Teorema 3.8 Para todo x ≥ 1, se cumple que ϑ(x) < 2x log 2. En particular,
ϑ(x) = O(x).
3.2. Las funciones de Chebyshev 63
Esto no basta para probar el teorema de los números primos, pero sí para
obtener un hecho notable que fue conjeturado por Joseph Bertrand en 1845 y
demostrado por Chebyshev en 1852. La prueba clásica de Euclides sobre la
existencia de infinitos primos consiste esencialmente en observar que, dado un
número cualquiera n, tiene que existir un primo n < p ≤ n! + 1, pero esta cota
del menor primo mayor que n es exageradamente elevada y puede mejorarse
mucho:
Teorema 3.9 (Postulado de Bertrand) Para cada número natural n ≥ 1,
existe al menos un primo n < p ≤ 2n. Equivalentemente, la sucesión de los
primos cumple pn+1 ≤ 2pn .
Demostración: Si
(2n)! Y
N= = pkp , (3.5)
(n!)2
p≤2n
Supongamos que existe un n > 29 = 512 para el que no existe ningún primo
p en las condiciones del enunciado. Sea N el dado por (3.5) y sea p | N , de
modo que p ≤ 2n, luego por hipótesis p ≤ n.
Supongamos que 2n/3 < p ≤ n, con lo que 2p ≤ 2n < 3p, luego también
p2 > 4n2 /9 > 2n. Entonces
kp = E[2n/p] − 2E[n/p] = 2 − 2 = 0,
lo que contradice que p | N . Por lo tanto, tiene que ser p ≤ 2n/3, luego, por el
teorema 3.8, tenemos que
P P 4
log p ≤ log p = ϑ(2n/3) ≤ n log 2.
p|N p≤2n/3 3
luego √
P P P
log N ≤ log p + kp log p ≤ log p + 2n log 2n
kp =1 kp ≥2 p|N
4 √
≤ n log n + 2n log 2n.
3
Por otra parte, N es el mayor de los 2n términos del desarrollo de (1 + 1)2n ,
luego 22n ≤ 2nN , luego
4 √
2n log 2 ≤ log 2n + log N ≤ n log 2 + (1 + 2n) log 2n,
3
de donde √
2n log 2 ≤ 3(1 + 2n) log 2n. (3.7)
Estamos suponiendo que n > 512, luego
log(n/512)
α= >0
10 log 2
y además
n
210(1+α) = 210 210α = 210 = 2n.
29
Así (3.7) se convierte en 210(1+α) ≤ 30(25+5α + 1)(1 + α), con lo que
25α ≤ 30 · 2−5 (1 + 2−5−5α )(1 + α) < (1 − 2−5 )(1 + 2−5 )(1 + α) < 1 + α,
pero
25α = e5α log 2 > 1 + 5α log 2 > 1 + α.
Tenemos así una contradicción, que prueba el postulado de Bertrand para todo
n > 512. Para primos menores también se cumple, pues basta considerar los
primos:
Cada uno es menor que el doble del anterior, por lo que entre n y 2n hay siempre
uno de ellos, para todo n < 631.
1.4
1.2
La gráfica representa la función pn+1 /pn ,
1.0
y muestra que la cota 2 que proporciona el 0.8
postulado de Bertrand no es muy fina. 0.6
0.4
Seguidamente obtenemos más estimacio- 0.2
nes de las funciones de Chebyshev:
0 200 400 600 800 1000
Por consiguiente,
log d(n) = log 2r = r log 2 = π(pr ) log 2
log n log 2 log n log 2 (1 − ) log n log 2
≥ ≥ ≥
log pr log log n − log A log log n
para todo n suficientemente grande, ya que la última desigualdad equivale a
(1 + ) log A
log log n ≥ − ,
(notemos que log A < 0, pues ϑ(x) < x).
66 Capítulo 3. La distribución de los números primos
1
P
Ahora bien, la serie de la izquierda no es más que n, donde A es el
n∈A
conjunto de los números naturales no divisibles entre primos menores que N .
Si llamamos P al producto de los primos menores o iguales que N tenemos que
nP + 1 ∈ A para todo n, luego
∞
X 1
< +∞,
n=1
nP +1
pero por otra parte es claro que esta serie diverge, por comparación con la serie
∞
1
P
n . Por lo tanto la serie de partida no puede converger.
n=1
Ahora necesitamos un resultado técnico:
Teorema 3.13 Si h > 0, entonces
logh (x/n) = O(x).
P
n≤x
Por lo tanto
E[x] ∞
x x
logh u logh u
X Z Z Z
logh (x/n) ≤ logh (x/n) dt = x du < x du.
n=2 1 1 u2 1 u2
En particular,
P P
log n = E[x] log x − log(x/n) = E[x] log x + O(x) = x log x + O(x).
n≤x n≤x
E[x/pm ] log p =
P P P P
log n = j(E[x], p) log p = E[x/n]Λ(n).
n≤x n≤x pm ≤x n≤x
luego Xx X
Λ(n) = log n + O(x) = x log x + O(x).
n
n≤x n≤x
Teorema 3.14
X Λ(n)
= log x + O(1).
n
n≤x
La serie converge por [An 8.62], pues puede mayorarse por la suma de 1/n3/2 .
log p
Ahora aplicamos 2.20 con cp = p y f (x) = 1/ log x. Así
X log p
C(x) = = log x + τ (x),
p
p≤x
donde τ (x) está acotada por el teorema anterior. Así obtenemos que
X1 Z x Z x Z x
C(x) C(t) τ (x) dt τ (t)
= + 2 dt = 1 + log x + + 2 dt.
p log x 2 t log t 2 t log t 2 t log t
p≤x
68 Capítulo 3. La distribución de los números primos
donde Z +∞
τ (t)
M = 1 − log log 2 + dt
2 t log2 t
y Z +∞
τ (x) τ (t)
R(x) = − dt = O(1/ log x).
log x x t log2 t
En definitiva: X1
= log log x + M + O(1/ log x).
p
p≤x
luego también
log π(x) log log x
+ −→ 1.
log x log x
El segundo sumando tiende a 0, con lo que
log π(x)
lím = 1.
x→+∞ log x
n log n
lím = 1,
n pn
es decir, pn ∼ n log n.
La diferencia ostensible entre las gráficas de pn /n y log n se debe a que, como
hemos señalado en la introducción, la convergencia a 1 de π(x)xlog x es muy lenta.
El teorema siguiente es un refinamiento asintótico del postulado de Bertrand:
Teorema 3.19 Para todo > 0, existe un x0 > 0 tal que, para todo x ≥ x0
existe un número primo x < p ≤ (1 + )x.
luego
π((1 + )x)
lím = 1 + .
x→+∞ π(x)
Por lo tanto, para todo x suficientemente grande, π(x) < π((1 + )x), luego
existe un primo x < p ≤ (1 + )x.
He aquí una consecuencia curiosa:
70 Capítulo 3. La distribución de los números primos
La tabla siguiente muestra el menor primo que empieza por 1, el menor primo
que empieza por 2, el menor primo que empieza por 3, y así sucesivamente hasta
el menor primo que empieza por 100.
11 2 3 41 5 61 7 83 97 101
11 127 13 149 151 163 17 181 19 2 003
211 223 23 241 251 263 271 281 29 307
31 3 203 331 347 353 367 37 383 397 401
41 421 43 443 457 461 47 487 491 503
5 101 521 53 541 557 563 571 587 59 601
61 6 203 631 641 653 661 67 683 691 701
71 727 73 743 751 761 773 787 79 809
811 821 83 8 419 853 863 877 881 89 907
911 929 937 941 953 967 97 983 991 1009
Teorema 3.20 Para cada número natural N > 0, existe un primo cuyas pri-
meras cifras decimales son las de N .
ω(n) = r, Ω(n) = e1 + · · · + er .
Demostración: Llamemos
P P P P
S1 = ω(n) = 1= E[x/p],
n≤x n≤x p|n p≤x
donde hay que entender que pm | x significa que hay un sumando para cada
potencia pm que divide a x, con m ≥ 1. Por consiguiente
E[x/pm+1 ].
P
S2 − S1 =
pm+1 ≤x
Así pues,
X 1
= M 0 − M + o(1)
pm+1
pm+1 ≤x
Teorema 3.23 Las funciones ω(n) y Ω(n) tienen orden normal log log n. Más
precisamente, para todo δ > 0, casi todos los números naturales cumplen
1
|ω(n) − log log n| ≤ (log log n) 2 +δ ,
luego, dado > 0, existe un n0 a partir del cual el miembro derecho es < , y
así, el conjunto de números que no cumple
ω(n)
log log n − 1 <
y (3.10) es o(x), y los que sí que cumplen estas dos condiciones cumplen también
|ω(n) − log log n| ≤ |ω(n) − log log x + log log x − log log n|
1 δ 1 δ 1
≤ (log log x) 2 + 2 + 1 ≤ (log log n + 1) 2 + 2 + 1 ≤ (log log n) 2 +δ ,
donde para que se cumpla la última desigualdad tenemos que exigir también
que n sea suficientemente grande.
Así pues, el conjunto de los números n ≤ x que no cumplen la desigualdad
del enunciado es o(x), que es lo que hay que probar. Notemos también que esta
reducción vale igualmente para Ω.
Consideremos todos los pares ordenados (p, q) de divisores primos de n tales
que p 6= q. En total hay
ω(n)2 − ω(n) =
P P P
1= 1− 1,
pq|n pq|n p2 |n
p6=q
donde en el penúltimo sumatorio hay que entender que es sobre los pares orde-
nados (p, q) tales que pq | n. Por lo tanto,
X X X x X x
ω(n)2 − ω(n) = E[ ] − E[ 2 ],
pq 2
p
n≤x n≤x pq≤x p ≤x
donde el penúltimo sumatorio es sobre los pares (p, q) tales que pq ≤ x. Ahora
X x X x X 1
E[ 2
]≤ 2
≤x = O(x).
p p p
p2
p2 ≤x 2
p ≤x
74 Capítulo 3. La distribución de los números primos
pues al quitar las partes enteras estamos sumando a lo sumo 2E[x] números
(porque cada pq = qp con p 6= q cuenta dos veces) menores que 1, luego el
término añadido es O(x). Usando el teorema anterior llegamos a que
X X 1
ω(n)2 = x + O(x log log x)
pq
n≤x pq≤x
pues al elevar al cuadrado el primer√ sumatorio, obtenemos una suma sobre todos
los pares (p, q) tales que p, q ≤ x, que están entre los que cumplen pq ≤ x,
y al elevar al cuadrado el último sumatorio obtenemos la suma sobre todos los
pares (p, q) tales que p, q ≤ x, los cuales incluyen a los que cumplen pq ≤ x. Los
dos extremos de la última fórmula son
luego X
ω(n)2 = x log2 log x + O(x log log x).
n≤x
A su vez,
(ω(n) − log log x)2 = ω(x)2 − 2 log log x ω(x) + E[x] log2 log x
P P P
n≤x n≤x n≤x
= x log2 log x + O(x log log x) − 2 log log x(x log log x + O(x))+
(x + O(1)) log2 log x
= x log2 log x − 2x log2 log x + x log2 log x + O(x log log x) = O(x log log x).
Supongamos finalmente que la cantidad de números n ≤ x que no cumplen
la condición (3.9) no fuera o(x). Esto significa que para todo > 0 existen
valores de x arbitrariamente grandes tales que hay más de x números n ≤ x
que cumplen
1
|ω(n) − log log x| > (log log x) 2 +δ ,
luego
(ω(n) − log log x)2 ≥ x(log log x)1+2δ ,
P
n≤x
es o(x), luego también lo es la cantidad de los que cumplen alguna de las dos.
Por el contrario, los n ≤ x que no cumplen ninguna, cumplen también
≤ (log log x)1/2 + (log log x)1/2+δ/2 < (log log x)1/2+δ
si x es suficientemente grande. Así pues, la cantidad de números n ≤ x que no
cumplen esto es o(x), que es lo que había que probar.
Así pues, tanto el orden medio como el orden normal de ω(n) y Ω(n) es
log log n. Para otras funciones aritméticas no tiene por qué darse la igualdad.
Teorema 3.24 La función log d(n) tiene orden normal log 2 log log n.
d(n) = (1 + e1 ) · · · (1 + er ),
y 2 ≤ 1 + ei ≤ 2ei , luego
y a su vez
ω(n) log 2 ≤ log d(n) ≤ Ω(n) log 2,
luego
ω(n) log d(n) Ω(n)
≤ ≤ .
log log n log 2 log log n log log n
Dado > 0, los dos extremos de las desigualdades anteriores distan de 1 menos
que para casi todo n, luego lo mismo le sucede al término central.
76 Capítulo 3. La distribución de los números primos
Explícitamente, esto significa que, para todo > 0, casi todos los números
naturales cumplen
(log n)(1−) log 2 = 2(1−) log log n ≤ d(n) ≤ 2(1+) log log n = (log n)(1+) log 2 .
Cabe señalar que estas cotas son bastante inferiores al valor medio de la función
d(n) que, de acuerdo con el teorema 2.23, es log x + 2γ − 1.
30
20
15
log x + 2γ − 1
10
Por lo tanto,
∞
ζ(2s) X λ(n)
= , para σ > 1.
ζ(s) n=1
ns
3.5. Series de Dirichlet 79
1
log ζ(s) = log + O(s − 1). (3.13)
s−1
Teorema 3.25 Se cumple que
P
M =γ+ (log(1 − 1/p) + 1/p) .
p
y la serie
X 1 X 1
<
p
2p(p − 1) n
2n(n − 1)
Por lo tanto:
Z +∞ Z 2
R(t) log log t + M
g(δ) + γ − M + log δ = δ dt − δ dt.
2 t1+δ 1 t1+δ
√
1/ δ
Ahora, llamando T = e y teniendo en cuenta que |R(t)| ≤ A/ log t, tenemos
que Z +∞ Z T Z +∞
R(t) dt dt
δ
1+δ
dt ≤ δA 1+δ log t
+ δA 1+δ log t
2 t 2 t T t
Z T Z +∞
δA dt δA dt δA log T A
≤ δ + 1+δ
= δ + δ
2 log 2 2 t log T T t 2 log 2 T log T
δA √ √
< √ + A δ = 2A δ.
δ
Por otra parte,
Z 2 Z 2
log log t + M | log log t| + M
δ
1+δ
dt ≤ δ dt
1 t 1 t
(notemos que la integral converge, pues se puede calcular haciendo el cambio
u = log t). Concluimos que
lím g(δ) + log δ = M − γ.
δ→0
Pero por (3.13) también tenemos que log ζ(1 + δ) + log δ → 0. Por consiguiente,
F (δ) = g(δ) + log δ − (log ζ(1 + δ) + log δ) → M − γ,
luego F (0) = 0, es decir,
X 1
1
log 1 − + =M −γ
p
p p
Por consiguiente:
X
1
1
X
1
lím log log x + γ + log 1 − + − = 0,
x→∞ p p p
p≤x p≤x
o también:
X 1
lím log log x + γ + log 1 − = 0,
x→∞ p
p≤x
y aplicando la exponencial,
Y 1
lím log x eγ 1− = 1.
x→+∞ p
p≤x
σ(p)/p = 1 − 1/p ≈ 1,
σ(n)
,
n log log n
y en ella podemos ver que está acotada.
Teorema 3.27
σ(n) φ(n) log log n
lím = eγ , lím = e−γ ,
n n log log n n n
donde σ es la función suma de divisores, φ es la función de Euler y γ es la
constante de Euler.
3.5. Series de Dirichlet 83
pa+1 − 1
σ(pa )φ(pa ) = (p − 1)pa−1 = p2a − pa−1 < (pa )2 .
p−1
Llamemos
σ(n) φ(n)eγ log log n
f1 (n) = , f2 (n) = .
neγ log log n n
Hay que probar que lím f1 = 1 y lím f2 = 1, para lo cual basta encontrar
n n
funciones F1 (t), F2 (t) que tiendan a 1 en +∞ de modo que
1
f1 (n) ≤ , f2 (n) ≥ F1 (log n)
F1 (log n)
para todo n ≥ 3, y
1
f1 (nj ) ≥ F2 (j), f2 (nj ) ≤ ,
F2 (j)
Y 1
γ
lím e log t 1− = 1,
t→+∞ p
p≤t
Tomemos ahora
pj ,
Q
nj =
p≤ej
de modo que log nj = jϑ(ej ) ≤ Ajej , para cierta constante A > 0, por 3.10.
Por lo tanto
log log nj ≤ A0 + j + log j.
Ahora:
Y 1
Y 1
1
1− > 1− = .
pj+1 p
pj+1 ζ(j + 1)
p≤ej
Por lo tanto,
e−γ Y 1
≥ = F2 (j),
ζ(j + 1)(A0 + j + log j) j
1 − p−1
p≤e
e−γ Y 1
lím = 1,
j→+∞ log ej
j
1 − p−1
p≤e
luego
j
lím F2 (j) = lím = 1.
j→+∞ j→+∞ A0 + j + log j
Vemos que algunos valores superan el límite superior. No hay nada de ex-
traño en ello, pero podemos afirmar que la sucesión de todos los valores por
encima del límite superior (si es infinita) converge a dicho límite superior. Los
puntos en el rango mostrado en la figura que sobrepasan el límite superior son
60, 72, 84, 120, 180, 240, 360, 720, 840, 2 520, 5 040.
El caso es que no se conoce ninguno mayor que 7! = 5 040.
Por el desarrollo en serie de Dirichlet de log ζ(s) tenemos que, para σ > 1,
∞ ∞
X Λ(n) −s log n X Λ(n)
log |ζ(s)| = Re e = cos(τ log n).
n=2
log n n=2
log n nσ
Por lo tanto,
log |ζ(σ)3 ζ(σ + τ i)4 ζ(σ + 2τ i)| = 3 log |ζ(σ)| + 4 log |ζ(σ + τ i)| + log |ζ(σ + 2τ i)|
∞
X Λ(n)
= (3 + 4 log(τ log n) + cos(2τ log n)) ≥ 0.
n=2
log n nσ
Por lo tanto,
ζ(σ + τ i) 4
3
((σ − 1)ζ(σ))
|ζ(σ + 2τ i)| ≥ 1 ,
σ−1 σ−1
para σ > 1 y cualquier τ . Si fuera ζ(1 + τ i) = 0, entonces
ζ(σ + τ i) ζ(1 + τ i + σ) − ζ(1 + τ i)
lím = lím = ζ 0 (1 + τ i)
σ→1 σ−1 σ→1 σ−1
y el miembro izquierdo de la desigualdad precedente tendería a
y al comparar esta serie con Φ(2s) concluimos que de hecho converge (a una
función holomorfa) en el semiplano σ > 1/2. Por lo tanto, Φ(s) se prolonga
analíticamente hasta la función meromorfa
ζ 0 (s) X log p
Φ(s) = − −
ζ(s) p
p2s − ps
en el semiplano σ > 1/2. Sus polos son los de la derivada logarítmica ζ 0 (s)/ζ(s),
que por [VC 3.20] son todos simples, y se corresponden con los ceros y polos de
ζ(s). Como ζ(s) tiene un polo simple en s = 1, resulta que Φ(s) tiene también un
polo simple en s = 1 con residuo 1. Los polos restantes de Φ(s) se corresponden
con los ceros de ζ(s), luego por el teorema anterior no hay ninguno más en la
recta σ = 1.
Ahora necesitamos un resultado técnico que traduce propiedades analíticas
de una serie de Dirichlet en una propiedad de la sucesión de sus coeficientes:
Teorema 3.31 Sea
∞
X an
f (s) =
n=1
ns
una serie de Dirichlet con coeficientes an ≥ 0 y convergente en un semiplano
σ > σ0 ≥ 0. Sea P
φ(x) = an .
log n≤x
Entonces Z +∞
f (s) = s φ(x)e−sx dx, para σ > σ0 .
0
k
P
donde h(k) = an . Sea σ > σ1 > σ0 > 0. Entonces
n=1
k Z log(k+1)
X an h(k) 1
s
=s φ(x)e−sx dx + σ1 (k + 1)s−σ1
. (3.14)
n=1
n 0 (k + 1)
Claramente
k
h(k) X an 1 1 k
0≤ ≤ ≤ f (σ1 ), (k + 1)s−σ1 = (k + 1)σ−σ1 → 0,
(k + 1)σ1 n=1
n σ1
Teorema 3.32 Sea φ : [0, +∞[ −→ R una función acotada, integrable en cada
intervalo [0, x], y supongamos que la función
Z +∞
g(s) = φ(x)e−sx dx,
0
Demostración: Llamemos
Z t
gt (s) = φ(x)e−sx dx.
0
Por el teorema [VC 1.23] tenemos2 que gt es una función entera. Claramente
Z +∞
φ(x) dx = lím gt (0),
0 t→+∞
las dos variables, cosa que aquí no tenemos garantizada, porque φ no tiene por qué ser continua,
pero la prueba de [VC 1.23] se basa en [An 8.57], y este teorema admite que el integrando
tenga un factor integrable acotado no necesariamente continuo, que dependa únicamente de x,
por lo que la prueba de [VC 1.23] se aplica igualmente a este caso.
3.6. Prueba del teorema de los números primos 89
s2
h(s) = (g(s) − gt (s))est 1 + 2 .
R
Así:
s2
Z
1 st ds
g(0) − gt (0) = h(0) = (g(s) − gt (s))e 1+ 2 .
2πi γ R s
Si C es una cota de φ, en el semiplano σ > 0 se cumple que
+∞ +∞
Ce−σt
Z Z
−sx
|e−sx | dx =
|g(s) − gt (s)| = φ(x)e dx ≤ C
,
t t σ
s2 1
st 2σ
e 1+ 2 = eσt 2 ,
R s R
donde hemos usado que |1 + s2 /|s|2 | = 2σ/|s|, como se comprueba sin más que
sustituir s = a + bi. Por lo tanto,
Z
1 h(s) 1 2C C
2πi ds ≤ 2
πR = .
γ1 s 2π R R
coincide con la integral de la misma función sobre el arco γ3 (u) = Reui , para
u ∈ [π/2, 3π/2], pues la integral sobre γ1 ∪ γ2 , al igual que la integral sobre
γ1 ∪ γ3 , es el valor del integrando en 0. Sobre γ3∗ tenemos igual que antes que
Z t Z t t
Ce−σt
Z
|gt (s)| = φ(x)e−sx dx ≤ C |e−sx | dx ≤ C e−σx dx =
,
0 0 −∞ σ
Tomamos una cota C 0 de g(s)(1+s2 /R2 )/s sobre γ2∗ y, por compacidad, podemos
tomar t0 > 0 tal que si t > x0 entonces,
2π
|est | <
2L(γ2 )C 0
para todo s ∈ γ2∗ . Esto hace que la integral sea menor que /2 y, en definitiva,
que |g(0) − gt (0)| < , para todo t ≥ t0 .
Ahora ya podemos probar un teorema cuya conclusión no tiene nada que ver
con series de Dirichlet: bajo ciertas hipótesis sobre una serie de Dirichlet, nos
asegura que una función que por hipótesis está acotada, de hecho converge a 1:
Teorema 3.33 Consideremos una serie de Dirichlet
∞
X an
f (s) =
n=1
ns
con coeficientes an ≥ 0 tales que exista una constante A > 0 de modo la función
de las sumas finitas P
θ(x) = an
n≤x
Demostración: La función Φ(x) del teorema 3.30 cumple las hipótesis del
teorema anterior. Notemos que la sucesión de sumas finitas de sus coeficientes
es precisamente ϑ(x), y la hipótesis ϑ(x) ≤ Ax es el teorema 3.8. La conclusión
es precisamente que ϑ(x) ∼ x.
Por consiguiente, todas las afirmaciones del teorema 3.7 son ciertas, y queda
demostrado el teorema de los números primos.
Probamos ahora una variante de 3.33 para sucesiones de números complejos
cualesquiera que usaremos en la sección siguiente.
92 Capítulo 3. La distribución de los números primos
Teorema 3.35 Sea f (s) una serie de Dirichlet en las condiciones del teorema
anterior y sea g(s) una serie de Dirichlet con coeficientes bn ∈ C convergente
en el semiplano σ > 1 y prolongable analíticamente a la recta σ = 1 salvo quizá
a s = 1, donde tiene un polo simple con residuo α (entendiendo que α = 0 si
P Supongamos que existe una constante C tal que |bn | ≤ Can y
no hay tal polo).
sea ψ(x) = bn . Entonces
n≤x
ψ(x)
lím = α.
x→+∞ x
Demostración: Supongamos primero que los coeficientes de g son reales
(y por lo tanto α también). Podemos tomar C > α. Entonces la función
f ∗ = (Cf + g)/(C + α) cumple las mismas hipótesis que f , luego el teorema
anterior nos da que
o equivalentemente,
ψf (x) ψg (x)
C lím + lím = C + α,
x→+∞ x x→+∞ x
donde el segundo límite existe porque existe el primero y el de la suma. Como
el primer límite es 1, la conclusión es inmediata.
En el caso general consideramos la serie
∞
X b̄n
ḡ(s) = s
.
n=1
n
Vemos que están repartidos de forma casi uniforme. Sólo podrían haber
estado mejor repartidos si uno de los acabados en 3 hubiera acabado en 1 o
en 9. Si consideramos intervalos de primos mayores, el número de primos con
cada terminación se aproxima más al 25%:
Hasta Total 1 3 7 9
100 23 5 (22%) 7 (30%) 6 (26%) 5 (22%)
1 000 166 40 (24%) 42 (25%) 46 (28%) 38 (23%)
10 000 1 227 306 (24.9%) 310 (25.3%) 308 (25.1%) 303 (24.7%)
100 000 9 590 2 387 (24.9%) 2 402 (25.1%) 2 411 (25.1%) 2 390 (24.9%)
En esta sección usaremos las técnicas que hemos empleado en la prueba del
teorema de los números primos para demostrar algo más fuerte:
y al integrar queda αψ = 0.
Como la dimensión de CG es |G|, concluimos que los caracteres forman una
base, con lo que toda función f ∈ CG se expresa como combinación lineal de los
caracteres de G.
Sea A un conjunto de primos no vacío y sea Ax = {p ∈ A | p ≤ x}. Diremos
que una aplicación λ : A −→ G está equidistribuida si para todo carácter χ de G
se cumple P
χ(λ(p)) Z
p∈Ax
lím = χ(g) dg. (3.15)
x→+∞ |Ax | G
πk (x) 1
lím = ,
x→+∞ π(x) − c φ(m)
donde c es el número de primos que dividen a m, pero es claro que dicha c puede
eliminarse.
Así pues, basta probar que todo carácter no principal χ módulo m cumple
que el límite de (3.15) es igual a 0.
3.7. Primos en progresiones aritméticas 95
Por el teorema [VC 5.27] esta función es holomorfa sobre toda la recta σ = 1.
Si separamos los términos con n ≥ 2 obtenemos una serie mayorada en módulo
por
X log p
,
pns
p,n≥2
y en la demostración del teorema 3.30 hemos probado que esta serie converge
en el semiplano σ > 1/2. Por consiguiente, la serie restante,
X log p χ(p)
p
ps
La conclusión es que X
log p χ(p) = o(x),
p≤x
donde los números B2k son los números de Bernoulli definidos en [VC 4.9] y
determinados recurrentemente en [VC 4.34]. Tras [VC 5.31] razonamos que esta
función holomorfa admite una prolongación analítica a una función meromorfa
en C con un único polo en el punto s = 1, que es un polo simple con residuo 1.
El teorema [VC 5.32] nos da el valor de esta prolongación en 0 y sobre los
números negativos:
Bn+1
ζ(−n) = − .
n+1
En particular la función dseta se anula en los enteros negativos pares, mientras
que ζ(0) = −1/2.
97
98 Capítulo 4. La función dseta de Riemann I
= π(2π)s−1 .
cos(πs/2)Π(s − 1)
-6
-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6
1 2
-1 1 2 3 4
-1 20 40 60 80 100
-2
-2
-3
100 Capítulo 4. La función dseta de Riemann I
Las gráficas muestran que la curva pasa muchas veces por el punto 0, es
decir, que la función dseta tiene ceros no triviales. Más adelante probaremos
que esto es realmente así.
x
t − E[t]
Z
E[x] s s
= s + − −s dt
x s − 1 (s − 1)xs−1 1 ts+1
x
x − E[x] t − E[t]
Z
s 1
= − s
− −s dt.
s−1 x (s − 1)xs−1 1 ts+1
Ahora bien, la última integral, para un x fijo, define una función entera,
luego por el principio de prolongación analítica, la igualdad
Z x
X 1 s x − E[x] 1 t − E[t]
= − − − s dt (4.1)
ns s−1 xs (s − 1)xs−1 1 ts+1
n≤x
Como esta función es integrable en [1, +∞[, el teorema [VC 1.24] nos da que
la integral de la primera igualdad del enunciado es una función holomorfa en
el semiplano σ > 0, luego por el principio de prolongación analítica no sólo es
válida cuando σ > 1, sino también cuando σ > 0. La segunda fórmula se obtiene
restando (4.1) de la primera.
6 10
4
5
2
0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0
-2
-5
-4
ζ(σ + 0i), k = 1
4
5
0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0
-2
-5
-4
ζ(σ + 0i), k = 2
3.0
2.5 15
2.0 10
5
1.5
Re ζ(σ + 6πi), k = 2
15
0.5 1.0 1.5 2.0 2.5 3.0
-0.2 10
5
-0.4
Im ζ(σ + 6πi), k = 2
Enseguida daremos una forma alternativa de acotar el error con más preci-
sión, pero antes vamos a extraer las primeras consecuencias teóricas del teorema
4.1:
4.1. Aproximación de la función dseta 103
15
3
10
2
5
1
5 10 15 20
-5
5 10 15 20
-1 -10
Re ζ(1/2 + τ I), k = 10
1 10
5
5 10 15 20
-1 5 10 15 20
-5
-2
-10
-3
-15
Im ζ(1/2 + τ I), k = 10
Teorema 4.2 Si 0 < s < 1, entonces ζ(s) < 0. Por consiguiente, los únicos
ceros reales de la función dseta son sus ceros triviales.
Aunque no nos va a ser necesario más adelante, vamos mejorar la cota del
error dada por la fórmula (4.2). Para ello calcularemos las colas de la serie de
Dirichlet por un método similar al empleado en el teorema [VC 3.18] para sumar
series mediante el teorema de los residuos.
|z −s | = |e−s log z | = |e−(σ+τ i)(log |z|+i arg z) | = e−σ log |z|+τ arg z = |z|−σ eτ arg z .
Por lo tanto,
Z
1 cot πz π −σ |τ |π/2
dz ≤ 2 (r + 1/2)M (r + 1/2) e
2i z s
Cr
π
M e|τ |π/2 (r + 1/2)1−σ ,
=
2
que ciertamente converge a 0 si σ > 1. Concluimos que
1 x+∞i cot πz
X 1 Z
ζ(s) − = − dz.
n<x
ns 2i x−∞i zs
Hemos probado que la integral existe como límite simultáneo de sus dos
extremos de integración, pero el límite es también finito si sólo hacemos tender
a infinito uno de los extremos, porque |z −s cot πz| ≤ M e|τ |π/2 |z|−σ , luego
Z x+T i Z T
cot πz i cot πz
dz = dy
x zs 0 zs
y el módulo del integrando está acotado por cy −σ (con c constante) que (siempre
para σ > 1) es integrable en [1, +∞[, luego la función es integrable en [1, +∞[
y, por continuidad, también en [0, +∞[. Por consiguiente podemos separar:
1 x 1 x+∞i cot πz
X 1 Z Z
cot πz
ζ(s) − = − dz − dz.
n<x
ns 2i x−∞i z s 2i x zs
e igualmente
Z x+∞i
1 1
s
dz = .
x z (s − 1)xs−1
Por lo tanto,
1 x cot πz − i
X 1 Z
1
ζ(s) − s
− s−1
= − dz
n<x
n (s − 1)x 2i x−∞i zs
Z x+∞i
1 cot πz + i
− dz, (4.3)
2i x zs
106 Capítulo 4. La función dseta de Riemann I
|z −s | = |z|−σ eτ arg z ≤ x−σ e|τ || arg z| ≤ x−σ e|τ | arctan(|y|/x) ≤ x−σ e|τ ||y|/x .
1 Es fácil ver que, si x ≥ 0, se cumple que arctan x ≤ x. Basta tener en cuenta que la
luego Z x
x−σ
1 cot πz − i
dz ≤
2π − |τ |/x .
2i
x−∞i zs
La segunda integral de (4.3) se acota análogamente y así llegamos a la fórmula
del enunciado.
Las gráficas siguientes muestran las acotaciones de la parte real de ζ(σ +6πi)
con k = 3 y ζ(1/2 + τ i) con k = 10, respectivamente, que proporcionan (4.2) y
el teorema anterior. En ambos casos, las más ajustadas son las correspondientes
al teorema anterior.
3.5
15
3.0
10
2.5
2.0 5
1.5
5 10 15 20
1.0
-5
0.5
-10
0.0 0.5 1.0 1.5 2.0 2.5 3.0
1 2 · 4.5−1/2
≤ −1/2
+ ≈ 0.5627.
6π4.5 2π − 6π/4.5
Por otra parte,
4
X 1
Re σ+6πi
= 1 + 2−σ cos(6π log 2) + 3−σ cos(6π log 3) + 4−σ cos(6π log 4)
n=1
n
uniformemente en τ .
En otras palabras, tomando σ0 suficientemente grande podemos asegurar
que en todo el semiplano σ > σ0 la función ζ(s) se aproxime a 1 tanto como
queramos. Podríamos decir que en los semiplanos σ > σ0 con σ0 grande es
donde “menos interesante” es la función dseta.
En particular, ζ(s) está acotada en cualquier semiplano σ > σ0 , con σ0 > 1.
Más precisamente, un hecho elemental es que, para todo s en el semiplano σ > 1,
se cumple que
∞ ∞
X 1 X 1
|ζ(s)| = ≤ = ζ(σ),
ns n=1 nσ
n=1
(uα1 , . . . , uαn ),
Por lo tanto,
N N ∞
X 1 X 1 X 1
σ
cos(t log n) ≥ cos(2π/q) σ
> cos(2π/q)ζ(σ) − σ
.
n=1
n n=1
n n
n=N +1
110 Capítulo 4. La función dseta de Riemann I
Por consiguiente
∞
X 1
|ζ(σ + it)| ≥ cos(2π/q)ζ(σ) − 2 .
nσ
n=N +1
Concluimos que
cos(2π/q) − 2N 1−σ ≥ 1 − .
Así pues, ahora podemos asegurar que ζ(s) no está acotada, por ejemplo,
sobre el cuadrante σ > 1, τ > 1, pues el supremo de |ζ(s)| sobre cada semirrecta
σ = σ0 en dicho cuadrante sigue siendo ζ(σ0 ), que tiende a +∞ cuando σ0 tiende
a 1.
Esto nos plantea el problema opuesto, es decir, estimar el orden de creci-
miento de ζ(s) en un cuadrante como σ > 1, τ ≥ 2 (la cota sobre τ es irrelevante,
con tal de que excluya un entorno del polo de la función dseta). El teorema 4.1
nos da la solución:
|s − 1| ≥ Im(s − 1) = τ,
se reduce a Z +∞
X 1 1 1 dt
|ζ(s)| ≤ σ
+ σ−1 + σ + |s| σ+1
n τx x x t
n≤x
X 1 1 1 σ+τ 1
< + σ−1 + σ + .
nσ τx x σ xσ
n≤x
En particular, si σ ≥ 1, τ ≥ 1, x ≥ 1
X 1 Z x
1 1 1+τ 1 τ τ
|ζ(s)| < + + + ≤4+ dt + = 4 + log x + ,
n τ x x 1 t x x
n≤x
x1−δ 3τ
≤ + x1−δ + δ .
1−δ δx
En particular, para x = τ , queda
τ 1−δ 3
|ζ(s)| ≤ + τ 1−δ + τ 1−δ = O(τ 1−δ ),
1−δ δ
donde la constante de O depende de δ.
1
χ(s) = π(2π)s−1 .
cos(πs/2)Π(s − 1)
En particular,
σ−1/2
2π
|χ(s)| = (1 + O(1/τ )).
τ
∞ ∞
!
1 X cr+1 −r X cr+2 −r 1
(c − ) τ + τ
2 r=0 r + 1 r=0
r+2 τ
pero τ −1 varía en ]0, 1/c[ y podemos ver las dos series como series de potencias
actuando sobre τ −1 , luego están acotadas en cualquier intervalo ligeramente
menor, de modo que la expresión anterior está acotada por A/τ , para todo
τ > c + 1, por ejemplo, y usando que la expresión inicial tiene que estar acotada
en el compacto [1 − c, 1 + c] × [δ, c + 1], podemos aumentar la constante para
que valga en todo el conjunto S.
4.2. El crecimiento de la función dseta 113
1 2 2
= iπ(σ+iτ )/2 = iσπ/2−πτ /2
cos(πs/2) e + e−iπ(σ+iτ )/2 e + e−iπσ/2+πτ /2
2eiπσ/2−πτ /2
= = 2eiπσ/2−πτ /2 (1 + O(1/τ )),
eiπσ−πτ + 1
pues
τ e−πτ
1 0
τ − 1 = → = 0.
eiπσ−πτ +1 |eiπσ−πτ + 1| 1
Notemos que el denominador admite una cota inferior independiente de σ, por
lo que la constante en O(1/τ ) es independiente de σ. En total,
√ π 1 π
χ(s) = π(2π)s−1 2πτ σ+iτ −1/2 ei 2 (σ− 2 )− 2 τ −iτ 2eiπσ/2−πτ /2 (1 + O(1/τ ))
σ+iτ −1/2
2π
= ei(τ +π/4) (1 + O(1/τ )).
τ
1
≤ (2π)−δ−1/2 ζ(1 + δ)τ 1/2−σ O(1) = O(τ 2 −σ ).
Demostración: Por [AA 2.18] sabemos que log |f (s)| es una función harmó-
nica donde f no se anula (y tiende a −∞ en los ceros de f ). También lo es k(σ)τ ,
donde
σ2 − σ σ − σ1
k(σ) = c1 + c2 ,
σ2 − σ1 σ2 − σ1
(pues k(σ)τ = aστ + bτ ), luego, dado > 0, también lo es
2c1 log τ
g(σ1 + iτ ) ≤ log A1 − 2c1 log τ − τ = log A1 − τ ( + ),
τ
y el miembro derecho tiende a −∞ cuando τ tiende a +∞.
Así pues, g(s) está acotada sobre la semirrecta {σ1 } × [τ0 , +∞[, y lo mismo
sucede con la semirrecta {σ2 } × [τ0 , +∞[. Sea M una cota común que sea
también una cota superior de log |f (s)| en el intervalo [σ1 , σ2 ] × {τ0 }, con lo que
también lo es de g.
Por otra parte,
(c − k0 ) log τ
= log A − τ ( + ),
τ
donde k0 es el mínimo de k(σ) en el intervalo [σ1 , σ2 ], y nuevamente el miembro
derecho tiende a −∞, uniformemente en σ, luego podemos encontrar un T > τ0
tal que si Im s = T se cumple que g(s) ≤ M en el segmento [σ1 , σ2 ] × {T }.
Así, la función g está acotada por M en toda la frontera del rectángulo
[σ1 , σ2 ] × [τ0 , T ], luego por el principio del máximo [AA 2.5] (aplicado al rectán-
gulo menos un disco alrededor de cada cero de f en el que g tome valores
menores que M ) concluimos que g(s) ≤ M en todo el rectángulo y, como T
es arbitrariamente grande, de hecho esto vale para toda la semi banda S. En
particular,
luego
|f (σ + iτ )| ≤ M τ k(σ) eτ ,
y, como M no depende de , concluimos que |f (σ+iτ )| ≤ M τ k(σ) . Esto es lo que
había que probar, pues si σ = (1 − λ)σ1 + λσ2 , entonces k(σ) = (1 − λ)c1 + λc2 .
Como consecuencia:
Por [An 4.19] sabemos que toda función convexa es continua, luego tenemos
que la función de Lindelöf es continua. Más aún, por [An 4.18] sabemos que µ
está por encima de su recta tangente en cualquier punto donde sea derivable.
En resumen, ahora podemos afirmar que la
gráfica de µ(σ) es como indica la figura. Cuando @µ(σ)
σ ≤ 0 se cumple que µ(σ) = 1/2 − σ (lo hemos @
1/2 − σ @
probado para σ < 0, pero por continuidad tam- HH
bién tiene que valer para σ = 0).
@ (1 − σ)/2
@HH
Para σ ≥ 1 es µ(σ) = 0 (y de nuevo el caso de @ HH
σ = 1 lo deducimos por continuidad). 0 1/2 1
En la banda crítica 0 < σ < 1 no sabemos cuánto vale exactamente µ(σ),
pero por la definición de convexidad el valor tiene que estar por debajo del
segmento que une los puntos (0, µ(0)) y (1, µ(1)), luego µ(σ) ≤ (1 − σ)/2, y
por otra parte µ(σ) debe estar por encima de las prolongaciones de las rectas
(tangentes) y = 1/2 − σ e y = 0. En otras palabras, la gráfica de µ(σ) debe
estar en el triángulo que muestra la figura.
En particular, tenemos, por ejemplo, que 0 ≤ µ(1/2) ≤ 1/4. En 6.9 proba-
remos que la cota 1/4 puede ser mejorada.
Hemos visto que, para σ > 1, se cumple que µ(σ) = 0 porque ζ(σ + τ i) está
acotada. En cambio, del hecho de que µ(1) = 0 no podemos deducir que la
función ζ(1 + τ i) esté acotada. De hecho, no lo está:
4.3 La función xi
Recordemos que la función xi es la dada por
|s|
| log(π −s/2 (s − 1))| ≤ log π + | log(s − 1)| = O(|s| log |s|),
2
pues q
log2 |s| + arg2 s
s
| log(s − 1)| arg2 s
= = 1+ ,
log |s| log |s| log2 |s|
y el argumento está acotado. Por la fórmula de Stirling [VC 4.44],
√
s+1 s s
log Π(s/2) = log 2π + log − + µ(s/2) = O(|s| log |s|),
2 2 2
pues
|s + 1| | log(s/2)| 1 1 | log |s| − log 2 + i arg(s/2)|
≤ +
2|s| log |s| 2 |s| log |s|
está claramente acotado y, por otra parte, por [VC 4.43], µ(s/2) = O(|s|−1 ).
Por último, por el teorema 4.7 tenemos que
1
|ζ(s)| ≤ cτ 1/2 = ec+ 2 log |s| = O(|s| log |s|).
lím ζ(s) = 1.
s→+∞
1+
Mξ (r) ≤ ecr log r < ecr ,
donde en el último paso hay que aumentar la constante c. Por lo tanto, el orden
de ξ es ≤ 1. Si fuera menor que 1, existiría una constante tal que Mξ (r) ≤ ecr ,
en contra de la nota precedente.
El teorema de Hadamard [VC 4.27] nos da entonces que existen constantes
A y B tales que
∞
Y s
ξ(s) = eA+Bs 1− es/ρn , (4.5)
n=0
ρ n
donde {ρn }∞n=0 es una enumeración de los ceros de ξ(s) repetidos según su mul-
tiplicidad y ordenados de modo que su módulo sea creciente, donde en principio
no excluimos que la función xi tenga un número finito de ceros o incluso que no
tenga ninguno, en cuyo caso el producto será finito o vacío (en cuyo caso hay
que entender que vale 1). Sin embargo, ahora podemos probar que de hecho
ξ(s) tiene infinitos ceros. Más aún:
Teorema 4.15 La función ξ(s) tiene infinitos ceros y, si {ρn }∞n=0 es una enu-
meración en la que cada uno se repite tantas veces como indica su multiplicidad
y de modo que la sucesión de los módulos sea creciente, entonces la serie
∞
X 1
n=0
|ρn |λ
luego
1 − s es/ρn ≤ ec|s|/|ρn | .
ρn
Por consiguiente, ∞
Y s
s/ρn
1− e ≤ ecα|s| ,
ρ
n=0 n
y también es claro que eA+Bs ≤ ec|s| , luego en total tendríamos que |ξ(s)| ≤ ec|s| ,
en contradicción con la nota posterior al teorema 4.13.
Así pues, el número de ceros es infinito y la serie de los inversos de sus módu-
los diverge. La convergencia de las series con exponentes λ > 1 es consecuencia
directa del teorema [VC 4.25], según el cual el exponente de convergencia de la
sucesión de los ceros tiene que ser ≤ 1 (y el argumento precedente prueba que
de hecho es exactamente 1).
Así pues, hemos probado que la función dseta de Riemann tiene infinitos
ceros no triviales.
Nota En lo sucesivo entenderemos que la variable ρ como índice de una suma o
producto recorre los ceros de la función ξ repetidos tantas veces como indica su
multiplicidad y ordenados de modo que la sucesión de sus módulos sea creciente.
1 Y s
s/2 −(γ/2+1)s
ξ(s) = (4π) e 1− es/ρ .
2 ρ
ρ
ξ 0 (s) X 1
1
=B+ + . (4.6)
ξ(s) ρ
s−ρ ρ
luego
∞
1 Π0 (1/2)
γ X 1 1 γ
− = + − = + log 2 − 1,
2 Π(1/2) 2 2k + 1 2k 2
k=1
por la expresión para log 2 que se sigue de la serie de Taylor del logaritmo.3
1
Por otra parte, según 4.3 tenemos que ζ(s) = s−1 + γ + g(s), con g(1) = 0,
luego ζ 0 (s) = − (s−1)
1 0
2 + g (s), luego
luego
ζ 0 (s)
1
lím + = γ.
s→1 ζ(s) s+1
En total √ √
γ γ
B = log + log 2 − 1 − γ = log 4π − − 1,
π+
2 2
lo que nos da la fórmula del enunciado.
ζ 0 (0) √ 1 Π0 (0)
= 1 + log π − +B
ζ(0) 2 Π(0)
√ γ √ γ
= 1 + log π+ + log 2 π − − 1 = log 2π,
2 2
luego ζ 0 (0) = − 21 log 2π.
Otra consecuencia:
pues las sumas parciales de esta serie son la subsucesión formada por los términos
pares de la serie anterior. Ahora bien,
1 1 ρ + ρ̄ 2 Re ρ 2
+ = = ≤ 2,
ρ ρ̄ |ρ|2 |ρ|2 |ρ|
y la última serie es convergente. Lo mismo vale para los sumandos con 1 − ρ y
1 − ρ̄, luego podemos descomponer la serie en suma de dos series de términos
positivos:
X 2 Re(1 − ρ) X 2 Re ρ
2
+ = −2B.
|1 − ρ| |ρ|2
Im ρ>0 Im ρ>0
y a su vez
∞ ∞
a2n (it)2n = (−1)n a2n t2n ,
P P
Ξ(t) = ξ(1/2 + it) =
n=0 n=0
lo que prueba que, tal y como habíamos afirmado, Ξ(t) es real cuando t ∈ R.
Además Ξ(t) = Ξ(−t).
Observemos que los ceros reales de Ξ(t) son las partes imaginarias de los
ceros de ξ(s) sobre la recta crítica σ = 1/2. He aquí la gráfica de Ξ(t) y una
ampliación en la que se ve claramente que ξ tiene un cero con parte real σ = 1/2
y parte imaginaria 14 < τ < 15.
0.5
0.005
0.4
0.004
0.3 0.003
0.2 0.002
0.001
0.1
13 14 15 16 17 18 19
5 10 15 -0.001
N (T )
lím = +∞,
T →+∞ T
pues
N (T ) 1 T 1 log T
= log + + O( ).
T 2π 2π 2π T
Como Φ(s) = Φ(1 − s) y Φ(s̄) = Φ(s) (esto es cierto para ξ, luego también
para Φ), la integral toma el mismo valor en los cuatro cuadrantes del la frontera
de R, de modo que
Φ0 (s) Φ0 (s)
Z Z
1 1
Im ds = Im ds,
4π R Φ(s) π γ Φ(s)
Nos falta estimar la última integral. Para ello observamos que Re ζ(s) no se
anula sobre γ1 , pues
∞ ∞ ∞
X e−it log n X cos(t log n) X 1
Re ζ(2 + it) = 1 + Re 2
= 1 + 2
≥ 1 − 2
n=2
n n=2
n n=2
n
π2
= 2 − ζ(2) = 2 − > 1/4.
6
Por lo tanto, si Re ζ(s) se anula en algún punto de γ, lo hace concretamente en
puntos de γ2 , es decir, en puntos de la forma σ + iT , con 1/2 ≤ σ ≤ 2. Además,
σ es entonces un cero de la función
1
g(s) = (ζ(s + iT ) + ζ(s − iT )),
2
que es holomorfa en C salvo en 1 ± iT .
Así, los ceros de Re ζ(s) sobre γ2 se corresponden con 2 + Ti
los ceros de una función holomorfa en un intervalo cerrado,
luego son un número finito. Digamos que g tiene m ceros en
[1/2, 2]. Vamos a aplicar el teorema [VC 2.19] a la función g
y a los discos D(2, 3/2) ⊂ D(2, 7/4). El menor de ellos con-
tiene al menos los m ceros de g que estamos considerando. 2
Como T > 3, tenemos que g es holomorfa en el disco mayor.
Por 4.7 con δ = 1/4 tenemos que, si |s − 2| = 7/4, entonces5
1
|g(s)| < c(|τ + T |3/4 + |τ − T |3/4 ) < c(2 + T )3/4 .
2
5 Por comodidad vamos a adoptar el convenio de que la letra c representará una constante
arbitraria que puede ir modificándose en cada paso. Así, en las dos desigualdades siguientes
la misma letra c representa a dos constantes distintas.
126 Capítulo 4. La función dseta de Riemann I
Por otra parte, hemos visto que g(2) = Re ζ(2 + iT ) > 1/4, luego
m
7
≤ 4c(2 + T )3/4 < T,
6
para T suficientemente grande. Por lo tanto m ≤ c log T , para T suficientemente
grande.
Por otra parte, γ queda dividida en m + 1 arcos consecutivos en los que
Re ζ(s) tiene signo constante. Si γj : [a, b] −→ C es uno de estos arcos, entonces
ζ 0 (s)
Z Z
1
Im ds = Im dz.
γj ζ(s) γj ◦ζ z
N (T + h) − N (T − h) = O(log T ),
Demostración: Llamemos
t t t
f (t) = log − ,
2π 2π 2π
de modo que
1 t
f 0 (t) =
log .
2π 2π
Por el teorema del valor medio, existe un T < α < T + h tal que
f (T + h) − f (T ) = f 0 (α),
luego
1 α 1 T +h
N (T +h)−N (T ) = log +O(log T ) ≤ log +O(log T ) = O(log T ).
2π 2π 2π 2π
Igualmente se prueba que N (T ) − N (T − h) = O(log T ), y la conclusión es
inmediata.
Más adelante necesitaremos esta estimación:
Teorema 4.21 X 1
= O(log2 T ).
Im ρ
0<Im ρ≤T
Ya en 1859 Riemann afirmó que era “muy probable” que todos los ceros
no triviales de la función dseta estuvieran sobre la recta crítica, y por ello tal
afirmación se conoce como la hipótesis de Riemann. Hasta el momento nadie la
ha demostrado ni refutado. Al margen del interés que tiene en sí misma, veremos
que tiene consecuencias destacables sobre la distribución de los números primos.
Π0 (z) 1 1
= log z + 1 + − 1 + µ0 (z) = log z + + µ0 (z),
Π(z) 2z 2z
donde
+∞
u − E[u] − 1/2
Z
µ(z) = − du,
0 z+u
luego
+∞
u − E[u] − 1/2
Z
µ0 (z) = du,
0 (z + u)2
y así, para x = Re z > 0,
Z +∞
1 +∞
Z
0 1/2 1 1
|µ (z)| ≤ 2 2
du ≤ 2
du = .
0 (x + u) + y 2 0 (x + u) 2x
130 Capítulo 4. La función dseta de Riemann I
En particular, si σ ≤ −1,
0
Π (−s) 1 1 π
Π(−s) ≤ | log(−s)| + 2 + 2 = log |s| + 2 + 1 < c log |s|,
En total, 0
ζ (s) π
ζ(s) < log(2π) + 2 c + c log |s| + c < c log |s|.
donde ρ recorre los ceros no triviales de la función dseta que cumplen la condi-
ción indicada.
Demostración: Observemos que la función del enunciado no se hace nunca
infinita. Podría pensarse que lo es cuando s = ρ es un cero no trivial de la función
dseta, pero hay que entender que en el sumatorio ρ aparece repetido tantas veces
como indica su multiplicidad, es decir, que la suma contiene a o(ζ, ρ)/(s − ρ),
pero esto es justamente la parte singular de la serie de Laurent de ζ 0 (s)/ζ(s)
alrededor de ρ, luego ésta se cancela con parte del sumatorio y la expresión toma
un valor finito.
4.4. La fórmula de Riemann-von Mangoldt 131
Tenemos que
ζ(s) = π s/2 Π−1 (s/2)(s − 1)−1 ξ(s) =
1 Y s
s −1 −1
(2π) Π (s/2)(s − 1) 1− es/ρ ,
2 ρ
ρ
luego
X 1 1
+ < c log τ.
| Im ρ−τ |<1 s + 3 − ρ ρ
Ahora que hemos separado los ceros que podrían coincidir con s, pasamos a
estimar la suma con s en lugar de s + 3. Para ello usamos 4.20:
X
1 1 X 1 1
+ − + =
Im ρ≥τ +1 s − ρ ρ s + 3 − ρ ρ
Im ρ≥τ +1
X X
1 1 3
= − =
Im ρ≥τ +1 s − ρ s + 3 − ρ Im ρ≥τ +1 (s − ρ)(s + 3 − ρ)
X 3 X 1
≤ =3
|τ − Im ρ||τ − Im ρ| (τ − Im ρ)2
Im ρ≥τ +1 Im ρ≥τ +1
∞ ∞
X X 1 X X 1
=3 ≤ 3
n=1 τ +n≤Im ρ<τ +n+1
(τ − Im ρ)2 n=1
n2
τ +n≤Im ρ<τ +n+1
∞
X log(τ + n) X log 2τ X log 2n
<c 2
≤c 2
+c
n=1
n n n>τ
n2
n≤τ
∞ ∞ ∞
X log 2 Xlog 2n X 1
≤c 2
+ c 2
+ c log τ 2
< c log τ.
n=1
n n=1
n n=1
n
Similarmente,
X X
1 1 1 1
+ − + < c log τ,
Im ρ≤τ −1 s − ρ ρ s + 3 − ρ ρ
Im ρ≤τ −1
luego
X
1 1
+ < c log τ.
| Im ρ−τ |≥1 s − ρ ρ
ζ 0 (s) X 1 1 1 Π0 (s/2)
− = log 2π − −
ζ(s) s−ρ s − 1 2 Π(s/2)
| Im ρ−τ |<1
X 1 1 X 1
+ + + .
s−ρ ρ ρ
| Im ρ−τ |≥1 | Im ρ−τ |<1
vez una función holomorfa en un entorno de ρ0 (la parte del sumatorio corres-
pondiente a otros posibles ceros). Así, ambos miembros son funciones continuas
de δ en un entorno de 0, luego la igualdad vale también para δ = 0, es decir,
para s = ρ0 .
En cuanto al miembro derecho, los sumandos del último término cumplen
1
≤ 1 ≤1
ρ τ − 1
π 1 σ
− arctan ≤ c.
4 2 2
Ahora ya podemos continuar nuestro estudio de ζ 0 (s)/ζ(s):
Teorema 4.24 Para cada T ≥ 2 existe T < τ < T + 1 con la propiedad de que
todo cero no trivial de la función dseta cumple
| Im ρ − τ |−1 ≤ c log τ
Por 4.20, el número de sumandos es O(log τ ), y cada uno de ellos tiene módulo
luego
ζ 0 (s)
= O(log2 τ ).
ζ(s)
Obviamente, tomando conjugados se concluye que ζ 0 (s)/ζ(s) = O(log2 |τ |)
en la banda −1 ≤ σ ≤ 2 y |τ | ≥ 2 siempre que la parte imaginaria sea opuesta
a una de las que cumplen el teorema anterior. Respecto al semiplano σ ≥ 2 no
hay nada que probar, pues en él la derivada logarítmica está acotada.
Terminamos con otra aplicación del teorema 4.23:
Teorema 4.25 La hipótesis de Lindelöf equivale a que, para todo σ > 1/2, se
cumple que
N (σ, T + 1) − N (σ, T ) = o(log T ).
donde ρ recorre los ceros de ζ(s) que cumplen la condición indicada. Por hi-
pótesis, dado > 0, tenemos que ζ(σ + iτ ) = O(τ /2 ) para σ = 12 + 4δ y para
σ = 27 − 4δ , y en la prueba del teorema de Lindelöf se ve que, para cualquier σ
intermedio, se cumple que
|ζ(σ + iτ )| ≤ M τ /2 ,
3 δ
2 − 4
X r
N log 3 δ
≤ log ≤
2 − 2
|ρ − 2 − iT |
|ρ−2−iT |<r1
2 + Ti
X r
log = o(log T ).
|ρ − 2 − iT |
|ρ−2−iT |<r
1
P
= o(log(T + k)) = o(log T ).
k=−2
136 Capítulo 4. La función dseta de Riemann I
ζ 0 (s) X 1 ζ 0 (s) X 1 X 1 X 1
− = − + − .
ζ(s) s−ρ ζ(s) s−ρ s−ρ s−ρ
ρ∈D | Im ρ−τ |<1 | Im ρ−τ |<1 | Im ρ−τ |≥1
ρ∈D
/ ρ∈D
contradicción. Por lo tanto, todos los sumandos de los dos últimos sumatorios
cumplen |s − ρ| ≥ δ.
Los ceros con | Im ρ − τ | < 1 cumplen también | Im ρ − T | < 1 + 3/2, luego
su número es O(log T ), por el teorema 4.20. Por lo tanto, el penúltimo suma-
torio es O(δ −1 log T ). Lo mismo vale para el último sumatorio, pues los ceros
de D cumplen | Im ρ − T | < 3/2, luego también son O(log T ) (de hecho, he-
mos visto que por hipótesis son o(log T ), pero en este punto de la prueba no
podemos aprovechar realmente la hipótesis del teorema, ya que en cualquier
caso sólo podemos afirmar que los dos últimos sumatorios son, conjuntamente,
O(δ −1 log T )).
Por último, teniendo en cuenta el teorema 4.23, para σ ≤ 2 tenemos que
3
En resumen, tenemos que, sobre la circunferencia |s − 2 − iT | = 2 − 2δ, la
función holomorfa F está acotada por M (δ, T ) = (c/δ) log T .
Consideremos ahora un punto de la circunferencia |s − 2 − iT | = 1/2. Aquí
la situación es más simple, pues la circunferencia está contenida en el semiplano
σ > 1.25, luego ζ 0 (s)/ζ(s) = O(1) y |s − ρ| ≥ 1/2, luego todos los sumandos
que aparecen en el sumatorio que define a F están acotados por 2. Ahora sí que
podemos aprovechar que el número de sumandos es o(log T ), para concluir que
el sumatorio es o(log T ).
Así pues, sobre la circunferencia |s − 2 − iT | = 1/2 tenemos que F está
acotada por M0 (T ) = o(log T ).
Con esto podemos aplicar el teorema de las tres circunferencias [VC 2.22]
al anillo D(0; 21 , 32 − 2δ). La conclusión es que si |s − 2 − iT | = r, donde
1/2 ≤ r ≤ 3/2 − 3δ, entonces
X 1
− (log(2 + iT − ρ) − log( + 3δ + iT − ρ)).
2
ρ∈D
Aquí observamos que log ζ(2 + iT ) = O(1), pues log ζ admite un desarrollo en
serie de Dirichlet absolutamente convergente en σ > 1. Tomando partes reales
queda que
1 X 1
log |ζ( + 3δ + iT )| = o(log T ) + (log |2 + iT − ρ| − log | + 3δ + iT − ρ|).
2 2
ρ∈D
luego
|ζ( 21 + 3δ + iT )|
lím = 0.
T →+∞ T
Esto significa que µ( 12 + 3δ) = 0 y, por continuidad, µ(1/2) = 0.
Capítulo V
139
140 Capítulo 5. La función dseta y los números primos
sería una aproximación de π(x) mucho mejor que Il(x). La figura de la izquierda
muestra las gráficas de π(x), Il(x) y R(x). La que queda claramente por arriba
es Il(x). La figura de la izquierda muestra los valores de R(x) − π(x). Vemos
que para x ≤ 10 000 el error no excede las 6 unidades.
30
25
20
15
10
20 40 60 80 100
Esta función coincide con ψ(x) en los puntos en los que ésta es continua, es
decir, salvo si x es potencia de primo, en cuyo caso ψ0 (x) = ψ(x) − 21 Λ(x) es el
punto medio del salto finito que presenta ψ en x.
Con esta salvedad, la fórmula explícita de von Mangoldt es la siguiente:
Teorema 5.1 Si1 x ≥ 2, se cumple que
X xρ 1
1
ψ0 (x) = x − − log 2π − log 1 − 2 ,
ρ
ρ 2 x
que es la serie análoga al primer término, pero para los ceros triviales de la
función dseta, en lugar de los no triviales.
El hecho de que la serie que aparece en la fórmula de von Mangoldt no
sea absolutamente convergente introduce dificultades en la demostración, las
cuales desaparecen en la prueba de otra fórmula explícita para una función más
sencilla, a saber, la integral
Z x X
ψ1 (x) = ψ(t) dt = Λ(n)(x − n).
1 n≤x
Las dos series que aparecen en esta fórmula son absolutamente convergentes.
La primera porque ρ+1
x xRe ρ+1 x2
ρ(ρ + 1) = |ρ||ρ + 1| ≤ |ρ|2 ,
ys c
≤ y ,
|f (s)| =
k+1
s(s + 1) · · · (s + k) |t|
Por otra parte, sobre C ∗ tenemos que |s+j| ≥ |s|−j ≥ R −k ≥ R −R/2 = R/2,
así como que |y s | = y σ ≤ y c , y aquí es crucial que y ≥ 1. Por lo tanto
k+1 c
2k+1 y c
Z
f (s) ds ≤ 2πR 2 y
k+1
< 2π .
C
R Tk
Los residuos los podemos calcular con el teorema [VC 3.16]: consideramos
k
Q
la función P (s) = (s + r), con lo que
r=0
y −j
Res(f, −j) = 0 ,
P (−j)
Q
ahora bien, P (s) = (s + j) (s + r), luego
r6=j
j−1 k
P 0 (−j) =
Q Q Q
(r − j) = (r − j) (r − j) =
r6=j r=0 r=j+1
j−1 k
(−1)j (k − j) = (−1)j j!(r − j)!
Q Q
(j − r)
r=0 r=j+1
Por lo tanto,
k k k
X X y −j 1 X k 1
Res(f, −j) = j
= (−1/y)j = (1 − 1/y)k .
j=0 j=0
(−1) j!(k − j)! k! j=0 j k!
1 Z 2±T i xs+1 ζ 0 (s) 3c log2 T x3 c log2 (2X)x3
ds ≤ ≤ .
2πi −1±T i s(s + 1) ζ(s) 2πT 2 X2
Todas las cotas tienden a 0 cuando X tiende a +∞, luego la integral sobre R
tiende a la integral sobre la recta σ = 2, que por el teorema anterior es ψ1 (x).
El teorema de los residuos nos da entonces que
P ∞
P
ψ1 (x) = − Res(f, −1) − Res(f, 0) − Res(f, 1) − Res(f, ρ) − Res(f, −2n),
ρ n=1
ζ 0 (−1)
Res(f, −1) = − .
ζ(−1)
ζ 0 (0)
Res(f, 0) = x = x log 2π
ζ(0)
x log2 (xT ) x
+O( + mín{1, } log x),
T T hxi
donde hxi > 0 representa la distancia de x a la potencia de primo más próxima.
De aquí se obtiene 5.1 sin más que hacer que T tienda a +∞. Necesitamos
una variante de 5.3:
y mín{1, πT | 1log y| } si y 6= 1,
c
|I(y, T ) − δ(y)| < c
πT si y = 1.
Observemos que −T i
1 Z −X−iT y s T y −X T
ds ≤ < ,
2πi −X+iT s πX πX
1 Z c−iT y s 1 Z c y u−iT
ds = du
2πi −X−iT s 2πi −X u − iT
Z c u
1 y 1
≤ du = (y c − y −X ),
2π −X T 2πT log y
e igualmente
1 Z −X+iT y s 1
ds ≤ (y c − y −X ).
2πi c+iT s 2πT log y
luego
1 Z c+iT y s 1 p yc
ds − 1 ≤ 2π c2 + T 2 √ = yc .
2πi c−iT s 2π
c2 + T 2
5.1. Fórmulas explícitas 147
1 T 1 T /c 1 1 +∞ 1
Z Z Z
c 1
= du = dv = − dv,
π 0 c2 + u 2 π 0 1 + v2 2 π T /c 1 + v 2
luego
1 Z c+iT 1 1 1 +∞ 1
Z
1 +∞ 1
Z
c
ds − ≤ dv < dv = .
2πi c−iT s 2 π T /c 1 + v 2 π T /c v 2 πT
Ahora definimos
c+iT
ζ 0 (s) xs
Z
1
J(x, T ) = − ds.
2πi c−iT ζ(s) s
Teorema 5.7 Si x > 0, c > 1 y T > 0, se cumple que
∞
X 1 c
|J(x, T ) − ψ0 (x)| < Λ(n)(x/n)c mín{1, } + Λ(x),
T | log(x/n)| T
n=1(6=x)
xc x log2 x x
|J(x, T ) − ψ0 (x)| = O( + + mín{1, } log x).
T (c − 1) T T hxi
c 3 log x Ax log2 x
Λ(x) ≤ ≤ ,
T T T
pues si Λ(x) 6= 0, es que x = pk y Λ(x) = log p ≤ log x.
Ahora vamos a descomponer en varios sumandos la serie que aparece en el
enunciado del teorema anterior. Para todo T suficientemente grande, tenemos
que
X 1 X Λ(n) 1
Λ(n)(x/n)c mín{1, } ≤ xc mín{1, }
T | log(x/n)| nc T log(4/3)
n≤(3/4)x n≤(3/4)x
∞
xc X Λ(n) xc ζ 0 (c) Axc
≤ = − ≤ ,
T log(4/3) n=1 nc T log(4/3) ζ(c) T (c − 1)
donde en la última desigualdad usamos que
ζ 0 (s) 1 1 + (s − 1)f (s) A
− = + f (s) = ≤ ,
ζ(s) s−1 s−1 s−1
donde f (s) es una función holomorfa en un entorno del intervalo [1, 3], por lo
que 1 + (s − 1)f (s) está acotada superiormente en dicho intervalo.
Similarmente, si n ≥ (5/4)x tenemos que | log(x/n)| ≥ | log(4/5)|, y llegamos
a la misma conclusión.
A continuación consideramos los sumandos para los que (3/4)x < n < x. Si
x1 es la mayor potencia de primo menor que x, podemos suponer que se cumple
(3/4)x < x1 < x, pues en caso contrario todos los sumandos indicados son nulos.
Para n = x1 tenemos que
x − x1 x − x1
log(x/n) = − log 1 − ≥ ,
x x
luego el sumando correspondiente en la serie es
1 x
Λ(n)(x/n)c mín{1, } ≤ Λ(x1 )(4/3)3 mín{1, }
T | log(x/n)| T (x − x1 )
x
= A mín{1, } log x
T hxi
(pues x1 = pk y Λ(x1 ) = log p ≤ log x1 < log x).
2 Como en este contexto la c tiene un significado concreto, en lo sucesivo usaremos la letra A
Finalmente consideramos los términos para x < n < (5/4)x, para lo cual
llamamos x2 a la menor potencia de primo mayor que x, y podemos suponer
que x < x2 < (5/4)x, o de lo contrario todos los términos son nulos. Para
n = x2 tenemos que
x2 − x x2 − x
| log(x/n)| = − log 1 − ≥ ,
x2 x2
luego el sumando correspondiente en la serie es
1 x
Λ(n)(x/n)c mín{1, } ≤ Λ(x2 ) mín{1, }
T | log(x/n)| T (x2 − x)
x
= A mín{1, } log x.
T hxi
Los demás términos corresponden a índices n = x2 + m, con 0 < m < x/4, luego
m m
| log(x/n)| = − log(x/n) > − log(x2 /n) = − log − ≥ .
x2 + m x2 + m
Por lo tanto,
x2 + m 2x + x/4 Ax
| log(x/n)|−1 ≤ ≤ ≤ ,
m m m
donde hemos usado el postulado de Bertrand 3.9, que nos asegura que existe un
primo (luego una potencia de primo) ≤ 2x. A su vez,
X 1
Λ(n)(x/n)c mín{1, }≤
T | log(x/n)|
x2 <n<(5/4)x
X Ax Ax log x X 1 Ax log2 x
Λ(x1 + m) ≤ ≤ .
Tm T m T
0<m<x/4 0<m<x/4
Reuniendo todas las cotas que hemos obtenido, tenemos la conclusión buscada.
mx1/m
π(x1/m ) ≤ c .
log x
5.2. Estimación del error en el teorema de los números primos 153
Por lo tanto,
x x
π(t) − Il(t)
Z Z
Il(t)
= π(x) log x − dt − dt.
2 t 2 t
Integrando por partes:
Z x
Il(t)
dt = Il(x) log x − x + 2,
2 t
X 1
X 1 1 E[x]
= n − + ,
ns ns (n + 1)s E[x]s
n≤x n≤x−1
y los restos
σ
ψ(x) ψ(x)/x x E[x] 1
= , E[x]s = E[x]σ−1
E[x]s xσ−1 E[x]
5.2. Estimación del error en el teorema de los números primos 155
x log2 (xT ) x
+O( + mín{1, } log x),
T T hxi
Al suponer que x es entero resulta que hxi ≥ 1, luego
Así pues,
x log2 (xT )
ψ0 (x) = x + O(xα log2 T + ).
T
Tomamos concretamente T = x1−α , y así
Teorema 5.13 Existe una constante c > 0 tal que no existen ceros de la función
dseta en la región
c
1− ≤ σ.
log(|τ | + 2)
Por lo tanto,
∞
ζ 0 (s) X Λ(n)
− Re = cos(τ log n).
ζ(s) n=1
nσ
Y en consecuencia
ζ 0 (σ) ζ 0 (σ + τ i) ζ 0 (σ + 2τ i)
−3 − 4 Re − Re ≥ 0. (5.3)
ζ(σ) ζ(σ + τ i) ζ(σ + 2τ i)
Como −ζ 0 (s)/ζ(s) tiene un polo simple en 1 con residuo 1, se cumple que
ζ 0 (σ) 1
− ≤ + O(1).
ζ(σ) σ−1
Combinando dos igualdades obtenidas en la prueba del teorema 4.16 tenemos
que
ζ 0 (s) 1 Π0 (s/2) X 1
1 1
− = + − + +c
ζ(s) 2 Π(s/2) s−1 ρ
s−ρ ρ
En la prueba del teorema 4.22 hemos visto que
Π0 (s/2) s 1
= log + + O(1/σ),
Π(s/2) 2 s
luego, si 1 ≤ σ ≤ 2, τ ≥ 2,
0
Π (s/2)
Π(s/2) ≤ log |s/2| + π/2 + 1 + O(1) = O(log τ ).
Así pues,
ζ 0 (s)
X 1 1
− Re =− Re + + O(log τ ).
ζ(s) ρ
s−ρ ρ
Los términos de la serie son no negativos, pues si ρ = β + γi,
1 σ−β 1 β
Re = ≥ 0, Re = 2 ≥ 0.
s−ρ |s − ρ|2 ρ |ρ|
Por lo tanto,
ζ 0 (σ + 2τ i)
− Re ≤ O(log τ ).
ζ(σ + 2τ i)
Fijemos ahora un cero ρ y tomemos s con 1 < σ ≤ 2 tal que τ = Im ρ ≥ 2. Si
eliminamos todos los términos de la suma menos el correspondiente a ρ resulta
que
ζ 0 (σ + τ i) 1
− Re ≤− + O(log τ ).
ζ(σ + τ i) σ−β
158 Capítulo 5. La función dseta y los números primos
Ahora sustituimos en (5.3) las cotas que hemos obtenido, y así obtenemos
que
3 4
0≤ − + O(log τ ).
σ−1 σ−β
Así pues, existe una constante c (independiente de s o de ρ) tal que
4 3
≤ + c log τ,
σ−β σ−1
donde β es la parte real de un cero no trivial, τ es su parte imaginaria y σ es
cualquier número real 1 < σ ≤ 2. En realidad, si σ ≥ 2 también tenemos que
4 4 3 1 3 3
≤ = + ≤ +1≤ + c log τ,
σ−β σ−1 σ−1 σ−1 σ−1 σ−1
si tomamos c ≥ 1/ log 2. Así pues,
4
σ−β ≥ ,
3/(σ − 1) + c log τ
luego
4 1 − c(σ − 1) log τ
1−β ≥ − (σ − 1) = .
3/(σ − 1) + c log τ 3/(σ − 1) + c log τ
Ahora elegimos concretamente σ = 1 + 1/(2c log τ ), para asegurar que el
numerador sea positivo. Así obtenemos que
1/2
1−β ≥ ,
5c log τ
y esto tiene que cumplirlo la parte real de cualquier cero no trivial con parte
imaginaria positiva. Cambiando la constante podemos pedir que
c
1− ≤β
log(τ + 2)
y, como los conjugados de los ceros también son ceros, lo mismo vale si τ < 0,
con |τ | en lugar de τ .
muestra la gráfica.
5.2. Estimación del error en el teorema de los números primos 159
donde
x log2 (xT ) x x log2 (xT )
R(x, T ) = O( + mín{1, } log x) = O( ),
T T hxi T
donde la última igualdad se debe a que si x es un número natural, entonces
hxi ≥ 1. Por el teorema anterior existe una constante tal que todo cero ρ tal
que | Im ρ| ≤ T cumple
c c
Re ρ < 1 − ≤1− .
log | Im ρ| log T
Por lo tanto,
c log x
|xρ | = xRe ρ = eRe ρ log x < xe− log T ,
luego, por 4.21,
ρ
X
x c log x
< xe− log T
X 1 c log x X 1
= 2xe− log T
| Im ρ|<T ρ
|ρ| Im ρ
| Im ρ|<T 0<Im ρ<T
c log x
= O(x log2 T e− log T ).
Llevando todo a la fórmula explícita,
x log2 (xT )
c log x
ψ(x) = ψ0 (x) + O(log x) = x + O(log x + x log2 T e− log T + ).
T
√ para todo número natural x ≥ 2 y todo T ≥ 2. Ahora tomamos x ≥ 3
Esto vale
y T = e log x > 2. Así obtenemos que
√
√
−c log x x(log x + log x)2 √
−c0 log x
ψ(x) = x+O(log x+x log x e + √ ) = x+O(xe ),
e log x
donde c0 es cualquier constante que cumpla c0 < 1, c0 < c. En efecto,
√
x log x e−c log x log x t
−c
√
0 log x
= (c−c0 )√log x = (c−c0 )t → 0,
xe ) e e
√
x(log x + log x)2 (t2 + t)2
(1−c
√
0 ) log x
= (1−c0 )t → 0.
xe e
160 Capítulo 5. La función dseta y los números primos
Notemos que hasta ahora sabíamos que ψ(x)/x convergía a 1, pero el teorema
anterior nos dice, más concretamente, que
ψ(x) √
= 1 + O(e−c log x ),
x
luego nos da información sobre la velocidad a que ψ(x)/x se acerca a 1. Obser-
vemos también que √
xe−c log x
lím = 0,
x→+∞ x/ log x
√
por lo que la estimación O(xe−c log x ) que proporciona el teorema siguiente es
mejor que la estimación O(x/ log x) que es la que, en principio, tenemos para
π(x) − Il(x).
Teorema 5.15 Existe una constante c > 0 tal que
√
π(x) = Il(x) + O(xe−c log x
).
Demostración: Por el teorema 5.10 tenemos que
√ √
ϑ(x) − x = O(xe−c log x
+ x1/2 log2 x) = O(xe−c log x
),
de donde, a su vez, de nuevo por 5.10,
√ √
x
c0 xe−c log x te−c log t √
Z
∗
π (x) − Il(x) = + c0 + c0 2 dt = O(xe−c log x ),
log x 2 t log t
pues el integrando está acotado, luego la integral es O(x). Por último
√
π(x) − Il(x) = π(x) − π ∗ (x) + π ∗ (x) − Il(x) = O(x1/2 ) + O(xe−c log x
)
√
−c log x
= O(xe ).
Más adelante probaremos (véanse las observaciones tras el teorema 6.16)
que, en realidad, el teorema anterior es cierto para toda constante c. La única
salvedad es que, cuanto mayor sea la c que elijamos, mayor tendrá que ser tam-
bién la constante implícita en O o, alternativamente, obtendremos una acotación
válida únicamente para x suficientemente grande. Por ejemplo, puede probarse
que √
|π(x) − Il(x)| < 0.1e−0.3 log x para x ≥ 59.
Si aumentamos el 0.1 hasta 0.7 la desigualdad vale para todo x, pero, a cambio
de incluir unos pocos números, la hacemos menos precisa para números mayores.
30
120
25
100
20
80
15
60
10 40
5 20
La figura muestra las funciones π(x) e Il(x) (en el centro) y las cotas superior
e inferior para π(x) que proporciona la desigualdad anterior.
5.2. Estimación del error en el teorema de los números primos 161
Veamos algunas aplicaciones del teorema de los números primos con la esti-
mación del error.
Teorema 5.16 Para todo n ≥ 0 y todo x > 2, se cumple que
n
x X k! x
π(x) = k
+ O( n+1 ).
log x log x log x
k=0
Una variante del teorema de los números primos Las gráficas muestran
que Il(x) aproxima mucho mejor a π(x) que x/ log x, pero la situación se puede
mejorar un poco.
x
log x−1 Il(x)
π(x)
150
x
log x
100
50
La figura muestra las gráficas de π(x), Il(x), x/ log x y también x/(log x−1).
Vemos que Il(x) sigue estando más cerca de π(x) (sus gráficas se superponen),
pero la función modificada x/(log x − 1) se acerca más a π(x) que la original
x/ log x.
162 Capítulo 5. La función dseta y los números primos
Esto ya lo había observado Gauss, que conjeturó que, de entre las funciones
x
,
log x − a
la que más se ajusta a π(x) es la correspondiente a a = 1. Ahora podemos
probar que tenía razón. Para ello definimos A(x) como la función dada por
x
π(x) =
log x − A(x)
o, equivalentemente
x
A(x) = log x − .
π(x)
Vamos a probar que
1
A(x) = 1 + O( ),
log x
lo que quiere decir que, cuanto mayor es x, el número que habría que restarle a
log x para obtener el valor exacto de π(x) se parece más a 1. Tenemos que
x log x log2 x
(A(x) − 1) log x = log2 x − log x − = log2 x − log x − ,
π(x) π(x) log x/x
y el teorema anterior nos proporciona la estimación
π(x) log x 1 2 1
=1+ + + R(x), R(x) = O( ).
x log x log2 x log2 x
Así,
log2 x
(A(x) − 1) log x = log2 x − log x − 1
1+ log x + log22 x + R(x)
1 + log2 x R(x) − 2
log x − log x R(x)
= 1 2 → 1.
1+ log x + log2 x
+ R(x)
En la introducción dejamos planteado que la suma de los primos ≤ x es
asintóticamente equivalente a Il(x2 ). Ahora podemos demostrar este hecho, e
incluso algo más general:
Teorema 5.17 Si k ≥ 1, se cumple que
X xk+1
pk ∼ ∼ Il(xk+1 ).
(k + 1) log x
p≤x
Por lo tanto,
x x
tk
X Z Z
pk = (π(x) − Il(x))xk − k (π(t) − Il(t))tk−1 dt + dt.
2 2 log t
p≤x
x
√ √ xk+1
Z
0 k+1 −c log x 0
≤cx e + kc tk dt = c0 xk+1 e−c log x
+ kc0
2 k+1
√
= O(xk+1 e−c log x
).
Por consiguiente,
P k
p − Il(xk+1 ) xk+1 √
k+1 −c log x
(k+1) log x 0 x e log x
p≤x
≤ c ≤ c0 c√log x ,
Il(xk+1 ) Il(xk+1 ) xk+1
(k+1) log x
e
Teorema 5.18 Sea η : [0, +∞[ −→ ]0, 1/2] una función decreciente de clase C 1
y supongamos que lím η 0 (t) = 0 y que 1/η(t) = O(log t). Si ζ(s) no tiene ceros
t→+∞
que cumplan 1 − η(|τ |) < σ y 0 < α < 1, entonces
ζ 0 (s)
= O(log2 |τ |)
ζ(s)
ζ 0 (s) 1
− = + f (s),
ζ(s) s−1
s0 = 1 + αH + iT s0
1
1 − (1 + α)H + i(T ± R) ∈ Ω± ,
2
es decir, que cumple
1
1 − η(|T | + R) < 1 − (1 + α)H,
2
como η es decreciente, esto valdrá para para todos los puntos con parte real
mayor y parte imaginaria menor (en valor absoluto), luego para todos los puntos
del círculo de radio R, y también es claro que si |T | es suficientemente grande la
circunferencia de radio R < 1 no cortará al eje real. Ahora bien, por el teorema
del valor medio, existe un |T | < t < |T | + R tal que
luego
1 1
1 − (1 + α)H − 1 + η(R + |T |) = − (1 + α)H + η(|T |) + Rη 0 (t)
2 2
1 1 1−α 1
= − (1 + α)H + H + (1 + 3α)Hη 0 (t) = ( + (1 + 3α)η 0 (t))H > 0
2 2 2 2
si |T | es suficientemente grande, ya que η 0 es negativa, pero tiende a 0.
Así pues, log ζ(s) es holomorfa en el disco |s − s0 | < R. En dicho disco
σ > 1/2 y |T | − 1 < |τ | < |T | + 1. Por 4.7, para δ = 1/2, tenemos que
por (4.4), teniendo en cuenta que a la serie le falta el primer término. Así
estamos en condiciones de aplicar el teorema [VC 2.17] para k = 1, según el
cual, si |s − s0 | ≤ r,
0
ζ (s) 2R 4(1 + 3α) 1
ζ(s) ≤ (R − r)2 (log |T | − Re log ζ(s0 )) < (1 − α)2 H (log |T | + αH ).
166 Capítulo 5. La función dseta y los números primos
luego las dos últimas integrales tienden a 0 cuando T tiende a ∞. Esto implica
que
x2 xs+1 ζ 0 (s)
Z
1
ψ1 (x) = − ds.
2 2πi φ s(s + 1) ζ(s)
3 Usamos A para representar constantes distintas.
5.3. Regiones sin ceros y el error en el teorema de los números primos 167
log2 (|t| + 2)
Z +∞ −αη(|t|)
ψ1 (x) 1 x
x2 − 2 ≤ A dt
−∞ |φ(t)(φ(t) + 1)|
+∞
|t|α log2 (|t| + 2)
Z
=A e−αη(|t|) log x−α log |t| dt
−∞ |φ(t)(φ(t) + 1)|
1
|t|α log2 (|t| + 2) log2 (t + 2)
Z Z +∞
≤ Ae−ω(x) dt + 2 2−α
dt = O(e−ω(x) ),
−1 (1/2)(3/2) 1 t
donde la última integral converge porque δ = 2 − α > 1, por lo que, para
0
cualquier 1 < δ 0 < δ, el integrando es claramente O(1/tδ ), y la convergencia de
Z +∞
dt
1 tδ 0
y
ω(x2 ) ≤ η(t1 ) log x2 + log t1 = ω(x1 ) + η(t1 )(log x2 − log x1 )
< ω(x1 ) + log x2 − log x1 .
Como ω(1) = 0 y log 1 − ω(1) = 0, concluimos además que
Veamos ahora cómo obtener una estimación del error de ψ(x) ∼ x a partir
de la que acabamos de obtener para ψ1 (x) y a su vez de aquí pasaremos a otra
para π(x):
168 Capítulo 5. La función dseta y los números primos
Teorema 5.20 Sea η : [0, +∞[ −→ ]0, 1/2] una función decreciente de clase C 1
y supongamos que lím η 0 (t) = 0 y que 1/η(t) = O(log t). Si ζ(s) no tiene ceros
t→+∞
que cumplan 1 − η(|τ |) < σ y 0 < α < 1, entonces
1 1
ψ(x) = x + O(xe− 2 αω(x) ), π(x) = Il(x) + O(xe− 2 αω(x) ),
1 x 1 x+h
Z Z
ψ(t) dt ≤ ψ(x) ≤ ψ(t) dt.
h x−h h x
los términos de los extremos son
ψ1 (x ± h) − ψ1 (x) (x ± h)2 /2 − x2 /2
=
±h ±h
1 1
+ O((x − h)2 e−αω(x−h) ) + O((x + h)2 e−αω(x+h) )
h h
h e−αω(x/2)
=x±+ O(x2 ),
2 h
donde hemos aplicado el teorema anterior y, para la última igualdad, hemos
usado que
2
(x ± h)2 2h h
2
≤1+ + ≤ 3,
x x x
así como que x/2 < x − h < x < x + h y, como ω es creciente
e−αω(x±h) ≤ e−αω(x/2) .
Por lo tanto,
h e−αω(x)
|ψ(x) − x| ≤ + O(x2 ).
2 h
1
Tomando h(x) = 12 xe− 2 αω(x) obtenemos la primera afirmación del enunciado.
Para probar la segunda parte, el teorema 5.10 nos da que
Z x
∗ ψ(x) − x 2 ψ(t) − t
π (x) − Il(x) = + + dt
log x log 2 2 t log2 t
Por la parte ya probada del teorema,
ψ(x) − x − 21 αω(x)
log x ≤ cxe
5.3. Regiones sin ceros y el error en el teorema de los números primos 169
1 1
y, como xe− 2 αω(x) > xe− 2 α log x = x1−α/2 > x1/2 > 1, eligiendo la constante c
suficientemente grande, también
2 1
< cxe− 2 αω(x) .
log 2
Además,
Z x Z x α (log t−ω(t)) Z x
ψ(t) − t e2 α
(log x−ω(x)) dt
2 dt ≤ c1 α/2
dt ≤ e 2
α/2
2 t log t 2 t 2 t
− 12 αω(x) α/2−1 1 1
xe x −
(1 − α/2)x1−α/2 (1 − α/2)21−α/2
− 12 αω(x) 1 1
= xe −
(1 − α/2)x2−α (1 − α/2)21−α/2 x1−α/2
1
= O(xe− 2 αω(x) ).
Concluimos que
1
π ∗ (x) − Il(x) = O(xe− 2 αω(x) ).
Por otra parte, el teorema 5.10 nos da también que π ∗ (x) − π(x) = O(x1/2 ),
de donde
1 1
π(x) − Il(x) = O(xe− 2 αω(x) ) + O(x1/2 ) = O(xe− 2 αω(x) ),
1
pues hemos visto que x1/2 /xe− 2 αω(x) < 1.
c log x
+ log t
log(t + 2)
√
para t ≥ 1 cumple ω(x) ≥ k log x, para cierto k > 0. En efecto, esto equivale
a que √
c log x log t
+√ ≥ k > 0.
log(t + 2) log x
Si 1 ≤ t ≤ 2 tenemos que
√ √
c log x log t c log 2
+√ ≥ .
log(t + 2) log x log 4
4 Alternativamente, si nos fijamos en los pocos pasos de las pruebas precedentes en los que
interviene η 0 , vemos que en realidad el teorema anterior sigue siendo válido si admitimos que
η no sea derivable en un número finito de puntos, con tal de que η 0 esté acotada donde está
definida. Así, en lugar de reducir la constante c, podemos usar η ∗ (t) = mín{1/2, η(t)}.
170 Capítulo 5. La función dseta y los números primos
Para t ≥ 2, existe una constante a > 0 tal que log(t + 2) ≤ a log t, con lo que
√ √ r
c log x log t c log x log t c
+√ ≥ +√ ≥2 ,
log(t + 2) log x a log t log x a
√
√ ver que la función f (u) = cu + 1/u alcanza su mínimo en u = 1/ c,
pues es fácil
y éste es 2 c.
Por consiguiente, el teorema anterior, con α = 1/2, por ejemplo, nos da que
que es la misma estimación que nos proporciona el teorema 5.15 (aunque obte-
nida ahora de forma más laboriosa).
Si suponemos la hipótesis de Riemann, entonces la función η(t) = 1/2 tam-
bién cumple las condiciones del teorema anterior, y ω(x) = 21 log x, con lo que
obtenemos
para todo 0 < α < 1, que es una estimación peor que (5.2).
Así pues, el teorema anterior no nos proporciona ninguna información que
no supiéramos ya. Naturalmente, su interés estriba en que puede aplicarse a
teoremas que proporcionen regiones sin ceros más amplias que la del teorema
5.13, o más complicadas que la que proporciona la hipótesis de Riemann, para
las que no hay argumentos directos que nos proporcionen una estimación del
error del teorema de los números primos más fácilmente que aplicando el teorema
anterior.
En la sección siguiente veremos un ejemplo del interés que tiene obtener
estimaciones mejores del error.
se cumple que
kxθ
π(x + kxθ ) − π(x) ∼ ,
log x
luego pn+1 − pn = O(pθn ).
pn+1 − pn = O(p1/2+
n ), para todo > 0.
1/2+
En realidad, esto equivale a que pn+1 − pn = o(pn ), pues
pn+1 − pn pn+1 − pn 1
1/2+
= 1/2+/2 1/2+/2
,
pn pn pn
0.08 0.8
0.06 0.6
0.04 0.4
0.02 0.2
0 200 400 600 800 1000 0 200 400 600 800 1000
172 Capítulo 5. La función dseta y los números primos
luego
N (T )
T ≤ cT
b/2−1
logB T.
En efecto, por una parte, ψ(x) = ψ0 (s)+O(log x), lo que nos permite cambiar
ψ0 por ψ. Por otra parte,
T log 2π
2 log 2π ≤ = O(1),
x log x log2 x
T 1 1 1 2 log x
2 2 log 1 − ≤ 2 (log(x2 − 1) + log x2 ) ≤ = O(1),
x log x x 2 2 log x log2 x
2
x log2 (xT ) log x2
T
≤ = 4,
x log2 x T log x
T x 1
2 mín{1, } log x ≤ = O(1).
x log x T hxi log x
5.4. Diferencias entre primos 173
luego
ψ(x + h) − ψ(x) X x log2 x
=1+O xRe ρ−1 + O( ).
h Th
| Im ρ|<T
Por otra parte, si {ρn }Nn=1 es una enumeración de los ceros de la función
dseta que cumplen | Im ρn | < T , repetidos según su orden, tenemos que
X X Z Re ρ
Re ρ−1 −1
(x −x )= xu−1 log x du
| Im ρ|<T | Im ρ|<T 0
N Z
X Re ρn N Z
X 1
= xu−1 log x du = χ[0,Re ρn ] (u)xu−1 log x du
n=1 0 n=1 0
Z N
1X Z 1 X
u−1
= χ[0,Re ρn ] (u)x log x du = xu−1 log x du
0 n=1 0 | Im ρ|<T
Re ρ≥u
Z 1
=2 N (u, T )xu−1 log x du.
0
Así pues,
X Z 1
Re ρ−1 −1
x = 2x N (T ) + 2 N (u, T )xu−1 log x du.
| Im ρ|<T 0
X Z 1−η(T )
xRe ρ−1 = 2x−1 N (T ) + 2 N (u, T )xu−1 log x du.
| Im ρ|<T 0
174 Capítulo 5. La función dseta y los números primos
1
1−θ <α< ,
b + B/A
de modo que
1 B
>b+ ≥b
α A
y así δ > 0 y
y así ψ(x + kxθ ) − ψ(x) ∼ kxθ (omitimos la k del término de error porque sólo
modifica la constante implícita, que ya dependía de k). Por otra parte,
X X
ψ(x + kxθ ) − ψ(x) = log p + log p.
x<p≤x+kxθ x<pm ≤x+kxθ
m≥2
5.5. Orden de crecimiento y localización de ceros 175
Para cada primo p tal que p2 ≤ x+kxθ , el número de veces que log p aparece
en el segundo sumatorio es menor que el máximo m que cumple pm ≤ x + kxθ ,
luego
log(x + kxθ )
X X
ψ(x + kxθ ) − ψ(x) −
log p≤ log p
θ
2 θ
log p
x<p≤x+x p ≤x+kx
X p
≤ log((k + 1)x) ≤ (k + 1)x log((k + 1)x) = O(x1/2 log x),
p2 ≤(k+1)x
luego
X
ψ(x + kxθ ) − ψ(x) = log p + O(x1/2 log x)
x<p≤x+kxθ
X X
= log x + log(p/x) + O(x1/2 log x)
x<p≤x+kxθ x<p≤x+kxθ
x + kxθ
= (π(x + kxθ ) − π(x)) log x + kxθ log( ) + O(x1/2 log x)
x
luego
k
|ψ(x+kxθ )−ψ(x)−(π(x+kxθ )−π(x)) log x| ≤ kxθ log(1+ )+O(x1/2 log x).
x1−θ
A su vez,
φ(t)
= o(eφ(t) ) y ζ(s) = O(eφ(τ ) )
θ(t)
para todo s que cumpla 1 − θ(τ ) ≤ σ ≤ 2, τ ≥ t0 , entonces existen un c > 0 y
un t1 ≥ t0 tal que ζ(s) no tiene ceros en la región
θ(2τ + 1)
1−c ≤ σ, τ ≥ t1 .
φ(2τ + 1)
Por ejemplo, si tomamos θ(t) = 1/2 y φ(t) = log(t), la hipótesis sobre la
función dseta es que ζ(s) = o(τ ), para 1/2 ≤ σ ≤ 2, lo cual se cumple por el
teorema 4.7. La conclusión es que ζ(s) no tiene ceros en la región
c
1− ≤ σ, τ ≥ t1
log(2τ + 1)
y, como
log(2τ + 1)
= O(1),
log(τ + 2)
esto equivale a que ζ(s) no tenga ceros en una región de la forma
c
1− ≤ σ,
log(|τ | + 2)
con lo que recuperamos el teorema 5.13. Naturalmente, la finalidad de este
teorema es mejorar 5.13 a partir de estimaciones mejores del crecimiento de la
función dseta.
Necesitamos un resultado previo sobre funciones holomorfas:
Teorema 5.23 Sea f una función holomorfa en un disco D(s0 , r) y tal que
f (s0 ) 6= 0 y |f (s)| < |f (s0 )|eM en el disco, para cierto M > 0. Supongamos
además que f no se anula en el semicírculo formado por los puntos de D(s0 , r)
con Re s > Re s0 . Entonces
− Re(f 0 (s0 )/f (s0 )) < 4M/r.
Si además f tiene un cero ρ0 en el segmento abierto de extremos s0 − r/2 y s0 ,
entonces
4M 1
− Re(f 0 (s0 )/f (s0 )) < − .
r s0 − ρ0
Demostración: Cambiando f por f (s)/f (s0 ) podemos suponer sin pérdida
de generalidad que f (s0 ) = 1. Igualmente, cambiando f por f (s − s0 ) podemos
suponer que s0 = 0. Consideremos la función
f (s)
g(s) = Q ,
1 − ρs
ρ
donde ρ recorre los ceros de f en el disco |s| ≤ r/2 (que son un número finito y
por hipótesis son todos no nulos), repetidos tantas veces como indica su orden.
5.5. Orden de crecimiento y localización de ceros 177
de donde
X 1 4M
− Re f 0 (0) − Re < ,
ρ
ρ r
o equivalentemente,
4M X 1
− Re f 0 (0) < + Re .
r ρ
ρ
En particular
ζ 0 (σ0 ) 5/4
− < .
ζ(σ0 ) σ0 − 1
Sea s0 = σ0 + bi, s00 = σ0 + 2bi, r = θ(2b + 1) ≤ 1. Entonces los discos
|s − s0 | ≤ r y |s − s00 | ≤ r están contenidos en la región
1 − θ(τ ) ≤ σ ≤ 2, τ ≥ τ0 .
178 Capítulo 5. La función dseta y los números primos
σ ≥ σ0 − r ≥ 1 − θ(2b + 1) ≥ 1 − θ(τ ).
∞
P µ(n)
Como 1/ζ(s) = ns y |µ(n)| ≤ 1, para 1 < σ ≤ 2
n=1
tenemos que bi
1 2
ζ(s) ≤ ζ(σ) ≤ σ − 1 .
En particular
1 2 φ(2b+1)
ζ(s0 ) ≤ σ0 − 1 ≤ 2e .
0 1 σ0
Igualmente
1 φ(2b+1)
ζ(s0 ) ≤ 2e .
0
Veamos ahora que existe una constante A1 > 0 tal que, si |s − s0 | ≤ r,
ζ(s) A1 φ(2b+1)
ζ(s0 ) < e .
Todo esto es válido para cualquier σ0 que cumpla (5.5). Observemos ahora
que, si elegimos t1 suficientemente grande, podemos tomar, concretamente,
1 θ(2b + 1)
σ0 = 1 + .
40A2 φ(2b + 1)
En efecto, se trata de comprobar que con esta elección se cumple
1 θ(2b + 1)
e−φ(2b+1) ≤ <δ
40A2 φ(2b + 1)
si b > t1 . Ahora bien, por hipótesis el término central tiende a 0, luego podemos
asegurar que sea < δ, y la otra desigualdad equivale a
φ(2b+1)
θ(2b+1) 1
< ,
eφ(2b+1) 40A2
y por hipótesis el miembro izquierdo tiende a 0. Con esta elección queda que
θ(2b + 1)
1−a≥ ,
1240A2 φ(2b + 1)
luego a + bi está en la región del enunciado. Esto es así bajo la hipótesis de que
a > σ0 − r/2, que hemos usado en (5.7). Pero si a ≤ σ0 − r/2, la misma elección
de σ0 nos da que
r 1 θ(2b + 1) θ(2b + 1)
a ≤ σ0 − =1+ − ,
2 40A2 φ(2b + 1) 2
180 Capítulo 5. La función dseta y los números primos
luego
20A2 − 1 θ(2b + 1)
= ,
40A2 φ(2b + 1)
y de nuevo tenemos la conclusión.
Capítulo VI
181
182 Capítulo 6. La función dseta de Riemann II
por lo que χ(s)χ(1 − s) = 1. Por otra parte, el teorema 4.8 nos da que
por lo que
X 1 X 1 X 1 X 1
ζ(s) − − χ(s) = ζ(s̄) − − χ(s̄)
ns n1−s ns̄ n1−s̄
n≤x n≤y n≤x n≤y
Teorema 6.3 Sea m un número natural, sea 0 < < 2π y sea C la unión de
las curvas C1 , C2 , C3 dadas por
Entonces, la función
z s−1 e−mz
Z
I(s, m) = dz
C ez − 1
es entera, no depende de la elección de y
X 1 e−πis Π(−s)
ζ(s) = + I(s, m).
ns 2πi
n≤m
Como en [VC 5.30], hay que entender que z s−1 se calcula con el logaritmo
con parte imaginaria nula sobre C1 , con parte imaginaria en [0, 2π] sobre C2 y
con parte imaginaria 2π sobre C3 .
Demostración: Como en [VC 5.30], la integral sobre C2 define una función
entera. Sobre C1 y C2 basta probar que, cuando |s − 1| ≤ M , el integrando está
acotado por una función integrable independiente de s. Ahora bien, sobre C1
tenemos que z = x > y |z s−1 | = xσ−1 ≤ xM , mientras que sobre C3 es
luego
+∞ +∞
xs−1 e−mx z s−1 e−mz xs−1 e−mx
Z Z Z
I(s, m) = − dx+ dz+e2πis dx
ex − 1 |z|= ez − 1 ex − 1
+∞
xs−1 e−mx z s−1 e−mz
Z Z
= (e2πis − 1) dx + dz.
ex − 1 |z|= ez − 1
Veamos seguidamente que, si σ > 1, la última integral tiende a 0 cuando
→ 0. En efecto, la función g(z) = e−mz /(ez − 1) es holomorfa en D(0, 2π)
6.1. La ecuación funcional aproximada 185
salvo en z = 0, donde tiene un polo simple, por lo que zg(z) está acotada. Por
lo tanto,
Z Z
z s−1 e−mz 2π (s−1)(log +iθ)
iθ iθ
dz = e g(e )ie dθ
z
|z|= e − 1
0
Z 2π
≤ cσ−1 e−τ θ dθ ≤ cσ−1 2πe2π|τ | ,
0
que ciertamente tiende a 0 con . Así pues, por el teorema anterior,
Z +∞ s−1 −mx
x e X 1
I(s, m) = (e2πis − 1) dx = (e 2πis
− 1)Π(s − 1) ζ(s) − .
0 ex − 1 ns
n≤m
Equivalentemente,
X 1 I(s, m) e−πis I(s, m)
ζ(s) = = =
ns (e2πis − 1)Π(s − 1) (eπis − e−πis )Π(s − 1)
n≤m
e−πis I(s, m)
= ,
2i(sen πs)Π(s − 1)
y así [VC 4.40] nos da la fórmula del enunciado, en principio para σ > 1, pero
por el principio de prolongación analítica la igualdad vale para todo s.
El paso siguiente es una variante de [VC 5.35]:
Teorema 6.4 Fijados 0 ≤ c ≤ 1/2, δ > 0, 0 < r < 2πcδ, r ≤ π/2, tomamos
s ∈ C tal que τ > δ y dos números reales x, y > δ. Llamamos m = E[x],
q = E[y], η = 2πy. Consideramos la curva φ = φ1 ∪ φ2 ∪ φ3 ∪ φ4 descrita en
la figura, donde φ2 es el segmento que une cη + (1 + c)ηi con −cη + (1 − c)ηi
salvo si éste corta a un disco D(2kπi, r), necesariamente con k = q o k = q + 1,
en cuyo caso rodeamos el disco con un arco de circunferencia de radio r, el
adecuado para que 2qπi quede por debajo del arco, pero 2(q + 1)πi quede por
encima. Entonces
e−πis Π(−s)
X 1 X 1 Z s−1 −mz
z e
ζ(s) = + χ(s) + dz,
ns n1−s 2πi φ ez − 1
n≤x n≤y
donde z s−1 = e(s−1) log z se calcula con el logaritmo de z con parte imaginaria
en ]0, 2π[.
x x
si x > log√2 (pues entonces e − 1 > e /2). Para x ≥ |α| se cumple que
|φj (x)| ≤ 2x. Por consiguiente, si llamamos f (z) al integrando, tenemos que
ηi
−cη + (1 − c)ηi
4πi
2πi
−2πi
−4πi
z s−1 e−mz
Z
+ dz.
φ ez − 1
El teorema [VC 3.16] nos permite calcular los residuos:
Res(f, 2nπi) + Res(f, −2nπi) = (2nπi)s−1 + (−2nπi)s−1
= (2nπ)s−1 (eπi(s−1)/2 + e3πi(s−1)/2 ) = eπis (2nπ)s−1 (−ie−πis/2 + ieπis/2 )
πs
= −2eπis (2nπ)s−1 sen .
2
188 Capítulo 6. La función dseta de Riemann II
Por consiguiente,
q
X
e−πis Π(−s) (Res(f, 2nπi) + Res(f, −2nπi))
n=1
πs X 1 X 1
= −2(2π)s−1 Π(−s) sen 1−s
= −χ(s) .
2 n n1−s
n≤q n≤q
|s|2 σ2 σ02
= + 1 ≤ + 1 = A2 ,
τ2 τ2 δ2
luego
1/2−σ
|s| |s|
≤ ≤ A.
τ τ
Por consiguiente,
e−πis Π(−s) = O(τ 1/2−σ eπτ /2 ).
Ahora nos ocupamos de la integral, pero antes conviene hacer una observa-
ción general, y es que, si A > 0, e−Aτ = O(y −1 ).
En efecto,
ye−Aτ = ye−A2πxy ≤ ye−A2πδy
y la última expresión tiende a 0 con y, luego está acotada.
Para cada número complejo z que no esté en el semieje real positivo, sobren-
tenderemos que sus coordenadas polares son z = ρeiθ , con 0 < θ < 2π.
Supongamos en primer lugar que z = φ1 (u) = u + (1 + c)ηi. Entonces
ρ ≥ (1 + c)η ≥ η y, si u ≥ log 2,
|ez − 1| ≥ eu − 1 ≥ eu /2.
6.1. La ecuación funcional aproximada 189
Además,
(1 + c)η
θ = arctan .
u
Como
d (1 + c)η u (1 + c)η 1 (1 + c)η 1
arctan + =− 2 2 2
+ >− 2 2
+ = 0,
du u η u + (1 + c) η η (1 + c) η η
y u ≥ cη, tenemos que
(1 + c)η u 1+c π c π
θ = arctan + ≥ arctan + c = + c − arctan = + A,
u η c 2 1+c 2
pues
Z c/(1+c) Z c/(1+c)
c dp dp
arctan = < = c.
1+c 0 1 + p2 0 (1 − p)2
Notemos que A < c ≤ 1/2. Teniendo en cuenta además que m + 1 ≥ x = τ /η,
vemos que
s−1 −mz
z e σ−1 −τ (π/2+A−u/η)−τ u/η
= Kη σ−1 e−τ (π/2+A) .
ez − 1 ≤ Kη e
1
≤ Kπη σ e−τ (π/2+A) + Kη σ−1 e−xπη ≤
x
1 K σ −τ (π/2+A)
Kπη σ e−τ (π/2+A) + Kη σ−1 e−τ π ≤ η e ≤
δ 2πδ 2
K σ −τ (π/2+A)
y e ,
δ2
pues, teniendo en cuenta que A ≤ 1/2 < π/2,
η σ−1 e−πτ 1 1
σ −τ (π/2+A)
= e−τ (π/2−A) ≤ .
η e η 2πδ
190 Capítulo 6. La función dseta de Riemann II
Por lo tanto,
|z s−1 e−mz | = ρσ−1 e−τ θ emcη ≤ η σ−1 e−τ (π/2+A+c) eτ c = η σ−1 e−τ (π/2+A) .
luego
Z (1−c)η
z s−1 e−mz
Z
σ−1 σ−1 −τ (π/2+A)
dz ≤ η ≤ Aη e du
φ3 ez − 1
−(2q+1)π
luego θ > 5π/4. Por otra parte, el mínimo de ρ se alcanza cuando u = 0, con lo
que ρ ≥ (2q + 1)π ≥ (q + 1)π > yπ = η/2. Por otra parte,
luego
s−1 −mz
z e σ−1 −τ θ −mu
≤ (η/2)σ−1 e−5τ π/4 e−mu ≤ 2η σ−1 e−5τ π/4 e−mu .
ez − 1 ≤ ρ e e
Si m ≥ 1, Z +∞
z s−1 e−mz
Z
σ−1 −5τ π/4
e−mu du
dz ≤ Aη e
φ4 ez − 1 −cη
σ−1 mcη−5τ π/4
≤ Aη e ≤ Aη σ−1 eτ c−5πτ /4 ≤ 2πAy σ−1 eτ c−5πτ /4 ,
donde hemos usado que mη = E[x]2πy ≤ 2πxy = τ . De aquí, a su vez,
≤ Aη σ e−5τ π/4 + Aη σ−1 e−5τ π/4 e−πη ≤ Aη σ e−5τ π/4 + Aη σ−1 e−5τ π/4 e−π2πxy
= Aη σ e−5τ π/4 + Aη σ−1 e−5τ π/4 e−πτ = Aη σ e−5τ π/4 + Aη σ−1 e−9τ π/4
= Aη σ e−5τ π/4 ≤ A2πy σ e−5πτ /4 ,
y nuevamente,
e−πis Π(−s) z s−1 e−mz
Z
dz = O(τ 1/2−σ y σ e−3πτ /4 ) = O(τ 1/2−σ y σ−1 ).
2πi φ4 ez − 1
Sólo falta considerar el caso en que z está sobre φ∗2 . Llamamos φ21 al seg-
mento √ √
φ21 (λ) = ηi + λeπi/4 , con − 2cη ≤ λ ≤ 2cη,
de modo que φ2 puede coincidir con φ21 , o bien es φ21 menos un segmento
central que es sustituido por un arco de circunferencia φ22 .
Nos ocupamos en primer lugar de φ21 . Si z = φ21 (λ), entonces
pues entonces la parte imaginaria del logaritmo está en ]−π/2, π/2[, y al su-
marle π obtenemos una parte imaginaria en ]0, 2π[. Ahora bien,
√ √ √
λ −iπ/4 λ 2 2 2 1
Re(1 + e )=1+ ≥1− = > 0.
η η 2 2 2 2
El logaritmo con parte imaginaria en ]0, 2π[ es el que necesitamos para cal-
cular la exponencial compleja z s−1 del integrando. Concretamente:
π λ −πi/4
z s−1 = e(s−1)( 2 i+log η+log(1+ η e ))
,
luego
λ −πi/4
|z s−1 | = η σ−1 e−π/2 eRe((s−1) log(1+ η e ))
.
El desarrollo de Taylor del logaritmo nos da que
λ −πi/4 λ λ2 λ3
log(1 + e ) = e−iπ/4 − 2 e−iπ/2 + 3 h0 (λ/η),
η η 2η η
√ √
para cierta función continua h0 : [− 22 , 22 ] −→ C (en realidad está definida
en ]−1, 1[, pero la restringimos a un intervalo compacto menor para tenerla
acotada). Por lo tanto,
√ √
λ −iπ/4 λ 2 2λ λ2
Re((s − 1) log(1 + e )) = (σ − 1) + τ − 2τ
η η 2 2 η 2η
λ3 λ3
+(σ − 1) Re h0 (λ/η) − τ Im h0 (λ/η).
η3 η3
Tanto σ − 1 como λ/η y h0 están acotados, luego concluimos que
√
λ2
(− π
2+
2 λ
2 η − 2η 2 )τ +(λ/η)3 h(λ/η)τ
|z s−1 | ≤ Aη σ−1 e ,
√ √
2 2
donde h : [− 2 , 2 ] −→ R es la parte real de h0 .
Por [VC 5.33] sabemos que existe una constante A > 0 tal que |ez −1|−1 ≤ A
en la región Ω que resulta de quitarle al plano complejo el disco abierto de centro
cada punto 2nπi y radio r. Tomemos un punto z ∈ Ω y distingamos varios casos:
• Si 0 ≤ Re z ≤ π/2, entonces
−mz+xz
e (x−m) Re z
≤ Aeπ/2 .
ez − 1 ≤ Ae
6.1. La ecuación funcional aproximada 193
• Si Re z < 0, entonces
−mz+xz
e e(x−m) Re z
ez − 1 = |ez − 1| ≤ A.
λ3 λ2
h(λ/η) 3
≤ 2,
η 4η
luego
λ2 λ3 λ2
− 2
τ + τ h(λ/η) 3 ≤ − 2 τ.
2η η 4η
De este modo
Z cη√2
z s−1 e−mz
Z
λ2
σ−1 −τ π/2 − 4η 2τ
z
dz ≤ Aη e √ e dλ ≤
φ21 e −1 −cη 2
Z +∞ τ 2
− 4η 2λ
Aη σ−1 e−τ π/2 e dλ
−∞
+∞
2Aη σ
Z
2
= √ e−τ π/2 e−t dt = O(η σ τ −1/2 e−τ π/2 ).
τ −∞
Ahora nos ocupamos del posible arco φ22 de centro 2kπi y radio r, cuya
longitud está acotada por 2πr, luego sólo tenemos que preocuparnos de acotar
el integrando. Recordemos que k = q o bien k = q + 1, pero que, en cualquier
caso, k 6= 0. Consideremos un punto
ri iθ
z = φ22 (θ) = 2kπi + reiθ = 2kπi(1 − e ).
2kπ
Entonces
π ri iθ
log z = i + log 2kπ + log(1 − e ),
2 2kπ
194 Capítulo 6. La función dseta de Riemann II
π π (σ − 1)ri iθ τ r iθ
(s − 1) log z = (σ − 1) i − τ + (s − 1) log 2kπ − e + e
2 2 2kπ 2kπ
(σ − 1)r2 2θi ri iθ τ r2 i ri iθ
+ 2 2
e h( e ) + 2 2 e2θi h( e ).
4k π 2kπ k 4π 2kπ
Es claro que el primer término y el sexto están acotados. Como hemos visto
al principio de la sección, no perdemos generalidad si suponemos que x ≤ y, y
vamos a ver que, bajo esta hipótesis, también está acotado el séptimo término.
Sólo hay que probar que lo está τ /k 2 . Ahora bien:
• Si y ≥ 1, tenemos que
2
2πy 2
τ 2πxy q+
≤ ≤ = 2π ≤ 8π.
k2 q2 q2 q
• Si δ ≤ y < 1, entonces
τ 2πy 2
≤ ≤ 2π.
k2 1
Así pues:
π τ r iθ
(s − 1) log z = −τ + (s − 1) log 2kπ + e + O(1).
2 2kπ
Por otra parte,
iθ iθ
e−mz = e−m(2kπi+re )
= e−mre .
Por consiguiente,
iθ τ
−m)−τ π
z s−1 e−mz = ere ( 2kπ 2 +(s−1) log 2kπ+O(1) .
Puesto que (ez − 1)−1 está acotada sobre φ∗22 , concluimos que, sobre este arco,
z s−1 e−mz
= O((2kπ)σ−1 e−τ π/2 ).
ez − 1
Más aún, k ≤ y + 1 ≤ 2y si y ≥ 1 y k = 1 ≤ y/δ en otro caso, luego en cualquier
caso k = O(y) y llegamos a que
z s−1 e−mz
Z
z −1
dz = O(η σ−1 e−τ π/2 ).
∗
φ22 e
En particular
X 1
|ζ(1/2 + τ i)| ≤ c + O(|τ |−1/4 )
√ n1/2
n≤ |τ |/2π
Ahora bien,
Z √|τ |/2π
X 1 1 p
≤1+ √ du ≤ 4 |τ |/2π = O(|τ |1/4 ),
√ n1/2 1 u
n≤ |τ |/2π
Teorema 6.5 Sea {bn } una sucesión decreciente de números reales ≥ 0, sea
n
P
sn = am , donde los am son números complejos, entendiendo que sM = 0.
m=M +1
Entonces N
X
an bn ≤ bM máx |sn |.
M <n≤N
n=M +1
luego
X N N
X −1
an bn ≤ |sn |(bn − bn+1 ) + |sN |bN
n=M +1 n=M +1
−1
N
!
X
≤ (bn − bn+1 ) + bN máx |sn | = bM +1 máx |sn | ≤ bM máx |sn |.
M <n≤N M <n≤N M <n≤N
n=M +1
La función ∞
1 r
P
t rz
ft (z) = e r=k+1
donde cada crk es una forma de grado m (es decir, un polinomio cuyos monomios
tienen todos grado m) con coeficientes positivos. Ciertamente, para m = 1 es
cr1 = t/r. Si vale para m, entonces
∞
X ∞
1 r m+1 X
∞
X t r
∞
X r−m(k+1)
X t r
t z = crm (t)z r z = cr−u,m (t) z
r r u
r=k+1 r=m(k+1) r=k+1 r=(m+1)(k+1) u=k+1
y, como cada cr−u,m (t) es una forma de grado m con coeficientes positivos, es
claro que lo mismo vale para el coeficiente de z r , pero ahora con grado m + 1.
A su vez,
∞
t
P 1 r ∞ ∞ ∞ ∞
rz 1 X 1 r m X 1
X X
e r=k+1
= t z = crm (t)z r
m=0
m! r m=0
m!
r=k+1 r=m(k+1)
∞ E[r/(k+1)]
X X crm (t)
= zr ,
r=0 m=0
m!
E[r/(k+1)]
P crm (t)
luego cr (t) = m! es un polinomio con coeficientes positivos.
m=0
Se cumple que
Xb Xb−a b−a
−it log n −it log(a+n) −it log a X −it log(1+ na )
e = e = e e =
n=a+1 n=1 n=1
k ∞
(−1)r−1 nr (−1)r−1 nr
b−a −it
X P P
r ar −it r ar
e r=1 r=k+1 ,
n=1
donde usamos que n/a ≤ (b − a)/a ≤ 1/2, luego la serie de Taylor de log(1 + z)
converge en n/a.
En términos del desarrollo precedente tenemos que
k
b b−a P (−1)r−1 nr X∞ n l
X X −it r
e−it log n =
r a
e r=1 (−1)l cl (it)
a
n=a+1 n=1 l=0
∞ b−a k
(−1)r−1 nr
X (−1)l cl (it) X P
−it r a r
l
= l
n e r=1 .
l=0 a n=1
Ahora usamos el teorema de Abel [VC 5.2], según el cual
b−a k
P (−1)r−1 nr
X −it r ar
l
ne r=1 =
n=1
b−a−1 u k b−a k
P (−1)r−1 nr P (−1)r−1 nr
X X −it r ar
X −it r ar
(ul − (u + 1)l ) e r=1 + (b − a)l e r=1 .
u=1 n=1 n=1
Por lo tanto,
b−a k
X l −it P (−1)r narr b−a−1
r−1
X
n e r=1 ≤
((u + 1)l − ul )M + (b − a)l M
n=1 u=1
Pero, como cl (t) tiene coeficientes positivos, se cumple que |cl (it)| ≤ cl (|it|),
luego
b ∞ l
X
−it log n
X b−a
e ≤ 2M cl (t)
a
n=a+1 l=0
∞ r ∞
P 1 (b−a)
P (b−a)r
t r ar t ar ((b−a)/a)k+1
= 2M e r=k+1
≤ 2M e r=k+1
≤ 2M et (1−(b−a)/a) .
200 Capítulo 6. La función dseta de Riemann II
A su vez
b−a k+1 1 −1/(k+1) k+1
t a t 2t 1
≤ =
1 − b−a
a
1/2 2k
k
y e1/2 < 2.
Esto nos lleva a su vez a estudiar la acotación de expresiones de la forma
m
X
S= e2πiP (n) ,
n=1
Hacemos el cambio n0 = n − r:
m n−1 m−1 mín{r+m,m}
X X 2 X X 2
|S|2 = e2πi(2αnr−αr +βr)
= e2πi(2αnr−αr +βr)
El primer factor tiene módulo 1 y el segundo se puede acotar como hemos visto.
Ahora obtenemos una cota para el caso general:
m−1
X mín{rX
1 +m,m}
Llamemos
mín{rX
1 +m,m} mín{rX
1 +m,m}
k−1
S1 = S1 (r1 ) = e2πi(P (n)−P (n−r1 )) = e2πi(αkr1 n +···)
.
n=máx{r1 +1,1} n=máx{r1 +1,1}
En estos términos,
m−1 m−1
!1−2/K m−1
!2/K
X X X
|S|2 ≤ |S1 (r1 )| ≤ 1 |S1 |K/2 ,
r1 =−m+1 r1 =−m+1 r1 =−m+1
m−1
!2/K
X
∗
|S|2 ≤ (2m)1−2/K mK/2 + |S1 |K/2 ,
r1 =−m+1
202 Capítulo 6. La función dseta de Riemann II
1 1
< 4a1− K t (k+1)K log1/K t,
pues esto equivale a que 1/4K−1 < 1 < log t. Es claro entonces que se cumple
la desigualdad del enunciado con c = 4.
Así pues, podemos suponer que
1 − k+1
1 1 1 1 1
at − at− k+1 = at− k+1 > 2.
2 4 2
Además, por (6.2) y (6.4),
b−a a
0≤N ≤ ≤ < 4t1/(k+1) , (6.5)
m m
mN ≤ b − a < mN + m.
Supongamos ahora que N ≤ 1, con lo que b − a < 2m, luego, por (6.4),
b
X 1 1
n ≤ b − a < 2m ≤ at− k+1 < at− k+1 log(k−1)/K a log1/K t,
−it
n=a+1
donde usamos que a > 4, t > 3, y de aquí se sigue la desigualdad del enunciado.
Por lo tanto, podemos suponer que N ≥ 2, con lo que
b−a
m≤ . (6.6)
2
En tal caso descomponemos
b
X a+m
X a+2m
X b
X
n−it = n−it + n−it + · · · + n−it
n=a+1 n=a+1 n=a+m+1 n=a+mN +1
= Σ1 + Σ2 + · · · ΣN +1 , (6.7)
donde en la última igualdad ha de entenderse que cada Σj es, por definición, el
sumando correspondiente del término anterior. Llamamos vj = a + jm, para
1 ≤ j ≤ N + 1, con lo que, por (6.5), (6.2) y (6.6),
a < a + m = v1 ≤ vj ≤ vN +1 = a + (N + 1)m ≤ a + m + b − a
b−a 5a
= b + m ≤ 2a + m ≤ 2a + ≤ < 3a. (6.8)
2 2
Por otra parte, por (6.4),
m m 1 1
< ≤ t− k+1 .
vj a 2
Así el teorema 6.6 nos da que |Σj | ≤ 4M . Por otra parte, según 6.7,
m−1
X t(k − 1)! r1 · · · rk−1 −1
|Sj |K ≤ 22K mK−1 + 2K mK−k mín{m, | sen( )| }.
r1 ,...,rk−1 =1
2(a + jm)k
A su vez,
1/K
m−1
X t(k − 1)! r1 · · · rk−1 −1
|Σj | ≤ 24 m1−1/K + 23 m1−k/K mín{m, | sen( )| } ,
r1 ,...,rk−1 =1
2(a + jm)k
≤ 24 (N + 1)m1−1/K + (6.9)
1/K
N m−1
1 X X t(k − 1)! r1 · · · rk−1 −1
23 m1−k/K (N + 1)1− K mín{m, | sen( )| } ,
j=0 r1 ,...,rk−1 =1
2(a + jm)k
Observemos que
k
1 1 vj+1 − vjk (vj+1 /vj )k − 1
− = =
vjk k
vj+1 vjk vj+1
k k
vj+1
1 + log N 2 2 1
≤m+2 < m + (1 + log 4t1/2 ) < m + (3 + log t)
D D D 2
7
< m + log t,
D
donde en la última desigualdad usamos que t > 3. Igualmente se concluye que
X 7
mín{m, | sen θj |−1 } < m + log t.
D
θj ∈[(g−1/2)π,gπ]
26 · 3k+1 − k+1
k 22k+8 − k+1
k 210k − k+1
k
≤ t ak log t < t ak log t ≤ t ak log t,
R R R
Por (6.10), el número de intervalos I que pueden contener algún θj es a lo
sumo
2
1 + tk!Ra−k < 1 + tk!Ra−k .
π
Por lo tanto,
N
X 210k − k+1
k
mín{m, | sen θj |−1 } < (1 + tk!Ra−k ) t ak log t
j=0
R
210k − k+1
k 1
= t ak log t + 210k k! t k+1 log t. (6.11)
R
Ahora:
m−1 m−1
!k−1
X 1 X 1
= ≤ (1 + log(m − 1))k−1
r1 ,...,rk−1
r · · · rk−1
=1 1 r=1
r
m−1
X 1 k
10k − k+1 k 10k 1
≤ 2 t a log t + 2 k! t k+1 log t
r1 ,...,rk−1 =1
r1 · · · rk−1
k k−1 1
≤ (2 log a)k−1 210k t− k+1 ak log t + ak−1 t− k+1 210k k! t k+1 log t
k 2−k
< 211K logk−1 a t− k+1 ak log t + 211K t k+1 log t ak−1 ,
k u−2 k−1
22 < 22 = 2K .
Q
donde hemos usado que k ≤ K y que k! ≤
u=2
A su vez:
N
X m−1
X 1/K
mín{m, | sen(θj )|−1 } <
j=0 r1 ,...,rk−1 =1
k 2−k
211 (log(k−1)/K a t− (k+1)K ak/K + t (k+1)K a(k−1)/K ) log1/K t,
6.2. El método de Hardy-Littlewood 207
donde hemos usado de nuevo que (véase la nota al pie de la página 204) que
(A + B)1/K ≤ A1/K + B 1/K . Llamamos V a la última expresión que hemos
obtenido. Con ella podemos volver a (6.9):
Xb 1
−it
n ≤ 24 (N + 1)m1−1/K + 23 m1−k/K (N + 1)1− K V.
n=a+1
1
+217 a1−1/K t (k+1)K ) log1/K t
1 1 1 k−1
< 218 a1− K t (k+1)K + at− (k+1)K log K a log1/K t.
y combinando esto con la fórmula (6.1) que nos proporciona la ecuación funcional
aproximada, concluimos que
En definitiva:
Teorema 6.10 Dado δ > 0, en la región del plano complejo determinada por
las condiciones 21 + δ ≤ σ ≤ 1, τ ≥ 3 se cumple que
X 1
ζ(s) = + O(1).
2
ns
n≤τ
E[τ 2 ]1/2−δ √ +∞
E[τ 2 ]1/2−δ √
Z
dt 1
≤ + 2τ ≤ + 2τ
τ E[τ 2 ] t3/2+δ τ (1/2 + δ)E[τ 2 ]1/2+δ
√ √ √
E[τ 2 ]1/2 1 2 2 2
≤ δ
+ 2δ
≤ δ + = O(1),
τ 9 (1/2 + δ)τ 9 (1/2 + δ)9δ
√ p √
donde hemos usado que |s|/τ ≤ 11 + τ 2 /τ = 1/τ 2 + 1 ≤ 2, y también que
√ p √
E[τ 2 ]1/2 /τ ≤ τ 2 + 1/τ = 1 + 1/τ 2 ≤ 2.
El caso k = 2 del teorema 6.8 basta para probar esta versión, válida en una
estrecha banda en el borde de la banda crítica:
23
Teorema 6.11 En la región del plano complejo determinada por 24 ≤ σ ≤ 1,
τ ≥ 3, se cumple que
X 1
ζ(s) = s
+ O(1)
√ n
n≤ τ
En efecto:
donde r es el único número natural que cumple 2r E[τ 1/2 ] < τ 2 ≤ 2r+1 E[τ 1/2 ].
Si 2h E[τ 1/2 ] ≤ τ 2 , por la observación tras el teorema 6.5,
h+1 1/2
2 E[τ ]
X 1 1
s
≤ hσ 1/2 ]σ
M,
n
n=2h E[τ 1/2 ]+1 2 E[τ
donde
m
X 1 h
2 E[τ 2 ] + 1 ≤ m ≤ 2h+1 E[τ 2 ] .
M = máx
niτ
n=2h E[τ 1/2 ]+1
con lo que
210 Capítulo 6. La función dseta de Riemann II
h+1 1/2
2 E[τ ]
X 1
≤ c(2−h(σ−1/2) E[τ 1/2 ]−(σ−1/2) τ 1/6
n s
(2h E[τ 1/2 ])1−σ τ −1/6 log1/2 (2h E[τ 1/2 ])) log1/2 τ
n=2h E[τ 1/2 ]+1
+
≤ c(E[τ 1/2 ]−(σ−1/2) τ 1/6 + τ 2−2σ τ −1/6 log1/2 τ 2 ) log1/2 τ
≤ c(τ −σ/2+1/4 τ 1/6 + τ 2−2σ τ −1/6 log1/2 τ ) log1/2 τ
≤ c(τ −1/16 log1/2 τ + τ −1/12 log τ ) = O(1/ log τ ).
τ2 τ2 3
h + 1 ≤ log ≤ log = log τ,
E[τ 1/2 ] τ 1/2 2
Notemos que 23/24 ≤ 31/32 = 1 − 1/26−1 ≤ 1 − 1/R, por lo que bajo las
hipótesis actuales se cumplen las del teorema anterior.
A su vez, para probar esto demostraremos primero que si k ≥ 3, K = 2k−1 ,
1
1 − 4K ≤ σ ≤ 1, τ > 3, entonces
X 1 1
= O(τ − 84K log2 τ ), (6.15)
ns
τ 2/(k+2) <n≤τ 2/(k+1)
donde
Xm
n−iτ 2h E[τ 2/(k+2) ] + 1 ≤ m ≤ 2h+1 E[τ 2/(k+2) ] .
M = máx
n=2h E[τ 2/(k+2) ]+1
1 k−1
+(2h E[τ 2/(k+2) ])1−σ τ − (k+1)K log K (2h E[τ 2/(k+2) ])) log1/K τ
3 1 2−2σ 1 k−1
≤ c((2h E[τ 2/(k+2) ])− 4K τ (k+1)K + τ k+1 τ − (k+1)K log K τ 1/2 ) log1/K τ
3 1 1 1
≤ c(τ − 2K(k+2) τ (k+1)K log1/K τ + τ 2K(k+1) τ − (k+1)K logk τ ).
Ahora usamos que
3 1 1−k 1 k−1
− + = =−
2K(k + 2) (k + 1)K 2K(k + 1)(k + 2) 2K(k + 2) k + 1
1 1 1
=−
<− ,
4kK 2 8kK
1 1 1 1
− =− <− ,
2K(k + 1) (k + 1)K 2K(k + 1) 8kK
luego
h+1 2/(k+2)
2 E[τ ]
X 1 1 1
≤ c(τ − 8kK log1/K τ + τ − 8kK logk/K τ )
n s
n=2h E[τ 2/(k+2) ]+1
1
= O(τ − 8kK log τ ).
Como el número de sumandos es a lo sumo log τ , obtenemos (6.15). Pasamos a
probar (6.14).
212 Capítulo 6. La función dseta de Riemann II
luego
− log 2
X 1 1 (log τ )
= O(τ − 4 log log τ log2 τ )
s
τ 2/(k+2) <n≤τ 2/(k+1) n
1−log 2
1 (log τ )
= O(e− 4 log log τ log2 τ ) = O(1/ log τ ).
Para comprobar la última igualdad basta ver que
α
1 (log τ )
lím e− 4 log log τ log3 τ = 0,
τ →+∞
x3
lím 1 xα
= 0.
x→+∞
e 4 log x
x3 x3 1
1 xα
≤ 10
= 9
.
e 4 log x P 1 xkα
P 1 1
+ 1 x10α−3
4k k! logk x 4k k! x3−kα logk x 410 10! log10 x
k=0 k=0
1 1
σ ≥1− >1− ,
2r R
donde R = 2r−1 . Se cumplen, pues, todas las hipótesis del teorema anterior.
Por consiguiente:
E[τ 2/r ]
X 1
|ζ(s)| ≤ + O(1). (6.16)
n=1
nσ
Ahora distinguimos dos casos:
log log τ
Caso 1 Si se cumple 1 − σ ≤ log τ , entonces
E[τ 2/r ]
X 1 4 log log τ log τ
σ
≤ elog τ log log τ log τ 2 log τ 2/r < 4e4
n=1
n log log τ
1 log τ
< 4e4 τ 4(1−σ)/ log 1−σ ,
log log τ
donde hemos acotado
E[τ 2/r ] E[τ 2/r ] Z n Z τ 2/r
X 1 X dx dx
≤1+ ≤1+ = 1 + log τ 2/r < 2 log τ 2/r ,
n=1
n n=2 n−1 x 1 x
log log τ
Caso 2 Si se cumple 1 − σ > log τ , entonces
1
log < log log τ − log log log τ < log log τ,
1−σ
luego " #
1
log 1−σ 1 1 1
r=E ≥ E[log ] > log
log 2 1−σ 2 1−σ
y
4(1−σ)
E[τ 2/r ] τ 2/r 1
τ 2/r(1−σ) τ log 1−σ log τ
Z
X 1 dx
< = < .
n=1
nσ 0 x σ 1−σ log log τ
Así pues, en ambos casos tenemos que
E[τ 2/r ]
X 1 1 log τ
< Cτ 4(1−σ)/ log 1−σ ,
n=1
nσ log log τ
ζ(s) = O(logA τ ).
(log log τ )2
Demostración: Notemos que 1 − log τ > 63/64, luego si σ ≤ 1 el
teorema anterior nos da que
4(log log τ )2
1 log τ log
log τ log τ
ζ(s) = O τ 4(1−σ)/ log 1−σ = O e (log log τ )2
log log τ log log τ
log5 τ
log τ
= O e4 log log τ = O( ) = O(log5 τ ).
log log τ log log τ
Para σ ≥ 1 se cumple que ζ(s) = O(log τ ), por el teorema 4.7.
6.2. El método de Hardy-Littlewood 215
para 1/2 ≤ σ ≤ 1.
Esto tiene interés a causa del teorema de Hoheisel, 5.21, pues, teniendo en
cuenta también el teorema 6.15, así podemos concluir que las hipótesis de 5.21
se satisfacen con b = 2 + 4c y B = 5 (y la A dada por el teorema 6.15), luego
llegamos a que todo θ que cumpla
1 + 4c + 5/A
<θ<1
2 + 4c + 5/A
satisface la conclusión del teorema. El mero hecho de que exista una constante
0 < θ < 1 que cumpla el teorema ya es algo notable en sí mismo, pero en el
capítulo siguiente demostraremos (véanse las observaciones tras el teorema 7.14)
que, en realidad, el teorema 6.15 se cumple para cualquier valor de A, luego
haciendo tender A a +∞ obtenemos una cota para θ que depende exclusivamente
de c, a saber:
Por hipótesis, Z σ2
lím |f (t ± iT )|2 dt = 0,
T →+∞ σ0
por lo que
Z b Z +∞
2
lím |f (γ2 (t))| dt = |f (σ2 + it)|2 dt = J(σ2 ),
T →+∞ a −∞
e igualmente
Z b
lím |f (γ1 (t))|2 dt = J(σ1 ),
T →+∞ a
lo que nos da la conclusión para σ0 .
218 Capítulo 6. La función dseta de Riemann II
0
σ2 −σ0 0 −σ
σ0 1
J(σ1 ) σ2 −σ1 J(σ2 ) σ2 −σ1 .
Pero la conclusión es trivialmente cierta para σ1 y σ2 , luego tenemos el teorema
probado para todos los valores σ = σ1 + (σ1 − σ1 ) 2mn , con 0 ≤ m ≤ 2n . Como
J(σ) es una función continua, esto implica que se cumple para todo σ.
La parte más técnica de la prueba del teorema de Ingham es el resultado
siguiente:
P µ(n)
Teorema 6.20 Sea MX (s) = ns , donde µ es la función de Möbius, y sea
n<X
pues X 2δ = e2δ log X > 16 (2δ log X)3 , por la serie de Taylor de la exponencial.
Por otra parte observamos que, si λ > 1, se cumple que
1 1
1 < log λ + < log λ + √ . (6.18)
λ λ
(basta ver que log λ − 1/λ − 1 es creciente, porque tiene derivada positiva).
220 Capítulo 6. La función dseta de Riemann II
∞
!2
X d(n) X d(m)d(n) 1
1+δ
+ 1/2 log(n/m) n1+δ
n=1
n (mn)
1≤m<n
Z +∞
4 d(m)d(n)
X 1+δ
= ζ (1 + δ) + dx
(mn) 1/2 log(n/m)
n x2+δ
1≤m<n
Z +∞
1+δ X d(m)d(n)
= ζ 4 (1 + δ) + χ
1/2 log(n/m) [n,+∞[
(x) dx
1 x2+δ (mn)
1≤m<n
Z +∞
4 1+δ X d(m)d(n)
= ζ (1 + δ) + dx
1 x2+δ (mn) 1/2 log(n/m)
m<n≤x
Z +∞
1+δ
4
< ζ (1 + δ) + c log3 x dx,
1 x1+δ
de nuevo por el teorema 2.24. Ahora, por una parte,
1 δ+1 2
ζ(1 + δ) ≤ 1 + = ≤ ,
δ δ δ
y por otra parte, como x−δ/2 log3 x está acotada,
Z +∞ Z +∞
1+δ 3 1 2A A
1+δ
A log x dx ≤ c 1+δ/2
dx ≤ ≤ 4.
1 x 1 x δ δ
En definitiva,
X d(m)d(n) A
1+δ
< 4.
(mn) log(n/m) δ
X≤m<n
Ahora obtendremos una estimación similar para 1/2 + iτ . Para ello obser-
vamos que de (x − y)2 ≥ 0 se sigue que (x + y)2 ≤ 2(x2 + y 2 ), luego
Como antes,
Z T Z T
2
X µ(m)µ(n)
|MX (1/2 + iτ )| dτ = dτ
0 0 m,n<X
m1/2+iτ n1/2−iτ
X µ(m)µ(n) Z T m iτ
= dτ ≤
m,n<X
(mn)1/2 0 n
X 1 Z T
X 1
T +4 cos log(m/n)τ dτ ≤
n
n<X m<n<X
(mn)1/2 0
X 1 X 1
T +4 .
n<X
n
m<n<X
(mn)1/2 log(m/n)
1 λ1/2
<1+ ,
log λ λ−1
luego
1 (mn)1/2
<1+ ,
log(m/n) m−n
luego
Z T X 1
X 1 1
|MX (1/2 + iτ )|2 dτ ≤ T +4 1/2
+
0 n (mn) n−m
n<X m<n<X
X X 2
≤ T log X + 4 ≤ T log X + 8X log X.
n
m<X n<X
1
T ≤ T 1+2c ≤ T 1+2c log X ≤ AT 2c (T + X) log X,
log 2
X 1 Z X
1 √ √
|fX (1/2 + iτ )| ≤ A( √ + 1) ≤ A( √ dx + 2) = A(2 X + 2) ≤ A X,
n 1 x
n<X
con lo que
Z T
|fX (1/2 + iτ )|2 dτ ≤ AT X ≤ AT 2c X ≤ AT 2c (T + X) log X,
0
T
I1/2 (T ) = O(T 2c (T + X) log X), I1+δ (T ) = O(( + 1)δ −4 ).
X
Consideramos ahora
s−1
φ(s) = fX (s), t > 3/π.
s cos(s/2t)
1 −|τ |/t σ 1 1
| cos(s/2t)|2 e−|τ |/t = e cos + (1 + e−2|τ |/t ) → (6.19)
2 t 4 4
cuando τ tiende a ±∞. Así pues,
1
= O(e−|τ |/t ).
| cos(s/2t)|2
6.3. El teorema de Ingham 223
Por 4.7, para σ ≥ 1/2 y |τ | ≥ 3 tenemos que ζ(s) = O(|τ |1/2 ), luego
s − 1 s − 1 s − 1 1/2 1/2
s fX (s) = s ζ(s)MX (s) − s ≤ A|τ | X + 2 ≤ AX|τ | .
lím |φ(σ + iτ )| = 0
|τ |→+∞
En total,
J1/2 ≤ At2c (t + X) log X.
Así podemos aplicar el teorema 6.19, según el cual, para cada 1/2 ≤ σ ≤ 1+δ,
tenemos que
σ−1/2
t −4
1/2+δ 1+δ−σ
Jσ ≤ A 1 + δ At2c (t + X) log X 1/2+δ ≤
X
1−2σ 4c(1+δ−σ)
X 2δ+1 t 2δ+1 (X + t) máx{Aδ −4 , A log X}.
Por otra parte, en dicha banda, para τ ≥ 1, en virtud de (6.19) tenemos que
1
Tomamos t = T , δ = A/ log(T + X), donde A = 4 log 3, para que se cumpla
0 < δ ≤ 1/4. Así,
1−2σ
X 2δ+1 ≤ X −(1−2δ)(2σ−1) ≤ X −(2σ−1)+2δ ≤ e2A X −(2σ−1) .
|x| |x|
Z Z
2 2
dxdy = lím 2 2
dxdy = 4r.
Br (0) x + y n A(0,1/n,r) x + y
para σ ≥ 2 y, por (2.5) sabemos que d(n) = o(n1/4 ), luego existe un X0 > 1 tal
que si X ≥ X0 , entonces, d(n) ≤ n1/4 y, usando el teorema 4.1,
X n1/4 X 1 X 1
|fX (s)| ≤ 2
= 7/4
= ζ(7/4) − ≤
n≥X
n
n≥X
n n<X
n7/4
X 1 1 1 4 1
ζ(7/4) − + ≤ + + 7/4 = O(X −3/4 ),
n≤X
n7/4 X 7/4 X 7/4 3X 3/4 X
Dados X > X0 y T > 4, elegimos T1 y T2 tales que 3 < T1 < 4 y T < T2 < T +1
y de modo que h no se anule en los segmentos 1/2 ≤ σ ≤ 2, τ = T1 , T2 .
Para 1/2 ≤ u ≤ 2, llamamos Ru al rectángulo de vértices T2 i
u+iT1 , 2+iT1 , 2+iT2 , u+iT2 . A lo sumo, hX puede anularse
en su lado izquierdo. Si no es así, tenemos que Ru
h0X (s)
Z
T1 i
ds = 2πiNh (u, T1 , T2 ).
Ru hX (s)
!
2 2+iT1 2+iT2 u+iT2
h0X (s) h0X (s) h0X (s)
Z Z Z Z
= ds + ds + ds du
σ0 u+iT1 hX (s) 2+iT1 hX (s) 2+iT2 hX (s)
2 u+iT1
h0X (s)
Z Z
+ dsdu.
σ0 u+iT2 hX (s)
Notemos que la integral sobre Ru no está definida para un número finito de
valores de u, pero esto es irrelevante. Lo mismo sucede con la integral interior
del último término. Desarrollándola tenemos que
Z 2 Z u+iT1 0 Z 2 Z T2 0
hX (s) hX (u + it)
dsdu = −i dtdu.
σ0 u+iT2 hX (s) σ0 T1 hX (u + it)
T2 2
h0X (u + it)
Z Z
= −i dudt,
T1 σ0 hX (u + it)
donde la integral interior no está definida para un número finito de valores de t.
Tomando partes imaginarias llegamos a que
Z 2 Z T2 Z 2 0
hX (u + it)
2π Nh (u, T1 , T2 ) du = − Re dudt +
σ0 T1 σ0 hX (u + it)
lo que nos permite aplicar el teorema anterior para acotar la integral del pri-
mer logaritmo. Para el segundo logaritmo no podemos usar el teorema, pero,
teniendo en cuenta (6.20), vemos que
− log |hX (2 + it)| ≤ − log(1 − |fX (2 + it)|2 ) ≤ 2|fX (2 + it)|2 < 1/X.
228 Capítulo 6. La función dseta de Riemann II
Por consiguiente, teniendo en cuenta que T2 < T +1 y que T2 −T1 ≤ T +1−3 < T ,
T2 2
h0X (u + it)
Z Z
− Re dudt <
T1 σ0 hX (u + it)
(T + 1)4c(1−σ0 ) T
A (T + X + 1) log4 (T + X + 1) + . (6.21)
X 2σ0 −1 X
Ahora estimaremos las otras tres integrales. La más sencilla es
2 2+iT2
h0X (s)
Z Z
Im dsdu ≤ 2(Im log hX (2 + iT2 ) − Im log hX (2 + iT1 )),
σ0 2+iT1 hX (s)
Ahora bien,
1 1
|Hj (s)| = | (hX (s + iTj ) + hX (s̄ + iTj ))| ≤ (|hX (s + iTj )| + |hX (s̄ + iTj )|)
2 2
y si |s − 2| = 7/4, entonces Re(s ± iTj ) = Re s ≥ 1/4 y
−2 ≤ Im s, Im s̄ ≤ 2,
luego
1 ≤ Im(s + iTj ), Im(s̄ + iTj ) ≤ T + 3,
y por lo tanto
mj
7
≤ máx{|hX (s)| | σ ≥ 1/4, 1 ≤ τ ≤ T + 3}.
6
6.3. El teorema de Ingham 229
Ahora bien,
|hX (s)| ≤ 1 + |fX (s)|2 ≤ 2 + |ζ(s)MX (s)|2 .
Por 4.7 tenemos que |ζ(s)| = O(τ 3/4 ) = O(T 3/4 ) = O((T + X)3/4 ), y
X 1 Z X
dx
|MX (s)| ≤ 1/4
≤ 1 + 1/4
= O(X 3/4 ) = O((T + X)3/4 ),
n<X
n 1 x
Concluimos que
Z 2
Nh (u, T1 , T2 ) du ≤ AT 4c(1−σ0 ) (T X 1−2σ0 + X 2(1−σ0 ) ) log4 (T + X),
σ0
Por otra parte, todos los ceros de la función dseta son ceros de h, luego,
llamando N (σ, T1 , T2 ) = N (σ, T2 ) − N (σ, T1 ), para 0 < δ < 1, tenemos que
Z 2 Z σ0 +δ
Nh (u, T1 , T2 ) du ≥ N (u, T1 , T2 ) du ≥ δN (σ0 + δ, T1 , T2 ).
σ0 σ0
Como la función dseta no tiene ceros con parte imaginaria ≤ T1 < 4, resulta
que N (σ, T ) ≤ N (σ, T1 , T2 ), luego, si 1/2 + δ ≤ σ ≤ 1, tomando σ0 = σ − δ,
tenemos que
A 4c(1−σ+δ)
N (σ, T ) ≤ T (T X 1−2σ+2δ + X 2(1−σ+δ) ) log4 (T + X). (6.22)
δ
230 Capítulo 6. La función dseta de Riemann II
El método de Vinogradov
x1 + · · · + xl = y1 + · · · + yl
x21 + ··· + x2l = y12 + · · · + yl2
............ ............
xk1 + · · · + xkl = y1k + · · · + ylk
tales que 1 ≤ xj , yj ≤ q.
Un hecho fundamental es que se cumple algo más general que lo que exige
la definición:
l l
(a + bxj )h = (a + byj )h ,
P P
h = 1, . . . , k,
j=1 j=1
tales que 1 ≤ xj , yj ≤ q.
231
232 Capítulo 7. El método de Vinogradov
que equivale a
h l h l
X h X X h X
r h−r
a b xh−r
j = ar bh−r yjh−r .
r=0
r j=1 r=0
r j=1
por lo que x0j , yj0 son una solución del sistema de 7.2 cambiando a por −ac, luego
el número de soluciones es igualmente Jlk (q).
El segundo hecho fundamental se deriva de la observación elemental de que,
si c es entero, entonces
Z 1
2πicx 1 si c = 0,
e dx =
0 0 si c 6= 0.
f k (x; a0 , . . . , ak ) = ak xk + · · · + a1 x + a0 ,
y a su vez
X
S k (q, a; a0 , . . . , ak ) = e2πif (n) . (7.1)
a<n≤a+q
El teorema siguiente puede interpretarse como que Jlk (q) es el valor medio
de |S k (q, a)|2l cuando las variables a1 , . . . , al varían sobre el cubo unitario. Por
ello la expresión Jlk (q) se conoce como “valor medio de Vinogradov”:
7.1. La conjetura de Vinogradov 233
Por lo tanto,
Z k Z 1
X Y h h h h
|S k (q, a)|2l da1 · · · dak = e2πi(m1 +···+ml −n1 +···−nl )ah dah
[0,1]k a<m1 ,...,ml ≤a+q h=0 0
a<n1 ,...,nl ≤a+q
1 + 4 + 4 = 9 = 2 + 2 + 5, 12 + 42 + 42 = 33 = 22 + 22 + 52 .
es q 2l , ya que existe una única solución para cada valor que demos arbitraria-
mente a las variables xj , yj (dentro del rango permitido). Así pues,
Podemos precisar cuál de las dos cotas inferiores es mejor. Basta observar
que
(2l)−k q 2l−k(k+1)/2
= (2l)−k q l−k(k+1)/2 ,
ql
luego si l ≤ k(k + 1)/2 el denominador es mayor, mientras que si l > k(k + 1)/2
es mayor el numerador siempre que q ≥ lk . Por lo tanto, para q ≥ lk , podemos
desglosar la cota inferior así:3
l
q si l ≤ k(k + 1)/2,
Jlk (q) ≥
(2l)−k q 2l−k(k+1)/2 si l > k(k + 1)/2.
Vinogradov conjeturó en 1935 que esta estimación por defecto del creci-
miento del valor medio Jlk (q) está muy cerca de ser una estimación por exceso.
Concretamente, conjeturó que, para todo > 0, cuando q tiende a ∞ se cumple
que
O(q l+ ) si l ≤ k(k + 1)/2,
Jlk (q) =
O(q 2l−k(k+1)/2+ ) si l ≥ k(k + 1)/2.
Vinogradov demostró una versión débil de su conjetura, y es costumbre
referirse a ella y a cualquiera de las mejoras subsiguientes como “teorema del
valor medio de Vinogradov”.
La conjetura ha sido demostrada recientemente (en 2015) por Bourgain,
Demeter y Guth mediante técnicas del análisis harmónico. Aquí no vamos a
entrar en la demostración de este hecho, sino que nos limitaremos a probar una
versión débil que será suficiente para nuestros fines. No obstante, terminaremos
esta sección con algunas observaciones elementales sobre la conjetura.
pues
l/lk l/lk
O(q lk + )l/lk O(q lk + ) O(q lk + )
1
= = = O(1).
q l+ q lk +lk /l q lk + q (lk /l−1)
x1 + x2 − y3 = (y1 − x3 ) + (y2 − x3 ) + x3
En efecto, eh−i está formado por todos los monomios M de grado h − i con
variables de grado 1. Cuando uno de estos monomios se multiplica por una
potencia Xji , obtenemos uno de los monomios de rih si Xj no está en M , o bien
h
uno de los monomios de ri+1 si Xj está en M . Además, así se obtienen todos
h h
los monomios de ri + ri+1 sin repeticiones.
Para i = h tenemos trivialmente:
ph e0 = ph = rhh ,
luego
Jlk (q) ≤ k!q 2l−k .
Esto ya es un refinamiento de (7.2) para el caso k ≥ l, pero necesitamos mucho
más. Para ello probamos primero un resultado sobre congruencias:
x1 + · · · + xk ≡ m1 (mód p)
·····················
xk1 + · · · + xkk ≡ mk (mód pk )
T ≤ k!M k pk(k−1)/2 .
Demostración: Dada una solución del sistema, es claro que podemos desa-
rrollar
xj = xj1 + pxj2 + · · · + pk xj,k+1 ,
donde 0 ≤ xjs ≤ p − 1, para s = 1, . . . , k y 0 ≤ xj,k+1 ≤ M − 1. Por hipótesis
xj1 6= xj 0 1 , para j 6= j 0 , y además
Como x11 , . . . , xk1 son distintos dos a dos, a lo sumo uno de ellos es 0.
Podemos suponer que x11 , . . . , xk−1,1 son no nulos, y entonces
x11 ··· xk−1,1
.. ..
Y
= x11 · · · xk−1,1 (xj1 − xi1 ) 6= 0,
. .
k−1
xk−1
x ··· i<j
11 k−1,1
240 Capítulo 7. El método de Vinogradov
luego cada valor para xk2 determina una única solución del sistema x12 , . . . , xk2 ,
luego en total tiene T2 = p soluciones. A continuación consideramos el sistema
Entonces
2l−k(k+1)/2
Jlk (q) < 3l2k q12k p1 k
Jl−k (q1 ),
0 ≤ xj , yj ≤ p1 q1 − 1,
para h = 1, . . . , k, con 0 ≤ mj , nj ≤ p1 − 1 y 0 ≤ xj , yj ≤ q1 − 1.
7.2. El teorema del valor medio de Vinogradov 241
Diremos que una solución del sistema (7.4) es de primera clase si los con-
juntos {m1 , . . . , ml } y {n1 , . . . , nl } tienen al menos k elementos cada uno, y en
caso contrario diremos que es de segunda clase. Así, podemos descomponer
Jlk (p1 q1 ) = J1 + J2 ,
Por otra parte, el número de l-tuplas de números menores que p1 es pl1 , luego
pares de l-tuplas con el primer par restringido a l-tuplas con a lo sumo k elemen-
tos distintos es a lo sumo k l pl+k−1
1 /2, y el número de sumandos del integrando
anterior es a lo sumo k l p1l+k−1 . Aplicamos la desigualdad de Hölder:
P
U (m1 ) · · · U (ml )U (n1 ) · · · U (nl ) ≤
mj ,nj
P 1/2l P 1/2l
|U (m1 )|2l ··· |U (nl )|2l .
mj ,nj mj ,nj
1 −1
pP
Ahora acotamos |U (mj )|2l , |U (nj )|2l ≤ |U (m)|2l , con lo que
m=0
P
U (m1 ) · · · U (ml )U (n1 ) · · · U (nl ) ≤
mj ,nj
1 −1
pP 1 −1
pP
k l p1l+k−1 |U (m)|2l ≤ k l q12k pl+k−1
1 |U (m)|2l−2k ,
m=0 m=0
Por consiguiente
Z 1 −1
pX
J2 ≤ k l q12k pl+k−1
1 |U (m)|2l−2k da1 · · · dak
[0,1]k m=0
1 −1 Z
pX
= k l q12k p1l+k−1 |U (m)|2l−2k da1 · · · dak ,
m=0 [0,1]k
formas distintas para obtener (con repeticiones) todas las l-tuplas con al menos k
números distintos. Así pues,
2 Z
l X
J1 ≤ U (m1 ) · · · U (ml )U (n1 ) · · · U (nl ) da1 · · · dak
k [0,1]k m ,n j j
2 Z
l X X
= U (m1 ) · · · U (ml ) U (n1 ) · · · U (nl ) da1 · · · dak
k [0,1]k m n
j j
2 Z 2
l X
= U (m1 ) · · · U (ml ) da1 · · · dak ,
k
[0,1]k m j
A su vez,
2 Z 2 X 2
l X
J1 ≤ U (m1 ) · · · U (mk ) U (mk+1 ) · · · U (ml ) da1 · · · dak ,
k
[0,1]k m j m j
donde en el primer sumatorio mj recorre las k-tuplas con todas sus componentes
distintas dos a dos y en el segundo las l−k-tuplas arbitrarias. Por la desigualdad
de Hölder:
P 2 P
U (mk+1 ) · · · U (ml ) = U (mk+1 ) · · · U (ml )U (nk+1 ) · · · U (nl ) ≤
mj mj ,nj
P 1/(2l−2k) P 1/(2l−2k)
|U (mk+1 )|2l−2k ··· |U (nl )|2l−2k
mj ,nj mj ,nj
1 −1
pP 1/(2l−2k) 1 −1
pP 1/(2l−2k)
= p12l−2k−1 |U (m)|2l−2k · · · p12l−2k−1 |U (m)|2l−2k
m=0 m=0
1 −1
pP
= p2l−2k−1
1 |U (m)|2l−2k .
m=0
Por lo tanto,
2 Z 1 −1
2 pX
l 2l−2k−1
X
J1 ≤ p1 U (m1 ) · · · U (mk ) |U (m)|2l−2k da1 · · · dak .
k
[0,1]k m j m=0
k
donde, recordemos, Jl−k (q1 ; c) es el número de soluciones del sistema de ecuacio-
nes (7.3), que es nulo si algún |ch | ≥ lq1h . Además hemos probado la desigualdad
k k
Jl−k (q1 ; c) ≤ Jl−k (q1 ).
Finalmente, llamemos T al número de soluciones del sistema de congruencias
(m1 − m + p1 x1 )h + · · · + (mk − m + p1 xk )h
Ahora vamos a acotar T usando el teorema 7.6. Esto requiere que p1 > k,
que ya sabemos que se cumple. Por el mismo argumento empleado en la prueba
del teorema 7.2, vemos que T es también el número de soluciones del sistema
de congruencias
(m1 − m + p1 + p1 x1 )h + · · · + (mk − m + p1 + p1 xk )h ≡
Además los x0j no son congruentes módulo p1 . Así pues, fijados m, nj , yj , existen
k(k−1)/2
a lo sumo k!M k p1 soluciones posibles del sistema
x0h 0h
1 + · · · + xk ≡
M = E[q1 p−k+1
1 ] + 1 ≤ 2E[q1 p−k+1
1 ] ≤ 2q1 p−k+1
1 .
7.2. El teorema del valor medio de Vinogradov 245
Por lo tanto
k+1−k(k−1)/2 2k+1−k(k+1)/2
T ≤ k!2k q 2k p1 = k!2k q 2k p1
luego
2k+1−k(k−1)/2
Nl ≤ k!2k q12k p1 k
Jl−k (q1 )
y a su vez
2
l 2l−k(k+1)/2 k 2l−k(k+1)/2 k l2k 2k
J1 ≤ k!2k q12k p1 Jl−k (q1 ) = q12k p1 Jl−k (q1 )
k k!
2l−k(k+1)/2
≤ 2l2k q12k p1 k
Jl−k (q1 ).
Recapitulando,
2l−k(k+1)/2
Jlk (q) < Jlk (p1 q1 ) = J1 + J2 ≤ 3l2k q12k p1 k
Jl−k (q1 ).
Entonces
Jlk (q) ≤ (3l2k )r 34lr−k(k+1)r/2 q 2l−k(k+1)/2+δr ,
donde r
k(k + 1) 1
δr = 1− .
2 k
q q1 qr−1 q r r
p1 · · · pr < 2 2 ···2 = 2r < 2r qq −(1−1/k) < 3r q 1−(1−1/k) .
q1 q2 qr qr
k 2(l−kr)
A esto añadimos la cota trivial Jl−kr (qr ) ≤ qr , con lo que
r
Jlk (q) < (3l2k )r 3(2l−k(k+1)/2)r q (1−(1−1/k) )(2l−k(k+1)/2) 2l
qr . (7.5)
Usando que
r
qr2l < 32lr q (1−1/k) 2l
,
se llega a la desigualdad del enunciado.
En realidad necesitaremos únicamente esta versión más débil del teorema
anterior, según el cual la conjetura de Vinogradov se cumple para = 1/2 y
todo l suficientemente grande (respecto de un k prefijado):
7.2. El teorema del valor medio de Vinogradov 247
k ≥ 2, l ≥ k 2 + 4k 2 log k, q ≥ 2.
3k 2 − k + 1
= + 8k 2 log k > 0,
2
la conclusión es inmediata.
Supongamos ahora que q ≥ 2−k (2k)2k . Definimos la sucesión qu como en la
prueba del teorema anterior. Notemos que, mientras qu−1 ≥ 2−k (2k)2k ≥ 43 se
1/k
cumple que qu−1 ≥ 2k 2 y es posible definir pu ≥ k 2 y, a partir de él, qu , que
cumplirá qu ≤ qu−1 , pues
qu−1 qu−1
qu − 1 ≤ ≤ ≤ qu−1 − 1,
p1 2
ya que la última desigualdad equivale a qu−1 ≥ 2. Distinguimos dos casos:
1. Existe un número natural 1 ≤ r ≤ 4k log k tal que
4k log k log 3 + 8k 2 log k log l + 8lk log k log 3 + 2lk log 2 + 4lk log k
≤ ck log k + clk log2 k + clk log k + clk + clk log k ≤ clk log2 k,
luego se cumple también la desigualdad del enunciado.
248 Capítulo 7. El método de Vinogradov
luego
Gx < (N − 1)aδ + δ + 1 ≤ (N − 1)aδ + aδ + 1 = N aδ + 1.
Hemos probado esto bajo la hipótesis de que Gx > 0, pero es trivialmente cierto
si Gx = 0.
Ahora dividimos el intervalo [−W δ, W δ] en E[2W +1] subintervalos [xi , xi+1 ]
de longitud menor que δ. Si D(φ(n)) ≤ W δ, existe un entero h para el que
φ(n) − h ∈ [−W δ, W δ], luego existe un i tal que φ(n) − h ∈ [xi , xi+1 ], luego
luego
xi + h ≤ φ(n) < xi + h + δ.
Así pues, n está contado en Gxi , luego el número total de tales n es a lo sumo
el número de intervalos E[2W + 1] por la cota N aδ + 1 al número posible de
tales n para cada intervalo, lo que nos da la conclusión del teorema.
Demostración: Llamemos
PP
+Q
C= e2πiF (n)
n=P +1
Ahora
q PP
+Q
q P +Q−q+m
q
P P
2πiF (n)
P
2πiF (n)
P
|qC| = e ≤
e + q
m=1 n=P +1 m=1 n=P +1+m m=1
P +Q−q 1/2l
≤ Q1−1/2l |T (n)|2l + q2 .
P
(7.7)
n=P +1
Equivalentemente,
k
X F r) (n)
F (m + n) − F (n) = mr + 2λθq k+1 ,
r=1
r!
Por lo tanto,
Z
2l
m(Ωn )|T (n)| = |T (n)|2l da1 · · · dak ≤
Ωn
Z
22l |S(q)|2l da1 · · · dak + m(Ωn )(4πkλq k+2 )2l .
Ωn
4 En la primera desigualdad usamos que |eix1 − eix2 | ≤ |x1 − x2 |. Geométricamente, esto
equivale a que la longitud de una cuerda en una circunferencia es menor o igual que su arco.
Analíticamente equivale a que |eix − 1| ≤ |x|, que a su vez equivale a que cos x − 1 + x2 /2 ≥ 0,
lo cual se prueba fácilmente viendo que la función tiene un mínimo en 0.
5 Concretamente, usamos que
1 1 1 1
(x + y)n = ( 2x + 2y)n ≤ (2x)n + (2y)n ≤ (2x)n + (2y)n .
2 2 2 2
252 Capítulo 7. El método de Vinogradov
luego
Z
2l 2l k(k+1)/2 k+1 −k
|T (n)| ≤ 2 q (λq ) |S(q)|2l da1 · · · dak + (4πkλq k+2 )2l .
Ωn
pues Ωn puede cubrirse por un número finito de trasladados del cubo unitario I k ,
luego la integral sobre Ωn es la suma de las integrales sobre los trasladados
(g + I k ) ∩ Ωn , que coinciden con las integrales sobre I k ∩ (−g + Ωn ), cuya suma
es la integral sobre Ω0n . Así pues,
Z
|T (n)|2l ≤ 22l q k(k+1)/2 (λq k+1 )−k |S(q)|2l χΩ0n da1 · · · dak + (4πkλq k+2 )2l ,
[0,1]k
y a su vez
P +Q−q
X Z P +Q−q
|T (n)|2l ≤ 22l q k(k+1)/2 (λq k+1 )−k |S(q)|2l
P
χΩ0n da1 · · · dak
n=P +1 [0,1]k n=P +1
P +Q−q
Por consiguiente, si a ∈ Ω0n0 , entonces
P
χΩ0n (a) es menor o igual que el
n=P +1
número de números n en [P + 1, P + Q − q] tales que D(Ak (n) − Ak (n0 )) ≤ λq.
Vamos a probar que este número es menor o igual que 4kq, para lo cual
aplicamos el teorema anterior con
φ(n) = Ak (n) − Ak (n0 ), M = P + 1, N = Q − q > 1, a = 2, δ = λ(k + 1),
W = q/(k + 1) (si Q − q = 1 el intervalo es {P + 1} y la conclusión es trivial).
Observemos que
1 k) F k+1) (x)
φ(n + 1) − φ(n) = (F (n + 1) − F k) (n)) = ,
k! k!
para cierto n < x < n + 1. La hipótesis del teorema implica que F k+1) no se
anula, luego no puede cambiar de signo en el intervalo que estamos considerando.
Observemos además que ni las hipótesis ni la conclusión del teorema que estamos
probando se alteran si cambiamos F por −F , por lo que no perdemos generalidad
si suponemos que el signo es positivo. Esto garantiza que φ cumple las hipótesis
del teorema anterior. La conclusión es que el número de números n que estamos
estudiando es menor o igual que
2q 2q q
(2λ(k + 1)Q + 1)( + 1) ≤ (2k + 3)( + 1) ≤ 3k + q < 4kq,
k+1 k+1 4
donde hemos usado que Q ≤ 1/λ y luego que k ≥ 7.
Así pues,
P +Q−q
|T (n)|2l ≤ 22l q k(k+1)/2 (λq k+1 )−k 4kqJlk (q) + (Q − q)(4πkλq k+2 )2l ,
P
n=P +1
+4Qπkλq k+1 + q
2
1/2l
≤ e(c/2)k log k 2Q1−1/2l q 3/2−k(k+1) 4kλ−k + 4Qπkλq k+1 + q.
254 Capítulo 7. El método de Vinogradov
k = E[log τ / log a] + 1.
con lo que
ak+1
Q< ≤ Q2 . (7.9)
τ
Consideramos la función F : [a, 2a] −→ R dada por
τ log x
F (x) = − .
2π
Así
k!τ
F k+1) (x) = (−1)k+1 .
2πxk+1
En [a, 2a] se cumple que
k+1)
τ F (x) τ
≤
≤ . (7.10)
2π(k + 1)(2a)k+1 (k + 1)! 2π(k + 1)ak+1
el extremo derecho del intervalo es el punto en el que el miembro central se reduce a la mitad.
Ése es el valor de λ para el intervalo siguiente, y, como el miembro izquierdo se obtiene del
derecho dividiéndolo entre 2k+1 , el número de intervalos necesarios para superar b será a lo
sumo k + 1.
256 Capítulo 7. El método de Vinogradov
Por la observación tras el teorema 6.5, para 0 < σ < 1 se cumple que
X
1 1 2 2
2Kkeck log k a1−ρ = O(keck log k a1−ρ−σ ),
s
≤ σ
(7.11)
a<n≤b n E[a]
2 log τ 2 log1/4 τ
k< <
log a A (log log τ )3/4
2 log τ
k< < log1/2 τ,
log a
1
luego log k < 2 log log τ y, usando de nuevo (7.13),
1/3
1 1 1
k log k < (log τ log log τ )1/4 < log a
A A A
1 1
= log1/3 a < 4/3 1/3 log1/3 τ,
A4/3 A 6
luego eligiendo A suficientemente grande se cumple (7.12).
En estas condiciones (7.11) nos da que
X 1
log τ 1−σ ck log2 k− 70klog a
2 log k
=O a e
ns log a
a<n≤b
3 −1) log a
log τ 1−σ (70c
=O a e 70k2 log k .
log a
258 Capítulo 7. El método de Vinogradov
y así
−B log a(log log τ )1/2
X 1 1−σ 1/2 τ
= O log τ a e log , (7.15)
ns
a<n≤b
1/4
τ (log log τ )5/4 log1/2 τ 1/4 5/4
≤ eA1 A2 log 1/2
≤ eA3 log τ (log log τ ) , (7.17)
A2 (log log τ )
7.3. Aplicación a la función dseta 259
1/2
≤ ≤ elog τ ≤ elog τ (log log τ ) ,
A2 (log log τ ) A2
1
A(log τ log log τ )3/4 < log(2u α) < log τ.
6
La segunda desigualdad se cumple si r ≥ 12, pues entonces
2 1
log(2u α) < log τ ≤ log τ,
r 6
mientras que la primera se cumple si A < A1 , pues entonces
O log τ e log1/2 τ .
260 Capítulo 7. El método de Vinogradov
se cumple que
1/4
τ (log log τ )5/4
ζ(s) = O(eA3 log ).
Observamos que
A
1− 3/4
≤ σ, τ ≥ τ0 .
log τ (log log τ )3/4
7.3. Aplicación a la función dseta 261
Observaciones Una versión más débil del teorema anterior afirma que, para
todo > 0, existe un τ1 tal que ζ(s) no tiene ceros en la región
1
σ >1− 3/4+
, τ ≥ τ1 .
log τ
En efecto, basta observar que
1 A
1− >1−
log3/4+ τ log3/4 τ (log log τ )3/4
para todo τ suficientemente grande, ya que esto equivale a que
A su vez, el teorema 5.20 nos proporciona una estimación mejor del error en
el teorema de los números primos. Concretamente, ahora podemos tomar
A
η(t) = 3/4
log t (log log t)3/4
(extendida con un valor constante para t < τ1 ). La conclusión es:
262 Capítulo 7. El método de Vinogradov
pues el miembro izquierdo, como función de t, decrece hacia −∞, luego pasa
por 0 una única vez. Entonces, para todo t ≥ 1, se cumple que
ω(x) ≥ log t0 ,
pero
A log x
log t0 = 3/4
,
log t0 (log log t0 )3/4
luego
(A log x)4/7
log t0 = ,
(log log t0 )3/7
luego
3 4 4
log log t0 < log log t0 + log log log t0 = log log x + log A < c log log x,
7 7 7
para x suficientemente grande, luego
(A log x)4/7
ω(x) ≥ ≥ c log4/7 x (log log x)−3/7 .
(log log t0 )3/7
Números compuestos
60
10
48
36
24
12
5 6
263
264 Capítulo 8. Números compuestos
e e2 + 3 e5 52
d(4n/5) = d(2e2 +1 3e3 5e5 −1 · · · prpr ) = d(n) ≥ d(n) > d(n),
e2 + 1 e5 + 1 33
contradicción.
Ahora consideramos
pr+1 n e −1 e2 e pr 32
d( ) = d(2e2 −1 3e3 · · · prpr pr+1 ) = 2d(n) ≥ 2d(n) = d(n),
2pr e2 + 1 epr +1 43
1 Aquí usamos que la sucesión n/(n + 1) = 1 − 1/(n + 1) es creciente.
8.1. Números altamente compuestos 265
pero
pr+1 n
< n,
2pr
ya que esto equivale a que pr+1 < 2pr , y esto es el postulado de Bertrand
(teorema 3.9).
El teorema anterior se expresa de forma más conveniente en términos de
primoriales.
N = P1 # · · · Pk #,
N= 2 · 3 · · · · · · · · · P1
2 · 3 · · · · · · P2
(8.1)
·········
2 · 3 · · · Pk
N = P1 # · · · Pk #,
Hay que tener presente que todas estas condiciones son necesarias, pero no
suficientes, para que un número N sea altamente compuesto. Podemos dar
algunas más:
e
Teorema 8.5 Si N = 2e2 · · · P1 P1 es un número altamente compuesto, entonces
ep
pi i < P16 , para todo primo pi ≤ P1 .
Nota Hemos probado que la condición N 0 < N se cumple cuando P1 > 4 096,
pero es fácil comprobar con un ordenador que también se cumple3 para los
primos en el rango 37 ≤ P1 < 4 096, por lo que el teorema vale, de hecho,
para P1 ≥ 37. Más aún, una comprobación directa sobre los números altamente
compuestos con P1 < 37 muestra que el teorema se cumple de hecho para
P1 ≥ 11, es decir, a partir de N = 55 440 = 24 · 32 · 5 · 7 · 11.
Similarmente podemos probar que, los números altamente compuestos sufi-
cientemente grandes tienen divisores primos con exponente 2:
Teorema
√ 8.8 Si √ N es un número altamente compuesto suficientemente grande
y 2 P1 < p < 4 P1 , entonces ep = 2.
√
Demostración: Supongamos que p > 2 P1 y que ep ≥ 3. Entonces,
con la misma notación de la prueba del teorema anterior, esta vez tomamos
N 0 = N P1∗ P1∗∗ /(p1 p2 p3 p4 ), que cumple N 0 < N cuando N es suficientemente
grande. Además
4
d(N 0 ) 4ep1 ep2 ep3 ep4 3
= ≥4 > 1,
d(N ) (ep1 + 1)(ep2 + 1)(ep3 + 1)(ep4 + 1) 4
contradicción. √
Sea ahora p el menor primo para el que ep = 1 y supongamos que p < 4 P1 .
Tomando N suficientemente grande, podemos √ asegurar que haya al menos 11
primos comprendidos estrictamente entre 4 P1 y P1 .
En efecto, por el teorema de los números primos, dado > 0, tomando P1
suficientemente grande podemos exigir que
√ √
π(P1 ) log P1 π(4 P1 ) log 4 P1
1−< , √ < 1 + ,
P1 4 P1
luego √
P1 p 4 P1
π(P1 ) > (1 − ) , π(4 P1 ) < (1 + ) √ .
log P1 log 4 P1
Por lo tanto,
√
p P1 4 P1
π(P1 ) − π(4 P1 ) > (1 − ) − (1 + ) √
log P1 log 4 P1
√ √
P1 4 P1 (1 − )P1 − (1 + )8 P1
≥ (1 − ) − (1 + ) √ = ,
log P1 log P1 log P1
y aplicando la regla de L’Hôpital vemos que el último término tiende a +∞
cuando P1 tiende a +∞. √ Por lo tanto, si P1 es suficientemente grande, podemos
exigir que π(P1 ) − π(4 P1 ) > 12.
3 Más concretamente, se trata de comprobar que si 5 ≤ P1 < 4 096 (hay 560 primos en
estas condiciones) y p es el menor primo p > P1 /2, entonces P1∗ P1∗∗ /p1 p2 p3 > 1 siempre que
P1 ≥ 37.
268 Capítulo 8. Números compuestos
La elección√de P1 garantiza que los 11 primos son distintos entre sí. Ahora
pi+t ≤ 2t+2 P1 y pj+t ≥ P1 /2t , luego
√ 7
N0 P1 K
≤K = √ <1
N P14 P1
cuando P1 es suficientemente grande. Por otra parte,
7
d(N 0 ) 3 1
= > 1,
d(N ) 2 24
pues todos los primos considerados tienen exponente 1 en N , luego los del nume-
rador dan lugar a factores 3/2 y los del denominador a factores 2. Así tenemos
nuevamente una contradicción.
Nota Las pruebas de los teoremas anteriores pueden generalizarse para probar
que un número altamente compuesto suficientemente grande es divisible entre
primos con cualquier exponente en un conjunto finito prefijado.
Hemos señalado que 7! es altamente compuesto, y es fácil ver que lo mismo
sucede con todos los factoriales precedentes, pero no con los posteriores:
Esto equivale a su vez a que (e2 −3)(e13 +2) > (e2 +1)(e13 +1), que se simplifica
hasta e2 > 4e13 + 7. Según (3.1) tenemos que
E[log
P2 n] E[log
P2 n]
e2 = E[n/2m ] > n/2m − 1 = n − n/2E[log2 n] − E[log2 n]
m=1 m=1
≥ n − 2 − log2 n,
donde hemos usado que n/2E[log2 n] ≤ n/2log2 n = 1 y que E[logn 2] ≤ log2 n + 1.
Similarmente,
E[log
P13 n] ∞
E[n/13m ] ≤ n/13m = n/12.
P
e13 =
m=1 m=1
8.2. Números altamente compuestos superiores 269
1.0
0.5
0 20 40 60 80 100.
luego n
d(n) ≤ d(N ) < d(N ).
N
e
Supongamos ahora N = 2e2 · · · prpr es altamente compuesto superior (res-
pecto de ) y consideremos un primo p | N . Sea N 0 = N/p, de modo que
d(N 0 ) d(N )
≤ . (8.2)
N 0 N
Usando que d es multiplicativa, esto se simplifica a
ep ep + 1
≤ ,
p −1)
(e p p ep
que equivale a
1
p ≤ 1 +
ep
o también a
log((ep + 1)/ep )
≤ (8.3)
log p
Similarmente, si p es un primo cualquiera y consideramos N 0 = N p, la
condición
d(N 0 ) d(N )
0
≤
N N
se traduce ahora en que
ep + 2 ep + 1
(e +1)
≤ ep ,
p p p
de donde
ep + 2
p ≥
ep + 1
o, equivalentemente,
log((ep + 2)/(ep + 1))
≥ . (8.4)
log p
Ahora necesitamos distinguir casos según que puedan darse o no las igual-
dades en las fórmulas precedentes. Para ello conviene introducir la notación
siguiente:
Ep = {F (p, e) | e = 1, 2, 3, . . .} .
S
Llamamos E = Ep .
p
8.2. Números altamente compuestos superiores 271
Teorema 8.12 Para cada número real > 0 tal que ∈/ E, la función d(n)/n
alcanza su máximo en un único número NQ
, que es, por definición, altamente
compuesto superior. Concretamente N = pep () , con
p
1
ep () = E .
p − 1
Demostración: Sea N un número donde d(n)/n alcance su valor máximo
(en principio puede haber más de uno). Si p divide a N con exponente ep , hemos
probado que
F (p, ep + 1) < < F (p, ep ), (8.5)
es decir,
log((ep + 2)/(ep + 1)) log((ep + 1)/ep )
<<
log p log p
o también
ep + 2 ep + 1
< p < ,
ep + 1 ep
de donde se llega fácilmente hasta
1 1
− 1 < ep < .
p −1 p −1
Esto implica en particular que (p − 1)−1 no es entero, y que, necesariamente,
e = E[(p − 1)−1 ], luego N está completamente determinado por .
Notemos ahora que, como la función F (p, e) decrece hacia 0 tanto si au-
menta p como si aumenta e, para cada δ > 0 hay un número finito de elementos
en E mayores que δ, lo que se traduce en que E puede ordenarse en una sucesión
decreciente4
1 > 2 > 3 > · · ·
convergente a 0. Definimos 0 = +∞.
Así, si > 0 no está en E, existe un único i tal que i < < i−1 , y se
cumple que cualquier par de números en ]i , i−1 [ tienen por debajo y por
encima exactamente los mismos elementos de E. A su vez, esto implica que
ambos cumplen las desigualdades (8.5) con los mismos valores de p y ep , por lo
que N es el mismo para ambos.
calculando F (2, 1), F (2, 2), F (2, 3), . . . hasta obtener un valor menor que δ, luego calculamos
F (3, 1), F (3, 2), . . . hasta obtener de nuevo un valor menor que δ, y seguimos así hasta que
F (p, 1) < δ. Por último ordenamos todos los números obtenidos.
272 Capítulo 8. Números compuestos
La tabla siguiente contiene los i > 0.3 junto con sus números altamente
compuestos superiores asociados. No obstante, cabe señalar que se ha impuesto
el convenio de no considerar al 1 como número altamente compuesto superior
(aunque sí que se considera altamente compuesto).
F (q, eq + 1) ≤ i ≤ F (q, eq ),
En particular, vemos que Ni−1 y Ni son divisibles entre los mismos pri-
mos con los mismos exponentes, salvo en el caso de q, que divide a Ni−1 con
exponente e − 1 y a Ni con exponente e. Por consiguiente, Ni = Ni−1 q.
Si partimos de (8.2) aplicado a N = Ni , N 0 = Ni /q = p, entonces (8.3) es la
igualdad i = F (p, eq ), luego (8.2) también es una igualdad, luego concluimos
que Ni−1 y Ni son ambos máximos de la función d(n)/ni .
Nos falta considerar la posibilidad de que un mismo i esté en varios conjun-
tos Ep . Observemos que i ∈ Ep equivale a que existe un e tal que i = F (p, e),
lo cual equivale a su vez a que pi = 1 + 1/e.
Si i ∈ Ep ∩ Eq , tenemos que pi y q i son números racionales. La conjetura
de las cuatro exponenciales (véase la página 12) implica que esto sólo es posible
si i es entero, lo cual a su vez sería absurdo porque 1 + 1/e no puede ser una
potencia de primo.
Así pues, la conjetura de las cuatro exponenciales implica que los conjuntos
Ep son disjuntos dos a dos, y el teorema anterior recoge todas las posibilidades
de generación de números altamente compuestos superiores.
Si no nos apoyamos en la conjetura, el teorema 1.7 nos asegura al menos
que un mismo i puede pertenecer a lo sumo a dos conjuntos Ep y Eq . En tal
caso, los mismos razonamientos del teorema anterior implican que la función
d(n)/ni alcanza su máximo en cuatro números, que son Ni−1 , Ni = Ni−1 pq,
Ni−1 p, Ni−1 q, con lo que hay dos números altamente compuestos superiores que
no aparecen en la sucesión de los Ni .
En particular vemos que una consecuencia de la conjetura de las cuatro ex-
ponenciales es que el cociente de dos números altamente compuestos superiores
consecutivos es un número primo. Sin la conjetura, sólo podemos afirmar que
los cocientes Ni /Ni−1 son primos o productos de dos primos, y en el segundo
caso (para los primos p < q), hay dos números altamente compuestos compren-
didos entre Ni−1 y Ni , a saber, Ni−1 < Ni−1 p < Ni−1 q < Ni . Sin embargo, es
muy poco probable que pueda darse este caso.
Una última observación es que si al calcular los i recordamos con qué valores
F (p, e) se genera cada uno, el cálculo de los Ni se reduce a ir multiplicando los
primos p. Por ejemplo, 1 = F (2, 1) ∈ E2 , 2 = F (3, 1) ∈ E3 , 3 = F (2, 2) ∈ E2 ,
4 = F (5, 1) ∈ E5 , etc., lo que se traduce en que la sucesión de los Ni es
1, 2, 2 · 3, 2 · 3 · 2, 2 · 3 · 2 · 5, ...
1.5
1.0
0.5
12, 18, 20, 24, 30, 36, 40, 42, 48, 54, 56, 60, 66, 70, 72, 78, 80, 84, 88, 90, 96, 100 . . .
Parece que los números abundantes tengan que ser pares, pero no es así. El
primer número abundante impar es 945. Otra conjetura falsa sobre los núme-
ros superabundantes es que son múltiplos de la suma de sus cifras. El primer
contraejemplo es 149 602 080 797 769 600.
Pese a que, como hemos indicado, los números superabundantes no tienen
por qué ser altamente compuestos ni viceversa, ambas clases de números com-
parten muchas propiedades. Por ejemplo:
q eq +2 − 1 q ep +2 − 1 pep +1 − 1
≥ ≥ ,
q eq +2 − q q ep +2 − q pep +1 − p
con lo que tenemos una contradicción.
Supongamos ahora que er ≥ 2. Si r = 1, es decir, si n = 2e , consideramos
n = n3/4 < n, y al simplificar σ(n0 )/n0 < σ(n)/n obtenemos
0
Operando llegamos a que 2e2 +1 < 79/7 ≈ 11.3, luego e3 ≤ 2, lo que a su vez
implica que e2 = e3 = 2 y n = 36. Ahora supongamos que r > 2, luego q ≥ 3,
y acotamos
1 q−1 p−1 1 1
1+ < 1+ 3 1+ 3 = 1+ 1+
t q −q p −p q(q + 1) p(p + 1)
1 1
< 1+ 2 1+ 2 .
q p
Por otro lado, la desigualdad
1 1 1
1+ 2 1+ 2 ≤1+
q p t
equivale a
1 1 1 1
+ 2+ 2 2 < ,
q2 p q p t
pero
1 1 1 2q 2 + 1 1 1
2
+ 2+ 2 2 ≤ , < ,
q p q p q4 4q t
la última desigualdad por el postulado de Bertrand. Luego si
2q 2 + 1 1
3
≤
q 4
e
Teorema 8.16 Si n = 2e2 · 3e3 · · · prpr es superabundante y 1 ≤ i < j ≤ r,
entonces epj toma uno de los valores k−1, k, k+1, donde k = E[epi log pi / log pj ].
ep epi +1
< pi i (pi pj − pi ) + pm
i < pi (pj − 1) + pm
i .
epj +1 epi +m
pj (pm−1
i − 1) + pj < pi (pj − 1) + pm−1
i .
lo que equivale a
2 m−1
pm
i (pj − 1) ≤ pj (pi − 1).
Si llamamos P = pm−1
i y p = pi , tenemos que P ≥ p ≥ 2 y P < pj < pP , y se
trata de probar que la función
F (t) = (P − 1)t2 − pP t + pP
F (P ) = (P − 1)P 2 − pP 2 + pP = P (P − 1)(P − p) ≥ 0.
σ(2n)/2n 2e2 +2 − 1
= e2 +2 > 1,
σ(n)/n 2 −2
luego el menor número superabundante mayor que n es menor o igual que 2n.
pues el primer factor está acotado. Por lo tanto, la sucesión σ(n)/n1+ alcanza
un valor máximo.
1.5
1.0
0.5
0 20 40 60 80 100
para p ≥ 2 y e > 0, con el convenio de que F (p, 0) = +∞. Más aún, observemos
que
log(1 + p(pp−1
e −1) )
1
log(1 + p+p2 +···+p e)
F (p, e) = = .
log p log p
La última expresión muestra que la función p 7→ F (p, e) es estrictamente de-
creciente, y una comprobación rutinaria nos da que biyecta el intervalo ]1, +∞[
con ]0, +∞[. También es claro que la sucesión e 7→ F (p, e) decrece hacia 0, para
e = 0, 1, 2, 3 . . .
Ep = {F (p, e) | e ≥ 1},
S
y llamamos E = Ep . Como F (p, e) decrece cuando crece cualquiera de sus
p
argumentos, es claro que los elementos de E forman una sucesión decreciente
1 > 2 > 3 > . . . convergente a 0. Definimos además 0 = +∞.
Observemos que si ∈ Ep entonces
1
p = 1 + ∈ Q,
p + p2 + · · · + pe
p1+ −1 p1+ −1
log p −1 log p −1
− 1 < ep − 1 < .
log p log p
Vemos que coinciden con los primeros números altamente compuestos supe-
riores, sin embargo, la coincidencia termina en N15 , pues el decimosexto número
colosalmente abundante es
N16 = 25 · 33 · 52 · 7 · 11 · 13 · 17 · 19 · 23,
26 · 33 · 52 · 7 · 11 · 13 · 17 · 19.
Bibliografía
[1] Alaoglu, L., y Erdős, P., On Highly Composite and Similar Numbers.
Transactions of the American Mathematical Society, 56 (3) (1944) 448–
469.
[6] Edwards, H.M. Riemann’s zeta function. Academic Press, New York, 1975.
[13] Ramanujan, S. Highly composite numbers. Proc. London Math. Soc. 14 (2)
(1915) 347–409
283
284 BIBLIOGRAFÍA
integral logarítmica, xi
285









