关于圆周率的难懂程序，懂pi算法的高手请来看看,！

[DreamSeeker]

#include <stdio.h>
int a=10000,b,c=2800,d,e,f[2801],g;
int 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);
}
这是一个三行就可以计算圆周率到999位的程序，可是小弟怎么也看不懂，他用的是什么算法？具体怎么实现的？如果想计算更多位应该怎么改？请看得懂的兄弟详细解释一下，不胜感谢！！

[mathe]

http://www.csdn.net/expert/TopicView.asp?id=335556