奇妙的π已經(jīng)精算到100萬億位，后面的秘密還有多少？

圓周率（簡稱π）是數(shù)學(xué)中最神秘又熟悉的數(shù)字之一。無論是在建筑、工程，還是航天、物理，π都發(fā)揮著至關(guān)重要的作用。簡單的說，圓周率就是一個(gè)圓的周長與其直徑之比?？雌饋硪玫竭@個(gè)數(shù)字很簡單，只要量一下一個(gè)圓的周長，用這個(gè)周長除以直徑，圓周率就得到了。

但困難的是，這個(gè)圓用什么尺子才能夠?qū)⒅荛L量的十分精確呢？從古到今，幾千年過去了，這個(gè)傷腦筋的難題也沒有完全得到解決。而得到的結(jié)果就是，圓周率是一個(gè)無理數(shù)，意味著它的小數(shù)部分無限多，永無止境，且沒有任何循環(huán)規(guī)律。

那么這個(gè)小數(shù)位數(shù)到底有多少呢？沒人知道，因?yàn)槠袢藗冎恢朗菬o限的，無限就是永無止境的。幾千年前，古人們就試圖用手工的方法來測量出精確的圓周長，他們用尺規(guī)和幾何方法手工計(jì)算π，得到了一個(gè)大致的數(shù)據(jù)，最終精確到了小數(shù)點(diǎn)后九位數(shù)，即3.141615926。這是咱們老祖宗祖沖之得到的，可以自豪一下。

今天，隨著科學(xué)的進(jìn)步，超級計(jì)算機(jī)將π推算到100萬億位，但也只能說明它還是個(gè)近似值，依然沒有窮盡。100萬億位是什么概念？就是一個(gè)人如果從出生就開始算數(shù)，不吃不喝24小時(shí)不停的算，平均每秒鐘算兩個(gè)數(shù)的話，要算158萬多年。也就是說，一個(gè)人從原始時(shí)代猿人老祖宗開始，世世代代接著算，算到現(xiàn)在也還沒算完。

那么，圓周率真的有這么多位嗎，人類開始時(shí)如何得到的呢？我們一起來了解。

古人的智慧：用割圓術(shù)計(jì)算π

無論是西方還是東方，早在遠(yuǎn)古時(shí)代就有那么一群智者對圓周率有了興趣，并開始研究。東西合璧殊途同歸，得到了圓周率的近似值。具體說來，有如下一些典型代表：

阿基米德的“多邊形逼近法”

最早系統(tǒng)計(jì)算π的數(shù)學(xué)家之一是古希臘的阿基米德（約公元前287—前212年）。他的想法很簡單：

1. 畫一個(gè)圓，然后在圓里面畫一個(gè)正六邊形，再在外面畫一個(gè)正六邊形。
2. 計(jì)算這兩個(gè)六邊形的周長，就可以得到一個(gè)π的上下限。
3. 然后，把六邊形的邊數(shù)加倍，變成十二邊形，再算一次周長。
4. 繼續(xù)加倍，變成二十四邊形、四十八邊形……，直到逼近圓的真實(shí)周長。

阿基米德用這種方法，把π的值估算在3.1408 到 3.1429 之間，這在沒有計(jì)算工具的時(shí)代已經(jīng)是非常驚人的成就！

祖沖之：計(jì)算到全球領(lǐng)先1000年的精度

祖沖之是南北朝時(shí)期（公元5世紀(jì)）中國古代數(shù)學(xué)家，他改進(jìn)了割圓術(shù)，將一個(gè)正元分割出24576個(gè)邊，通過對這24576邊形的測量計(jì)算，得到了比阿基米德的精度更高的π近似值：3.1415926到3.1415927，這個(gè)精度在世界范圍內(nèi)領(lǐng)先了1000年！

但問題是，割圓術(shù)計(jì)算π非常慢。如果想要更精確的π，就得畫邊數(shù)更多的多邊形，計(jì)算量成倍增長，最終會變得難以承受。所以，祖沖之之后的1000年間，π的精度再也難以提升。到了近代，數(shù)學(xué)家們開始尋找更高效的計(jì)算方法，讓π的計(jì)算速度大幅提升。

數(shù)學(xué)公式的力量：不用畫圖，也能精確計(jì)算π，且比割圓術(shù)快很多

古人計(jì)算π是用形狀逼近圓，而現(xiàn)代數(shù)學(xué)家則用公式直接計(jì)算π，這種方法比割圓術(shù)快得多，且更精準(zhǔn)。比較有代表性的公式有：

18世紀(jì)：馬青公式（Machin's Formula）

1706年，英國數(shù)學(xué)家約翰·馬青（John Machin）提出了一個(gè)快速計(jì)算π的公式：

π=16tan??1(1/5)?4tan??1(1/239)\pi = 16 \tan^{-1} (1/5) - 4 \tan^{-1} (1/239)π=16tan?1(1/5)?4tan?1(1/239)

這個(gè)公式可以讓數(shù)學(xué)家用級數(shù)展開的方式快速計(jì)算π，而不需要使用古代的割圓術(shù)方法。從18世紀(jì)到19世紀(jì)，數(shù)學(xué)家們不斷改進(jìn)這類公式，比如：

• 約翰遜（John W. Wrench）等人手工計(jì)算π到808位（1946年），當(dāng)時(shí)是世界紀(jì)錄。
• 手工計(jì)算一次 π 的100位數(shù)值，大約需要數(shù)個(gè)月的時(shí)間，計(jì)算錯誤率較高。高斯-勒讓德算法（Gauss-Legendre Algorithm）——倍增精度
• 先選兩個(gè)數(shù)，一個(gè)代表圓的半徑，另一個(gè)代表多邊形的周長。
• 不斷用數(shù)學(xué)公式調(diào)整這兩個(gè)數(shù)，使它們越來越接近真實(shí)的π值。
• 這個(gè)方法每計(jì)算一次，π的精度就能翻倍，計(jì)算速度比割圓術(shù)快得多。但在現(xiàn)代計(jì)算機(jī)問世之前，即便有了更好的公式，對π計(jì)算的精度和速度比割圓術(shù)快了許多，傳統(tǒng)的手工計(jì)算依然是很復(fù)雜緩慢的。一直到1948年，英國的弗格森和美國的倫奇共同發(fā)表了π的808位小數(shù)值，成為人工計(jì)算圓周率值的最高紀(jì)錄 ?，F(xiàn)代計(jì)算機(jī)問世，對π的計(jì)算才有了質(zhì)的飛躍。

到了1970年代，數(shù)學(xué)家發(fā)現(xiàn)了一種指數(shù)級提高精度的方法：

• 先選兩個(gè)數(shù)，一個(gè)代表圓的半徑，另一個(gè)代表多邊形的周長。
• 不斷用數(shù)學(xué)公式調(diào)整這兩個(gè)數(shù)，使它們越來越接近真實(shí)的π值。
• 這個(gè)方法每計(jì)算一次，π的精度就能翻倍，計(jì)算速度比割圓術(shù)快得多。

但在現(xiàn)代計(jì)算機(jī)問世之前，即便有了更好的公式，對π計(jì)算的精度和速度比割圓術(shù)快了許多，傳統(tǒng)的手工計(jì)算依然是很復(fù)雜緩慢的。一直到1948年，英國的弗格森和美國的倫奇共同發(fā)表了π的808位小數(shù)值，成為人工計(jì)算圓周率值的最高紀(jì)錄 。現(xiàn)代計(jì)算機(jī)問世，對π的計(jì)算才有了質(zhì)的飛躍。

**計(jì)算機(jī)計(jì)算π：****1949年，第一次使用計(jì)算機(jī)計(jì)算π，**效率突飛猛進(jìn)

1949年，美國人約翰·馮·諾依曼（John von Neumann）和尼古拉斯·梅特羅波利斯（Nicholas Metropolis）首次使用電子計(jì)算機(jī)計(jì)算π。他們使用的計(jì)算機(jī)是ENIAC（電子數(shù)字積分計(jì)算機(jī)），計(jì)算到了2037位，打破了人類歷史上所有的手工計(jì)算紀(jì)錄。

計(jì)算效率提升對比：

• 手工計(jì)算808位 π 需要幾個(gè)月，但ENIAC 用70小時(shí)就算出了2037位，速度提升了數(shù)百倍。
• ENIAC每秒可計(jì)算5000次加法，遠(yuǎn)超人工手算。這次計(jì)算標(biāo)志著計(jì)算機(jī)在數(shù)學(xué)計(jì)算中的首次突破，也讓π的計(jì)算邁入了現(xiàn)代計(jì)算時(shí)代。

1960年代 - 1980年代：計(jì)算機(jī) π 計(jì)算突破百萬位

隨著計(jì)算機(jī)技術(shù)的發(fā)展，數(shù)學(xué)家開始使用更加高效的算法，比如高斯-勒讓德算法（Gauss-Legendre Algorithm）。這使得π的計(jì)算速度指數(shù)級增長，計(jì)算效率提升對比：從1949年 ENIAC 計(jì)算2037位需要70小時(shí)，到1982年 CRAY-1 計(jì)算800萬位僅需5小時(shí)，速度提升了10萬倍以上。

1990年代至今：超級計(jì)算機(jī)計(jì)算π

1987年，數(shù)學(xué)家楚德諾夫斯基兄弟提出了一種超快的計(jì)算方法：

1π=12∑k=0∞(?1)k(6k)!(545140134k+13591409)(3k)!(k!)3(640320)3k+3/2\frac{1}{\pi} = 12 \sum_{k=0}^{\infty} \frac{(-1)^k (6k)! (545140134k + 13591409)}{(3k)! (k!)^3 (640320)^{3k+3/2}}π1=12k=0∑∞(3k)!(k!)3(640320)3k+3/2(?1)k(6k)!(545140134k+13591409)

這個(gè)公式雖然看起來很復(fù)雜，但它的優(yōu)勢是——計(jì)算速度極快！

• 每次計(jì)算，π的精度可以增加100萬位！

• 這是目前計(jì)算機(jī)計(jì)算π時(shí)最常用的方法。進(jìn)入21世紀(jì)后，數(shù)學(xué)家開始使用更加高效的算法（比如楚德諾夫斯基算法）和超級計(jì)算機(jī)，π的計(jì)算速度和精度進(jìn)一步提升。手工計(jì)算1000位 π 需要數(shù)月，計(jì)算機(jī)可以在秒級完成。隨著數(shù)學(xué)公式的優(yōu)化和計(jì)算機(jī)的不斷升級，計(jì)算速度越來越快，計(jì)算效率提升對比：

• 1949年 ENIAC 計(jì)算2037位 π 需要70小時(shí)，2022年超級計(jì)算機(jī)計(jì)算100萬億位僅需157天。

• 計(jì)算位數(shù)從2000位增長到100萬億位，增長了50億倍，計(jì)算時(shí)間卻只增長了1000倍左右，說明計(jì)算機(jī)的計(jì)算能力大幅提升?，F(xiàn)代的π值，計(jì)算機(jī)已經(jīng)完全取代了手工計(jì)算。截至2023年，瑞士的研究團(tuán)隊(duì)已經(jīng)用超級計(jì)算機(jī)將π計(jì)算到100萬億（101?）位，打破了歷史記錄。今天，我們計(jì)算π的主要瓶頸已經(jīng)不再是數(shù)學(xué)公式，而是計(jì)算機(jī)硬件的速度和存儲容量。未來，隨著量子計(jì)算的發(fā)展，π的計(jì)算可能會變得更加高效。

計(jì)算這么多位的π值到底有什么用呢？

其實(shí)，日常生活中，我們用的π通常不會超過3.1416，就連NASA計(jì)算航天器軌道時(shí)，也只用到15位（3.14159265358979）。但是，計(jì)算超高精度的π仍然有很多意義：

• 測試超級計(jì)算機(jī)的性能：計(jì)算π需要大量計(jì)算資源，能測試計(jì)算機(jī)的處理能力和穩(wěn)定性。
• 數(shù)學(xué)研究：數(shù)學(xué)家想知道，π的無窮小數(shù)部分是否真的完全隨機(jī)，或者是否有隱藏的數(shù)學(xué)規(guī)律。
• 科學(xué)工程：在某些精密科學(xué)研究（比如量子計(jì)算、黑洞模擬）中，需要超高精度的π值。盡管我們可能永遠(yuǎn)不需要全部的100萬億位π，但計(jì)算它的過程，推動了數(shù)學(xué)、計(jì)算機(jī)科學(xué)和工程學(xué)的發(fā)展。

π的精確值計(jì)算仍在繼續(xù)，未來沒有止境

因?yàn)棣惺莻€(gè)無理數(shù)，因此是小數(shù)點(diǎn)是不循環(huán)無限多的，這就決定了π值的精確值求索永無止境。從阿基米德的割圓術(shù)到超級計(jì)算機(jī)的楚德諾夫斯基算法，人類對π的追求已經(jīng)持續(xù)了2000多年。今天，π的計(jì)算精度已經(jīng)遠(yuǎn)超實(shí)際需求，但科學(xué)家們?nèi)匀辉诓粩嗵魬?zhàn)極限，不僅僅是為了計(jì)算π，而是為了推動數(shù)學(xué)和計(jì)算機(jī)科學(xué)的發(fā)展。

由此，π不僅僅是一個(gè)數(shù)學(xué)常數(shù)，它是人類智慧和科技進(jìn)步的象征。未來，人類將會用更先進(jìn)的方法，把π計(jì)算得更遠(yuǎn)，或許能揭開它更深層的秘密！

時(shí)空通訊原創(chuàng)文章，請尊重作者版權(quán)，感謝閱讀。