python10:圆周率的计算(蒙特卡罗法)

【问题描述】

求解圆周率可以采用蒙特卡罗方法,在一个正方形中撒点,根据在 1/4 圆内点的数量占总撒点数的比例计算圆周率值。

请以 123 作为随机数种子,获得用户输入的撒点数量,编写程序输出圆周率的值,保留小数点后 6 位。

【输入形式】
【输出形式】
【样例输入】

输入:"1024"

【样例输出】

输出:"3.218750"

【样例说明】
【评分标准】

扇形面积/正方形面积: ss/sz=pi/4
用扇形中点的个数模拟扇形的面积(正方形同理):ss=ps/n,sz=n
所以:pi=4*ss/sz=4*(ps/pz)

# ss/sz=pi/4
# 用扇形中点的个数模拟扇形的面积(正方形同理):ss=ps/n,sz=n
# pi=4*ss/sz=4*(ps/pz)

import math
import random

random.seed(123)

n = eval(input())  # 点的总个数
ps = 0  # 用于统计在扇形中的点个数

for i in range(n):
    # 随机点的坐标(x,y)
    x = random.random()  # 获得一个0~1之间的随机小数
    y = random.random()

    # 点到圆心的距离:d<=1,点在扇形区域内;
    d = math.sqrt(x * x + y * y)
    if d <= 1:
        ps += 1

pi = 4 * (ps / n)
print('%.6f'%pi)

 

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值