计算圆周率近似值

利用圆与外切正方形的面积之比来近似计算 PI 的值

扇形面积在正方形中所占的比例为 K

而 K 值可通过随机撒点 所占的比模拟近似求得

因为 rand() 函数产生的随机数范围是 0~ 32767 ,RAND_MAX 代表最大的随机数 故
下面有 0 <= X <= 1; 0<= Y <= 1;
当xx + xx <=1时 落在扇形里;

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<time.h>

//生成0~1 之间的随机数; 
double frandom()
{
	//srand(time(NULL)); 
	return (double) rand() / RAND_MAX; 
}
int main()
{
	double x,y,p;
	int i,count,n;
	
	srand((unsigned int)time(NULL)) ; //设置种子
	 
	scanf("%d",&n);
	
	//模拟随机投点n 次;
	count =0; 
	 
	for(int i=0;i<n;i++)
	{
		x= frandom();
		y= frandom();
	//	printf("%f   %f\n",x,y);
		if(sqrt(x*x + y*y)<=1) // 若点位于 扇形内;
			count ++;
	}
	p= 4.0 * count / n; //计算 PI 的近似值 
	printf("%.4f\n",p); //输出 PI 的近似值 
	
	return 0;
 } 

srand() 函数可产生 0~32767 的随机数 srand() 用来设置 rand() 产生随机数种子
在调用 rand()函数前 如果没有用 srand() 函数设置随机数种子,随机数种子 将自动设成相同值 1;
将会导致 rand ()所产生的随机数都一样,所以在调用 rand()函数时 先利用 srand()设好随机数种子。
通常用 time(NULL) 的返回值做种子 ,time(NULL)的返回值是 国际标准时间公元1970 年1月1日零点 以来经过的秒数。

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值