1 Motivations
Many problems from combinatorial optimization, and in particular from scheduling, present a high degree of neutrality. Such a property means that a lot of different solutions have the same fitness value. This is a critical situation for local search techniques (although this is not the only one), since it becomes difficult to find a way to reach optimal solutions. However, up to now, this neutrality property has been underexploited to design efficient search methods. Barnett [1]
proposes a heuristic (the
Netcrawler process), adapted to neutral landscapes, that consists of a random neutral walk with a mutation mode adapted to local neutrality. The persequence mutation rate is optimized to jump from one neutral network to another. Stewart [2] proposes an Extrema Selection for evolutionary optimization in order to find good solutions in a neutral search space. The selection aims at accelerating the evolution during the search process once most solutions from the population have reached the same level of performance. To each solution is assigned an endogenous performance during the selection step to explore the search space area with the same performance more largely, and to reach solutions with better fitness values. Verel et al. [3] propose a new approach, the Scuba search, which has been tested on MaxSAT problems and on NKlandscapes with neutrality. This local search exploits the evolvability (ability of random variations to reach improving solutions) of neutral networks. On neutral plateaus, the search is guided by the evolvability of solutions. When there is no neighboring solution with a higher evolvability on the plateau, a fitness “jump” is performed.Knowing that the permutation flowshop scheduling problem (PFSP) of minimizing the makespan is hard in the general case [4], a large number of metaheuristics have been proposed so far for its resolution. Recently, 25 methods have been tested and their performance has been compared by Ruiz and Maroto [5]. The Iterated Local Search of Stützle [6] has been shown to be one of the most efficient algorithm to solve Taillard’s FSP instances [7]. More recently, Ruiz and Stützle [8] proposed a simple variant of Stützle’s ILS based on greedy mechanisms. The perturbation is made of a destruction and of a construction phase: jobs are removed from the solution and then reinserted in order to get a new configuration that yields the best possible fitness value. Given that the PFSP is known to have a high neutrality with respect to the insertion neighborhood operator [9], we can assume that a lot of moves which do not change the fitness value are allowed with such a perturbation strategy. Furthermore, the acceptance criterion used in [8] has the particularity to accept equivalent, and then neutral, solutions. We argue that this neutrality property could be used more explicitly, and in a more simple way.
This paper will not bring to a new heuristic which produces some new bestknown solutions for a given scheduling problem. The goal of this preliminary work is to give a minimal, and yet efficient approach based on Iterated Local Search (ILS), which exploits the neutrality property of the search space in a novel way. The approach explicitly balances the search between the exploitation of the plateaus in the landscape, and the exploration of new parts of the search space. Two main questions are addressed in this paper: () What are the performances of this neutralitybased approach on solving a difficult scheduling problem where neutrality arises? () What are the costs and the benefits of such an exploitation? A Neutralitybased Iterated Local Search (NILS) that performs neutral walks along the search is proposed. The performances and the dynamics of NILS are deeply analyzed on the PFSP.
The paper is organized as follows. Section 2 is dedicated to the flowshop scheduling problem, to the required definitions of neutral fitness landscapes, and to the main principles of iterated local search. Section 3 presents the Neutralitybased Iterated Local Search (NILS) proposed in the paper. Section 4 is devoted to the analysis of the NILS efficiency to solve the PFSP. Finally, the last section concludes the paper and gives suggestions for further research.
2 Background
2.1 The Permutation Flowshop Scheduling Problem
The Flowshop Scheduling Problem is one of the most investigated scheduling problem from the literature. The problem consists in scheduling jobs on machines . Machines are critical resources, i.e. two jobs cannot be assigned to the same machine at the same time. A job is composed of tasks , where is the task of , requiring machine . A processing time is associated with each task . We here focus on a Permutation Flowshop Scheduling Problem (PFSP), where the operating sequences of the jobs are identical for every machine. As a consequence, a feasible solution can be represented by a permutation of size (determining the position of each job within the ordered sequence), and the size of the search space is then .
In this study, we consider the minimization of the makespan, i.e. the total completion time, as the objective function. Let be the completion date of task , the makespan can be defined as follows:
2.2 Neighborhood and Local Search
The design of local search algorithms requires a proper definition of the neighborhood structure for the problem under consideration. A neighborhood structure is a mapping function that assigns a set of solutions to any feasible solution . is called the neighborhood of , and a solution is called a neighbor of . A neighbor of solution results of the application of a move operator performing a small perturbation to solution . This neighborhood is one of a key issue for the local search efficiency. A solution is a local optimum iff no neighbor has a better fitness value: , .
For the PFSP, we consider the insertion operator. This operator is known to be one of the bestperforming neighborhood structure for the PFSP [6, 10]. It can be defined as follows. Let us consider an arbitrary solution, represented here by a permutation of size (the number of jobs). A job located at position is inserted at position . The jobs located between positions and are shifted. The number of neighbors per solution is then [6], where stands for the size of the permutation.
2.3 Neutral Fitness Landscape
A fitness landscape [11] can be defined by a triplet , where is a set of admissible solutions (i.e. the search space), , a neighborhood structure, and is a fitness function that can be pictured as the height of the corresponding solutions.
A neutral neighbor is a neighboring solution having the same fitness value, and the set of neutral neighbors of a solution is then . The neutral degree of a given solution is the number of neutral solutions in its neighborhood. A fitness landscape is said to be neutral if there are many solutions with a high neutral degree. A neutral fitness landscape can be pictured by a landscape with many plateaus. A plateau is a set of pairwise neighboring solutions with the same fitness values: . A portal in a plateau is a solution that has at least one neighbor with a better fitness value, i.e. a lower fitness value in a minimization context. Notice that a local optimum can still have some neutral neighbors. A plateau which contains a local optimum is called a Local Optima Plateau (LOP).
The average or the distribution of neutral degrees over the landscape may be used to test the level of neutrality of a problem. This measure plays an important role in the dynamics of local search algorithms [12, 13]. When the fitness landscape is neutral, the main features of the landscape can be described by its LOP, that may be sampled by neutral walks. A neutral walk from to is a sequence of solutions belonging to where and and for all , . Hence, to escape from a LOP, a heuristic method has to perform a neutral walk in order to find a portal.
2.4 Neutrality in the PFSP
In scheduling problems, and in particular in flowshop scheduling problems, it is well known that several solutions may have the same fitness value, i.e. the neutral degree of feasible solutions is not null. This has been confirmed in a recent work that analyzes the neutrality of Taillard’s PFSP instances [9]. It has been shown that the neutrality increases when the number of jobs and/or the number of machines increase. Moreover, this study revealed that very few local optima have no neutral neighbor and that the local optima plateaus are numerous and quite large. Experiments highlighted that most random neutral walks on a LOP are able to find a portal, i.e. numerous portals exist. This indicates that exploring a plateau allows to lead to a more interesting solution. These comments lead us to make proposals on the way to exploit this neutrality in order to guide the search more efficiently. More precisely, we will study how this neutrality can be exploited within an iterated local search algorithm.
2.5 Iterated Local Search
Iterated Local Search (ILS) [14] is a simple and powerful heuristic methodology that applies a local move iteratively to a single solution. In order to design an ILS, four elements have to be defined: () an initialization strategy, () a (basic) local search, () a perturbation, and () an acceptance criterion. The initialization strategy generates one solution from the search space. The local search must lead the current solution to a local optimum. The perturbation process allows to modify a solution in a different way that the neighborhood relation, in order to jump over the landscape. The acceptance criterion determines whether the local optimum should be kept for further treatments. Let us remark that both the perturbation process and the acceptance criterion could be based on an history of solutions found during the search. Thus, during the search, if the current optimum is not satisfying, the search process is able to restart from a better local optimum, found in a previous iteration. The search is stopped once a termination condition is satisfied.
The general scheme of ILS, given in [14] is as follows: () after generating an initial solution, the local search finds a local optimum, until the termination condition is met, () the current solution is perturbed, () the local search is applied to find another local optimum, () it is accepted or not regarding the acceptation criterion. The best solution found along the search is returned as output of the algorithm.
3 Neutralitybased Iterated Local Search
In this section, we propose an ILS based on neutrality, called Neutralitybased Iterated Local Search (NILS). We first give the local search and the perturbation strategy used at every NILS iteration, and then we discuss the general principles of NILS.
3.1 Local Search: FirstImproving HillClimbing
NILS is based on the iteration of a FirstImproving HillClimbing (FIHC) algorithm, where the current solution is replaced by the first encountered neighbor that strictly improves it. Each neighbor is explored only once, and the neighborhood is evaluated in a random order. FIHC stops on a local optimum.
3.2 Perturbation: Neutral Walkbased Perturbation
In general, there exists two possible ways of escaping from a local optimum in a neutral fitness landscape: either performing neutral moves on the Local Optima Plateau (LOP) until finding a portal, or performing a kick move which is a ‘large step’ move. When a neutral move is applied, it is assumed that the exploitation of the neutral properties helps to find a better solution. On the contrary, when a kick move is applied, it supposes that portals are rare, and that the exploration of another part of the search space is more promising.
The proposed Neutral Walkbased Perturbation (NWP), given in Algorithm 1, deals with this tradeoff between exploitation and exploration of LOP. First, NWP performs a random neutral walk on a LOP. The maximum number of allowed neutral steps on the LOP is tuned by a parameter denoted by (Maximal Number of Steps). Along the neutral walk, as soon as a better solution (a portal) is found, the neutral walk is stopped and the current solution is replaced by this neighboring solution. Otherwise, if the neutral walk does not find any portal, the solution is kicked. The kick move corresponds to a large modification of the solution. Like in FIHC, each neighbor is explored only once, and the neighborhood is evaluated in a random order.
3.3 NILS: a Neutralitybased ILS
NILS is based on the FIHC local search and the NWP perturbation scheme. Its acceptance criterion always accepts the current solution. After the initialization of the solution and the first execution of the FIHC, NILS iterates two phases: () a neutral phase which performs neutral moves on LOP, () a strictly improving phase. The neutrality of the problem is taken into account during phase (). Indeed, the neutral walk on the LOP is used to cross a large part of the search space. When the density of portals on plateaus is low, the parameter avoids to spread on the plateau unnecessarily. In such a case, a “restart” technique is used: the solution is kicked to escape from its neighborhood. Phase () strictly improves the current solution using the neighborhood operator under consideration. The neutral phase allows to exploit the neutrality of the search space from the local optima by visiting the corresponding LOP. Thus, a LOP is considered as a large frontier until a better local optimum. When a portal is found, the FIHC algorithm is executed in order to find this local optimum. However, when the frontier is too large, such a local optimum can be difficult to reach. Thus, the NILS algorithm restarts its search from another part of the search space, expecting that the next LOP would be easier to cross. As a consequence, the tradeoff between exploitation and exploration of the plateaus is directly tuned by the single NILS parameter: the value.
4 Neutralitybased Iterated Local Search for the Permutation Flowshop Scheduling Problem
4.1 Experimental Design
Experiments are driven using a set of benchmark instances originally proposed by Taillard [7] for the flowshop scheduling problem, and widely used in the literature [6, 10]. We investigate different values for the number of machines and for the number of jobs . The processing time of job and machine
is an integer value generated randomly, according to a uniform distribution
. For each problem size (), 10 instances are available. Note that, as mentioned on the Taillard’s website^{1}^{1}1http://mistic.heigvd.ch/taillard/problemes.dir/ordonnancement.dir/ordonnancement.html, very few instances with machines have been solved to optimality. For  and machine instances, optimal solutions have been found, requiring for some of them a very long computational time. Hence, the number of machines seems to be very determinant in the problem difficulty. That is the reason why the results of the paper will be exposed separately for each number of machines.The performances of the NILS algorithm are experimented on the first Taillard’s instance of each size, given by the number of jobs () and the number of machines (). For each size, several values are investigated:

For ,

For ,

For ,

For ,

For ,
NILS is tested times for each value and each problem size. Let us remark that when , NILS corresponds to an ILS algorithm without any neutrality consideration. Hence, we are able to compare an ILS that uses neutrality to a classical one. As the distribution of the
final fitness values is not symmetrical, the average and the standard deviation are less statistically meaningful. Therefore, for each problem size and each
value, the median and the boxplot of the 30 fitness values of solutions found are presented. For the other statistical quantities, the average value and the standard deviation of the 30 executions are computed.In order to let any configuration enough time to converge, the NILS termination condition is given in terms of a maximal number of evaluations, set to . The insertion operator (Section 2.2) defines the neighborhood used in the first improving hillclimbing (FIHC) and in the random neutral walk (NWP). In the NILS perturbation scheme, the solution is kicked when the neutral walk does not find any portal. In this study, the kick move corresponds to 3 randomly chosen exchange moves. As empirically shown in [10], this kick is instancedependent. However, as our work does not attempt to set the best parameters for the kick move, we choose a proper value that obtained reasonably good performances for all instances.
4.2 Experimental Results

Literature  NILS  NILS  NILS  

best  median  
5  20  1278  1278  1278  1278 
50  2724  2724  2724  2724  
100  5493  5493  5493  5493  
10  20  1582  1582  1582  1582 
50  2991  3025  3042  3025  
100  5770  5770  5784  5771  
200  10862  10872  10885  10872  
20  20  2297  2297  2297  2297 
50  3847  3894  3925  3917  
100  6202  6271  6375.5  6315.5  
200  11181  11303  11405.5  11334.5  
500  26059  26235  26401.5  26335.5 
This section examines the benefit of taking neutrality into account in the search process. We denote by NILS, the NILS algorithm with the following parameter setting: . Results obtained with and are compared with each other. Table 1 presents, for each problem size ( ), the fitness value of the best solution found in comparison to the best known solution from the literature. Let us remind that a minimization problem is under consideration. The median of the 30 fitness values found by NILS and NILS with the larger value (NILS) are also given. For machine instances, the best solution is reached by every NILS configuration. This means that using the neutral walk is as interesting as the classical ILS model. The same conclusion can be made for all job instances. For the remaining instances, the performance is quite far from the best known results from the literature. However, the median value of NILS is better than the one of NILS. For example, for and , NILS clearly outperforms NILS. The best results are reached when neutral walks on plateaus are allowed, which validates the use of neutrality.
When the results are compared with the degree of neutrality, we can give additional conclusions. For machine instances, the average degree of neutrality of solutions is quite low (around in average) [9]. However, the performance of NILS is better. This could be surprising at first sight. In fact, for those instances, the degree of neutrality is sufficient for the neutral random walks to cross a wide search space with the same fitness level. For the instances with a higher degree of neutrality, the density of portals on the plateaus decreases, and a pure random walk on a plateau is less efficient to find a portal quickly.
4.3 Influence of the Parameter
Since random neutral walks on a local optima plateau seem to be efficient in finding improving solutions, we can wonder whether the neutral walk should be large or not. The parameter corresponds to the maximal number of neutral steps allowed to move on the LOP. In this section, we study the performance using different values that depend on the number of jobs and so, that is related to the neighborhood size. Figure 1 shows the boxplot of fitness values for instances up to jobs and machines. For machine instances (top), the median fitness value starts by decreasing with the value, and then stabilizes: the results are equivalent for large values. As they are not distributed normally, the MannWhitney statistical test was used to check this tendency. The pairwise tests validate that the results are not statistically different for highest values. Increasing the value of the parameter does not deteriorate the performance. Therefore, it appears that this parameter can be set to a large value. For machine instances (bottom), best performances are always found by a NILS algorithm that exploit the LOP. The same happens for and (Figure 2 (a)). However, for and , the performance decreases with the value. This remark is confirmed by the MannWhitney statistical test. Therefore, in some cases, an overexploitation of the plateaus can become too costly. Here, we see that the tradeoff between neutrality exploitation and search space exploration should be considered carefully as it seems to exist a moresuited value for some instance sizes.
(a)  (b) 
(c)  (d) 
The performances after evaluations are better with the maximum number of neutral steps. Figure 2 (b) represents the median of the 30 fitness values according to different numbers of evaluations for and . For each possible number of evaluations, the median performance is always higher when the value is large. Moreover, the larger the problem size, the better the fitness value found using neutrality along the search (see Figure 2 (b)). Whatever the fitness level, from to , using neutrality clearly improves the performance.
4.4 Benefits and Costs of Neutral Moves
In order to deeply analyze the tradeoff between exploitation and exploration, tuned by the
parameter, we here estimate the benefit and the cost of neutral moves on the plateaus. Figure
2(c) computes the number of portals reached by the neutral walks for different values for and . When the value increases, the number of reached portals increases, even if the slope decreases. This indicates that the few additional portals reached with a larger value are interesting ones, since they lead to better solutions (Figure 1).The cost of the neutral moves is the total number of evaluations when the neutral walk fails to find a portal. Figure 2 (d) gives this number for and . Even if this cost increases with the value, we can notice that the slope decreases. This indicates that a larger neutral walk keeps interesting performance. Nevertheless, the slope of the cost is higher than the one of the benefit. It suggests that a tradeoff value can be reached in favor of a kickmove exploration for very large values.
5 Conclusion and Future Works
In this paper, we exploited the property of neutrality in order to design an efficient local search algorithm. Indeed, numerous combinatorial optimization problems contain a significant number of equivalent neighboring solutions. A Neutralitybased Iterated Local Search (NILS), that allows neutral walks to move on the plateaus, is proposed. The performance of the NILS algorithm has been experimented on the permutation flowshop scheduling problem by regarding the maximum length of neutral walks. As revealed in a previous study [9], the fitness landscape associated with this problem has a high degree of neutrality. Our experimental analysis shows that neutral walks allow to find improving solutions, and that the longer the neutral walk, the better the solution found. NILS is able to take advantage of applying neutral moves on plateaus, without any prohibitive additional cost. This can be explained by the relatively high density of portals over plateaus which can be reached by a random neutral walk.
In their iterated greedy algorithm proposed to solve the flowshop problem investigated in the paper, Ruiz et al. [8] combine a steepest descent local search with a destruction and a construction phases for perturbation, and an acceptance criterion based on a temperature. Let us remark that the authors suggest that the temperature level has not a significant influence on the overall performance. However, for the problem at hand, a solution has a non negligible number of neutral neighbors and the local optima plateaus are quite large. Thus, the destruction and construction phases can easily lead to a solution with the same fitness value. This perturbation can be compared to the NILS perturbation between a random neutral walk and a kickmove exploration. In the NILS algorithm, this tradeoff is explicitly tuned by a single parameter.
Considering those remarks together with the neutral exploitation of NILS, it suggests that the good performance of Ruiz et al.’s algorithm [8]
is probably due to the degree of neutrality of the flowshop scheduling problem. NILS uses this property more explicitly, with a lower implementation cost and less parameters.
This first version of the NILS algorithm explores the neighborhood in a random order, and both the components of hillclimbing and of random neutral walk are based on a firstimproving strategy. Future works will be dedicated to design efficient guiding methods based on the evolvability of solutions, in order to reach the portals more quickly. Moreover, a similar work will allow to better understand the dynamics of the NILS algorithm on other combinatorial optimization problems where a high degree of neutrality arises. Hopefully, this will allow us to propose an adaptive control of the single parameter of the NILS algorithm: the maximum length of the neutral walk.
References

[1]
Barnett, L.:
Netcrawling  optimal evolutionary search with neutral networks.
In: Congress on Evolutionary Computation CEC2001, IEEE Press (2001) 30–37
 [2] Stewart, T.: Extrema selection: Accelerated evolution on neutral networks. In: Congress on Evolutionary Computation CEC2001, IEEE Press (2001) 25–29
 [3] Verel, S., Collard, P., Clergue, M.: Scuba Search : when selection meets innovation. In: Congress on Evolutionary Computation CEC2004, IEEE Press (2004) 924 – 931
 [4] Johnson, S.M.: Optimal two and threestage production schedules with setup times included. Naval Research Logistics Quarterly 1 (1954) 61–68
 [5] Ruiz, R., Maroto, C.: A comprehensive review and evaluation of permutation flowshop heuristics. European Journal of Operational Research 165(2) (2005) 479–494
 [6] Stützle, T.: Applying iterated local search to the permutation flow shop problem. Technical Report AIDA9804, FG Intellektik, TU Darmstadt (1998)
 [7] Taillard, E.: Benchmarks for basic scheduling problems. European Journal of Operational Research 64 (1993) 278–285
 [8] Ruiz, R., Stützle, T.: A simple and effective iterated greedy algorithm for the permutation flowshop scheduling problem. European Journal of Operational Research 177(3) (2007) 2033–2049
 [9] Marmion, M.E., Dhaenens, C., Jourdan, L., Liefooghe, A., Verel, S.: On the neutrality of flowshop scheduling fitness landscapes. In: Learning and Intelligent OptimizatioN (LION 5). LNCS, SpringerVerlag (2011)
 [10] Ruiz, R., Maroto, C.: A comprehensive review and evaluation of permutation flowshop heuristics. European Journal of Operational Research 165(2) (2005) 479–494
 [11] Stadler, P.F. In: Towards a theory of landscapes. Volume 461. Springer Berlin / Heidelberg (1995) 78–163
 [12] Wilke, C.O.: Adaptative evolution on neutral networks. Bull. Math. Biol 63 (2001) 715–730
 [13] Verel, S., Collard, P., Tomassini, M., Vanneschi, L.: Fitness landscape of the cellular automata majority problem: view from the “Olympus”. Theor. Comp. Sci. 378 (2007) 54–77
 [14] Lourenço, H.R., Martin, O., Stützle, T.: Iterated local search. In: Handbook of Metaheuristics. Volume 57 of International Series in Operations Research & Management Science. Kluwer Academic Publishers (2002) 321–353
Comments
There are no comments yet.