Dado prime $p=8n+1$. Entonces
$$x^2-2py^2 = -1\tag1$$
es que no solucionable ,
$$p_1= 17, 73, 89, 97, 193, 233, 241, 257, 281, 337, 353, 401, 433, 449, 577, 593,601, 617, 641,\dots$$
pero es solucionable ,
$$p_2= 41, 113, 137, 313, 409, 457, 521, 569, 761, 809, 857, 953, 1129, 1201, 1321, 1601,\dots$$
Comparar a la de los números primos de la forma $p=u^2+32v^2$ (A105389):
$$p_3 = 41, 113, 137, \color{brown}{257}, 313, \color{brown}{337}, \color{brown}{353}, 409, 457, 521, 569, \color{brown}{577}, \color{brown}{593}, 761, 809, 857, \color{brown}{881}, 953, \dots$$
(También, $p_3$ tiene clase número $h(-p)$ divisible por 8.)
P: Es $p_2$ un subconjunto de a $p_3$?
(Lo que es equivalente, para $p=8n+1$, es cierto que una condición necesaria pero no suficiente tal que $(1)$ es solucionable, es que $p = u^2+32v^2$?)
He comprobado que todos solucionable $p = 8n+1 \leq 18089$ tiene la forma $u^2+32v^2$, pero no sé si todo solucionable $p$ tienen esa forma.
$\color{blue}{Edit}$: (En respuesta a Jagy la respuesta.) Los números primos de la forma $p=u^2+64v^2$ (A014754),
$$p_4 = 73, 89, 113, 233, 257, 281, 337, 353, 577, 593, 601, 617, 881, 937, 1033, 1049, 1097\dots$$
pero ni $p_1$ ni $p_2$ es un subconjunto de a $p_4$. Sin embargo, los números primos de la forma $p=u^2+64v^2=16n+9$,
$$p_5 = 73, 89, 233, 281, 601, 617, 937, 1033, 1049, 1097,\dots$$
como resultado de Dirichlet, es irresoluble por lo que es un subconjunto de a $p_1$.