pub fn crt_inv(xs: &[u16], ps: &[u16]) -> u128
Compute the value x given a list of CRT primes and residues.