Sea $p$ un primo impar. Sabemos que el polinomio $x^{p-1}-1$ se descompone en factores lineales módulo $p$. Si $p$ es de la forma $4k+1$ entonces podemos escribir $$x^{p-1}-1=x^{4k}-1=(x^{2k}+1)(x^{2k}-1).$$ El teorema de Lagrange nos dice que cualquier congruencia de polinomios de grado $n$ módulo $p$ tiene a lo sumo $n$ soluciones. Por lo tanto, podemos deducir de esta factorización que $-1$ es un residuo cuadrático módulo $p$. De manera similar, si $p$ es de la forma $3k+1$ podemos escribir $4(x^{p-1}-1)=4(x^{3k}-1)=(x^k-1)((2x^{k}+1)^2+3)$ y deducir que $-3$ es un residuo cuadrático módulo $p.
¿Podemos demostrar de esta manera que $-2$ es un residuo cuadrático módulo $p$ si $p$ es de la forma $8k+1$ o $8k+3$?
Nota que estoy interesado solo en este método específico. Sé cómo probar esto usando diferentes medios.