monte.cpp

蒙特卡罗求圆周率算法

#include <stdlib.h>
#include <iostream>
using namespace std;
void main()
{
	double m = 0;

	for(int i=0;i<100000;i++)
	{
		double ans1=(double)rand()/RAND_MAX; 
		double ans2=(double)rand()/RAND_MAX; 
		if((ans1*ans1+ans2*ans2)<1)m++;		
	}
	cout<<m/100000*4<<endl;
}