Dejemos que $x$ representan los dos primeros dígitos y $y$ representan los dos segundos dígitos. Así,
$100x+y=x^2+y^2$
$x(100-x)=y(y-1)$
El lado derecho está parejo. Módulo $5$ el lado izquierdo es $\in\{0,1,2\}$ y el lado derecho es $\in\{0,1,4\}$ . Así, ambas partes terminan con $0$ o $6$ . Además, el lado izquierdo debe ser $\in\{0,1\}\bmod3$ y el lado derecho $\in\{0,2\}\bmod3$ por lo que ambos lados son múltiplos de $3$ .
Así que necesitamos $x$ sea un múltiplo de $3$ o uno más que tal, $x(100-x)$ para terminar en $0$ o $6$ (forzando $x$ para terminar en uno de los dígitos $0,2,8$ ), y $4[x(100-x)]+1=(2y-1)^2$ para ser un cuadrado perfecto. Nos propusimos probar los valores elegibles de dos dígitos de $x$ de $10$ hasta $50$ (nueve ensayos), sabiendo que cualquier solución para $x$ en este rango irá acompañada de una solución $100-x$ por simetría del producto $x(100-x)$ en el lado izquierdo. Los bloques son casos en los que no conseguimos el cuadrado perfecto necesario para $4[x(100-x)]+1=(2y-1)^2$ :
$x=10, x(100-x)=900, 4×900+1=3601$
$x=12, x(100-x)=1056, 4×1056+1=4225=65^2$
$x=18, x(100-x)=1476, 4×1476+1=5905$
$x=22, x(100-x)=1716, 4×1716+1=6865$
$x=28, x(100-x)=2016, 4×2016+1=8065$
$x=30, x(100-x)=2100, 4×2100+1=8401$
$x=40, x(100-x)=2400, 4×2400+1=9601$
$x=42, x(100-x)=2436, 4×2436+1=9745$
$x=48, x(100-x)=2496, 4×2496+1=9985$
Así que sólo $x=12$ trabaja entre los números de abajo $50$ de la que $100-x=88$ también funcionará entre los números superiores a $50$ y para este par de soluciones vemos $2y-1=65, y=33$ . Por lo tanto,
$1233=12^2+33^2$
$8833=88^2+33^2$