一个圆周率的算法

#include <windows.h>
#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);  
   return 0;
}

在网上看到的.