Diapositiva 1

1
Números primos y Benford generalizado
Bartolo Luque Lucas Lacasa ETSI
Aeronáuticos Dpto. Matemática Aplicada y
Estadística Universidad Politécnica de
Madrid Octavio Miramontes Instituto de
Física UNAM
Madrid 23 de Febrero de 2006
2
Números primos
"primo" "de base"
No parece que haya una regla que parezca gobernar
la sucesión de los números primos.
Para comprobar si un número x es primo, basta con
probar que ningún número y menor o igual a vx es
divisor de x.
God may not play dice with the universe, but
something strange is going on with the prime
numbers. Paul Erdös
3
It is evident that the primes are randomly
distributed but, unfortunately we don't know
what 'random' means. R.C. Vaughan
4
The Counting Prime Function
5
La distribución de números primos parece ser
aleatoria. Por ejemplo, hay probablemente
infinitos primos gemelos y existen gaps
arbitrariamente largos entre primos.
6
Sin embargo, la función p(x) exhibe un
sorprendente "buen comportamiento".
7
"For me, the smoothness with which this curve
climbs is one of the most astonishing facts in
mathematics." Don Zagier, "The first 50 million
primes" Mathematical Intelligencer, 0 (1977)
1-19
8
Si contamos el número de enteros pares inferiores
a n, tenemos que hay exactamente n/2 si n es
par y (n1)/2 si n es impar De modo que cuando
n tiende a infinito la densidad de números pares
tiende a 1/2. Decimos entonces que los números
pares tienen una densidad límite igual a 1/2.
Exactamente los mismo pasa para los impares. Y
los múltiplos de 3 tienen densidad límite 1/3.
Los números no múltiplos de 3 tienen densidad
límite 2/3. Y en general los múltiplos de r gt 0
tendrán densidad límite 1/r y los no múltiplos
de r (1 - 1/r).
9
Dificultades de la densidad límite
(1) Un conjunto dado de números no tiene por qué
tener siempre una densidad límite. Consideremos
como ejemplo el conjunto D 1, 4,5,6,7,
16,17,18,...,31, 64,65,... obtenido
reagrupando todos los enteros m comprendidos
entre 22N (incluido) y 22N1 (excluido) para N
0, 1, 2, 3, ... Están incluidos 20, 21) 22,
23) 24, 25) 26, 27)... Están excluidos (21,
22 (23, 24 (25, 26... La densidad de los
elementos de D inferiores a un número m será
10
Dificultades de la densidad límite
Pero esta densidad no posee límite para m 22N
es 3/8 y para m 22N1 es 5/8. Oscila
continuamente...
11
Otra aparente dificultad es que ciertos conjuntos
puedan tener densidad límite nula. No se trata de
una paradoja y es fácil encontrar ejemplos sean
los cuadrados 0, 1, 4, 9, 16, ... , n2. El
número de cuadrados inferiores a un número m es
exactamente el número de enteros inferiores a
vm, número que es inferior a vm 1. La densidad
de cuadrados es inferior entonces a (vm 1)/m
que tiende a cero cuando m tiende a
infinito. Decimos que un conjunto infinito que
tiene una densidad límite nula es un conjunto que
se rarifica. Y es el caso de la densidad de los
números primos. Se trata de un conjunto que
aunque es infinito se hace cada vez más y más
tenue, y su densidad límite se hace nula.
12
Teorema de rarefacción de Legendre
El conjunto de los números primos admite una
densidad límite nula.
13
Teorema de rarefacción de Hadamard-Poussin
Existe una función simple f(x) tal que
Los dos teoremas de rarefacción anteriores nos
dicen que si existe una función tal, tendrá que
cumplir que f(x)/x tiene límite cero cuando x
tiende a infinito y la suma de las inversas de
f(x) será infinita.
14
C.F. Gauss was the first to note, towards the
middle of the19th century, a certain regularity
in the series, in that the average distance
between consecutive primes is about Ln(n) 1.
1 D. Zagier, The first 50 million prime
numbers, Math. Intell. (1977) 719.
In 1963 Ulam found that when arranging the
integers in a grid, the primes form a spiral 2.
2 M.L. Stein, S.M. Ulam, M.B. Wells, A visual
display of some properties of the distribution of
primes, Am. Math. Mon. 71 (1964) 516520.
15
A la edad de 15 años Gauss, observando los
valores de tablas de primos, propuso que
alrededor de un número x la proporción de primos
era aproximadamente 1/Ln x. Que la densidad de
los números primos cercanos a x es alrededor de
1/Ln x. O sea que tomando un número alrededor de
x la probabilidad de que sea primo es 1/Ln
x. Por ejemplo, alrededor del número 1000 habrá
1/ Ln 1000 14,7 de números primos, alrededor
de 1000000 7,23 Esto condujo a Gauss a
proponer la función Li(x) como aproximación de
?(x).
16
Versión estadística de la distribución de números
primos
Sucesión de urnas con i bolas cada una.
De las i bolas i / Ln i bolas son negras y i -
i / Ln i bolas son rojas.
Posición de los primos
17
El teorema de los números primos
En 1896, de la Valee Poussin y Hadamard probaron
simultáneamente lo que se había sospechado
durante mucho tiempo, el teorema de los números
primos
El número de primos que no excede a x es
asintótico a x/log x. En otras palabras, la
probabilidad de que un número x escogido al azar
sea primo es 1/log x.
18
El teorema nos dice que x/log x es, hasta cierto
punto, una buena aproximación a p(x) . Al decir
que "a(x) es asintótico a b(x)" o "a(x) b(x)"
decimos que el límite de a(x)/b(x) es 1 cuando x
tiende a infinito. Pero, observemos que
a(x)  b(x) no significa que a(x) - b(x) sea
pequeño.
19
The logarithmic integral function Li(x)
Zagier en su artículo dice al respecto "within
the accuracy of our picture, the two coincide
exactly."
20
Li(x) no tiene primitiva analítica simple. Pero
podemos aproximar Li(x) por una función gracias a
un teorema de análisis que nos dice que dos
funciones para las cuales sus derivadas se
comportan de la misma manera en el infinito
implica que las funciones mismas se comportan
igual en el infinito. Busquemos entonces una
función cuya derivada se comporte en el infinito
como la derivada de Li(x) que es de hecho 1/Ln x.
Por ejemplo x/Lnx ya que su derivada es 1/Lnx
- 1/Ln2x se comporta como 1/Lnx en el infinito.
21
Supongamos que la densidad de primos f(x) va como
Su derivada es
Luego no funciona bien..... Sin embargo si
tomamos como densidad
22
Esta función actúa como x/Lnx
porque su derivada se comporta como 1/Ln x.
23
Supongamos
24
Tchebychev
25
Sea L la llamada constante de Liouville-Erdös.
Puede demostrarse que todo primo es
26
B. Riemann
Ueber die Anzahl der Primzahlem unter einter
gegebenen Grösse (1859)
27
La función de Riemann R(x)
Los números de Möbius se definen como cero
cuando n es divisible por un cuadrado y como
(-1)k en caso contrario. Donde k es el número de
distintos factores primos de n.
28
Recordemos que p(x) no es una función escalonada.
Sin embargo la función Riemann R(x) es una
aproximación suave. Podemos interpretar a R(x)
como una densidad media de los números primos. E
interpretar entonces a la diferencia R(x) - p(x)
como fluctuaciones locales.
29
Are Prime Numbers Regularly Ordered? Z. Gamba, J.
Hernando and L. Romanelli Physics Letters A 145,
no. 2,3 (2 April 1990), 106-108.
Consideran R(x) - p(x) como una "señal" y
calculan sus exponentes de Lyapunov para decidir
si el "mecanismo" que la crea es caótico o no.
Los autores concluyen "...a regular pattern
describing the prime number distribution cannot
be found. Also, from a physical point of view,
we can say that any physical system whose
dynamics is unknown but isomorphic to the prime
number distribution has a chaotic behaviour."
30
La función zeta ?(s)
Euler la llamó función zeta en 1737. Consideró
que s era un real mayor que 1.
31
Repitamos la operación para el siguiente primo 3.
32
Producto de Euler para la función zeta.
33
La función zeta ?(s) de Riemann
Bernhard Riemann hacia 1859 generalizó la función
zeta a números s x iy complejos. Aquí vemos
una representación gráfica del módulo de la
función z de Riemman ? (s). Obsérvese el polo
en s 1.
34
Aquí vemos una representación gráfica del módulo
de la inversa de la función z de Riemman 1/?
(s). De este modo podemos ver fácilmente los
ceros de la función z como polos. Los ceros
parece que vayan paralelos y cercanos al eje
imaginario.
35
Hipótesis de Riemann
Los pares negativos (-2, -4, -6, etc.) son ceros
de la función zeta, los llamados ceros triviales.
La hipótesis de Riemman afirma que todos los
demás ceros, llamados no triviales, tienen parte
real igual a ½. Es decir, que son de la forma ½
iy.
Grafica de y frente al módulo
36
Los 10.000 primeros millones de ceros de la
función zeta están en la línea crítica ½ (Mayo
2002).
37
La función zeta de Riemman está profundamente
conectada con la distribución de los números
primos. Helge von Koch probó en 1901 que la
hipótesis de Riemann es equivalente a
La distribución de ceros no triviales de la
función zeta de Riemann sobre la línea crítica
es "la dual" de la distribución de los números
primos.
38
El teorema de los números primos (II)
De forma más precisa, el teorema de los números
primos dice que
para alguna constante positiva a y
El error está de hecho íntimamente conectado con
la hipótesis de Riemann.
39
La distribución de ceros no triviales de la
función zeta de Riemann sobre la línea crítica
es "la dual" de la distribución de los números
primos.
La función zeta de Riemman está profundamente
conectada con la distribución de los números
primos. Helge von Koch probó en 1901 que la
hipótesis de Riemann es equivalente a
En otras palabras la hipótesis de Riemann nos
dice que Li(x) es una buena aproximación.
40
(No Transcript)
41
(No Transcript)
42
DISCUSSION ON BENFORDS LAW AND ITS
APPLICATION LI ZHIPENG, CONG LIN AND WANG HUAJIA,
Oct 2004
3.3.2 The Prime-number Series The prime number
series is rather uniform below 100000, with the
probability of each possible first significant
digit being between 12.5 and 10.4. Moreover,
using the upper and lower bounds of function pi
from the prime number theorem, it can be shown
that the prime number sequence approximates a
uniform distribution.
43
Si los números primos están distribuidos "al
azar", entonces para un número suficientemente
grande de ellos, la distribución de sus primeros
dígitos significativos (1,2,39) debería ser
igual a 100/9 11,1111111111.
44
(No Transcript)
45
Existe alguna suerte de misteriosa dinámica bajo
la distribución de los números primos? Es la
"forma" familiar de la secuencia de primos el
resultado de algún tipo de dinámica o proceso
evolutivo? I. V. Volovich
Versión estadística de la distribución de
números primos? dinámica
46
Normalicemos P(d)
47
(No Transcript)
48
Estamos trabajando en el intervalo 1, 10D1),
es decir en D 1 décadas. Normalicemos ahora
la función densidad de probabilidad de primos
De modo que la función densidad de primos es
Comprobemos que esta densidad nos da directamente
una ley de Benford generalizada normalizada
49
(No Transcript)
50
Estamos trabajando en el intervalo 1, 10D1),
es decir en D 1 décadas. Supongamos que la
función es
De modo que la función densidad de primos es
Comprobemos que esta densidad nos da directamente
una ley de Benford generalizada normalizada
51
(No Transcript)
52
Qué significa esta densidad de probabilidad?
que un número x tomado del intervalo aparecerá
con probabilidad f(x) en mi lista de números que
cumple la ley de Benford generalizada. Si
recorriera el intervalo número a número y lanzara
un dado cargado con la anterior probabilidad
para determinar si el número es o no primo,
cuál sería el número medio de primos?
53
Si pasamos de décadas y hablamos de N, estamos
diciendo que
Tomando como dependencia para alfa
54
Ecuación auto-consistente
Supongamos
55
Ecuación auto-consistente
56
(No Transcript)
57
Ecuación auto-consistente
Supongamos
58
(No Transcript)
59
Pierre Dusart Dusart99 demostró que si x gt 598
entonces
(The upper bound holds for all x gt 1.) This
gives a tight bound for larger x. Note x/log x lt
pi(x) for x gt 10.
P. Dusart, "The kth prime is greater than k(ln
kln ln k-1) for kgt 2," Math. Comp., 68225
(January 1999) 411--415.  MR 99d11133 (Abstract
available)
60
(No Transcript)
61
(No Transcript)
62
(No Transcript)
63
(No Transcript)
64
(No Transcript)
65
(No Transcript)
66
(No Transcript)
67
(No Transcript)
68
La función zeta ?(z)
Euler la llamó función zeta en 1737. Consideró
que z era un real mayor que 1.
69
Repitamos la operación para el siguiente primo 3.
70
Producto de Euler para la función zeta.
71
La función zeta ?(s) de Riemann
Bernhard Riemann hacia 1859 generalizó la función
zeta a números s x iy complejos. Aquí vemos
una representación gráfica del módulo de la
función z de Riemman ? (s). Obsérvese el polo
en s 1.
72
Aquí vemos una representación gráfica del módulo
de la inversa de la función z de Riemman 1/?
(s). De este modo podemos ver fácilmente los
ceros de la función z como polos. Los ceros
parece que vayan paralelos y cercanos al eje
imaginario.
73
Hipótesis de Riemann
Los pares negativos (-2, -4, -6, etc.) son ceros
de la función zeta, los llamados ceros triviales.
La hipótesis de Riemman afirma que todos los
demás ceros, llamados no triviales, tienen parte
real igual a ½. Es decir, que son de la forma ½
iy.
Grafica de y frente al módulo
74
Los 10.000 primeros millones de ceros de la
función zeta están en la línea crítica ½ (Mayo
2002).
75
(No Transcript)
76
(No Transcript)
77
Normalicemos P(d)
78
(No Transcript)
79
Estamos trabajando en el intervalo 1, 10D1),
es decir en D 1 décadas. Normalicemos ahora
la función densidad de probabilidad de ceros
De modo que la función densidad de ceros es
Como antes esta densidad nos da directamente una
ley de Benford generalizada normalizada.
80
(No Transcript)
81
Si pasamos de décadas y hablamos de N, estamos
diciendo que
Tomando como dependencia para alfa
82
(No Transcript)