Lektion
Die Theorie — Primzahlsieb & Ulam-Spirale
Das Feld druckt zwei verschiedene Arten von Zahl, und es lohnt sich, sie zu trennen, bevor man irgendetwas anderes liest. π(N) ist eine Anzahl — wie viele Primzahlen das Sieb tatsächlich stehen ließ, 25 bei N = 100 und 95 bei N = 500, exakt und nicht gerundet. Die beiden Zeilen darunter sind Schätzungen ebendieser Anzahl. Der Primzahlsatz behauptet, dass eine von ihnen im Grenzwert das Verhältnis richtig trifft; er behauptet nicht, dass sie die Zahl richtig trifft. Das sind sehr verschiedene Versprechen, und in der Fehlerspalte wird der Unterschied sichtbar.
Was die einzelnen Symbole bedeuten
N- die Obergrenze, bis zu der gesiebt wird, und die einzige Eingabe. Der Regler reicht bis 10000.
π(N)- die Anzahl der Primzahlen, die N nicht überschreiten. Gezählt, nicht geschätzt — es sind schlicht so viele Felder, wie das Sieb nicht durchgestrichen hat.
N / ln N- die einfachste Schätzung dieser Anzahl. Bei jedem N, das dieser Regler erreicht, liegt sie unter der Wahrheit.
Li(N)- das Integrallogarithmus
∫₂ᴺ dt/ln t, eine schärfere Schätzung. Bei jedem N, das dieser Regler erreicht, liegt sie über der Wahrheit.
Woher die Formel kommt
- Das Sieb fragt nie, ob eine Zahl prim ist. Es beginnt bei 2, streicht jedes Vielfache von 2, geht zur nächsten noch stehenden Zahl, streicht jedes Vielfache davon und wiederholt das. Primalität wird nicht geprüft; sie ist das, was übrig bleibt.
- Man muss nur Vielfache von Primzahlen bis
√Nstreichen. Ist eine Zahl n ≤ N zusammengesetzt, so ist n = a·b mit a ≤ b, also a ≤ √n ≤ √N — jede zusammengesetzte Zahl hat einen Teiler höchstens in Höhe der Wurzel und wurde daher schon gestrichen, als dieser Teiler an der Reihe war. Bis 100 zu sieben braucht nur die Durchgänge für 2, 3, 5 und 7. - Zähle die Überlebenden, und du hast
π(N)exakt. Deshalb ist die oberste Zeile eine Tatsache und die beiden darunter sind Meinungen: das Sieb liefert die Anzahl als Nebenprodukt seines Durchlaufs. - Nun vergleiche. Der Primzahlsatz sagt
π(N) · ln N / N → 1für wachsendes N. Beachte, was er nicht sagt: ein gegen 1 strebendes Verhältnis erlaubt es dem prozentualen Fehler, sehr lange groß zu bleiben — und die nächste Zeile zeigt genau das.
So liest du, was du siehst
Lies die beiden Fehlerwerte als Wettrennen und verändere N. Bei N = 100 führt die grobe Schätzung: 13.1% gegen Lis 16.3%. Gehe auf N = 200, und es kehrt sich um — 17.9% gegen 6.9% — und es kehrt sich nie wieder zurück. Das eigentlich Bemerkenswerte ist aber, was N/ln N nicht tut. Über zwei Größenordnungen lautet sein Fehler 14.8%, 13.1%, 15.3%, 13.8%, 12.2%: er schwankt und verbessert sich kaum. Li geht über dieselbe Spanne von 16.3% auf 2.8%. Beide Schätzungen erfüllen den Primzahlsatz; nur eine davon taugt in den Größenordnungen, die du sehen kannst. Ebenso beständig sind die Vorzeichen — N/ln N liegt hier bei jedem N unter der Anzahl, Li über ihr.
- Setzt voraus
- Dass N klein genug ist, um vollständig gesiebt zu werden: die Anzahl ist exakt, weil jede Zahl bis N tatsächlich vorgehalten und gestrichen wird — deshalb endet die Obergrenze bei 10000 und nicht bei 10¹⁰. Die Fehlerprozente beziehen sich auf π(N), messen also die Schätzungen und nie die Anzahl.
- Versagt, wenn
- Li(N) > π(N) gilt bei jedem N, das diese Seite erreichen kann, und sieht damit aus wie ein Gesetz. Es ist keines. Littlewood bewies 1914, dass die Differenz unendlich oft das Vorzeichen wechselt, es also N gibt, bei denen Li zu wenig zählt — und im Jahrhundert seither hat niemand auch nur ein einziges vorgezeigt. Bays und Hudson drückten den ersten Wechsel 1999 unter etwa 1.4×10³¹⁶, und enger ist er seither nicht geworden. Diese Seite zeigt dir also ein Muster, das nachweislich nicht allgemein gilt, und keine Reglerstellung erreicht das Gegenbeispiel. Diese Lücke zwischen dem Sichtbaren und dem Wahren ist die ehrliche Gestalt dieses Themas.
Aufgaben vollständig gelöst
-
Sieben der 25 Primzahlen unter 100 von Hand 5 Schritte
Es gibt 25 Primzahlen unter 100. Sieben Sie sie von Hand aus — das macht weniger Arbeit, als man vermuten würde — und testen Sie dann den Primzahlsatz an einer Stichprobe, die weit zu klein für ihn ist.
-
Die Ersparnis des Siebs ist der Grund, warum es einen eigenen Namen verdient. Wenn n ≤ 100 zusammengesetzt ist, zerfällt es als ab, und der kleinere Faktor kann √100 = 10 nicht überschreiten. Das Streichen der Vielfachen jeder Primzahl bis 10 entfernt somit jede zusammengesetzte Zahl: Vier Primzahlen erledigen die gesamte Aufgabe.
-
Streichen Sie die Vielfachen von 2, 3, 5 und 7, jeweils beginnend bei deren Quadrat, da alles darunter bereits gestrichen ist. Fünfundzwanzig Zahlen bleiben übrig.
-
Der Primzahlsatz besagt, dass die Anzahl asymptotisch N/ln N beträgt. Bei N = 100 ist das 21,7.
-
Das ist um 13 % zu niedrig, was ein asymptotisches Ergebnis bei 100 durchaus sein darf. Umgekehrt gelesen ist es hier dennoch nützlich: Die Dichte der Primzahlen nahe N ist 1/ln N, sodass etwa 22 % der Zahlen nahe 100 Primzahlen sind, gegenüber 7,2 % nahe einer Million. Primzahlen dünnen sich logarithmisch aus, was in der Tat sehr langsam geschieht.
-
Die Lücken stimmen überein. Die mittlere Lücke unter 100 beträgt 100/25 = 4, und die größte ist 8 — die Folge von 89 bis 97. Das Doppelte des Durchschnitts und nicht schlechter.
Antwort
Das Sieb liefert π(100) = 25 und das Werkzeug gibt den Schätzwert mit 21,7 an, 13,1 % zu niedrig. Dieser Fehler ist kein Versagen des Satzes; er macht die Konvergenzgeschwindigkeit des Satzes sichtbar, und diese Konvergenz ist bekanntermaßen uneilig — erhöht man N auf 10 000, um vier Größenordnungen nach oben, liegt die Schätzung immer noch im zweistelligen Bereich zu niedrig. Was auf jeder Skala gilt, ist die Aussage über die Dichte. Eine Zahl von 4,6 nahe 100 ist prim, eine von 13,8 nahe einer Million, und weil der Logarithmus so langsam wächst, gehen die Primzahlen nie aus und werden nicht einmal wirklich selten.
-
-
Die offensichtliche Schätzung für 8 Zwillingspaare unter 100 6 Schritte
Die Tafel zählt 25 Primzahlen unter 100 und 8 Primzahlzwillinge. Für die erste Zahl gibt es eine berühmte Schätzung, die auf 13 % herankommt. Versuche die naheliegende Schätzung für die zweite — und sieh zu, wie sie um einen Faktor danebenliegt, der einen Namen hat.
-
Nimm die gedruckte Anzahl und die gedruckte Schätzung als Ausgangspunkt: nahe N ist eine Zahl mit etwa 1/ln N Wahrscheinlichkeit prim, bei N = 100 also rund 0,2171.
-
Nun nimm an, n und n+2 seien unabhängig. Ist jede mit Wahrscheinlichkeit 1/ln N prim, so sind es beide mit dem Quadrat davon.
-
Multipliziere mit N und vergleiche mit dem Werkzeug. Die Schätzung sagt 4,72 Zwillingspaare unter 100; das Sieb fand 8. Das ist nicht knapp daneben — das sind 40 % zu wenig.
-
Die Unabhängigkeit ist der falsche Schritt, und eine einzige Primzahl zeigt warum. Nimm eine ungerade Primzahl p: ein zufälliges n wird von p in einem von p Fällen erschlagen, das Paar aber immer dann, wenn n oder n+2 durch p teilbar ist — zwei Reste von p. Die Überlebensrate ist also (p−2)/p und nicht das Quadrat von (p−1)/p, das die Unabhängigkeit unterstellt.
-
Multipliziere diese Korrektur über alle ungeraden Primzahlen, und sie konvergiert gegen eine Konstante. Verdoppelt ist sie 1,3203, und angewandt hebt sie die Schätzung auf 6,23.
-
Bei N = 100 immer noch zu wenig gegenüber 8 — die Schätzung ist asymptotisch, und 100 ist eine kleine Zahl. Sie konvergiert: 156 Paare unter 10⁴, 6917 unter 10⁶.
Antwort
Die naive Zahl ist 4,72, die korrigierte 6,23, und die Wahrheit ist 8 — und der Korrekturfaktor 1,3203 ist ein unendliches Produkt über die Primzahlen. Diese Konstante ist der Preis dafür, Unabhängigkeit anzunehmen, wo keine herrscht, und sie hat die Gestalt fast jeder schweren Frage des Fachs: die Primzahlen sind zufällig genug, dass eine Heuristik funktioniert, und strukturiert genug, dass sie eine Korrektur braucht, die niemand von Grund auf herleiten kann. Die Primzahlzwillingsvermutung besagt, dass dieser Schätzung nie die Paare ausgehen, und sie ist bis heute offen.
-
Quellen (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.