九州算术数学论坛's Archiver

Math 发表于 2007-7-16 18:19

一个很牛的计算Л(PI)的c程序

[code]#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);  
}  [/code]

计算结果:
[free]314159265358979323846264338327950288
419716939937510582097494459230781640
628620899862803482534211706798214808
651328230664709384460955058223172535
940812848111745028410270193852110555
964462294895493038196442881097566593
344612847564823378678316527120190914
564856692346034861045432664821339360
726024914127372458700660631558817488
152092096282925409171536436789259036
001133053054882046652138414695194151
160943305727036575959195309218611738
193261179310511854807446237996274956
735188575272489122793818301194912983
367336244065664308602139494639522473
719070217986094370277053921717629317
675238467481846766940513200056812714
526356082778577134275778960917363717
872146844090122495343014654958537105
079227968925892354201995611212902196
086403441815981362977477130996051870
721134999999837297804995105973173281
60963185  [/free]

jznhljg 发表于 2008-4-1 20:21

飞燕?........

页: [1]

Powered by Discuz! Archiver 7.0.0  © 2001-2007 Comsenz Inc.