// //////////PREPROCESSOR////////// // #include #include #include #include // //////////CLASSES////////// // class Type { protected: std::default_random_engine generator; public: virtual void simulate() = 0; virtual void getResult() = 0; }; class Laplace : public Type { private: int *possibility; int it, n, size; public: Laplace(int s, int n) :size(s), n(n) { unsigned seed = static_cast(time(nullptr)); generator.seed(seed); possibility = new int[size]; for (it = 0; it < size; ++it) { possibility[it] = 0; } } ~Laplace() { delete[] possibility; } void simulate() { std::uniform_int_distribution distribution(0, size - 1); for (int it = n; it > 0; --it) { ++possibility[distribution(generator)]; } } void getResult() { std::cout << "Result: Hits: Out of: Average:" << std::endl; for (it = 0; it < size; ++it) { std::cout << std::setw(4) << it + 1 << std::setw(6) << possibility[it] << std::setw(7) << n << std::setw(12) << (double)n/(double)size << std::endl; } } }; class Bernoulli : public Type { private: int length, result = 0; double probability; public: Bernoulli(int n, double p) :length(n), probability(p) { unsigned seed = static_cast(time(nullptr)); generator.seed(seed); } void simulate() { std::bernoulli_distribution distribution(probability); for (int it = length; it > 0; --it) { if (distribution(generator)) ++result; } } void getResult() { std::cout << "Hits on result: " << result << std::endl; std::cout << "Other hits: " << length - result << std::endl; } }; // //////////MAIN FUNCTION////////// // int main() { Type *t; int choice; std::cout << "What kind of experiment do you want to simulate?\n1. Laplace\n2. Bernoulli" << std::endl; try { std::cin >> choice; switch (choice) { case 1: { int results, length; std::cout << "Please enter the number of possible results:" << std::endl; std::cin >> results; if (results <= 0) { throw "The experiment cant be simulated for that amount of results."; } else if (std::cin.bad() || std::cin.fail()) { throw "Something went wrong. Please try again!"; } std::cout << "How often should the experiment run?" << std::endl; std::cin >> length; if (length <= 0) { throw "The experiment cant be simulated for that amount of runs."; } else if (std::cin.bad() || std::cin.fail()) { throw "Something went wrong. Please try again!"; } Laplace l(results, length); t = &l; t->simulate(); std::cout << "Here are the results:" << std::endl; t->getResult(); } break; case 2: { int length; double probability; std::cout << "How often should the experiment run?" << std::endl; std::cin >> length; if (length <= 0) { throw "The experiment cant be simulated for that amount of runs."; } else if (std::cin.bad() || std::cin.fail()) { throw "Something went wrong. Please try again."; } std::cout << "What is the probability for the result to be hit?" << std::endl; std::cin >> probability; if (probability < 0.0 || probability > 1.0) { throw "Invalid probability."; } else if (std::cin.bad() || std::cin.fail()) { throw "Something went wrong. Please try again."; } Bernoulli b(length, probability); t = &b; t->simulate(); std::cout << "Here are the results:" << std::endl; t->getResult(); } break; default: { throw "Invalid mode. Please try again."; } } } //try catch (const char* msg) { std::cerr << "Error: " << msg << std::endl; } return 0; }

Random Experiment Simulator