The polynomial whose roots are all real

Suppose $p(x)=a_nx^n+a_{n-1}x^{n-1}+\cdots+a_0\in \mathbb{R}[x]$ is a polynomial whose roots are all real where $a_n=1$. We want to show that

The polynomial $g(x)=\sum_{i=0}^{n} \binom{n}{i}a_ix^i\in \mathbb{R}[x]$ has only real roots.

This problem is the Exercise 5 of the first section of chapter 5 in the page 257 of the book "Algebra" by T.T.Moh.

I have tried a couple of days to solve it out, but fail in end.

Thanks.

Perhaps there's an easier way, but here's a sketch using some general facts:

For any polynomial $f$ it follows from Rolle's theorem that the derivative $f'$ has at most as many non-real zeros as $f$ itself. (If $f$ is of degree $n$ and has $k$ real zeros (and hence $n-k$ nonreal zeros), then $f'$ is of degree $n-1$ and has at least $k-1$ real zeros, since $f'$ has at least one zero between any two adjacent zeros of $f$; hence $f'$ has at most $(n-1)-(k-1) = n-k$ nonreal zeros.)

More generally (for any real number $\alpha$), $f'-\alpha f$ has at most as many non-real zeros as $f$ itself; this can be shown by counting zeros of $e^{-\alpha x} f(x)$ and its derivative using Rolle's theorem as above. (If $\alpha \neq 0$, then $f'-\alpha f$ is of degree $n$ instead of $n-1$, so one also has to take into account that $e^{-\alpha x} f(x) \to 0$ at either $+\infty$ or $-\infty$ in order to find at least one more real zero of its derivative.)

In other words, the operator $D-\alpha$ (where $D=d/dx$ and $\alpha\in\mathbf{R}$) doesn't increase the number of non-real zeros of a polynomial that it acts upon. Thus, if $p(x)=c (x-\alpha_1) \dots (x-\alpha_n)$ has only real zeros, then the operator $p(D)=c (D-\alpha_1) \dots (D-\alpha_n)$ doesn't increase the number of non-real zeros; in particular, it preserves the property of having only real zeros.

Another operation which preserves the property of having only real zeros is "reversing the coefficients", i.e., mapping $f(x) = \sum_0^n a_k x^k$ to $x^n f(1/x) = \sum_0^n a_k x^{n-k}$. (This is because the zeros of the new polynomial are exactly the reciprocals of the zeros of the original polynomial, except that the root $x=0$ may appear or disappear or change its multiplicity.)

Now let $p$ be your polynomial with only real zeros. Then the reversed polynomial $q(x)=\sum_0^n a_k x^{n-k}$ has only real zeros too, so the operator $q(D)$ doesn't increase the number of non-real zeros. Apply this operator to the polynomial $x^n$ (which has no non-real zeros). Then the resulting polynomial $r(x)=q(D)x^n=\sum_0^n b_k D^{n-k} x^n = n! \sum_0^n \frac{b_k}{k!} x^k$ has no non-real zeros. The conclusion is that we can divide the coefficients by $k!$ without ruining the property of having only real zeros. Now reverse $r(x)$ and perform our newly-discovered operation again; this has the effect of further dividing the coefficients by $(n-k)!$ while preserving the reality of the zeros, and then we're done.