A novel hybrid algorithm for generalized traveling salesman problems in different environments

A swap sequence-based particle swarm optimization (SSPSO) technique and genetic algorithm (GA) are used in tandem to develop a hybrid algorithm to solve generalized traveling salesman problem. Local search algorithm K-Opt is occasionally used to move any stagnant solution. Here, SSPSO is used to find the sequence of groups of a solution in which a tour to be made and cities from different groups of the sequence are selected using GA. The K-Opt algorithm (for $$K=3$$ K=3 ) is used periodically for a predefined number of iterations to improve the quality of the solutions. The algorithm is capable of solving the problem in crisp as well as in imprecise environment. For this purpose, a general fitness evaluation rule for the solutions is proposed. The efficiency of the algorithm is tested in crisp environment using different size benchmark problems available in TSPLIB. In crisp environment, the algorithm gives 100% success rate for problems up to considerably large sizes. Imprecise problems are generated from crisp problems randomly using a rule and are solved using the proposed approach. The obtained results are discussed. Moreover it is observed that the proposed algorithm finds multiple optimal paths, when they exist, both for the crisp problems and their fuzzy variations.


Introduction
The traveling salesmen problem (TSP) is one of the standard combinatorial discrete optimization problems.The generalized traveling salesman problem (GTSP) is an extension of TSP and is a NP-hard problem.The GTSP has been introduced by [14], Saksena [31], and Srivastava [33] in the context of computer record balancing and of visit sequencing through welfare agencies since the 1960s.The problem consists of a set of n cities where travel cost c i j between (i-th city and j-th city) any two cities is known.The n cities are divided into several groups, i.e., g 1 , g 2 , . . ., g m .Here, m is the numbers of groups.A city may belong to one or more than one group.The objective is to find the possible minimum-cost Hamiltonian cycle (tour) through the groups.The GTSP has various real-world applications such as mail delivery [19], welfare agency routing ( [31], material flow system design [19], vehicle routing [19], and computer file sequencing [14].The existing algorithms for GTSP are mainly based on dynamic programming techniques [8,14,18,26,31,33], which transfer GTSP into TSP.Some heuristic and meta-heuristic algorithms have been proposed for solving GTSP: a random key genetic algorithm [34], a memetic algorithm [13], an efficient composite heuristic [30], reinforcing ant colony system [27], etc. Yang et al. [37] present a new heuristic method-based ant colony optimization (ACO) for solving GTSP.Since 1995, particle swram optimization (PSO) [16] has been proven to succeed in continuous optimization problems, and many works have been done effectively in this area.PSO can be used to solve GTSP also.Using the concept of the swap operator and swap sequence, and redefining some operators of PSO on the basis of them, Shi et al. [32] present another heuristic method based on PSO for solving TSP as well as GTSP.On the other hand, the genetic algorithm (GA) was introduced by Holland, inspired by Darwin's theory in 1970s [15].The idea behind GA is to model the natural evolution by using genetic inheritance together with Darwin's theory.Nowadays, GA is one of the important heuristic algorithms for NP-hard combinatorial optimization problems.Bontouxa et al. [3] propose a memetic algorithm using GA for solving GTSP.Zhao and Zhu [40] present an innovative GA-based algorithm for solving GTSP.
In all the above literature, it is implicitly assumed that the traveling cost from one city to another is fixed, i.e., crisp in nature.Traveling cost from one city to another depends on the conveyance used for traveling.It varies slightly depending on the availability of the conveyance, condition of the road, etc., though its value normally lies in an interval.So, these costs are imprecise in nature.There are different estimation approaches for imprecise data.Among them, stochastic estimation [22], fuzzy estimation [5,12,17] and rough estimation [24] have drawn more attention.For reallife problem like GTSP, fuzzy estimation of travel costs is more appropriate.It is less error prone, as these estimations are based on experts' opinion.Also, not much past data are required for such an estimation.Due to this reason, it is better to model the costs of a GTSP as fuzzy numbers.Though there are some research works on TSP incorporating fuzzy costs [4,5,17], none has solved GTSP in a fuzzy environment.
From the above discussion, it is clear that there are some lacunas in the existing literature of GTSP which are summarized below: • In the existing literature, GTSPs are considered only in a crisp environment.• Existing algorithms for GTSP are not capable of solving GTSPs in an imprecise environment.• The presence of multiple paths of any GTSP is overlooked/ignored by the researchers.
To remove the above-mentioned shortcomings in this research paper, a GTSP is considered whose costs are fuzzy in nature.To solve this problem, an algorithm is proposed which is capable of solving a GTSP in crisp environment as well as in fuzzy environment.The proposed algorithm is a hybridization of SSPSO, GA and K-Opt.Initially, a sequence of groups in which a tour is to be made is selected randomly and for each sequence a city from each group is selected using GA.After the first iteration, SSPSO is used to rearrange the sequence of groups for improvement of the solution and cities from different groups of the sequence selected using GA.In this GA, two types of crossover operations-single-point crossover and multipoint crossover-are performed depending on some random functions.For mutation operation, an adaptive mutation function is used that continues to perform mutation operation until no better movement is found in a predefined number of iterations.A strong local search algorithm K-Opt (K = 3) is used periodically for a predefined number of iterations to improve the quality of the solutions.This predefined number of iterations may vary with the size of the problem.The algorithm is capable of solving the problem in crisp as well as in imprecise environment.A general fitness evaluation scheme is proposed for the purpose.As a result, the algorithm can be used to solve GTSP in fuzzy, rough and stochastic environments.The efficiency of the algorithm is tested in crisp environment using different size benchmark problems available in TSPLIB [29].In crisp environment, the algorithm gives 100% success rate for problems up to considerably large sizes.Imprecise problems are generated from crisp problems randomly using a proposed rule.Imprecise problems are solved and the obtained results are discussed.
The rest of the paper is organized as follows: in Sect.2, problem definition and some mathematical prerequisites are presented.In Sect.3, some features of SSPSO are discussed.The features of GA are discussed in Sect. 4. The K-Opt algorithm is presented in Sect. 5.The proposed algorithm is presented in Sect.6.The experimental results are discussed in Sect.7. A brief conclusion is drawn in Sect.8.At the end, the reference list is presented.

Problem definition and mathematical prerequisites
The GTSP can be described as the problem of seeking a special Hamiltonian cycle with the lowest cycle cost in a complete weighted graph.Let G = {V, E, C} be a complete weighted graph where , . . ., n}} are a set of cities, set of edges between cities, and the set of travel cost between two cities, respectively.The symbols v i , e i j , c i j represent the i-th city, the edge connecting cities v i and v j , and the travel cost corresponding to edge e i j , respectively.The city set V is partitioned into m possibly intersecting groups g 1 , g 2 , . . ., g m with |g j | ≥ 1 and V = m j=1 g j , where |g j | is the number of elements in group |g j |.The special Hamiltonian cycle is required to pass through all of the groups, but not all of the cities differing from that of TSP.In real-life problem, c i j represents the cost of travel between city i and city j, or distance between city i and city j, or time required to travel from city i to city j, etc.There are two different kinds of GTSP under the above-mentioned framework of the special Hamiltonian cycle [8,19]: (1) the cycle passes through exactly one city in each group (Fig. 1) [14,31,33] and (2) the cycle passes through at least one city in each group (Fig. 2) [18,26].The first kind of GTSP is known as E-GTSP, where E stands for equality [8].In this paper only the E-GTSP, i.e., the first case, is discussed and is still called GTSP for convenience.The problem is to find a travel schedule for a salesman who likes to visit exactly one city of each group (g 1 , g 2 , . . ., g m ) and comes back to the starting city with minimum tour cost.Here, it is assumed that groups of cities are known, i.e., the salesman will not determine the group.
Fuzzy number A fuzzy number ã is a fuzzy set in (set of real numbers) characterized by a membership function, μ ã (x), which is both normal (i.e., ∃ at least one point x ∈ s.t.μ ã (x) = 1) and convex [39].

Triangular fuzzy number (TFN)
A TFN ã = (a 1 , a 2 , a 3 ) has three parameters a 1 , a 2 , a 3 , where a 1 < a 2 < a 3 , whose membership function μ ã (x) is given by Possibility (Pos), necessity (Nes) measure Let ã and b be two fuzzy numbers with membership functions μ ã and μ b, respectively.Then taking the degree of uncertainty as the semantics of fuzzy number: where is any one of the relations >, <, =, ≤, ≥.Analogously, if b is a crisp number, say b, then On the other hand, the necessity measure of an event ã b is a dual of the possibility measure.The grade of necessity of an event is the grade of impossibility of the opposite event and is defined as where ã b represents the complement of the event ã b.Credibility (Cr) measure Using (3) and ( 4), the credibility measure of an event ã b is denoted by Cr( ã b) and is defined as [21] Cr( ã b) = and is given by [17] Cr( Fuzzy GTSP It is stated earlier that in a real-life GTSP, c i j represents the cost of travel between city i and city j, or distance between city i and city j, or time to travel from city i to city j, etc.Though the distance between two cities is normally fixed, the time to travel from a city to another is always imprecise in nature in the fuzzy sense.The same property holds for traveling cost from a city to another.The distance of travel from a city to another also sometimes varies due to bad roadways or different/alternative ways.So in the general problem, it is better to consider c i j as fuzzy number ci j .In this research work, ci j is considered as a TFN (c i j1 , c i j2 , c i j3 ).GTSP with imprecise ci j is named fuzzy GTSP.
3 Swap sequence-based particle swarm optimization for GTSP PSOs are exhaustive search algorithms based on the emergent motion of a flock of birds searching for food [6,16] and has been extensively used/modified to solve complex decisionmaking problems in different fields of science and technology [1,7,9,11,12,28].A PSO normally starts with a set of potential solutions (called swarm) of the decision-making problem under consideration.Individual solutions are called particles and location of food is the optimal solution.In simple terms, the particles are flown through a multi-dimensional search space, where the position of each particle is adjusted according to its own experience and that of its neighbors.Each particle i has a position vector (X i (t)), a velocity vector (V i (t)), the position at which the best fitness (X pbesti (t)) has been encountered by the particle so far, and the best position of all particles (X gbest (t)) in the current generation t.In generation (t + 1), the position and velocity of the particle are changed to X i (t + 1) and V i (t + 1) using the following rules: The parameters c 1 and c 2 are set to constant values, which are normally taken as 2, r 1 and r 2 are two random values uniformly distributed over [0, 1], and w(0 < w < 1) is the inertia weight which controls the influence of previous velocity on the new velocity.It is mainly used to solve continuous optimization problems.It is also used to solve traveling salesman problems, where swap sequence and swap operations are used to find the velocity of a particle and its updating [20,35,38].A PSO that uses swap sequence and swap operation is called SSPSO.As discussed in §2, in a GTSP a potential solution is represented by a sequence of cities and one and only one city is selected from a cluster.SSPSO is used to find the sequence of groups from which the cities are to be selected for a solution.In SSPSO, swap operations on different groups are used to update the sequence of groups of a solution.A swap sequence represents a sequence of swap operations used to transform the sequence of groups from a solution to another solution.The basic operations of SSPSO are briefly presented below: Swap operator Consider a normal solution sequence of group in GTSP with n cities, and m groups X = (x 1 , x 2 , . . ., x m ), where x i ∈ {g 1 , g 2 , . . ., g m } and each g i are distinct.
Here, swap operator, SO(i, j) is defined as exchange of group x i and x j in solution sequence X .Then we define X = X + SO(i, j) as a new sequence on the operating operator SO(i, j) on X .So, the plus sign '+' above has its new meaning.We explain with a concrete exam-ple: suppose there is a GTSP problem with six groups, and X = (x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ) = (g 1 , g 3 , g 5 , g 2 , g 4 , g 6 ) be a sequence.Let the swap operator be SO (2,4), then , groups of position 2 and position 4 are exchanged.
Swap sequence A swap sequence SS is made up of one or more swap operators.Let SS = (SO 1 , SO 2 , . . ., SO p ), where SO 1 , SO 2 , . . ., SO p are swap operators.Swap sequence acting on a sequence of groups of a solution means all the swap operators of the swap sequence act on the sequence of groups in order.This can be described by the following formula: Different swap sequences acting on the same sequence of groups may produce the same new sequence.All these swap sequences are named an equivalent set of swap sequences.
In the equivalent set, the sequence which has the least swap operator is called basic swap sequence of the set or basic swap sequence (BSS) in short.
Several swap sequences can be merged into a new swap sequence.Here, the operator ⊕ is defined as merging two swap sequences into a new swap sequence.Suppose there are two swap sequences, SS1 and SS2, which act on one sequence X in order, namely SS1 first and SS2 second, a new sequence of groups X is obtained.Let there be another swap sequence SS acting on the same sequence of groups X and gets the sequence of groups X , then SS is called merging of SS1 and SS2 and is represented as: Here, SS and SS1 ⊕ SS2 are in the same equivalent set.
The construction of basic swap sequence Suppose there are two solutions, A and B, and our task is to construct a basic swap sequence SS which can act on B to get solution A.
The transformation of the particle updating formulas For solving the GTSP formula, (7) and (8) of PSO have to be transformed using swap sequences and swap operations as follows: Here, r 1 , r 2 are random numbers between 0 and 1. Velocity ) should be maintained with the probability of r 1 , i.e., each swap operator in BSS (X pbesti (t) − X i (t)) should be selected with probability r 1 .The same meaning is for the expression r 2 (X gbest (t) − X i (t)).From here, it is seen that the bigger the value of r 1 , the greater is the influence of X pbesti (t) for more swap operators.

Genetic algorithm for GTSP
GA are exhaustive search algorithms based on the mechanics of natural selection and genesis (crossover, mutation, etc.) and have been introduced by Holland [15] inspired by Darwin's theory in the 1970s.The idea behind GA is to model the natural evolution by using genetic inheritance together with Darwin's theory.Basic operations of GA are selection, crossover and mutation.In natural genesis, it is known that chromosomes are the main carriers of hereditary information from parent to offspring and that genes, which present hereditary factors, are lined up on chromosomes.At the time of reproduction, crossover and mutation take place among the chromosomes of parents.In this way, hereditary factors of parents are mixed up and carried to their offsprings.Again, Darwinian principle states that only the fittest animals can survive in nature.So, a pair of fittest parent normally reproduces a better offspring.The same phenomenon is followed to create a genetic algorithm for an optimization problem.A GA normally starts with a set of potential solutions (called initial population) of the decision-making problem under consideration.Individual solutions are called chromosomes and the fitness of each chromosome is evaluated by the evaluation process.Solutions are selected from the population randomly depending on their fitness by the selection process for the mating pool.Crossover and mutation operations happen among these selected solutions to get a new set of solutions and it continues until terminating conditions are encountered.
It has already been stated that to solve GTSP, PSO is used to determine the sequence of groups in which a sales-man should travel.GA is used to find appropriate cities from different groups in a sequence of groups determined by PSO, so that the total cost spent by the salesman is a minimum for that particular sequence.Consider a GTSP with n cities v 1 , v 2 , . . ., v n and m groups g 1 , g 2 , . . ., g m .Let X = (x 1 , x 2 , . . ., x m ), where x i ∈ {g 1 , g 2 , . . ., g m } be a sequence of groups for a solution which is determined by PSO in a particular iteration.The major steps for the selection of cities from different groups by GA for such a sequence are discussed below.
Initialization Initially, a set of N -solutions where y i j ∈ x j for j = 1, 2, . . ., m and are randomly selected from x j .
Evaluation The fitness of a solution Y i is taken as the ratio of the number of solutions of the population dominated by Y i and the population size N and is represented by f (Y i ).Let the total travel cost spent by a salesman for a solution Y i be C(Y i ) and for another solution Selection All solutions of the population do not take place in the evolution process (crossover and mutation).Mainly, solutions with higher fitness are selected for the process.Various types of selection techniques are available to select solutions from the initial population for the matting pool [10,23].Here, roulette wheel (RW) selection technique [23]is used to select N solutions from the population for the mating pool depending on their fitness.The following are the steps of this selection: (i) Find the total fitness of the population F = N i=1 f (Y i ).(ii) Calculate the probability of selection p i of each solution Y i by the formula p i = f (Y i ) F .(iii) Calculate the cumulative probability cp i for each solution Y i by the formula cp i = i j=1 p j .(iv) Generate a random number r in the interval (0, 1).

(vi) Repeat steps (iv) and (v) N times to select N solutions
for the mating pool.One solution may be selected more than once.CY 1 , CY 2 , created by two-point crossover process, on the parents are given below [10](c.f., Fig. 4): Mutation The mutation process takes place on some selected solutions from the mating pool.Selection takes place according to the probability of mutation.The steps for this purpose are presented below: (i) Selection for mutation For each solution of N pop , generate a random number r from the interval (0, 1).If r < pm, then the solution is taken for mutation, where pm is the probability of mutation.(ii) Mutation process To mutate a selected solution PY = ( py 1 , py 2 , . . ., py m ), select a random integer number k in the range [1, m].Then py k is replaced by a randomly selected city from its group to get a new mutated solution.This process is repeated for a finite number of iterations or until an improvement of the selected solution is made.

K-Opt operation for GTSP
K-Opt is a local search algorithm [2] mainly used for TSP which is based on the exchange of K parts (sub-tours) and their reverses (reverse sub-tours) of a tour (path) of the TSP under consideration to find a better tour.In the PSO part of the algorithm, it is tacitly used to find the optimal sequence of groups in which a salesman should select a particular city from a group in his tour to minimize the cost.If it is assumed that groups are connected by edges, then, while breaking (removes) K edges of a sequence of groups in a tour, there are (K − 1)!2 K −1 ways to reconnect it (including the initial sequence) to form a valid sequence.Each new combination gives a new sequence of groups.We find the corresponding optimal tour using GA.Among these sequences, one may produce a better tour than the tour produced by the original sequence and can be taken as an improvement.In the case of the 2-Opt algorithm, we remove two edges from the sequence and reconnect all combinations of sub-sequences and their reverses (Fig. 5).Continue this process until no 2-Opt improvements can be found.Similarly in the case of 3-Opt, breaking three edges in a sequence there are in total eight cases of reconnection (Fig. 6).If a sequence is 3-optimal, it is also 2-optimal [2].Continue break (remove) edges from the sequence, i.e., K = 1, 2, 3, . . ., n and get new algorithms, such as 2-Opt, 3-Opt, 4-Opt and so on.But increase of K increases the time complexity.Due to this, the 3-Opt operation is used and it is found that it acts better than the 2-Opt operation for large size GTSPs.In the proposed method, the 3-Opt operation is periodically used in a predefined number of iterations to improve the quality of the solutions in the PSO part of the algorithm and at the end of the algorithm.

Proposed algorithm for GTSP
A Hybrid algorithm based on SSPSO and GA algorithm with a local search (K-Opt) for GTSP Consider a GTSP with n cities {v 1 , v 2 , . . ., v n } and m groups {g 1 , g 2 , . . ., g m }.
The group g k contains n k cities, i.e., m k=1 n k = n.In the algorithm, the i-th solution is represented by a triplet {X i , Y i , C(Y i )}, where X i = (x i1 , x i2 , . . ., x im ), represents the sequence of groups in which travel is to be made, i.e., x i j ∈ {g 1 , g 2 , . . ., g m } and each x i j are distinct.Y i = (y i1 , y i2 , . . ., y im ) represents the actual path corresponding to X i , i.e., y i j ∈ x i j for j = 1, 2, . . ., m. C(Y i ) is the path cost of Y i .PSO is used to find the optimal sequence (X i ) of groups of a solution and GA is used to select cities from different groups to find the optimal path Y i for X i .Here, SS is the swarm size, MaxGen is the positive integer that represents the maximum number of iterations for PSO and t is the iteration counter.for k = 1 to SS do Determine Vk(t) using equation ( 9).Determine Xk(t) using equation (10).

GA algorithm for selection of cities from different groups to find the optimal path corresponding to a sequence of groups X k (t)
In the algorithm, c is the iteration counter, Maxgen is the maximum number of iterations, pc is the probability of crossover, and pm is the probability of mutation.
iii.Generate Ipop = {Y1, Y2, ..., YN }. iv.Evaluate path length and fitness of each solution.v. Repeat a. Select N solutions from Ipop, for mating pool using Roulette-wheel selection process (One solution may be selected more than once).Let this set be Npop.b.Select solutions from Npop, for crossover and mutation depending on pc and pm respectively.c.Make crossover on selected solutions for crossover.d.Make mutation on selected solutions for mutation.e. Set Ipop = Npop.f. c = c + 1. Until(c> maxgen).vi.Return the path with smallest tour cost.vii.End Algorithm

K-Opt operation on sequence of groups for improvement of sequence of groups X k (t)
The detailed algorithm of K-Opt operation for K = 3 is presented below.In the algorithm, a one-dimensional array of size m, X tem k (t), is used to represent a temporary sequence of group corresponding to X k (t) in iteration t.Y tem k (t) is used to represent the actual path corresponding to X tem k (t).X ki (t) and X r ki (t), i = 1, 2, 3 are one-dimensional arrays used to represent the sub-sequence of group and revers_sub-sequence of group of the original sequence of groups X k (t).Maxit3 is the maximum number of iterations.
for i = 1 to Maxit3 do • Remove three edges (randomly selected) from the sequence of group X k (t), it makes three sub-sequences of group X ki (t), i = 1, 2, 3. • Reverses of the contents of these sub-sequence of groups are called revers_sub-sequence, represented as X r ki (t), i = 1, 2, 3, i.e., X r k1 (t) = revers_sub-sequence of groupX k1 (t), X r k2 (t) = revers_sub-sequence of X k2 (t), X r k3 (t) = revers_sub-sequence of X k3 (t).• Now, combining the sub-sequences of groups {X k1 (t), X k2 (t), X k3 (t)}, {X r k1 (t), X r k2 (t), X r k3 (t)}, a new sequence of groups can be formed in the following eight combinations: Create a new sequence of group from the combination and let it be Xtem k (t).Find Ytem k (t) using GA algorithm ( §6.1). Determine

Experimental results
All computational experiments are conducted with Dev C++ 5.8.3, core i3 CPU @ 2.10 GH z , Windows 8.1 Operating System and 4 GB RAM.The performance of the proposed algorithm is tested using different size standard GTSPs from TSPLIB.Each problem algorithm is tested by running the program five times for different seeds of random number generator and the best solution is selected.The average cost of these solutions and percentage of relative error (Error(%)) according to the cost of optimal solution are calculated.The percentage of relative Error(%) is calculated using the following equation.

Error (%)
= average cost of solutions − cost of the optimal solution cost of the optimal solution × 100. ( The results obtained by the proposed algorithm for 21 different test problems from TSPLIB are presented in Table 1 (Fig. 7).
In Table 1, the 1st column stands for the problem name, 2nd (n) for the number of cities, 3rd (m) for the number of groups in the GTSP, 4th for the cost of optimal solution for each problem [29], 5th for the cost of the best solution obtained by the proposed method, 6th-10th for the obtained results in different runs of the algorithm using different seeds of random number generator for each problem, 11th and 12th columns provide some statistical information about the problem, such as average cost of the obtained solutions and percentage of relative error of the average cost of the solutions.The last column presents the average computational time in seconds.It is found from Table 1 that the solution produced by the algorithm for each of the considered problems is the same as its optimal solutions [29].As a result, the obtained average cost of solutions of each problem is the same as the cost of the optimal solution of the corresponding problem.The relative error of each problem is exactly zero, i.e., the proposed algorithm obtained exact solution for each problem considered for testing.As GTSP belongs to the NP-hard problem group, it is not possible to find the time complexity of the algorithm used for finding its solution.But a graph is drawn using different size problems and corresponding computational time required for finding the optimal solution (cf.Fig. 8).As expected, it is clear from Fig. 8 that computational time increases with problem size.It should be noted that computational time depends on several factors like processor speed, operating system, size of catch memory, different software running (in back-ground) in the machine, etc.
Again, solutions (paths and corresponding costs) obtained in different runs of the algorithm using different seeds of random number generator are analyzed and it is observed that for many problems there are multiple optimal paths.The different optimal paths for some problems are listed in Table 2.
Fuzzy GTSP For fuzzy GTSP, no standard test problems are available in the literature.So, here, crisp GTSP problems from TSPLIB are used to generate fuzzy GTSP problems.Only TFNs are used to represent costs of a fuzzy GTSP.Let c i j be the cost of travel from city i to city j for a crisp GTSP from TSPLIB.For the corresponding fuzzy GTSP, fuzzy cost of travel ci j from city i to city j is considered as ci j = (c i j1 , c i j2 , c i j3 ), where c i j2 = c i j , c i j1 = c i j − R1, c i j3 = d i j + R2, and R1 and R2 are randomly generated in the interval (0, R × c i j /100).For 5% fuzziness R is taken as 5, for 10% fuzziness R is taken as 10, etc. From each crisp GTSP, three different fuzzy GTSPs are generated considering R = 5, R = 10 and R = 15, respectively.After generating a problem, it is solved using the proposed algorithm for fuzzy GTSP and the results obtained are presented in Table 3.It is found from Table 3 that for small size problems and for small amount of fuzziness (5%), a mid value of optimum cost as well as optimal path is same as the corresponding crisp GTSP.While fuzziness increases, the optimum cost as well as optimum path varies gradually.All these observations agree with reality.The last column of Table 3 presents the computational time (s) for the problems with 5 % fuzziness.For fuzzy GTSP, also it is found that for some problems there are multiple optimal paths.For multiple paths, the mid values of the fuzzy objective (TFN) for all the optimal paths of a problem are the same.The other two components (left and right limits) differ by a very small amount (cf.Table 4).
Table 5 represents a comparison of some computational results of the proposed algorithm with other existing algorithms in the literature.In Table 5, the results of PSO [32], GCGA [36] and GA [34] are taken from the corresponding reference paper and the results of the proposed method are taken from a single run of the algorithm for the corresponding problem.From Table 5, it is clear that the proposed approach is better compared to some other existing approaches in the literature with respect to the accuracy (Error(%)).For the test problems like 14ST 70, 20K RO D100, 22P R107, 25P R124, and 40D198, other algorithms also produce optimal solution, but with respect to the statical information(error) the proposed algorithm is better compared to some other algorithms in the literature.On the other hand, GCGA [36] provide results in least computational time, but its success rate is not 100%.Compared with GA [34], the proposed algorithm does not always provide optimum results in lesser computational time.In some problems, our algorithm takes less computational time and in some other problems GA [34] takes less time.But it should be noted that computational time depends on several factors such as processor speed, operating system, size of catch memory, and different software running (in the background) in the machine.The PSO [32] and GCGA [36] produce optimal results in some problems, but they are not capable of producing optimal solution in all runs of the algorithms for different test problems (at least up to a moderate size).For this reason, in this paper PSO and GA are combined together to create a hybrid algorithm to improve the efficiency of the algorithm.Probability of crossover ( pc) and probability of mutation ( pm) are two important parameters of GA, and the performance of GA on a particular problem mostly depends on them.Though these parameters are problem dependent, three test problems (relatively large sizes) are used to check the performance (with respect to computation time) of GA with respect to these parameters and a comparative study table is presented in Table 6.The result is obtained for a particular seed of random number generator.From the table, it is observed that the algorithm takes minimum run time to produce optimal solution of the test problems for pc = 0.8 and pm = 0.6.Due to this reason, these values of pc and pm are taken for obtaining the results of other problems also.A graphical representation of the requirement of time to produce optimal solution due to different pc and pm for the two test problems-25PR124 and 29PR144-are presented in Fig. 7.
Table 7 represents the results obtained by the proposed method for different test problems using 2-Opt and 3-Opt operations in the algorithm.It is observed that, for small size problems like 11E I L51 and 14ST 70, both the approaches produce the same solution as the optimal one.For large size problems, the algorithm produces different solutions using 2-Opt and 3-Opt.Problems for which optimal solutions are obtained by the algorithms are presented in boldface in Table 7.It is clear from Table 7 that for all the problems, the algorithm with 3-Opt produces better result than that using 2-Opt.So in the proposed algorithm, 3-Opt operation is used.

Conclusion
For the first time combining the features of SSPSO technique and GA, a hybrid algorithm is developed to solve GTSP in different environments.SSPSO is used to find the sequence of groups in a solution and for each sequence cities from different groups are selected using GA.Here K-Opt algorithm (for K = 3) is used periodically in a predefined number of iterations in the PSO portion of the algorithm to improve the quality of the solutions.It is tacitly used to find the optimal sequence of groups in which a salesman should select a particular city from a group in his tour to minimize the tour cost.
Existing algorithms for GTSP are not suitable for GTSP in an imprecise environment.The proposed algorithm is capable of solving the problem in crisp as well as in imprecise environment.The efficiency of the algorithm is tested in crisp environment using different size benchmark problems available in TSPLIB [29].In crisp environment, the algorithm gives 100% success rate for problems up to considerably large problem sizes.As the results of fuzzy GTSP and crisp GTSP are the same for relatively small size problems with small amount of fuzziness, it can be concluded that the algorithm is efficient for solving fuzzy GTSP.The algorithm can be used to solve GTSPs in different other imprecise environments, like stochastic environment, rough environment, etc.The only difference is that a comparison criteria of rough objectives have to be developed using a valid comparison operator like trust measure [24].In stochastic environment, random objectives have to be be compared using a valid comparison technique like the chance constraints approach [25].
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

1 2 [ 5 ) 1
Pos( ã b) + Nes( ã b)].(Lemma If ã and b are two fuzzy numbers, then the credibility measure of the event ã < b is denoted by Cr( ã < b) Lemma1).It is a valid fuzzy comparison, since cr ( Ã < B) + cr ( Ã ≥ B) = 1.The comparison of objective values for the particles of PSO is made in the same manner.
(iv) Upgradation of parents If the path cost of CY 1 is smaller than that of PY 1 , then replace PY 1 with CY 1 .Similarly, if the path cost of CY 2 is smaller than that of PY 2 , then replace PY 2 with CY 2 .

Table 1
Computational results obtained by different runs of the algorithm due to considered test problems from TSPLIB

Table 2
Different optimal paths of some test problems from TSPLIB obtained by the algorithm

Table 3
Results obtained for different fuzzy TSPs generated from some problems of TSPLIB by the proposed algorithm

Table 4
Different paths obtained for different fuzzy TSPs generated from some problems of TSPLIB by the proposed algorithm

Table 5
Comparison of results obtained by the proposed algorithm with respect to some other existing algorithms in the literature

Table 6
Time required (in s) to find the optimal solution of different test problems from TSPLIB using different values of pm and pc