正态分布随机数示例

#include"stdafx.h"#include<random>const double eps = 1e-6;const int N = 1000000;int main(){double mu1 = 0.0, sigma1 = 1.0, mu2 = 0.0, sigma2 = 1.0;std::default_random_engine e;std::normal_distribution<double> x(mu1, sigma1);std::normal_distribution<double> y(mu2, sigma2);int ans = 0;for (int i = 0; i < N; i++){double x1 = x(e);double y1 = y(e);double heart = (x1*x1 + y1*y1 - 1)*(x1*x1 + y1*y1 - 1)*(x1*x1 + y1*y1 - 1) - x1*x1*y1*y1*y1;if (heart < 0.0)ans++;}printf("%.1lf\n", ans*1.0 / N);system("pause");return 0;}