國外有位程式設計師童鞋想做一個輸出 π 值的程式,嘗試過很多種演算法。不過不管用哪種演算法,他最後都有一個困惑:如何才能確定計算的 π 值是正確的呢?
檢視下文之前,大家可以先想想。來看看 π 值計算世界紀錄保持者 Mysticial 的回覆。
問題(來自Ishan Sharma):
我曾試過用各種不同的方法,來做一個按序輸出顯示 π 值的程式。我試過泰勒級數法,但是它收斂地非常慢(這是之後,當我將我的結果和網上的值對比時發現的)。總之,我正試著找出更好的演算法。
所以,當我在寫這個程式時總遇到問題,不管使用什麼演算法,我都有一個困惑:我如何才能知道我計算的數字是正確的呢?
最高票數的回答(來自Mysticial):
由於我是計算 π 值的世界紀錄保持者(10 萬億位),我來談一下我的拙見:
除非你是想創造一個新的世界記錄,否則通常的解決辦法就是對比已知值驗證計算值。所以,這很簡單。
其實,我就有一個列出了數字片段的網頁。它就是為了驗證計算值。http://www.numberworld.org/digits/Pi/
但是,當你進入世界記錄的範圍,就沒有已知值讓你對比驗證了。
通常,驗證計算值的標準方法是用另一種演算法再計算一次。所以,只要有一個計算錯誤,數字就不會相同。
這經常會花兩倍多的時間(因為第二種演算法一般會慢一點)。但這是驗證計算值的唯一方法,只要你涉及了從沒被計算過的未知領域或是一個新的世界記錄,就只能這麼做。
在超級計算機保持世界記錄的時候,有兩種不同的常見AGM演算法:
它們都是時間複雜度為 O(N log(N)^2) 的演算法,而且相當簡單。但是,如今的演算法有些不同了。在最近的3個世界記錄中,我們只通過最快的公式(Chudnovsky 公式)進行一次運算,而不是兩次。
這個演算法比較難,但是它比AGM演算法快很多。
然後我們用BBP公式進行數字提取,來驗證這個二進位制數字。
使用這個公式,你可以計算任何二進位制數字,並且不用計算它之前的所有數字。因此,它用來驗證最後少數的二進位制計算值。所以,它比一個完整的計算要快得多。
它的優點是:
1.只需要進行一次複雜的運算。
它的缺點是:
- 需要實現Bailey–Borwein–Plouffe (BBP)公式。
- 需要多一步來驗證:基數從二進位制到十進位制的轉換。
我曾列出過一些細節來說明為什麼驗證最後幾個數字就能證明所有數字都是正確的。但這是顯而易見的,因為任何計算錯誤都會傳遞到最後一個數字。
現在,最後一步(驗證轉換)實際上是非常重要的。一位前世界紀錄保持者曾呼籲我們不要這麼做,因為,我一開始並沒有給出它詳細的過程。
所以,現在我把我部落格的這個片段貼出來:
N = #所需的小數位數
p = 64位素數
用10進位制演算法計算A,用二進位制演算法算B。
如果 A=B,那麼由於“極大可能性”,轉換是正確的。
深入閱讀,可檢視我的這篇博文《Pi – 5 萬億位》。