1p := (x^2 - 3*x*y + y^2)^4*(3*x - 7*y + 2)^5; 2q := (x^2 - 3*x*y + y^2)^3*(3*x - 7*y - 2)^6; 3qqq:=gcd(q,p); 4 5factor(qqq); 6 7