pub fn binomial_coefficient(n: usize, k: usize) -> u64
n choose k, computed multiplicatively so it stays exact for the small values derivative formulas need.
n choose k