In this talk, we are concerned with long-time behavior of Euler-Maruyama schemes associated with regime-switching diffusion processes. The existence and uniqueness of numerical invariant measures are addressed (i) for regime-switching diffusion processes with finite state spaces by the Perron-Frobenius theorem if the ``averaging condition" holds, and, for the case of reversible Markov chain, via the principal eigenvalue approach provided that the principal eigenvalue is positive; (ii) for regime-switching diffusion processes with countable state spaces by means of a finite partition method and an M-Matrix theory. We also reveal that numerical invariant measures converge in the Wasserstein metric to the underlying ones.