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