很牛的圆周率算法(800位)

圆周率:     3.1415926535897932384626433832795028841971693993

 

#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;
}