The properties of hot hadronic matter are of great importance to the studies of heavy-ion collisions, cosmology and compact star formation. I briefly outline the current methods in use in the lattice simulations of QCD thermodynamics at zero and nonzero density. I discuss the most recent results for the QCD phase transition, critical behavior and the equation of state.