蒙特卡洛方法求圆周率

题目:

尝试使用蒙特卡罗法计算圆周率(π)的值。原理如下:

mengtekaluo.jpg

如图:边长为1的正方形内部有一个半径为1的四分之一个圆。
现往该正方形内随机投点,数量足够多的情况下,落入圆内的点m与落入整个
外切正方形的点n的数量比值大概为:  
m/n = (1/4)pi ,然后就可以得到π的值。
让点(x,y)投在整个矩形中,x与y的取值范围为(0≤x≤1, 0≤y≤1)
注意
1.请使用random库中的uniform来生成随机数。
2.使用运算符x*x求x的平方。
3.使用x**(1/2)求x的平方根。

输入格式:

在一行中给出种子和投点次数,中间以空格隔开。

输出格式:

对每一组输入,在一行中输出pi的值,并保留五位小数。

输入样例:

在这里给出一组输入。例如:

2 10000

输出样例:

在这里给出相应的输出。例如:

3.13560

 代码展示:

from random import *
a,b=map(int,input().split())
seed(a)
c=0
for i in range(b):
    x=uniform(0,1)
    y=uniform(0,1)
    ##取值范围是0到1
    if x**2+y**2<=1:
        c+=1
pi=float(4*(c/b))
print("%.5lf"%pi)

  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值