二项式分布的随机数

Question

我需要从二项式分布中快速生成大量随机数，以适应截然不同的试验规模（但是，大多数试验规模很小）。我希望不必手动编写算法（参见，例如 11 月份的相关讨论），因为我是一名新手程序员，不喜欢重新发明轮子。 Boost 似乎不提供二项分布变量的生成器，但 TR1 和 GSL 即可。是否有充分的理由选择其中之一，或者我根据自己的情况写一些定制的东西更好？我不知道这是否有意义，但我将在整个程序中交替从均匀分布和二项式分布生成数字，并且我希望它们共享相同的种子并最大限度地减少开销。我希望得到一些关于我应该考虑的建议或例子。

Answer 1

Boost 1.43 似乎支持 二项分布。您可以使用 boost::variate_generator 将随机源连接到类型
您想要从中采样的分布。

所以你的代码可能看起来像这样（免责声明：未经测试！）：

boost::mt19937 rng;                 // produces randomness out of thin air
                                    // see pseudo-random number generators
const int n = 20;
const double p = 0.5;
boost::binomial<> my_binomial(n,p);      // binomial distribution with n=20, p=0.5
                                         // see random number distributions
boost::variate_generator<boost::mt19937&, boost::binomial<> >
         next_value(rng, my_binomial);     // glues randomness with mapping
int x = next_value();                      // simulate flipping a fair coin 20 times

Answer 2

您误解了 Boost 模型 - 您选择随机数生成器类型，然后选择一个分布，将 RNG 生成的值分散到该分布上。 这个答案中有一个非常简单的例子，它使用均匀分布，但其他分布使用相同的基本模式 - 生成器和分布完全解耦。