How can I prove that if $g(X)\in \mathbb Q[X]$ and $\zeta\in\mathbb C\backslash \mathbb R$, therefore $$g(\zeta)\in\mathbb R\implies g(\zeta)=h(\zeta+\bar\zeta)$$ for a certain polynomial $h(X)\in\mathbb Q[X]$ ?
Let $$g(X)=a_0+a_1X+...+a_nX^n,$$ I know that $$g(\zeta)\in\mathbb R\implies g(\zeta)=\overline{g(\zeta)}\underset{a_i\in\mathbb Q}{=}g(\bar\zeta)$$
therefore $$a_1(\zeta-\bar\zeta)+a_2(\zeta^2-\bar\zeta^2)+...+a_n(\zeta^n-\bar\zeta^n)=0,$$ and thus $$\Im(\zeta^k-\bar\zeta^k)=0,$$ for all $k$, but I don't know how to continue.