#include<iostream> #include<cstdio> #include<algorithm> using namespace std; int a=10000,b,c=2800,d,e,f[2801],g; int 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); return 0; } 文章来源: 很牛的圆周率算法(800位)