求π的神级外星人程序(800位)

今天在书上偶尔看到了一个能够精确求出π的前800位的神级程序，令人惊奇的是此程序的主体只有3行。

此程序的作者不详，大家都叫外星人程序。首先膜拜一下。。。。

以下为源程序：

#include <iostream>using namespace std;int main() {long a = 10000, b = 0, c = 2800, d, e = 0, f[2801], g;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); return 0;}

我看网上有人对此程序分析过，有一个相对来说比较好理解的程序（其实也不好理解，菜鸟一个。。。）：

#include <iostream>using namespace std;long a=10000,b,c=2800,d,e,f[2801];int main() {for(b = 0; b < c; b++)f[b]=a/5;for(;c > 0; c -= 14) {d = 0;for(b=c;b >0;b--){d = d*b + f[b]*a;f[b] = d % (2*b-1);d /= 2*b-1;}printf("%.4d", e+d/a);e = d%a;}return 0;}

注意：以上的输出格式为"%.4d"，否则会出现细微的错误。（具体原因不详）
如果想要详细了解求解π的原理及解释，请参考：求圆周率π一万位程序分析