The direct simulation Monte Carlo method for the Boltzmann equation is modified by an additional displacement in the advection process and an enhanced collision rate in order to obtain the exact hard sphere equation of state at all densities. This leads to consistent thermodynamic and transport properties in the low density (Boltzmann) regime. At higher densities transport properties are comparable to the predictions of the Enskog model. The algorithm is faster than molecular dynamics at low and moderate densities and readily run on a parallel architecture