一个很牛的计算pi的c程序

#include   
#include   

long 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);  
}

计算结果： 
3141592653589793238462643383279502884197169399375105820974944592307816406
2862089986280348253421170679821480865132823066470938446095505822317253594
0812848111745028410270193852110555964462294895493038196442881097566593344
6128475648233786783165271201909145648566923460348610454326648213393607260
2491412737245870066063155881748815209209628292540917153643678925903600113
3053054882046652138414695194151160943305727036575959195309218611738193261
1793105118548074462379962749567351885752724891227938183011949129833673362
4406566430860213949463952247371907021798609437027705392171762931767523846
7481846766940513200056812714526356082778577134275778960917363717872146844
0901224953430146549585371050792279689258923542019956112129021960864034418
1598136297747713099605187072113499999983729780499510597317328160963185