How many lattice paths are there from $(0,0)$ to $(2n,2n)$ that do not go through one of the points $(2i-1,2i-1)$ for $i=1,\dots,n$?
My idea is to count the number of total lattice paths one can take from $(0,0)$ to $(2n,2n)$. There are ${4n \choose 2n}$ such paths. Then subtract the number of paths that are not valid. In counting these, I reasoned that we must avoid the "odd points" inside the grid with height and width of $2n$. I counted the number of paths that take these of points to be ${4 \choose 2}^{n-1}{2 \choose 1}{2 \choose 1}$ with the reasoning that from $(0,0)$ to $(1,1)$, there are ${2 \choose 1}$ paths, similarly for $(2n-1,2n-1)$ to $(2n,2n)$. Now, there are a total of $n-1$ "odd points" we consider and the number of paths from say $(1,1)$ to $(3,3)$ is ${4 \choose 2}$, we consider $n-1$ such scenarios. But in comparing my result, it is wrong, I seem to be undercounting the number of invalid paths that I need to subtract from the total paths.
Edit: The result is expected to be the Catalan numbers of the form $C_{2n+1}$.
Edit 2: I've reworked the problem to make the first couple of terms match $C_{2n+1}$, by removing from the total number of lattice paths the invalid paths (a sum of all the possible cases by which we choose how many and which odd points our invalid path has gone through). It seems to be some recursive function, any ideas how to express this recursively?
Let's call your sequence of valid, $(2i-1,2i-1)$ avoiding, paths $\langle a_k \rangle$. Then it has a recurrence, starting with $a_0=1$ shown below:
$$a_k=\sum_{i=1}^{k}2c_{2i-1}a_{k-i}\tag{1}$$
Since any path has some first point where it touches the diagonal at an even point $(2i,2i)$ there are $2c_{2i-1}$ catalan paths to this point from $(0,0)$ (one set of $c_{2i-1}$ paths above the diagonal and one below) then $a_{k-i}$ valid paths from $(2i,2i)$ to $(2k,2k)$. Hence the are $2c_{2i-1}a_{k-i}$ valid paths whose first intersection with the diagonal is $(2i,2i)$. Summing over all possible first diagonal intersection points $i=1,\ldots,k$ gives $(1)$.
Call the generating function for $\langle a_k \rangle$ $f(z)=\sum_{j\ge 0}a_jz^j$ and the catalan number generating function $C(z)$, then the odd catalan numbers have generating function
$$C_o(z)=\frac{1}{2}z^{1/2}(C(z^{1/2})+C(-z^{1/2}))=\sum_{j\ge 1}c_{2j-1}z^j\tag{2}$$
and even catalan number generating function
$$C_e(z)=\frac{1}{2}(C(z^{1/2})+C(-z^{1/2}))=\sum_{j\ge 0}c_{2j}z^j\tag{3}$$
then $(1)$ can be represented by the generating function relation
$$1+2C_o(z)f(z)=f(z)$$
so that
$$f(z)=(1-2C_o(z))^{-1}\tag{4}$$
It is known that the catalan number generating function is
$$C(z)=\frac{1}{2z}(1-\sqrt{1-4z})=\sum_{j\ge 0}c_jz^j$$
so $(2)$ and $(3)$ become
$$C_o(z)=\frac{1}{4}\left(2-\sqrt{1-4z^{1/2}}-\sqrt{1+4z^{1/2}}\right)\tag{2*}$$
$$C_e(z)=\frac{1}{-4z^{1/2}}\left(\sqrt{1-4z^{1/2}}-\sqrt{1+4z^{1/2}}\right)\tag{3*}$$
Now putting $(2\text{*})$ in $(4)$ gives
$$f(z)=\frac{2}{\sqrt{1-4z^{1/2}}+\sqrt{1+4z^{1/2}}}$$
then multiplying top and bottom by $\sqrt{1-4z^{1/2}}-\sqrt{1+4z^{1/2}}$ gives
$$\begin{align}f(z)&=\frac{2\left(\sqrt{1-4z^{1/2}}-\sqrt{1+4z^{1/2}}\right)}{(1-4z^{1/2})-(1+4z^{1/2})}\\[2ex] &=\frac{1}{-4z^{1/2}}\left(\sqrt{1-4z^{1/2}}-\sqrt{1+4z^{1/2}}\right)\\[2ex] &=C_e(z)\end{align}$$
hence $a_k=c_{2k}$ are the even catalan numbers
$$1,2,14,132,1430,\ldots$$