How to prove $\int_{-\infty}^{\infty}\binom{n}{x}^2\binom{2n+x}{x} dx=\binom{2n}{n}^2$

Inspired by this question (How prove this $\sum_{k=0}^n \binom{n}{k} \binom{(p-1)n}{k} \binom{pn+k}{k} = \binom{pn}{n}^2$) and the identity given there for the easiest case $p=2$: $$\sum_{k=0}^{n}\binom{n}{k}^2\binom{2n+k}{k}=\binom{2n}{n}^2$$

I have checked the following identity numerically and it seems to be true:

$$\int_{-\infty}^{\infty}\binom{n}{x}^2\binom{2n+x}{x} dx=\binom{2n}{n}^2$$.

Any ideas how this could be proved?

This question is related to my (answered) question https://mathoverflow.net/questions/255758/integral-of-power-of-binomials-equal-to-sum-of-power-of-binomials but I thought it would be more appropriate to post it here, on Stack Exchange, because it is no real research question.

• Most likely a beta function type problem? Dec 22, 2016 at 21:54
• @Nemo it could have involved the beta function as a ratio of gamma functions. Dec 24, 2016 at 12:47

1 st proof. By elementary algebra $$\binom{n}{x}^2\binom{2n+x}{x}=\frac{(n!)^2}{(2n)!}\frac{\prod_{k=1}^{2n}(x+k)}{\prod_{k=1}^{n}(x-k)^2}\frac{\sin^2 \pi x}{\pi^2 x^2}.$$ Now one can apply residue theorem to the function $\frac{(n!)^2}{(2n)!}\frac{\prod_{k=1}^{2n}(z+k)}{\prod_{k=1}^{n}(z-k)^2}\frac{1-e^{2\pi i z}}{2\pi^2 z^2}$ which has only simple poles at $z=m,~0\le m\le n$ and asymptotics $O(1/z^2),~z\to\infty$ in the upper half plane. The contour of integration is the interval $(-R,R)$ with small semicircles near the poles, and closed by a semicircle of radius $R$ in the upper half plane.
If $m,n\in\mathbb{N},~0<\beta<1$ $$f(x)=\binom{m}{\beta x}\binom{n}{\beta x}\binom{m+n+\beta x}{m+n}$$ then for arbitrary $s\in\mathbb{C}$ $$\beta\sum_{k=-\infty}^\infty f(s+k)=\beta\int_{-\infty}^\infty f(x)dx=\binom{m+n}{n}^2.$$
An outline of the proof is given below. First by observing that $$f(x)=\frac{m!n!}{(m+n)!}\frac{\prod_{k=1}^{m+n}(\beta x+k)}{\prod_{k=1}^{m}(\beta x-k)\prod_{k=1}^{n}(\beta x-k)}\frac{\sin^2\pi\beta x}{\pi^2\beta^2 x^2},$$one can see that $f(x)$ does not have any poles and $f(z)\sim e^{2\pi\beta |\text{Im}z|}/|z|^2$, $z\to\infty$. Now consider the Fourier transform $$\int_{-\infty}^\infty f(x)e^{iax}dx.$$ Contour integration shows that if $|a|>2\pi\beta$ this integral is $0$. This means that $f(x)$ is band limited, i.e. its Fourier spectrum is limited to the band $|\omega|<2\pi\beta$. Now Poisson summation gives $$\sum_{k=-\infty}^\infty f(s+k)=\int_{-\infty}^\infty f(x)dx,\quad 2\pi\beta<2\pi.$$ This means that the sum on the left does not depend on $s$ and $\beta$. So one can put $s=0$, $\beta=1$ and use the result to compute this sum.$\Box$