  Let $(R, \mathfrak{m})$ be a complete discrete valuation ring with the finite
residue field $R/\mathfrak{m} = \mathbb{F}_{q}$. Given a monic polynomial $P(t)
\in R[t]$ whose reduction modulo $\mathfrak{m}$ gives an irreducible polynomial
$\bar{P}(t) \in \mathbb{F}_{q}[t]$, we initiate the investigation of the
distribution of $\mathrm{coker}(P(A))$, where $A \in \mathrm{Mat}_{n}(R)$ is
randomly chosen with respect to the Haar probability measure on the additive
group $\mathrm{Mat}_{n}(R)$ of $n \times n$ $R$-matrices. One of our main
results generalizes two results of Friedman and Washington. Our other results
are related to the distribution of the $\bar{P}$-part of a random matrix
$\bar{A} \in \mathrm{Mat}_{n}(\mathbb{F}_{q})$ with respect to the uniform
distribution, and one of them generalizes a result of Fulman. We heuristically
relate our results to a celebrated conjecture of Cohen and Lenstra, which
predicts that given an odd prime $p$, any finite abelian $p$-group (i.e.,
$\mathbb{Z}_{p}$-module) $H$ occurs as the $p$-part of the class group of a
random imaginary quadratic field extension of $\mathbb{Q}$ with a probability
inversely proportional to $|\mathrm{Aut}_{\mathbb{Z}}(H)|$. We review three
different heuristics for the conjecture of Cohen and Lenstra, and they are all
related to special cases of our main conjecture, which we prove as our main
theorems. For proofs, we use some concrete combinatorial connections between
$\mathrm{Mat}_{n}(R)$ and $\mathrm{Mat}_{n}(\mathbb{F}_{q})$ to translate our
problems about a Haar-random matrix in $\mathrm{Mat}_{n}(R)$ into problems
about a random matrix in $\mathrm{Mat}_{n}(\mathbb{F}_{q})$ with respect to the
uniform distribution.
