如何查找可四舍五入为目标浮点数的最短十进制表示

问题本质 #

你遇到的是IEEE 754浮点数的固有精度限制：单精度float只有24位二进制有效位（23位显式尾数+1位隐含前导1），换算为十进制大约是7~8位有效精度，超出这个长度的十进制数位不会影响最终的浮点值，都会被舍入到最近的可表示浮点值，所以你修改尾数末尾的数位不会改变变量的实际存储值。

查找最短唯一十进制表示的方法 #

直接使用标准库函数（最简便） #

C99及以上标准的格式化输出函数已经内置了对应的能力，你只需要使用对应精度的%g/%e格式符即可：

• 对于float类型，使用%.9g作为格式符：

  printf("%.9g", FLT_MIN);
  // 输出结果为 1.17549435e-38，共9位有效数字
  

• 对于double类型，使用%.17g作为格式符。

这里的9和17是经过计算的固定值：对于有b位二进制有效位的浮点数，保证十进制串转换回浮点数时唯一对应原值的最短有效位数为 ⌈b * log10(2)⌉ + 1，float对应24位二进制有效位，计算得9；double对应53位二进制有效位，计算得17。

你可以验证：将1.17549435e-38f赋值给float变量后再用%.100e打印，输出和直接打印FLT_MIN的结果完全一致，没有任何区别。

手动实现逻辑 #

如果需要自己实现查找逻辑，可以按照以下步骤执行：

• 确定目标浮点数的相邻可表示值：找到比目标值小的最大可表示浮点数prev，和比目标值大的最小可表示浮点数next
• 计算舍入区间：在默认的就近舍入（ties to even）模式下，所有落入开区间 ((prev + target)/2, (target + next)/2) 的十进制数，都会被舍入为目标浮点数
• 在上述区间中搜索有效十进制位数最少的数值，就是你需要的最短表示

目前工业界常用的实现算法是Dragon4、Grisu以及优化版的Ryu算法，这些算法都可以快速输出最短唯一的浮点十进制表示，不需要依赖大整数运算，性能很高。

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