I have prove that, for any given positive integer $p,$ parametric solution of the Diophantine equation
$$\frac{1}{p}=\frac{1}{x}+\frac{1}{y}$$
can be written in the form $x=ac(a+b),y=bc(a+b),$ where $p=abc.$
Proof
Let
$\frac{1}{p}=\frac{1}{x}+\frac{1}{y} ,x,y∈Z^+.$
Then $x+y=t$ and $xy=pt$ for some $t∈Z^+.$
Now the quadratic equation $z^2-tz+pt=0$ has two integer roots $x,y.$
Discriminant of this equation can be written as $Δ_(x,y)=t^2-4pt=q^2, q∈Z^+.$
The quadratic equation $t^2-4pt-q^2=0$ gives the value of $t.$
$Δ_t=16p^2+4q^2=4r^2,r∈Z^+.$
$4p^2+q^2=r^2,r∈Z^+.$
This equation is of the form of Pythagoras equation.
Therefore $p=abc,q=(a^2-b^2 )c$ and $r=(a^2+b^2 )c$ where $a,b,c$ are parameters.
Backward substitution gives that $t=(a+b)^2 c.$
Hence we can obtain that
$x=ac(a+b),y=bc(a+b).$
Then I was try to find the general parametric solution of the Diophantine equation
$$\frac{1}{p}=\frac{1}{x}+\frac{1}{y}+\frac{1}{z} ,x,y,z∈Z^+.$$
I have found some particular solutions like,
$$\frac{1}{n}=\frac{1}{n+2}+\frac{1}{n(n+1)}+\frac{1}{(n+1)(n+2)}$$
$$\frac{1}{n}=\frac{1}{n+1}+\frac{1}{n(n+2)}+\frac{1}{n(n+1)(n+2)}$$
$$\frac{1}{n}=\frac{1}{n+1}+\frac{1}{(n+1)^2} +\frac{1}{n(n+1)^2 }$$
$$\frac{1}{n}=\frac{1}{n+1}+\frac{1}{n(2n+1)}+\frac{1}{(n+1)(2n+1)}$$
$$\frac{1}{n}=\frac{1}{n+1}+\frac{1}{(n^2+n+1)}+\frac{1}{n(n+1)(n^2+n+1)}.$$
But still I have no idea about how to attack the general one.
Here I have three questions.
1) Is there any different proof for general solution of first equation than my proof ?
2) Is there any general parametric solution for the second Diophantine equation ?
3) Is there any reference for these type of Diophantine equations ?
Got it. Your equation is $$ xy = px + py, $$ $$ xy - px - py = 0, $$ $$ xy - px - py + p^2 = p^2, $$ $$ (x-p)(y-p) = p^2. $$ Apparently this observation occurs at Number of solution for $xy +yz + zx = N$
All solutions are given by finding a divisor $w$ of $p^2,$ with triple $$ \color{magenta}{ \left( p, \; \; p + w, \; \; p + \frac{p^2}{w} \; \right).} $$ If $w < p$ these are in order, if $w=p$ it is just $(p,2p,2p),$ if $w > p$ it is a repeat but out of order. So, the total number of solutions is $$ \frac{1 + d(p^2)}{2}, $$ where $d(n)$ is the number of positive divisors of $n.$
Note that the primitive triples, $\gcd(p,x,y),$ come when my $w$ is $1$ or some other square, so $p^2/w$ is also a square, in addition we require $\gcd(w,p^2/w)= 1$; for example $(6,10,15)$ with $w=4$ and $p^2/w = 9.$
OR $$ (30,31,930); \; \; (30,34,255); \; \; (30,39,130); \; \; (30,55,66). $$
$p$ up to $30.$
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
All primitive solutions are given by finding a divisor $w$ of $p$ such that $\gcd(w,p/w) = 1$ with triple $$ \color{magenta}{ \left( p, \; \; p + w^2, \; \; p + \frac{p^2}{w^2} \; \right).} $$ To keep them ordered we also choose $w \leq \sqrt p.$ If $p$ is a square in the first place, larger than $1,$ then $w=\sqrt p$ does not ever give a primitive solution anyway, that just gives $(p,2p,2p).$
Here are just the primitive ones for $p \leq 30$ and then $p=210.$
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=