We investigate the thermalization of high-energy particles injected from the perturbative decay of inflaton during the pre-thermal phase of reheating in detail. In general, thermalization takes a relatively long time in a low-temperature plasma; therefore, the instantaneous thermalization approximation is not justified, even for the reheating of the Standard Model (SM) sector. We consider a pure Yang-Mills (YM) theory as an approximation of the SM sector or a possible dark sector, considering the Landau-Pomeranchuk-Migdal effect, a quantum interference effect in a finite temperature plasma. We perform the first numerical calculation to solve the time evolution of the system, including the redshift due to the expansion of the Universe, and show the details of the temperature evolution near the maximum and the behavior of the quasi-attractors at later times. The maximal temperature Tmax and time scale tmax are determined quantitatively, such as Tmax ≃ 0.05 × ΓIMPI2/mI32/5mI\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\left({\\Gamma}_I{M}_{\ extrm{PI}}^2/{m}_I^3\\right)}^{2/5}{m}_I $$\\end{document} and tmax ≃ 2 × 103 × ΓIMPI2/mI3−3/5mI−1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\left({\\Gamma}_I{M}_{\ extrm{PI}}^2/{m}_I^3\\right)}^{-3/5}{m}_I^{-1} $$\\end{document} in the SM-like system, where mI and ΓI are the mass and decay rate of inflaton. We also provide a similar formula for pure SU(N) and SO(N) YM theories for general values of N and coupling constant α, including Tmax ∝ α4/5 and tmax ∝ N−2α−16/5 behaviors and their numerical coefficients. The thermalization occurs in a finite time scale, resulting in a lower maximal temperature of the Universe after inflation than that under the instantaneous thermalization approximation.