Я использую SymPy для разложения многочленов по комплексным числам.
import sympy as sp
x, y = sp.symbols('x y')
z = sp.expand( (x - 1) * (y - 1) )
factorization_1 = sp.factor(z)
factorization_2 = sp.factor(z, extension = [sp.I])
Здесь z содержит выражение xy - x - y + 1, а factorization_1 действительно имеет ожидаемое значение (x - 1)(y - 1), но по какой-то причине factorization_2 имеет значение x - 1.
Откуда взялся коэффициент (y - 1) go?