Computing the shortest elementary flux modes in genome-scale metabolic networks