是否存在比1.0更适合表示实数1+1e-16的double类型数值？

不存在比1.0更接近实数1+1e-16的double类型数值，这种不精确性是浮点数本身的存储机制特性导致的，并非实现方式的错误。

核心原因：IEEE 754双精度浮点数的精度限制 #

double类型遵循IEEE 754标准，是64位浮点数，其中包含52位尾数（加上隐含的1位，总共53位有效二进制精度），对应的十进制有效精度约为15-17位。

在1.0这个数值点上，相邻两个可表示的double值之间的间隔（称为ulp，单位最后位置）是2^-52，换算成十进制大约是2.22e-16。而你计算的1+1e-16中，1e-16的数值小于这个ulp的一半，意味着它离1.0的距离比离下一个可表示的double值（1 + 2^-52）更近。按照IEEE 754的默认舍入规则（舍入到最近的可表示值，若距离相等则舍入到偶数尾数值），最终结果会被舍入为1.0，不存在其他更接近真实值的double数值可选。

示例验证 #

In [64]: 1+1e-16
Out[64]: 1.0

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