Assume $z_1,z_2,z_3 \in\Bbb C$ satisfy
$$(*)\quad\quad\quad\,\,\,\, z_1^2 \,+\, z_2^2 \,+\, z_3^2 \,\,=\,\, z_1z_2 \,+\, z_2z_3 \,+\, z_3z_1.$$
It is easy to see, by opening all parentheses, that $(*)$ holds iff
$$\,\,\,(**)\quad\,\, (z_1\,-\,z_2)^2 \,+\, (z_2\,-\,z_3)^2 \,+\, (z_3\,-\,z_1)^2 \,\,=\,\, 0.$$
First assume $z_j=z_k,$ for some $j\not=k.$ Then $(**)$ implies $z_1=z_2=z_3.$ This situation we accept as a degenerate equilateral triangle. Otherwise, the three $z_j$ are distinct points and we have a proper triangle.
The property of being equilateral is invariant under any rotation about the origin, and under any translation. It is easy to see that the property $(**)$ is also so invariant, in the first case via $(z_je^{i\theta}\,-\,z_ke^{i\theta})^2\,=\,e^{i2\theta}(z_j\,-\,z_k)^2,$ and in the second case via $((z_j\,+\,\zeta )\,-\,(z_k\,+\,\zeta ))^2\,=\,(z_j\,-\,z_k)^2.$
Therefore we may assume $z_1=0$ and $z_2=\xi>0.$ Our given $(*)$ is the quadratic $\,\xi^2\,+\,z_3^2\,=\,\xi z_3.$ This quadratic admits two solutions: $\,\xi(\frac{1}{2}\,+\,i\frac{\sqrt 3}{2})\,$ and $\,\xi(\frac{1}{2}\,-\,i\frac{\sqrt 3}{2}).$
The triangle is equilateral.