R语言中factorial()计算52阶乘结果异常，手动计算与函数结果不符

你没操作出错，问题出在R默认数值类型的精度限制上：

• R里默认用双精度浮点数存储数值，这种类型只有53位有效数字的精度，当数值超过2^53（约9e15）后，无法精确表示所有整数。
• 手动连乘1*2*…*52时，从21!开始，乘积就超过了2^53的范围，后续每一步乘法都会丢失精度，最终结果偏离准确值。
• 而factorial()函数内部不是直接做逐次整数乘法，它调用伽马函数gamma(n+1)计算，这个函数针对大数值做了精度优化，能返回更接近准确值的结果（但本质还是浮点数近似，精度损失更小）。

验证示例 #

# 查看双精度浮点数能精确表示的最大整数
2^53
# [1] 9007199254740992

# 20!小于2^53，手动连乘和factorial结果一致
factorial(20) == 1*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19*20
# [1] TRUE

# 21!超过2^53，手动连乘开始出现误差
factorial(21) == 1*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19*20*21
# [1] FALSE

解决方法：用任意精度整数计算 #

如果需要精确的大整数阶乘，可使用gmp包的factorialZ()函数，它支持任意精度的整数运算：

library(gmp)
# 计算精确的52!
exact_fact52 <- factorialZ(52)
exact_fact52
# 输出准确的52阶乘整数结果

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