随机算法之求π

最新推荐文章于 2024-04-07 19:31:14 发布

weixin_34217711

最新推荐文章于 2024-04-07 19:31:14 发布

阅读量1k

点赞数 1

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/520xiuge/p/6377871.html

用随机方法求π,具体方法如下:
易想到用 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\n",4.0*count/n);
	}
	return 0;
}

参照:http://www.acmerblog.com/random-%cf%80-2512.html

转载于:https://www.cnblogs.com/520xiuge/p/6377871.html