  A set of nonnegative matrices $\mathcal{M}=\{M_1, M_2, \ldots, M_k\}$ is
called primitive if there exist indices $i_1, i_2, \ldots, i_m$ such that
$M_{i_1} M_{i_2} \ldots M_{i_m}$ is positive (i.e. has all its entries $>0$).
The length of the shortest such product is called the exponent of
$\mathcal{M}$. The concept of primitive sets of matrices comes up in a number
of problems within control theory, non-homogeneous Markov chains, automata
theory etc. Recently, connections between synchronizing automata and primitive
sets of matrices were established. In the present paper, we significantly
strengthen these links by providing equivalence results, both in terms of
combinatorial characterization, and computational aspects. We study the maximal
exponent among all primitive sets of $n \times n$ matrices, which we denote by
$\exp(n)$. We prove that $\lim_{n\rightarrow\infty} \tfrac{\log \exp(n)}{n} =
\tfrac{\log 3}{3}$, and moreover, we establish that this bound leads to a
resolution of the \v{C}ern\'{y} problem for carefully synchronizing automata.
We also study the set of matrices with no zero rows and columns, denoted by
$\mathcal{NZ}$, due to its intriguing connections to the \v{C}ern\'{y}
conjecture and the recent generalization of Perron-Frobenius theory for this
class. We characterize computational complexity of different problems related
to the exponent of $\mathcal{NZ}$ matrix sets, and present a quadratic bound on
the exponents of sets belonging to a special subclass. Namely, we show that the
exponent of a set of matrices having total support is bounded by $2n^2 -5n +5$.
