The equation $x_1x_2+y_1y_2=n$ is well-studied (Ingham, Heath-Brown, Deshouillers & Iwaniec, Ismoilov) because it arises in an additive divisor problem. The number of solutions in positive integers is $n(c_2\log^2n+c_1\log n+c_0)+O(n^{1-\delta})$ where $c_i$'s are some explicit arithmetic functions of $n$.
The same tools allow to prove similar formula for the number of solutions of the equation $ax_1x_2+by_1y_2=n$. In this case $c_i$'s will depend on $n,a$ and $b$.
Is it possible to find this result in a literature?