We show the formulation and the result of test calculation of the overlap of excites states in the quasiparticle random-phase approximation for calculating the nuclear matrix elements of the neutrino-less double beta decay. Our method uses the ground states of the QRPA explicitly. The feasibility of that calculation is demonstrated, and the effectiveness of the truncation used in our method is shown.