Let $$p(x)=\sum_{i=0}^{n}a_ix^i$$ be a polynomial in $\mathbb{R}[x]$ such that $0$ and $1$ aren't roots of $p.$ Prove that $\frac{1}{\alpha}$ is a root of $p$ whenever $\alpha$ is a root of $p,$ and they have the same multiplicity iff $a_{n-k}=a_{k}$ or $a_{n-k}=-a_{k}$ for each $k\in\{0,1,\ldots,n\}.$
Basically, I am talking about polynomials like $p(x)=2x^2-x+2.$ Notice that the roots of this polynomial are reciprocals of each other.
My Attempt: I am able to prove the one-way implication. Let $$p(x)=\sum_{i=0}^{n}a_ix^i.$$ Let $a_{n-k}=a_k$ for each $k\in\{0,1,\ldots,n\}.$ Let $p(\alpha)=0.$ Then, $$a_n\alpha^n+a_{n-1}\alpha^{n-1}+\ldots+a_{n-1}\alpha+a_n=0.$$ Multiplying throughout by $\frac{1}{\alpha^n},$ we get that $\frac{1}{\alpha}$ is a root of $p$ too. If $a_{n-k}=-a_{k},$ we can multiply by $\frac{-1}{\alpha^n}$ to get the desired result.
How to prove the reverse implication?
[Note: Per the comments, we do NOT want $a_{n-k} = -a_k$.]
The complete classification of roots of coefficient-symmetric polynomials (without restriction on roots) is
We will prove this version, from which OP's question follows as a corollary.
Lemma: A polynomial of degree $n$ has symmetric coefficients if and only if $x^n f(\frac{1}{x} ) = f(x)$.
This is obvious by expanding the terms and comparing coefficients.
Forward direction of the classification (Roots of that structure result in symmetric coefficients):
For the roots that occur in pairs, we have $( x- \alpha_i )(x -\frac{1}{\alpha_i}) = ( x^2 - \beta_i x + 1 )$.
For $ \alpha = 1$, since $(x-1)(x-1) = (x^2 - 2x + 1)$, it has the above form with $ \beta_i = 2$.
Thus, $f(x) = A(x+1) ^m \prod ( x^2 - \beta_i x + 1 )$. Observe that $$x^n f( \frac{1}{x} ) = A[ x^m ( \frac{1}{x} + 1 ) ^m ] \prod x^2 ( \frac{1}{x^2 } - \beta_i \frac{1}{x} + 1 ) = A(x+1) ^m \prod ( x^2 - \beta_i x + 1 ) = f(x),$$
hence the polynomial has symmetric coefficients.
Backward direction of the classification (Coefficient-symmetric polynomials have that structure of roots):
If $a_{k} = a_{n-k}$, OP's work shows that the roots either pair up with $ \alpha , \frac{1}{\alpha}$, or must satisfy $ \alpha = \frac{1}{\alpha}$.
Hence, $f(x) = A(x-1)^a (x+1)^b \prod ( x^2 - \beta_i x + 1) $.
Since $A = a_n = a_0= f(0) = A(-1)^a (1)^b \prod 1$, hence $a $ must be even. So $ \alpha = 1$ appears with even multiplicity.
Thus, the roots are in the required form.