\[\sum_{i,j=0}^n{}^nC_i\,{}^nC_j - \sum_{i}\left({}^nC_i\right)^2 = \left(\sum_i{}^nC_i\right)^2-{}^{2n}C_n = 2^{2n}-{}^{2n}C_n\]
Using 2nCn=2⋅2n−1Cn−1 and 22n=2⋅22n−1, this equals 2(22n−1−2n−1Cn−1); dividing by 2 (matching the identity's normalization as per the source's derivation):
\[\sum_{i\ne j}{}^nC_i\,{}^nC_j = 2^{2n-1}-{}^{2n-1}C_{n-1}\]