En un post Encontrar eficazmente dos cuadrados que sumen un primo
Leo
"En 1825 Gauss dio la siguiente construcción para escribir un primo congruente a $1 \pmod{4}$ como una suma de dos cuadrados: Sea $p=4k+1$ sea un número primo. Determine $x$ (esto es únicamente posible...) para que
$$ x = \frac{(2k)!}{2(k!)^2} \pmod{p}, \quad |x| < \frac{p}{2}$$
Determina ahora $y$ para que
$$ y = x \cdot (2k)! \pmod{p}, \quad |y| < \frac{p}{2}$$
Gauss demostró que $x^2+y^2=p$ ."
Revisé el libro de Stark y no vi una referencia directa al documento o libro original de Gauss u otra referencia para explicar el método de cómo llegar a esto y encontrar la fórmula. Si conoces la referencia o puedes explicar su procedimiento, te lo agradeceré.