正态分布随机数示例

#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;
}