求一个用c语言算圆周率的程序

#include
long a=10000,b,c=2800,d,e,f[2801],g;
void main()
{
for(;b-c;) f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf("%04d",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
}
这个问题我问过.
能算800位,我给的网址是我问的,回答里有解释,和能算1000位的优化算法,很好的.