经典PI求解C程序

#include <stdio.h>
#include <stdlib.h>
long a=10000, b, c=2800, d, e, f[2801], g;

void main()
{
 for(;b<2800;) 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);
 scanf("%s");
}