1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
|
// 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;
}
|