题目：小明向他的女朋友仙仙求婚，在求婚戒指上刻了一个大大的爱心。仙仙看到爱心想考验一下小明，出了一道题。方程 (x2+y2−1)2−x2y2=0 能画出一个美丽的爱心，现在给定一个点 (x,y) ，其中， X∼N((μ1,σ21) ， Y∼N((μ2,σ22) ，这个点在这个爱心里面的概率是多少，精确到0.1（小数点后保留一位小数，比如0.1, 0.2等）
#include <iostream> #include <math.h> using namespace std; /*生成标准正态分布随机数*/ double gaussRand(){ static double V1,V2,S; static int phase=0; double X; if(phase==0){ do{ double U1=(double)rand()/RAND_MAX; double U2=(double)rand()/RAND_MAX; V1=2*U1-1; V2=2*U2-1; S=V1*V1+V2*V2; }while(S>=1||S==0); X=V1*sqrt(-2*log(S)/S); } else X=V2*sqrt(-2*log(S)/S); phase=1-phase; return X; } /*将随机值从正态分布N(0,1)转化到高斯分布N(mu,sigma)*/ double gaussValue(double val,double mu,double sigma) { return mu+sigma*val; } /*判断随机值是否在曲线内*/ bool isInside(double x,double y) { double x2=x*x; double y2=y*y; double temp=x2+y2-1; return (temp*temp-x2*y2)>0?false:true; } /*模拟*/ double learnCurve(double mu1,double sigma1,double mu2,double sigma2) { int n=10000;//模拟次数 int count=0; int i=n; while(i-->0){ double x=gaussValue(gaussRand(),mu1,sigma1); double y=gaussValue(gaussRand(),mu2,sigma2); if(isInside(x,y)) count++; } return (double) count/n; } int main() { double mu1,mu2,sigma1,sigma2; cin>>mu1>>sigma1>>mu2>>sigma2; double res=learnCurve(mu1,sigma1,mu2,sigma2); printf("%.1f\n",res); return 0; }参考：
http://blog.csdn.net/shijing_0214/article/details/60465542 http://www.cnblogs.com/tsingke/p/6194737.html
