Podemos observar que el problema se puede expresar de esta manera: dado un $n \times n$ matriz $A$ con entradas sobre $\mathbb{F}_p$ (el campo finito de orden $p$ ), ¿cuál es la probabilidad de que la matriz tenga el determinante $1$ ?
En primer lugar, contemos cuántas matrices definidas de este modo tienen un determinante distinto de $0$ . Para que una matriz como esta tenga determinante no $0$ la primera columna puede ser cualquier vector distinto del vector nulo, por lo que tenemos $p^n - 1$ la segunda columna puede ser cualquier vector distinto de la primera columna multiplicado por un escalar, por lo que tenemos $p^n - p$ opciones...
Procediendo de esta manera tenemos que el número que buscamos es $$ r = (p^n - 1)(p^n - p)(p^n - p^2)\cdots (p^n - p^{n - 1})$$ En otras palabras, hemos encontrado $|GL(n, \mathbb{F}_p)| = r$ .
Ahora bien, puesto que $GL(n, \mathbb{F}_p)$ es un grupo bajo multiplicación y $\text{det} : GL(n, \mathbb{F}_p) \longrightarrow \mathbb{F}_p$ es un homomorfismo de grupo cuyo núcleo es el grupo de matrices con determinante $1$ y alcance $\mathbb{F}_p^*$ tenemos que el grupo de matrices con determinante $1$ es un subgrupo de $GL(n, \mathbb{F}_p)$ con índice $p - 1$ . Así que este subgrupo tiene $\frac{r}{p -1}$ elementos.
A partir de ahí tenemos $$q_n = \frac{r}{p^{n^2}(p -1)} = \frac{p^{n^2}\left(1 - \frac{1}{p^{n}}\right)\cdots \left(1 -\frac{1}{p}\right)}{p^{n^2}(p - 1)} = \frac{\left(1 - \frac{1}{p^n}\right)\cdots \left(1 - \frac{1}{p}\right)}{p - 1}$$ .
Lamentablemente no encuentro la forma de calcular $\lim_{n \to \infty} q_n$ Espero que alguien pueda ayudarnos en los comentarios o en otras respuestas.