一个很牛的求圆周率的6行代码

int 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);
} 

很早以前找到过这段代码，感觉蛮牛的，有兴趣的可以研究下。