Monte Carlo algorithm for hot phonons in polar semiconductors