Question

我正在尝试使用boost在C ++ - 11中生成随机数。我已经看到一些关于这个的好帖子，我能够生成随机数，将我的所有代码放在main函数中，就像在这个网页（http://www.sitmo.com/article/generating-random-numbers-in-c-with-boost/）中所建议的那样。

但是，我无法封装可由不同功能调用的随机数生成器（RNG）。我希望这个封装的RNG接收数字范围和随机种子作为参数。我尝试按照这篇文章的建议（Boost random number generator）来做到这一点。

我实现了以下代码。但两者都失败了：随机种子的变化和价值范围的定义。

#include <boost/random/variate_generator.hpp>
#include <boost/generator_iterator.hpp>
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/normal_distribution.hpp>

typedef boost::mt19937                     ENG;    // Mersenne Twister
typedef boost::normal_distribution<double> DIST;   // Normal Distribution
typedef boost::variate_generator<ENG,DIST> GEN;    // Variate generator

class RNG {

  private:
     ENG eng;
     DIST dist;

  public:
     GEN gen;
     RNG(double min,double max,int seed)
     : dist(min,max), gen(eng,dist)
     {eng.seed(seed); }
};

这是我的主要功能。

int main(int argc, char **argv){
    if(argc == 2){
        int myseed = atoi(argv[1]);
        double minimum = 1.;
        double maximum = 10.;
        int N = 20;

        RNG myrn(minimum,maximum,myseed);

        for (int i = 0; i < N; i++){
            std::cout << myrn.gen() << std::endl;
        }
            return 0;
    }
    exit(1);
}

此实现独立于我使用的种子生成相同的随机数序列。此外，随机数不在最小值和最大值之间的范围内。你们有谁知道如何解决这两个问题吗？

Answer 1

boost::variate_generator按值获取其构造函数参数，因此当您在mem-initializer列表中传递默认构造的mt19937时，它会复制它。在构造函数体内播种引擎对该副本没有影响。要解决此问题，请将模板参数类型更改为引用。

typedef boost::variate_generator<ENG&,DIST> GEN;
//                                  ^

接下来，我认为你不想要normal_distribution，它的构造函数参数不是分布产生的值范围的最小值和最大值，而是分布的均值和标准差。

您可能正在寻找的分布是uniform_real_distribution。

typedef boost::random::uniform_real_distribution<> DIST;

完成这些修复后，您的代码应该按照您的意愿运行。

Live demo

使用C ++ 11编译器甚至不需要Boost来执行此操作，<random>标头提供了您需要的所有组件。唯一显着的区别是缺少variate_generator类，而是在您想要生成随机数时使用引擎作为参数调用分布对象。

#include <random>

typedef std::mt19937                     ENG;    // Mersenne Twister
typedef std::uniform_real_distribution<> DIST;   // Uniform Distribution

class RNG {

  private:
     ENG eng;
     DIST dist;

  public:
     DIST::result_type gen() { return dist(eng); }
     RNG(double min,double max,int seed)
     : dist(min,max)
     {eng.seed(seed); }
};

Live demo

封装的随机数生成器在C ++ - 11中使用boost

1 个答案: