pi.cpp

一个不错的算法,计算圆周率的PI,快到100位了...好东东大家分享,不过程序虽短,很难看明白的

#include "stdio.h"
void main()
{
    long a=10000,b=0,c=2800,d=0,e=0,f[2801],g=0;
    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 );
 }