开张纪念

有人说阶乘计算很慢，确实，如果按照定义去计算，那么多次乘法是很慢的。

但是事实上我们有现成的工具可用，比如伽玛函数，斯特林公式。

其中gmp库中就有现成的gamma函数，所以大数的阶乘并不慢。

import gmpy2
from gmpy2 import mpz
def stirling(z):
        return gmpy2.gamma(z+1)
import time
t=time.clock()
print(stirling(mpz(30000000)))
t=time.clock()-t
print(t)

　　结果：

2.0996305395095615e+211284807
0.000302

速度很快，但是更大的数就不行了，因为受制于gmp浮点数的上限