Jordan--Landau theorem for matrices over finite fields

Abstract: Given a positive integer $r$ and a prime power $q$, we estimate the probability that the characteristic polynomial $f_{A}(t)$ of a random matrix $A$ in $\mathrm{GL}_{n}(\mathbb{F}_{q})$ is square-free with $r$ (monic) irreducible factors when $n$ is large. We also estimate the analogous probability that $f_{A}(t)$ has $r$ irreducible factors counting with multiplicity. In either case, the main term $(\log n)^{r-1}((r-1)!n)^{-1}$ and the error term $O((\log n)^{r-2}n^{-1})$, whose implied constant only depends on $r$ but not on $q$ nor $n$, coincide with the probability that a random permutation on $n$ letters is a product of $r$ disjoint cycles. The main ingredient of our proof is a recursion argument due to S. D. Cohen, which was previously used to estimate the probability that a random degree $n$ monic polynomial in $\mathbb{F}_{q}[t]$ is square-free with $r$ irreducible factors and the analogous probability that the polynomial has $r$ irreducible factors counting with multiplicity. We obtain our result by carefully modifying Cohen's recursion argument in the matrix setting, using Reiner's theorem that counts the number of $n \times n$ matrices with a fixed characteristic polynomial over $\mathbb{F}_{q}$.

\section{Introduction}

\hspace{3mm} The purpose of this paper is to give some concrete arithmetic applications of the following fact: for any finitely many distinct $d_{1}, \dots, d_{r} \in \bZ_{\geq 1}$, the distribution of numbers of degree $d_{1}, \dots, d_{r}$ irreducible factors of the characteristic polynomial $f_{A}(t)$ of a matrix $A$ chosen uniformly at random from the set $\Mat_{n}(\bF_{q})$ of $n \times n$ matrices over $\bF_{q}$ converges to the distribution of numbers of length $d_{1}, \dots, d_{r}$ cycles of a random $\sigma \in S_{n}$, as $q$ goes to infinity.

\

\begin{prop}\label{main} Let $d_{1}, \dots, d_{r} \in \bZ_{\geq 1}$ be distinct and $k_{1}, \dots, k_{r} \in \bZ_{\geq 0}$ not necessarily distinct. We have

$$\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k_{j} \text{ degree } d_{j} \\
\text{irreducible factors for } 1 \leq j \leq r \\
\text{counting with multiplicity}
\end{array}\right)
= \Prob_{\sigma \in S_{n}}\left(\begin{array}{c}
\sigma \text{ has } k_{j} \text{ disjoint} \\
d_{j}\text{-cycles for } 1 \leq j \leq r
\end{array}\right).$$

\

where $f_{A}(t)$ is the characteristic polynomial of $A \in \Mat_{n}(\bF_{q})$.\footnote{The event on the right-hand side means that when we write $\sigma$ as a unique product of disjoint cycles, there are $k_{j}$ cycles of length $d_{j}$ in the product for $1 \leq j \leq r$. The special case $r = 1$ is indirectly mentioned in ``Conclusion'' section of \cite{Sto}.}
\end{prop}

\

\begin{rmk} After the previous version of this paper, we have been noted that Proposition \ref{main} is known due to Hansen and Schumutz (\cite{HS}, Theorem 1.1). Their result is for $\GL_{n}(\bF_{q})$ rather than $\Mat_{n}(\bF_{q})$, but these two are easily interchangeable. Indeed, the last line of their proof (but not the statement), together with a well-known fact regarding some statistics of degree $n$ monic polynomials in $\bF_{q}[t]$ for large $q$, implies Proposition \ref{main}. The proof of Hansen and Schumutz uses Stong's work \cite{Sto}, on which our proof also depends, and they even compute some asymptotic errors with respect to $q$ and $r$. Since our work only considers the limit in $q$, our proof is simpler, so we decided to keep it in the paper.

\

\hspace{3mm} In the following subsections, we discuss two applications of Proposition \ref{main}. The first one, Theorem \ref{p-adic}, considers the distribution of the cokernel of a random $n \times n$ matrix over $\bZ_{p}$, the ring of $p$-adic integers, with respect to the Haar (probability) measure on $\Mat_{n}(\bZ_{p}) = {\bZ_{p}}^{n^{2}}$, for large $p$ and $n$. According to our best knowledge and the consultations with the experts we have talked to, Theorem \ref{p-adic} is new. The second one, Theorem \ref{JL}, is a matrix version of Landau's theorem in number theory that estimates the number of irreducible factors of a random characteristic polynomial for large $n$ when $q \ra \infty$. Just as Proposition \ref{main} can be deduced from the work of Hansen and Schumutz \cite{HS}, Theorem \ref{JL} can also be deduced from their work if we apply another result of S. Cohen in \cite{Coh}. The interesting difference in our exposition is that we instead used Jordan's result in \cite{Jor} about the statistics of a random permutation in $S_{n}$ for large $n$ rather than the statistics of degree $n$ monic polynomials in $\bF_{q}[t]$ for large $q$ and $n$ in order to deduce Theorem \ref{JL}.
\end{rmk}

\

\subsection{The distribution of cokernels of Haar $\bZ_{p}$-random matrices with large $p$} Let $\bZ_{p}$ be the ring of $p$-adic integers, namely the inverse limit of the system of projection maps $\cdots \tra \bZ/(p^{3}) \tra \bZ/(p^{2}) \tra \bZ/(p)$.

\

\begin{thm}\label{p-adic} Fix any distinct $d_{1}, \dots, d_{r} \in \bZ_{\geq 1}$. For any $d \in \bZ_{\geq 1}$ denote by $\mc{M}(d, p)$ the set of all monic irreducible polynomials of degree $d$ in $\bF_{p}[t]$. Then the limit

$$\lim_{p \ra \infty} \Prob_{A \in \Mat_{n}(\bZ_{p})}
\left(\begin{array}{c}
\coker(P(A)) = 0 \\
\text{ for all } \bar{P} \in \bigsqcup_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right)$$

\

exists, where for $P(A)$, we use any lift of $\bar{P}(t) \in \bF_{p}[t]$ in $\bZ_{p}[t]$ under the reduction modulo $p$ and the probability is according to the Haar measure on $\Mat_{n}(\bZ_{p})$. Moreover, we have

$$\lim_{n \ra \infty} \lim_{p \ra \infty} \Prob_{A \in \Mat_{n}(\bZ_{p})}
\left(\begin{array}{c}
\coker(P(A)) = 0 \\
\text{ for all } \bar{P} \in \bigsqcup_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right)
= e^{-1/d_{1}}  \cdots e^{-1/d_{r}}.$$
\end{thm}

\

\begin{rmk} Taking the two limits in the particular order matters in Theorem \ref{p-adic}. For instance, we are unsure whether we can change the order of these limits, which may lead to an interesting consequence regarding the Cohen-Lenstra measure introduced in \cite{CL}. It is worth noting that our proof can be easily modified so that one may replace $\bZ_{p}$ with $\bF_{q}\llb t \rrb$, reduce modulo $(t)$ in place of $p$, and let $q \ra \infty$ instead of $p \ra \infty$ in Theorem \ref{p-adic}, although we will just focus on $\bZ_{p}$. The reader does not need to have much expertise of the Haar measure on $\Mat_{n}(\bZ_{p})$ to follow our proof. The only property of the Haar measure on $\Mat_{n}(\bZ_{p})$ we will use is that each fiber under the mod $p$ projection map $\Mat_{n}(\bZ_{p}) \tra \Mat_{n}(\bF_{p})$ has the same measure, which is necessarily $1/p^{n^{2}}$.
\end{rmk}

\

\hspace{3mm} We conjecture that the large $p$ limit in Theorem \ref{p-adic} has to do with independent Poisson random variables with means $1/d_{1}, \dots, 1/d_{r}$ in more generality. However, we have little empirical evidence, so it would be interesting to see any more examples that either support or disprove this conjecture.

\

\begin{conj} Fix any distinct $d_{1}, \dots, d_{r} \in \bZ_{\geq 1}$ and not necessarily distinct $k_{1}, \dots, k_{r} \in \bZ_{\geq 0}$. Then the limit

$$\lim_{p \ra \infty} \Prob_{A \in \Mat_{n}(\bZ_{p})}
\left(\begin{array}{c}
|\coker(P_{j}(A))| = p^{d_{j}k_{j}} \\
\text{ for all } (\bar{P}_{1}, \dots, \bar{P}_{r}) \in \prod_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right)$$

\

must exist. Moreover, we must have

$$\lim_{n \ra \infty}\lim_{p \ra \infty} \Prob_{A \in \Mat_{n}(\bZ_{p})}
\left(\begin{array}{c}
|\coker(P_{j}(A))| = p^{d_{j}k_{j}} \\
\text{ for all } (\bar{P}_{1}, \dots, \bar{P}_{r}) \in \prod_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right)
= \frac{e^{-1/d_{1}}(1/d_{1})^{k_{1}}}{k_{1}!}  \cdots \frac{e^{-1/d_{r}}(1/d_{r})^{k_{r}}}{k_{r}!}.$$
\end{conj}

\

\subsection{Jordan-Landau theorem for matrices over finite fields} It is a folklore that the cycle decomposition of a random permutation is analogous to the prime decomposition of a random positive integer. For instance, Granville (Section 2.1 of \cite{Gra}) notes that given $k \in \bZ_{\geq 1}$, a theorem of Jordan (p.161 of \cite{Jor}) says

$$\Prob_{\sigma \in S_{n}}(\sigma \text{ has } k \text{ disjoint cycles}) \sim \frac{(\log n)^{k-1}}{(k-1)! n}$$

\

for large $n \in \bZ_{\geq 1}$, while a theorem of Landau (p.211 of \cite{Lan} or Theorem 437 on p.491 of \cite{HW}) says

\begin{align*}
\Prob_{1 \leq N \leq x}\left(\begin{array}{c}
N \text{ has } k \text{ distinct} \\
\text{prime factors}
\end{array}\right) &\sim \Prob_{1 \leq N \leq x}\left(\begin{array}{c}
N \text{ has } k \text{ prime factors} \\
\text{counting with multiplicity}
\end{array}\right)\\
&\sim \frac{(\log \log x)^{k-1}}{(k-1)!\log x}
\end{align*}

\

for large $x \in \mathbb{R}_{\geq 1}$. The first probability is given uniformly at random from the set $S_{n}$ of permutations on $n$ letters. The second probability is given by choosing the integer $N$ uniformly at random from the set $\{1, 2, \dots, \lfloor x \rfloor\}$. We wrote ``$f(x) \sim g(x)$ for large $x$'' to mean that $f(x)/g(x) \ra 1$ as $x \ra \infty$. This is a generalization of the Prime Number Theorem (i.e., the case $k = 1$).

\

\hspace{3mm} By the Prime Number Theorem, for large $x$, we may expect one prime in every interval of length $\log x$. For large $n$, we expect one cycle in $n$ permutations (since there are $(n-1)!$ $n$-cycles while there are far less cycles of shorter lengths), so Jordan's result is analogous to Landau's result. We take this analogy further. In particular, summing the identity given in Proposition \ref{main} over all tuples $(d_{1}, \dots, d_{r}) \in (\bZ_{\geq 1})^{r}$ and $(k_{1}, \dots, k_{r}) \in (\bZ_{\geq 0})^{r}$ such that

\begin{itemize}
	\item $r \in \bZ_{\geq 0}$,
	\item $n = \sum_{j=1}^{r}k_{j}d_{j}$, and
	\item $k = \sum_{j=1}^{r} k_{j}$,
\end{itemize}

\

we have

$$\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right)
= \Prob_{\sigma \in S_{n}}(\sigma \text{ has } k \text{ disjoint cycles}),$$

\

which is asymptotically $(\log n)^{k-1}/((k-1)!n)$ for large $n$, as noted above. Using the facts that almost all matrices $A \in \Mat_{n}(\bF_{q})$ have square-free characteristic polynomials and almost all monic polynomials of degree $n$ are square-free when $q \ra \infty$ (see Section \ref{Landau}), we may obtain the following.

\

\begin{thm}[Jordan-Landau theorem for $\bF_{q}$-matrices]\label{JL} Given $k \in \bZ_{\geq 1}$, for large $n \in \mathbb{Z}_{\geq 1}$, we have

\begin{align*}
\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ distinct} \\
\text{irreducible factors}
\end{array}\right) &= \lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right) \\
&\sim \frac{(\log n)^{k-1}}{(k-1)! n}.
\end{align*}
\end{thm}

\

\begin{rmk} There is an analogous result of Theorem \ref{JL}, where the sample space is the set $\bA^{n}(\bF_{q})$ of all degree $n$ monic polynomials in $\bF_{q}[t]$ due to S. Cohen \cite{Coh}. It says given $k \in \bZ_{\geq 1}$, for large $n \in \mathbb{Z}_{\geq 1}$, we have

\begin{align*}
\lim_{q \ra \infty}\Prob_{f \in \bA^{n}(\bF_{q})}
\left(\begin{array}{c}
f(t) \text{ has } k \text{ distinct} \\
\text{irreducible factors}
\end{array}\right) &\sim \lim_{q \ra \infty}\Prob_{f \in \bA^{n}(\bF_{q})}
\left(\begin{array}{c}
f(t) \text{ has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right) \\
&\sim \frac{(\log n)^{k-1}}{(k-1)! n}.
\end{align*}
\end{rmk}

\

With these two results in mind, we note that the push-forward measure $\Mat_{n}(\bF_{q}) = \bA^{n^{2}}(\bF_{q}) \tra \bA^{n}(\bF_{q})$ given by taking the characteristic polynomial of a matrix does not change the probability of the event we compute in $\bA^{n}(\bF_{q})$, when we take $q \ra \infty$. This can be also taken as another example to be added in the list of similarly behaving combinatorial examples in \cite{ABT}, a survey paper by Arratia, Barbour, and Tavar\'e.

\

\subsection{Fulman's generalization and further directions} Recently, we have also realized that Fulman generalized Proposition \ref{main} to certain linear algebraic groups $G$ over $\bF_{q}$, where $G = \GL_{n}$ is the desired special case for Proposition \ref{main} (\cite{Ful99}, Theorem 15), where $S_{n}$ occurs as the Weyl group of $G(\ol{\bF_{q}}) = \GL_{n}(\ol{\bF_{q}})$. Fulman's argument needs to assume that the probability that a random element in $G(\bF_{q}) \subset G(\ol{\bF_{q}})$ is regular semisimple converges to $1$, as $q \ra \infty$, which we previously thought as a mere heuristics. There is a statement in Fulman's paper (\cite{Ful99}, Theorem 5 and Theorem 7) that proves a similar property related to this hypothesis for $G = \GL_{n}$, but it takes $n \ra \infty$ before we can let $q \ra \infty$. However, for $g \in \GL_{n}(\bF_{q})$ saying that $g$ is regular semisimple turns out to mean that the characteristic polynomial of $g$ is square-free in $\bF_{q}[t]$, so using a similar argument as in Lemma \ref{sfmat}, we see that $G = \GL_{n}$ satisfies Fulman's hypothesis because the number of degree $n$ monic square-free polynomials $f(t) \in \bF_{q}[t]$ such that $f(0) \neq 0$ is precisely

$$q^{n} - 2q^{n-1} + 2q^{n-2} - \cdots + (-1)^{n-1}2q + (-1)^{n},$$

\

which replaces Lemma \ref{sf} in the argument. We are not sure whether the above count is explicitly written in literature, but this is certainly well-known (e.g., more general phenomena are dealt in various works such as \cite{BE}, \cite{Kim}, and \cite{VW}). More specifically, one can obtain this count by counting the $\bF_{q}$-points of the unordered configuration space $\Conf^{n}(\bA^{1} \sm \{0\})$ of $n$ points on the affine line minus the origin over $\bF_{q}$. The generating function for such count $|\Conf^{n}(\bA^{1} \sm \{0\})(\bF_{q})|$ is given by

$$\sum_{n=0}^{\infty} |\Conf^{n}(\bA^{1} \sm \{0\})(\bF_{q})| t^{n} = \frac{Z_{\bA^{1} \setminus \{0\}}(t)}{Z_{\bA^{1} \setminus \{0\}}(t^{2})} = \frac{1 - qt^{2}}{(1 + t)(1 - qt)},$$

\

where $Z_{X}(t)$ means the zeta series of an algebraic variety over $\bF_{q}$, and expanding the right-hand side, one may obtain the desired count. We are not sure whether this hypothesis may or may not be satisfied by other linear algebraic groups $G$ over $\bF_{q}$. The reader must note that the reasons that this works for $G = \GL_{n}$ are as follows:

\be
	\item we are aware of the size of each fiber of the map $\Mat_{n}(\bF_{q}) = \bA^{n^{2}}(\bF_{q}) \rightarrow \bA^{n}(\bF_{q})$ given by taking the characteristic polynomials and
	\item the preimage of the set of degree $n$ monic polynomials with nonzero constant terms under this map is precisely $\GL_{n}(\bF_{q})$.
\ee

Nevertheless, it would be extremely interesting if one can identify which algebraic groups $G$, other than $\GL_{n}$, over $\bF_{q}$ produce the sets $G(\bF_{q})$ of $\bF_{q}$-points that can replace $\Mat_{n}(\bF_{q})$ with the Weyl groups of $G(\ol{\bF_{q}})$, replacing $S_{n}$ in Proposition \ref{main}. If such algebraic groups are given in a sequence $G_{n}$ in $n \in \mathbb{Z}_{\geq 1}$, then we can hope that studying the corresponding sequence $W_{n}$ of Weyl groups asymptotically in $n$ may produce analogous statements for Theorem \ref{p-adic} and Theorem \ref{JL}.

\

\subsection{Organization for the rest} In Section \ref{Landau}, we will show that Proposition \ref{main} implies Theorem \ref{JL}. In Section \ref{Haar}, we will show that Proposition \ref{main} implies Theorem \ref{p-adic}. We will provide a crucial formula due to Stong in Section \ref{Stong}, but our proof is slightly different from the original one, as we use direct counting of Young diagrams. Stong's formula will be used in proving Proposition \ref{main} in Section \ref{mainproof}. Finally, in Section \ref{permproof}, we will provide another proof of Lemma \ref{perm}, an influential result of Shepp and Lloyd, which states that the number of length $d$ cycles of a random permutation in $S_{n}$ asymptotically follows the Poisson distribution with mean $1/d$ when $n$ is large, and the probability is given independently for finitely many different choices of $d$. It is interesting that our proof is entirely combinatorial, while the original proof was probabilistic. Since the result is quite popular in literature, our alternative proof of Lemma \ref{perm} may be known to experts, although we were unable to locate it.

\

\subsection{Acknowledgment} We thank Yifeng Huang, Nathan Kaplan, and Jeff Lagarias for helpful conversations. G. Cheong was supported by NSF grant DMS-1162181 and the Korea Institute for Advanced Study for his visits to the institution regarding this research. M. Yu was supported by a KIAS Individual Grant (SP075201) via the Center for Mathematical Challenges at Korea Institute for Advanced Study. He was also supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2020R1C1C1A01007604).

\

\section{Proposition \ref{main} implies Theorem \ref{JL}}\label{Landau}

\hspace{3mm} We already know that Proposition \ref{main} implies that

$$\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right) \sim \frac{(\log n)^{k-1}}{(k-1)! n}$$

\

for large $n$ in the statement of Theorem \ref{JL}. We need to justify the statement about distinct irreducible factors. We need two lemmas, well-known to the experts, to provide a complete exposition.

\

\begin{lem}\label{sf} Let $\bA^{n}(\bF_{q})$ be the set of all monic polynomials of degree $n$ in $\bF_{q}[t]$. For $n \geq 2$, we have

$$|\{f \in \bA^{n}(\bF_{q}): f(t) \text{ is square-free in } \bF_{q}[t]\}| = q^{n}(1 - q^{-1}).$$

\

In particular, we have

$$\lim_{q \ra \infty}\Prob_{f \in \bA^{n}(\bF_{q})}(f(t) \text{ is square-free in } \bF_{q}[t]) = 1.$$
\end{lem}

\

\begin{proof} This is a well-known fact in literature. For a proof, see Lemma 3.4 of \cite{CWZ}. A more combinatorial proof using partitions is also available (e.g., Proposition 5.9 in \cite{VW} with $a = 2$).
\end{proof}

\

\begin{lem}\label{sfmat} For any $n \in \bZ_{\geq 1}$, we have

$$\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}(f_{A}(t) \text{ is square-free in } \bF_{q}[t]) = 1.$$
\end{lem}

\

\begin{proof} Fix any monic square-free polynomial

$$f(t) = P_{1}(t) \cdots P_{r}(t),$$

\

where $P_{i}(t)$ are distinct monic irreducible polynomials in $\bF_{q}[t]$. By Theorem 2 in \cite{Rei}, the number of $A \in \Mat_{n}(\bF_{q})$ such that $f_{A}(t) = f(t)$ is $$q^{n^{2} - n} \cdot \frac{(1 - q^{-1})(1 - q^{-2}) \cdots (1 - q^{-n})}{(1 - q^{-\deg(P_{1})}) (1 - q^{-\deg(P_{2})}) \cdots (1 - q^{-\deg(P_{r})})} \geq q^{n^{2} - n} \cdot \frac{(1 - q^{-1})(1 - q^{-2}) \cdots (1 - q^{-n})}{(1 - q^{-n})^{r}}.$$

\

For $n = 1$, the map $A \mapsto f_{A}(t)$ is bijective, and thus letting $q \ra \infty$, we are done. Foe $n \geq 2$, by Lemma \ref{sf}, we have

$$\Prob_{A \in \Mat_{n}(\bF_{q})}(f_{A}(t) \text{ is square-free in } \bF_{q}[t]) \geq (1 - q^{-1}) \cdot \frac{(1 - q^{-1})(1 - q^{-2}) \cdots (1 - q^{-n})}{(1 - q^{-n})^{r}},$$

\

so letting $q \ra \infty$ gives the result.
\end{proof}

\

\begin{proof}[Proof that Proposition \ref{main} implies Theorem \ref{JL}] Again, we only need to show the first identity in the statement. By Lemma \ref{sfmat}, we have

\begin{align*}
\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ distinct} \\
\text{irreducible factors}
\end{array}\right) &= \lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ is square-free in } \bF_{q}[t] \text{ and}\\
\text{has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right) \\
&= \lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}
\left(\begin{array}{c}
f_{A}(t) \text{ has } k \text{ irreducible factors} \\
\text{counting with multiplicity}
\end{array}\right),
\end{align*}

\

because the number of monic non-square-free polynomials of degree $n$ is negligible as $q \ra \infty$.
\end{proof}

\

\section{Proposition \ref{main} implies Theorem \ref{p-adic}}\label{Haar}

\subsection{Useful lemmas} The following result, due to Shepp and Lloyd, will be crucial in applying Proposition \ref{main} to obtain Theorem \ref{p-adic}. In \cite{SL}, Shepp and Lloyd obtained the result by showing that the the characteristic function of the distribution of the numbers of cycles of fixed lengths $d_{1}, \dots, d_{r}$ of a random permutation converges to the characteristic function of the distribution given by the independent Poisson random variables with the means $1/d_{1}, \dots, 1/d_{r}$ and then applying L\'evy's theorem. We will provide our own combinatorial proof of this result in Section \ref{permproof}.

\

\begin{notn}
For any permutation $\sigma \in S_{n}$, we denote by $m_{d}(\sigma)$ the number of $d$-cycles in the cycle decomposition of $\sigma$.
\end{notn}

\

\begin{lem}[cf. p.342 in \cite{SL}]\label{perm} Fix $r \in \bZ_{\geq 0}$. Given distinct $d_{1}, \dots, d_{r} \in \bZ_{\geq 1}$ and not necessarily distinct $k_{1}, \dots, k_{r} \in \bZ_{\geq 0}$, we have

$$\lim_{n \ra \infty}\Prob_{\sigma \in S_{n}}(m_{d_{j}}(\sigma) = k_{j} \text{ for } 1 \leq j \leq r) = \frac{e^{-1/d_{1}}(1/d_{1})^{k_{1}}}{k_{1}!} \cdots \frac{e^{-1/d_{r}}(1/d_{r})^{k_{r}}}{k_{r}!},$$

\

which means that the number of cycles of length $d_{1}, \dots, d_{r}$ of a random permutation of $n$ letters are asymptotically given by independent Poisson random variables with means $1/d_{1}, \dots, 1/d_{r}$ when $n$ is large.
\end{lem}

\

\hspace{3mm} We need the following lemma to deduce Theorem \ref{p-adic} from Proposition \ref{main}. This lemma contains some ideas used for proving Theorem C of \cite{CH}, but the proof here is simpler, since we only need a special case of their argument.

\

\begin{lem}\label{red} Let $P_{1}, \dots, P_{r} \in \bZ_{p}[t]$ be any monic polynomials whose images in $\bF_{p}[t]$ under the mod $p$ reduction are distinct irreducible polynomials. Then

$$\Prob_{A \in \Mat_{n}(\bZ_{p})}\left(\begin{array}{c}
\coker(P_{j}(A)) = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right) = \Prob_{A \in \Mat_{n}(\bF_{p})}\left(\begin{array}{c}
A[\bar{P}_{j}^{\infty}] = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right),$$

\

where $A[\bar{P}_{j}^{\infty}]$ means the $\bar{P}_{j}$-part of the $\bF_{p}[t]$-module structure given by the matrix action $A \acts \bF_{p}^{n}$.
\end{lem}

\begin{proof} First, note that for any $A \in \Mat_{n}(\bF_{p})$, we have $\dim_{\bF_{p}}\ker(\bar{P}_{j}(A)) = \dim_{\bF_{p}}\coker(\bar{P}_{j}(A))$ and $\ker(\bar{P}_{j}(A)) = A[\bar{P}_{j}^{\infty}]/\bar{P}_{j}(t)A[\bar{P}_{j}^{\infty}]$, as $\bF_{p}[t]$-modules where the action of $t$ is given by multiplying the matrix $A$ on the left. Hence, the following are equivalent:

\

\begin{itemize}
	\item $A[\bar{P}_{j}^{\infty}] = 0$;
	\item $\ker(\bar{P}_{j}(A)) = 0$;
	\item $\coker(\bar{P}_{j}(A)) = 0$,
\end{itemize}

\

so it is enough to show that

$$\Prob_{A \in \Mat_{n}(\bZ_{p})}\left(\begin{array}{c}
\coker(P_{j}(A)) = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right) = \Prob_{A \in \Mat_{n}(\bF_{p})}\left(\begin{array}{c}
\coker(\bar{P}_{j}(A)) = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right).$$

\

For any $B \in \Mat_{n}(\bZ_{p})$, the right exactness of $(-) \otimes_{\bZ_{p}} \bF_{p}$ implies that $$\coker(B) \otimes_{\bZ_{p}} \bF_{p} \simeq \coker(\bar{B}),$$ where $\bar{B} \in \Mat_{n}(\bF_{p})$ is the reduction of $B$ mod $p$. Hence, by Nakayama's lemma, we have $\coker(B) = 0$ if and only if $\coker(\bar{B}) = 0$. Thus, if we consider the mod $p$ reduction map $\pi : \Mat_{n}(\bZ_{p}) \tra \Mat_{n}(\bF_{p})$, we have $$\pi^{-1}\left\{\begin{array}{c}
A \in \Mat_{n}(\bF_{p}) : \coker(\bar{P}_{j}(A)) = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right\} = \left\{\begin{array}{c}
A \in \Mat_{n}(\bZ_{p}) : \coker(P_{j}(A)) = 0 \\
\text{ for } 1 \leq j \leq r
\end{array}\right\}.$$

\

Since each fiber under $\pi$ has the same measure $1/p^{n^{2}}$, this finishes the proof.
\end{proof}

\

\subsection{Proof of Theorem \ref{p-adic} given Proposition \ref{main}} We now prove Theorem \ref{p-adic}, assuming Proposition \ref{main}. Proposition \ref{main} will be proven in Section \ref{mainproof}.

\

\begin{proof}[Proof of Theorem \ref{p-adic} given Proposition \ref{main}] Let $\mc{M}(q, d)$ be the set of degree $d$ monic irreducible polynomials in $\bF_{q}[t]$ so that $M(q, d) = |\mc{M}(q, d)|$. By Lemma \ref{red}, we have

\begin{align*}\Prob_{A \in \Mat_{n}(\bZ_{p})}
\left(\begin{array}{c}
\coker(P(A)) = 0 \\
\text{ for all } P \in \bigsqcup_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right) &= \Prob_{A \in \Mat_{n}(\bF_{p})}
\left(\begin{array}{c}
A[\bar{P}_{j}^{\infty}] = 0 \\
\text{ for all } P \in \bigsqcup_{j=1}^{r}\mc{M}(p, d_{j})
\end{array}\right) \\
&= \Prob_{A \in \Mat_{n}(\bF_{p})}
\left(\begin{array}{c}
f_{A}(t) \text{ has no irreducible factors of} \\
\text{degrees } d_{1}, \dots, d_{r}
\end{array}\right).
\end{align*}

\

Hence, we may apply Proposition \ref{main} and Lemma \ref{perm} to finish the proof.
\end{proof}

\

\section{Stong's formula}\label{Stong}

\subsection{Set-up} Let $R$ be a PID, and say $\mf{m}$ is a maximal ideal of $R$ such that $R/\mf{m} = \bF_{q}$. Denote by $\mc{P}$ the set of all partitions including the empty partition $\es$. Then any finite length (or equivalently, finite size) $\mf{m}^{\infty}$-torsion $R$-module is isomorphic to

$$H_{\mf{m},\ld}^{R} := R/\mf{m}^{\ld_{1}} \op \cdots \op R/\mf{m}^{\ld_{l}},$$

\

for a unique partition $\ld = (\ld_{1}, \dots, \ld_{l}) \in \mc{P}$. We will always assume $\ld_{1} \geq \cdots \geq \ld_{l} > 0$, and we write $|\ld| := \ld_{1} + \cdots + \ld_{l}$. The case $l = 0$ will correspond to the empty partition $\ld = \es$. Note that

$$|H_{\mf{m},\ld}^{R}| = q^{\ld_{1} + \cdots + \ld_{l}} = q^{|\ld|}.$$

\

We will write $H_{\mf{m},\ld} := H_{\mf{m},\ld}^{R}$ if the meaning of $R$ is evident from the context. Denote by $m_{d}(\ld)$ the number of parts of size $d$ in $\ld$, and define

$$n(\ld) := 0 \cdot \ld_{1} + 1 \cdot \ld_{2} + 2 \cdot \ld_{3} + \cdots + (l-1) \cdot \ld_{l}.$$

\

\hspace{3mm} The number $|\Aut_{R}(H_{\mf{m},\ld}^{R})|$ of automorphisms of $H_{\mf{m},\ld}^{R}$ can be computed by noting the fact that

$$|\Aut_{R}(H_{\mf{m},\ld}^{R})| = |\Aut_{R_{\mf{m}}}(H_{\mf{m}R_{\mf{m}},\ld}^{R_{\mf{m}}})|$$

\

and the following formula due to Macdonald.

\

\begin{lem}[(1.6) on p.181 in \cite{Mac}]\label{Mac} Let $(R, \mf{m})$ be a DVR (discrete valuation ring) with the finite residue field $R/\mf{m} = \bF_{q}$. Then we have

$$|\Aut_{R}(H_{\mf{m},\ld})| = q^{|\ld|+2n(\ld)} \prod_{d = 1}^{\infty}\prod_{i=1}^{m_{d}(\ld)}(1 - q^{-i}).$$
\end{lem}

\

\hspace{3mm} To our best knowledge, the following result is due to Stong in \cite{Sto}, but we rephrase Stong's result in a more general setting and provide a slightly different proof. We will essentially go through some key ideas in Lemma 6 in \cite{Sto}, which relies on the Fine-Herstein theorem (i.e., the number of $n \times n$ nilpotent matrices over $\bF_{q}$ is $q^{n(n-1)}$), but unlike the reference, we will avoid the use of M\"obius inversion along with exponentiation, logarithm, differentiation and integration of power series by counting partitions (i.e., Young diagrams) instead.

\

\begin{lem}[cf. Proposition 19 in \cite{Sto}]\label{key} Let $R$ be any Dedekind domain with a maximal ideal $\mf{m}$ such that $R/\mf{m} = \bF_{q}$. Then

$$\sum_{M \in \Mod_{R_{\mf{m}}}^{<\infty}} \frac{u^{q\dim_{\bF_q}(M)}}{|\Aut_{R}(M)|} = \prod_{i=1}^{\infty}\frac{1}{1 - q^{-i}u^{q}},$$

\

where $u$ is a complex variable with $|u| < q^{1/q}$. Equivalently, we have

$$\sum_{\ld \in \mc{P}} \frac{y^{|\ld|}}{|\Aut_{R}(H_{\mf{m},\ld})|} = \prod_{i=1}^{\infty}\frac{1}{1 - q^{-i}y},$$

\

where $y$ is a complex variable with $|y| < q$.
\end{lem}

\begin{proof} As explained in the beginning of this section, finite $\mf{m}^{\infty}$-torsion $R$-modules are parametrized by partitions. Hence, taking $y = u^{q}$, we see that the two given statements are equivalent. We will prove the second statement. Applying Lemma \ref{Mac}, this reduces our problem to the case where $R = \bF_{q}[t]$ and $\mf{m} = (t)$ by replacing $R_{\mf{m}}$ with $\bF_{q}[t]_{(t)}$ or $\bF_{q}\llb t \rrb$. This reduction lets us rewrite the left-hand side of the desired identity as

$$\sum_{n = 0}^{\infty}\sum_{\ld \vdash n}\frac{y^{n}}{|\Stab_{\GL_{n}(\bF_{q})}(J_{0,\ld})|},$$

\

where $J_{0,\ld}$ is the Jordan canonical form given by the $0$-Jordan blocks of whose sizes are equal to the parts of the partition $\ld$, and the action $\GL_{n}(\bF_{q}) \acts \Mat_{n}(\bF_{q})$ is given by the conjugation. By the orbit-stabilizer theorem, we have

$$\frac{1}{|\Stab_{\GL_{n}(\bF_{q})}(J_{0,\ld})|} = \frac{|\GL_{n}(\bF_{q}) \cdot J_{0,\ld}|}{|\GL_{n}(\bF_{q})|}.$$

\

Recall that $n \times n$ matrices all of whose eigenvalues are $0$ are precisely nilpotent matrices. Thus, we have

\begin{align*}
\sum_{n = 0}^{\infty}\sum_{\ld \vdash n}\frac{y^{n}}{|\Stab_{\GL_{n}(\bF_{q})}(J_{0,\ld})|} &= \sum_{n = 0}^{\infty}\frac{|\Nil_{n}(\bF_{q})|y^{n}}{|\GL_{n}(\bF_{q})|} \\
&= 1 + \sum_{n = 1}^{\infty}\frac{q^{n(n-1)}y^{n}}{(q^{n} - 1)(q^{n} - q) \cdots (q^{n} - q^{n-1})} \\
&= 1 + \sum_{n = 1}^{\infty}\frac{(q^{-1}y)^{n}}{(1 - q^{-1})(1 - q^{-2}) \cdots (1 - q^{-n})} \\
&= \sum_{n = 0}^{\infty}(q^{-1}y)^{n} \left(\sum_{j_{1}=0}^{\infty} q^{-j_{1}}\right) \left(\sum_{j_{2}=0}^{\infty} q^{-2j_{2}}\right) \cdots \left(\sum_{j_{n}=0}^{\infty} q^{-nj_{n}}\right) \\
&= \sum_{n = 0}^{\infty}(q^{-1}y)^{n} \sum_{j_{1}=0}^{\infty} \sum_{j_{2}=0}^{\infty} \cdots \sum_{j_{n}=0}^{\infty} q^{-j_{1} - 2j_{2} - \cdots - nj_{n}} \\
&= \prod_{i=0}^{\infty} (1 + q^{-i}(q^{-1}y) + q^{-2i}(q^{-1}y)^{2} + \cdots) \\
&= \prod_{i=0}^{\infty} \frac{1}{1 - q^{-(i+1)}y} \\
&= \prod_{i=1}^{\infty} \frac{1}{1 - q^{-i}y},
\end{align*}

\

where for the second equality we used the elementary identity $|\GL_{n}(\bF_{q})| = (q^{n} - 1)(q^{n} - q) \cdots (q^{n} - q^{n-1})$, and Fine-Herstein theorem (e.g., Theorem 1 of \cite{FH}), which gives the number $|\Nil_{n}(\bF_{q})|$ of nilpotent matrices in $\Mat_{n}(\bF_{q})$:

$$|\Nil_{n}(\bF_{q})| = q^{n(n-1)}.$$

\

For the sixth equality, we have used the fact that the coefficient of $Y^{n}$ of the following product

$$\prod_{i=0}^{\infty}(1 + X^{i}Y + X^{2i}Y^{2} + \cdots)$$

\

is equal to

$$\sum_{j_{1}=0}^{\infty} \sum_{j_{2}=0}^{\infty} \cdots \sum_{j_{n}=0}^{\infty} X^{j_{1} + 2j_{2} + \cdots + nj_{n}}.$$

\

where $X, Y$ are considered to be complex numbers varying within the open unit disk centered at $0$ (i.e., $|X|, |Y| < 1$) so that we can take $X = q^{-1}$ and $Y = q^{-1}y$ with $|y| < q$ in our proof for the sixth equality in the chain of equalities above. Indeed, when we expand the given product, we have

\begin{align*}
\prod_{i=0}^{\infty}(1 + X^{i}Y + X^{2i}Y^{2} + \cdots) &= \sum_{
m_{0}, m_{1}, m_{2}, \dots \geq 0} X^{0 \cdot m_{0} + 1 \cdot m_{1} + 2 \cdot m_{2} + \cdots} Y^{m_{0} + m_{1} + m_{2} + \cdots}\\
 &= \sum_{n=0}^{\infty} Y^{n} \sum_{\substack{
m_{0}, m_{1}, m_{2}, \dots \geq 0 \\
m_{0} + m_{1} + m_{2} + \cdots  = n}} X^{m_{1} + 2m_{2} + \cdots},
\end{align*}

\

so it is enough to show that

$$\sum_{\substack{
m_{0}, m_{1}, m_{2}, \dots \geq 0 \\
m_{0} + m_{1} + m_{2} + \cdots  = n}} X^{m_{1} + 2m_{2} + \cdots} = \sum_{j_{1}=0}^{\infty} \sum_{j_{2}=0}^{\infty} \cdots \sum_{j_{n}=0}^{\infty} X^{j_{1} + 2j_{2} + \cdots + nj_{n}}.$$

\

Note that we have a bijection

$$\{(m_{0}, m_{1}, m_{2}, \dots) \in \bZ_{\geq 0}^{\infty} : m_{0} + m_{1} + m_{2} + \cdots = n \} \lra \{(m_{1}, m_{2}, \dots) \in \bZ_{\geq 0}^{\infty} : m_{1} + m_{2} + \cdots \leq n \}.$$

\

given by $(m_{0}, m_{1}, m_{2}, \dots) \mapsto (m_{1}, m_{2}, \dots)$. This reduces our problem to the following:

$$\sum_{\substack{
m_{1}, m_{2}, \dots \geq 0 \\
m_{1} + m_{2} + \cdots  \leq n}} X^{m_{1} + 2m_{2} + \cdots} = \sum_{j_{1}=0}^{\infty} \sum_{j_{2}=0}^{\infty} \cdots \sum_{j_{n}=0}^{\infty} X^{j_{1} + 2j_{2} + \cdots + nj_{n}}.$$

\

If $n = 0$, both sides are $1$, so let $n \geq 1$. Let $A_{n}$ be the set of partitions whose number of parts counting with multiplicity (i.e., the row lengths of Young diagrams) is $\leq n$ and $B_{n}$ the set of partitions whose parts (i.e., the column lengths of Young diagrams) are $\leq n$. Then we have a bijection $A_{n} \lra B_{n}$ given by taking conjugation, so in particular, we have

$$\sum_{\ld \in A_{n}}X^{|\ld|} = \sum_{\ld \in B_{n}}X^{|\ld|}.$$

\

Now, noting that

$$\sum_{\ld \in A_{n}}X^{|\ld|} = \sum_{\substack{
m_{1}, m_{2}, \dots \geq 0 \\
m_{1} + m_{2} + \cdots  \leq n}} X^{m_{1} + 2m_{2} + \cdots}$$

\

and

$$\sum_{\ld \in B_{n}}X^{|\ld|} = \sum_{j_{1}=0}^{\infty} \sum_{j_{2}=0}^{\infty} \cdots \sum_{j_{n}=0}^{\infty} X^{j_{1} + 2j_{2} + \cdots + nj_{n}},$$

\

we finish the proof.
\end{proof}

\

\section{Proof of Proposition \ref{main}}\label{mainproof}

\hspace{3mm} For our proof of Proposition \ref{main}, we need to analyze polynomials that encode some information about random matrices in $\Mat_{n}(\bF_{q})$ and random permutations in $S_{n}$. Such a polynomial is called a ``cycle index''.

\

\subsection{Cycle index} Given any permutation group $G \leqslant S_{n}$, we define the \textbf{cycle index} of $G$ as the following polynomial:

$$\mc{Z}(G, \bs{x}) := \frac{1}{|G|}\sum_{g \in G}x_{1}^{m_{1}(g)} \cdots x_{n}^{m_{n}(g)},$$

\

where we recall that $m_{d}(g)$ means the number of $d$-cycles in the cycle decomposition of $g$. Cycle indices were used to solve various enumeration problems, and to our best knowledge, these were invented independently by Redfield \cite{Red} and P\'olya \cite{Pol}. (P\'olya's paper is in German, but there is an English translation \cite{PR} by Read). We will use the cycle index $\mc{Z}(S_{n}, \bs{x})$ of the full symmetric group $S_{n}$. Note that for any $\sigma, \tau \in S_{n}$, we note that $\sigma$ and $\tau$ are conjugate to each other in $S_{n}$ if and only if $x_{1}^{m_{1}(\sigma)} \cdots x_{n}^{m_{n}(\sigma)} = x_{1}^{m_{1}(\tau)} \cdots x_{n}^{m_{n}(\tau)}$, so $\mc{Z}(S_{n}, \bs{x})$ captures some information about the conjugate action $S_{n} \acts S_{n}$.

\

\hspace{3mm} We will also make use of the following polynomial:

$$\mc{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) := \frac{1}{|\GL_{n}(\bF_{q})|}\sum_{A \in \Mat_{n}(\bF_{q})}\prod_{P \in |\bA^{1}_{\bF_{q}}|}x_{P, \mu_{P}(A)},$$

\

where $\mu_{P}(A) = (\mu_{1}, \dots, \mu_{l})$ is the partition defined by the $P$-part of the $\bF_{q}[t]$-module structure defined by the matrix multiplication $A \acts \bF_{q}^{n}$:

$$A[P^{\infty}] \simeq \bF_{q}[t]/(P)^{\mu_{1}} \op \cdots \op \bF_{q}[t]/(P)^{\mu_{l}}.$$

\

We denoted by $|\bA^{1}_{\bF_{q}}|$ the set of all monic irreducible polynomials in $\bF_{q}[t]$. For $P \in |\bA^{1}_{\bF_{q}}|$ and a nonempty partition $\ld$, the notation $x_{P,\ld}$ means a formal variable associated to the pair $(P, \ld)$. For the empty partition, we define $x_{P,\es} := 1$. We call $\mc{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ the \textbf{cycle index} of the conjugate action $\GL_{n}(\bF_{q}) \acts \Mat_{n}(\bF_{q})$. The terminology makes sense because each monomial $\prod_{P \in |\bA^{1}_{\bF_{q}}|}x_{P, \mu_{P}(A)}$ is characterized by an orbit of such action. That is, two matrices $A$ and $B$ are in the same orbit if and only if $\prod_{P \in |\bA^{1}_{\bF_{q}}|}x_{P, \mu_{P}(A)} = \prod_{P \in |\bA^{1}_{\bF_{q}}|}x_{P, \mu_{P}(B)}$. The notion of $\mc{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ is introduced by Stong \cite{Sto}, generalizing a similar definition for the conjugate action $\GL_{n}(\bF_{q}) \acts \GL_{n}(\bF_{q})$ introduced by Kung \cite{Kun}. We will use the following factorization results for the generating functions for the cycle indicies.

\

\begin{lem}[p.14 of \cite{PR}]\label{fac1} In the ring $\bQ [\bs{x}] \llb u \rrb$ of formal power series in $u$ over the polynomial ring $\bQ [\bs{x}] = \bQ[x_{1},x_{2}, \dots]$, we have

$$\sum_{n=0}^{\infty}\mc{Z}(S_{n}, \bs{x}) u^{n} = \prod_{d=1}^{\infty}e^{x_{d}u^{d}/d}.$$
\end{lem}

\begin{proof} Given $\ld \vdash n$ (i.e., a partition of $n$), the number of permutations in $S_{n}$ with cycle type $\ld$ is

$$\frac{n!}{m_{1}(\ld)!1^{m_{1}(\ld)} \cdots m_{n}(\ld)!n^{m_{n}(\ld)}},$$

\

so

$$\mc{Z}(S_{n}, \bs{x}) = \sum_{\ld \vdash n}\frac{x_1^{m_{1}(\ld)} \cdots x_{n}^{m_{n}(\ld)}}{m_{1}(\ld)!1^{m_{1}(\ld)} \cdots m_{n}(\ld)!n^{m_{n}(\ld)}}.$$

\

Thus, we have

\begin{align*}
\sum_{n=0}^{\infty}\mc{Z}(S_{n}, \bs{x})u^{n} &= \sum_{n=0}^{\infty}\sum_{\ld \vdash n}\frac{x_1^{m_{1}(\ld)} \cdots x_{n}^{m_{n}(\ld)} u^{1 \cdot m_{1}(\ld)} \cdots u^{n \cdot m_{n}(\ld)}}{m_{1}(\ld)!1^{m_{1}(\ld)} \cdots m_{n}(\ld)!n^{m_{n}(\ld)}} \\
&= \sum_{\ld \in \mc{P}} \prod_{d=1}^{\infty} \frac{ (x_{d}u^{d}/d)^{m_{d}(\ld)} }{m_{d}(\ld)!} \\
&= \prod_{d=1}^{\infty}\sum_{m=0}^{\infty}\frac{(x_{d}u^{d}/d)^{m}}{m!} \\
&= \prod_{d=1}^{\infty}e^{x_{d}u^{d}/d},
\end{align*}

\

as desired.
\end{proof}

\

\begin{lem}[Lemma 1 (2) in \cite{Sto}]\label{fac2} In the ring $\bQ [\bs{x}] \llb u \rrb$ of formal power series in $u$ over the polynomial ring $\bQ [\bs{x}] = \bQ[x_{P,\ld}]_{P \in |\bA^{1}_{\bF_{q}}|, \ld \in \mc{P} \sm \{\es\}}$, we have

$$\sum_{n=0}^{\infty}\mc{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^n = \prod_{P \in |\bA^{1}_{\bF_{q}}|} \sum_{\ld \in \mc{P}} \frac{x_{P,\ld}u^{|\ld|\deg(P)}}{|\Aut_{\bF_{q}[t]}(H_{P, \ld})|},$$

\

where $H_{P, \ld} := H_{(P),\ld}^{\bF_{q}[t]}$, following the notation defined in the beginning of Section \ref{Stong}.
\end{lem}

\

\begin{rmk} Lemma \ref{fac2} says that on the left-hand side, the coefficient of $u^{n}$ is given by

$$\sum_{|\ld^{(1)}|\deg(P_{1}) + \cdots + |\ld^{(s)}|\deg(P_{s}) = n}\frac{x_{P_{1}, \ld^{(1)}} \cdots x_{P_{s}, \ld^{(s)}}}{|\Aut_{\bF_{q}[t]}(H_{P_{1}, \ld^{(1)}})| \cdots |\Aut_{\bF_{q}[t]}(H_{P_{s}, \ld^{(s)}})|},$$

\

where the sum is over distinct $P_{1}, \dots, P_{s} \in |\bA^{1}_{\bF_{q}}|$ and not necessarily distinct $\ld^{(1)}, \dots, \ld^{(s)}$ with the stipulated condition. This is a finite sum, so it makes sense to evaluate any complex numbers in the variables $x_{P, \ld}$ for $\ld \neq \es$ of the infinite product in the statement of Lemma \ref{fac2} and get an identity in $\bC \llb u \rrb$. (Recall that $x_{P, \es} = 1$.)
\end{rmk}

\

\subsection{Connection between two cycle indices} We define $\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ to be the polynomial given by specializing $x_{P,\ld} = x_{\deg(P)}^{|\ld|}$ in the cycle index $\mc{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$, which also makes sense for the empty partition because $x_{P,\es} = 1 = x_{\deg(P)}^{0}$. More explicitly, we have

\begin{align*}
\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) &= \frac{1}{|\GL_{n}(\bF_{q})|}\sum_{A \in \Mat_{n}(\bF_{q})}\prod_{P \in |\bA^{1}_{\bF_{q}}|}x_{\deg(P)}^{|\mu_{P}(A)|} \\
&= \frac{1}{|\GL_{n}(\bF_{q})|}\sum_{A \in \Mat_{n}(\bF_{q})}x_{1}^{m_{1}(A)} \cdots x_{n}^{m_{n}(A)},
\end{align*}

\

where $m_{d}(A)$ is the number of degree $d$ monic irreducible polynomials, counting with multiplicity, in $\bF_{q}[t]$ dividing the characteristic polynomial $f_{A}(t)$ of $A$. Note that only $x_{1}, \dots, x_{n}$ will occur in $\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q})$ because the degree of $f_A(t)$ is $n$ for $A \in \Mat_{n}(\bF_{q})$, so $f_A(t)$ is not divisible by any irreducible polynomial of degree $> n$. In the concluding remark of \cite{Sto}, Stong essentially observed that there is a close relationship between $\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ and $\mc{Z}(S_{n}, \bs{x})$ when $q$ is large. We rigorously formulate what he might have meant.

\

\begin{lem}[cf. ``Conclusion'' in \cite{Sto}]\label{conv} We have

$$\lim_{q \ra \infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) = \mc{Z}(S_{n}, \bs{x}),$$

\

meaning that the coefficients of $\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ converge to the coefficients of $\mc{Z}(S_{n}, \bs{x})$ as $q \ra \infty$.
\end{lem}

\begin{proof} Let $x_{d} \in \bC$ such that $|x_{d}| \leq 1$ for any $d \in \bZ_{\geq 1}$. Applying Lemma \ref{fac2} and then Lemma \ref{key}, we have

\begin{align*}
\sum_{n=0}^{\infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^{n} &= \prod_{P \in |\bA^{1}_{\bF_{q}}|} \sum_{\ld \in \mc{P}} \frac{(x_{\deg(P)}u^{\deg(P)})^{|\ld|}}{|\Aut_{\bF_{q}[t]}(H_{P, \ld})|} \\
&= \prod_{P \in |\bA^{1}_{\bF_{q}}|} \prod_{i=1}^{\infty}(1 - x_{\deg(P)}(q^{-i}u)^{\deg(P)})^{-1} \\
&= \prod_{d=1}^{\infty}\prod_{i=1}^{\infty}(1 - x_{d}(q^{-i}u)^{d})^{-M(q,d)},
\end{align*}

\

for $|u| < 1$, where $M(q, d)$ is the number of monic irreducible polynomials in $\bF_{q}[t]$ with degree $d$. Since $|x_{d}| \leq 1$ for all $d \geq 1$ so that we have

$$|\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})|^{1/n} \leq \left(\prod_{i=1}^{n}\frac{1}{1 - q^{-i}}\right)^{1/n} \leq \frac{1}{1 - q^{-1}},$$

\

the radius of convergence of the power series $\sum_{n=0}^{\infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^{n}$ in $u$ is at least $1 - q^{-1} > 0$. Since our statement only involves finitely many $x_{d}$, we may set all but finitely many of them to be $0$, say $x_{d} = 0$ for all $d > m$ for fixed $m \in \bZ_{\geq 1}$. This lets us have

$$\sum_{n=0}^{\infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^{n} = \prod_{d=1}^{m}\prod_{i=1}^{\infty}(1 - x_{d}(q^{-i}u)^{d})^{-M(q,d)}.$$

\

In what follows, we will use the fact that

$$\lim_{q \ra \infty}\frac{M(q,d)}{q^{d}/d} = 1,$$

\

for any $d \geq 1$, which can be found as Theorem 2.2 in \cite{Ros}. (Note that for $M(q, 1) = q$, so we do not need to take the limit to see this for $d = 1$.) Now, for $|u| < 1 - q^{-1}$, we have

\begin{align*}
\lim_{q \ra \infty}\sum_{n=0}^{\infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^{n} &= \prod_{d=1}^{m}\lim_{q \ra \infty}\prod_{i=1}^{\infty}(1 - x_{d}(q^{-i}u)^{d})^{-M(q,d)} \\
&= \prod_{d=1}^{m}\lim_{q \ra \infty}\prod_{i=1}^{\infty}\left((1 - x_{d}(q^{-i}u)^{d})^{-q^{d}/d}\right)^{M(q,d)d/q^{d}} \\
&= \prod_{d=1}^{m}\prod_{i=1}^{\infty}\lim_{q \ra \infty}\left((1 - x_{d}(q^{-i}u)^{d})^{-q^{d}/d}\right)^{M(q,d)d/q^{d}} \\
&= \prod_{d=1}^{\infty}\lim_{q \ra \infty}(1 - x_{d}(q^{-1}u)^{d})^{-q^{d}/d} \\
&= \prod_{d=1}^{\infty}\left(\lim_{q \ra \infty}\left(1 - \frac{x_{d}u^{d}}{q^{d}}\right)^{- q^{d}/(x_{d}u^{d})}\right)^{(x_{d}u^{d})/d} \\
&= \prod_{d=1}^{\infty}e^{x_{d}u^{d}/d}.
\end{align*}

\

Applying Lemma \ref{fac1} to the last expression, we see that

$$\lim_{q \ra \infty}f_{q}(u) = f(u),$$

\

where

$$f_{q}(u) := \sum_{n=0}^{\infty}\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})u^{n}$$

\

and

$$f(u) := \sum_{n=0}^{\infty}\mc{Z}(S_{n}, \bs{x})u^{n}$$

\

are holomorphic functions in $u$ defined for $|u| < 1-q^{-1}$. Take any $0 < \epsilon < 1-q^{-1}$ and write $C_{\eps} := \{u \in \bC : |u| = \eps\}$. By the Cauchy integral formula, we have

$$\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) = \frac{1}{2\pi i}\varointctrclockwise_{z \in C_{\eps}}\frac{f_{q}(z)}{z^{n+1}}dz$$

\

and

$$\mc{Z}(S_{n}, \bs{x}) = \frac{1}{2\pi i}\varointctrclockwise_{z \in C_{\eps}}\frac{f(z)}{z^{n+1}}dz.$$

\

Thus, by the Dominated Convergence Theorem, we have

\begin{align*}
\lim_{q \ra\infty} \bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) &= \lim_{q \ra \infty} \frac{1}{2\pi i}\varointctrclockwise_{z \in C_{\eps}}\frac{f_{q}(z)}{z^{n+1}}dz \\
&= \frac{1}{2\pi i}\varointctrclockwise_{z \in C_{\eps}}\frac{f(z)}{z^{n+1}}dz \\
&= \mc{Z}(S_{n}, \bs{x}).
\end{align*}

\

Since $x_{1}, \dots, x_{m}$ are arbitrary with the restriction $|x_{d}| \leq 1$, this is enough to prove the desired statement about the convergence of the coefficients in $x_{1}, \dots, x_{m}$ by the Cauchy integral formula for several variables (e.g., Theorem 2.1.3 of \cite{Fie}) with the Dominate Convergence Theorem.
\end{proof}

\

\subsection{Proof of Proposition \ref{main}} Again, Lemma \ref{conv} says that the coefficients of the polynomial

$$\bs{Z}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) = \frac{1}{|\GL_{n}(\bF_{q})|}\sum_{A \in \Mat_{n}(\bF_{q})}x_{1}^{m_{1}(A)} \cdots x_{n}^{m_{n}(A)} \in \bQ[x_{1}, \dots, x_{n}]$$

\

converge to the coefficients of the polynomial

$$\mc{Z}(S_{n}, \bs{x}) = \frac{1}{|S_{n}|}\sum_{\sigma \in S_{n}}x_{1}^{m_{1}(\sigma)} \cdots x_{n}^{m_{n}(\sigma)} \in \bQ[x_{1}, \dots, x_{n}]$$

\

as $q \ra \infty$, where the convergences are happening in $\bR$ or $\bC$, although the limits still lie in $\bQ$. This implies that, setting $x_{d,0} := 1$ and letting $x_{d,m}$ be formal for any $d, m \in \bZ_{\geq 1}$, the coefficients of

$$\hat{\bs{Z}}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) := \frac{1}{|\GL_{n}(\bF_{q})|}\sum_{A \in \Mat_{n}(\bF_{q})}x_{1, m_{1}(A)} \cdots x_{n, m_{n}(A)} \in \bQ
\left[\begin{array}{c}
x_{1,1}, \dots, x_{1,n}, \\
\cdots, \cdots, \cdots \\
x_{n,1}, \dots, x_{n,n}
\end{array}\right]$$

\

converge to the coefficients of

$$\hat{\mc{Z}}(S_{n}, \bs{x}) := \frac{1}{|S_{n}|}\sum_{\sigma \in S_{n}}x_{1, m_{1}(\sigma)} \cdots x_{n, m_{n}(\sigma)} \in \bQ
\left[\begin{array}{c}
x_{1,1}, \dots, x_{1,n}, \\
\cdots, \cdots, \cdots \\
x_{n,1}, \dots, x_{n,n}
\end{array}\right].$$

\

To refer back to this fact, we record this as follows.

\

\begin{lem}\label{conv'} As $q \ra \infty$, the coefficients of the polynomial $\hat{\bs{Z}}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x})$ converge to those of $\hat{\mc{Z}}(S_{n}, \bs{x})$. In particular, for any evaluation of the variables $x_{d,m}$ in $\bC$, we have

$$\lim_{q \ra \infty}\hat{\bs{Z}}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) = \hat{\mc{Z}}(S_{n}, \bs{x}).$$
\end{lem}

\

\hspace{3mm} We are now ready to prove Proposition \ref{main}:

\begin{proof}[Proof of Proposition \ref{main}] Recall the statements of Proposition \ref{main}, since we will use the notations in them. We will also use the notations given in this section. If we evaluate

\begin{itemize}
	\item $x_{d_{1},k_{1}} = \cdots = x_{d_{r},k_{r}} = 1$,
	\item $x_{d_{j},m} = 0$ for all $m \neq 0, k_{j}$ for $1 \leq j \leq r$, and
	\item $x_{d,m} = 1$ for all the other cases from above,
\end{itemize}

\

then we have

$$\hat{\bs{Z}}([\Mat_{n}/\GL_{n}](\bF_{q}), \bs{x}) = \frac{\left|\left\{\begin{array}{c}
A \in \Mat_{n}(\bF_{q}) : f_{A}(t) \text{ has } k_{j} \text{ irreducible}  \\ \text{factors of degree } d_{j} \text{ for } 1 \leq j \leq r
\end{array}
\right\}\right|}{|\GL_{n}(\bF_{q})|},$$

\

where we count the factors with multiplicity, and

$$\hat{\mc{Z}}(S_{n}, \bs{x}) = \frac{\left|\left\{\begin{array}{c}
\sigma \in S_{n} : \sigma \text{ has } k_{j} \text{ cycles}  \\ \text{of length } d_{j} \text{ for } 1 \leq j \leq r
\end{array}
\right\}\right|}{|S_{n}|}.$$

\

Thus, applying Lemma \ref{conv'} and noting that

$$\lim_{q \ra \infty}\frac{|\Mat_{n}(\bF_{q})|}{|\GL_{n}(\bF_{q})|} = \lim_{q \ra \infty}(1 - q^{-1})^{-1} \cdots (1 - q^{-n})^{-1} = 1,$$

\

we have

$$\lim_{q \ra \infty}\Prob_{A \in \Mat_{n}(\bF_{q})}\left(\begin{array}{c}
A \in \Mat_{n}(\bF_{q}) : f_{A}(t) \text{ has } k_{j}  \\ \text{irreducible factors of} \\
\text{degree } d_{j} \text{ for } 1 \leq j \leq r
\end{array}
\right) = \Prob_{\sigma \in S_{n}}\left(\begin{array}{c}
\sigma \in S_{n} : \sigma \text{ has } k_{j} \text{ cycles}  \\ \text{of length } d_{j} \text{ for } 1 \leq j \leq r
\end{array}
\right),$$

\

as desired.
\end{proof}

\

\section{Combinatorial proof of Lemma \ref{perm}}\label{permproof}

\hspace{3mm} In this section, we give a proof of Lemma \ref{perm}, a result due to Shepp and Lloyd in \cite{SL}, originally proven by computing the characteristic functions of the distributions. Unlike their proof, we will directly compute the desired probability with a more combinatorial method. Most of our argument will be encoded in the following lemma.

\

\begin{lem}\label{issue} Take

\begin{itemize}
	\item $x_{d_{j},m} = 0$ with $1 \leq j \leq s$ and $m \neq 0$,
	\item $x_{d_{j},m} = 0$ with $s + 1 \leq j \leq r$ and $m \neq 0, k_{j}$, and
	\item $x_{d,m} = 1$ for any other ones not on the above list.
\end{itemize}

Then $$\sum_{n=0}^{\infty} \hat{\mc{Z}}(S_{n}, \bs{x})u^{n} = \left( \frac{e^{-u^{d_1}/d_{1}} \cdots e^{-u^{d_r}/d_{r}}}{1 - u} \right) \cdot \left( \prod_{j=s+1}^{r} \left(1 + \frac{(u^{d_j}/d_{j})^{k_{j}}}{k_{j}!} \right) \right) \in \bC \llb u \rrb.$$
\end{lem}

\

\begin{rmk}\label{Tau} Before proving Lemma \ref{issue}, we first see how it implies Lemma \ref{perm}.  We will make use of the following useful observation: for any $f(u) = c_{0} + c_{1}u + c_{2}u^{2} + \cdots \in \bC \llb u \rrb$, if $f(1) = c_{0} + c_{1} + c_{2} + \cdots$ exists, then the limit of the coefficient sequence of the power series

$$\frac{f(u)}{1 - u} = b_{0} + b_{1}u + b_{2}u^{2} + \cdots$$

\

is given by

$$\lim_{n \ra \infty} b_{n} = f(1),$$

\

because $(1 - u)^{-1} = 1 + u + u^{2} + \cdots$ so that $b_{n} = c_{0} + c_{1} + \cdots + c_{n}$, which is the coefficient of $u^{n}$ for the power series $f(u)(1 - u)^{-1}$.
\end{rmk}

\

\begin{proof}[Proof of Lemma \ref{perm}] Using Lemma \ref{issue} with $s = r$, we have

$$\sum_{n=0}^{\infty}\Prob_{\sigma \in S_{n}}\left( m_{d_{j}}(\sigma) = 0 \text{ for } 1 \leq j \leq r \right) u^{n} = \frac{e^{-u^{d_1}/d_{1}} \cdots e^{-u^{d_r}/d_{r}}}{1 - u}.$$

\

Using Remark \ref{Tau}, this implies that

$$\lim_{n \ra \infty}\Prob_{\sigma \in S_{n}}\left(m_{d_{j}}(\sigma) = 0 \text{ for } 1 \leq j \leq r\right) = e^{-1/d_{1}} \cdots e^{-1/d_{r}},$$

\

as claimed. For the general case, we work with the induction on the number of nonzero integers among $k_{1}, \dots, k_{r}$. The base case is when such number is $0$, which is exactly what we have proved above. Suppose that the result is true when such number is $0, 1, \dots, s-1$, where $s - 1 \geq 0$. Then we assume that there are $s$ nonzero elements among $k_{1}, \dots, k_{r}$, so permuting them if necessary, say

$$k_{1} = \cdots = k_{s} = 0,$$

\

while

$$k_{s+1}, \dots, k_{r} \neq 0.$$

\

Applying Lemma \ref{issue}, we get

\begin{align*}
\sum_{n=0}^{\infty}&\Prob_{\sigma \in S_{n}}\left(
\begin{array}{c}
m_{d_{j}}(\sigma) = 0 \text{ for } 1 \leq j \leq s, \text{ and } \\
m_{d_{s+1}}(\sigma) \in \{0, k_{s+1}\}, \dots, m_{d_{r}}(\sigma) \in \{0, k_{r}\}
\end{array}
\right) u^{n} \\
&= \left( \frac{e^{-u^{d_1}/d_{1}} \cdots e^{-u^{d_r}/d_{r}}}{1 - u} \right) \cdot \left( \prod_{j=s+1}^{r} \left(1 + \frac{(u^{d_j}/d_{j})^{k_{j}}}{k_{j}!} \right) \right),
\end{align*}

\

so

\begin{align*}
\sum_{\eps_{s+1} \in \{0, k_{s+1}\}} \cdots \sum_{\eps_{r} \in \{0, k_{r}\}} & \lim_{n \ra \infty}\Prob_{\sigma \in S_{n}}\left(
\begin{array}{c}
m_{d_{1}}(\sigma) = \cdots = m_{d_{s}}(\sigma) = 0, \\
m_{d_{s+1}}(\sigma) = \eps_{s+1}, \dots, m_{d_{r}}(\sigma) = \eps_{r}
\end{array}
\right) \\
&= e^{-1/d_{1}} \cdots e^{-1/d_{r}} \prod_{j=s+1}^{r} \left(1 + \frac{(1/d_{j})^{k_{j}}}{k_{j}!} \right)
\end{align*}

\

Applying the induction hypothesis after expanding the right-hand side of the above identity implies that

$$\Prob_{\sigma \in S_{n}}\left(
\begin{array}{c}
m_{d_{1}}(\sigma) = \cdots = m_{d_{s}}(\sigma) = 0, \\
m_{d_{s+1}}(\sigma) = k_{s+1}, \dots, m_{d_{r}}(\sigma) = k_{r}
\end{array} \right) = e^{-1/d_{1}} \cdots e^{-1/d_{r}} \frac{(1/d_{s+1})^{k_{s+1}}}{k_{s+1}!} \cdots \frac{(1/d_{r})^{k_{r}}}{k_{r}!},$$

\

which finishes the proof.
\end{proof}

\

\hspace{3mm} Finally, we provide a proof of Lemma \ref{issue}.

\

\begin{proof}[Proof of Lemma \ref{issue}] For now, let us not evaluate any variables in $\bs{x} = (x_{d,m})$. Recall that
$$\hat{\mc{Z}}(S_{n}, \bs{x}) = \frac{1}{|S_{n}|}\sum_{\sigma \in S_{n}}x_{1, m_{1}(\sigma)} \cdots x_{n, m_{n}(\sigma)},$$

\

where $x_{d,0} = 1$ and $x_{d,m}$ are defined to be formal variables for $d, m \geq 1$. We observe that

$$\hat{\mc{Z}}(S_{n}, \bs{x}) = \sum_{\ld \vdash n}\frac{x_{1,m_{1}(\ld)} \cdots x_{n,m_{n}(\ld)}}{m_{1}(\ld)!1^{m_{1}(\ld)} \cdots m_{n}(\ld)!n^{m_{n}(\ld)}},$$

\

so

$$\sum_{n=0}^{\infty}\hat{\mc{Z}}(S_{n}, \bs{x})u^{n} = \prod_{d=1}^{\infty}\sum_{m=0}^{\infty}\frac{x_{d,m}u^{dm}}{m!d^{m}} = \prod_{d=1}^{\infty}\sum_{m=0}^{\infty}\frac{x_{d,m}(u^d/d)^{m}}{m!}.$$

\

Note that taking all $x_{d,m} = 1$ in the identity, we get

$$\frac{1}{1 - u} = 1 + u + u^{2} + \cdots + = \prod_{d=1}^{\infty}\sum_{m=0}^{\infty}\frac{(u^d/d)^{m}}{m!}.$$

\

Hence, with the given evaluation in the variables $\bs{x} = (x_{d,m})$, we get

\begin{align*}
\sum_{n=0}^{\infty}\hat{\mc{Z}}(S_{n}, \bs{x}) u^{n} &= \left( \prod_{d=1}^{\infty}\sum_{m=0}^{\infty}\frac{(u^d/d)^{m}}{m!} \right) \cdot \left( \prod_{j=1}^{r} \left( \sum_{m=0}^{\infty}\frac{(u^{d_j}/d_{j})^{m}}{m!} \right)^{-1} \right) \cdot \left( \prod_{j=s+1}^{r} \left(1 + \frac{(u^{d_j}/d_{j})^{k_{j}}}{k_{j}!} \right) \right) \\
&= \left( \frac{e^{-u^{d_1}/d_{1}} \cdots e^{-u^{d_r}/d_{r}}}{1 - u} \right) \cdot \left( \prod_{j=s+1}^{r} \left(1 + \frac{(u^{d_j}/d_{j})^{k_{j}}}{k_{j}!} \right) \right).
\end{align*}

\

Since the identities makes sense as those of holomorphic functions in $u \in \bC$ with $|u| < 1$, this finishes the proof.
\end{proof}

\

\newpage

Combinatorics (math.CO); Number Theory (math.NT)