The vehicle routing problem (VRP) is a challenging combinatorial optimization problem. This research focuses on the problem under which a manufacturer needs to outsource materials from other suppliers and to ship the materials back to the company. Heterogeneous vehicles are available to ship the materials, and each vehicle has a limited loading capacity and a limited travelling distance. The purpose of this research is to study a multiple vehicle routing problem (MVRP) with soft time window and heterogeneous vehicles. Two models, using mixed integer programming (MIP) and genetic algorithm (GA), are developed to solve the problem. The MIP model is first constructed to minimize the total transportation cost, which includes the assignment cost, travelling cost, and the tardiness cost, for the manufacturer. The optimal solution can present multiple vehicle routing and the loading size of each vehicle in each period. The GA is next applied to solve the problem so that a near-optimal solution can be obtained when the problem is too difficult to be solved using the MIP. A case of a food manufacturing company is used to examine the practicality of the proposed MIP model and the GA model. The results show that the MIP model can obtain the optimal solution under a short computational time when the scale of the problem is small. When the problem becomes non-deterministic polynomial hard (NP-hard), the MIP model cannot find the optimal solution. On the other hand, the GA model can obtain a near-optimal solution within a reasonable amount of computational time. This paper is related to several important topics of the Symmetry journal in the areas of mathematics and computer science theory and methods. In the area of mathematics, the theories of linear and non-linear algebraic structures and information technology are adopted. In the area of computer science, theory and methods, and metaheuristics are applied.Unit of measure: min.Symmetry 2018, 10, 650 9 of 20 Table 6 lists the problem configurations for the three cases, including the number of periods, number of suppliers, number of vehicles, size of the transportation network, number of variables for the MIP, and number of constraints for the MIP. In the first case, there are five periods, the number of suppliers is five, the number of vehicles is three, the transportation network is 6 × 6, and the numbers of variables and constraints for the MIP model are 931 and 1171, respectively. In the second case, there are seven periods, the number of suppliers is nine, the number of vehicles is four, the transportation network is 10 × 10, and the numbers of variables and constraints for the MIP model are 3903 and 3977, respectively. In the third case, there are nine periods, the number of suppliers is 12, the number of vehicles is five, the transportation network is 13 × 13, and the numbers of variables and constraints for the MIP model are 9758 and 9919, respectively.