This answer does a good job at explaining that if $I$ is primary monomial ideal in $k[x_1, \dots, x_n]$, then $I = (x_{i_1}^{a_1}, \ldots, x_{i_m}^{a_m}, m_1, \ldots, m_k)$ where $m_1, \ldots, m_k$ only involve variables among $x_{i_1}, \ldots, x_{i_m}$.
However, it doesn't show that if $I = (x_{i_1}^{a_1}, \ldots, x_{i_m}^{a_m}, m_1, \ldots, m_k)$ as above, then $I$ is primary.
How do we prove this?
One thing we know is that if $I = (x_{i_1}^{a_1}, \ldots, x_{i_m}^{a_m}, m_1, \ldots, m_k)$ where $m_1, \ldots, m_k$ only involve variables among $x_{i_1}, \ldots, x_{i_m}$, then ${\rm rad}(I)=(x_{i_1}, \ldots, x_{i_m})$ which is prime.