经典算法求PI

Ayala

蒲丰投针,其实这算法本质是个求定积分的算法

1. 
   
2. #include <windows.h>
   
3. #define N 10000
   
4. #define C 20
   
5. #define R 19
   
6. 
   
7. __int64 a=1000000, b=0, c=N*R, d, e, f[N*R+1], g;   
   
8. main()   
   
9. {   
   
10. for (;b-c;) f[b++]=C*N;
    
11. for(;d=0,g=c*2;c-=R,printf("%.6d",e+d/a),e=d%a)   
    
12. for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);   
    
13. system("pause");
    
14. return 0;   
    
15. }   

复制代码

0xAA55

无力吐槽abcdefg，不过代码还真简短

千鹤诗音

我今天在http://outofmemory.cn/code-snipp ... c-algorithm-PI-tell 上也看到一个求pi的算法了呢。不过好长，囧。

乘简

N久前看到过，3行代码计算PI的东西，你这个有点长了。。。

Ayala

乘简 发表于 2017-5-4 12:45
N久前看到过，3行代码计算PI的东西，你这个有点长了。。。

当然有点长 因为对原代码进行了词义解释 方便计算精度的扩展

0xAA55

乘简 发表于 2017-5-4 12:45
N久前看到过，3行代码计算PI的东西，你这个有点长了。。。

代码长度不能决定其优秀程度

乘简

原来代码是一样的。。。

乘简

0xAA55 发表于 2017-5-5 07:12
代码长度不能决定其优秀程度

不好意思，发错了。。。

乘简

0xAA55 发表于 2017-5-5 07:12
代码长度不能决定其优秀程度

不好意思，是我看错了。。。

(⊙o⊙)

这就是精髓的

watermelon

蒙特卡洛方法与数值模拟第一课就是讲的蒲丰式投针实验，

yzw92

发帖辛苦，谢谢分享