CN112001678A - Method for reducing tobacco distribution mileage - Google Patents
Method for reducing tobacco distribution mileage Download PDFInfo
- Publication number
- CN112001678A CN112001678A CN202010870443.1A CN202010870443A CN112001678A CN 112001678 A CN112001678 A CN 112001678A CN 202010870443 A CN202010870443 A CN 202010870443A CN 112001678 A CN112001678 A CN 112001678A
- Authority
- CN
- China
- Prior art keywords
- solution
- delivery
- algorithm
- new
- current
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 241000208125 Nicotiana Species 0.000 title claims abstract description 21
- 235000002637 Nicotiana tabacum Nutrition 0.000 title claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims description 72
- 238000002922 simulated annealing Methods 0.000 claims description 36
- 238000010845 search algorithm Methods 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 4
- 230000008707 rearrangement Effects 0.000 abstract 1
- 235000019504 cigarettes Nutrition 0.000 description 11
- 238000005457 optimization Methods 0.000 description 9
- 238000011161 development Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000003442 weekly effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000003064 k means clustering Methods 0.000 description 3
- 238000010187 selection method Methods 0.000 description 3
- 239000000779 smoke Substances 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000000003 hoof Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 235000019505 tobacco product Nutrition 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
- G06Q10/083—Shipping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Entrepreneurship & Innovation (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method for reducing tobacco delivery mileage, which comprises the following steps: step S1, adjusting the delivery structure, namely acquiring the number of vehicle delivery times in a week and acquiring the number of delivery vehicles; step S2, constructing a minimum delivery unit of the client; step S3, rearrangement of the delivery route by the minimum delivery unit of the customer, the present invention can effectively reduce the repeat delivery distance.
Description
Technical Field
The invention relates to the technical field of tobacco company distribution, in particular to a method for reducing tobacco distribution mileage.
Background
In recent years, a tobacco business system of the whole province guides the high-quality development hoof to stably advance by a new development concept, systematically performs source opening, submergence digging and throttling, becomes an important way for continuously improving the operation benefit and promoting the operation, quality and efficiency improvement of enterprises, and as an industry logistics enterprise of important links of cost reduction and efficiency improvement in the tobacco business system, the development idea of provincial bureaus should be implemented in the cigarette distribution main industry and the new non-cigarette distribution state, the delivery link with the largest cost occupation ratio is focused, the artificial intelligence algorithm is introduced, the route planning is more reasonable, the effective integration of operation resources is promoted, the total mileage of non-cigarette integrated distribution of cigarettes is continuously reduced, the distribution cost is saved, and the high-quality development of the enterprises with the good quality is realized.
The existing logistics transportation capacity can only meet a part of non-smoke distribution requirements; namely, third-party logistics distribution influences the gross profit rate of non-tobacco products, and according to the win-win goal, how to effectively control the total logistics distribution mileage and simultaneously considering the development of the non-tobacco distribution is a new challenge of cigarette non-tobacco integration lower hem in front of logistics companies. The original delivery route arrangement of the tobacco company has the big problems that the total delivery mileage is relatively long, the delivery quantity is relatively small, and how to reduce the total delivery mileage becomes a problem to be solved under the comprehensive development of cigarette non-smoke integrated delivery work.
Disclosure of Invention
In order to overcome the problems, the invention aims to provide a method for reducing the tobacco delivery mileage, which effectively reduces the delivery mileage and improves the delivery efficiency.
The invention is realized by adopting the following scheme: a method of reducing tobacco delivery mileage, the method comprising the steps of:
step S1, adjusting the delivery structure, namely acquiring the number of vehicle delivery times in a week and acquiring the number of delivery vehicles;
step S2, constructing a minimum delivery unit of the client;
step S3, the delivery route is rearranged by the customer minimum delivery unit.
Further, the step S2 is further specifically: dividing the customers into a plurality of minimum distribution units according to the similarity degree of the longitude and the latitude through a clustering algorithm, and merging the similar customers to form respective minimum distribution units; set the customer set to be deliveredV={v1,v2,...,vnThe longitude and latitude attribute of the client is xiConstructing k clustering subsets { c1,c2,...,ck-minimizing the client differentiation in k subsets of clusters; namely, it isWherein,for clustering subset { c1,c2,...,ckMiddle cluster cj(j ═ 1, 2.. times, k) cluster center, viRepresenting a cluster cjAll clients in; the selection method of the cluster center is as follows:wherein, | cjI is the number of customers in cluster j, xjDenotes cjLongitude and latitude data of the middle client;
the clustering algorithm takes the Euclidean distance as the similarity of the client nodes, so the Euclidean distance d (x) of any two client nodesi,xj) Expressed as:
further, the step S3 is further specifically: step 31: obtaining a minimum delivery unit of a client;
step 32: setting initialization parameters, i.e. setting the initial temperature T0Termination temperature Tf;
Step 33: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the simulated annealing algorithm is obtained after continuous adjustment;
step 34: by local search to obtain a new solution, i.e. by making a one for the current optimum pointRandomly varied to generate a new solution XkAcquiring an objective function f (X) of a new simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the objective function0) Let k be k +1, k being the number of iterations;
step 35: judging whether the new solution can become the current optimal solution, namely if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 36: judging whether the simulated annealing algorithm is finished or not, and if the current temperature is more than the termination temperature TfReturning to step 34, otherwise, outputting the current solution as the optimal solution, and ending the algorithm.
Further, the method further comprises: step S4, optimizing the delivery route through a tabu search algorithm; the step S4 further includes: based on the minimum delivery units of the customers, defining the positions of the customer nodes as the average longitude and latitude information of each minimum delivery unit, judging whether the distance between any two customer nodes is smaller than a preset Q value or not through a proximity KNN algorithm, if so, performing two neighborhood searching operations, and recording an accessed neighborhood solution; if not, abandoning the pair of nodes and carrying out the same operation on other nodes.
Further, the step S4 is further specifically: step 41: obtaining a minimum delivery unit of a client;
step 42: setting initialization parameters, setting initial temperature T0Termination temperature TfAnd the value of G in the KNN algorithm; discarding the farthest nodes which do not conform to the KNN algorithm through a preset G value, and reducing the range of neighborhood searching;
step 43: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the hybrid simulated annealing algorithm is obtained after continuous adjustment;
step 44: by searching locally to obtain a new solution, i.e. at the current solution X0Respectively carrying out local search operations of an exchange method and an interpolation method in the neighborhood of the target to generate a new feasible solution X, obtaining an objective function f (X) of a new simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the objective function0) Let k be k +1, k being the number of iterations;
step 45: judging whether the new solution can become the current optimal solution, if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 46: judging whether all neighborhood searching is finished, namely under the current temperature T, if not, turning to the step 47, otherwise, turning to the step 44;
step 47: judging whether the hybrid simulated annealing algorithm is finished or not, and if the current temperature is more than the termination temperature TfGo back to step 44, otherwise output the current solution as the optimal solution, and end the algorithm.
The invention has the beneficial effects that: the bottleneck problem of manual line planning which is troubled for a long time is solved, a cigarette non-cigarette integrated 4D distribution new mode based on the ant colony working principle is preliminarily explored and formed, intelligent optimization model improvement methods such as a hybrid simulation annealing algorithm covering the minimum delivery unit and based on tabu search are built, the delivery mileage is effectively reduced, the delivery efficiency is improved, the characteristics of strong reproducibility and generalizability are achieved in the industry, and certain effect is achieved in the aspects of operation cost control and informatization construction.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
Fig. 2 is a detailed flowchart of the present invention for optimizing the delivery route by a tabu search algorithm.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1, the present invention provides a method for reducing the mileage of tobacco distribution, the method comprising the steps of:
step S1, adjusting the delivery structure, namely acquiring the number of vehicle delivery times in a week and acquiring the number of delivery vehicles;
step S2, constructing a minimum delivery unit of the client;
the step S2 further includes: dividing the customers into a plurality of minimum distribution units according to the similarity degree of the longitude and the latitude through a clustering algorithm, and merging the similar customers to form respective minimum distribution units; setting the distributed client set as V ═ V1,v2,...,vnThe longitude and latitude attribute of the client is xiConstructing k clustering subsets { c1,c2,...,ck-minimizing the client differentiation in k subsets of clusters; namely, it isWherein,for clustering subset { c1,c2,...,ckMiddle cluster cj(j ═ 1, 2.. times, k) cluster center, viRepresenting a cluster cjAll clients in; the selection method of the cluster center is as follows:wherein, | cjI is the number of customers in cluster j, xjDenotes cjLongitude and latitude data of the middle client;
the clustering algorithm takes the Euclidean distance as the similarity of the client nodes, so the Euclidean distance d (x) of any two client nodesi,xj) Expressed as:
step S3, the delivery route is rearranged by the customer minimum delivery unit.
The step S3 further includes: step 31: obtaining a minimum delivery unit of a client;
step 32: setting initialization parametersNumber, i.e. setting the initial temperature T0Termination temperature Tf;
Step 33: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the simulated annealing algorithm is obtained after continuous adjustment;
step 34: obtaining a new solution by local search, i.e. randomly changing the current optimum point to generate a new solution XkAcquiring an objective function f (X) of a new simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the objective function0) Let k be k +1, k being the number of iterations;
step 35: judging whether the new solution can become the current optimal solution, namely if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 36: judging whether the simulated annealing algorithm is finished or not, and if the current temperature is more than the termination temperature TfReturning to step 34, otherwise, outputting the current solution as the optimal solution, and ending the algorithm.
Referring to fig. 2, the method further includes: step S4, optimizing the delivery route through a tabu search algorithm; the step S4 further includes: based on the minimum delivery units of the customers, defining the positions of the customer nodes as the average longitude and latitude information of each minimum delivery unit, judging whether the distance between any two customer nodes is smaller than a preset Q value or not through a proximity KNN algorithm, if so, performing two neighborhood searching operations, and recording an accessed neighborhood solution; if not, abandoning the pair of nodes and carrying out the same operation on other nodes.
The step S4 further includes: step 41: obtaining a minimum delivery unit of a client;
step 42: setting initialization parameters, setting initial temperature T0Termination temperature TfTo do so byAnd the G value in the KNN algorithm; discarding the farthest nodes which do not conform to the KNN algorithm through a preset G value, and reducing the range of neighborhood searching;
step 43: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the hybrid simulated annealing algorithm is obtained after continuous adjustment;
step 44: by searching locally to obtain a new solution, i.e. at the current solution X0Respectively carrying out local search operations of an exchange method and an interpolation method in the neighborhood of the target to generate a new feasible solution X, obtaining an objective function f (X) of a new simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the objective function0) Let k be k +1, k being the number of iterations;
step 45: judging whether the new solution can become the current optimal solution, if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 46: judging whether all neighborhood searching is finished, namely under the current temperature T, if not, turning to the step 47, otherwise, turning to the step 44;
step 47: and judging whether the hybrid simulated annealing algorithm is finished, if the current temperature is higher than the termination temperature Tf, returning to the step 44, otherwise, outputting the current solution as the optimal solution, and finishing the algorithm.
The invention is further illustrated below with reference to a specific embodiment:
structure for implementing and adjusting delivery of countermeasure
(1) Calculating the number of laps of the turnover vehicle
In 2016 + 2018, the cigarette sales increase is relatively small, and the influence on the vehicle loading rate is mainly the number of trips. The average weekly sales is calculated by calling the total amount of the year 2018 in the whole year, and the weekly trips meeting the 75% loading rate index are calculated by taking the full load of 95 pieces and the provincial load rate index of 75% as basic parameters.
Calculating the formula: the number of vehicle-out passes is 75% of the average round-trip-rate/vehicle loading rate/95 of the number of single-vehicle full-load pieces. As shown in Table 1 below
TABLE 1
It is found by calculation that the weekly vehicle trips are controlled within 187 trips when the annual vehicle loading rate is 75%.
(2) Measuring and calculating distribution vehicle
After the delivery service data is accurately analyzed, the number of vehicles required to be delivered is determined by taking 187 weeks as reference conditions, the number of the vehicle passes and the proportion of clients on the day are respectively calculated according to 26, 27 and 28 vehicles, and the calculation results are as follows:
calculating the formula: the departure times are 4 days (tuesday to friday) on the day line, 5 days (tuesday to saturday) on the day line, 1 day (saturday) on the day delivery line
Customer proportion on the day: current day line, number of delivery customers (press 65 customers) day/total number of customers
Through data analysis, of the vehicle plans around the weekly trip 187, the daily delivery ratio of 27 vehicle plans and 26 vehicle plans is close to 30%, and the daily delivery customer ratio of 28 vehicle plans is only 22%. So the proposal of 28 vehicles is eliminated.
Through measurement and calculation, 27 distribution vehicles are adopted (namely one vehicle is reduced on the basis of the original delivery vehicle), the average daily average number of users of a single vehicle is increased by 4, the single vehicle distribution mileage is increased by about 8 kilometers, the total single vehicle distribution time is estimated to be increased by about 0.5 hour, the total average single vehicle distribution time is estimated to be 7.4 hours, the distribution time of a C25 line is only met, and the working time is 8.22 hours after 8 hours.
By adopting 26 delivery vehicles (namely, two vehicles are reduced on the basis of the original delivery vehicles), the average daily average number of users of a single vehicle is increased by 8, the delivery mileage of the single vehicle is increased by about 16 kilometers, the total delivery time of the single vehicle is expected to be increased by about 1 hour, the total delivery time of the single vehicle is expected to be 7.9 hours, and the working time of 13 delivery lines exceeds 8 hours.
By contrast, it is considered relatively reasonable to use 27 cars for delivery.
And implementing a second strategy: customer minimization unit
And dividing the customers into a plurality of minimum distribution units according to the similarity of the longitude and the latitude through a clustering algorithm, and merging the customers close to the city compartment customer service part to form respective minimum distribution units.
The K-means clustering algorithm for the pu tian tobacco stream can be described as: given a set V ═ V { V } containing 1420 city block customers1,v2,...,vn1420, where the latitude and longitude attribute of the client is xiConstructing k clustering subsets { c1,c2,...,ckAnd (4) minimizing the customer differentiation in the k clustering subsets. The objective function of the K-means clustering algorithm is shown as formula (9), and the right side of the objective function is also the error sum of squares function of the class.
Namely, it isWherein,for clustering subset { c1,c2,...,ckMiddle cluster cj(j ═ 1, 2.. times, k) cluster center, vjRepresenting a cluster cjAll clients in; the selection method of the cluster center is formula (10):wherein, | cjI is the number of customers in cluster j, xjDenotes cjLongitude and latitude data of the middle client;
the K-means clustering algorithm of the Pu Tian tobacco takes the Euclidean distance as the similarity of the client nodes, so the Euclidean distance of any two client nodesFormula d (x)i,xj) Expressed as:
the division of the minimum distribution unit of the urban area service part is schematically shown in the following table 2.
TABLE 2 minimum delivery Unit schematic of urban service department
The purpose of the customer base being divided into a number of minimum delivery units is to reduce the pressure on the optimization of the line, the number of minimum delivery units should not be too large. Meanwhile, if the number of customers in the minimum delivery unit is too large, the accuracy of the route optimization is affected, and therefore, the number of the minimum delivery unit is not too small.
And (3) implementing a strategy III: rearranging a delivery route:
the route planning of the present invention is done manually. The artificial line optimization can be reluctantly completed in mountain service departments with sparser customer distribution, but for urban service departments with dense customers and large demand, the artificial line optimization takes time and has poor effect. The invention is based on the actual situation of the Pu Tian tobacco material flow, applies a simulated annealing algorithm, and comprises the following specific steps:
step 31: obtaining a minimum delivery unit of a client;
step 32: setting initialization parameters, i.e. setting the initial temperature T0Termination temperature Tf;
Step 33: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining the objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the simulated annealing algorithm is obtained after continuous adjustment;
step 34: by local search to obtain new solutions, i.e. for the current optimumThe point is randomly changed to generate a new solution XkAcquiring a new target function f (X) of the simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the target function0) Let k be k +1, k being the number of iterations;
step 35: judging whether the new solution can become the current optimal solution, namely if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 36: judging whether the simulated annealing algorithm is finished or not, and if the current temperature is more than the termination temperature TfReturning to step 34, otherwise, outputting the current solution as the optimal solution, and ending the algorithm.
Through the intelligent line optimization of the subject group, nearly twenty thousand cigarettes are distributed by the service department of the urban compartment area every week, and account for 80.8 percent of all orders of the urban compartment area. After optimization, the number of customers 289 is sent to the current day, which accounts for 20.3% of the total number of customers in the urban area, and is nearly twice of that of the manual optimization line, so that the customer experience is greatly improved.
And implementing a strategy of four: taboo search optimized delivery line
According to the method, the distances from the client nodes to all the client nodes in the training set are calculated, K adjacent client nodes closest to the current client node are selected, and a KNN algorithm is adopted and combined with two neighborhood search algorithms, so that the calculation complexity of local search is reduced, the complexity of a hybrid simulated annealing algorithm is reduced, and the operation efficiency of the algorithm is improved.
Based on the minimum delivery units proposed above, the client node location is defined as the average longitude and latitude information of each minimum delivery unit, and there are 223 client nodes taking the urban area as an example. Judging whether the distance between any two nodes is smaller than a preset G value or not by a neighbor KNN algorithm for any two nodes, if so, performing two neighborhood searching operations, and recording an accessed neighborhood solution; if not, abandoning the pair of nodes and carrying out the same operation on other nodes. And performing secondary promotion on the basis of the re-simulation algorithm, and establishing a hybrid simulated annealing algorithm based on tabu search, wherein the specific steps are as follows, and a flow chart is shown in fig. 2.
The step S4 further includes: step 41: obtaining a minimum delivery unit of a client;
step 42: setting initialization parameters, setting initial temperature T0Termination temperature TfAnd the value of G in the KNN algorithm; discarding the farthest nodes which do not conform to the KNN algorithm through a preset G value, and reducing the range of neighborhood searching;
step 43: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Calculating the objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the hybrid simulated annealing algorithm is obtained after continuous adjustment;
step 44: by searching locally to obtain a new solution, i.e. at the current solution X0Respectively carrying out local search operations of an exchange method and an interpolation method in the neighborhood of the target to generate a new feasible solution X, calculating a new target function f (X) of the simulated annealing algorithm, and calculating the increment delta f (f) (X) -f (X) of the target function0) Let k be k +1, k being the number of iterations;
step 45: judging whether the new solution can become the current optimal solution, if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution according to the probability p ═ exp (-delta f/T);
step 46: judging whether all neighborhood searches are finished, namely, under the current temperature T, if all neighborhood searches in a tabu table in the KNN algorithm are not finished, turning to the step 47, and otherwise, turning to the step 44;
step 47: and judging whether the hybrid simulated annealing algorithm is finished, if the current temperature is higher than the termination temperature Tf, returning to the step 44, otherwise, outputting the current solution as the optimal solution, and finishing the algorithm.
In a word, the bottleneck problem of manual line planning which is troubled for a long time is solved, a new cigarette non-smoke integration type 4D distribution mode based on the ant colony working principle is preliminarily explored and formed, intelligent optimization model improvement methods such as a mixed simulated annealing algorithm covering the minimum delivery unit of delivery and based on tabu search are built, the delivery mileage is effectively reduced, the delivery efficiency is improved, the characteristics of strong reproducibility and popularization are achieved in the industry, and certain effect is achieved in the aspects of operation cost control and information construction.
The above description is only a preferred embodiment of the present invention, and all equivalent changes and modifications made in accordance with the claims of the present invention should be covered by the present invention.
Claims (5)
1. A method of reducing tobacco delivery mileage, comprising: the method comprises the following steps:
step S1, adjusting the delivery structure, namely acquiring the number of vehicle delivery times in a week and acquiring the number of delivery vehicles;
step S2, constructing a minimum delivery unit of the client;
step S3, the delivery route is rearranged by the customer minimum delivery unit.
2. A method of reducing tobacco delivery mileage as recited in claim 1, wherein: the step S2 further includes: dividing the customers into a plurality of minimum distribution units according to the similarity degree of the longitude and the latitude through a clustering algorithm, and merging the similar customers to form respective minimum distribution units; setting the distributed client set as V ═ V1,v2,...,vnThe longitude and latitude attribute of the client is xiConstructing k clustering subsets { c1,c2,...,ck-minimizing the client differentiation in k subsets of clusters; namely, it isWherein,for clustering subset { c1,c2,...,ckMiddle cluster cj(j ═ 1, 2.. times, k) cluster center, viRepresenting a cluster cjAll clients in; the cluster center is selected byThe formula:wherein, | cjI is the number of customers in cluster j, xjDenotes cjLongitude and latitude data of the middle client;
3. a method of reducing tobacco delivery mileage as recited in claim 1, wherein: the step S3 further includes: step 31: obtaining a minimum delivery unit of a client;
step 32: setting initialization parameters, i.e. setting the initial temperature T0Termination temperature Tf;
Step 33: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the simulated annealing algorithm is obtained after continuous adjustment;
step 34: obtaining a new solution by local search, i.e. randomly changing the current optimum point to generate a new solution XkAcquiring a new target function f (X) of the simulated annealing algorithm, and acquiring the increment delta f (f) (X) -f (X) of the target function0) Let k be k +1, k being the number of iterations;
step 35: judging whether the new solution can become the current optimal solution, namely if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 36: judging whether the simulated annealing algorithm is finished or not, and if the current temperature is more than or equal toTermination temperature TfReturning to step 34, otherwise, outputting the current solution as the optimal solution, and ending the algorithm.
4. A method of reducing tobacco delivery mileage as set forth in claim 2, wherein: the method further comprises the following steps: step S4, optimizing the delivery route through a tabu search algorithm; the step S4 further includes: based on the minimum delivery units of the customers, defining the positions of the customer nodes as the average longitude and latitude information of each minimum delivery unit, judging whether the distance between any two customer nodes is smaller than a preset Q value or not through a proximity KNN algorithm, if so, performing two neighborhood searching operations, and recording an accessed neighborhood solution; if not, abandoning the pair of nodes and carrying out the same operation on other nodes.
5. A method of reducing tobacco delivery mileage as set forth in claim 4, wherein: the step S4 further includes: step 41: obtaining a minimum delivery unit of a client;
step 42: setting initialization parameters, setting initial temperature T0Termination temperature TfAnd the value of G in the KNN algorithm;
discarding the farthest nodes which do not conform to the KNN algorithm through a preset G value, and reducing the range of neighborhood searching;
step 43: generating an initial solution X0I.e. to solve for X initially0As the current optimal solution Xbest=X0Obtaining an objective function f (X) of the simulated annealing algorithm0) (ii) a In the step, the customer node in the sub-path with the maximum order quantity is inserted into the sub-path with the minimum customer order quantity, and an initial feasible solution of the hybrid simulated annealing algorithm is obtained after continuous adjustment;
step 44: by searching locally to obtain a new solution, i.e. at the current solution X0Respectively carrying out local search operations of an exchange method and an interpolation method in the neighborhood of the target to generate a new feasible solution X, obtaining an objective function f (X) of a new simulated annealing algorithm, and obtaining an increment delta f (f) (X) -f (X) of the objective function0) Let k be k +1, k being the number of iterationsCounting;
step 45: judging whether the new solution can become the current optimal solution, if delta f is less than 0, accepting the new solution X to become the current optimal solution, XbestX; otherwise, receiving X as a new current optimal solution by taking the probability p as exp (-delta f/T), wherein T is the temperature T corresponding to the optimal solution X;
step 46: judging whether all neighborhood searching is finished, namely under the current temperature T, if not, turning to the step 47, otherwise, turning to the step 44;
step 47: judging whether the hybrid simulated annealing algorithm is finished or not, and if the current temperature is more than the termination temperature TfGo back to step 44, otherwise output the current solution as the optimal solution, and end the algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010870443.1A CN112001678B (en) | 2020-08-26 | 2020-08-26 | Method for reducing tobacco delivery mileage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010870443.1A CN112001678B (en) | 2020-08-26 | 2020-08-26 | Method for reducing tobacco delivery mileage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112001678A true CN112001678A (en) | 2020-11-27 |
CN112001678B CN112001678B (en) | 2022-07-19 |
Family
ID=73471922
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010870443.1A Active CN112001678B (en) | 2020-08-26 | 2020-08-26 | Method for reducing tobacco delivery mileage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112001678B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000112925A (en) * | 1998-10-07 | 2000-04-21 | Ibm Japan Ltd | Method/system for generating delivery route for miniaturizing number of routes |
JP2010111452A (en) * | 2008-11-04 | 2010-05-20 | Internatl Business Mach Corp <Ibm> | Method, device, and program for obtaining delivery route of a plurality of vehicles |
CN102880798A (en) * | 2012-09-20 | 2013-01-16 | 浪潮电子信息产业股份有限公司 | Variable neighborhood search algorithm for solving multi depot vehicle routing problem with time windows |
CN103383756A (en) * | 2013-07-22 | 2013-11-06 | 浙江省烟草公司绍兴市公司 | Planning method for tobacco logistics distribution routes |
CN107392539A (en) * | 2017-07-23 | 2017-11-24 | 浙江省烟草公司温州市公司 | A kind of cigarette delivery management system |
CN108764754A (en) * | 2018-06-08 | 2018-11-06 | 浪潮软件股份有限公司 | A kind of tobacco distribution vehicle dispatch system based on geographic area grouping |
US20190079975A1 (en) * | 2017-09-11 | 2019-03-14 | Hefei University Of Technology | Scheduling method and system based on hybrid variable neighborhood search and gravitational search algorithm |
US20200017117A1 (en) * | 2018-07-14 | 2020-01-16 | Stephen Milton | Vehicle-data analytics |
CN110930091A (en) * | 2019-11-05 | 2020-03-27 | 哈尔滨工业大学(威海) | Express delivery terminal network point optimization integration method based on neighborhood search simulated annealing algorithm |
CN111178440A (en) * | 2019-12-31 | 2020-05-19 | 哈工汇智(深圳)科技有限公司 | Path optimization method for efficiently reducing total transport cost of AGV |
-
2020
- 2020-08-26 CN CN202010870443.1A patent/CN112001678B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000112925A (en) * | 1998-10-07 | 2000-04-21 | Ibm Japan Ltd | Method/system for generating delivery route for miniaturizing number of routes |
JP2010111452A (en) * | 2008-11-04 | 2010-05-20 | Internatl Business Mach Corp <Ibm> | Method, device, and program for obtaining delivery route of a plurality of vehicles |
CN102880798A (en) * | 2012-09-20 | 2013-01-16 | 浪潮电子信息产业股份有限公司 | Variable neighborhood search algorithm for solving multi depot vehicle routing problem with time windows |
CN103383756A (en) * | 2013-07-22 | 2013-11-06 | 浙江省烟草公司绍兴市公司 | Planning method for tobacco logistics distribution routes |
CN107392539A (en) * | 2017-07-23 | 2017-11-24 | 浙江省烟草公司温州市公司 | A kind of cigarette delivery management system |
US20190079975A1 (en) * | 2017-09-11 | 2019-03-14 | Hefei University Of Technology | Scheduling method and system based on hybrid variable neighborhood search and gravitational search algorithm |
CN108764754A (en) * | 2018-06-08 | 2018-11-06 | 浪潮软件股份有限公司 | A kind of tobacco distribution vehicle dispatch system based on geographic area grouping |
US20200017117A1 (en) * | 2018-07-14 | 2020-01-16 | Stephen Milton | Vehicle-data analytics |
CN110930091A (en) * | 2019-11-05 | 2020-03-27 | 哈尔滨工业大学(威海) | Express delivery terminal network point optimization integration method based on neighborhood search simulated annealing algorithm |
CN111178440A (en) * | 2019-12-31 | 2020-05-19 | 哈工汇智(深圳)科技有限公司 | Path optimization method for efficiently reducing total transport cost of AGV |
Also Published As
Publication number | Publication date |
---|---|
CN112001678B (en) | 2022-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106980912B (en) | Multipoint distribution line planning method and system | |
Zhu et al. | Joint transportation and charging scheduling in public vehicle systems—A game theoretic approach | |
CN109117993B (en) | Processing method for optimizing vehicle path | |
Rong et al. | Mining efficient taxi operation strategies from large scale geo-location data | |
Yang et al. | Dynamic pricing for integrated energy-traffic systems from a cyber-physical-human perspective | |
Ren et al. | Improved ant colony optimization for the vehicle routing problem with split pickup and split delivery | |
CN111091329B (en) | Semi-open type vehicle path optimization method for multi-vehicle-type transportation of hazardous chemicals | |
Asghari et al. | Environmental and social implications of incorporating carpooling service on a customized bus system | |
CN113222275A (en) | Vehicle path optimization method considering space-time distance under time-varying road network | |
CN111667114A (en) | Intelligent route planning recommendation method based on space-time big data fusion | |
Meng et al. | Customer‐Oriented Vehicle Routing Problem with Environment Consideration: Two‐Phase Optimization Approach and Heuristic Solution | |
Wang et al. | Collaborative multicenter reverse logistics network design with dynamic customer demands | |
CN116415882A (en) | Instant distribution order distribution system for collaborative distribution of riders and unmanned vehicles | |
JP7291252B2 (en) | Processing route information | |
Hua et al. | Optimality-guaranteed algorithms on the dynamic shared-taxi problem | |
CN114418213A (en) | Urban electric vehicle scheduling method and system based on deep reinforcement learning | |
CN112001678B (en) | Method for reducing tobacco delivery mileage | |
CN112613701B (en) | Finished cigarette logistics scheduling method | |
CN113344267A (en) | Logistics network resource allocation optimization method based on cooperation | |
Lai et al. | Data-driven flexible vehicle scheduling and route optimization | |
CN116882883A (en) | Power grid material distribution path planning method utilizing big data | |
Lai et al. | Utility-based matching of vehicles and hybrid requests on rider demand responsive systems | |
Davatgari et al. | Electric vehicle supply equipment location and capacity allocation for fixed-route networks | |
Sağlam et al. | Prediction of parking space availability using arima and neural networks | |
Tang et al. | Forward and reverse courier logistics distribution path planning based on multi-objective optimization genetic algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230404 Address after: 351100 No. 185, Licheng South Avenue, Chengxiang District, Putian City, Fujian Province Patentee after: Fujian tobacco company Putian Co. Address before: 351131 No. 1459, Shaolin South Street, West Tianwei Town, Licheng District, Putian City, Fujian Province Patentee before: PUTIAN TOBACCO LOGISTICS Co.,Ltd. |