如何提升π计算精度？k取29/30时结果为0的问题排查

AI 小北

2026-8-18

问题分析与修正

核心问题

你遇到的精度不足、k=29/30时输出0的问题，根源在于两点：

类型不匹配与精度浪费
outside被定义为int类型，虽然k≤30时2^k在32位int范围内，但pow(2,k)返回的是double，赋值给int属于不必要的类型转换。后续大整数乘极小浮点数时，结合浮点数精度限制，容易触发舍入为0的情况。
double类型精度上限
double仅能提供约16位十进制有效数字，当k≥29时，经过28/29次sqrt(2+inside)计算后，inside的值会无限接近2，此时2 - inside的差值会小到被double舍入为0，最终sqrt(2-inside)结果为0，导致乘积为0。

修正方案

1. 调整变量类型与计算精度

将inside和outside改为long double类型（通常提供约19位十进制精度），同时使用对应精度的数学函数（powl、sqrtl），避免过早触发精度丢失：

2. 修正后的完整代码

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

using namespace std;

int main() {
    int k;
    long double inside = 0;
    long double outside = 0;

    cout << "Please enter in a value for k..." << endl;
    cin >> k;

    while (k < 1 || k > 30) {
        cout << "Invalid --- k must be >=1 and <=30" << endl;
        cout << "Please enter in a value for k..." << endl;
        cin >> k;
    }

    outside = powl(2, k);
    
    for (int loopNum = 1; loopNum < k; loopNum++) {
        inside = sqrtl(2 + inside);
    }

    inside = sqrtl(2 - inside);
    long double ApproxPI = outside * inside;

    cout << fixed << setprecision(20);
    cout << "Approximation of PI = " << static_cast<double>(ApproxPI) << endl;

    return 0;
}