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

王朝delphi·作者佚名  2006-01-09
窄屏简体版  字體: |||超大  

近来,在网上发现一个能精确计算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 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语言非常灵活,写出的代码的确很精练,不过他的编程风格确实不值得提倡,太难懂了,非C语言高手很难搞懂他写的什么。不过如果有谁也能写出更精练的Delphi代码也请贴出,供大家学习。我写的只是为了翻译,基本上是一一对应的,没有做更多的优化。

 
 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
 
 
© 2005- 王朝網路 版權所有 導航