蒙特卡羅方法
利用隨機(jī)點(diǎn)分布，讓隨機(jī)點(diǎn)足夠的多，用指定區(qū)域隨機(jī)點(diǎn)數(shù)與總區(qū)域隨機(jī)點(diǎn)數(shù)做比值來求取指定區(qū)域面積

求解圓周率，在一個(gè)正方形中撒點(diǎn)，根據(jù)在1/4圓內(nèi)點(diǎn)的數(shù)量占總?cè)鳇c(diǎn)數(shù)的比例計(jì)算圓周率值。

代碼如下：

#CalPi.py from random import random from time import perf_counter DARTS = 1000*1000 #當(dāng)前在區(qū)域中點(diǎn)的總數(shù)量 hits = 0.0 #目前在圓的內(nèi)部的點(diǎn)的數(shù)量 start = perf_counter() for i in range (1,DARTS+1):x,y = random(),random()dist = pow(x**2+y**2,0.5)if dist <= 1.0:hits = hits+1 pi = 4*(hits/DARTS) print("圓周率的值是：{}".format(pi)) print("運(yùn)行時(shí)間是：{:.5f}s".format(perf_counter()-start))

現(xiàn)以123作為隨機(jī)數(shù)種子，獲得用戶輸入的撒點(diǎn)數(shù)量，編寫程序輸出圓周率的值，保留小數(shù)點(diǎn)后6位。

代碼如下：

from random import random,seed seed(123) darts = eval(input()) hits = 0.0 for i in range (darts):x,y = random(),random()d = pow(x**2+y**2,0.5)if d <= 1.0:hits = hits+1 pi = 4*(hits/darts) print("圓周率的值為:{:.6f}".format(pi))

總結(jié)

以上是生活随笔為你收集整理的【Python】蒙特卡罗方法计算圆周率及给定随机数种子的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。