PI-撒币算法

首先构造一个单位正方形和一个四分之一圆，

然后假设你有一堆硬币

你开始随机对上述构造的正方形，撒币！！！！！

当然这个硬币可能在圆里，也可能在圆外

只要你的硬币够多，那么你的硬币将构成1/4圆。

通过计数其中落入内切圆的硬币的个数，有：

如果一共投入a个硬币，其中有b个落入圆中，则只要硬币均匀，假定圆周的半径为r，则：

你撒币撒得越多，那么pi就越准确。

所以以上算法就是蒙特拉罗算法。

代码：

#random返回随机生成的一个实数，它在[0,1)范围内
from random import random
import math

sabi = 1200000
count = 0
for i in range(1,sabi):
    x = random()
    y = random()
    distance = math.sqrt(pow(x,2)+pow(y,2))
    if distance <= 1.0:
        count += 1
    pi = 4*(count/sabi)
print(pi)

结果：

3.1420666666666666