C语言大数浮点数除法计算结果异常问题咨询

原因分析 #

问题核心是float类型的精度限制：

• float仅拥有24位有效二进制位（含1位隐含位），对应约7-8位十进制有效数字。86997171是8位十进制数，已经触及float的精度上限，转换为float时会被近似存储，无法保留精确值。
• 实际转换后，86997171对应的float值是86997168（这是float能精确表示的最接近原数的整数），用这个数除以48，结果刚好是1812441.0，所以ceil计算后还是1812441，和预期的1812442不符。
• 而3/2的数值较小，float能精确表示1.5，因此计算结果符合预期。

解决办法 #

有两种可靠的解决思路：

1. 改用double类型 #

double拥有53位有效二进制位，对应15-17位十进制有效数字，完全能精确表示86997171这类数，计算时不会丢失精度。修改后的代码如下：

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

int main(void)
{
    // 86997171 / 48 = 1812441.0625
    int max = ceil((double) 86997171 / 48);
    double max2 = ((double) 86997171)/ 48;
    printf("max = %i, max2 = %lf\n", max, max2);
    int max3 = ceil((double) 3 / 2);
    double max4 = ((double) 3) / 2;
    printf("max3 = %i, max4 = %lf\n", max3, max4);
}

运行后会输出正确结果：

max = 1812442, max2 = 1812441.062500
max3 = 2, max4 = 1.500000

2. 使用整数运算实现ceil #

对于正整数除法的ceil需求，完全可以用整数运算避免浮点误差，公式为(被除数 + 除数 - 1) / 除数。示例代码：

#include <stdio.h>

int main(void)
{
    int a = 86997171;
    int b = 48;
    int max = (a + b - 1) / b;
    printf("max = %i\n", max);
    
    int a2 = 3;
    int b2 = 2;
    int max3 = (a2 + b2 -1)/b2;
    printf("max3 = %i\n", max3);
}

这种方法既没有浮点误差，运算效率也更高。

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