Let $w = \dfrac{z+i}{z-i}$, then $w^{2016} = 1$, hence we solve for $w$ first. $$|w^{2016}| = |1+ 0i| = 1 \Rightarrow |w| = 1$$ $$\text{arg}(w^{2016}) = 0+2k\pi \Rightarrow \text{arg}(w) = \dfrac{k\pi}{1008}$$
Hence $w = e^{k\pi i/1008}$
And then $$w = \dfrac{z+i}{z-i} \Rightarrow z = \dfrac{-i(1+w)}{1-w}$$
We know that $w = e^{k\pi i/1008} = \cos \dfrac{k\pi}{1008} + i\sin \dfrac{k\pi}{2016}$
Denote $\dfrac{k\pi}{1008} = a$, now we expand and we have $$z = \dfrac{-i(1+\cos a + i\sin b)}{1- \cos a - i\sin a} = \dfrac{-i-i\cos a + \sin a}{1- \cos a - i\sin a}$$
Now recall we can use the denominator of $z$ and times it by the denominator's conjugate $$z = \dfrac{(-i-i\cos a + \sin a)(1-\cos a + i\sin a)}{(1- \cos a - i\sin a)(1-\cos a + i\sin a)} = \dfrac{\sin a}{1- \cos a} \in \mathbb{R}$$
I have came up with the above solution and i think it should be right, but then i read the answer and it says that $z$ is real if and only if $z = \overline{z}$. Now i dont dispute that but they claim that $\overline{z} = \dfrac{(-i)(-1-\overline{b})}{1-\overline{b}}$ where $b$ is just $e^{k\pi i/1008}$. I am confused over one thing, which is how they managed to know what $\overline{z}$ is. For me my understanding of conjugate of $z$ is changing the sign of $z$'s imaginary part. Which i clearly cant see any short cut?