Hay una versión un poco más general de este resultado: si un campo finito tiene un número impar $q$ de elementos, entonces la ecuación $x^2+1=0$ tiene una solución en el campo si y sólo si $q \equiv1\pmod4 $ . En su caso, por supuesto, puede tomar $ \mathbb Z/p \mathbb Z$ como un campo finito, en cuyo caso $q=p$ .
Aquí hay una simple prueba. Agrupa todos los elementos no cero de tu campo en paquetes $\{a,-a,a^{-1},-a^{-1}\}$ . Es fácil ver que si se hubiera generado el paquete a partir de cualquiera de sus tres elementos que no sean $a$ seguiría dando el mismo paquete. En otras palabras, uno tiene particionado el $q-1$ elementos no cero en paquetes. Sin embargo, no todos los paquetes tienen cuatro elementos distintos: aunque no puede suceder que $a=-a$ (ya que $q$ es impar), puede suceder que $a=a^{-1}$ o que $a=-a^{-1}$ en ambos casos el paquete tiene dos elementos en lugar de cuatro.
Ahora $a=a^{-1}$ sucede si y sólo si $a^2=1$ y desde que $a^2-1=(a-1)(a+1)$ esto ocurre precisamente para $a=1$ y $a=-1$ dando un paquete de este tipo, sin importar el campo. Para el otro tipo $a=-a^{-1}$ esto sucede si y sólo si $a^2=-1$ . Eso no tiene por qué ocurrir, pero si ocurre, hay precisamente dos soluciones a esta ecuación cuadrática, que juntas definen un solo paquete de este segundo tipo. Ahora bien, como todos los paquetes restantes son de tamaño $4$ y todos los tamaños de paquetes suman $q-1$ se ve fácilmente que un paquete del segundo tipo existe si y sólo si $q-1$ es un múltiplo de $4$ .