Shock waves play an important role in turbulent astrophysical media by compressing the gas and dissipating the turbulent energy into the thermal energy. Here, we study shocks in magnetohydrodynamic turbulence using high-resolution simulations. Turbulent Mach numbers of M-turb = 0.5-7 and initial magnetic fields of plasma beta beta(0) = 0.1-10 are considered, targeting turbulences in interstellar and intracluster media. Specifically, we present the statistics of fast and slow shocks, such as the distribution of shock Mach numbers (M-s) and the energy dissipation at shocks, based on refined methodologies for their quantifications. While most shocks form with low M-s, strong shocks follow exponentially decreasing distributions of M-s. More shocks appear for larger M-turb and larger beta(0). Fast shock populations dominate over slow shocks if beta(0) >> 1, but substantial populations of slow shocks develop in the cases of beta less than or similar to 1, i.e., strong background fields. The shock dissipation of turbulent energy occurs preferentially at fast shocks with M-s less than or similar to of a few to several, and the dissipation at strong shocks shows exponentially decreasing functions of M-s. The energy dissipation at shocks, normalized to the energy injection, epsilon(shock)/epsilon(inj), is estimated to be in the range of similar to 0.1-0.5, except for the case of M-turb = 0.5 and beta(0) = 0.1, where the shock dissipation is negligible. The fraction decreases with M-turb; it is close to similar to 0.4-0.6 for M-turb = 0.5, while it is similar to 0.1-0.25 for M-turb = 7. The rest of the turbulent energy is expected to dissipate through the turbulent cascade. Our work will add insights into the interpretations of physical processes in turbulent interstellar and intracluster media.