Integrated Tomato Picking and Distribution Scheduling Based on Maturity

: Agriproducts have the characteristics of short lifespan and quality decay due to the maturity factor. With the development of e-commerce, high timelines and quality have become a new pursuit for agriproduct online retailing. To satisfy the new demands of customers, reducing the time from receiving orders to distribution and improving agriproduct quality are signiﬁcantly needed advancements. In this study, we focus on the joint optimization of the fulﬁllment of online tomato orders that integrates picking and distribution simultaneously within the context of the farm-to-door model. A tomato maturity model with a ﬁrmness indicator is proposed ﬁrstly. Then, we incorporate the tomato maturity model function into the integrated picking and distribution schedule and formulate a multiple-vehicle routing problem with time windows. Next, to solve the model, an improved genetic algorithm (the sweep-adaptive genetic algorithm, S-AGA) is addressed. Finally, we prove the validity of the proposed model and the superiority of S-AGA with di ﬀ erent numerical experiments. The results show that signiﬁcant improvements are obtained in the overall tomato supply chain e ﬃ ciency and quality. For instance, tomato quality and customer satisfaction increased by 5% when considering the joint optimization, and the order processing speed increased over 90% compared with traditional GA. This study could provide scientiﬁc tomato picking and distribution scheduling to satisfy the multiple requirements of consumers and improve agricultural and logistics sustainability.


Introduction
With the development of e-commerce platforms, the farm-to-door model has received critical attention in recent years. It can not only avoid the redundant part of distribution effectively but can also reduce the loss in the handling of agriproducts during the process. Merchants gain more profit by providing a value-added service for agriproducts through picking, packing, warehousing, and transportation [1], and consumers get their orders faster in the meanwhile. However, the quality of agriproducts will decline after picking due to the property of maturation. Time will increase quality decay after picking and distribution, affecting consumer satisfaction. Fresh agriproducts, represented by tomatoes, have a short lifespan and quality decay nature. Respiration affects the food flavor, and the agriproducts even lose commercial value during storage and transportation [2], resulting in huge losses [3]. Therefore, scheduling a proper tomato picking and distribution scheme so as to reduce quality decay, improve the sustainability of the tomato supply chain, and satisfy consumers' multiple requirements has been an essential problem for tomato trading in the farm-to-door context.
In Section 4, we present a sweep-adaptive genetic algorithm to solve the proposed model. In Section 5, we report the numerical results and the comparison of computer performance between the genetic algorithm (GA) and the sweep-adaptive genetic algorithm (S-AGA). In Section 6, we discuss the conclusions and future works.

Tomato Maturity Model
Maturity is an important indicator to measure tomato quality [15], so, in this study, we propose a tomato maturity model to be the measurement index of quality. Nowadays, the measurement of tomato maturity and quality mainly relies on artificial sensory assessment [16], physicochemical analysis [17], and nondestructive testing technology [18]. Tomato maturity is a complex physiological and biochemical process, during which a series of changes will take place. The main characteristic is the deterioration of firmness. Therefore, firmness is an important indicator of fruit maturity [19]. According to the study and measurement of Holt [20], Patrick, and Marie [21] and Zsom et al. [22] on the nondestructive testing of tomato firmness, we can take firmness as a critical indicator to measure tomato maturity. Although some studies have added maturity to the research on logistics distribution [5,23,24], most of them only considered the quality change of products in the process of distribution and ignored the change in the process of large-scale orders and long-time picking. Moreover, the rate of ripeness also influences the picking process. The picking process and distribution process influence and restrict each other. Therefore, this study considers integrated tomato picking and distribution scheduling based on maturity.
The tomato maturity can be divided into six stages according to six ripening periods [25], and tomatoes mature naturally after being picked during the six stages. Though the duration of tomato transformation from one stage to the next is different, its surface characteristics and physicochemical properties have a certain trend [15]. Therefore, tomato maturity can be expressed by a certain character of tomatoes.
Firmness is an important indicator to measure tomato maturity, and it shows a very obvious change with time [19]. In this study, partial least-square is used to establish the maturity model. The Dalian Traditional No.1 tomato is selected as the experimental material and measured at different maturity stages. We stored the samples in a constant temperature and humidity box at 25 • C and employed the professional US FTC TMS-Pro Texture analyzer with a cylindrical probe to get the data of tomato firmness. Six representative data points were selected in each stage. We have chosen the maximum data points (Data 1), the middle data points (Data 2), and the minimum data points (Data 3) of each stage, respectively, to represent the firmness of tomato maturity and the changes of firmness on different days were obtained, as shown with the dotted line in Figure 1. It can be seen that firmness is negatively correlated with time. The SPSS statistics 19.0 package is applied for data analysis.
Shiue [26] pointed out that the deterioration rate of perishable items with a lifespan is usually expressed in exponential form. Therefore, this study uses the exponential form to carry out the maturity model. The tomato maturity model is shown in Equation (1). Here, t denotes the time since picking; R(t) denotes tomato firmness at different maturity stages. Visualize the fitted function, as shown in the solid curve in Figure 1, being close to the variation trend of the data curve. Numerical examination is shown in Table 1. In the model, t is selected to produce an adjusted goodness of fit (R-square), R 2 = 0.977. The significance test (p = 0.000) and the F-test are both passed (F = 171.914 < 224.583). It can be concluded that there is an extremely significant correlation between time and firmness. Each specific range of firmness corresponds to a different tomato maturity stage, guiding the practical tomato picking process. In this study, we use the principle of proximity to classify tomato maturity based on firmness. For example, Maturity 8 corresponds to the tomato firmness interval (31,35].

Model Formulation
We considered an integrated tomato picking and distribution schedule based on maturity and modeled this problem as a multiple-vehicle routing problem with time windows (m-VRPTW). The objective is to minimize the total cost, including distribution cost and discipline cost.

Problem Description
Farm-to-door distribution is a superior model to reduce the logistics process and costs because only suppliers and customers participate in the whole agriproduct supply chain. However, in such a context, the supply chain will meet a lot of challenges on account of repaid delivery and high customer requirements. The farm-to-door distribution process of agriproduct online retailing contains order processing, picking, packaging, vehicle assignment, and distribution. In this study, we illustrate the picking, vehicle assignment, and distribution to formulate an integrated tomato picking and distribution schedule based on maturity, as shown in Figure 2. It contains one tomato production base, one packaging area, various consumers located in different places, and multiple types of vehicles. The tomato production base receives numerous orders = { , , … , }, including online and offline orders within the specified scope. After order processing, tomatoes are picked, packed, and delivered directly to customers = {1,2, … , } with the time required windows (ET i ,LT i ) , using vehicle = { , , … , | | } of = { , , … , | | } . Optimal vehicle routings are programmed after considering tomato maturity R(t) and time windows (ET i ,LT i ). The tomato maturity function, order details including coordinates, demand, and the type and load of vehicles in the production base are known.

Model Formulation
We considered an integrated tomato picking and distribution schedule based on maturity and modeled this problem as a multiple-vehicle routing problem with time windows (m-VRPTW). The objective is to minimize the total cost, including distribution cost and discipline cost.

Problem Description
Farm-to-door distribution is a superior model to reduce the logistics process and costs because only suppliers and customers participate in the whole agriproduct supply chain. However, in such a context, the supply chain will meet a lot of challenges on account of repaid delivery and high customer requirements. The farm-to-door distribution process of agriproduct online retailing contains order processing, picking, packaging, vehicle assignment, and distribution. In this study, we illustrate the picking, vehicle assignment, and distribution to formulate an integrated tomato picking and distribution schedule based on maturity, as shown in Figure 2. It contains one tomato production base, one packaging area, various consumers located in different places, and multiple types of vehicles. The tomato production base receives numerous orders O = {O 1 , O 2 , . . . , O i }, including online and offline orders within the specified scope. After order processing, tomatoes are picked, packed, and delivered directly to customers N = {1, 2, . . . , d} with the time required windows (ET i , LT i ), using vehicle K = k 1 , k 2 , . . . , k |k| of M = m 1 , m 2 , . . . , m |m| . Optimal vehicle routings are programmed after considering tomato maturity R(t) and time windows (ET i , LT i ). The tomato maturity function, order details including coordinates, demand, and the type and load of vehicles in the production base are known.
To clarify the scope of the application of this paper, some assumptions are proposed for the model system: (1) the orders are delivered by the same vehicle along one route, and each customer is served only once; (2) we ignore the customer service time and order processing time; (3) the variable shipping cost is proportional to the driving distance; (4) the traffic condition is good, no blocking or situation will Sustainability 2020, 12, 7934 5 of 17 be met; (5) each vehicle delivers orders from the tomato production base and returns to it immediately after distribution; (6) the rule of tomato picking is that one vehicle picks another vehicle after it has been picked, and each vehicle will be used as transport at once after it has been picked. To clarify the scope of the application of this paper, some assumptions are proposed for the model system: (1) the orders are delivered by the same vehicle along one route, and each customer is served only once; (2) we ignore the customer service time and order processing time; (3) the variable shipping cost is proportional to the driving distance; (4) the traffic condition is good, no blocking or situation will be met; (5) each vehicle delivers orders from the tomato production base and returns to it immediately after distribution; (6) the rule of tomato picking is that one vehicle picks another vehicle after it has been picked, and each vehicle will be used as transport at once after it has been picked.
The main notations are shown in Table 2. The decisions to be taken in this model are as follows: (1) how many of each type of vehicle should be selected; (2) the optimal picking sequence and corresponding tomato picking maturity; (3) which order is delivered by which type of vehicle to the consumers and the optimal served sequence, with the objective of minimizing the total cost, including distribution cost and discipline cost.   The main notations are shown in Table 2. The decisions to be taken in this model are as follows: (1) how many of each type of vehicle should be selected; (2) the optimal picking sequence and corresponding tomato picking maturity; (3) which order is delivered by which type of vehicle to the consumers and the optimal served sequence, with the objective of minimizing the total cost, including distribution cost and discipline cost.

Objective Function
In the integrated tomato picking and distribution model, we take two different objectives into consideration, i.e., the distribution cost and the discipline cost generated by time windows.
(1) Distribution Cost Let F m k denote the fixed cost of vehicle k of m, and x m ojk is a binary decision variable, taking the value of 1 if vehicle k of m delivers products from the tomato production base to customer j, and 0 otherwise. Let c m ij denote the unit shipping cost from node i (tomato production base or customers) to node j (tomato production base or customers) using vehicle m. d ij is the distance between node i and node j.
(2) Discipline Cost We consider time windows in this study to satisfy customer demand. Discipline cost can be divided into two parts. One is the part where vehicles arrive before time windows and the other is after time windows. Both of the two parts have different impacts on customers, hence the use of two parameters, α and β, to denote the different extent of the time. The calculation is given in

Constraints
We formulated a multiple-vehicle routing problem with time windows, and the constraints are given as follows: Sustainability 2020, 12, 7934 Constraint (4) limits the load of each vehicle, which cannot exceed its rated capacity. Constraint (5) indicates that the number of vehicles selected cannot exceed the capacity of the base. Constraint (6) denotes that subloops are eliminated. Constraint (7) limits each vehicle departing and returning to the tomato production base to, at most, one time. Constraint (8) points out flow conservation, which means once the vehicles arrive at the required location and finish the service, it will leave immediately. Constraints (9) and (10) indicate that there is only one vehicle on each route. Constraint (11) denotes that each customer is served only by one vehicle and only once. Constraint (12) constrains the arc traffic of vehicles. Constraint (13) points out the time relationship between one node and the subsequent nodes. Constraint (13) notes that the time when the vehicle reaches the next node is later than the time when it reaches the current node, and M is a positive number with a large value. Constraint (14) indicates the final time when all the tomato orders that the vehicle needs to distribute are picked. Constraint (15) denotes that the delivery can only start after finishing all the picking. Constraint (16) calculates the distribution time to the customer, ignoring service time and order processing time. Constraint (17) indicates that the tomato maturity delivered to customers is optimal. Constraints (18)- (20) define the conditions on the decision variables.

Programming Model
According to the discussion in Sections 3.2 and 3.3, the total objective function is to minimize obj 1 and obj 2 by subjecting them to Constraints (4)- (20). The integrated tomato picking and distribution scheduling model can be expressed as min ob j1 + ob j2 (21) Note that Model (21) is a multiple-vehicle routing problem with time windows (m-VRPTW), and it has proven to be an NP-hard problem [27]. If the tomato production base only has one type of vehicle, Model (21) is equivalent to the VRPTW problem. If Constraints (13)-(18) are removed and ET i = 0, LT i = M, and only one type of vehicle is available, Model (21) is equivalent to the VRP problem.

Model Solution Based on Sweep-Adaptive Genetic Algorithm
The genetic algorithm (GA), as a heuristic algorithm, is often used to solve VRP problems [28][29][30]. It simulates the evolution mechanism in nature to search for the optimal solution. GA starts from the initial population and continuously changes the population chromosomes by mutual iteration to seek the approximate optimal solution [31]. It has some disadvantages, such as high time complexity and slow convergence speed. However, the sweep algorithm can divide problems into several parts for parallel operation, which effectively improves the search speed. The self-adaptive algorithm plays an important role in improving convergence accuracy and accelerating the convergence speed. Based on the advantages of the above algorithms, this study proposes a sweep-adaptive genetic algorithm to solve the integrated tomato picking and distribution model.
The sweep-adaptive genetic algorithm (S-AGA) is a hybrid heuristic algorithm based on the genetic algorithm. This study formulates corresponding criteria and forms multiple effective initial solutions according to constraints. Then, the optimal solution is reached through selection, adaptive crossover, variation, and termination conditions, minimizing the total distribution and discipline costs. The pseudocode of the algorithm is shown in Algorithm 1 below.
Mutation for the current individual end if end while } until (Fitness less than T f or generation above G)

Code and Population Initialization
The feasible solution is expressed by the chromosome in GA. Additionlly, we use the integer coding method to encode the chromosome; 0 represents the tomato production base, and the remaining nodes are added with 1 based on the original serial number. For example, chromosome S = 174,216,351. As shown in chromosome S, there are two routes, including Customers 7, 4, and 2 served in Route 1 and Customers 6, 3, and 5 served in Route 2. They are all served successively.
The initial population is the origin of evolution, which has an impact on the operation performance of GA so that the population size should be carefully selected. In this paper, we divide the population into several parallel subpopulations. Each subpopulation represents the feasible solution of each vehicle, being a part of the optimal solution of the optimization model. There are many ways of initialization, such as chaotic initialization and random initialization. Random initialization is selected in this study.

Constraint Treatment
Traditional GA has a high time-complexity. It is not suitable for large-scale studies. The sweep algorithm is used in this study to reduce the difficulty of model solving, preprocessing data, and meeting the constraints of load and number of vehicles via grouping first and processing later. It is shown in Algorithm 2 from Step 1 to Step 5, translating the VRP problem into several simple TSP problems.

Algorithm 2: Steps of the sweep algorithm
Step1: Sort vehicle capacity from large to small. Q m denotes vehicle capacity of m. Sorted capacity is as follows: Q 1 , Q 2 , · · · , Q m ; Step2: Establish a coordinate system with the tomato production base as the origin. Centered on the production base and bounded by the abscissa, the ray sweeps all customer nodes. i denotes the customers according to the sweep sequence. d i denotes the customer demand. Then, calculate the interpolation ∆ im = Q m − d i between the remaining capacity of each vehicle and the customer nodes; Step3: Continue to sweep the customer nodes until the total demand of the sweep customer nodes is larger than the capacity of vehicle m, that is ∆ im < 0; Step4: Record the former customers meeting constraints ∆ im < 0 and ∆ (i−1)m < 0. Select the minimum positive value or zero value so we can get the corresponding customers. Allocate customer node i and its former nodes to the corresponding type of vehicles through constraint ∆ (i−1)m > 0. If this type of vehicle is allocated completely, select the second small positive value; Step5: Take customer i as the new starting point and return to Step2 for the next round of sweep.

Fitness and Selection
Fitness can be used to evaluate the quality of a chromosome. Chromosomes of high quality are easy to survive and exist on a large scale, while those of poor quality will exist on a small scale or even be extinct. The objective function is to minimize the total cost. To get the optimal solution, we use the reciprocal of the objective function to denote the fitness, as shown in Equation (23). Z i indicates the objective function value of chromosome i.
The selection calculation is to determine the probability of being selected according to the size of the fitness. Higher fitness will be better selected. The next generation of surplus individuals is calculated by the proportional selection method.

Crossover and Mutation
Crossover and mutation operation can generate new individuals and increase population diversity. In this study, a sequential crossover is used, and a self-adaptive crossover is added to improve population diversity [32], as shown in Equation (24). The crossover operation of this paper is order crossover. For example, Chromosome A = 2 8 4 10 5 1 7 3 6 9 or Chromosome B = 5 6 7 1 10 2 8 3 9 4. First, select a gene sequence with the same position in Chromosomes A and B randomly, such as the 3 to 6 genes in Chromosomes A (_ _4 10 5 1 _ _ _ _) and B (_ _ 7 1 10 2 _ _ _ _). Then, keep the proposed gene sequence unchanged, exchange the remaining gene sequences of Chromosomes A and B in the original order, that is, the new Chromosomes A1 = 6 7 4 10 5 1 2 8 3 9 and B1 = 8 4 7 1 10 2 5 3 6 9.
where f and f denote the larger fitness of two individuals by crossover or mutation operations; k 1 , k 2 , k 3 , k 4 are constants; f max indicates the largest fitness of the population; f avg denotes the average fitness of the population.

Termination Condition
The termination condition of this paper achieves a predetermined set of evolutionary generations when iterations reach 1000 or the fitness value between successive generations is less than the threshold value. Otherwise, the above steps will be repeated until the termination conditions are met. When evolution is stopped, select the chromosome with the highest fitness in each subpopulation as the corresponding optimal scheme. All the schemes are combined together to get the overall scheme design of the problem. The overall process is shown in Figure 3.  (25) where f and f ' denote the larger fitness of two individuals by crossover or mutation operations; k 1 , k 2 , k 3 , k 4 are constants; f max indicates the largest fitness of the population; f avg denotes the average fitness of the population.

Termination Condition
The termination condition of this paper achieves a predetermined set of evolutionary generations when iterations reach 1000 or the fitness value between successive generations is less than the threshold value. Otherwise, the above steps will be repeated until the termination conditions are met. When evolution is stopped, select the chromosome with the highest fitness in each subpopulation as the corresponding optimal scheme. All the schemes are combined together to get the overall scheme design of the problem. The overall process is shown in Figure 3.

Numerical Experiments and Simulation Analysis
To verify the effectiveness of the model, this study used S-AGA and traditional GA for a simulation operation aimed at different scales and compared the results. In this section, we adopt MATLAB R2017b as an MINLP solver to solve the integrated tomato picking and distribution model. In Section 5.1, a small-scale numerical experiment with one tomato production base and 20 customers is carried out to prove the efficiency and feasibility of the proposed model. In Section 5.2, we address 11 different large-scale numerical experiments to verify the superiority of S-AGA compared with traditional GA.
The tomato production base has two types of vehicles. The tomato picking speed of the production base is 1 t/h, and the change of tomato maturity after picking is consistent with the proposed tomato maturity model. According to practical experience, the optimal tomato maturity is 9, corresponding to firmness 29N. Considering daily life, customers are more likely to receive the agriproducts early. Therefore, we set α = 0 and β = 20 CNY/h. At the same time, we set related parameters of algorithm S-AGA, including the initial population of 50, maximum iterations of 500, and self-adaptive crossover and mutation parameters k 1 = 1.0, k 2 = 0.5, k 3 = 0.8, k 4 = 0.5. Table 3 demonstrates the details of 20 orders. Table 4 gives the vehicle information.

Numerical Results and Management Significance
The basic numerical experiment is tested 10 times, and the optimal scheduling of the tomato supply chain is reported in Table 5, including vehicle selection, tomato picking sequence, picking maturity, and firmness, as well as the optimal route schedule. Tomato picking maturity can be obtained by the corresponding table of firmness and maturity shown in Figure 1. The firmness interval (31,35) corresponds to Maturity 8, while the firmness interval (27,31) corresponds to Maturity 9. We take the order of the table column as the picking order. The 20 customers are represented on the distance axis in the form of a hollow circular in Figure 4. Blue lines denote the route of Vehicle A, while red lines denote the route of Vehicle B, depicted in Figure 4a. The evolution graph of S-AGA is depicted in Figure 4b. It can be seen that the proposed model can provide a complete scheme of the tomato supply chain, improving the overall operational efficiency of logistics enterprises and the product quality service.
Cost (CNY) 6134.98 4892.05 the distance axis in the form of a hollow circular in Figure 4. Blue lines denote the route of Vehicle A, while red lines denote the route of Vehicle B, depicted in Figure 4a. The evolution graph of S-AGA is depicted in Figure 4b. It can be seen that the proposed model can provide a complete scheme of the tomato supply chain, improving the overall operational efficiency of logistics enterprises and the product quality service. In Figure 4a, we obviously find that the route of Vehicle B is not the shortest because we should consider time windows of customers and weigh the cost of logistics against the discipline cost of tomatoes not being delivered within the request time windows in real life. If we pay more logistics costs than discipline costs, we may ignore customer satisfaction to minimize the whole cost. However, if we pay more discipline costs compared with logistics costs, enterprises would like to take the longer way to meet time windows and customer satisfaction. In fact, in order to minimize the overall tomato supply chain cost, we should consider a reasonable number and type of vehicles and routes.
6134.98 4892.05  In Figure 4a, we obviously find that the route of Vehicle B is not the shortest because we should consider time windows of customers and weigh the cost of logistics against the discipline cost of tomatoes not being delivered within the request time windows in real life. If we pay more logistics costs than discipline costs, we may ignore customer satisfaction to minimize the whole cost. However, if we pay more discipline costs compared with logistics costs, enterprises would like to take the longer way to meet time windows and customer satisfaction. In fact, in order to minimize the overall tomato supply chain cost, we should consider a reasonable number and type of vehicles and routes.
We ignore the time of tomato picking and only consider the deterioration in the process of transportation. The new optimal scheme of the tomato supply chain is reported in Table 6, and new optimal routes and the evolution graph are described in Figure 5. It can be seen that part of the firmness and maturity of tomatoes change when they are delivered.
The total cost is reduced by 1%, while customer satisfaction based on tomato quality decreases from 100% to 95%. The example in this section is on a small scale; the change of tomato maturity is not significant for short-distance distribution. However, in modern society, with convenient transportation, delivery areas of suppliers are larger. Moreover, for some agricultural products like juicy peach, whose maturity changes faster, the quality change of agricultural products will be more obvious, and the superiority of this model will be further highlighted.

Large-Scale Example Setting and Performance Analysis of Algorithm
In this section, we use 11 different large-scale numerical experiments to demonstrate the advantage and efficiency of S-AGA. Customer size increases from 50 to 1000. The computational performance results of S-AGA and GA are reported via CPU time and iteration number.
We generate N customers randomly. Furthermore, among them, 70% are online customers, while 30% are offline customers. The online customer demand D i follows the uniform distribution U (0.005, 0.015), while the offline customer demand D i follows the uniform distribution U (0.6, 2). The earliest arrival time ET i is generated by U (0, 24), and the latest arrival time LT i is generated by ET i +U (15,50).
We reported the computational results of GA and S-AGA in Table 7. It can be seen that with the increase in customer scale, GA and S-AGA show an increasing trend in CPU time and iterations. However, the time gap of S-AGA is much faster than that of GA, which is more than 90%. Comparing Table 6 with Table 5, we can find that the cost of both schemes is almost the same. The total cost is reduced by 1%, while customer satisfaction based on tomato quality decreases from 100% to 95%. The example in this section is on a small scale; the change of tomato maturity is not significant for short-distance distribution. However, in modern society, with convenient transportation, delivery areas of suppliers are larger. Moreover, for some agricultural products like juicy peach, whose maturity changes faster, the quality change of agricultural products will be more obvious, and the superiority of this model will be further highlighted.

Large-Scale Example Setting and Performance Analysis of Algorithm
In this section, we use 11 different large-scale numerical experiments to demonstrate the advantage and efficiency of S-AGA. Customer size increases from 50 to 1000. The computational performance results of S-AGA and GA are reported via CPU time and iteration number.
We generate N customers randomly. Furthermore, among them, 70% are online customers, while 30% are offline customers. The online customer demand D i follows the uniform distribution U (0.005, 0.015), while the offline customer demand D i follows the uniform distribution U (0.6, 2). The earliest arrival time ET i is generated by U (0, 24), and the latest arrival time LT i is generated by ET i +U (15,50).
We reported the computational results of GA and S-AGA in Table 7. It can be seen that with the increase in customer scale, GA and S-AGA show an increasing trend in CPU time and iterations. However, the time gap of S-AGA is much faster than that of GA, which is more than 90%. Additionally, S-AGA shows more advantages of iterations, reducing more than 50% compared with GA.  Figure 6a demonstrates the comparison of CPU time between S-AGA and GA. It shows that S-AGA is significantly superior to GA in CPU time. Figure 6b shows the trend of iterations of S-AGA and GA. From Case 1 to Case 7, the iteration gaps become smaller, while from Case 8 to Case 9, the iteration gaps get larger compared with Case 7. However, we find that the iterations of GA are all larger than that of S-AGA, and the largest iteration gap between S-AGA and GA is 67%. Generally, the performance of S-AGA proposed in this paper is better than the traditional GA in solving integrated tomato picking and distribution. Moreover, we conclude that S-AGA has good convergence and faster convergence speed, and it is suitable for providing schemes quickly when dealing with large-scale customers. Additionally, S-AGA shows more advantages of iterations, reducing more than 50% compared with GA.  Figure 6a demonstrates the comparison of CPU time between S-AGA and GA. It shows that S-AGA is significantly superior to GA in CPU time. Figure 6b shows the trend of iterations of S-AGA and GA. From Case 1 to Case 7, the iteration gaps become smaller, while from Case 8 to Case 9, the iteration gaps get larger compared with Case 7. However, we find that the iterations of GA are all larger than that of S-AGA, and the largest iteration gap between S-AGA and GA is 67%. Generally, the performance of S-AGA proposed in this paper is better than the traditional GA in solving integrated tomato picking and distribution. Moreover, we conclude that S-AGA has good convergence and faster convergence speed, and it is suitable for providing schemes quickly when dealing with large-scale customers.

Conclusions
This paper studied an integrated tomato picking and distribution schedule based on tomato maturity in the context of a farm-to-door model. We presented a tomato maturity model with firmness and incorporated it into the proposed integrated tomato picking and distribution model, with the objectives of minimizing the total cost in terms of obtaining optimal tomato quality and satisfaction for customers. In the proposed model, we addressed an improved algorithm, S-AGA, to solve the model and get the schedule of the tomato supply chain. In addition, numerical experiments were conducted to demonstrate the efficiency and effectiveness of S-AGA in solving the proposed

Conclusions
This paper studied an integrated tomato picking and distribution schedule based on tomato maturity in the context of a farm-to-door model. We presented a tomato maturity model with firmness and incorporated it into the proposed integrated tomato picking and distribution model, with the objectives of minimizing the total cost in terms of obtaining optimal tomato quality and satisfaction for customers. In the proposed model, we addressed an improved algorithm, S-AGA, to solve the model and get the schedule of the tomato supply chain. In addition, numerical experiments were conducted to demonstrate the efficiency and effectiveness of S-AGA in solving the proposed model compared with traditional GA. The simulation result of a small-scale example shows that we can provide tomato quality and improve by over 5% customer satisfaction without additional cost. Moreover, the result of a large-scale example demonstrates that both CPU time and iteration of S-AGA are superior to GA. Therefore, with the analysis of the example, the designed algorithm and proposed model established in this paper are reasonable. The promotion of this model is also conducive to realizing the personalized customization of customers and improving the sustainability of agriculture and logistics. Users can choose different tomato maturities to store, sell, or eat, meeting the needs of diversified customers.
There are specific areas in this study that can be researched further. This study only considered a single tomato production base; hence multiple bases could be considered as a prospect. Moreover, tomato maturity only considered one indicator, so we could consider multiple indicators to propose in the tomato maturity model. Then, another prospect could incorporate more uncertainty factors into the integrated tomato picking and distribution schedule, such as the uncertainties of transportation and the processing of tomato picking. Finally, our research is based on the context of the farm-to-door model, so another model and related factors could be taken into consideration, such as warehouse storage time. We leave these problems for future work.