Next Article in Journal
Special Issue on Advanced Learning Methods for Complex Data
Previous Article in Journal
A Comparison of Word Embeddings and N-gram Models for DBpedia Type and Invalid Entity Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Traveling-Salesman-Problem Algorithm Based on Simulated Annealing and Gene-Expression Programming

1
Global Energy Interconnection Research Institute Co., Ltd., Beijing 102200, China
2
Power Systems Artificial Intelligence Joint Laboratory of SGCC, Beijing 102200, China
3
Institute of Advanced Technology, Nanjing University of Posts and Telecommunications, Nanjing 230001, China
4
Electric Power Research Institute of State Grid Shanghai Municipal Electric Power Company, Shanghai 200437, China
*
Author to whom correspondence should be addressed.
Information 2019, 10(1), 7; https://doi.org/10.3390/info10010007
Submission received: 22 November 2018 / Revised: 17 December 2018 / Accepted: 21 December 2018 / Published: 25 December 2018

Abstract

:
The traveling-salesman problem can be regarded as an NP-hard problem. To better solve the best solution, many heuristic algorithms, such as simulated annealing, ant-colony optimization, tabu search, and genetic algorithm, were used. However, these algorithms either are easy to fall into local optimization or have low or poor convergence performance. This paper proposes a new algorithm based on simulated annealing and gene-expression programming to better solve the problem. In the algorithm, we use simulated annealing to increase the diversity of the Gene Expression Programming (GEP) population and improve the ability of global search. The comparative experiments results, using six benchmark instances, show that the proposed algorithm outperforms other well-known heuristic algorithms in terms of the best solution, the worst solution, the running time of the algorithm, the rate of difference between the best solution and the known optimal solution, and the convergent speed of algorithms.

1. Introduction

With the development of Big Data technology, graphics computing is applied to all kinds of fields. For example, in power-inspection route planning [1,2], all task points and equipment can be regarded as nodes in the graph, and path length between each task point or equipment can be regarded as an edge; thereby, a road graph on power inspection is formed. In another example, an entire social network [3,4,5] is a graph in which the user is the vertex of the graph, and the relationship between the user and another user is regarded as the edge. As can be seen from the above examples, as a storage or presentation form in a specific application, how to mine the knowledge contained in a graph by quickly traversing the vertices in the graph is a research hotspot in the field of graphics computing.
Traversing all the vertices in the graph and repeating is an important application direction in graph-vertex traversal. Traversing all the vertices in the graph is a classical traveling-salesman problem (TSP) in the shortest time. It is well known that TSP is an NP-hard problem [6,7]. To solve the problem, besides traditional backtracking algorithms, branch-and-bound algorithms, and greedy algorithms, we mainly use heuristic search algorithms to optimize. These heuristic algorithms include the simulated annealing algorithm [8,9,10,11,12], tabu-search algorithm [13,14,15], ant-colony optimization algorithm [16,17,18], and genetic algorithm [19,20,21,22,23,24]. However, the optimization time of the simulated annealing algorithm is very long, the tabu-search algorithm has a strong dependence on the initial solution and can only be serialized, genetic algorithm is easy to fall into local optimization and has poor convergence performance, and the ant-colony algorithm also has long optimization time and is also prone to falling into local optimization [25].
Inspired by the characteristics of genetic expression, Candida Ferreira first proposed the concept of gene-expression programming (GEP) in 2001 [26,27,28]. It is a new member of the evolutionary-computing family based on genetic algorithms and genetic programming with a strong-function mining ability and high-function mining efficiency. Compared with traditional genetic algorithms and genetic programming, GEP has the following advantages [29,30,31]: (1) it uses simple encoding to solve a complex problem; (2) its convergent speed is improved by 2–4 times in comparison to genetic algorithms and genetic programming; and (3) it is easy for it to overcome combination exploration and premature phenomena. This paper proposes a traveling-salesman problem based on simulated annealing and GEP (TSP-SAGEP) in order to traverse all cities in the shortest time. The major contributions of our work are listed as follows:
(1)
To better improve the ability to global search, we applied a simulated annealing algorithm to the genetic operation of GEP, and propose an improved GEP algorithm based on simulated annealing (IGEP-SA).
(2)
In order to solve traveling-salesman problems like graph vertex traversal optimization, we present a traveling-salesman optimization algorithm on the basis of IGEP-SA.
(3)
Experimental results show that the proposed algorithm outperforms traditional algorithms in terms of the best solution, the worst solution, the running time of the algorithm, the rate of difference between the best solution and the known optimal solution, and the convergent speed of the algorithms.
The remainder of this paper is organized as follows. Section 2 discusses the related works. Section 3.1 introduces the improved GEP algorithm based on simulated annealing. Section 4 focuses on the traveling-salesman problem based on simulated annealing and GEP. Section 5 conducts detailed comparative experiments. Finally, conclusions are given in Section 6.

2. Related Work

Zhan et al. [8] proposed a new simulated annealing algorithm, called a list-based simulated annealing algorithm, in order to solve the traveling-salesman problem. Experimental results indicated that the proposed algorithm in this paper had competitive performance compared with the other algorithms. Yu et al. [9] presented a hybrid simulated annealing algorithm based on tabu search to solve the traveling-salesman problem. Experimental results demonstrated that the proposed algorithms improved accuracy and efficiency. Absalom et al. [10] presented a simulated annealing algorithm based on symbiotic-organism search in order to better solve the traveling-salesman problem. Comparative results showed that the proposed algorithm had advantages in terms of convergence, average execution time, and percentage deviations. Behnck et al. [11] applied a modified simulated annealing algorithm to SUAVs’ path, planning to find an optimal path that included most of points. Zhao et al. [12] used an improved simulated annealing algorithm to solve the traveling-salesman problem. Results indicated that the proposed algorithm in this paper could obtain the best solutions for most TSPLIB benchmarks. Xu et al. [13] proposed a new tabu-search algorithm based on the evolutionary and ant-colony algorithms to effectively solve the traveling-salesman problem. Archetti et al. [14] proposed an integer linear programming model based on the tabu-search algorithm to solve the vehicle-routing problem. Simulated experiments illustrated that the proposed algorithm had good performance compared with other algorithms. Chiang et al. [15] used a quantum-inspired tabu-search algorithm to optimize the traveling-salesman problem. Experimental results verified effectiveness of the proposed algorithm. Yang et al. [18] proposed an improved ant-colony optimization algorithm based on swarm intelligence to obtain the best solution for the traveling-salesman problem. Comparative results showed that the algorithm had higher accuracy and efficiency compared with classic optimization algorithms. He et al. [16] put forward an improved ant-colony algorithm based on new crossover and mutation operations to find the best solution for the traveling-salesman problem. Simulation results demonstrated that the proposed algorithm had advantages in terms of stability and optimization capacity. Mohsen et al. [17] proposed a hybrid ant-colony optimization algorithm based on the simulated annealing algorithm to solve the traveling-salesman problem. El-Samak et al. [19] applied the affinity propagation clustering technique to optimize the genetic algorithm for finding the best solutions to the traveling-salesman problem. Rani et al. [20] proposed the improved genetic algorithm by a roulette-wheel selection operator with different crossovers and mutation rates to solve the traveling-salesman problem. Experimental results demonstrated that the algorithm in this paper is better than the existing crossover operator in other algorithms. Deng et al. [21] presented the improved genetic algorithm by using k-means to generate a new initial-population strategy for solving the traveling-salesman problem. Wang et al. [22] solved the traveling-salesman problem by merging two optimization strategies into the traditional genetic algorithm. Comparative results showed that the proposed algorithm in this paper could find better solutions than a traditional optimal algorithm. Wang et al. [23] solved the traveling-salesman problem by means of a multioffspring genetic algorithm according to biological evolutionary and mathematical ecological theory. Compared with the traditional genetic algorithm, the proposed algorithm in this paper was faster and decreased the evolutionary time of the best solution. Changdar et al. [24] introduced a multiobjective genetic algorithm into solving a multiobjective solid traveling-salesman problem. In order to compare performance of the genetic algorithm and the simulated annealing ant-colony optimization, Haroun et al. [25] performed a comparative experiment to evaluate the performance of the three algorithms in terms of computational time and the shortest distance in solving the traveling-salesman problem.
GEP [26] is also a kind of heuristic algorithm. Many researchers have already applied GEP to combinatorial optimization. Ferreira et al. [32] applied a new chromosomal organization based on multigenes in GEP to solve the combinatorial optimization problem. Results demonstrated that the algorithm had high efficiency and accuracy. Sabar et al. [33] used gene-expression programming to generate the acceptance criterion for solving an combinational optimization problem. Comparative results showed that the solution was better than the other results obtained from the existing algorithms.

3. Improved GEP Algorithm Based on Simulated Annealing

Although the GEP algorithm can quickly find the solution for combinatorial optimization problems, such as the traveling-salesman problem, GEP is still easy to fall into the local optimum. In order to better solve the problem, this paper proposes an improved GEP algorithm based on simulated annealing (IGEP-SA) to avoid GEP falling into a local optimum.

3.1. Simulated Annealing

Simulated annealing (SA) is an effective and general form of optimization and is based on simulating the annealing of solids [34,35]. It is very useful in finding global optima in the presence of local-optima large numbers. The simulated-annealing algorithm starts from a higher temperature, which is called the initial temperature. When the temperature gradually decreases, the solution of the algorithm tends to be stable. However, the solution may be a local optimal solution. Then in SA, such a local optimal solution is jumped with a certain probability to find the global optimal solution of the objective function. The flow chart of the simulated-annealing algorithm is shown in Figure 1.
From Figure 1, we can see that the Metropolis rule is very important for SA to find the optimal solution. Suppose that state x o l d becomes a state x new when the system is subject to some disturbance. Then, in the Metropolis rule, the energy of the system also changes from E ( x o l d ) to E ( x new ) , and the acceptance probability of the system changing from state x o l d to state x new is p :
p = 1 if E ( x n e w ) < E ( x old ) exp ( E ( x new ) E ( x old ) T ) if E ( x n e w ) E ( x old )
The steps of the simulated-annealing algorithm are as follows (Algorithm 1):
Algorithm 1: SA
Input:    initial temperature T 0 , minimum temperature T m i n , maximum number of iteration K m ax ,
probability of temperature drop ρ ;
Output: the optimal solution x best ;
1. Generating an initial solution x 0 ;
2.  x best x 0 ;
3. Computing the value of the objective function f ( x 0 ) and f ( x b e s t ) ;
4.  T i T 0 ;
5. while T i > T m i n do
6.      Δ f f ( x new ) f ( x best ) ;
7.     if Δ f < 0 then
8.          x best x new ;
9.     end if
10.    if Δ f 0 then
11.         p e Δ f T ;
12.        if c random [ 0 , 1 ] p then
13.              x best x new ;
14.        else
15.              x best x best ;
16.        end if
17.    end if
18.     i i + 1 ;
19.     T i ρ × T i ;
20. end while
21. Return x best ;

3.2. IGEP-SA

Because the simulated-annealing algorithm has the advantage of jumping out of the local optimum, gene-expression programming was improved by using the simulated-annealing algorithm in this paper. Firstly, this paper used simulated annealing to dynamically generate a GEP population. Then, we applied the simulated-annealing algorithm to the mutation operator of GEP to improve the ability of global search.

3.2.1. New Population Generation Based on Simulated Annealing

In order to improve the diversity of the GEP population, this paper proposes new population-generation algorithm based on SA (NPG-SA). The steps of NPG-SA are shown as follows (Algorithm 2):
Algorithm 2: NPG-SA
Input:     current temperature T, number of population P o p S i z e , fitness value of previous
generation population f ( o l d ) , and fitness value of population after genetic operations f ( n e w ) ;
Output: new population N e w P o p ;
1.  i 1 ;
2. while i < PopSize do
3.        Δ f f ( new ) [ i ] f ( old ) [ i ] ;
4.        p e Δ f T ;
5.       if c random [ 0 , 1 ] p then
6.            NewPop . add ( i ) ;
7.       end if
8.        i i + 1 ;
9. end while
10. Return N e w P o p ;

3.2.2. New Mutation Operator Based on Simulated Annealing

We know that it is very important for GEP to have a diverse population. A mutation operator is the most efficient way. By mutation operator, the best individual of the population in GEP can be found. However, the existing mutation operator in GEP is implemented according to a certain probability. The way depends on the initialization of mutation probability so that the mutation operator in GEP is strongly subjective. This paper proposes a new mutation operator based on SA (NMO-SA). The steps of NMO-SA are shown in Algorithm 3.
Algorithm 3: NMO-SA
Input:     current temperature T, number of population P o p S i z e , fitness value of the previous
generation population f ( o l d ) , and fitness value of population after genetic operations f ( n e w ) ;
Output: new population N e w P o p ;
1.  i 1 ;
2.  j SelectPo i n t ( o l d P o p ) ;
3. while i < PopSize do
4.        Δ f f ( new ) [ i ] f ( old ) [ i ] ;
5.        q e Δ f T ;
6.        p 0 . 001 ; / / initial mutation rate
7.       if c random [ 0 , 1 ] q then
8.            N e w P o p Mutation ( o l d P o p , j , c ) ;
9.       else
10.          N e w P o p Mutation ( o l d P o p , j , p ) ;
11.     end if;
12.       j SelectPo i n t ( N e w P o p ) ;
13.       i i + 1 ;
14. end while
15. Return N e w P o p ;

3.2.3. Description of IGEP-SA

To better improve the evolutional efficiency of traditional GEP and avoid falling into the local optimum, we applied Algorithms 2 and 3 to IGEP-SA. The steps of IGEP-SA are shown in Algorithm 4.
Algorithm 4: IGEP-SA
Input:     current temperature T, number of population P o p S i z e , fitness value of the previous
generation population f ( o l d ) , fitness value of population after genetic operations f ( n e w ) , probability
of temperature drop ρ , M a x G e n , P t , and P r ;
Output: new population N e w P o p ;
1.  P o p I n i t P o p ( P o p S i z e ) ;
2.  f ( o l d ) E v a l u a t e f i t n e s s ( P o p ) ;
3. while i < MaxGen do
4.        f ( n e w ) E v a l u a t e f i t n e s s ( P o p ) ;
5.        P o p NPG SA ( T , P o p S i z e , f ( o l d ) , f ( n e w ) ) ;
6.        P o p NMO SA ( T , P o p S i z e , f ( o l d ) , f ( n e w ) ) ;
7.        P o p I S T r a n s p o s i t i o n ( P t , P o p , P o p S i z e ) ;
8.        P o p R I S T r a n s p o s i t i o n ( P t , P o p , P o p S i z e ) ;
9.        P o p G e n e T r a n s p o s i t i o n ( P t , P o p , P o p S i z e ) ;
10.       P o p O n e P o i n t R e c o m b i n a t i o n ( P r , P o p , P o p S i z e ) ;
11.       P o p T w o P o i n t R e c o m b i n a t i o n ( P r , P o p , P o p S i z e ) ;
12.       P o p G e n e R e c o m b i n a t i o n ( P r , P o p , P o p S i z e ) ;
13. end while
14. N e w P o p P o p ;
15. R e t u r n N e w P o p ;

4. Traveling-Salesman Problem Based on Simulated Annealing and GEP

The traveling-salesman problem is a kind of classical combinatorial optimization problem. To simplify the solution of the problem, in this paper, we suppose that the traveling-salesman problem is symmetric. We know that the traveling-salesman problem is an NP-hard problem. For the TSP with N cities, there are ( N 1 ) ! solutions. With the increase of N, the number of solutions for TSP grow exponentially. We propose a traveling-salesman problem based on simulated annealing and GEP (TSP-SAGEP) in order to better solve the problem.

4.1. TSP-SAGEP Code

The traveling-salesman problem is a kind of combinational optimization problem. For a combinational optimization problem, a chromosome of GEP is a multigene chromosome composed of one-element genes. However, in a multigene chromosome, one-element genes are very important because they can be organized in multigene families (MGFs). These MGFs are composed of clusters of related-gene encoding. An example is shown as follows:
Example 1.
The different cities in the traveling-salesperson problem can be encoded in an MGFs where each gene code represent a particular city visited by the salesman. Consider the simple chromosome below, composed of one MGF with nine cities:
The spatial structure of the chromosome, which is shown in Figure 2, is shown in Figure 3. This is the traveling route of the traveling-salesman problem, where B is both the starting and finishing city and shown in gray.

4.2. Fitness Function of the TSP-SAGEP

The fitness function is very important for GEP to solve the traveling-salesman problem. We know that the objective of the traveling-salesman problem needs to find the shortest route on the condition that we start and lastly finish the city. The fitness function helps algorithm move in the correct direction.
It is obvious that we cannot directly use tour length as a measure of fitness, as the shorter the tour is, the fitter the individual [26]. So, for each generation, fitness f i of the ith individual in the gth generation is evaluated by the equation:
f i = T g t i + 1
where T g is the length of the largest tour is encoded in the chromosomes of the current population, and t i is the length of the tour encoded in the ith individual. From Equation (1), we know that the fitness of the worst individual of the population is always equal to one.

4.3. Algorithm Description

In order to better solve graph vertex traversal optimization problems, this paper proposes a traveling-salesman problem based on TSP-SAGEP. The description of TSP-SAGEP is similar to Algorithm 4, but the differences between TSP-SAGEP and Algorithm 4 are as follows:
(1)
In TSP-SAGEP, the population is initialized according to Figure 1.
(2)
In TSP-SAGEP, the fitness value of the individual in the population is evaluated according to Equation (1).
(3)
TSP-SAGEP returns the best graph vertex traversal line.

5. Experiment and Analysis

5.1. Experimental Environment

To better explain the effectiveness and feasibility of the proposed algorithms, we conducted related experiments in a laboratory environment. All experiments were implemented in Java on an Intel Core-i7 PC and Windows 10. All experiments were done by using 6 TSP standard benchmark problems, with different lengths, from TSPLIB [36,37]. Experimental datasets are shown in Table 1. These datasets can be divided into three groups according to their dimension. The first group includes four instances that vary in dimension, between 26 and 70 cities. The second group includes four instances that vary in dimension, between 120 and 442 cities. The last group includes four instances that vary in dimension, between 666 and 1304 cities. Finally, after conducting the experiments 50 times for each instance, the best experimental results were obtained.
As a kind of heuristic algorithm, the performance of TSP-SAGEP is affected by the different values of the parameters. Table 2 provides the parameters that the proposed algorithm takes.
The length of the chromosome in TSP-SAGEP is especially set according to the dimensions in Table 1.

5.2. Experimental Analysis

Experiment 1: To verify the effectiveness of IGEP-SA, we compare dTSP-SAGEP with traditional GEP for all datasets in Table 1. Four measures that include the best solution, the worst solution, number of generations, and the average running time of two algorithms were used to evaluate two algorithms. Table 3 shows that TSP-SAGEP algorithm is better than the traditional GEP algorithm in the above four measures. Meanwhile, we observed the relationship between the number of runs and the tour length for all datasets in Table 1. Results are shown in Figure 4.
Table 3 shows that the TSP-SAGEP algorithm is superior over a traditional GEP in computational results for all datasets in Table 1. The best solution based on TSP-SAGEP is better than the solutions based on a traditional GEP for all test instances. In Experiment 1, we used the rate of difference to evaluate the advantages and disadvantages of experiment results based on TSP-SAGEP and traditional GEP. The rate of difference is defined in Equation (2).
ρ = P b e s t P o p t i m a l P b e s t
where ρ represents the rate of difference, P b e s t represents the best solution based on TSP-SAGEP, and P o p t i m a l represents the known optimal solution. Rate of difference ρ based on TSP-SAGEP and traditional GEP is 0, 0, 0, 0.00636, 0.033, and 0.0287 for fri26, st70 and gr120 instances. For pcb442, gr666 and rl1304 instances, respectively, and the rate of difference based on TSP-SAGEP and traditional GEP is 0.00065, 0.000207, 0.000616, 0.00615, 0.0535, and 0.1008, respectively. These results illustrate that the best solution based on TSP-SAGEP is better than traditional GEP. This is mainly because new population generation policy and mutation operator based on simulated annealing can improve the diversity of the GEP population and increase the probability of searching for a global solution. Figure 4 demonstrates that, when the number of runs increases from 1 to 10, tour length based on TSP-SAGEP is closer to the optimal tour length than traditional GEP for six instances. From Figure 4, we can also find that, in all 10 runs, the number of solving the best solution based on TSP-SAGEP is 7, 7, 7, 6, 5 and 5 for the fri26, st70, gr120, pcb442, gr666, and rl1304 datasets, respectively.
Experiment 2: In order to better demonstrate the advantage of the proposed algorithm in this paper, we compared TSP-SAGEP with Ezugwu2017 [10], Mohsen2016 [17], and Wang 2016 [23] in terms of the best solution, the worst solution, and rate of difference, which is obtained by the Equation (2). Results are shown in Table 4. Figure 5 shows the comparison of the number of obtaining the best solution based on five algorithms for six instances in all 10 runs. Figure 6 shows comparison of convergent speed between TSP-SAGEP and three other algorithms for six instances. In Experiment 2, we also used convergent speed to evaluate the performance of four algorithms. Convergent speed is defined as follows.
C s = P G T G
where C s represents the convergent speed of the algorithm, P G represents the current generation corresponding to the best solution, and T G represents the total number of generations. From Equation (3), we know that the smaller the C s value of the algorithm is, the faster the convergence speed of the algorithm.
Table 4 shows that TSP-SAGEP algorithm is superior over Ezugwu2017, Mohsen2016, and Wang 2016 in computational results for all datasets in Table 1. The best solution based on TSP-SAGEP is better than the solutions based on other three algorithms for all test instances. According to Equation (2), we can see that the rate of difference based on TSP-SAGEP is the same as that of the other three algorithms for the fri26, st70 and gr120 instances, while the rate of difference based on TSP-SAGEP is obviously less than the other three algorithms for the pcb442, gr666 and rl1304 instances. Experimental results illustrate that, in comparison to the other three algorithms, the best solution obtained by TSP-SAGEP is the closest to the optimal solution.
According to Figure 5, we know that the number of the best solutions obtained by TSP-SAGEP is obviously more than other algorithms for the pcb442, gr666, and rl1304 instances, and the number of the best solutions obtained by TSP-SAGEP is equal to the number of the best solutions obtained by other algorithms for the fri26 and st70 instances. Generally, for the pcb442, gr666, and rl1304 instances, the number of the best solutions obtained by TSP-SAGEP is improved by 28.57%, 42.86%, and 57.14%, respectively, in comparison to Ezugwu2017 and Wang2016; compared with Mohsen2016, the number of the best solutions obtained by TSP-SAGEP is improved by 28.57%, 42.86%, and 57.14%, respectively. This illustrates that TSP-SAGEP easily solves the best solution on the condition that the number of runs is the same. Figure 6 indicates that convergent speed of TSP-SAGEP is superior over the other three algorithms for six instances. Specifically, compared with Ezugwu2017, Mohsen2016, and Wang2016, the convergent speed of TSP-SAGEP is improved by 75%, 62.5%, and 57.14% for the fri26, st70, and gr120 instances, respectively; for pcb442, the convergent speed of TSP-SAGEP is improved by 62.5%, 58.9%, and 50%, in comparison to Ezugwu2017, Mohsen2016, and Wang2016, respectively. However, when the number of cities increases to 666 and 1304, the convergent speed of four algorithms is almost the same. This is mainly because that the increase in the number of cities cause the increase of solution scale so that it is difficult for the above four algorithms to solve the best solution in the 1000 generations.

6. Conclusions

In this paper, we introduce a traveling salesman problem based on TSP-SAGEP in order to traverse all vertices of a graph in the shortest time. To verify the performance of the proposed algorithm, two experiments were conducted for six benchmark instances obtained from the TSPLIB, and the experimental results of TSP-SAGEP were compared with three other heuristic algorithms. Results show that TSP-SAGEP outperformed the other heuristic algorithms proposed in [10,17,23,26]. Specifically, compared with traditional GEP, Ezugwu2017, Mohsen2016, and Wang2016, TSP-SAGEP had the best solution quality, low error between the best solution, a known optimal solution for all instances, and high convergent speed. In the future, we will enhance the proposed algorithm in this paper by using other GEP operations in order to better solve the problem of large-scale graph vertex traverse. We will also introduce the proposed algorithm into the asymmetric application of graph vertex traverse.

Author Contributions

A.-H.Z. and S.D. contributed to the conception of the study; A.-H.Z., L.-P.Z., B.H. and Y.S. designed and performed the experiments; H.-B.Q. and S.P. analyzed the data; and A.-H.Z. wrote the manuscript.

Funding

This research was funded by State Grid Company Research Project grant number 5455HJ160005 and the National Natural Science Foundation of China grant number 51507084, 61363037.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, Y.; Li, Z.-W. Application of Ant Colony Algorithm to Power Cable Patrol Route Planning. Comput. Syst. Appl. 2015, 5, 135–139. [Google Scholar]
  2. Wang, J.; Wang, X.; Chen, Y.; Cai, G. Intelligent patrolling robot addressing program based on graph method. Autom. Electr. Power Syst. 2007, 9, 78–81. [Google Scholar]
  3. Robins, G.; Pattison, P.; Kalish, Y.; Lusher, D. An introduction to exponential random graph (p*) models for social networks. Soc. Netw. 2007, 29, 173–191. [Google Scholar] [CrossRef]
  4. Soussi, R.; Aufaure, M.-A.; Baazaoui, H. Towards social network extraction using a graph database. In Proceedings of the 2010 Second International Conference on Advances in Databases Knowledge and Data Applications (DBKDA), Menuires, France, 2–6 June 2010; pp. 28–34. [Google Scholar]
  5. Cattuto, C.; Quaggiotto, M.; Panisson, A.; Averbuch, A. Time-varying social networks in a graph database: A Neo4j use case. In Proceedings of the First International Workshop on Graph Data Management Experiences and Systems, New York, NY, USA, 22–27 June 2013; pp. 1–6. [Google Scholar]
  6. Arora, S. Polynomial time approximation schemes for Euclidean traveling salesman and other geometric problems. J. ACM 1998, 45, 753–782. [Google Scholar] [CrossRef] [Green Version]
  7. Cook, W.; Espinoza, D.; Goycoolea, M. Computing with Domino-Parity Inequalities for the TSP; Georgia Institute of Technology, Industrial and Systems Engineering: Atlanta, GA, USA, 2007; Volume 19, pp. 356–365. [Google Scholar]
  8. Zhan, S.; Lin, J.; Zhang, Z.; Zhong, Y. List-based simulated annealing algorithm for traveling salesman problem. Comput. Intell. Neurosci. 2016, 2016, 8–18. [Google Scholar] [CrossRef] [PubMed]
  9. Lin, Y.; Bian, Z.; Liu, X. Developing a dynamic neighborhood structure for an adaptive hybrid simulated annealing–tabu search algorithm to solve the symmetrical traveling salesman problem. Appl. Soft Comput. 2016, 49, 937–952. [Google Scholar] [CrossRef]
  10. Ezugwu, A.E.S.; Adewumi, A.O.; Frîncu, M.E. Simulated annealing based symbiotic organisms search optimization algorithm for traveling salesman problem. Expert Syst. Appl. 2017, 77, 189–210. [Google Scholar] [CrossRef]
  11. Behnck, L.P.; Doering, D.; Pereira, C.E.; Rettberg, A. A modified simulated annealing algorithm for SUAVs path planning. IFAC-PapersOnLine 2015, 48, 63–68. [Google Scholar] [CrossRef]
  12. Zhao, D.; Xiong, W.; Shu, Z. Simulated annealing with a hybrid local search for solving the traveling salesman problem. J. Comput. Theor. Nanosci. 2015, 12, 1165–1169. [Google Scholar] [CrossRef]
  13. Xu, D.; Weise, T.; Wu, Y.; Lassig, J.; Chiong, R. An investigation of hybrid tabu search for the traveling salesman problem. Bio-Inspired Comput. Theor. Appl. 2015, 562, 523–537. [Google Scholar]
  14. Archetti, C.; Guastaroba, G.; Speranza, M.G. An ILP-refined tabu search for the directed profitable rural postman problem. Discret. Appl. Math. 2014, 163, 3–16. [Google Scholar] [CrossRef]
  15. Chiang, H.-P.; Chou, Y.-H.; Chiu, C.-H.; Kuo, S.-Y.; Huang, Y.-M. A quantum-inspired Tabu search algorithm for solving combinatorial optimization problems. Soft Comput. 2014, 18, 1771–1781. [Google Scholar] [CrossRef]
  16. Min, H.; Dazhi, P.; Song, Y. An improved hybrid ant colony algorithm and its application in solving TSP. In Proceedings of the 2014 IEEE 7th Joint International Information Technology and Artificial Intelligence Conference (ITAIC), Chongqing, China, 20–21 December 2014; pp. 1534–1551. [Google Scholar]
  17. Mohsen, A.M. Annealing ant colony optimization with mutation operator for solving tsp. Comput. Intell. Neurosci. 2016, 2016, 1–13. [Google Scholar] [CrossRef] [PubMed]
  18. Yang, J.; Ding, R.; Zhang, Y.; Cong, M.; Wang, F.; Tang, G. An improved ant colony optimization (I-ACO) method for the quasi travelling salesman problem (Quasi-TSP). Int. J. Geogr. Inf. Sci. 2015, 29, 1534–1551. [Google Scholar] [CrossRef]
  19. El-Samak, A.F.; Ashour, W. Optimization of traveling salesman problem using affinity propagation clustering and genetic algorithm. J. Artif. Intell. Soft Comput. Res. 2015, 5, 239–245. [Google Scholar] [CrossRef]
  20. Rani, K.; Kumar, V. Solving travelling salesman problem using genetic algorithm based on heuristic crossover and mutation operator. Int. J. Res. Eng. Technol. 2014, 3, 27–34. [Google Scholar]
  21. Deng, Y.; Liu, Y.; Zhou, D. An improved genetic algorithm with initial population strategy for symmetric TSP. Math. Probl. Eng. 2015, 2015, 1–6. [Google Scholar] [CrossRef]
  22. Wang, Y. The hybrid genetic algorithm with two local optimization strategies for traveling salesman problem. Comput. Ind. Eng. 2014, 70, 124–133. [Google Scholar] [CrossRef]
  23. Wang, J.; Ersoy, O.K.; He, M.; Wang, F. Multi-offspring genetic algorithm and its application to the traveling salesman problem. App. Soft Comput. 2016, 43, 415–423. [Google Scholar] [CrossRef]
  24. Changdar, C.; Mahapatra, G.S.; Pal, R.K. An efficient genetic algorithm for multi-objective solid travelling salesman problem under fuzziness. Swarm Evol. Comput. 2014, 15, 27–37. [Google Scholar] [CrossRef]
  25. Haroun, S.A.; Jamal, B.; Hicham, E.H. A Performance Comparison of GA and ACO Applied to TSP. Int. J. Comput. Appl. 2015, 117, 28–35. [Google Scholar] [CrossRef] [Green Version]
  26. Ferreira, C. Gene Expression Programming in problem Solving; Springer: Berlin, Germany, 2002; pp. 635–653. [Google Scholar]
  27. Ferreira, C. Automatically Defined Functions in Gene Expression Programming; Springer: Berlin, Germany, 2006; pp. 21–56. [Google Scholar]
  28. Ferreira, C. Genetic representation and genetic neutrality in gene expression programming. Adv. Complex Syst. 2002, 5, 389–408. [Google Scholar] [CrossRef]
  29. Ferreira, C. Gene Expression Programming: Mathematical Modeling by an Artificial Intelligence; Springer: Berlin, Germany, 2006. [Google Scholar]
  30. Wu, J.; Li, T.Y.; Jiang, Y.; Li, Z.L.; Liu, Y.Y. Gene Expression Programming Based on Diversified Development Strategy. J. Jilin Univ. Inf. Sci. Ed. 2010, 28, 396–403. [Google Scholar]
  31. Tang, C.; Zhang, T.; Zuo, J.; Rui, W.; Jia, X. Knowledge discovery based on Gene Expression Programming-History, achievements and future directions. Comput. Appl. 2004, 10, 7–10. [Google Scholar]
  32. Ferreira, C. Combinatorial Optimization by Gene Expression Programming: Inversion Revisited; Springer: Berlin, Germany, 2002; pp. 160–174. [Google Scholar]
  33. Sabar, N.R.; Ayob, M.; Kendall, G.; Qu, R. A dynamic multiarmed bandit-gene expression programming hyper-heuristic for combinatorial optimization problems. IEEE Trans. Cybern. 2015, 45, 217–228. [Google Scholar] [CrossRef] [PubMed]
  34. Van Laarhoven, P.J.M.; Aarts, E.H.L. Simulated Annealing; Springer: Berlin, Germany, 1987; pp. 7–15. [Google Scholar]
  35. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by simulated annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed]
  36. Reinelt, G. TSPLIB-A traveling salesman problem library. ORSA J. Comput. 1991, 3, 376–384. [Google Scholar] [CrossRef]
  37. Discrete and Combinatorial Optimization—IWR, Heidelberg. Available online: http://comopt.ifi.uni-heidelberg.de/software/TSPLIB95/tsp/ (accessed on 21 November 2018).
Figure 1. A flowchart of the simulated-annealing algorithm.
Figure 1. A flowchart of the simulated-annealing algorithm.
Information 10 00007 g001
Figure 2. Chromosome of a multigene family.
Figure 2. Chromosome of a multigene family.
Information 10 00007 g002
Figure 3. Expression of the above chromosome.
Figure 3. Expression of the above chromosome.
Information 10 00007 g003
Figure 4. Comparison of tour length based on TSP-SAGEP and GEP for six datasets with an increase in the number of runs.
Figure 4. Comparison of tour length based on TSP-SAGEP and GEP for six datasets with an increase in the number of runs.
Information 10 00007 g004
Figure 5. Comparison of the number of the best solutions obtained by four algorithms for six instances in all ten runs.
Figure 5. Comparison of the number of the best solutions obtained by four algorithms for six instances in all ten runs.
Information 10 00007 g005
Figure 6. Comparison of convergent speed among four algorithms for six instances.
Figure 6. Comparison of convergent speed among four algorithms for six instances.
Information 10 00007 g006
Table 1. Datasets used in experiments.
Table 1. Datasets used in experiments.
InstanceCitiesOptimal Solution
fri2626937
st7070675
gr1201206942
pcb44244250,778
gr666666294,358
rl13041304252,948
Table 2. Parameters of traveling-salesman problem based on simulated annealing and gene-expression programming (TSP-SAGEP).
Table 2. Parameters of traveling-salesman problem based on simulated annealing and gene-expression programming (TSP-SAGEP).
Parameter ItemValue
Initial temperature100
Minimum temperature 10 8
The probability of temperature drop0.98
Population size500
The number of runs10
The number of generations1000
IS Transposition rate0.3
RIS Transposition rate0.3
Gene Transposition rate0.1
One-point recombination rate0.3
Two-point recombination rate0.3
Gene recombination rate0.1
Table 3. Comparison of four measures between TSP-SAGEP and a traditional GEP algorithm.
Table 3. Comparison of four measures between TSP-SAGEP and a traditional GEP algorithm.
InstancesOptimalAlgorithmBestWorstAverageGenerationsTime (s)
fri26937TSP-SAGEP9379719411507.0671
  GEP94311329618008.1629
st70675TSP-SAGEP67569967715012.1135
  GEP69882271080014.0931
gr1206942TSP-SAGEP69427406699520019.4612
  GEP714710,035762490022.3573
pcb44250,778TSP-SAGEP50,81152,14750,87830051.0964
  GEP51,09254,11351,67690057.1428
gr666294,358TSP-SAGEP294,419305,036295,54299079.7853
  GEP310,982410,368328,45999085.1204
rl1304252,948TSP-SAGEP253,104271,582254,699990160.2458
  GEP281,296334,871290,611990165.5572
Table 4. Comparison of three measures between four algorithms.
Table 4. Comparison of three measures between four algorithms.
InstancesOptimalAlgorithmBestWorstAverageRateTime (s)
fri26937TSP-SAGEP93797194107.0671
  Ezugwu2017937120495108.5247
  Mohsen2016937129995508.5916
  Wang2016937115194808.3103
st70675TSP-SAGEP675699677012.1135
  Ezugwu2017675846701014.9045
  Mohsen2016675907715015.1037
  Wang2016675831695014.0985
gr1206942TSP-SAGEP694274066995019.4612
  Ezugwu2017694211,2377786025.3581
  Mohsen2016694211,3547801027.0433
  Wang2016694211,0087655023.6709
pcb44250,778TSP-SAGEP50,81152,14750,8780.0006551.0964
  Ezugwu201751,10755,72351,9580.0064461.0876
  Mohsen201651,14356,09852,0440.0071463.9011
  Wang201651,09755,00651,8280.0062459.4571
gr666294,358TSP-SAGEP294,419305,036295,5420.0002179.7853
  Ezugwu2017311,855417,702331,0240.0561190.0014
  Mohsen2016312,096420,809332,5280.0568492.0745
  Wang2016311,003411,984329,1990.0535288.0163
rl1304252,948TSP-SAGEP253,104271,58225,46990.00062160.2458
  Ezugwu2017288,013340,015297,8130.12175169.0031
  Mohsen2016290,057341,298299,3050.12798170.1165
  Wang2016286,799335,990294,6370.11803167.9084

Share and Cite

MDPI and ACS Style

Zhou, A.-H.; Zhu, L.-P.; Hu, B.; Deng, S.; Song, Y.; Qiu, H.; Pan, S. Traveling-Salesman-Problem Algorithm Based on Simulated Annealing and Gene-Expression Programming. Information 2019, 10, 7. https://doi.org/10.3390/info10010007

AMA Style

Zhou A-H, Zhu L-P, Hu B, Deng S, Song Y, Qiu H, Pan S. Traveling-Salesman-Problem Algorithm Based on Simulated Annealing and Gene-Expression Programming. Information. 2019; 10(1):7. https://doi.org/10.3390/info10010007

Chicago/Turabian Style

Zhou, Ai-Hua, Li-Peng Zhu, Bin Hu, Song Deng, Yan Song, Hongbin Qiu, and Sen Pan. 2019. "Traveling-Salesman-Problem Algorithm Based on Simulated Annealing and Gene-Expression Programming" Information 10, no. 1: 7. https://doi.org/10.3390/info10010007

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop