Priming the Mersenne twister PRNG

Priming the Mersenne twister PRNG

本文关键字:PRNG twister Mersenne the Priming      更新时间:2023-10-16

似乎围绕MT19937的使用有一些神话,特别是,一旦播种了'某些'一定数量的'一定数量的'生成器生成的位置,应忽略以至于只有尽可能接近伪随机性。

我看到的代码示例如下:

boost::mt19937::result_type seed = 1234567; //taken from some entropy pool etc
boost::mt19937 prng(seed);
boost::uniform_int<unsigned int> dist(0,1000);
boost::variate_generator<boost::mt19937&,boost::uniform_int<unsigned int> > generator(prng,dist);
unsigned int skip = 10000;
while (skip--)
{
   generator();
}
//now begin using for real.
....

我的问题是:

  1. 这是神话还是有些事实?

  2. 如果它可行,应该忽略多少位?作为我看到的数字
    似乎是任意的

在第一个评论中引用的论文,默瑟恩·扭动的初始化,不仅是某人,而且他是该论文的两个合着者之一。

根据Boost文档,使用单个32位整数(4个字节)作为该发电机的种子的问题是,发电机的内部状态为2496字节。这么小的种子需要一段时间才能传播到发电机的其余状态,尤其是因为Twister并不意味着在密码上安全。

,这并不奇怪。

要解决您对需要运行发电机一段时间以开始的担忧,您需要替代(和显式)构造函数。

template<typename SeedSeq> explicit mersenne_twister_engine(SeedSeq &);

这是第三条评论的精神,您可以用比单个整数更长的东西初始化。序列提供的序列来自某些发电机。要使用熵池,请将发电机从熵池写成适配器,根据需要从池中返回值。