如何在控制台完整输出__float128类型的所有小数位？

问题分析 #

你的代码存在两个核心问题：

1. 阶乘计算用long int存储，会快速溢出（n≥20时就会超出范围），导致e的计算精度完全丧失；
2. C++标准库的cout没有原生支持__float128类型的输出，直接输出会编译失败。

修正方案 #

1. 修复阶乘计算的溢出问题

将阶乘函数的内部变量和返回值全部改为__float128类型，避免整数溢出：

__float128 fact(unsigned int n) {
    if (n == 0 || n == 1) return 1.0Q; // Q后缀标识__float128字面量
    __float128 fac = 1.0Q;
    for (unsigned int x = n; x > 1; x--) {
        fac *= x;
    }
    return fac;
}

2. 实现__float128的高精度输出

GCC配套的quadmath库提供了quadmath_snprintf函数，可将__float128格式化为字符串，再通过cout输出，完整保留精度。

完整可运行代码 #

#include <iostream>
#include <quadmath.h>

__float128 fact(unsigned int n) {
    if (n == 0 || n == 1) return 1.0Q;
    __float128 fac = 1.0Q;
    for (unsigned int x = n; x > 1; x--) {
        fac *= x;
    }
    return fac;
}

int main() {
    __float128 e = 0.0Q;
    // __float128约有34位有效数字，循环到x=30已足够收敛到精度上限
    for (unsigned int x = 0; x <= 30; x++) {
        e += 1.0Q / fact(x);
    }

    // 格式化__float128为字符串，保留全部有效数字
    char buf[256];
    quadmath_snprintf(buf, sizeof(buf), "%.34Qf", e);
    std::cout << "高精度常数e：" << buf << std::endl;

    return 0;
}

编译说明 #

编译时需要链接quadmath库，使用以下命令：

g++ -o e_calc e_calc.cpp -lquadmath

内容的提问来源于stack exchange，提问作者cum-bria