求π的神级外星人程序(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"，否则会出现细微的错误。（具体原因不详）
如果想要详细了解求解π的原理及解释，请参考：求圆周率π一万位程序分析