gcd( a^4 + 2 a^2 b + b^2 - 2 a^2 c - 2 b c + c^2, 4 a^3 + 4 a b - 4 a c )