概率算法 蒙特卡罗算法

/***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;
}