Sequential computation of elementary modes and minimal cut sets in genome-scale metabolic networks using alternate integer linear programming