fn lchoose(n: f64, k: f64) -> f64
Log of the binomial coefficient: lchoose(n, k) = lgamma(n+1) - lgamma(k+1) - lgamma(n-k+1).