概率算法 蒙特卡罗算法

/***Date:2014.12.11***/

//概率算法思想：统计学思路。

//基本过程：1）将问题转化为应的容易计算面积的几何图形S，问题结果对应几何图形中的某一部分S1；

////////////2)向几何图形中随机撒点；

////////////3）统计几何图形S、S1中的点数，根据二者面积关系以及二者中的点数来计算得到结果；

////////////4）判断结果是否达到需要精度结果，若符合则输出，否则继续（2）步骤。

//蒙特卡罗（Monte Carlo）算法

//均匀撒点：利用随机函数实现，产生(0,1)之间随机坐标值[x,y];

//区域判断：通过计算 x^2 + y^2 <=1 实现。

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#include <time.h>

double MontePI(int n)

{

double PI,x,y;

int i,sum;

sum = 0;

srand((unsigned)time(NULL));//64位要强转一下,要不然虽可执行，但不安全，PS，我电脑VS都64位

for (i = 1;i < n;i ++)

{

x=(double)rand()/RAND_MAX;

y=(double)rand()/RAND_MAX;

if (x*x + y*y <= 1)

{

sum ++;

}

}

PI = 4.0 * sum /n;

return PI;

}

int main()

{

int n;

double PI;

printf("蒙特卡罗概率算法计算 π：\n");

printf("输入点的数量(数量越多，结果越接近真实值)：");

scanf_s("%d",&n);

PI = MontePI(n);

printf("π =  %f\n",PI);

system("pause");

return 0;

}