For any positive integer $n$, I conjecture that the polynomial
$$P(x)=\sum_{i=0}^{\lfloor n/2 \rfloor} (-1)^i\dfrac{(2i)!}{2^i\cdot i!}\dbinom{n}{2i}x^{n-2i}$$
has $n$ real roots.
I've checked some polynomials in Mathematica and I strongly believe this result is true.
Does anyone have any reference or proof for this result?
Any help would be appreciated.
Let's look at a small example. The first few such polynomials are:
We can compute the total number of real roots of $P_4$ via Sturm's theorem. The details are given at the linked Wikipedia article. The summary of how it works is that we apply the Euclidean algorithm to $P_4$ and $P_4'$, and remember the remainders we get at each step, with a suitable choice of sign. This is called the "canonical Sturm chain" for $P_4$. Then the the number of roots of $P_4$ in $(a,b)$ can be obtained by taking the difference between the number of sign changes we see in the canonical Sturm chain at $a$, and the number we see at $b$.
Explicitly, the polynomials we end up considering are:
Since the leading coefficient of each of these polynomials is positive and each has degree $1$ smaller than the one before, they alternate in sign near $-\infty$, so the sign sequence is $+-+-+$: there are $4$ changes of sign. Similarly, there are $0$ changes of sign near $\infty$ as all the polynomials are positive there. So Sturm's theorem tells us that $P_4$ has $4$ real roots.
The same argument works in general. We will establish the following two identities:
\begin{align} P_n'&=nP_{n-1}\\ x P_{n-1}-P_n&=(n-1)P_{n-2} \end{align} Together, these imply that the canonical Sturm chain for $P_n$ consists of the polynomials $P_k$ for $k \leq n$, each scaled by some positive constant multiple. Since each $P_n$ is monic, this in turn means that the canonical Sturm chain has $n$ sign changes at $-\infty$ and $0$ sign changes at $\infty$, and so $P_n$ has $n$ real roots.
Each of these identities can be established by direct calculation. We adopt the convention that any coefficients which have negative factorials in the denominator are $0$; this allows us to forget about the upper limit of the sum. So: $$P_n(x)=\sum_{i=0}^\infty (-1)^i \frac{(2i)!}{i!2^i}\binom{n}{i}x^{2n-i}=\sum_{i=0}^\infty (-1)^i \frac{n!}{(n-2i)!i!2^k}x^{2n-i}$$
With this in mind, we have: \begin{align} P_n'(x)&=\sum_{i=0}^\infty (-1)^i \frac{n!}{(n-2i)!i!2^i}(n-2i)x^{n-1-2i}\\ &=\sum_{i=0}^\infty (-1)^i \frac{n!}{(n-1-2i)!i!2^i}x^{n-1-2i}\\ &=n\sum_{i=0}^\infty (-1)^i \frac{(n-1)!}{(n-1-2i)!i!2^i}x^{n-1-2i}\\ &=nP_{n-1}(x) \end{align} and: \begin{align}P_n(x)-xP_{n-1}(x)&=\sum_{i=0}^\infty (-1)^i \left[ \frac{n!}{(n-2i)!i!2^i}-\frac{(n-1)!}{(n-1-2i)!i!2^i}\right]x^{n-2i}\\ &=\sum_{i=0}^\infty(-1)^i\frac{(n-1)!}{(n-1-2i)!i!2^i}\left(\frac{n}{n-2i}-1\right)x^{n-2i}\\ &=\sum_{i=1}^\infty (-1)^i \frac{(n-1)!}{(n-1-2i)!i!2^i}\frac{2i}{n-2i}x^{n-2i}\\ &=\sum_{i=1}^\infty(-1)^i \frac{(n-1)!}{(n-2i)!(i-1)!2^{i-1}}x^{n-2i}\\ &=(n-1)\sum_{i=1}^\infty (-1)^i \frac{(n-2)!}{(n-2i)!(i-1)!2^{i-1}}x^{n-2i}\\ &=-(n-1)\sum_{j=0}^\infty (-1)^j \frac{(n-2)!}{(n-2-2j)!j!2^j}x^{n-2-2j}\\ &=-(n-1)P_{n-2}(x) \end{align} which completes the proof.
(I don't know a whole lot about this, but Wikipedia tells me that polynomials which obey the identity $P_n'=nP_{n-1}$ form an Appell sequence, which means there may be a simpler or more intuitive way of thinking about all of this in terms of umbral calculus...)