R：大数字错误的算术

Question

当R执行简单算术的大数字时，我在R中遇到了一个问题。有人可以解释发生的事情以及如何工作吗？

> a <- 51569
> b <- a + 6924851946518374722
> b
[1] 6924851946518425600 # problem is visible here already
> c <- b - 6924851946518374722
> c
[1] 51200

因此，添加大数字时的错误为469。为什么？

Answer 1

R使用双精度浮点来表示值。精度取决于您的计算机。在我的Intel Machine上，带有IEEE双打，我获得了52位精确度：

> options(digits=22)
> (0:1) + 2^52
[1] 4503599627370496 4503599627370497
> (0:1) + 2^53
[1] 9007199254740992 9007199254740992

Answer 2

r有一个包装GMP的软件包，GNU多个精度库：

cran：ran> cran：r GMP手册

library(gmp)
a <- as.bigz("51569")
b <- a + as.bigz("6924851946518374722")
c <- b - as.bigz("6924851946518374722")
c

Big Integer ('bigz') :
[1] 51569

正如评论中指出的那样，报价是必要的。

as.bigz(6924851946518374722) # Does not work as intended
Big Integer ('bigz') :
[1] 6924851946518374400

因为否则r将参数转换为双重，在转换为bigz之前丢失精度。

还请注意，R尚未支持64位整数，例如

[1] 4294967296L
Warning message:
non-integer value 4294967296 qualified with L; using numeric value