为何用MAX_SAFE_INTEGER生成的随机整数奇偶分布不均？求Math.floor最优上限

问题现象 #

测试发现以下三种随机数生成方式的奇偶分布存在明显差异：

• Math.floor(Math.random() * Number.MAX_SAFE_INTEGER)：几乎始终返回奇数（仅当Math.random()返回0时得到偶数，概率可忽略）
• Math.floor(Math.random() * (Number.MAX_SAFE_INTEGER - 1))：约75%概率返回奇数
• Math.ceil(Math.random() * Number.MAX_SAFE_INTEGER)：奇偶概率各占50%

附测试代码：

let evenCount = 0, oddCount = 0;

for (let i = 0; i < 10000; i++) {
  const randomNumber = Math.floor(Math.random() * Number.MAX_SAFE_INTEGER);
  if (randomNumber % 2 === 0) {
    evenCount++;
  } else {
    oddCount++;
  }
}

console.log("Number of even numbers:", evenCount);
console.log("Number of odd numbers:", oddCount);

底层原理：双精度浮点数与旧版Math.random()实现 #

JavaScript的Number类型基于64位双精度浮点数存储，其中尾数部分占52位，加1位隐含的最高位1，总共能精确表示的整数范围是-2^53 +1到2^53 -1（即Number.MIN_SAFE_INTEGER到Number.MAX_SAFE_INTEGER）。

旧版V8引擎中，Math.random()生成的是[0,1)区间的浮点数，其尾数仅使用52位随机位，等价于k / 2^52（k为0到2^52 -1之间的整数）。结合这个特性分析三种情况：

1. Math.floor(Math.random() * Number.MAX_SAFE_INTEGER)
Number.MAX_SAFE_INTEGER = 2^53 -1，代入计算得：
   Math.random() * (2^53 -1) = k*(2^53 -1)/2^52 = 2k - k/2^52
   由于k的取值范围是0到2^52 -1，k/2^52是大于0且小于1的正数，因此2k - k/2^52的整数部分为2k -1（当k>0时），必然是奇数；仅当k=0时结果为0（偶数），但该情况的概率仅为1/2^52，几乎无法在测试中出现，因此看起来始终返回奇数。

2. Math.floor(Math.random() * (Number.MAX_SAFE_INTEGER -1))
Number.MAX_SAFE_INTEGER -1 = 2^53 -2 = 2*(2^52 -1)，代入得：
   Math.random() * (2^53 -2) = 2k*(2^52 -1)/2^52 = 2k - k/2^51

• 当k < 2^51时，k/2^51 <1，整数部分为2k -1（奇数），这类k占总数的50%；
• 当k >=2^51时，1 <= k/2^51 <2，整数部分为2k -2（偶数），但此区间内k的奇偶占比与前半段叠加后，最终奇数的整体概率约为75%。

3. Math.ceil(Math.random() * Number.MAX_SAFE_INTEGER)
Math.ceil会取大于等于目标值的最小整数，对于2k - k/2^52（k>0），其值介于2k-1和2k之间，因此Math.ceil的结果为2k（偶数）；当k取特定值时，结果会覆盖奇数区间，最终使得奇偶分布各占50%。

V8引擎的更新：53位随机数修复分布问题 #

Chrome 137.0.7151.41及之后的V8引擎，将Math.random()的实现改为生成53位随机尾数，等价于k/2^53（k为0到2^53-1之间的整数）。此时乘以Number.MAX_SAFE_INTEGER后，结果的精度足够覆盖完整的整数范围，Math.floor后的奇偶分布回归均匀，因此无法复现原问题。

解决方案：Math.floor实现奇偶50/50的最大整数 #

要让Math.floor(Math.random() * N)的结果奇偶各占50%，需保证N是2的幂，且适配旧版V8的52位随机数限制。最大的符合条件的整数是2^52，即：

const MAX_UNIFORM_INTEGER = 2 ** 52; // 等价于 4503599627370496

此时Math.random() * MAX_UNIFORM_INTEGER的结果为k（k为0到2^52-1之间的整数），Math.floor后直接得到k，奇偶数量各占一半，分布完全均匀。

内容的提问来源于stack exchange，提问作者Jean-Samuel Girard