蒙特卡罗算法串行计算

蒙特卡罗算法是我在算法学习中最喜欢的算法之一，不知道什么原因，就是很欣赏，那种感觉就像你喜欢在雨中散步一样，没有什么能给你的，但是你享受到雨中的浪漫或者优雅。我一直觉得，像蒙特卡罗和拉斯维加斯这样的算法，在解决一些近视计算问题上，能提供给程序员很快捷和巧妙的方法。

下面给出用蒙特卡洛计算圆周率的程序，很简单，但是，要说明的是它的并行算法，看《蒙特卡罗算法并行计算》
#include
#include
#include
#include

using namespace std;

#define pear 50000//扔的“黄豆”数目
int isInCir(double x,double y)
{
if ((x*x+y*y) <= 1)
{
return 1;//在1/4园内的话，返1
}
return 0;
}

int main()
{
int i;
int num = 0;//在园内的豆豆数目
double x,y;
double mianji;
srand((unsigned)time(NULL));
for (i=0;i
{
x = rand()%RAND_MAX/float(RAND_MAX);
y = rand()%RAND_MAX/float(RAND_MAX);

if (isInCir(x,y) == 1)
{
num++;
}
}
cout<<"豆豆落在园内的个数："<<num<<endl;
cout<<"最大随机数："<<RAND_MAX<<endl;
mianji = (num*1.0)/pear;
//printf("Monte Carlo算法模拟出的半径为一的1/4圆的面积为%f\n",mianji);
cout<<"Monte Carlo算法模拟出的半径为一的1/4圆的面积为"<<mianji<<endl;
cout<<"计算出的pi值为："<<4*mianji<<endl;
return 0;
}