R语言计算2.0^64.0-1.0输出异常，是否超出最大数值？

问题重现 #

运行以下R代码：

> sprintf("%f",2.0^64.0-1.0)
[1] "18446744073709551616.000000"

实际预期结果应为 18446744073709551615，这个偏差不是因为超出R语言的最大数值，而是双精度浮点数的精度限制导致的。

核心原因 #

• R默认的数值类型是双精度浮点数（double），遵循IEEE 754标准，这种类型的有效精度只有53位二进制位。
• 2^64-1 是一个64位全1的整数，已经超出了双精度浮点数能精确表示的整数范围（双精度最多能精确表示2^53以内的整数）。当计算这个值时，它会被舍入到最近的可精确表示的浮点数，也就是2^64，最终输出结果就变成了18446744073709551616。

解决方法 #

如果需要精确处理这类大整数，可以用以下两种方式：

• 使用bit64包的64位整数类型：

  library(bit64)
  # 直接构造目标大整数
  as.integer64("18446744073709551615")
  

• 使用gmp包进行大整数运算：

  library(gmp)
  # 精确计算2^64 - 1
  as.bigz(2)^64 - 1
  

内容的提问来源于stack exchange，提问作者Jean-Luc Dupouey