This question is inspired by the cute answer to Order of general- and special linear groups over finite fields.
The formula $$ |\mathrm{GL}_n(\mathbb F_q)|=q^{\frac{n(n-1)}2}(q-1)(q^2-1)\cdots(q^n-1). $$ is obtained in that answer by simply but cleverly counting the number of linearly independent $n$-tuples of vectors in $\mathbb F_q^n$.
On the other hand, we know a lot about the structure of the group $\mathrm{GL}_n(\mathbb F_q)$: we can choose (in many ways) a maximal torus; let us take the one consisting of all invertible diagonal matrices, which is a subgroup of order $(q-1)^n$. We can next locate the unipotent radical of the corresponding Borel subgroup which in our case is the subgroup of all upper triangular matrices with $1$s along the main diagonal, thus has order $q^{\frac{n(n-1)}2}$. This accounts for $q^{\frac{n(n-1)}2}(q-1)^n$ elements.
How to account for the remaining factors $\frac{q^2-1}{q-1}$, $\frac{q^3-1}{q-1}$, ..., $\frac{q^n-1}{q-1}$? Do they also correspond to some subgroups that can be named, or maybe some explicitly describable conjugacy classes?