计算圆周率pi，C++里用拉马努金公式求解

拉马努金的公式：

 

使用的macos xcode软件编译：

 代码如下：

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
long long fact(int n) { //递归函数，阶乘的递归
    if (n <= 1) {
        return 1;
    } else {
        return n * fact(n - 1);
    }
}

int main() {
    int x = 10; //次数
    for (int i = 0; i < x; ++i) {
        long double pi = 0.0;

        for (int k = 0; k < i; ++k) {
            long double term = static_cast<long double>(fact(4 * k)) * (1103 + 26390 * k) / (pow(fact(k), 4) * pow(396, 4 * k)); //pow函数计算幂
            pi += term;
        }

        pi *= (2 * sqrt(2)) / 9801;
        pi = 1 / pi;
        
        cout << "第" <<i + 1<<"项: " << setprecision(20) << pi << endl;
    }

    return 0;
}

运行结果：

​​​​​​​ 

 

圆周率前几十位为：3.14159 26535 89793 23846 26433；上述代码从第四项开始就不能更精确了，可能是long double只能精确到小数点后15位。