【图】精确计算PI小数点后800位小数

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

　　///////////////////////////////////////////////////////////
　　// C语言源程序功能：精确计算PI小数点后799位小数
　　// 转自 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);
　　// }
　　///////////////////////////////////////////////////////////
　　// 版权所有，转载请注明出处
　　// 该C函数的 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 c0 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 b0 then d:=d*b;
　　    until b=0;
　　    c:=c-14;
　　    Result:= Result + format('%.4d',[e+ d div a]);
　　    e:=d mod a;
　　 end;
　　end;