We extend to QCD an efficient method for lattice gauge theory with dynamica
l fermions. Once the eigenvalues of the Dirac operator and the density of s
tates of pure gluonic configurations at a set of plaquette energies (propor
tional to the gauge action) are computed, thermodynamical quantities derivi
ng from the partition function can be obtained for arbitrary flavor number,
quark masses and wide range of coupling constants, without additional comp
utational cost. Results for the chiral condensate and gauge action are pres
ented on the 10(4) lattice at flavor number N-f = 0, 1, 2, 3, 4 and many qu
ark masses and coupling constants. New results in the chiral limit for the
gauge action and its correlation with the chiral condensate, which are usef
ul for analyzing the QCD chiral phase structure, are also provided.