I am looking to factor a diophantine equation. However, I am unable to figure out the steps in order to reach the simplification:
Find all integral solutions to the equation:
$(x^2 + 1)(y^2 + 1) + 2(x − y)(1 − xy) = 4(1 + xy)$
However, I am unsure as to how this is equivalent to:
$x^2y^2 − 2xy +1+ x^2 + y^2 − 2xy + 2(x − y)(1 − xy) = 4$.
Please let me know the exact steps to achieve that answer.