Experimented ST-monad version of divModPolynomial, but it turns out to be inefficient (see division-bench.html)