首页 > 代码库 > 一种计算e的方法

一种计算e的方法

原理:平均e个(0,1)之间的随机数之和会大于1.
原因:n个数之和小于1的概率是1/n!,则n个数之和大于1的概率则是1-1/n!;
恰好n个数之后大于1的概率,等于n个数之后大于1的概率减去n-1个数之和大于1的概率,即
(1-1/n!)
- (1-1/(n-1)!) = (n-1)/n!
则n的期望为
sum( (n-1)/n! * n ) = sum( (n-1)/(n-1)! )
=
sum(n/n!)
又因为sum(n/n!)=e,所以平均e个(0,1)之间随机数之和大于1.
代码如下:可惜收敛较慢,效率较低。
#!/usr/bin/python
‘‘‘A
method to calculate e
‘‘‘
__author__ =
"hankjin<hankjin@163.com>"
import random,math
def
calce(iter=1000000):
    ‘‘‘ calculate math.e with <iter>
iterations.
    ‘‘‘
    res=0
    for i in range(iter):
       
tmp=0.0
        while tmp<1:
            res+=1
           
tmp+=random.random()
    return float(res)/iter
if
__name__==‘__main__‘:
    print(‘real e=‘, math.e)
    print(‘calculated
e=%s‘, calce())