Okay, so the expression just becomes xy(x^2-2y)+1=(x-1)(y-1)
Case 1: Suppose both x,y are positive. If x^2>2y the LHS grows faster than the RHS (xy>(x-1)(y-1) for positive x,y) so equality cannot be achieved. However if x^2<2y, then the second factor on lhs is negative. but if x and y are positive then rhs is non negatuve. hence equality can be achieved if xy (x^2-2y) is -1 which can be achieved via the solution (1,1)
Case 2: if x,y both zero, it is a solution (0,0)
Case 3: If one of x,y is greater than zero while the other is zero, then RHs is always negative while the LHS is 1
Hence the solution set is (1,1) and (0,0)