c++求圆周率pi的估算值

在C++中，求圆周率π的估算值可以使用多种算法，例如著名的Leibniz公式、Chudnovsky算法或者是蒙特卡罗方法等。这里我们简单举例使用最基本的迭代公式，如Bailey–Borwein–Plouffe (BBP) 公式，这个公式可以直接计算π的小数部分：

#include <iostream>
#include <cmath>

double piEstimate(int n_terms) {
    double numerator = 4;
    double denominator = 1;
    double term = 1;
    double pi = 0;

    for (int i = 0; i < n_terms; ++i) {
        pi += numerator / denominator;
        term *= -4 * i / (2 * i + 1);
        numerator *= (6 * i + 5);
        denominator *= (2 * i + 1) * (2 * i + 3);
    }
    
    return pi;
}

int main() {
    int terms = 10000; // 要求的精度，更多项会得到更精确的结果
    double estimatedPi = piEstimate(terms);
    std::cout << "Estimated value of π with " << terms << " terms: " << estimatedPi << std::endl;