如何生成[std::numeric_limits<double>::lowest(), max()]范围的均匀随机数

直接使用std::uniform_real_distribution<double>构造[std::numeric_limits<double>::lowest(), std::numeric_limits<double>::max()]的分布不可行，因为该分布要求上限与下限的差值不能超过numeric_limits<double>::max()，而这个区间的差值是2 * max()，明显超出限制。

根据你对“均匀”的不同需求，有两种实现方案：

方案1：每个有限double值等概率（位层面均匀） #

这种方案生成的每个合法有限double值（包括非规格化数）概率完全相同，适合需要遍历所有可能double数值的场景：

#include <random>
#include <cstdint>
#include <limits>

double generate_bit_uniform_double() {
    static std::mt19937_64 rng(std::random_device{}());
    
    // 生成符号位：0表示正数，1表示负数，各50%概率
    std::uniform_int_distribution<uint64_t> sign_dist(0, 1);
    uint64_t sign = sign_dist(rng) << 63;
    
    // 生成指数位：避开全1的情况（对应无穷/NaN），范围0~2046
    std::uniform_int_distribution<uint64_t> exp_dist(0, 0x7FE);
    uint64_t exponent = exp_dist(rng) << 52;
    
    // 生成尾数位：覆盖所有52位可能的取值
    std::uniform_int_distribution<uint64_t> mantissa_dist(0, 0xFFFFFFFFFFFFF);
    uint64_t mantissa = mantissa_dist(rng);
    
    // 组合所有位并转换为double
    uint64_t bits = sign | exponent | mantissa;
    return *reinterpret_cast<double*>(&bits);
}

注意事项： #

• 该方法排除了无穷大、NaN等非法值，只生成有限double值
• 数值分布并非线性均匀：靠近0的数值密度远高于±max()附近的数值，这是由double的存储特性决定的（小数值精度更高）

方案2：线性均匀分布（数值区间均匀） #

如果你需要的是数值范围上的线性均匀（即任意长度相同的子区间被选中的概率相同），可以通过缩放[0,2)区间的均匀数来实现：

#include <random>
#include <limits>

double generate_linear_uniform_double() {
    static std::mt19937 rng(std::random_device{}());
    static std::uniform_real_distribution<double> dist(0.0, 2.0);
    
    // 将[0,2)映射到(-max(), max())
    return (dist(rng) - 1.0) * std::numeric_limits<double>::max();
}

特点： #

• 生成的数值在(-max(), max())范围内线性均匀分布，符合你最初想用uniform_real_distribution达到的效果
• 理论上不包含端点max()和lowest()，但由于这两个值的概率为0，实际使用中可以忽略；若必须包含，可额外以极低概率返回这两个值

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