Posting my own answer to my question after a few months because:
- This answer is directly related to the source of the question.
- This answer shows the numerator is a power of 2 without finding the explicit closed formula for the value.
This question originally came from taking this question and asking myself the obvious generalization:
Find $g(x)$, a polynomial of degree $2n+1$, such that $g(x)+1$ is divisible by $(x-1)^{n+1}$ and $g(x)-1$ is divisible by $(x+1)^{n+1}$.
One answer is to note that $g'(x)$ must be divisible by $(1-x)^n$ and $(1+x)^n$, hence must be a multiple of $(1-x^2)^{n}$. Find the anti-derivative $G(x)$ of $(1-x^2)^n$ with $G(0)=0$, and then your answer is $g(x)=\frac{1}{G(1)}G(x)$.
Turns out, $G(1)$ is the value in my question above. Indeed, that was the origin of my question. I started computing these values and saw the pattern.
But one new answer to that original question gives a clear reason why there is a power of $2$ in the numerator of $G(1)$. I've adjusted that argument to handle the more general case:
Note that $$\begin{align}2^{2n+1} &= \left((1+x)+(1-x)\right)^{2n+1}\\
&=\sum_{k=0}^{n}\binom{2n+1}{k}\left((1+x)^k(1-x)^{2n+1-k} +(1+x)^{2n+1-k}(1-x)^{k}\right)\\
&=\sum_{j=0}^{n}\binom{2n+1}{n-j}\left((1+x)^{n-j}(1-x)^{n+1+j} + (1+x)^{n+1+j}(1-x)^{n-j}\right)
\end{align}$$
Dividing by $2^{2n}$ and rearranging, we get:
$$1-\frac{1}{2^{2n}}\sum_{j=0}^{n} \binom{2n+1}{n-j}(1+x)^{n+1+j}(1-x)^{n-j} = -1 + \frac{1}{2^{2n}}\sum_{j=0}^{n} \binom{2n+1}{n-j}(1+x)^{n-j}(1-x)^{n+1+j}$$
So if $g(x)=1-\frac{1}{2^{2n}}\sum_{j=0}^{n} \binom{2n+1}{n-j}(1+x)^{n+1+j}(1-x)^{n-j}$, then $g(x)-1$ is divisible by $(1+x)^{n+1}$ and $g(x)+1$ is divisible by $(x-1)^{n+1}$.
Now, the lead coefficient of $g(x)$ is $\frac{(-1)^n}{(2n+1)G(1)}$ in the original solution for $g(x)$, and it is $$\frac{1}{2^{2n}}\sum_{j=0}^{n}\binom{2n+1}{n-j}(-1)^{n-j+1}=\frac{M}{2^{2n}}$$ for an integer $M$ in this new solution.
So we only need to know that there is a unique solution for $g(x)$. That's relatively easy to do: If $g_1,g_2$ are solutions, then $g_1(x)-g_2(x)$ must be divisibly by both $(x-1)^{n+1}$ and $(x+1)^{n+1}$, and hence must be zero since it is of degree at most $2n+1.$
So this shows that $G(1)$ must have a numerator equal to a power of $2$.
Indeed, we get the general formula for $f(n)$ if we prove:
$$M=\sum_{j=0}^n (-1)^{j+1}\binom{2n+1}{j} = (-1)^n\binom{2n}{n}$$
Then we get that $f(n)=\frac{2^{2n}}{(2n+1)\binom{2n}{n}}$.