pub fn poly_divide<E: ExtElem>(p: &mut [E], z: E) -> E
In-place polynomial division.
Take the coefficients in P, and divide by (X - z) for some z, return the remainder.