在C ++中设置离散分布

时间:2015-07-01 06:02:04

标签: c++ random distribution

经过几个小时的努力解决这个问题,我找不到任何有关我的错误的解释。

我希望计算机选择0到120(含)之间的随机数(加权)。我有一个数组,interval [],它保存从0到120(含)的数字。我有另一个数组,weight [],它保存了选择数组中第i个元素的概率。我想为这些数据定义概率分布函数。

这是我尝试过的。我得到一个错误,说没有构造函数的实例匹配参数列表。

我的代码片段:

std::vector< int> weights(121);
for (int i = 0; i < 121; i++)
{
    weights[i] = (teamData[i]).S();
}
discrete_distribution<> dist(weights.begin(), weights.end());

2 个答案:

答案 0 :(得分:1)

试试这段代码,看看你是否遇到同样的问题。

 // piecewise_constant_distribution
#include <iostream>
#include <array>
#include <random>

int main()
{
  const int nrolls = 10000; // number of experiments
  const int nstars = 100;   // maximum number of stars to distribute

  std::default_random_engine generator;
  std::array<double,6> intervals {0.0, 2.0, 4.0, 6.0, 8.0, 10.0};
  std::array<double,5> weights {2.0, 1.0, 2.0, 1.0, 2.0};
  std::piecewise_constant_distribution<double>
    distribution (intervals.begin(),intervals.end(),weights.begin());

  int p[10]={};

  for (int i=0; i<nrolls; ++i) {
    int number = distribution(generator);
    ++p[number];
  }

  std::cout << "a piecewise_constant_distribution:" << std::endl;
  for (int i=0; i<10; ++i) {
    std::cout << i << "-" << i+1 << ": ";
    std::cout << std::string(p[i]*nstars/nrolls,'*') << std::endl;
  }
  return 0;
}

答案 1 :(得分:1)

来自linked page(强调我的)

  

std::piecewise_constant_distribution生成随机浮点数   数字,均匀分布在几个数字中   子区间 [b i ,b i + 1 ,每个都有自己的权重 w i < / em>的。这套   区间边界和权重集是这个的参数   分布。

它期望浮点权重和边界,并且比边界重量少一个。它也将输出0到120之间的整数,但是浮点数。

您正在传递整数权重和边界,因此无法编译。但即使你修复了它,你仍然会从中获得浮点值......

相反,正如您发现的那样,您需要disrete_distribution,您可以这样设置:(从链接页面文档修改)

#include <iostream>
#include <map>
#include <random>

int main()
{
    // Setup the random bits
    std::random_device rd;
    std::mt19937 gen(rd());

    // Setup the weights (in this case linearly weighted)
    std::vector<int> weights;
    for(int i=0; i<120; ++i) {
        weights.push_back(i);
    }

    // Create the distribution with those weights
    std::discrete_distribution<> d(weights.begin(), weights.end());

    // use the distribution and print the results.
    std::map<int, int> m;
    for(int n=0; n<10000; ++n) {
        ++m[d(gen)/10];
    }
    for(auto p : m) {
        std::cout << p.first*10 << " - "<<p.first*10+9 << " generated " << p.second << " times\n";
    }
}