- ex p = poly.collect(x);
- ex cont = p.lcoeff(x);
- for ( int i=p.degree(x)-1; i>=p.ldegree(x); --i ) {
- cont = gcd(cont, p.coeff(x,i));
- if ( cont == 1 ) break;
- }
- ex pp = expand(normal(p / cont));
+ ex unit, cont, pp;
+ poly.unitcontprim(x, unit, cont, pp);