随机算法之求π
用随机方法求π,具体方法如下:
易想到用 S(正方形内切圆 )与S(正方形）比值求得
我们首先对内切圆圆心为坐标原点建立坐标系, 考虑到随机函数rand()
产生0~RAND_MAX的随机数,所以产生的点均落在第一象限,
第一象限1/4圆和圆外 的点的个数比值即为 S(正方形内切圆 )与S(正方形）比值
#include <stdio.h>
#include <stdlib.h> //srand,rand,RAND_MAX
#include <time.h> //time
#include <math.h>
int main(){
//比较10次运算结果
for(int i = 0 ; i < 10 ; i ++ ) {
double x,y;
int n,count;
srand((unsigned int)time(NULL));//初始化伪随机数生成器
printf("产生随机数个数(尽量较大):");
scanf("%d",&n);
count = 0;
for(int j = 0 ; j < n ; j ++ ){
x = rand()/(double)RAND_MAX;//0.0~1.0之间
y = rand()/(double)RAND_MAX;
if(sqrt(1-x*x) >= y) //圆内
count ++;
}
printf("pai:%5.8f
",4.0*count/n);
}
return 0;
}
参照:http://www.acmerblog.com/random-%cf%80-2512.html