弱问：这几行代码看不懂

shelleycao

这个可以计算pi的值，可是这些10000,560000,56010等这些裸数字是什么意思啊。
这几个短短的代码的思路是什么啊。自己看不懂。//shy~~
# include "stdio.h"
long a=10000,b,c=56000,d,e,f[56010],g;
void main(){
    for(;b-cf[b++]=a/5;
    for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
        for(b=c;d+=f*a,f=d%--g,d/=g--,--b;d*=b);
}

yuxh

原帖由 "shelleycao" 发表：
这个可以计算pi的值，可是这些10000,560000,56010等这些裸数字是什么意思啊。
}

不知道原始公式是怎么样的。。。
我把程序转换了一下，这样比较清楚：

1. # include "stdio.h"
   
2. long a=10000,b,c=5600,d,e,f[5600];
   
3. void main(){
   
4.     for(b = 0; b < c; b++)
   
5.         f[b]=a/5;
   
6.     for(;c >; 0; c -= 14) {
   
7.         d = 0;
   
8.         for(b=c;b >; 0;b--){
   
9.             d = d*b + f[b]*a;
   
10.             f[b] = d % (2*b-1);
    
11.             d /= 2*b-1;
    
12.         }
    
13.         printf("%.4d", e+d/a);
    
14.         e = d%a;
    
15.     }
    
16. }
    

复制代码

只知道10000表示每次取4位精度，56000表示有效数字的位数是56000/14*4=16000位，56010是为了防止运算时数据溢出吧，其实不用，改成56000也可以
很经典的算法啊!数学家才能写得出这样的程序吧？

yuxh

明白了，公式是