圓周率是一個(gè)極其馳名的數(shù)。從有文字記載的歷史開(kāi)始，這個(gè)數(shù)就引進(jìn)了外行人和學(xué)者們的興趣。作為一個(gè)非常重要的常數(shù)，圓周率最早是出于解決有關(guān)圓的計(jì)算問(wèn)題。僅憑這一點(diǎn)，求出它的盡量準(zhǔn)確的近似值，就是一個(gè)極其迫切的問(wèn)題了。事實(shí)也是如此，幾千年來(lái)作為數(shù)學(xué)家們的奮斗目標(biāo)，古今中外一代一代的數(shù)學(xué)家為此獻(xiàn)出了自己的智慧和勞動(dòng)。回顧歷史，人類(lèi)對(duì) π 的認(rèn)識(shí)過(guò)程，反映了數(shù)學(xué)和計(jì)算技術(shù)發(fā)展情形的一個(gè)側(cè)面。 π 的研究，在一定程度上反映這個(gè)地區(qū)或時(shí)代的數(shù)學(xué)水平。德國(guó)數(shù)學(xué)史家康托說(shuō)："歷史上一個(gè)國(guó)家所算得的圓周率的準(zhǔn)確程度，可以作為衡量這個(gè)國(guó)家當(dāng)時(shí)數(shù)學(xué)發(fā)展水平的指標(biāo)。"直到19世紀(jì)初，求圓周率的值應(yīng)該說(shuō)是數(shù)學(xué)中的頭號(hào)難題。為求得圓周率的值，人類(lèi)走過(guò)了漫長(zhǎng)而曲折的道路，它的歷史是饒有趣味的。我們可以將這一計(jì)算歷程分為幾個(gè)階段。

　　通過(guò)實(shí)驗(yàn)對(duì) π 值進(jìn)行估算，這是計(jì)算 π 的的第一階段。這種對(duì) π 值的估算基本上都是以觀察或?qū)嶒?yàn)為根據(jù)，是基于對(duì)一個(gè)圓的周長(zhǎng)和直徑的實(shí)際測(cè)量而得出的。在古代世界，實(shí)際上長(zhǎng)期使用 π ＝3這個(gè)數(shù)值。最早見(jiàn)于文字記載的有基督教《圣經(jīng)》中的章節(jié)，其上取圓周率為3。這一段描述的事大約發(fā)生在公元前950年前后。其他如巴比倫、印度、中國(guó)等也長(zhǎng)期使用3這個(gè)粗略而簡(jiǎn)單實(shí)用的數(shù)值。在我國(guó)劉徽之前"圓徑一而周三"曾廣泛流傳。我國(guó)第一部《周髀算經(jīng)》中，就記載有圓"周三徑一"這一結(jié)論。在我國(guó)，木工師傅有兩句從古流傳下來(lái)的口訣：叫做："周三徑一，方五斜七"，意思是說(shuō)，直徑為1的圓，周長(zhǎng)大約是3，邊長(zhǎng)為5的正方形，對(duì)角線(xiàn)之長(zhǎng)約為7。這正反映了早期人們對(duì)圓周率 π 和√2 這兩個(gè)無(wú)理數(shù)的粗略估計(jì)。東漢時(shí)期官方還明文規(guī)定圓周率取3為計(jì)算面積的標(biāo)準(zhǔn)。后人稱(chēng)之為"古率"。

　　早期的人們還使用了其它的粗糙方法。如古埃及、古希臘人曾用谷粒擺在圓形上，以數(shù)粒數(shù)與方形對(duì)比的方法取得數(shù)值�；蛴脛蛑啬景邃彸蓤A形和方形以秤量對(duì)比取值……由此，得到圓周率的稍好些的值。如古埃及人應(yīng)用了約四千年的 4 (8/9)² = 3.1605。在印度，公元前六世紀(jì)，曾取 π= √10 = 3.162。在我國(guó)東、西漢之交，新朝王莽令劉歆制造量的容器??律嘉量斛。劉歆在制造標(biāo)準(zhǔn)容器的過(guò)程中就需要用到圓周率的值。為此，他大約也是通過(guò)做實(shí)驗(yàn)，得到一些關(guān)于圓周率的并不劃一的近似值�，F(xiàn)在根據(jù)銘文推算，其計(jì)算值分別取為3.1547，3.1992，3.1498，3.2031比徑一周三的古率已有所進(jìn)步。人類(lèi)的這種探索的結(jié)果，當(dāng)主要估計(jì)圓田面積時(shí)，對(duì)生產(chǎn)沒(méi)有太大影響，但以此來(lái)制造器皿或其它計(jì)算就不合適了。

　　憑直觀推測(cè)或?qū)嵨锒攘�，�?lái)計(jì)算 π 值的實(shí)驗(yàn)方法所得到的結(jié)果是相當(dāng)粗略的。

　　真正使圓周率計(jì)算建立在科學(xué)的基礎(chǔ)上，首先應(yīng)歸功于阿基米德。他是科學(xué)地研究這一常數(shù)的第一個(gè)人，是他首先提出了一種能夠借助數(shù)學(xué)過(guò)程而不是通過(guò)測(cè)量的、能夠把 π 的值精確到任意精度的方法。由此，開(kāi)創(chuàng)了圓周率計(jì)算的第二階段。

圓周長(zhǎng)大于內(nèi)接正四邊形而小于外切正四邊形，因此 2√2 ＜ π ＜ 4 。
當(dāng)然，這是一個(gè)差勁透頂?shù)睦�子。�?jù)說(shuō)阿基米德用到了正96邊形才算出他的值域。

　　阿基米德求圓周率的更精確近似值的方法，體現(xiàn)在他的一篇論文《圓的測(cè)定》之中。在這一書(shū)中，阿基米德第一次創(chuàng)用上、下界來(lái)確定 π 的近似值，他用幾何方法證明了"圓周長(zhǎng)與圓直徑之比小于 3+(1/7) 而大于 3 + (10/71) "，他還提供了誤差的估計(jì)。重要的是，這種方法從理論上而言，能夠求得圓周率的更準(zhǔn)確的值。到公元150年左右，希臘天文學(xué)家托勒密得出 π ＝3.1416，取得了自阿基米德以來(lái)的巨大進(jìn)步。

割圓術(shù)。不斷地利用勾股定理，來(lái)計(jì)算正N邊形的邊長(zhǎng)。

　　在我國(guó)，首先是由數(shù)學(xué)家劉徽得出較精確的圓周率。公元263年前后，劉徽提出著名的割圓術(shù)，得出 π ＝3.14，通常稱(chēng)為"徽率"，他指出這是不足近似值。雖然他提出割圓術(shù)的時(shí)間比阿基米德晚一些，但其方法確有著較阿基米德方法更美妙之處。割圓術(shù)僅用內(nèi)接正多邊形就確定出了圓周率的上、下界，比阿基米德用內(nèi)接同時(shí)又用外切正多邊形簡(jiǎn)捷得多。另外，有人認(rèn)為在割圓術(shù)中劉徽提供了一種絕妙的精加工辦法，以致于他將割到192邊形的幾個(gè)粗糙的近似值通過(guò)簡(jiǎn)單的加權(quán)平均，竟然獲得具有4位有效數(shù)字的圓周率 π ＝3927/1250 ＝3.1416。而這一結(jié)果，正如劉徽本人指出的，如果通過(guò)割圓計(jì)算得出這個(gè)結(jié)果，需要割到3072邊形。這種精加工方法的效果是奇妙的。這一神奇的精加工技術(shù)是割圓術(shù)中最為精彩的部分，令人遺憾的是，由于人們對(duì)它缺乏理解而被長(zhǎng)期埋沒(méi)了。

　　恐怕大家更加熟悉的是祖沖之所做出的貢獻(xiàn)吧。對(duì)此，《隋書(shū)?律歷志》有如下記載："宋末，南徐州從事祖沖之更開(kāi)密法。以圓徑一億為丈，圓周盈數(shù)三丈一尺四寸一分五厘九毫二秒七忽，?數(shù)三丈一尺四寸一分五厘九毫二秒六忽，正數(shù)在盈?二限之間。密率：圓徑一百一十三，圓周三百五十五。約率，圓徑七，周二十二。"

　　這一記錄指出，祖沖之關(guān)于圓周率的兩大貢獻(xiàn)。其一是求得圓周率

　　　　3.1415926 ＜ π ＜ 3.1415927

　　其二是，得到 π 的兩個(gè)近似分?jǐn)?shù)即：約率為22／7；密率為355／113。

　　他算出的 π 的8位可靠數(shù)字，不但在當(dāng)時(shí)是最精密的圓周率，而且保持世界記錄九百多年。以致于有數(shù)學(xué)史家提議將這一結(jié)果命名為"祖率"。

　　這一結(jié)果是如何獲得的呢？追根溯源，正是基于對(duì)劉徽割圓術(shù)的繼承與發(fā)展，祖沖之才能得到這一非凡的成果。因而當(dāng)我們稱(chēng)頌祖沖之的功績(jī)時(shí)，不要忘記他的成就的取得是因?yàn)樗�站在�?shù)學(xué)偉人劉徽的肩膀上的緣故。后人曾推算若要單純地通過(guò)計(jì)算圓內(nèi)接多邊形邊長(zhǎng)的話(huà)，得到這一結(jié)果，需要算到圓內(nèi)接正12288邊形，才能得到這樣精確度的值。祖沖之是否還使用了其它的巧妙辦法來(lái)簡(jiǎn)化計(jì)算呢？這已經(jīng)不得而知，因?yàn)橛涊d其研究成果的著作《綴術(shù)》早已失傳了。這在中國(guó)數(shù)學(xué)發(fā)展史上是一件極令人痛惜的事。

中國(guó)發(fā)行的祖沖之紀(jì)念郵票

　　祖沖之的這一研究成果享有世界聲譽(yù)：巴黎"發(fā)現(xiàn)宮"科學(xué)博物館的墻壁上著文介紹了祖沖之求得的圓周率，莫斯科大學(xué)禮堂的走廊上鑲嵌有祖沖之的大理石塑像，月球上有以祖沖之命名的環(huán)形山……

　　對(duì)于祖沖之的關(guān)于圓周率的第二點(diǎn)貢獻(xiàn)，即他選用兩個(gè)簡(jiǎn)單的分?jǐn)?shù)尤其是用密率來(lái)近似地表示 π 這一點(diǎn)，通常人們不會(huì)太注意。然而，實(shí)際上，后者在數(shù)學(xué)上有更重要的意義。

　　密率與 π 的近似程度很好，但形式上卻很簡(jiǎn)單，并且很優(yōu)美，只用到了數(shù)字1、3、5。數(shù)學(xué)史家梁宗巨教授驗(yàn)證出：分母小于16604的一切分?jǐn)?shù)中，沒(méi)有比密率更接近 π 的分?jǐn)?shù)。在國(guó)外，祖沖之死后一千多年，西方人才獲得這一結(jié)果。

　　可見(jiàn)，密率的提出是一件很不簡(jiǎn)單的事情。人們自然要追究他是采用什么辦法得到這一結(jié)果的呢？他是用什么辦法把圓周率從小數(shù)表示的近似值化為近似分?jǐn)?shù)的呢？這一問(wèn)題歷來(lái)為數(shù)學(xué)史家所關(guān)注。由于文獻(xiàn)的失傳，祖沖之的求法已不為人知。后人對(duì)此進(jìn)行了各種猜測(cè)。

　　讓我們先看看國(guó)外歷史上的工作，希望能夠提供出一些信息。

　　1573年，德國(guó)人奧托得出這一結(jié)果。他是用阿基米德成果22／7與托勒密的結(jié)果377／120用類(lèi)似于加成法"合成"的：(377-22) / (120-7) = 355/113。

　　1585年，荷蘭人安托尼茲用阿基米德的方法先求得：333/106 ＜ π ＜ 377/120，用兩者作為 π 的母近似值，分子、分母各取平均，通過(guò)加成法獲得結(jié)果：3 ((15+17)/(106+120) = 355/113。

　　兩個(gè)雖都得出了祖沖之密率，但使用方法都為偶合，無(wú)理由可言。

　　在日本，十七世紀(jì)關(guān)孝和重要著作《括要算法》卷四中求圓周率時(shí)創(chuàng)立零約術(shù)，其實(shí)質(zhì)就是用加成法來(lái)求近似分?jǐn)?shù)的方法。他以3、4作為母近似值，連續(xù)加成六次得到祖沖之約率，加成一百十二次得到密率。其學(xué)生對(duì)這種按部就班的笨辦法作了改進(jìn)，提出從相鄰的不足、過(guò)剩近似值就近加成的辦法，（實(shí)際上就是我們前面已經(jīng)提到的加成法）這樣從3、4出發(fā)，六次加成到約率，第七次出現(xiàn)25／8，就近與其緊鄰的22／7加成，得47／15，依次類(lèi)推，只要加成23次就得到密率。

　　錢(qián)宗琮先生在《中國(guó)算學(xué)史》（1931年）中提出祖沖之采用了我們前面提到的由何承天首創(chuàng)的"調(diào)日法"或稱(chēng)加權(quán)加成法。他設(shè)想了祖沖之求密率的過(guò)程：以徽率157／50，約率22／7為母近似值，并計(jì)算加成權(quán)數(shù)x=9，于是 (157 + 22×，9) / (50+7×9) = 355/113，一舉得到密率。錢(qián)先生說(shuō)："沖之在承天后，用其術(shù)以造密率，亦意中事耳。"

　　另一種推測(cè)是：使用連分?jǐn)?shù)法。

　　由于求二自然數(shù)的最大公約數(shù)的更相減損術(shù)遠(yuǎn)在《九章算術(shù)》成書(shū)時(shí)代已流行，所以借助這一工具求近似分?jǐn)?shù)應(yīng)該是比較自然的。于是有人提出祖沖之可能是在求得盈 二數(shù)之后，再使用這個(gè)工具，將3.14159265表示成連分?jǐn)?shù)，得到其漸近分?jǐn)?shù)：3，22／7，333／106，355／113，102573／32650…

　　最后，取精確度很高但分子分母都較小的355／113作為圓周率的近似值。至于上面圓周率漸近分?jǐn)?shù)的具體求法，這里略掉了。你不妨利用我們前面介紹的方法自己求求看。英國(guó)李約瑟博士持這一觀點(diǎn)。他在《中國(guó)科學(xué)技術(shù)史》卷三第19章幾何編中論祖沖之的密率說(shuō)："密率的分?jǐn)?shù)是一個(gè)連分?jǐn)?shù)漸近數(shù)，因此是一個(gè)非凡的成就。"

　　我國(guó)再回過(guò)頭來(lái)看一下國(guó)外所取得的成果。

　　1150年，印度數(shù)學(xué)家婆什迦羅第二計(jì)算出 π= 3927/1250 = 3.1416。1424年，中亞細(xì)亞地區(qū)的天文學(xué)家、數(shù)學(xué)家卡西著《圓周論》，計(jì)算了3×2²⁸＝805,306,368邊內(nèi)接與外切正多邊形的周長(zhǎng)，求出 π 值，他的結(jié)果是：

　　 π＝3.14159265358979325

　　有十七位準(zhǔn)確數(shù)字。這是國(guó)外第一次打破祖沖之的記錄。

　　16世紀(jì)的法國(guó)數(shù)學(xué)家韋達(dá)利用阿基米德的方法計(jì)算 π 近似值，用 6×2¹⁶正邊形，推算出精確到9位小數(shù)的 π 值。他所采用的仍然是阿基米德的方法，但韋達(dá)卻擁有比阿基米德更先進(jìn)的工具：十進(jìn)位置制。17世紀(jì)初，德國(guó)人魯?shù)婪蛴昧藥缀跻簧�的時(shí)間鉆研這個(gè)問(wèn)題。他也將新的十進(jìn)制與早的阿基米德方法結(jié)合起來(lái)，但他不是從正六邊形開(kāi)始并將其邊數(shù)翻番的，他是從正方形開(kāi)始的，一直推導(dǎo)出了有2⁶²條邊的正多邊形，約4,610,000,000,000,000,000邊形！這樣，算出小數(shù)35位。為了記念他的這一非凡成果，在德國(guó)圓周率 π 被稱(chēng)為"魯?shù)婪驍?shù)"。但是，用幾何方法求其值，計(jì)算量很大，這樣算下去，窮數(shù)學(xué)家一生也改進(jìn)不了多少。到魯?shù)婪蚩梢哉f(shuō)已經(jīng)登峰造極，古典方法已引導(dǎo)數(shù)學(xué)家們走得很遠(yuǎn)，再向前推進(jìn)，必須在方法上有所突破。

　　17世紀(jì)出現(xiàn)了數(shù)學(xué)分析，這銳利的工具使得許多初等數(shù)學(xué)束手無(wú)策的問(wèn)題迎刃而解。 π 的計(jì)算歷史也隨之進(jìn)入了一個(gè)新的階段。

　　這一時(shí)期人們開(kāi)始擺脫求多邊形周長(zhǎng)的繁難計(jì)算，利用無(wú)窮級(jí)數(shù)或無(wú)窮連乘積來(lái)算 π 。

　　1593年，韋達(dá)給出

 　　這一不尋常的公式是 π 的最早分析表達(dá)式。甚至在今天，這個(gè)公式的優(yōu)美也會(huì)令我們贊嘆不已。它表明僅僅借助數(shù)字2，通過(guò)一系列的加、乘、除和開(kāi)平方就可算出 π 值。

　　接著有多種表達(dá)式出現(xiàn)。如沃利斯1650年給出：


　　1706年，梅欽建立了一個(gè)重要的公式，現(xiàn)以他的名字命名：


　　再利用分析中的級(jí)數(shù)展開(kāi)，他算到小數(shù)后100位。

　　這樣的方法遠(yuǎn)比可憐的魯?shù)婪蛴么蟀肷鷷r(shí)間才摳出的35位小數(shù)的方法簡(jiǎn)便得多。顯然，級(jí)數(shù)方法宣告了古典方法的過(guò)時(shí)。此后，對(duì)于圓周率的計(jì)算像馬拉松式競(jìng)賽，紀(jì)錄一個(gè)接著一個(gè)：

　　1844年，達(dá)塞利用公式：


　　算到200位。

　　19世紀(jì)以后，類(lèi)似的公式不斷涌現(xiàn)， π 的位數(shù)也迅速增長(zhǎng)。1873年，謝克斯利用梅欽的一系列方法，級(jí)數(shù)公式將 π 算到小數(shù)后707位。為了得到這項(xiàng)空前的紀(jì)錄，他花費(fèi)了二十年的時(shí)間。他死后，人們將這凝聚著他畢生心血的數(shù)值，銘刻在他的墓碑上，以頌揚(yáng)他頑強(qiáng)的意志和堅(jiān)韌不拔的毅力。于是在他的墓碑上留下了他一生心血的結(jié)晶： π 的小數(shù)點(diǎn)后707位數(shù)值。這一驚人的結(jié)果成為此后74年的標(biāo)準(zhǔn)。此后半個(gè)世紀(jì)，人們對(duì)他的計(jì)算結(jié)果深信不疑，或者說(shuō)即便懷疑也沒(méi)有辦法來(lái)檢查它是否正確。以致于在1937年巴黎博覽會(huì)發(fā)現(xiàn)館的天井里，依然顯赫地刻著他求出的 π 值。

　　又過(guò)了若干年，數(shù)學(xué)家弗格森對(duì)他的計(jì)算結(jié)果產(chǎn)生了懷疑，其疑問(wèn)基于如下猜想：在 π 的數(shù)值中，盡管各數(shù)字排列沒(méi)有規(guī)律可循，但是各數(shù)碼出現(xiàn)的機(jī)會(huì)應(yīng)該相同。當(dāng)他對(duì)謝克斯的結(jié)果進(jìn)行統(tǒng)計(jì)時(shí)，發(fā)現(xiàn)各數(shù)字出現(xiàn)次數(shù)過(guò)于參差不齊。于是懷疑有誤。他使用了當(dāng)時(shí)所能找到的最先進(jìn)的計(jì)算工具，從1944年5月到1945年5月，算了整整一年。1946年，弗格森發(fā)現(xiàn)第528位是錯(cuò)的（應(yīng)為4，誤為5）。謝克斯的值中足足有一百多位全都報(bào)了銷(xiāo)，這把可憐的謝克斯和他的十五年浪費(fèi)了的光陰全部一筆勾銷(xiāo)了。

　　對(duì)此，有人曾嘲笑他說(shuō)：數(shù)學(xué)史在記錄了諸如阿基米德、費(fèi)馬等人的著作之余，也將會(huì)擠出那么一、二行的篇幅來(lái)記述1873年前謝克斯曾把 π 計(jì)算到小數(shù)707位這件事。這樣，他也許會(huì)覺(jué)得自己的生命沒(méi)有虛度。如果確實(shí)是這樣的話(huà)，他的目的達(dá)到了。

　　人們對(duì)這些在地球的各個(gè)角落里作出不懈努力的人感到不可理解，這可能是正常的。但是，對(duì)此做出的嘲笑卻是過(guò)于殘忍了。人的能力是不同的，我們無(wú)法要求每個(gè)人都成為費(fèi)馬、高斯那樣的人物。但成為不了偉大的數(shù)學(xué)家，并不意味著我們就不能為這個(gè)社會(huì)做出自己有限的貢獻(xiàn)。人各有其長(zhǎng)，作為一個(gè)精力充沛的計(jì)算者，謝克斯愿意獻(xiàn)出一生的大部分時(shí)光從事這項(xiàng)工作而別無(wú)報(bào)酬，并最終為世上的知識(shí)寶庫(kù)添了一小塊磚加了一個(gè)塊瓦。對(duì)此我們不應(yīng)為他的不懈努力而感染并從中得到一些啟發(fā)與教育嗎？

　　1948年1月弗格森和倫奇兩人共同發(fā)表有808位正確小數(shù)的 π 。這是人工計(jì)算 π 的最高記錄。

　　1946年，世界第一臺(tái)計(jì)算機(jī)ENIAC制造成功，標(biāo)志著人類(lèi)歷史邁入了電腦時(shí)代。電腦的出現(xiàn)導(dǎo)致了計(jì)算方面的根本革命。1949年，ENIAC根據(jù)梅欽公式計(jì)算到2035（一說(shuō)是2037）位小數(shù)，包括準(zhǔn)備和整理時(shí)間在內(nèi)僅用了70小時(shí)。計(jì)算機(jī)的發(fā)展一日千里，其記錄也就被頻頻打破。

ENIAC：一個(gè)時(shí)代的開(kāi)始

　　1973年，有人就把圓周率算到了小數(shù)點(diǎn)后100萬(wàn)位，并將結(jié)果印成一本二百頁(yè)厚的書(shū)，可謂世界上最枯燥無(wú)味的書(shū)了。1989年突破10億大關(guān)，1995年10月超過(guò)64億位。1999年9月30日，《文摘報(bào)》報(bào)道，日本東京大學(xué)教授金田康正已求到2061.5843億位的小數(shù)值。如果將這些數(shù)字打印在A4大小的復(fù)印紙上，令每頁(yè)印2萬(wàn)位數(shù)字，那么，這些紙摞起來(lái)將高達(dá)五六百米。來(lái)自最新的報(bào)道：金田康正利用一臺(tái)超級(jí)計(jì)算機(jī)，計(jì)算出圓周率小數(shù)點(diǎn)后一兆二千四百一十一億位數(shù)，改寫(xiě)了他本人兩年前創(chuàng)造的紀(jì)錄。據(jù)悉，金田教授與日立制作所的員工合作，利用目前計(jì)算能力居世界第二十六位的超級(jí)計(jì)算機(jī)，使用新的計(jì)算方法，耗時(shí)四百多個(gè)小時(shí)，才計(jì)算出新的數(shù)位，比他一九九九年九月計(jì)算出的小數(shù)點(diǎn)后二千六百一十一位提高了六倍。圓周率小數(shù)點(diǎn)后第一兆位數(shù)是二，第一兆二千四百一十一億位數(shù)為五。如果一秒鐘讀一位數(shù)，大約四萬(wàn)年后才能讀完。

　　不過(guò)，現(xiàn)在打破記錄，不管推進(jìn)到多少位，也不會(huì)令人感到特別的驚奇了。實(shí)際上，把 π 的數(shù)值算得過(guò)分精確，應(yīng)用意義并不大�，F(xiàn)代科技領(lǐng)域使用的 π 值，有十幾位已經(jīng)足夠。如果用魯?shù)婪虻?5位小數(shù)的 π 值計(jì)算一個(gè)能把太陽(yáng)系包圍起來(lái)的圓的周長(zhǎng)，誤差還不到質(zhì)子直徑的百萬(wàn)分之一。我們還可以引美國(guó)天文學(xué)家西蒙?紐克姆的話(huà)來(lái)說(shuō)明這種計(jì)算的實(shí)用價(jià)值：

　　"十位小數(shù)就足以使地球周界準(zhǔn)確到一英寸以?xún)?nèi)，三十位小數(shù)便能使整個(gè)可見(jiàn)宇宙的四周準(zhǔn)確到連最強(qiáng)大的顯微鏡都不能分辨的一個(gè)量。"

　　那么為什么數(shù)學(xué)家們還象登山運(yùn)動(dòng)員那樣，奮力向上攀登，一直求下去而不是停止對(duì) π 的探索呢？為什么其小數(shù)值有如此的魅力呢？

　　這其中大概免不了有人類(lèi)的好奇心與領(lǐng)先于人的心態(tài)作怪，但除此之外，還有許多其它原因。

奔騰與圓周率之間的奇妙關(guān)系……

　　1、它現(xiàn)在可以被人們用來(lái)測(cè)試或檢驗(yàn)超級(jí)計(jì)算機(jī)的各項(xiàng)性能，特別是運(yùn)算速度與計(jì)算過(guò)程的穩(wěn)定性。這對(duì)計(jì)算機(jī)本身的改進(jìn)至關(guān)重要。就在幾年前，當(dāng)Intel公司推出奔騰（Pentium）時(shí)，發(fā)現(xiàn)它有一點(diǎn)小問(wèn)題，這問(wèn)題正是通過(guò)運(yùn)行 π 的計(jì)算而找到的。這正是超高精度的 π 計(jì)算直到今天仍然有重要意義的原因之一。

　　2、 計(jì)算的方法和思路可以引發(fā)新的概念和思想。雖然計(jì)算機(jī)的計(jì)算速度超出任何人的想象，但畢竟還需要由數(shù)學(xué)家去編制程序，指導(dǎo)計(jì)算機(jī)正確運(yùn)算。實(shí)際上，確切地說(shuō)，當(dāng)我們把 π 的計(jì)算歷史劃分出一個(gè)電子計(jì)算機(jī)時(shí)期時(shí)，這并非意味著計(jì)算方法上的改進(jìn)，而只是計(jì)算工具有了一個(gè)大飛躍而已。因而如何改進(jìn)計(jì)算技術(shù)，研究出更好的計(jì)算公式，使公式收斂得更快、能極快地達(dá)到較大的精確度仍是數(shù)學(xué)家們面對(duì)的一個(gè)重要課題。在這方面，本世紀(jì)印度天才數(shù)學(xué)家拉馬努揚(yáng)得出了一些很好的結(jié)果。他發(fā)現(xiàn)了許多能夠迅速而精確地計(jì)算 π 近似值的公式。他的見(jiàn)解開(kāi)通了更有效地計(jì)算 π 近似值的思路。現(xiàn)在計(jì)算機(jī)計(jì)算 π 值的公式就是由他得到的。至于這位極富傳奇色彩的數(shù)學(xué)家的故事，在這本小書(shū)中我們不想多做介紹了。不過(guò)，我希望大家能夠明白 π 的故事講述的是人類(lèi)的勝利，而不是機(jī)器的勝利。

　　3、還有一個(gè)關(guān)于 π 的計(jì)算的問(wèn)題是：我們能否無(wú)限地繼續(xù)算下去？答案是：不行！根據(jù)朱達(dá)偌夫斯基的估計(jì)，我們最多算10⁷⁷位。雖然，現(xiàn)在我們離這一極限還相差很遠(yuǎn)很遠(yuǎn)，但這畢竟是一個(gè)界限。為了不受這一界限的約束，就需要從計(jì)算理論上有新的突破。前面我們所提到的計(jì)算，不管用什么公式都必須從頭算起，一旦前面的某一位出錯(cuò)，后面的數(shù)值完全沒(méi)有意義。還記得令人遺憾的謝克斯嗎？他就是歷史上最慘痛的教訓(xùn)。

　　4、于是，有人想能否計(jì)算時(shí)不從頭開(kāi)始，而是從半截開(kāi)始呢？這一根本性的想法就是尋找并行算法公式。1996年，圓周率的并行算法公式終于找到，但這是一個(gè)16進(jìn)位的公式，這樣很容易得出的1000億位的數(shù)值，只不過(guò)是16進(jìn)位的。是否有10進(jìn)位的并行計(jì)算公式，仍是未來(lái)數(shù)學(xué)的一大難題。

　　5、作為一個(gè)無(wú)窮數(shù)列，數(shù)學(xué)家感興趣的把 π 展開(kāi)到上億位，能夠提供充足的數(shù)據(jù)來(lái)驗(yàn)證人們所提出的某些理論問(wèn)題，可以發(fā)現(xiàn)許多迷人的性質(zhì)。如，在 π 的十進(jìn)展開(kāi)中，10個(gè)數(shù)字，哪些比較稀，哪些比較密？ π 的數(shù)字展開(kāi)中某些數(shù)字出現(xiàn)的頻率會(huì)比另一些高嗎？或許它們并非完全隨意？這樣的想法并非是無(wú)聊之舉。只有那些思想敏銳的人才會(huì)問(wèn)這種貌似簡(jiǎn)單，許多人司空見(jiàn)慣但卻不屑發(fā)問(wèn)的問(wèn)題。

　　6、數(shù)學(xué)家弗格森最早有過(guò)這種猜想：在 π 的數(shù)值式中各數(shù)碼出現(xiàn)的概率相同。正是他的這個(gè)猜想為發(fā)現(xiàn)和糾正向克斯計(jì)算 π 值的錯(cuò)誤立下了汗馬功勞。然而，猜想并不等于現(xiàn)實(shí)。弗格森想驗(yàn)證它，卻無(wú)能為力。后人也想驗(yàn)證它，也是苦于已知的 π 值的位數(shù)太少。甚至當(dāng)位數(shù)太少時(shí)，人們有理由對(duì)猜想的正確性做出懷疑。如，數(shù)字0的出現(xiàn)機(jī)會(huì)在開(kāi)始時(shí)就非常少。前50位中只有1個(gè)0，第一次出現(xiàn)在32位上�？墒�，這種現(xiàn)象隨著數(shù)據(jù)的增多，很快就改變了：100位以?xún)?nèi)有8個(gè)0；200位以?xún)?nèi)有19個(gè)0；……1000萬(wàn)位以?xún)?nèi)有999，440個(gè)0；……60億位以?xún)?nèi)有599，963，005個(gè)0，幾乎占1／10。

　　其他數(shù)字又如何呢？結(jié)果顯示，每一個(gè)都差不多是1／10，有的多一點(diǎn)，有的少一點(diǎn)。雖然有些偏差，但都在1／10000之內(nèi)。

　　7、人們還想知道： π 的數(shù)字展開(kāi)真的沒(méi)有一定的模式嗎？我們希望能夠在十進(jìn)制展開(kāi)式中通過(guò)研究數(shù)字的統(tǒng)計(jì)分布，尋找任何可能的模型??如果存在這種模型的話(huà)，迄今為止尚未發(fā)現(xiàn)有這種模型。同時(shí)我們還想了解： π 的展開(kāi)式中含有無(wú)窮的樣式變化嗎？或者說(shuō)，是否任何形式的數(shù)字排列都會(huì)出現(xiàn)呢？著名數(shù)學(xué)家希爾伯特在沒(méi)有發(fā)表的筆記本中曾提出下面的問(wèn)題： π 的十進(jìn)展開(kāi)中是否有10個(gè)9連在一起？以現(xiàn)在算到的60億位數(shù)字來(lái)看，已經(jīng)出現(xiàn)：連續(xù)6個(gè)9連在一起。希爾伯特的問(wèn)題答案似乎應(yīng)該是肯定的，看來(lái)任何數(shù)字的排列都應(yīng)該出現(xiàn)，只是什么時(shí)候出現(xiàn)而已。但這還需要更多 π 的數(shù)位的計(jì)算才能提供切實(shí)的證據(jù)。

　　8、在這方面，還有如下的統(tǒng)計(jì)結(jié)果：在60億數(shù)字中已出現(xiàn)連在一起的8個(gè)8；9個(gè)7；10個(gè)6；小數(shù)點(diǎn)后第710150位與3204765位開(kāi)始，均連續(xù)出現(xiàn)了七個(gè)3；小數(shù)點(diǎn)52638位起連續(xù)出現(xiàn)了14142135這八個(gè)數(shù)字，這恰是的前八位；小數(shù)點(diǎn)后第2747956位起，出現(xiàn)了有趣的數(shù)列876543210，遺憾的是前面缺個(gè)9；還有更有趣的數(shù)列123456789也出現(xiàn)了。

　　如果繼續(xù)算下去，看來(lái)各種類(lèi)型的數(shù)字列組合可能都會(huì)出現(xiàn)。

　　在1777年出版的《或然性算術(shù)實(shí)驗(yàn)》一書(shū)中，蒲豐提出了用實(shí)驗(yàn)方法計(jì)算 π 。這個(gè)實(shí)驗(yàn)方法的操作很簡(jiǎn)單：找一根粗細(xì)均勻，長(zhǎng)度為 d 的細(xì)針，并在一張白紙上畫(huà)上一組間距為 l 的平行線(xiàn)（方便起見(jiàn)，常取 l = d/2），然后一次又一次地將小針任意投擲在白紙上。這樣反復(fù)地投多次，數(shù)數(shù)針與任意平行線(xiàn)相交的次數(shù)，于是就可以得到 π 的近似值。因?yàn)槠沿S本人證明了針與任意平行線(xiàn)相交的概率為 p = 2l/πd 。利用這一公式，可以用概率方法得到圓周率的近似值。在一次實(shí)驗(yàn)中，他選取 l = d/2 ，然后投針2212次，其中針與平行線(xiàn)相交704次，這樣求得圓周率的近似值為 2212/704 = 3.142。當(dāng)實(shí)驗(yàn)中投的次數(shù)相當(dāng)多時(shí)，就可以得到 π 的更精確的值。

　　1850年，一位叫沃爾夫的人在投擲5000多次后，得到 π 的近似值為3.1596。目前宣稱(chēng)用這種方法得到最好結(jié)果的是意大利人拉茲瑞尼。在1901年，他重復(fù)這項(xiàng)實(shí)驗(yàn)，作了3408次投針，求得 π 的近似值為3.1415929，這個(gè)結(jié)果是如此準(zhǔn)確，以致于很多人懷疑其實(shí)驗(yàn)的真?zhèn)�。如美�?guó)猶他州奧格登的國(guó)立韋伯大學(xué)的L?巴杰就對(duì)此提出過(guò)有力的質(zhì)疑。

　　不過(guò)，蒲豐實(shí)驗(yàn)的重要性并非是為了求得比其它方法更精確的 π 值。蒲豐投針問(wèn)題的重要性在于它是第一個(gè)用幾何形式表達(dá)概率問(wèn)題的例子。計(jì)算 π 的這一方法，不但因其新穎，奇妙而讓人叫絕，而且它開(kāi)創(chuàng)了使用隨機(jī)數(shù)處理確定性數(shù)學(xué)問(wèn)題的先河，是用偶然性方法去解決確定性計(jì)算的前導(dǎo)。

　　在用概率方法計(jì)算 π 值中還要提到的是：R?查特在1904年發(fā)現(xiàn)，兩個(gè)隨意寫(xiě)出的數(shù)中，互素的概率為6／π²。1995年4月英國(guó)《自然》雜志刊登文章，介紹英國(guó)伯明翰市阿斯頓大學(xué)計(jì)算機(jī)科學(xué)與應(yīng)用數(shù)學(xué)系的羅伯特?馬修斯，如何利用夜空中亮星的分布來(lái)計(jì)算圓周率。馬修斯從100顆最亮的星星中隨意選取一對(duì)又一對(duì)進(jìn)行分析，計(jì)算它們位置之間的角距。他檢查了100萬(wàn)對(duì)因子，據(jù)此求得 π 的值約為3.12772。這個(gè)值與真值相對(duì)誤差不超過(guò)5％。

　　通過(guò)幾何、微積分、概率等廣泛的范圍和渠道發(fā)現(xiàn) π ，這充分顯示了數(shù)學(xué)方法的奇異美。 π 竟然與這么些表面看來(lái)風(fēng)馬牛不相及的試驗(yàn)，溝通在一起，這的確使人驚訝不已。