Another exponent problem. $GL(n,q)$ is the group of invertible $n\times n$ matrices over the finite field $GF(q)$, where $q$ is a prime power. I am trying to figure out the exponent of this group.
Ideas:
We'd be doin real good if we could find the set of element orders $\{o(A):A\in GL(n,q)\}$, so I am looking for things to help me figure this out.
The order of $GL(n,q)$ is $$|GL(n,q)|=\prod_{k=0}^{n-1}\left(q^n-q^k\right)$$ The set of prime factors of this number is a little unpredictable so I'm not sure whether I can use this.
If $q$ is prime, the maximum order of an element in $GL(n,q)$ is at most $p^n-1$ (and this bound is strict I think). Also I guess the maximal power of $p$ that occurs as an order of this group is $p^{\lceil \log_pn\rceil}$. (source)
I am aware of this, which I believe implies that $\text{Exp}(GL(2,p^n))=\text{lcm}\{p(p^{2n}-1),p-1\}$. I am trying to think of how to extend it to the general case but I'm not sure.
And I'm stuck; can't think of anything else that might help.
Let us write $G=GL(n,q)$ for your group.
By the theorem that states that the so called Block Jordan Canonical Form of matrices over finite fields exist, we know that every matrix in $G$ is conjugate to one which is a block matrix with blocks of the form $$\begin{pmatrix}C_f\\I&C_f\\&I&C_f\\&&\ddots&\ddots\end{pmatrix} \tag{$\star$}$$ with $f$ an irreducible polynomial over the field $\mathbb F_q$ of some degree $d$, $C_d$ its companion matrix and $I$ the identity matrix of size $d$; see, for example, this notes.
It follows that the exponent of $G$ is the LCM of the orders of all matrices of this form with total size at most $n$.
So let $A$ be a matrix of the form $(\star)$, with $r\times r$ blocks, $f\in\mathbb F_q[X]$ an irreducible polynomial of degree $d$, $C_f$ its companion matrix and $I=I_d$ the $d\times d$ identity matrix. If we let $D$ be the «diagonal part» and $N$ be the «lower triangular part» (so that $D$ is a block diagonal matrix with $r$ blocks equal to $C_f$ along the diagonal and $N$ is a block matrix with non-zero blocks equal to $I_d$ exactly where these appear in $A$), we have $A=D+N$, $DN=DN$ and $N^r=0$. It follows immediately from this and the binomial formula that $$A^s=D^s+\binom{s}{1}D^{s-1}N+\binom{s}{2}D^{s-2}N^2+\cdots+\binom{s}{r-1}D^{s-r+1}N^{r-1}$$ for all $s\geq0$. It is easy to see what each term in this sum looks like, and using that we can see that $A^s=I_{dr}$ if and only if $D^s=I_{dr}$ and $\binom{s}{i}D^{s-i}=0$ for all $i\in\{1,\dots,r-1\}$. This happens iff $$\text{$C^s=I_d$ and $\binom{s}{i}=0$ for all $i\in\{1,\dots,r-1\}$.}$$ Of course, the order of $A$ is the least positive number $s$ which satisfies this condition. Plainly, it depends only on $f$ and on $r$, so we can call it $o(f,r)$. Moreover, it is clear that $o(f,r)\mid o(f,r+1)$, so the exponent of $G$ is $$LCM\{o(f,\lfloor n/\deg(f)\rfloor:f\in\mathbb F_q[X],\text{$f$ irreducible}, \deg(f)\leq n\}$$
If $n=2$, it is not very difficult to compute this: I get $LCM(p(p^{2n}-1),p(p^n-1))$, which happily equals with the formula I gave in a comment to the question.