Let
$$\chi(\lambda) = \det(\lambda I_n - A ) = \lambda^n - e_1 \lambda^{n-1} + e_2 \lambda^{n-2} + \cdots + (-1)^n e_n$$
be the characteristic polynomial of $A$. By Cayley–Hamilton theorem (which is valid for any commutative ring), we have
$$\chi(A) = A^n - e_1 A^{n-1} + e_2 A^{n-2} + \cdots + (-1)^n e_n I_n = 0\tag{*1}$$
Let $p_k = \text{tr}A^k$ for $k = 1, \ldots, n$. The $e_k$ are elementary symmetric polynomials in the roots of $\chi(\lambda)$ and they are related to $p_k$ through the Newton identities (again valid for any commutative ring):
$$\begin{align}
e_1 = & p_1,\\
2 e_2 = & e_1 p_1 - p_2,\\
3 e_3 = & e_2 p_1 - e_1 p_2 + p_3,\\
4 e_4 = & e_3 p_1 - e_2 p_2 + e_1 p_3 - p_4,\\
\vdots
\end{align}$$
As a result,
$$p_1 = p_2 = \cdots = p_n = 0\quad\implies\quad e_1 = 2 e_2 = 3 e_3 = \cdots = n e_n = 0$$
Multiply $(*1)$ by $n!$ then give us $\;n!A^n = 0\;$.