pub fn multinomial<T>(k: &[T]) -> Twhere T: Integer<Output = T> + Clone + for<'a> Add<&'a T>,
Calculate the multinomial coefficient.