Existen $\displaystyle \sum_{i=0}^n \frac{N(n)}{N(i)N(n-i)}$ soluciones a su problema, donde $N(r)$ es el número de elementos en $GL_r(\mathbb{F}_p)$ .
Cualquier $X$ es invertible, y $X\in GL_n(\mathbb{F}_p)$ es una solución si y sólo si es conjugada en $GL_n(\mathbb{F}_p)$ a una matriz diagonal con la primera $i$ entradas iguales a $1$ y por último $n-i$ entradas iguales a $-1$ . Esto explica el número de sumandos. Ahora $g\in GL_n(\mathbb{F}_p)$ conmuta con dicha matriz diagonal si y sólo si preserva ambos eigenspaces. Por lo tanto, existen $\frac{N(n)}{N(i)N(n-i)}$ elementos en la clase de conjugación.
La fórmula puede hacerse aún más explícita sustituyendo la fórmula por $N(r)$ .