c++ - 存储用于蒙特卡罗模拟的 random_device 种子

Question

在 C++ 上，我使用 random_device 播种我的 64 位 mersenne twister，如下所示：

random_device rand_dev;
mt19937_64 mersenne_generator(rand_dev());

然后在循环内调用分布函数：

double rand_no = uniform_real_distribution<double>(0.0,1.0)(mersenne_generator);

换句话说，据我所知，使用 mt19937_64 的标准方式。我将它用于蒙特卡洛模拟，到目前为止效果很好。

但是，对于未来的模拟，我希望模拟结果能够基于存储的种子完全重现。这样运行：

mt19937_64 mersenne_generator(stored_rand_dev);

对于相同的参数集和给定的 stored_rand_dev 将产生相同的输出。

我的第一个想法是将 random_device 的结果存储为某种整数或字符串。但是，将 random_device 简单地类型转换为整数是行不通的，使用 to_string 也行不通。我想知道是否有另一种存储 random_device 的方法？最好采用可以占据 tsv 文件中单个单元格的形式。

或者，我想使用随机设备从 0 到 INT_MAX 选择一个无符号整数，并将其用作 mt19937 的种子。但是，我不确定这是否会影响 mt19937 与使用 random_device 生成的随机数的质量。任何人都可以就此提出建议吗？

标签： c++randomdevicemersenne-twister