圆周率的计算（C语言版，绝对精华）

#include<stdio.h>
long 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);
 scanf("%s");
 return 0;
}

太牛的程序了，算法了不起，哪位老兄可以解释一下，让小弟长长见识！