双精度数值舍入引发取余运算异常，如何解决？

你的问题核心是二进制浮点数的精度局限性：30e-9无法被二进制浮点数精确存储，实际存为29.9999997e-9，转成纳秒后得到29.9999997，模5的结果4.9999997被cout默认的输出精度四舍五入显示为5，但逻辑上你期望的是0（因为原始值应当是30，能被5整除）。

以下是几种实用的解决方法：

• 转整数处理（最可靠）
  既然你在处理纳秒级的时间值，直接将其转换为整数类型可以彻底避免浮点数误差。先对浮点数做四舍五入，再转成整数后取模：

  #include <iostream>
  #include <cmath>
  #include <cstdint>
  
  int main() {
      double s = 30e-9;
      int64_t ns = static_cast<int64_t>(round(s * 1e9)); // 四舍五入为整数纳秒
      int64_t ref_ns = 5;
      std::cout << (ns % ref_ns); // 输出0，符合预期
      return 0;
  }
  

• 修正浮点数模运算的结果
  如果必须保留浮点数运算，可以判断模运算的结果是否接近ref_ns（比如差值小于一个极小的阈值，如1e-6），如果是则返回0，否则返回原结果：

  #include <iostream>
  #include <cmath>
  #include <algorithm>
  
  double safe_fmod(double value, double mod) {
      double result = fmod(value, mod);
      // 若结果接近mod，视为余数0；若结果接近负数，调整为正余数
      if (std::abs(result - mod) < 1e-6) {
          return 0.0;
      } else if (result < 0 && std::abs(result + mod) < 1e-6) {
          return 0.0;
      }
      return result;
  }
  
  int main() {
      double s = 30e-9;
      double ns = s * 1e9;
      double ref_ns = 5;
      std::cout << safe_fmod(ns, ref_ns); // 输出0
      return 0;
  }
  

• 调整输出精度（仅用于调试）
  如果只是不想让cout把4.9999997显示成5，可以设置更高的输出精度，查看真实的浮点数结果：

  #include <iostream>
  #include <iomanip>
  #include <cmath>
  
  int main() {
      double s = 30e-9;
      double ns = s * 1e9;
      double ref_ns = 5;
      std::cout << std::fixed << std::setprecision(7) << fmod(ns, ref_ns); // 输出4.9999997
      return 0;
  }
  

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