<s>[INST] Evaluate $|z_1 + z_2 + z_3|$ on unit circle given $ \frac{ z_1^2}{z_2 z_3} + \frac{ z_2^2}{z_1 z_3} + \frac{ z_3^2}{z_1 z_2} +1 =0$ 
given three complex numbers $(z_1 , z_2 , z_3)$ lying on the unit circle and related by the equation $ \frac{ z_1^2}{z_2 z_3}  + \frac{ z_2^2}{z_1 z_3} + \frac{ z_3^2}{z_1 z_2} +1 =0$, find the sum of all possible values of $|z_1 + z_2 + z_3|$

I have no idea how to solve it systematically but with some guessing I found that $(1,1,-1) $ is a solution triplet for this constraint. And, if $(1,1,-1)$ is a solution then by the symmetry of how the equation is, it means that $(-1,1,1)$ and $(1,-1,1)$ are also solutions.
But, I'm not sure how to rule out/ find more solutions. Ideally, the answer I"m looking for is a general method to approach these kind of questions.
 [/INST] Let us define the angles $\alpha$ and $\beta$ as
\begin{eqnarray}
\alpha &=& \mathrm{Arg}\left[\frac{z_2}{z_1}\right] \, ,
\\
\beta &=& \mathrm{Arg}\left[\frac{z_3}{z_1}\right] \, .
\end{eqnarray}
Both the condition and the norm $|z_1 + z_2 + z_3|$ can be written in terms of only $\alpha$ and $\beta$ since
\begin{eqnarray}
0 &=& \frac{z^2_1}{z_2 z_3} + \frac{z^2_2}{z_3 z_1} + \frac{z^2_3}{z_1 z_2} + 1 \, 
\\
&=& \frac{1}{\frac{z_2}{z_1} \frac{z_3}{z_1}} + \frac{\left(\frac{z_2}{z_1}\right)^2}{\frac{z_3}{z_1}} + \frac{\left(\frac{z_3}{z_1}\right)^2}{\frac{z_2}{z_1}} + 1 \, ,
\end{eqnarray}
and
\begin{equation}
|z_1 + z_2 + z_3| = |z_1|\left|1 + \frac{z_2}{z_1} + \frac{z_3}{z_1}\right| = \left|1 + \frac{z_2}{z_1} + \frac{z_3}{z_1}\right| \, .
\end{equation}
This means that for each pair of values of $\alpha$ and $\beta$ that solve the condition above, $z_1$ (or, equivalently, its $\mathrm{Arg}$) still remains an unconstrained degree of freedom, which is probably considered a redundancy when computing the sum that the exercise asks for.
An additional source or redundancy is the fact that $\alpha$ and $\beta$ that solve the condition above are interchangeable, given the symmetry of the expression.
If we focus on only the ranges $0 \leq \alpha < 2 \pi$ and $0 \leq \beta < 2 \pi$, the pairs of values that solve the condition are
\begin{equation}
\left( \pi, \pi \right), \left( 0, \pi \right), \left( \frac{\pi}{3} , \frac{2\pi}{3} \right) , \left( \frac{\pi}{3} , \frac{5\pi}{3} \right), \left( \frac{4\pi}{3} , \frac{5\pi}{3} \right) \, ,
\end{equation}
where, again, each can be either $\alpha$ or $\beta$. The corresponding values of $|z_1 + z_2 + z_3|$ are
\begin{equation}
1, 1, 2, 2, 2.
\end{equation}
From here it's not completely clear how much redundancy should we consider in order to compute the sum, but the fundamental components are there. If we are expected to sum each possible value of the norm once then we get $1+2 = 3$.
</s>