A functional equation and its application to the characterization of gamma distributions

The functional equation $$ f\left(x\right)g\left(y\right)=p\left(x+y\right)q\left(\frac{x}{y} \right) $$ is investigated for almost all ${\left(x,\,y\right)\in\mathbb{R}^{2}_{+}}$ and for the measurable functions ${f,\,g,\,p,\,q:\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}}$ . This equation is related to the Lukács characterization of gamma distribution.