这是代码:
#include <stdio.h>
long f[2801];
int main()
{
long i = 0, c = 2800, d = 0, e = 0, g = 0;
for (i = 0; i < c; ++i)
f[i] = 2000;
for (;;) {
d = 0;
g = c * 2;
if (!g)
break;
i = c;
for(;;) {
d += f[i] * 10000;
--g;
f[i] = d % g;
d /= g;
--g;
--i;
if (!i) break;
d *= i;
}
printf("%.4ld",e+d/10000);
e = d % 10000;
c -= 14;
}
return 0;
}
我的问题是:这段代码如何计算具有高小数精度的pi以及它使用的数学公式是什么?
最佳答案 这是由荷兰CWI研究所的Dik T. Winter编写的PI程序的格式化副本.最初以两行或三行的模糊形式书写,Dik和其他几种变体根据数学系列的评估输出不同数量的PI(例如800,15,000等).
它是一类称为“套管算法”的程序,旨在输出特定数量的数字.您可以通过谷歌搜索Dik Winter和’spigot算法’了解更多信息.一些示例点击:
Computing Pi in C详细分析了未解决问题的算法.