有人能为我解释一下用这种方法计算圆周率的原理和三条for语句该如何理解吗

Question

#include"stdio.h"
#include"stdlib.h"

long a =10000, b = 0, c = 10000,d,e,f[10001],g;

int main(int argc, char const *argv[])
{
    for(; b != c; f[b++] = a/5);
    for(; d = 0,g = c*2,c>5000;c--,printf("%.4d",e+d/a),e = d%a)
    for(b = c;d += f[b]*a,f[b]=d % --g, d /= g--, --b; d*=b );

    system("pause");//这行和上面的第二个include是我自己加的。
    return 0;
}

我发现我对于for语句的了解，对于c语言的了解还很困乏。这是一个算圆周率的代码，我看不懂这是咋跑起来的。有人能为我解释一下用这种方法计算圆周率的原理和三条for语句该如何理解吗？
麻烦大家了。

Answer 1

回答不易，回答有用请你采纳，点击回答右侧采纳即可！
这段代码是用来计算圆周率的。圆周率是一个无理数，可以通过算法来近似计算。这里使用的是一种叫做"莫比乌斯级数"的方法。

第一个for语句是初始化数组f，使得f[b] = a/5 (其中b为循环变量，a为10000)。

第二个for语句是主要的循环，c从10000开始递减，直到5000。在每次循环中，通过计算d和e的值，输出e+d/a的值，并将e的值赋值为d mod a。

第三个for语句是在第二个for语句中计算d、f[b]和g的值。b从c开始，d += f[b]a，f[b]=d mod --g， d /= g--, --b， d=b。

这个程序的输出是近似的圆周率值。