Notice,
(xy−1)(yz−1)(zx−1)=x2y2z2−(x+y+z)xyz+xy+yz+zx−1.
The problem is equivalent to finding distinct positive integers x,y,z such that xyz∣(xy+yz+zx−1).
Assume without loss of generality that x>y>z.
Then xyz≤xy+yz+zx−1. If z=1, then xy≤xy+y+x−1, which is a contradiction since y is a positive integer.
(2) When z=2, 2xy∣(xy+2y+2x−1)⇒2xy≤xy+2y+2x−1⇒xy≤2y+2x−1.
Since z=2, we know y=3.
Thus, 6x∣(5x+5)⇒6x≤5x+5⇒x≤5.
Combining x>y=3, we know x=4,5.
Upon inspection, x=5 satisfies the condition.
Therefore, the positive integer solutions (x,y,z) that satisfy the condition are (2,3,5) and their permutations.