While proving some facts about matrix group operations on finite fields, I stumbled across the following question:
What is the order of the group of invertible $n\times n$ matrices over a finite field of prime order $p$?
The answer seems to be $$\left|\mathrm{GL}_n(\mathbb F_p)\right|=\prod_{k=1}^n\left(p^n-p^{k-1}\right)\text,$$ but I have not yet seen a satisfying proof: I took a look at this one, but found it to be somewhat sloppy. Is there a nice formal proof of this fact?