无法理解这个pi计算算法

Question

我在网站上看到了一个pi计算算法，它看起来像是：

#include <stdio.h>
int a[52514],b,c=52514,d,e,f=1e4,g,h;
main(){
 for(;b=c-=14;h=printf("%04d",e+d/f)){
    for(e=d%=f;g=--b*2;d/=g){
        d=d*b+f*(h?a[b]:f/5);
        a[b]=d%--g;}
   }
}

据说这段代码是基于这种扩展，但我不明白代码和扩展之间的相对关系。

  

pi = sigma（（i！）^ 2 * 2 ^（i + 1））/（2i + 1）！   （i = 0到无限）

你能帮我解释一下吗？谢谢。

Answer 1

pi+3=sigma( (m!)^2 * 2^m * m / (2*m)! ) (m=1 to infinite).

算法的S Pflouffe使用它。