The mean-field equations for the three-state Potts model $$ \mathcal{H}=-J\sum_{\langle ij \rangle}\delta_{\sigma_i\sigma_j}\:,\qquad\sigma_i=1,\,2,\,3 $$ is equivalent to $$ \mathcal{H}=-J\sum_{\langle ij \rangle}\vec{s}_i\cdot\vec{s}_j $$ with $$ \vec{s}_i= \begin{pmatrix} 1\\0 \end{pmatrix}, \begin{pmatrix} -1/2\\\sqrt{3}/2 \end{pmatrix}, \begin{pmatrix} -1/2\\-\sqrt{3}/2 \end{pmatrix}. $$
This is a problem from J. M. Yeomans book Statistical Mechanics of Phase Transitions and I was wondering how to show this for an arbitrary number of particles.
If we had just two particles it could easily be shown with a small table as follows:
$$ \begin{align} \text{State particle 1}&\quad |\quad\text{State particle 2}&\quad \quad\text{Potts}&\quad |\quad\text{Heisenberg}\\ 1&\quad |\quad 1&\quad 1&\quad |\quad 1\\ 1&\quad |\quad 2&\quad 0&\quad |\quad -1/2\\ 1&\quad |\quad 3&\quad 0&\quad |\quad -1/2\\ 2&\quad |\quad 1&\quad 0&\quad |\quad -1/2\\ 2&\quad |\quad 2&\quad 1&\quad |\quad 1\\ 2&\quad |\quad 3&\quad 0&\quad |\quad -1/2\\ 3&\quad |\quad 1&\quad 0&\quad |\quad -1/2\\ 3&\quad |\quad 2&\quad 0&\quad |\quad -1/2\\ 3&\quad |\quad 3&\quad 1&\quad |\quad 1\\ \end{align} $$
In both models we find the same resulting picture, which becomes more clear or obvious if we would use arrows or something, like when Potts equals to one then $\uparrow$, when Heisenberg is equal to one the $\uparrow$, when Potts equals 0 or Heisenberg $-1/2$ then $\downarrow$.
But I am not even sure if this is what Yeomans had in mind with this problem and how it could be expanded to $N$ particles instead of just two. I think a permutation might do it or that the expectation value of $\mathcal{H}$
$$ \langle\mathcal{H}\rangle=\frac{1}{Z}\sum_{configurations}\!\mathcal{H}\cdot e^{-\beta\mathcal{H}} $$
with $Z$ the partition function and $\beta=\frac{1}{k_BT}$ ($k_B$ the Boltzmann constant) should be the same for both models.
Or is there a more elegant / straightforward solution?