def _logpmf(self, x, n, p):
        k = floor(x)
        combiln = (gamln(n+1) - (gamln(k+1) + gamln(n-k+1)))
        return combiln + special.xlogy(k, p) + special.xlog1py(n-k, -p)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/binom_gen._logpmf/original