C++实现digital_root时0加6结果为5的异常问题排查

问题背景 #

• 开发场景：使用C++完成Codewars平台的数字根（digital root）算法练习题，测试用例要求计算167346的数字根，正确结果应为9，但代码首次运行就出现0加6计算结果为5的异常。

问题代码 #

int digital_root(int n, int rtn = 0) {
   double fract(0), intval(0), dn(n);
   dn = dn / 10.0;
   fract = modf(dn, &intval);

   std::cout << "rtn before : " << rtn << std::endl;
   std::cout << "fract : " << fract << std::endl;
   std::cout << "fract times 10 : " << fract*10 << std::endl;

   rtn+=(fract*10.0);

   std::cout << "rtn plus fract times 10 : " << rtn << std::endl;

   if (rtn < 10 && intval == 0)
     return rtn;
   else if (rtn > 10 && intval == 0)
     return digital_root(rtn, 0);

   return digital_root(intval, rtn);
}

运行日志 #

rtn before : 0
fract : 0.6
fract times 10 : 6
rtn plus fract times 10 : 5
rtn before : 5
fract : 0.4
fract times 10 : 4
rtn plus fract times 10 : 9
rtn before : 9
fract : 0.3
fract times 10 : 3
rtn plus fract times 10 : 12
rtn before : 12
fract : 0.7
fract times 10 : 7
rtn plus fract times 10 : 18
rtn before : 18
fract : 0.6
fract times 10 : 6
rtn plus fract times 10 : 24
rtn before : 24
fract : 0.1
fract times 10 : 1
rtn plus fract times 10 : 25
rtn before : 0
fract : 0.5
fract times 10 : 5
rtn plus fract times 10 : 5
rtn before : 5
fract : 0.2
fract times 10 : 2
rtn plus fract times 10 : 7

异常原因 #

该异常由浮点数精度误差+隐式类型转换截断共同导致：

• 二进制浮点数无法精确表示全部十进制小数，代码中计算得到的fract实际并非精确值0.6，而是无限接近0.6但略小于0.6的二进制近似值（典型值为0.59999999999999986）。
• cout默认仅输出6位有效数字，打印浮点数时会自动做四舍五入，因此日志中显示fract : 0.6、fract times 10 : 6，但fract*10的真实值是略小于6的近似值（典型值为5.999999999999998）。
• 将该浮点数结果累加到int类型变量rtn时，C++会执行浮点数到整数的隐式转换，转换规则为直接丢弃小数部分、不做四舍五入，5.99999999999经截断后得到整数5，最终出现0加6结果为5的异常。
• 代码还存在其他同类隐患：递归调用时将double类型的intval直接传入int类型形参，同样会触发精度截断问题，本质是错误使用浮点数处理纯整数运算场景。

修复方向：完全移除浮点数运算逻辑，改用整数取余%、整除/操作拆分每一位数字，可从根源避免这类精度问题；数字根本身存在O(1)时间复杂度的数学公式，无需递归即可计算。

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