Approximate constants of motion for vibrational spectra of many-oscillator systems with multiple anharmonic resonances