pub fn poly_mul(a: &[BigInt], b: &[BigInt], p: &BigInt) -> Vec<BigInt>
Polynomial product mod p (coefficients low-to-high).
p