Path integral Monte Carlo method in the expanded ensemble is used for calculation of the ratio of partition functions for different classes of permutations treating the problem of several interacting identical particles (fermions) in an external field. Wang-Landau algorithm is used for adjustment of balancing factors. For systems consisting of greater than two number of particles we propose an advanced variant of our approach which implies calculation of the ratio of positive and negative contributions to the partition function. Densities and energies of the sequence of excited states starting from the ground state for a system of non interacting quantum particles are calculated in turn, one by one, by means of considering systems with artificially excluded lowest energy levels and further obtaining of the ”ground state” of each next system constructed in this way. The idea of evaluation of densities of excited states for a quantum system with interparticle interaction by evaluating density difference betwee