To solve this integrated problem, we develop a method that combines Benders’ decomposition and column generation. We conduct computational experiments with real-world data from a major U.S. carrier.