Adenda. El cálculo en la otra respuesta no es del todo riguroso en el sentido de que la serie infinita que aparece no converge en una vecindad de cero para $z.$ Esto se puede arreglar utilizando un soporte Iverson diferente.
Supongamos que, como antes, tratamos de evaluar $$S = \sum_{r=0}^n {r+n-1\choose n-1} {3n-r\choose n}$$
que es $$S_2-S_1 = \sum_{r=0}^{2n} {r+n-1\choose n-1} {3n-r\choose n} - \sum_{r=n+1}^{2n} {r+n-1\choose n-1} {3n-r\choose n}.$$
Comience por evaluar $S_2.$
Poner $${3n-r\choose n} = \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{2n-r+1}} \frac{1}{(1-w)^{n+1}} \; dw.$$
y utilizar el siguiente soporte Iverson $$[[0\le r \le 2n]] = \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{z^{r}}{z^{2n+1}}\frac{1}{1-z} \; dz.$$
Esta segunda integral controla el rango de manera que podemos extender la suma hasta el infinito para obtener $$\frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{2n+1}} \frac{1}{(1-w)^{n+1}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{1}{z^{2n+1}} \frac{1}{1-z} \sum_{r=0}^\infty {r+n-1\choose n-1} z^r w^r \; dz \; dw.$$
Esto se simplifica a $$\frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{2n+1}} \frac{1}{(1-w)^{n+1}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{1}{z^{2n+1}} \frac{1}{1-z} \frac{1}{(1-wz)^n} \; dz \; dw.$$
Evaluamos la integral interna utilizando el hecho de que los residuos en los tres polos suman cero. El residuo en $z=0$ es la suma $S_2$ que estamos tratando de calcular. El residuo en $z=1$ rinde
$$-\frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{2n+1}} \frac{1}{(1-w)^{2n+1}} \; dw = -{2n+2n\choose 2n} = -{4n\choose 2n}.$$
Para el residuo en $z=1/w$ reescribir la integral de la siguiente manera: $$\frac{(-1)^n}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{3n+1}} \frac{1}{(1-w)^{n+1}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{1}{z^{2n+1}} \frac{1}{1-z} \frac{1}{(z-1/w)^n} \; dz \; dw.$$
Necesitamos una derivada que calculamos mediante la regla de Leibniz: $$\frac{1}{(n-1)!} \left(\frac{1}{z^{2n+1}} \frac{1}{1-z}\right)^{(n-1)} \\ = \frac{1}{(n-1)!} \sum_{q=0}^{n-1} {n-1\choose q} (-1)^q \frac{(2n+q)!}{(2n)! \times z^{2n+1+q}} \frac{(n-1-q)!}{(1-z)^{1+n-1-q}} \\ = \sum_{q=0}^{n-1} {2n+q\choose 2n} (-1)^q \frac{1}{z^{2n+1+q}} \frac{1}{(1-z)^{n-q}}.$$
Evaluar en $z=1/w$ para conseguir
$$\sum_{q=0}^{n-1} {2n+q\choose 2n} (-1)^q w^{2n+1+q} \frac{w^{n-q}}{(w-1)^{n-q}}.$$
Sustituya esto de nuevo en la integral en $w$ para obtener $$\sum_{q=0}^{n-1} {2n+q\choose 2n} (-1)^q \frac{(-1)^n}{2\pi i} \int_{|w|=\epsilon} \frac{1}{w^{3n+1}} \frac{1}{(1-w)^{n+1}} \frac{w^{3n+1}}{(w-1)^{n-q}} \; dw \\ = \sum_{q=0}^{n-1} {2n+q\choose 2n} \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{1}{(1-w)^{2n-q+1}} \; dw = 0.$$
Hemos demostrado que $$S_2 = {4n\choose 2n}.$$
Continuando con $S_1$ vemos que $$S_1 = \sum_{r=0}^{n-1} {r+2n\choose n-1} {2n-1-r\choose n} = \sum_{r=0}^{n-1} {r+2n\choose n-1} {2n-1-r\choose n-1-r}.$$
Para esta suma no se necesita el corchete de Iverson ya que la segunda binomial controla el rango a través de la siguiente integral:
$${2n-1-r\choose n-1-r} = \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n-1-r}}{w^{n-r}} \; dw.$$
Además, introduzca $${r+2n\choose n-1} = \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{r+2n}}{z^n} \; dz.$$
Esto da la integral $$\frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n-1}}{w^{n}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{2n}}{z^n} \sum_{r\ge 0} \frac{w^r (1+z)^r}{(1+w)^r} \; dz \; dw \\ = \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n-1}}{w^{n}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{2n}}{z^n} \frac{1}{1-w(1+z)/(1+w)} \; dz \; dw \\ = \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n}}{w^{n}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{2n}}{z^n} \frac{1}{1+w-w(1+z)} \; dz \; dw \\ = \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n}}{w^{n}} \frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{(1+z)^{2n}}{z^n} \frac{1}{1-wz} \; dz \; dw.$$
Extrayendo el residuo interior obtenemos $$\sum_{q=0}^{n-1} {2n\choose n-1-q} w^q$$ que da como resultado $$\sum_{q=0}^{n-1} {2n\choose n-1-q} \frac{1}{2\pi i} \int_{|w|=\epsilon} \frac{(1+w)^{2n}}{w^{n-q}} \; dw \\= \sum_{q=0}^{n-1} {2n\choose n-1-q} {2n\choose n-1-q}.$$
Esto es $$\sum_{q=0}^{n-1} {2n\choose q}^2$$ que puede ser evaluado por inspección como en la otra respuesta, y recogiendo las contribuciones de $S_2$ y $S_1$ volvemos a obtener
$$\frac{1}{2} \left({4n\choose 2n} + {2n\choose n}^2\right).$$
Observación. Para ser totalmente rigurosos también debemos demostrar que el residuo en infinito de
$$\frac{1}{2\pi i} \int_{|z|=\epsilon} \frac{1}{z^{2n+1}} \frac{1}{1-z} \frac{1}{(1-wz)^n} \; dz$$
es cero. Recordemos la fórmula del residuo en el infinito $$\mathrm{Res}_{z=\infty} h(z) = \mathrm{Res}_{z=0} \left[-\frac{1}{z^2} h\left(\frac{1}{z}\right)\right]$$
que en este caso da como resultado $$-\mathrm{Res}_{z=0} \frac{1}{z^2} z^{2n+1} \frac{1}{1-1/z} \frac{1}{(1-w/z)^n} \\ = -\mathrm{Res}_{z=0} z^{2n} \frac{1}{z-1} \frac{1}{(1-w/z)^n} \\ = -\mathrm{Res}_{z=0} z^{3n} \frac{1}{z-1} \frac{1}{(z-w)^n}$$ que es cero por inspección.