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

Question

我正在使用mpmath来编写一个二项式β函数，以计算极小的概率。

我需要的大多数函数都已经包含在 mpmath 中或易于重写。然而，可悲的是，scipy.misc.comb并非如此。我尝试查看源代码，但它似乎基于从我找不到其代码的文件中导入的一些binom函数。

如何在不使用 for/while 循环来计算阶乘的情况下重写 scipy 梳子函数（或编写我自己的函数）？

Answer 1

好的，所以显然mpmath确实有一个阶乘函数，具有看似直观的名称，mpmath.factorial -.-

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

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