\(\left(xy-1\right)|\left(x^3+x\right)\Rightarrow\left(xy-1\right)|x\left(x^2+1\right)\)mà \(\left(x,xy-1\right)=1\)nên \(\left(xy-1\right)|\left(x^2+1\right)\)
\(\Leftrightarrow\left(xy-1\right)|\left(x^2+1+xy-1\right)\Leftrightarrow\left(xy-1\right)|\left(x+y\right)\).
Đặt \(x+y=z\left(xy-1\right)\Leftrightarrow x+y+z=xyz\).
Không mất tính tổng quát, giả sử \(x\ge y\ge z\)thì \(xyz=x+y+z\le3x\Leftrightarrow3\ge yz\ge z^2\Rightarrow z=1\Rightarrow y\in\left\{1;2;3\right\}\).
Thử từng trường hợp của \(y\)chỉ thấy \(y=2\)có nghiệm \(x=3\)thỏa mãn.
Vậy phương trình ban đầu có các nghiệm là: \(\left(1,3\right),\left(1,2\right),\left(2,3\right),\left(2,1\right),\left(3,2\right),\left(3,1\right)\).