from numtools import fact # I didn't come up with this. fact(2 * 20) / (fact(20) ** 2)
def c(n, r): return fact(n) / (fact(r) * fact(n - r))