def logChoose(n, c): return lgamma(n+1) - lgamma(c+1) - lgamma(n-c+1)
def logChoose(n, c): return lgamma(n + 1) - lgamma(c + 1) - lgamma(n - c + 1)