C++梅森旋轉演算法生成亂數(mersenne_twister_engine)詳解

2020-07-16 10:04:29
mersenne_twister_engine 類別範本實現了梅森選擇演算法,它被這樣叫是因為週期長度是一個梅森素數。

梅森素數是 2n-1 形式的素數,因此 7 和 127 是梅森素數;當然,用在這個演算法中的梅森素數更大。這個引擎的應用非常廣泛,因為它可以生成非常高品質的序列,但存在速度相對較慢的缺點。這個演算法很複雜並且包含很多的引數,因此在此不作解釋。

梅森素數演算法範例

對於定義具體生成器的 mersenne_twister_engine 的範例有兩個型別別名。mtl9937 生成隨機的無符號 32 位整數,mtl9937_64 生成無符號的 64 位整數。mtl9937 亂數生成器的週期長度為 219937-1,它因此得名。

可以用和其他生成相同的方式使用它們:
std::random_device rd;
std::mtl9937_64 rng {rd()}; // Generates random 64-bit integers
std::uniform_real_distribution<long double> dist {-5.OL, 5.OL};
for(size_t i {}; i < 8; ++i)
    std::cout << std::setw(5)<< dist(rng)<<" "; //-2.57481 3.0546 -1.6438 2.14798 -3.84095 0.973843 -2.98971 -2.1067