使用求和公式近似π值的C++代码错误排查求助

1. 编译器警告与阶乘函数的问题 #

• 窄化转换警告：你写的factorial(2.0 * i + 1)中，2.0*i会生成double类型值，传入要求unsigned int参数的factorial函数时，会触发double转int的窄化警告。只需把2.0 * i改成整数运算2*i，就能得到整数类型的参数，消除警告。
• 整数溢出：factorial函数返回int类型，但阶乘增长极快——13的阶乘就超过了32位int的最大值（2147483647），会导致溢出，计算结果完全错误。应该把阶乘函数的返回类型改为double，或者直接在求和函数中递推计算每一项，避免单独计算大阶乘。

2. 求和公式的逻辑错误 #

• 项的计算错误：你使用的factorial(i)/factorial(2i+1)不符合π的近似求和公式（推测是基于arcsin(1)的泰勒展开，π=2*arcsin(1)）。正确的每一项应该是( (2i)! )/(4^i * (i!)^2 * (2i+1))，直接计算大阶乘容易溢出，更高效的方式是递推每一项：每一项等于前一项乘以(2i-1)/(2i*(2i+1))（i从1开始），完全不用单独计算阶乘。
• 返回值错误：你写的return (2*pi_half + 1 )逻辑错误，正确的应该是return 2 * pi_half——因为pi_half已经初始化为1.0（对应展开式的第一项），后续累加的是其余项，而π=2arcsin(1)=2(1 + Σ各项)。

3. 循环范围问题 #

循环条件i < m意味着输入m时，只会循环m-1次。如果m是你要计算的总项数（包括初始的1），这个逻辑是对的；如果m是额外要累加的项数，需要调整为i <= m。

修正后的代码（递推版，无溢出风险） #

#include <iostream>
using namespace std;

double calculatePi(unsigned int m) {
    double pi_half = 1.0;
    double current_term = 1.0; // 用于递推的初始项
    for (unsigned int i = 1; i < m; ++i) {
        // 递推计算当前项：前一项 * (2i-1)/(2i*(2i+1))
        current_term *= (2*i - 1.0) / (2*i * (2*i + 1.0));
        pi_half += current_term;
    }
    return 2 * pi_half;
}

int main() {
    unsigned int m;
    cout << "请输入计算所用的项数m: ";
    cin >> m;
    cout << "近似π值为: " << calculatePi(m) << endl;
    return 0;
}

保留阶乘函数的修正版（适合理解公式） #

如果一定要保留阶乘函数，需要修改返回类型为double，同时递推计算4^i避免溢出：

#include <iostream>
using namespace std;

double factorial(unsigned int f);
double sum2(unsigned int m);

double sum2(unsigned int m = 1) {
    double pi_half = 1.0;
    double power4 = 1.0; // 递推计算4^i，初始为4^0=1
    for (unsigned int i = 1; i < m; ++i) {
        power4 *= 4;
        double numerator = factorial(2*i);
        double denominator = power4 * factorial(i) * factorial(i) * (2*i + 1);
        pi_half += numerator / denominator;
    }
    return 2 * pi_half;
}

double factorial(unsigned int f) {
    double fact = 1.0;
    for (unsigned int i = f; i > 0; --i) {
        fact *= i;
    }
    return fact;
}

int main() {
    unsigned int m;
    cout << "请输入计算所用的项数m: ";
    cin >> m;
    cout << sum2(m) << endl;
    return 0;
}

内容的提问来源于stack exchange，提问作者lolamay