c – 用Chudnovsky算法计算Pi数

发布时间：2020-12-16 07:03:22 所属栏目：百科来源：网络整理

导读：所以这里有一些根据Chudnovsky算法计算Pi的C代码： #include iostream#include cmath#include iomanipint fact(int digit){ int result = 1; for (int i=digit; i1; i--) { result *= i; } return result;}int main(){ long double pi = 0.0; int precision

所以这里有一些根据Chudnovsky算法计算Pi的C代码：

#include <iostream>
#include <cmath>
#include <iomanip>

int fact(int digit)
{
    int result = 1;
    for (int i=digit; i>1; i--)
    {
        result *= i;
    }
    return result;
}

int main()
{
    long double pi = 0.0;
    int precision = 10;
    for ( int k=0; k<=precision; k++)
    {
        pi += ( pow(-1,k)*fact((6*k))*(13591409+545140134*k))/(fact(3*k)*pow(fact(k),3)*pow(640320,(3*k+1.5)));
    }
    pi = 1.0/(12*pi);
    std::cout<<std::setprecision(80)<<pi<<'n'<<M_PI<<'n';
}

问题是这段代码只返回了15位正确的PI编号,另一部分是垃圾…当精度变量设置为12或更大时,它返回-nan.任何人都可以帮助我吗？ TIA
编辑：Ubuntu 12.10 x64,gcc 4.7.2

解决方法

在fact()中使用double作为fact()返回类型以及“result”变量.

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!