[英]Setting up a Discrete Distribution in C++
经过数小时的努力解决此问题后,我找不到我的错误的任何解释。
我希望计算机选择0到120(含)之间的随机数(加权)。 我有一个数组interval [],其中保存着0到120(含)之间的数字。 我有另一个数组weights [],其中包含在数组中选择第ith个元素的概率。 我想为这些数据定义一个概率分布函数。
这是我尝试过的。 我收到一条错误消息,指出没有构造函数实例与参数列表匹配。
我的代码段:
std::vector< int> weights(121);
for (int i = 0; i < 121; i++)
{
weights[i] = (teamData[i]).S();
}
discrete_distribution<> dist(weights.begin(), weights.end());
尝试这段代码,看看是否遇到相同的问题。
// 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;
}
从您的链接页面 (重点是我的)
std::piecewise_constant_distribution
生成随机浮点数 ,这些浮点数均匀分布在几个子间隔[b i ,b i + 1 ]的每个子间隔中,每个子间隔都有自己的权重w i 。 间隔边界集和权重集是此分布的参数。
它期望浮点权重和边界,并且权重比边界小一。 它也不会输出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";
}
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.