【问题描述】
求解圆周率可以采用蒙特卡罗方法,在一个正方形中撒点,根据在 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)