La proposición. Para un entero positivo $p$, $\left(p^2+1\right)^2$ es divisible por algún número entero positivo $d\equiv 1\pmod{p}$ tal que $$d\notin\left\{1,p^2+1,\left(p^2+1\right)^2\right\}$$ si y sólo si
$$p=t_s(b)\text{ for some integers }b\geq 1\text{ and }s\geq 2\,,$$
donde
$$t_j(b):=\frac{\left(\frac{\left(b^2+2\right)+b\,\sqrt{b^2+4}}{2}\right)^j-\left(\frac{\left(b^2+2\right)-b\sqrt{b^2+4}}{2}\right)^j}{\sqrt{b^2+4}}$$
para todos los enteros positivos $b$ e $j$. Además, si $p$ es en esta forma $t_s(b)$, entonces hay al menos dos enteros positivos no de la forma $1$, $p^2+1$, e $\left(p^2+1\right)^2$ que dividen $\left(p^2+1\right)^2$, que son
$$t_s(b)\,t_{s-1}(b)+1\text{ and }\left(b^2+2\right)\,\left(t_s(b)\right)^2-t_s(b)\,t_{s-1}(b)+1\,.$$
Por último, el único primer número natural $p$ con esta propiedad es $p=3$.
Supongamos que $p$ e $d$ son enteros positivos tales que a$d\equiv 1\pmod{p}$ e $d$ divide $\left(p^2+1\right)^2$. Sin pérdida de generalidad, podemos suponer que la $d\leq p^2+1$; de lo contrario, reemplace $d$ por $\dfrac{\left(p^2+1\right)^2}{d}$. Debido a $d\equiv 1\pmod{p}$, tenemos
$$d=pk+1\text{ for some integer }k\text{ such that }0\leq k\leq p\,.$$
Como $d\mid \left(p^2+1\right)^2$, obtenemos que $d=pk+1$divide
$$\left(p^2+1\right)^2+(pk+1)(-2p^2+pk-1)=p^2\,\left(p-k\right)^2\,.$$
Desde $\gcd\left(pk+1,p^2\right)=1$, se deduce que
$$\frac{(p-k)^2}{pk+1}\in\mathbb{Z}\,.$$
Por supuesto, $k=0$ e $k=p$ son las posibles soluciones, pero que conducen a la trivial divisores $d=1$ e $d=p^2+1$.
De Zvi la respuesta en este enlace, todas las soluciones integrales $(p,k)$ con $p>k>0$ tomar la forma
$$(p,k)=\big(t_{j+1}(b),t_j(b)\big)\text{ for some positive integers }b\text{ and }j\,,$$
donde $\left\{t_j(b)\right\}_{j\in\mathbb{Z}_{\geq 0}}$ es una secuencia definida por las condiciones iniciales $t_0(b):=0$ e $t_1(b):=b$, junto con la recurrencia
$$t_{j}(b):=\left(b^2+2\right)\,t_{j-1}(b)-t_{j-2}(b)\text{ for all integers }j\geq 2\,.$$
Para un dado positivo intger $b$, tenga en cuenta que
$$t_j(b)=\frac{\left(\frac{\left(b^2+2\right)+b\,\sqrt{b^2+4}}{2}\right)^j-\left(\frac{\left(b^2+2\right)-b\sqrt{b^2+4}}{2}\right)^j}{\sqrt{b^2+4}}\text{ for all }j=0,1,2,\ldots\,.$$
A partir de aquí, llegamos a la conclusión de que todos los enteros positivos $p$ para que $\left(p^2+1\right)^2$ es divisible por algún número entero positivo $d\equiv 1\pmod{p}$ no de la forma $1$, $p^2+1$o $\left(p^2+1\right)^2$ debe satisfacer
$$p=t_{s}(b)\text{ for some integers }b>0\text{ and }s>1\,.$$
Entre $1$ a $100$, enteros positivos $p$ en dicho formulario son
$$t_2(1)=3\,,\,\,t_3(1)=8\,,\,\,t_2(2)=12\,,\,\,t_4(1)=21\,,\,\,t_2(3)=33\,,$$
$$t_5(1)=55\,,\,\,t_3(2)=70\,,\text{ and }t_4(1)=72\,.$$
Por ejemplo, con $p:=t_2(2)=12$, tenemos $\left(p^2+1\right)^2=145^2$, el cual es divisible por $$d:=5^2=25\equiv 1\pmod{p}\,,$$ which is not of the form $1$, $p^2+1$, or $\a la izquierda(p^2+1\right)^2$. In particular, if $p$ es el primer número natural, entonces debemos
$$p=t_2(1)=3\,,$$
donde $\left(p^2+1\right)^2=10^2=100$ tiene un divisor $d:=4\equiv 1\pmod{p}$ no de la forma $1$, $p^2+1$o $\left(p^2+1\right)^2$.