This problem derives from an expression of probability in random walk. I hope to
prove that \begin{equation*} \sum_{l=0}^{n} \binom{n}{l}^2 (x+y)^{2l} (x-y)^{2(n-l)} = \sum_{l=0}^{n} \binom{2l}{l} \binom{2(n-l)}{n-l} x^{2l}y^{2(n-l)} \end{equation*}
It seems quite right, I verified this equation from $n=1$ to $6$ without any idea to prove it.
This equality is found by @Sangchul Lee.
The RHS is the coefficient of $z^{2n}$ in the product between $$ \sum_{l\geq 0}\binom{2l}{l}(xz)^{2l}=\frac{1}{\sqrt{1-4x^2 z^2}}\qquad\text{and}\qquad \sum_{l\geq 0}\binom{2l}{l}(yz)^{2l}=\frac{1}{\sqrt{1-4y^2 z^2}} $$ i.e. the coefficient of $z^{2n}$ in $\frac{1}{\sqrt{(1-4x^2 z^2)(1-4y^2 z^2)}}$.
In the LHS we may recognize a Legendre polynomial: $$ P_n(t)=\frac{1}{2^n}\sum_{l=0}^{n}\binom{l}{n}^2(t+1)^l (t-1)^{n-l} $$ hence the given identity turns out to be a consequence of the following identity, about the generating function for Legendre polynomials: $$ \frac{1}{\sqrt{1-2xt+t^2}}=\sum_{n\geq 0}t^n P_n(x) $$ which can be proved through Rodrigues' formula or Bonnet's recursion formula.