Consider the polynomial with real coefficients $P(x)=x^6+a x^5+b x^4+c x^3+b x^2+a x+1$, and let $x_1, x_2, \ldots, x_6$ be its zeros. Prove that $$ \prod_{k=1}^6\left(x_k^2+1\right)=(2 a-c)^2 $$
By vieta, we know that $$-a=\sum_{cyc} x_1$$ and $$-c=\sum_{cyc} x_1 x_2 x_3$$ Also, $$P(x)=\prod_{i=1}^6(x-x_i)$$ How do I proceed from here?