用公式求π的近似值c++

在C++中，计算圆周率π的精确值通常涉及到数学库函数，因为π是一个无理数，无法通过简单的算术运算得到精确值。不过我们可以使用一些算法来求得其近似值，如著名的莱布尼茨级数、马赫林公式等。

例如，你可以使用Chudnovsky算法，它是一个快速收敛的级数，用于计算π的近似值。下面是使用C++的一个简单示例：

#include <iostream>
#include <cmath>
#include <complex>

double pi_chudnovsky(int n) {
    double pi = 0;
    for (int k = 0; k <= n; ++k) {
        complex<double> term(4, -8 * std::pow(-1., k));
        term /= factorial(k + 1);
        term *= pow(13591409 + 545140134 * std::pow(k, 3), 1./6);
        pi += term.real();
    }
    return pi;
}

long long factorial(int n) {
    if (n == 0 || n == 1)
        return 1;
    else
        return n * factorial(n - 1);
}

int main() {
    int precision = 100; // 指定求解精度
    double pi_approx = pi_chudnovsky(precision);
    std::cout << "π的近似值（精度到" << precision << "位）: " << pi_approx << std::endl;
    return 0;
}

这个程序会输出π的近似值，数值越大的n，计算结果就越接近真实值。注意这只是一个简化的例子，实际应用中可能需要考虑算法效率、浮点数舍入误差等问题。