Using Numerical Stochastic Perturbation Theory within three-dimensional pure SU(3) gauge theory, we estimate the last unknown renormalization constant that is needed for converting the vacuum energy density of this model from lattice regularization to the MS scheme. Making use of a previous non-perturbative lattice measurement of the plaquette expectation value in three dimensions, this allows us to approximate the first non-perturbative coefficient that appears in the weak-coupling expansion of hot QCD pressure.