多精度浮点`scipy.misc.comb`

Question

我正在使用mpmath编写一个适合计算极小概率的二项式beta函数。

我需要的大多数函数都已经包含在mpmath中，或者很容易重写。然而，不幸的是，scipy.misc.comb并非如此。我试着看了源代码，但它似乎是based on从a file导入的某个binom函数，在其中我找不到它的代码。

如何在不使用for/while循环计算阶乘的情况下重写scipy comb函数(或编写自己的函数)？

Answer 1

好吧，显然mpmath确实有一个阶乘函数，它的名字看起来很直观，mpmath.factorial -。

因此，我的查询的解决方案很简单：

def mp_comb(N,k):
    val= mpmath.factorial(N)/(mpmath.factorial(k)*mpmath.factorial(N-k))
    return val