$$\text{Buscar} \;\sum_{(a,b) \S}un
\qquad \text{donde} \qquad
S=\{(a,b)\in \mathbb Z^+ : (a-b)\sqrt{ab}=2016\}$$
Deje $(a,b) \in S$ y deje $x = \sqrt{ab}$. A continuación,$(a-b) = \dfrac{2016}{x}$.
Así
$$\left\{
\begin{array}{c}
a = \sqrt{ab + \frac 14(a-b)^2} + \frac 12(a-b) \\
b = \sqrt{ab + \frac 14(a-b)^2} - \frac 12(a-b) \\
\end{array}
\right\}
\implica
\left\{
\begin{array}{c}
a = \sqrt{x^2 + \dfrac{1008^2}{x^2}} + \dfrac{1008}{x} \\
b = \sqrt{x^2 + \dfrac{1008^2}{x^2}} - \dfrac{1008}{x} \\
\end{array}
\right\}$$
Necesitamos resolver $x^2 + y^2 = z^2$ donde $xy = 1008$.
Por ahora, debido a la simetría, podemos suponer $x < y$. A continuación, las posibilidades de $x$
$$x \in \{1, 2, 3, 4, 6, 7, 8, 9, 12, 14, 16, 18, 21, 24, 28\}$$
De trabajo mod 11.
Los residuos cuadráticos módulo $11$
\begin{array}{|c|cccccc|}
\hline
x & 0 & 1,10 & 2,9 & 3,8 & 4,7 & 5,6 \\
x^2 & 0 & 1 & 4 & 9 & 5 & 3 \\
\hline
\end{array}
Desde $xy = 1008 = 7 = 18 = 29 = 40 = 51 = 62 = 73 = 84 = 95 \pmod{11}$, las posibilidades de $x,y$ (que representan el $x, y$ simetría)
\begin{array}{|c|ccccc|}
\hline
x & 1 & 2 & 3 & 4 & 5 \\
y & 7 & 9 & 6 & 10 & 8 \\
x^2 + y^2 & 6 & 8 & 1 & 6 & 1 \\
\hline
\end{array}
Nos encontramos con $x \pmod{11} \in \{3,5,6,8\}$.
Este whittles la lista de abajo para
$$x \in \{3, 6, 8, 14, 16, 28\}$$
De trabajo mod 13.
Los residuos cuadráticos módulo $13$
\begin{array}{|c|cccccc|}
\hline
x & 0 & 1,12 & 2,11 & 3,10 & 4,9 & 5,8 & 6,7 \\
x^2 & 0 & 1 & 4 & 9 & 3 & 12 & 10 \\
\hline
\end{array}
Desde $xy = 1008 = 7 = 20 = 33 = 46 = 59 = 72 = 85 = 98 \pmod{13}$,
las posibilidades de $x,y$ (que representan el $x, y$ simetría)
\begin{array}{|c|ccccc|}
\hline
x & 1 & 2 & 3 & 4 & 6 & 8\\
y & 7 & 10 & 11 & 5 & 12 & 9\\
x^2+y^2 & 11 & 0 & 0 & 2 & 11 & 2 \\
\hline
\end{array}
Nos encontramos con $x \pmod{13} \in \{ 2, 3, 10, 11 \}$ Este whittles la lista de abajo para
$$x \in \{3, 16, 28\}$$
Calculamos el $x, y,$$z$.
\begin{array}{ccc}
x & y=\dfrac{1008}{x} & z=\sqrt{x^2+y^2} \\
\hline
3 & 336 & 336.0133926 \\
\color{red}{16} & \color{red}{63} & \color{red}{65} \\
28 & 36 & 45.607017 \\
\hline
\end{array}
Por lo tanto
\begin{array}{ccc|cc|c}
x & y=\frac{1008}{x}
& z=\sqrt{x^2+\left(\frac{1008}{x}\right)^2}
& a=z+y & b=z-y & (a-b)\sqrt{ab} \\
\hline
16 & 63 & 65 & 128 & 2 & 2016 \\
63 & 16 & 65 & 81 & 49 & 2016 \\
\hline
\end{array}
Por lo $\sum_{(a,b)\in S}a = 128 + 81 = 209$