蒙特卡洛法求圆周率100亿数据
代码[kou@hodge python]$ cat circle.py#!/usr/bin/env pythoncoding=utf-8import timeimport randomhits=0pi=0DARTS=100000*100000start=time.perf_counter()for i in range(DARTS):x,y=random.random...
·
代码
import time
import random
hits=0
pi=0
DARTS=100000*100000
start=time.perf_counter()
for i in range(DARTS):
x,y=random.random(),random.random()
dist=pow(x ** 2+y**2,0.5)
if dist <= 1.0:
hits+=1
pi=4*(hits/DARTS)
print("圆周率的值是{:.10f}".format(pi))
print("程序运行时间为{}s".format(time.perf_counter()-start))
100万个数据
100亿个数据
三个半小时运行时间
更多推荐
已为社区贡献3条内容
所有评论(0)