QCD equation of state in a virial expansion