一段计算派(∏)的代码

一段计算派(∏、圆周率) 的代码

    long a=10000,b,c=2800,d,e,f[2801],g;
    char r[1000];
    char t[6];
    memset(r,0,sizeof(r));
    b=d=e=g=0;
    for(;b-c;)
        {
        f[b++]=a/5;
        }
    for(;d=0,g=c*2;c-=14,sprintf(t,"%.4d/0",e+d/a),e=d%a,strcat(r,t))
        {
        for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
        }
       
运行结束后，数组r内容是PAI的前800位(31415926535......)，网上找到的但是我没看懂原理-_-!!!