Lição
A teoria — Crivo de primos e espiral de Ulam
O painel imprime dois tipos distintos de número, e vale a pena separá-los antes de ler mais alguma coisa. π(N) é uma contagem — quantos primos o crivo deixou realmente de pé, 25 em N = 100 e 95 em N = 500, exata e sem arredondamento. As duas linhas por baixo são estimativas dessa mesma contagem. O teorema dos números primos afirma que uma delas acerta o quociente no limite; não afirma que acerta o número. São promessas muito diferentes, e é na coluna do erro que a diferença aparece.
O que significa cada símbolo
N- o teto até ao qual se criva, e a única entrada. O cursor vai até 10000.
π(N)- a quantidade de primos que não excedem N. Contada, não estimada — são simplesmente as casas que o crivo não riscou.
N / ln N- a estimativa mais simples dessa contagem. Em todo o N que este cursor alcança fica abaixo da verdade.
Li(N)- o logaritmo integral
∫₂ᴺ dt/ln t, uma estimativa mais fina. Em todo o N que este cursor alcança fica acima da verdade.
De onde vem a fórmula
- O crivo nunca pergunta se um número é primo. Começa no 2 e risca todos os múltiplos de 2, passa ao número seguinte ainda de pé, risca todos os seus múltiplos e repete. A primalidade não é testada; é o que sobra.
- Basta riscar múltiplos de primos até
√N. Se um número n ≤ N é composto, então n = a·b com a ≤ b, logo a ≤ √n ≤ √N — todo o composto tem um divisor no máximo igual à raiz, e por isso já foi riscado quando chegou a vez desse divisor. Crivar até 100 exige apenas as passagens do 2, 3, 5 e 7. - Conta os sobreviventes e tens
π(N)exatamente. É por isso que a linha de cima é um facto e as duas de baixo são opiniões: o crivo produz a contagem como subproduto de terminar. - Agora compara. O teorema dos números primos diz que
π(N) · ln N / N → 1quando N cresce. Repara no que ele não diz: um quociente que tende para 1 permite que o erro percentual continue grande durante muitíssimo tempo — e a linha seguinte mostra exatamente isso.
Como ler o que vê
Lê os dois valores de erro como uma corrida e muda o N. Em N = 100 a estimativa grosseira vai à frente: 13.1% contra os 16.3% de Li. Vai até N = 200 e inverte-se — 17.9% contra 6.9% — e nunca mais se inverte. Mas o que interessa mesmo observar é aquilo que N/ln N não faz. Ao longo de duas ordens de grandeza o seu erro dá 14.8%, 13.1%, 15.3%, 13.8%, 12.2%: oscila e quase não melhora. Li, no mesmo intervalo, passa de 16.3% para 2.8%. Ambas as estimativas cumprem o teorema dos números primos; só uma serve nas escalas que consegues ver. Os sinais são igualmente constantes — aqui N/ln N fica abaixo da contagem em todo o N, e Li fica acima.
- Pressupõe
- Que N é suficientemente pequeno para ser crivado por inteiro: a contagem é exata porque cada número até N é mesmo guardado e riscado, e é por isso que o teto é 10000 e não 10¹⁰. As percentagens de erro são tomadas em relação a π(N), pelo que medem as estimativas e nunca a contagem.
- Falha quando
- Li(N) > π(N) em todo o N que esta página consegue alcançar, o que lhe dá ar de lei. Não é. Littlewood provou em 1914 que a diferença muda de sinal infinitas vezes, portanto existem valores de N em que Li fica aquém — e no século que passou ninguém exibiu um único. Bays e Hudson colocaram o primeiro cruzamento abaixo de cerca de 1.4×10³¹⁶ em 1999, e desde então ninguém o estreitou. Esta página mostra-te, então, um padrão que está provado não ser universal, e nenhuma posição do cursor alcança o contraexemplo. Essa distância entre o que se vê e o que é verdade é a forma honesta deste assunto.
Problemas resolvidos na íntegra
-
Crivar os 25 primos abaixo de 100 à mão 5 passos
Existem 25 números primos abaixo de 100. Passe-os pelo crivo à mão — dá menos trabalho do que imaginaria — e depois teste o teorema dos números primos numa amostra demasiado pequena para ele.
-
A poupança do crivo é a razão pela qual ele merece um nome. Se n ≤ 100 for composto, fatora-se como ab, e o fator menor não pode exceder √100 = 10. Portanto, eliminar os múltiplos de todos os primos até 10 remove todos os compostos existentes: quatro primos fazem todo o trabalho.
-
Risque os múltiplos de 2, 3, 5 e 7, cada um começando no seu próprio quadrado, porque tudo o que está abaixo disso já foi eliminado. Sobrevivem vinte e cinco números.
-
O teorema dos números primos diz que a contagem é assintoticamente N/ln N. Para N = 100 isso dá 21,7.
-
Está 13% abaixo, o que é aceitável para um resultado assintótico em 100. Lido ao contrário, continua a ser útil aqui: a densidade de primos perto de N é 1/ln N, pelo que cerca de 22% dos números perto de 100 são primos, contra 7,2% perto de um milhão. Os primos escasseiam de forma logarítmica, o que é de facto muito lentamente.
-
Os intervalos concordam. O intervalo médio abaixo de 100 é 100/25 = 4, e o maior é 8 — a sequência de 89 a 97. O dobro da média, e nada pior.
Resposta
O crivo dá π(100) = 25 e a ferramenta apresenta a estimativa em 21,7, 13,1% abaixo. Esse erro não é uma falha do teorema; é a taxa de convergência do teorema tornada visível, e a convergência é notoriamente pausada — aumente N para 10 000, quatro ordens de grandeza acima, e a estimativa continua abaixo por dois dígitos. O que se mantém em todas as escalas é a leitura da densidade. Um número em cada 4,6 perto de 100 é primo, um em cada 13,8 perto de um milhão, e como o logaritmo cresce tão devagar, os primos nunca se esgotam e nem sequer se tornam verdadeiramente raros.
-
-
A estimativa óbvia para 8 pares gémeos abaixo de 100 6 passos
O painel conta 25 primos abaixo de 100 e 8 pares de gémeos. O primeiro número tem uma estimativa célebre que chega aos 13 %. Experimente a estimativa óbvia para o segundo — e veja-a falhar por um fator que tem nome.
-
Tome a contagem e a estimativa impressas como ponto de partida: perto de N um número é primo com probabilidade cerca de 1/ln N, e em N = 100 isso dá uns 0,2171.
-
Suponha agora que n e n+2 são independentes. Se cada um é primo com probabilidade 1/ln N, ambos o são com o quadrado disso.
-
Multiplique por N e compare com a ferramenta. A estimativa dá 4,72 pares gémeos abaixo de 100; o crivo encontrou 8. Não fica um pouco aquém: fica 40 % aquém.
-
A independência é o passo falso, e basta um primo para o mostrar. Tome um primo ímpar p: um n ao acaso é eliminado por p uma vez em cada p, mas o par é eliminado sempre que n ou n+2 for divisível por p, o que são dois restos em p. A taxa de sobrevivência é (p−2)/p, e não o quadrado de (p−1)/p que a independência supunha.
-
Multiplique essa correção sobre todos os primos ímpares e ela converge para uma constante. O seu dobro é 1,3203, e aplicá-la eleva a estimativa a 6,23.
-
Continua abaixo de 8 em N = 100 — a estimativa é assintótica e 100 é um número pequeno. Converge: 156 pares previstos abaixo de 10⁴, 6917 abaixo de 10⁶.
Resposta
A contagem ingénua é 4,72, a corrigida 6,23, e a verdade é 8 — e o fator de correção 1,3203 é um produto infinito sobre os primos. Essa constante é o preço de supor independência onde não a há, e tem a forma de quase todas as perguntas difíceis da área: os primos são suficientemente aleatórios para que uma heurística funcione, e suficientemente estruturados para que ela precise de uma correção que ninguém sabe deduzir do nada. A conjetura dos primos gémeos afirma que a esta estimativa nunca faltam pares, e continua em aberto.
-
Referências (4)
- Why the page’s most convincing pattern is not a law — the difference changes sign infinitely often: J. E. Littlewood, "Sur la distribution des nombres premiers." Comptes Rendus de l’Académie des Sciences, Paris 158, 1869–1872, 1914. No DOI: the volume predates them.
- How far the first sign change has been pinned down, and where the 1.4×10³¹⁶ comes from: C. Bays & R. H. Hudson, "A new bound for the smallest x with π(x) > li(x)." Mathematics of Computation 69, 1285–1296, 1999.
- How far the Riemann hypothesis has actually been verified — note this is a bound on height, a different measure from a count of zeros: D. J. Platt & T. S. Trudgian, "The Riemann hypothesis is true up to 3·10^12." Bulletin of the London Mathematical Society 53, 792–797, 2021.
- The prize, and its official problem statement: Clay Mathematics Institute, Millennium Prize Problems — the Riemann Hypothesis.