Let Snβ be the set of permutations of [n] and let [n]n be the set of n-tuples in [n], that is, of maps f:[n]β[n]; both are nonempty and finite, by claims 3 and 4 of Finiteness of Cartesian Products, Tuple Sets, and Permutation Sets, and Snβ is a subset of [n]n. For fβ[n]n let Bfβ be the real nΓn matrix with (Bfβ)ijβ=Bf(i)jβ for all i,jβ[n], and put
Ξ±fβ=i=1βnβAif(i)β.
Expanding the product of the rows. Fix ΟβSnβ. By Product of Real Matrices,
(AB)iΟ(i)β=k=1βnβAikβBkΟ(i)β(iβ[n]),
so Generalized Distributivity: Expanding a Product of Finite Sums, applied with m=n to the family whose entry at (i,k) is AikβBkΟ(i)β, gives
i=1βnβ(AB)iΟ(i)β=fβ[n]nββΒ i=1βnβAif(i)βBf(i)Ο(i)β.
By claim 2 of Properties of Finite Products each inner product factors as
i=1βnβAif(i)βBf(i)Ο(i)β=Ξ±fβi=1βnβ(Bfβ)iΟ(i)β.
Multiplying by sgn(Ο) and using claim 4 of Properties of a Sum over a Finite Index Set,
sgn(Ο)i=1βnβ(AB)iΟ(i)β=fβ[n]nββΞ±fβsgn(Ο)i=1βnβ(Bfβ)iΟ(i)β.
Interchanging the two sums. Summing the last identity over ΟβSnβ and interchanging the two sums by claim 5 of Peeling, Splitting, and Interchange for Sums over a Finite Index Set, applied to the map on the product of Snβ and [n]n whose value at (Ο,f) is the summand above, and then taking the constant factor Ξ±fβ out of the inner sum by claim 4 of Properties of a Sum over a Finite Index Set,
det(AB)=fβ[n]nββΞ±fβΟβSnβββsgn(Ο)i=1βnβ(Bfβ)iΟ(i)β=fβ[n]nββΞ±fβdetBfβ,
the last equality by the form of the determinant recorded in Row Properties of the Determinant.
Discarding the non-injective terms. Let fβ[n]n and suppose f is not bijective. By claim 1 of An Injective Self-Map of a Finite Set is a Bijection it is then not injective, so there are p,qβ[n] with pξ =q and f(p)=f(q). Consequently (Bfβ)pjβ=Bf(p)jβ=Bf(q)jβ=(Bfβ)qjβ for every jβ[n], and claim 4 of Row Properties of the Determinant gives detBfβ=0, so the term at f vanishes by Zero Products and Elementary Identities in a Field. The remaining f are exactly the elements of Snβ, a nonempty subset of [n]n, so claim 4 of Peeling, Splitting, and Interchange for Sums over a Finite Index Set gives
det(AB)=fβSnβββΞ±fβdetBfβ.
Conclusion. For fβSnβ the matrix Bfβ is the row permutation of B by f, so claim 3 of Row Properties of the Determinant gives detBfβ=sgn(f)detB. Taking the constant factor detB out of the sum by claim 4 of Properties of a Sum over a Finite Index Set,
det(AB)=(fβSnβββsgn(f)i=1βnβAif(i)β)detB=detAΒ detB.