Color path-integral Monte Carlo simulations of quark–gluon plasma