Многоточность с плавающей точкой `scipy.misc.comb`

Question

Я использую mpmath написать биномиальную бета-функцию, подходящую для вычисления крайне малых вероятностей.

Большинство функций, которые мне нужны, уже включены в mpmath или легко переписываются. К сожалению, однако, scipy.misc.comb не является. Я пытался смотреть на исходный код, но, похоже, он основан на некоторых binom Функция импортирована из файла, в котором я не могу найти его код.

Как я мог бы переписать функцию гребенчатого гребня (или написать свою собственную), не используя циклы for/while для вычисления факториалов?

Answer 1

Ок, видимо mpmath имеет факториальную функцию, с обманчиво интуитивным названием, mpmath.factorial -.-

Таким образом, решение моего запроса просто:

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