Please check [my post](https://math.stackexchange.com/questions/3485228/the-condition-that-eulers-prime-generating-polynomial-is-a-composite-number).

It indicates that the number of divisors of $x^2+x+41$ is equal to the number of lattice points of $X^2+163Y^2-2(2x+1)Y-1=0$.

This formula is transformed in this way.

$$163X^2+163^2Y^2-2\cdot163(2x+1)Y=163$$
$$163X^2+\{163Y-(2x+1)\}^2-(2x+1)^2=163$$
$$163X^2+(163Y-2x-1)^2=4x^2+4x+164$$


$X':=163Y-2x-1,\ Y':=X$ and we divide both sides by 4,

$$\frac{{X'}^2+{163Y'}^2}{4}=x^2+x+41$$
$$N\left(\frac{X'+Y'\sqrt{-163}}{2}\right)=x^2+x+41$$

$\frac{X'+Y'\sqrt{-163}}{2}$ is an element of $\mathbb Q(\sqrt{-163})$.

This formula indicates the following:
The number of elements of $\mathbb Q(\sqrt{-163})$ with norm $x^2+x+41$ is linked to the number of divisors of $x^2+x+41$.


The phenomenon you are interested in is based on this fact.