We calculate the three-loop thermodynamic potential of QCD at finite temperature and chemical potential(s) using the hard-thermal-loop perturbation theory (HTLpt) reorganization of finite temperature and density QCD. The resulting analytic thermodynamic potential allows us to compute the pressure, energy density, and entropy density of the quark-gluon plasma. Using these we calculate the trace anomaly, speed of sound, and second-, fourth-, and sixth-order quark number susceptibilities. For all observables considered we find good agreement between our three-loop HTLpt calculations and available lattice data for temperatures above approximately 300 MeV.