蒙特卡羅演算法串列計算

蒙特卡羅演算法是我在演算法學習中最喜歡的演算法之一，不知道什麼原因，就是很欣賞，那種感覺就像你喜歡在雨中散步一樣，沒有什麼能給你的，但是你享受到雨中的浪漫或者優雅。我一直覺得，像蒙特卡羅和拉斯維加斯這樣的演算法，在解決一些近視計算問題上，能提供給程式員很快捷和巧妙的方法。

下面給出用蒙特卡洛計算圓周率的程式，很簡單，但是，要說明的是它的並行演算法，看《蒙特卡羅演算法並行計算》
#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;
}