pub fn checked_multinomial(n: u64, ni: &[u64]) -> Result<f64>
Computes the multinomial coefficient: n choose n1, n2, n3, ...
n choose n1, n2, n3, ...
If the elements in ni do not sum to n
ni
n