minor changes in modgcd() (univariate case)
modular GCD for polynomials in Z[x, y]
minor changes in modgcd() (bivariate case)
some improvements in modgcd() (bivariate case)
Implement Brown's GCD algorithm for polynomials in Z_p[x0, ..., x{k-1}]