近來，在網(wǎng)上發(fā)現(xiàn)一個能精確計算PI后面數(shù)百位的C程序，不敢獨享，把其翻譯成PASCAL版，供大家參考，該函數(shù)寫的非常精巧，僅用了三行代碼，就能精確計算出PI后800位小數(shù)。小弟不才，對數(shù)學(xué)造詣不深，不明白他算法的原理是什么，如果有懂的請不吝賜規(guī)教。下面將該函數(shù)的C版和PASCAL版一并貼出，供大家參考。

///////////////////////////////////////////////////////////
// C語言 源程序  功能：精確計算PI小數(shù)點后799位小數(shù)
// 轉(zhuǎn)自 http://www.PRogramfan.net/club/showbbs.asp?id=16434
//
//  #include <stdlib.h>
//  #include <stdio.h>
//  long a=10000,b,c=2800,d,e,f[2801],g;
//  main()
//  {
//    for(;b-c;) f[b++]=a/5;
//    for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
//      for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
//  }
///////////////////////////////////////////////////////////
// 版權(quán)所有，轉(zhuǎn)載請注明出處
// 該C函數(shù)的 Delphi 版
// 編譯： 黃濤  Hunter@Shentong.com.cn
//  2004-9-7
function CalcPI :string;
//long a=10000,b,c=2800,d,e,f[2801],g;
var a,b,c,d,e,g:integer;
    f :array[0..2800] of integer;
begin
  a:=10000;b:=0;c:=2800;e:=0;f[2800]:=0;
  //for(;b-c;) f[b++]=a/5;
  while (b-c)<>0 do begin f[b]:= a div 5; inc(b) end;
  //for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
  while c<>0 do begin
    d:=0; g:=c*2; b:=c;
    repeat                                 
      //for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
      d:=d+f[b]*a;dec(g);f[b]:=d mod g;
      d:=d div g;dec(g);dec(b);
      if b<>0 then d:=d*b;
    until b=0;
    c:=c-14;
    Result:= Result + format('%.4d',[e+ d div a]);
    e:=d mod a;
  end;
end;

     通過上面的代碼可以看出，C語言非常靈活，寫出的代碼的確很精練，不過他的編程風(fēng)格確實不值得提倡，太難懂了，非C語言高手很難搞懂他寫的什么。不過如果有誰也能寫出更精練的Delphi代碼也請貼出，供大家學(xué)習(xí)。我寫的只是為了翻譯，基本上是一一對應(yīng)的，沒有做更多的優(yōu)化。