Eulerian quasisymmetric functions