We derive a reduction formula that expresses the survival rate for the three-flavor neutrino oscillation using the two-flavor one to next-to-leading order when there is one resonance due to the matter effect. We numerically find that the next-to-leading reduction formula is extremely accurate and the improvement is relevant for the precision test of solar neutrino oscillation and the indirect measurment of CP violation in the leptonic sector. We also derive a reduction formula, which is slightly different from that previously obtained, in the case when there are two resonances. We numerically verify that this reduction formula is quite accurate and is valid for a wider parameter region than are those previously obtained.