C语言浮点运算下溢错误：f(x)=√(x²+1)-1适配大小数值实现问题

1. 猜测正误说明 #

你对浮点精度的理解存在偏差：float标注的6位精度指的是有效数字位数，而非小数点后6位。1e-7的有效数字仅有1位，本身远未到float的精度下限，不会直接被截断为0。
你遇到的本质是浮点计算的灾难性抵消问题：当x极小的时候，x²的量级远小于1，计算1 + x²时，double类型（53位二进制有效位，对应约15~17位十进制有效数字）会丢失x²的有效信息，导致sqrt(1 + x²)的结果直接等于1，最终减1得到0。这是算法逻辑导致的精度损失，和x本身的存储精度无关，所以你单纯切换long double也解决不了根本问题，只是推迟了问题出现的x量级。

2. 现有代码的其他问题 #

• 计算平方不要用pow(x,2)：pow是通用幂次计算函数，专门算平方直接写x*x精度更高、性能更好。
• 转换long double的操作无效：代码里的1.是double类型常量，计算时会把long double类型的g降为double计算，要使用long double常量需要写成1.L。
• 极大值近似逻辑有边界隐患：如果输入的x是负数，x-1的结果会完全错误，建议先取绝对值计算再处理符号。

3. 最优解决方案 #

不用换更高精度的浮点类型，直接对原式做分子有理化即可从算法层面消除抵消问题：
原表达式变形后为：
f(x) = x² / (sqrt(x² + 1) + 1)
转换后的表达式在x极小的时候，分母近似等于2，结果近似为x²/2，完全没有精度损失；x极大的时候，分母近似等于|x| + 1，结果近似等于|x|，和原表达式的极限完全一致，甚至可以不需要分支判断覆盖全量程。

4. 修改后的代码示例 #

#include <math.h>
#include <stdio.h>

double feval(double x) {
    double x_sq = x * x;
    // 直接用有理化后的表达式计算，全量程适用
    return x_sq / (sqrt(x_sq + 1.0) + 1.0);
}

int main(void)
{
    double x;
    printf("Input: ");
    scanf("%lf", &x);
    double b = feval(x);
    printf("%f\n", b);
    return 0;
}

如果确实需要保留极大值的近似分支优化，可以调整为：

double feval(double x) {
    const double MAX_THRESH = 1e299;
    if (fabs(x) > MAX_THRESH) {
        // 原表达式极大值极限为|x| - 1，可根据判题要求调整
        return x > 0 ? x - 1 : -x - 1;
    }
    double x_sq = x * x;
    return x_sq / (sqrt(x_sq + 1.0) + 1.0);
}

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