Journal of Industrial Engineering

Volume 2013 (2013), Article ID 347825, 5 pages

http://dx.doi.org/10.1155/2013/347825

Research Article

A Novel Metaheuristic for Travelling Salesman Problem

Vahid Zharfi and Abolfazl Mirzazadeh

Industrial Engineering Department, Faculty of Engineering, Kharazmi University, Tehran 31979-37551, Iran

Received 30 April 2013; Accepted 3 July 2013

Academic Editor: Alan Chan

Copyright © 2013 Vahid Zharfi and Abolfazl Mirzazadeh. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

One of the well-known combinatorial optimization problems is travelling salesman problem (TSP). This problem is in the fields of logistics, transportation, and distribution. TSP is among the NP-hard problems, and many different metaheuristics are used to solve this problem in an acceptable time especially when the number of cities is high. In this paper, a new meta-heuristic is proposed to solve TSP which is based on new insight into network routing problems.

1. Introduction

Although there are various classic methods for solving optimization problems, but they are not always able to solve real and applied optimization problems. It is generally believed that these problems are known as difficult and complicated problems. TSP is a typical example of a very hard combinatorial optimization problem. TSP can be modelled as an undirected weighted graph, such that cities are the graph’s vertices, paths are the graph’s edges, and a path’s distance is the edge’s length. It is a minimization problem starting and finishing at a specified vertex after having visited each other’s vertex exactly once. Metaheuristics can solve high-dimensional problems quickly. Heuristics may be classified into two families: specific heuristics and metaheuristics. Specific heuristics are tailored and designed to solve a specific problem and/or instance. Metaheuristics are general-purpose algorithms that can be applied to solve almost any optimization problem. They may be viewed as upper level general methodologies that can be used as a guiding strategy in designing underlying heuristics to solve specific optimization problems [1]. The word metaheuristic was first used by Glover [2] during introducing tabu search concepts.

Metaheuristics gained more popularity in few recent decades, and various algorithms are proposed to solve the problem such as Genetic Algorithm [3–5], Ant Colony [6–10], Particle Swarm [11], and Simulated Annealing [12, 13].

In this paper, a novel metaheuristic is proposed for solving TSP. This algorithm is used to foresee the best next city so as to gain shorter rout. The 2-Opt and 3-Opt algorithms are also used in local search. They are probably the most basic local search heuristics for the TSP.

The paper is organized as follows: Section 2 illustrates the background theory of suggested algorithm. In Section 3, formulating algorithms are presented. Algorithm improvement by using k-Opt is discussed in Section 4. The algorithm is employed into several standard TSP instances in the next section and the results are reported. According to the obtained results, it is noticeable that the algorithm achieves proper answers, quickly so that the answer of the first iteration will be a good and acceptable answer.

Finally, Section 6 concludes the paper.

2. Background Theory of the Algorithm

In this model, it is supposed that every node in the network is a vibration source and has a special weight which vibrates the edges, so that the nodes with the heavier weight release more vibration. This emitted force by a node distributes in edges that are connected to. Shared amount of each node is dependent on the edges distances. Shorter edges have more quota of this force. In other words, if the length of an edge connected to the node is shorter, less vibration falls down along the edge and it is vibrated more intensely. There is a simple network (Figure 1) to illustrate the theory. The searcher is located at the node in the following figure and wants to move to the next node after .

Figure 1: Routing in a network.

Signs which lead to are amount of vibrations that each edge indicates. Therefore, vibration of the edges can be considered as tendency that each neighbour node has in order to be selected by the searcher.

Node as a sample has relationship with nodes , and . The distance for each edge is determined, and the sum of the inverted distances equals . Weight of this node is . Therefore, the vibration that this node sends to the node equals:

Also, about nodes and , the amount of vibration which is sent to node is, respectively,

If all weights are considered equal to 1, the searcher receives a total vibration with an amount of in node . As a result, the searcher at the node selects with 14% possibility , b with 34% possibility , and with 52% possibility .

Movement from node in addition to direct route can be transferred to node of other routes such as “" and “.” If we also consider binary routes, calculations would be changed. Thus, by considering routes with maximum two edges, there are three different routes from to a (Table 1). Also, about nodes and , similar calculations will be performed (Tables 2 and 3).

Table 1: Emitted force by node “c” in the network.

Origin
Destination
Route
Route distance
Inversion of route distance
Final edge of the route
Amount of vibration in the final edge

c
a
c →b→a
3
0.33
b →a

c
a
c →a
5
0.20
c →a

c
a
c →d→a
6
0.16
d →a

c
b
c →a→b
7
0.14
a →b

c
b
c →b
1
1.00
c →b

c
d
c →a→d
8
0.12
a →d

c
d
c →d
3
0.33
c →d

Total

2.28

Table 2: Emitted force by node “b” in the network.

Origin
Destination
Route
Route distance
Inversion of route distance
Final edge of the route
Vibration force in the final edge

b
a
b →a
2
0.50
b →a

b
a
b →c→a
6
0.16
c →a

b
c
b →a→c
7
0.14
a →c

b
c
b →c
1
1.00
b →c

b
d
b →a→d
5
0.20
a →d

b
d
b →c→d
4
0.25
c →d

Total

2.25

Table 3: Emitted force by node “d” in the network.

Origin
Destination
Route
Route distance
Inversion of route distance
Final edge of the route
Vibration force in the final edge

d
a
d →c →a
8
0.12
c →a

d
a
d →a
3
0.33
d →a

d
b
d →a →b
5
0.20
a →b

d
b
d →c →b
4
0.25
c →b

d
c
d →a →c
8
0.12
a →c

d
c
d →c
3
0.33
d →c

Total

1.35

According to the tables, it is possible to obtain an amount of vibration in each edge. For example, edge “” receives a force with an amount of unit of , unit of , and a force with an amount of unit of node . Therefore, “” edge vibrates as follows:

It is supposed that all weights are 1. Therefore, the amount of vibration in “” is 0.25. Similarly, “” and “” are 0.37 and 0.31, respectively. It should be mentioned that the total amount of these three values is not necessarily equal to 1. In the investigated example final total amount of vibration force of the three edges which end with is 0.93.

Now, the searcher is faced with vibrant edges in node Therefore, it selects “” with 40% possibility and moves toward “,” selects “” with 27% possibility and moves toward “,” or also selects “” with 33% possibility and moves toward node “.” Total amount of these possibilities equals 1.

3. Model Formulation

In the previous section we got acquainted with the method of routing by the algorithm. As mentioned, edges with more vibration are selected with more possibility.

The vibration of a supposed edge “” is derived from two factors:

(i)direct force which is emitted from jth node;(ii)indirect forces that are emitted from other nodes and distribute up to “.” In other words, it is released from node to node after passing node :

In this relation is the distance of the route initiated from that ends in after passing . Through definition is sum of inverted distances of routes connected to which has up to two edges (see (5)); it is possible to gain the force which is applied from node to edge by:

so that represents weight of kth node. is the applied force for “” resulting from . Vibration of “” equals

The final result of the previous calculations is matrix in which represents vibration degree in the edge “.” It is necessary to use another element in order to increase ability of this matrix. To do so we can use vision power matrix.

In distance matrix, represents the length of an edge which is connector of nodes and . The vision matrix is defined as follows:

By combination of two matrixes and , it is possible to use both visual sense and mechanical sensor of the searcher for detection and recognition of the best route.

The result of performed calculations in this section is possibility matrix :

The parameters and are investigated from the obtained results. A proper amount of these parameters is dependent on nature and structure of the investigated problem. Adjustment of and has an important role in algorithm efficiency. Investigation of the obtained results shows that and are appropriate for most of the instances.

4. Algorithm Improvement by Local Search

The suggested algorithm searches to get to the best route in the network globally. But it makes less chance to envelop all answers. Consequently, local search is necessary in order to cover the solution space more entirely. A local search algorithm starts from a candidate solution and then iteratively moves to a neighbour solution.

In optimization, 2-Opt and 3-Opt are simple local search algorithms for solving the travelling salesman problem [14, 15]. In this paper, we propose to combine the 2-Opt, 3-Opt, and the suggested algorithm to improve the answers. The algorithm steps are as follows:

Initialize matrix

For (i < Iteration Max) do:

    While the searcher complete a solution

            Choose probabilistically the next state to move;

            Add that move to the tabu list;

    end

    Implementation of the 2-Opt and 3-Opt;

    If (local best solution better than global solution)

        Save local best solution as global solution;

    end

    Update matrix for each edge that the searcher chose

    (Increase elements of matrix related to the selected

    route up to percent);

end

5. Computational Results

In this section we are going to solve different standard examples by using the proposed algorithm. A large number of standard samples are available in public references. These references include famous samples for combinatorial optimization problems. One of the valid references in travelling salesman problems is TSPLIB [16].

5.1. Quality of Solutions

In order to estimate the algorithm function from absolute difference obtained by optimized answer, the following relation is used:

so that is the obtained answer and is the optimized answer.

5.2. Robustness

Metaheuristics must be able to act properly in an extended range of samples or problems with the same parameters. It is possible to adjust metaheuristic parameters for a collection of samples in a good way which has less efficiency for other samples.

Different problems of travelling salesman have been solved with different sizes and structures by the algorithm. For all instances, 9, 2, and 0.2% are allocated for , and , respectively (Table 4).

Table 4: Searching for optimized solution of different travelling salesman problems.

Instance
Best known
First sol.
Best sol.
Average sol.
Worst sol.
Error

ftv33
1286
1462
1286
1286
1286
0.0000000

ftv38
1530
1591
1530
1535
1546
0.0032679

swiss42
1273
1364
1273
1273
1273
0.0000000

ftv44
1613
1742
1615
1619
1625
0.0037197

ry48p
14422
15341
14466
14493
14507
0.0049230

berlin52
7542
7951
7542
7542
7542
0.0000000

ftv64
1839
2022
1909
1925
1940
0.0467645

eil76
538
559
538
538
538
0.0000000

eil101
629
717
630
638
650
0.0143084

The optimal path length of the algorithm’s iterations for berlin52 is shown in Figure 2.

Figure 2: The optimal path length of proposed algorithm for berlin52.

According to the obtained results, it is noticeable that the algorithm achieves proper answers quickly so that the answer of the first iteration will be a good and acceptable answer. Accordingly, although the suggested algorithm shows its success alone in this issue, it can be used as a powerful initial solution generator for other solution-based metaheuristics in order to use advantages of other algorithms and getting more close to the optimized answer.

6. Conclusion

In this paper we presented a new metaheuristic that has new insight into network problems. At the end of every iteration, the algorithm also tries to improve both the performance of the algorithm and the quality of solutions for the Tsp by using local search algorithms in the name of 2-Opt and 3-Opt. As a result, the proposed algorithm can be used for variety network problems due to its insight to these kinds of problems. Therefore, through extending this behaviour we can move toward solving different optimization problems by this algorithm in a way that each position of problem solutions represents a node and each algorithm step is passing from one state to another state. Future work will be conducted to improve the proposed algorithm and examine other local search algorithms to enhance the performance of the method for improving solutions. Additional improvements might lie on the combination of our approach with the other metaheuristics like genetic algorithm, tabu search, simulated annealing, and so on.

References

T. El-ghazali, Meta-Heuristic from Design to Implementation, Wiley, 2009.

F. Glover, “Future paths for integer programming and links to artificial intelligence,” Computers and Operations Research, vol. 13, no. 5, pp. 533–549, 1986.

S. Yuan, B. Skinner, S. Huang, and D. Liu, “A new crossover approach for solving the multiple travelling salesmen problem using genetic algorithms,” European Journal of Operational Research, vol. 228, no. 1, pp. 72–82, 2013.

Y. Nagata and D. Soler, “A new genetic algorithm for the asymmetric traveling salesman problem,” Expert Systems with Applications, vol. 39, no. 10, pp. 8947–8953, 2012.

W. Hui, “Comparison of several intelligent algorithms for solving TSP problem in industrial engineering,” Systems Engineering Procedia, vol. 4, pp. 226–2235, 2012.

S. M. Chen and C. Y. Chien, “Solving the traveling salesman problem based on the genetic simulated annealing ant colony system with particle swarm optimization techniques,” Expert Systems with Applications, vol. 38, no. 12, pp. 14439–14450, 2011.

J. Yang, X. Shi, M. Marchese, and Y. Liang, “Ant colony optimization method for generalized TSP problem,” Progress in Natural Science, vol. 18, no. 11, pp. 1417–1422, 2008.

M. Dorigo and L. M. Gambardella, “Ant colonies for the travelling salesman problem,” BioSystems, vol. 43, no. 2, pp. 73–81, 1997.

J. Bai, G. K. Yang, Y. W. Chen, L. S. Hu, and C. C. Pan, “A model induced max-min ant colony optimization for asymmetric traveling salesman problem,” Applied Soft Computing, vol. 13, no. 3, pp. 1365–1375, 2013.

A. Colorni, M. Dorigo, F. Maffioli, V. Maniezzo, G. Righini, and M. Trubian, “Heuristics from nature for hard combinatorial optimization problems,” International Transactions in Operational Research, vol. 3, no. 1, pp. 1–21, 1996.

X. H. Shi, Y. C. Liang, H. P. Lee, C. Lu, and Q. X. Wang, “Particle swarm optimization-based algorithms for TSP and generalized TSP,” Information Processing Letters, vol. 103, no. 5, pp. 169–176, 2007.

K. Meer, “Simulated annealing versus metropolis for a TSP instance,” Information Processing Letters, vol. 104, no. 6, pp. 216–219, 2007.

X. Geng, Z. Chen, W. Yang, D. Shi, and K. Zhao, “Solving the traveling salesman problem based on an adaptive simulated annealing algorithm with greedy search,” Applied Soft Computing Journal, vol. 11, no. 4, pp. 3680–3689, 2011.

C. Engels and B. Manthey, “Average-case approximation ratio of the 2-opt algorithm for the TSP,” Operations Research Letters, vol. 37, no. 2, pp. 83–84, 2009.

K. H. Hsieh, “Fourier descriptors for 2-Opt and 3-Opt heuristics for traveling salesman problem,” Journal of the Chinese Institute of Industrial Engineers, vol. 28, no. 3, pp. 237–246, 2011.

http://elib.zib.de/pub/Packages/mp-testdata/tsp/tsplib/tsplib.html