C语言计算高精度圆周率pi程序

这段C语言代码非常简单，原理自己分析吧，可以用于计算高精度的圆周率pi

#include <stdlib.h>
#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);  
}