The identity can be rewritten as
$$\sum_{\substack{|x|<p\\ 2|x}}\sum_{r|p^2-x^2}\left(\frac{-3}{r}\right)=p+2,$$
because for $x=0$ the inner sum is $1-1+1=1$. Writing $x=2c$, the identity becomes
$$\sum_{|c|<p/2}\,\sum_{r|p^2-4c^2}\left(\frac{-3}{r}\right)=p+2.$$
The inner sum counts the number of integral representations $p^2-4c^2=a^2+ab+b^2$ divided by $6$, hence the identity is equivalent to the statement that the number of integral representations of $p^2$ by the quadratic form $a^2+ab+b^2+4c^2$ equals $6(p+2)$. We shall verify this by Siegel's mass formula, as it appears in Über die analytische Theorie der quadratischen Formen, Ann. of Math. 36 (1935), 527-606.

The quadratic form $a^2+ab+b^2+4c^2$ has discriminant $12$, hence it is alone in its genus. So, by Siegel's mass formula, the number of representations can be calculated as a product of local densities:
$$r(p^2)=\alpha_\infty\alpha_2\alpha_3\alpha_p\prod_{q\nmid 6p}\alpha_q.$$
By Hilfssatz 26 and (71) and the line below (59) in Siegel's paper,
$$\alpha_\infty=\frac{p}{\sqrt{3}}\cdot\frac{\pi^{3/2}}{\Gamma(3/2)}=\frac{2\pi }{\sqrt{3}}p.$$
By Hilfssatz 13 in Siegel's paper,
$$\alpha_2=\frac{3}{2}\qquad\text{and}\qquad\alpha_3=\frac{4}{3}.$$
By Hilfssatz 16 in Siegel's paper,
$$\alpha_p=\left(1-p^{-2}\right)\left(1+\frac{p^{-1}}{1+p^{-1}}\right)=(1-p^{-1})(1+2p^{-1}).$$
Finally, by Hilfssatz 12 in Siegel's paper,
$$\prod_{q\nmid 6p}\alpha_q=\prod_{q\nmid 6p}(1+\chi(q)q^{-1})=\frac{2}{1-p^{-1}}\prod_{q\neq 3}(1+\chi(q)q^{-1}),$$
where $\chi$ denotes the nontrivial quadratic character modulo $3$. Therefore,
$$r(p^2)=(p+2)\frac{8\pi}{\sqrt{3}}\prod_{q\neq 3}(1+\chi(q)q^{-1}).$$
We can identify the product over $q\neq 3$ as
$$\prod_{q\neq 3}(1+\chi(q)q^{-1})=\prod_{q\neq 3}\frac{1-q^{-2}}{1-\chi(q)q^{-1}}=\frac{9}{8}\cdot\frac{6}{\pi^2}L(1,\chi)=\frac{3\sqrt{3}}{4\pi},$$
hence in the end
$$r(p^2)=(p+2)\frac{8\pi}{\sqrt{3}}\cdot\frac{3\sqrt{3}}{4\pi}=6(p+2).$$
The proof is complete.