Consideremos
$$ \frac{1}{x} - \frac{1}{y} = \frac{1}{z}, $$
y
$$ h = \textrm{gcd}(x,y,z). $$
Digamos que
$$ \begin{eqnarray} x &=& h x',\\ y &=& h y',\\ z &=& h z',\\ \end{eqnarray} $$
entonces obtenemos
$$ \frac{1}{x'} - \frac{1}{y'} = \frac{1}{z'}, $$
donde
$$ \textrm{gcd}(x',y',z') = 1. $$
Podemos escribir
$$ \Big(y'-x'\Big)z' = x'y'. $$
Miremos a impar e incluso:
$$ \begin{array}{ccccl} x' & y' & z'\\ \textrm{even} & \textrm{even} & \textrm{even} & \Rightarrow & \textrm{gcd}(x',y',z') > 1.\\ \textrm{odd} & \textrm{odd} & \textrm{odd} & \Rightarrow& \textrm{contradiction as $ xy $ is odd but $ (x-y)z $ is even}.\\ \textrm{odd} & \textrm{even} & \textrm{even}\\ \textrm{even} & \textrm{odd} & \textrm{even} \end{array} $$
Así que escribimos
$$ \begin{eqnarray} y' &=& x' + 2p + 1,\\ z' &=& 2q, \end{eqnarray} $$
cuando
$$ x'^2 + \Big(2p+1\Big)x' = \Big(2p+1\Big)2q, $$
así que
$$ x' = - \frac{2p+1}{2} + \frac{2p+1}{2} \sqrt{1 + \frac{8q}{2p+1}}, $$
lo que significa que
$$ q = \frac{1}{2} r \Big( r + 1 \Big) \Big( 2p + 1 \Big), $$
para que
$$ \begin{eqnarray} x' &=& r \Big( 2p + 1 \Big),\\ y' &=& \Big( r + 1 \Big) \Big( 2p + 1 \Big),\\ z' &=& r \Big( r + 1 \Big) \Big( 2p + 1 \Big).\\ \end{eqnarray} $$
Sin embargo, $\textrm{gcd}(x',y',z') = 1$ Por lo tanto $2p+1=1$ De ahí que
$$ \begin{eqnarray} x &=& h r,\\ y &=& h \Big( r + 1),\\ z &=& h r \Big( r + 1 \Big).\\ \end{eqnarray} $$
Está claro que
$$ \begin{eqnarray} hxyz &=& \left[ h^2 r \Big( r + 1 \Big) \right]^2,\\ h(y-x) &=& \Big[ h r \Big]^2,\\ \end{eqnarray} $$
por lo tanto, ambos $hxyz$ y $h(y-x)$ son cuadrados perfectos.