pub fn crt(x: u128, ps: &[u16]) -> Vec<u16>
Compute the CRT representation of x with respect to the primes ps.