精确计算PI小数点后800位小数

近来，在网上发现一个能精确计算PI后面数百位的C程序，不敢独享，把其翻译成PASCAL版，供大家参考，该函数写的非常精巧，仅用了三行代码，就能精确计算出PI后800位小数。小弟不才，对数学造诣不深，不明白他算法的原理是什么，如果有懂的请不吝赐规教。下面将该函数的C版和PASCAL版一并贴出，供大家参考。

///////////////////////////////////////////////////////////
// C语言源程序功能：精确计算PI小数点后799位小数
// 转自 http://www.programfan.net/club/showbbs.asp?id=16434
//
// #include

 1<stdlib.h>   
 2// #include <stdio.h>   
 3// long a=10000,b,c=2800,d,e,f[2801],g;   
 4// main()   
 5// {   
 6// for(;b-c;) f[b++]=a/5;   
 7// for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)   
 8// for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);   
 9// }   
10///////////////////////////////////////////////////////////   
11// 版权所有，转载请注明出处   
12// 该C函数的 Delphi 版   
13// 编译： 黄涛 [email protected]   
14// 2004-9-7   
15function CalcPI :string;   
16//long a=10000,b,c=2800,d,e,f[2801],g;   
17var a,b,c,d,e,g:integer;   
18f :array[0..2800] of integer;   
19begin   
20a:=10000;b:=0;c:=2800;e:=0;f[2800]:=0;   
21//for(;b-c;) f[b++]=a/5;   
22while (b-c)&lt;&gt;0 do begin f[b]:= a div 5; inc(b) end;   
23//for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)   
24while c&lt;&gt;0 do begin   
25d:=0; g:=c*2; b:=c;   
26repeat   
27//for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);   
28d:=d+f[b]*a;dec(g);f[b]:=d mod g;   
29d:=d div g;dec(g);dec(b);   
30if b&lt;&gt;0 then d:=d*b;   
31until b=0;   
32c:=c-14;   
33Result:= Result + format('%.4d',[e+ d div a]);   
34e:=d mod a;   
35end;   
36end; 
37
38通过上面的代码可以看出，C语言非常灵活，写出的代码的确很精练，不过他的编程风格确实不值得提倡，太难懂了，非C语言高手很难搞懂他写的什么。不过如果有谁也能写出更精练的Delphi代码也请贴出，供大家学习。我写的只是为了翻译，基本上是一一对应的，没有做更多的优化。</stdio.h></stdlib.h>