Úvod
V tomto článku budú popísané ďalšie súvislosti na predchádzajúci diel o empirickom preverení Riemannovej hypotézy. Krátke zhrnutie z predchádzajúceho dielu je, že ekvivaletné podmienky pre Riemannovú hypotézu zrejme platia, pričom v predchádzajúcom diely boli preverené ideálne čísla (teda tie, čo maximalizujú Grönwall vzťah) na bežnom pc do veľkosti zhruba 10 exp 6722000. Zároveň bola vypočítaná aj aproximácia z vlastnosti priebehu ideálnych čísiel.Referencia:
(1) RAMANUJAN, ROBIN, HIGHLY COMPOSITE NUMBERS, AND THE RIEMANN HYPOTHESIS
(2) http://math.colgate.edu/~integers/l33/l33.pdf
(3) Prime_number_theorem
![Screenshot - 11_ 1 003jpg](https://files.vlastnawebstranka.websupport.sk/69/24/69249aa2-81f3-4e8d-b054-f214d3c456dc.jpg)
Poďme sa teda pozrieť na ďalšie súvislosti empirického výpočtu z predchádzajúceho dielu. Výpočet bol vykonaný pre dve nosné sekvencie, teda: pre sekvenciu (1) a sekvenciu (3) - ideálne čísla. Odporúčam si pozrieť predchádzajúci článok, kde sú vysvetlené jednotlivé pojmy ako sekvencia (1), (3), ideálne čísla a podobne.
Sekvencia (1) - postupnosť čísiel 2,6,30,210..
![Screenshot - 11_ 1jpg](https://files.vlastnawebstranka.websupport.sk/87/c3/87c32594-188e-496e-b26f-eaf3ac5d6f9f.jpg)
obr.1 Tabuľka z pokročilého testovania sekvencie (1).
Sekvencia (3) - postupnosť ideálnych čísiel (highly composite numbers)
![Screenshot - 11_ 1 002jpg](https://files.vlastnawebstranka.websupport.sk/c4/07/c407d1b2-5a2a-44ac-81f4-ae80ab789ccc.jpg)
Obr.2 Tabuľa z pokročilého testovania sekvencie (3).
V prvom rade si všimneme veľkú podobnosť priebehu sekvencie (1) a (3). Tieto priebeh sú zrkadlovo otočené. Nie sú však dokonale symetrické. U sekvencie vzniká taktiež limitné úroveň na hodnote 1.082..Pri sekvencii (3) je to hodnota e ˆ gama.
Ramunujan nazval sekvneciu (3) vysoko- zložené čísla. Teda ideálne čísla, viď. predošlé diely - program sú highly composite number.
Vráťme sa však na chvíľu k prvo-číselnej vete:
![999613e5-b676-46b0-8198-968175f932a6jpg](https://files.vlastnawebstranka.websupport.sk/0f/ab/0fab0e7b-145e-4e89-b1bb-405bf41768ee.jpg)
Obr. 3 počet prvočísiel v N
![ea5c0b15-3c89-48e8-ac80-42356a535cfejpg](https://files.vlastnawebstranka.websupport.sk/4f/be/4fbe7f69-fbcb-4e61-b591-315fd3a9065f.jpg)
Z tabuliek však vidieť, že platí aj ďalšia vlastnosť (Autor: Robopol):
![Screenshot - 11_ 1 004jpg](https://files.vlastnawebstranka.websupport.sk/61/cc/61cc939d-6c2f-4ffd-89df-44a8b3689353.jpg)
![Screenshot - 11_ 1 005jpg](https://files.vlastnawebstranka.websupport.sk/79/cd/79cdc601-b1d3-48c6-bb05-fccd4cdd212d.jpg)
Ešte presnejšie vzťahy dostaneme tieto:
![Screenshot - 11_ 1 006jpg](https://files.vlastnawebstranka.websupport.sk/71/2f/712ff785-1035-44cf-b8f4-740a8bfeccab.jpg)
Mertensova veta:
Veľmi pekný vzťah (úzky súvis) medzi zeta funkciou (Riemannova hypotéza) a Mertensova veta, vzťah ukazuje napr. tu - odkaz.Lagarias theorem:
![Screenshot - 25_ 1jpg](https://files.vlastnawebstranka.websupport.sk/a0/cf/a0cf11cb-2b65-4fb1-a6b1-3c284e7259a4.jpg)
zdroj: referencia (2)
Odvodenie cez Lagarias theorem
Vo vzťahu vystupuje tzv. Hn - čo je harmonic series (harmonické čísla), ktoré sú dané vzťahom:![Screenshot - 25_ 1 002jpg](https://files.vlastnawebstranka.websupport.sk/5b/e4/5be4709c-c2e6-4abd-a4e0-20af0286e6a5.jpg)
Viac o Harmonický číslach nájdeme napr. tu: odkaz
Nakoľko táto rada pomaličky diverguje, teda výpočet pre veľké čísla by zabral veľa času a výpočtového výkonu pre harmonické čísla platí:
![Screenshot - 25_ 1 003jpg](https://files.vlastnawebstranka.websupport.sk/5f/f9/5ff9672c-fe51-4f64-94f9-8865eac49bcb.jpg)
Z toho potom dostaneme:
![Screenshot - 4_ 2jpg](https://files.vlastnawebstranka.websupport.sk/9d/34/9d34a877-ec3b-49e0-9a5d-59c60b93c71b.jpg)
Z tohto výpočtu plynie, že ak vysoko - zložené čísla (sekvencia (3)) majú sigma (N) menšiu ako udáva Lagariasov vzťah, potom by nebola nikdy prekročená hodnota e^gama. Kľúčové je teda preukázať, že vysoko- zložené čísla (ich sigma) neprekročí Lagariasov vzťah. To by bol dôkaz ekvivaletných podmienok pre Riemannovú hypotézu.
Môžeme teda skontrolovať ako to bude, pre veľmi veľké čísla z tabuľky, viď. vyššie. No, ale ak sa pozrieme na upravený Lagarias vzťah (ako sa hovorí : lenivosť je matka pokroku) zistíme, že sme dospeli ku Grönwall vzťahu, kde je len malá zmena:
Lagarias vzťah je vždy väčší o faktor: +ln(N)+gama a ešte je člen N ln(ln N) zväčšený o koeficienty gama a násobený e^ε. Pri veľmi veľkých číslach však bude ten rozdiel (Lagariass- Grönwall) sústredený v ln(N), pre výpočet sigma(n). Teda Lagarias vzťah nie je nutné preverovať, keďže hodnoty sú väčšie ako u Grönwall vzťahu, kde algoritmus preveril obrovské čísla, ktoré potvrdili ekvivalentné podmienky pre Riemannovu hypotézu do zhruba 10 exp 6722000.
Ideálne čísla - highly composite numbers
Vráťme sa k schéme ideálnych čísiel, pretože pre dôkaz je potrebné formulovať jednoznačný vzťah pre tieto vysoko - zložené čísla. V predchádzajúcom diely je vzťah pre výpočet sigma ideálnych čísiel, zároveň táto formula je uvedená aj v referencii (2):teda platí:
![highly_composite_numbersjpg](https://files.vlastnawebstranka.websupport.sk/30/88/3088a72a-fa91-4030-99d3-92376c874ad3.jpg)
Rýchly výpočet vo wolframe nám ukáže napr. pre k=100:
sigma (k=100)/N= odkaz
11.2676 je približne rovné ln(ln(N))*e^gama=11.323929 : N je vyrátané programom v pythone, viď. predošlý diel.
sigma(k=500)/N= odkaz
14.600 je približne rovné ln(ln(N))*e^gama=14.598277...
Mali by však platiť aj tieto vzťahy:
![Screenshot - 21_ 3jpg](https://files.vlastnawebstranka.websupport.sk/38/1f/381ffd39-43f5-4a24-b017-faca6fc1ff70.jpg)
Preverme tieto nové vzťahy malým programom v pythone:
![Screenshot - 31_ 1jpg](https://files.vlastnawebstranka.websupport.sk/23/21/23214f0b-d039-4ca3-b491-744da0d4b243.jpg)
obr.4 zdrojový kód programu
Verzia programu na test Reformulation of Riemann hypothesis:
Download súboru z GitHub: sigma-max_test.py
Výsledky testov pre veľmi veľké čísla ukazuje tabuľka:
![Screenshot - 4_ 2 002jpg](https://files.vlastnawebstranka.websupport.sk/b0/20/b02032f8-f602-4294-9fe0-0a37c1821dea.jpg)
Z testu teda nerovnosť platí a pre veľké čísla vzťah doslova kopíruje priebeh sigma(max). Rozdiel sa stále postupne zväčšuje, čo znamená, že v zmysle ekvivalentných podmienok pre platnosť Riemannovej hypotézy naopak má logaritmus o niečo málo navrch oproti sigma(N)/N. V zmysle Guy Robin empirický výpočet vývoja by e ^ ln(ln(N)) bol o máličko silnejší faktor ako Sigma(N)/N. To je ďalšia indícia, že ekvivalentné podmienky pre Riemannovú hypotézu zrejme platia. Ostáva teda už len analytický dôkaz, že sigma(N)/N je smerom do nekonečna vždy menšia ako e ^ gama* ln(ln(N)).
Hľadanie analytického dôkazu
Poďme skúsiť, čo nám derivácia týchto funkcii povie (reformulácia vzťahov Robopol). Pri sigma/N však ide o diskrétnu postupnosť. Z tohto dôvodu (derivácia nedáva zmysel) sa teda pozrieme na susedné body (susedné prvočísla) viď. nižšie. Zoberieme si počiatočný bod p(k) a pozrieme sa ako rastie hodnota pre nasledujúce p(k+1). Na obrázku č.6 je odvodenie.![Screenshot - 8_ 2jpg](https://files.vlastnawebstranka.websupport.sk/b1/4d/b14debc2-92b2-4522-896f-6bad0e58fa6a.jpg)
obr.6 Odvodenie rovnice pre susedné prvočísla
Dostali sme rovnicu, ktorá hovorí ako rastie funkcia pre ľubovolné p(k) - prvočíslo a jeho suseda. Keďže prvočísla vykazujú malé fluktuácie okolo vyhladenej krivky napr. v zmysle vzťahu pi(x)~Li(x), viď. prvo-číselná veta r.1849 nemôžeme očakávať, že ľavá strana rovnice bude pre všetky p(k) > ako pravá strana rovnice. To by však mohlo platiť pre vyhladenú krivku, ktorá aproximuje pi(x). Poďme to preskúmať.
V zmysle referencie (3) - Prime number theorem je vzťah:
![Screenshot - 8_ 2 002jpg](https://files.vlastnawebstranka.websupport.sk/4a/52/4a524f92-dbca-4b25-8936-2c4dcf083506.jpg)
Obr. 7 Max prime aproximation.
Vytvorime rovnicu, ktorá nám hovorí aká je x- ová vzdialenosť medzi dvoma prvočíslami (z vyhladenej krivky). No urobíme zjednodušenú variantu z rovnice x/log x, pretože pre uvedený vzťah (obr.7) dáva veľmi dlhé rovnice. Z tohto dôvodu urobíme malé zjednodušenie, kde sa princíp nemení.
![Screenshot - 24_ 2jpg](https://files.vlastnawebstranka.websupport.sk/ea/d1/ead1a3a3-b442-43ac-b94a-b62f42a38284.jpg)
Obr.8 Odvodené rovnice pre x/log(x). Poznámka: log(x) je prirodzený logaritmus.
Výsledky ručných testov vo Wolframe:
- Pre x/log(x) vychádza podmienka log(x+delta)>log(x)*(x+delta)/(x+delta-1) pri ručnom náhodnom testovaní splnená. Nenašiel som porušenie do čísiel, ktoré Vám Wolfram zadarmo vypočíta.
- Pre horné ohraničenie x/(log(x)-1.1) podmienka: log(x+delta)>log(x)*(x+delta)/(x+delta-1) pri ručnom testovaní zhruba od 10 ^ 6 nie je splnená, pravá strana rovnice je väčšia (v rozdiele na 6 desatinnom mieste a vyššie, záleží od veľkosti testovaného čísla).
- Pre aproximáciu pi(x): x/(log(x)-1), podmienka: log(x+delta)>log(x)*(x+delta)/(x+delta-1) pri ručnom náhodnom testovaní splnená. Nenašiel som porušenie do čísiel, ktoré Vám Wolfram zadarmo vypočíta.
V zmysle referencie (3) postačuje splniť nasledujúcu podmienku horného ohraničenia pi(x):
![Screenshot - 16_ 2 003jpg](https://files.vlastnawebstranka.websupport.sk/f5/30/f530924e-e26c-4a9d-9e00-6f4e49ce7d7a.jpg)
Teda nie je nutné preukázať to pre x/(log(x)-1.1). Ručné testovanie (vo Wolframe) ukázalo, že pre toto ohraničenie nie je splnená podmienka (viď. vyššie).
zvolme si epsilon=1.00005, potom otestujeme zopár čísiel a zistíme pre rovnicu pi(x)=x/(log(x)-1.00005):
napr. pre x=10^11 bude delta_x=25.3712
podmienka je splnená min delta_x=25.3284
Zopár som ich vyskúšal a podmienka je splnená do 10^12. Zdá sa, že aproximáciu pi(x): x/(log(x)-1.00005) spĺňa. Ručné testovanie rovníc sedí pre podmienku horného ohraničenia pi(x). Tým je teda empiricky do určitého čísla, ktoré Wolfram vypočíta podmienka splnená.
Empirické riešenie rovnice cez Python:
![Screenshot - 22_ 2jpg](https://files.vlastnawebstranka.websupport.sk/57/52/57529047-9718-4c32-a7e6-80ca635e129f.jpg)
Obr.9 Zdrojový kód numerického riešenia rovnice
Verzia programu na test numerického riešenia rovnice:
Download súboru z GitHub: numerical_solution_1.py
Výsledky testov od 10 ^2 do 10 ^ 14 ukazuje tabuľka:
![Screenshot - 22_ 2 002jpg](https://files.vlastnawebstranka.websupport.sk/e2/bf/e2bf7c9d-a3dd-4350-8167-57406b8ca547.jpg)
Obr.10 priebeh testovania: od 10^2 do 10 ^14
Zhodnotenie: Výsledky sú veľmi podobné ako vo Wolframe, výsledky rovnice však numerické metódy (bez úpravy premenných atď.) spoľahlivo ukazujú max do 10 ^12. Z priebehu však vidíme, že pre epsilon=1.00005 rozdiel (delta-delta_min) klesá smerom k nule. Z testovania nie je však možné povedať, či sa (delta-delta_min) dostane do záporných čísiel smerom k nekonečnu.
----pokračovanie v nasledujúcom článku---