Let $p$ be prime s.t. $p=a^2+4b^2, a\equiv 1\,mod\,4$.
Then, for $f=X^3+X$, $a_E(p)=2a$?
(Here, $a_E(p):=p+1-\#\{(x,y)\in\{0,...p-1\}^2|x^3+x\equiv y^2\,mod\,p\}$)
Let $p$ be prime s.t. $p=a^2+4b^2, a\equiv 1\,mod\,4$.
Then, for $f=X^3+X$, $a_E(p)=2a$?
(Here, $a_E(p):=p+1-\#\{(x,y)\in\{0,...p-1\}^2|x^3+x\equiv y^2\,mod\,p\}$)
Copyright © 2021 JogjaFile Inc.
For the direct reference i will cite a related result from [Koblitz, Introduction to elliptic curves and modular forms]. It is a Theorem from II, §2, The zeta function of $E_n$, pages 59, 60, where $E_n$ is the "congruent number problem" elliptic curve: $$ E_n\ :\ y^2 = x^3-n^2x\ . $$ The results needed to prove the result in the sequel are using the Gauss and Jacobi sums, defined for a multiplicative character $\chi:\Bbb F_q^\times\to\Bbb C^\times$ by $$ \begin{aligned} g(\chi) &= \sum_x\chi(x)\psi(x)\ ,\\ J(\chi_1,\chi_2) &= \sum_x\chi_1(x)\chi_2(1-x)\ . \end{aligned} $$ (Above $\psi$ is an additive character, the sums are over elements in the field.)
It is related, because the given curve $Y^2=X^3+X$ and the curve $E_2$ of equation $y^2=x^3-4x$ are isogenous. Now let us cite / state:
This is a generalization of the posted result.
The result for $p\equiv 3$ modulo four is simple to obtain, the polynomial $f$ given by $x^3+x$, respectively $x^3-n^2x$ being odd, so we can group points $\ne 0$ in $\Bbb F_p$ in pairs $\{x,-x\}$, and exactly one of $f(x)$, $f(-x)=-f(x)$ is a squaer.
The result for $p\equiv 1$ modulo four and $n=2$ leads to the posted question. (The $2$ is a square, so $\left(\frac2p\right)$ is $=1$.)
Comment:
Here are some numerical examples, supported by sage code for the first primes:
These explicit numerical examples are here, hoping that the questions gets more support.