如何将伪随机ULong转换为0≤x<1的统计随机双精度浮点数？

你的核心问题出在两个地方：一是混淆了ulong数值的字节表示和double的位模式字节表示，二是对IEEE 754 double的位结构理解有误。

错误原因分析 #

你之前的代码用BitConverter.ToDouble解析resultUL的字节数组，这是把resultUL的数值转换成double，而不是把resultUL当成double的位模式来解析。另外，你构造的resultUL对应的double位模式中，指数位为0（偏移后指数为0，实际指数为-1023），所以得到的是极小的非规格化数，完全不符合预期。

正确实现方案 #

要生成[0,1)区间的double，我们需要利用IEEE 754 double的位结构：

• 符号位（最高位）：0（保证正数）
• 指数位：设置为1023（偏移后指数，对应实际指数0，构造出的数范围是[1.0, 2.0)）
• 尾数位：用xoshiro256**生成的ulong的高52位（double尾数只有52位，将ulong右移12位即可得到）

然后通过减法将[1.0,2.0)映射到[0.0,1.0)，这个减法是线性变换，不会损失随机性（仅受限于double的精度，无法保留原ulong的全部64位随机性，这是double的固有特性）。

正确代码如下：

public static double DoubleFromRand(ulong rand)
{
    // 构造double的位模式：符号位0，指数位1023（0x3FF前缀），尾数为rand的高52位
    ulong bitPattern = (rand >> 12) | 0x3FF0000000000000UL;
    // 将位模式转换为double，得到[1.0, 2.0)区间的数值
    double temp = BitConverter.Int64BitsToDouble((long)bitPattern);
    // 映射到[0.0, 1.0)区间
    return temp - 1.0;
}

验证 #

当传入UInt64.MaxValue时，rand >> 12得到52位全1的数值，构造的位模式为0x3FFFFFFFFFFFFFFF，对应的double是2 - 1/(2^52)，减去1后得到1 - 1/(2^52)，这是小于1且最接近1的double值，符合你的预期。

补充说明 #

如果你完全不想使用减法，唯一的办法是构造指数位为1022的位模式，但这样得到的数范围是[0.5,1.0)，无法覆盖整个[0,1)区间，因此减法是实现完整区间映射的必要步骤，且不会破坏随机性的均匀分布。

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