计算派(圆周率)近似值的超小程序

开始看到时觉得很不可思议！

#include <stdio.h>

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