La prueba en la Wikipedia que se vincula en su intento anterior en esta pregunta es de hecho el uso de la misma definición de "determinante" que está utilizando, sólo que con diferente terminología.
Wikipedia utiliza la definición de determinante como la suma de todos los términos de la forma
$$\mathrm{sgn}(\tau)a_{1\tau(1)}\cdots a_{n\tau(n)}$$
donde $a_{ij}$ $(i,j)$ésima de la matriz, $\tau$ es una permutación de $\{1,2,\ldots,n\}$ (es decir, un bijection $\tau\colon\{1,2,\ldots,n\}\to\{1,2,\ldots,n\}$), y $\mathrm{sgn}(\tau)$ es el signo de $\tau$.
En su descripción, números de $i_1,\ldots,i_n$ debe ser, precisamente, los números de $\{1,2,\ldots,n\}$ en un poco de orden; dado un sumando
$$(-1)^ka_{1i_1}\cdots a_{ni_n}$$
como en su descripción, definimos $\tau\colon\{1,2,\ldots,n\}\to\{1,2,\ldots,n\}$$\tau(j) = i_j$. A continuación, $\tau$ es una permutación de $\{1,2,\ldots,n\}$, y podemos expresar el sumando como
$$(-1)^ka_{1\tau(1)}\cdots a_{n\tau(n)}.$$
Así que la única cosa que queda es mostrar que su factor de $(-1)^k$ es precisamente igual a $\mathrm{sgn}(\tau)$.
El signo de una permutación $\tau$ es igual a $1$ si $\tau$ puede ser escrito como el producto de un número par de transposiciones, y es igual a $-1$ si puede ser escrito como el producto de un número impar de transposiciones. Una "transposición" es precisamente un "intercambio de dos elementos".
Para cada "cambio" que se realiza para la lista de $i_1,i_2,\ldots,i_n$, considerar la transposición que corresponde a ese cambio. Realizar el intercambio es equivalente a componer $\tau$ con que la transposición. El hecho de que después de realizar la $k$ de estos intercambios de obtener la identidad indica que $\tau$ se puede escribir como un producto de $k$ transposiciones, y por lo tanto ese $\mathrm{sgn}(\tau)=(-1)^k$... exactamente el mismo que el factor de que usted tiene.
Por ejemplo, en su ejemplo,$a_{13}a_{21}a_{34}a_{42}$, la permutación $\tau$ asociado a este factor es, en discontinuo ciclo de la notación $\tau=(1,3,4,2)$, y en dos de la línea de la notación
$$\tau = \left(\begin{array}{cccc}
1&2&3&4\\
3&1&4&2
\end{array}\right).$$
El intercambio de la $1$ $3$ es equivalente a componer $\tau$ con la transposición de que los intercambios $1$ $3$ ($(1,3)$ en discontinuo ciclo de notación); a continuación, el intercambio de $2$$3$; y, finalmente, el intercambio de $3$$4$, con lo que consigue $(-1)^3 = -1$. Pero esto es lo mismo que decir que $\tau$ es la composición de los tres permutaciones que acabo de mencionar, es decir, que $(1,3,4,2) = (1,3)(2,3)(3,4)$, que es fácil de verificar. Así que la transposición de paridad impar, por lo tanto $\mathrm{sgn}(\tau)=-1$, lo que da exactamente el mismo resultado como el tuyo.
Desde las dos definiciones de determinante son de hecho idénticos, excepto por la terminología utilizada para describir los sumandos, la prueba en la Wikipedia es en el hecho de que "la prueba [...] con [su] definición".