an order effect of neighborhood structures in...

9
Research Article An Order Effect of Neighborhood Structures in Variable Neighborhood Search Algorithm for Minimizing the Makespan in an Identical Parallel Machine Scheduling Ibrahim Alharkan, 1 Khaled Bamatraf, 1 Mohammed A. Noman , 1 Husam Kaid , 1 Emad S. Abouel Nasr, 1,2 and Abdulaziz M. El-Tamimi 1 1 Industrial Engineering Department, College of Engineering, King Saud University, P.O. Box 800, Riyadh 11421, Saudi Arabia 2 Faculty of Engineering, Mechanical Engineering Department, Helwan University, Cairo 11732, Egypt Correspondence should be addressed to Mohammed A. Noman; [email protected] Received 11 September 2017; Revised 6 February 2018; Accepted 27 February 2018; Published 22 April 2018 Academic Editor: Ton D. Do Copyright © 2018 Ibrahim Alharkan et al. is 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. Variable neighborhood search (VNS) algorithm is proposed for scheduling identical parallel machine. e objective is to study the effect of adding a new neighborhood structure and changing the order of the neighborhood structures on minimizing the makespan. To enhance the quality of the final solution, a machine based encoding method and five neighborhood structures are used in VNS. Two initial solution methods which were used in two versions of improved VNS (IVNS) are employed, namely, longest processing time (LPT) initial solution, denoted as HIVNS, and random initial solution, denoted as RIVNS. e proposed versions are compared with LPT, simulated annealing (SA), genetic algorithm (GA), modified variable neighborhood search (MVNS), and improved variable neighborhood search (IVNS) algorithms from the literature. Computational results show that changing the order of neighborhood structures and adding a new neighborhood structure can yield a better solution in terms of average makespan. 1. Introduction Identical parallel machine scheduling (IPMS) with the objec- tive of minimizing the makespan is one of the combinational optimization problems. It is known to be NP-hard by Garey and Johnson [1] since it does not have a polynomial time algorithm. Exact algorithms such as branch and bound [2] and cutting plane algorithms [3] solve this type of IPM and find optimal solution for small size instances. As the problem size increases, the exact algorithms are inefficient and take much time to get a solution. at disadvantages bring a need for heuristics and meta- heuristics that give optimal or near optimal solution within a reasonable amount of time. Longest Processing Time Rule (LPT) proposed by Mokotoff [4] is the first heuristic applied in IPMS which has a tight worst case performance of bound of 4/3–1/3, where is the number of parallel machines. LPT is based on distributing jobs on machines according to maximum processing time and the remaining jobs go one by one to the least loaded machine until assigning all the jobs to the machines. e LPT heuristic performs well for makespan criteria but the solution obtained is oſten local optima. Later, Coffman et al. [5] proposed MULTIFIT algorithm that is based on techniques from bin-packing. Blackstone Jr. and Phillips [6] proposed a simple heuristic for improving LPT sequence by exchange jobs between processors to reduce makespan. Lee and Massey [7] combine two heuristics, LPT and MULTIFIT, to form a new one. e heuristic uses LPT heuristic as an initial solution for the MULTIFIT heuristic. e performance of the combined heuristic is better than LPT and the error bound is not worse than the MULTIFIT. Yue [8] proved the bound for MULTIFIT to be 13/11. Lee and Massey [9] extend the MULTIFIT algorithm and show that the error bound of implementing the algorithm is only 1/10. Garey and Johnson [1] proposed that 3-phase composite heuristic consists of constructive phase and two improvement phases with no preliminary sort of processing times. ey showed that their proposed heuristic is quicker than LPT. Ho and Hindawi Mathematical Problems in Engineering Volume 2018, Article ID 3586731, 8 pages https://doi.org/10.1155/2018/3586731

Upload: vuongbao

Post on 31-Jul-2018

214 views

Category:

Documents


0 download

TRANSCRIPT

Research ArticleAn Order Effect of Neighborhood Structures in VariableNeighborhood Search Algorithm for Minimizing the Makespanin an Identical Parallel Machine Scheduling

Ibrahim Alharkan1 Khaled Bamatraf1 Mohammed A Noman 1 Husam Kaid 1

Emad S Abouel Nasr12 and Abdulaziz M El-Tamimi1

1 Industrial Engineering Department College of Engineering King Saud University PO Box 800 Riyadh 11421 Saudi Arabia2Faculty of Engineering Mechanical Engineering Department Helwan University Cairo 11732 Egypt

Correspondence should be addressed to Mohammed A Noman mmohammed1ksuedusa

Received 11 September 2017 Revised 6 February 2018 Accepted 27 February 2018 Published 22 April 2018

Academic Editor Ton D Do

Copyright copy 2018 Ibrahim Alharkan et al This is an open access article distributed under the Creative Commons AttributionLicense which permits unrestricted use distribution and reproduction in any medium provided the original work is properlycited

Variable neighborhood search (VNS) algorithm is proposed for scheduling identical parallel machine The objective is to studythe effect of adding a new neighborhood structure and changing the order of the neighborhood structures on minimizing themakespan To enhance the quality of the final solution a machine based encoding method and five neighborhood structures areused inVNS Two initial solutionmethodswhichwere used in two versions of improvedVNS (IVNS) are employed namely longestprocessing time (LPT) initial solution denoted as HIVNS and random initial solution denoted as RIVNSThe proposed versionsare compared with LPT simulated annealing (SA) genetic algorithm (GA) modified variable neighborhood search (MVNS) andimproved variable neighborhood search (IVNS) algorithms from the literature Computational results show that changing the orderof neighborhood structures and adding a new neighborhood structure can yield a better solution in terms of average makespan

1 Introduction

Identical parallel machine scheduling (IPMS) with the objec-tive of minimizing the makespan is one of the combinationaloptimization problems It is known to be NP-hard by Gareyand Johnson [1] since it does not have a polynomial timealgorithm Exact algorithms such as branch and bound [2]and cutting plane algorithms [3] solve this type of IPM andfind optimal solution for small size instances As the problemsize increases the exact algorithms are inefficient and takemuch time to get a solution

That disadvantages bring a need for heuristics and meta-heuristics that give optimal or near optimal solution withina reasonable amount of time Longest Processing Time Rule(LPT) proposed by Mokotoff [4] is the first heuristic appliedin IPMS which has a tight worst case performance of boundof 43ndash13119898 where 119898 is the number of parallel machinesLPT is based on distributing jobs on machines according tomaximum processing time and the remaining jobs go one by

one to the least loaded machine until assigning all the jobs tothe machinesThe LPT heuristic performs well for makespancriteria but the solution obtained is often local optima LaterCoffman et al [5] proposed MULTIFIT algorithm that isbased on techniques from bin-packing Blackstone Jr andPhillips [6] proposed a simple heuristic for improving LPTsequence by exchange jobs between processors to reducemakespan Lee and Massey [7] combine two heuristics LPTand MULTIFIT to form a new one The heuristic uses LPTheuristic as an initial solution for the MULTIFIT heuristicTheperformance of the combined heuristic is better than LPTand the error bound is not worse than theMULTIFIT Yue [8]proved the bound for MULTIFIT to be 1311 Lee and Massey[9] extend the MULTIFIT algorithm and show that the errorbound of implementing the algorithm is only 110 Gareyand Johnson [1] proposed that 3-phase composite heuristicconsists of constructive phase and two improvement phaseswith no preliminary sort of processing times They showedthat their proposed heuristic is quicker than LPT Ho and

HindawiMathematical Problems in EngineeringVolume 2018 Article ID 3586731 8 pageshttpsdoiorg10115520183586731

2 Mathematical Problems in Engineering

Wong [10] introduce Two-Machine Optimal Schedulingwhich uses lexicographic searchTheir method performs bet-ter that LPTMULTIFIT andMULTIFIT extension algorithmand it takes less amount of CPU time than MULTIFIT andMULTIFIT extension algorithms

Riera et al [11] proposed two approximate algorithmsthat use LPT as an initial solution and compare themwith dynamic programming and MULTIFIT algorithmsAlgorithm 1 uses exchange between two jobs to improvethe makespan Algorithm 2 schedules a job such that thecompletion time and process time of the selected job arenear the bound Their second algorithm is compared withMULTIFIT algorithms and results showed similarity to theMULTIFIT algorithm but their algorithm reduces CPUtime with respect to MULTIFIT heuristic Cheng and Gen[12] applied memetic algorithm to minimize the maximumweighted absolute lateness on PMS and showed that it out-performs genetic algorithm and the conventional heuristicsGhomi and Ghazvini [13] proposed a pairwise interchangealgorithm and it gave near optimal solution in a short periodof timeMin andCheng [14] proposed a genetic algorithmGAusing machine code They showed that GA outperforms LPTand SA and is suitable for large scale IPMS problems Guptaand Ruiz-Torres [15] proposed a LISTFIT heuristic basedon bin-backing and list scheduling The LISTFIT generatean optimal or near optimal solution and outperforms LPTMULTIFIT and COMBINE heuristics Costa et al [16]proposed algorithm inspired by the immune systems ofvertebrate animals Lee et al [17] proposed a simulatedannealing (SA) approach for makespan minimization onIPMS It chooses LPT as an initial solution Computationalresults showed that the SAheuristic outperforms the LISTFITand pairwise interchange (PI) algorithms Moreover it isefficient for large scale problems Tang and Luo [18] proposea new ILS algorithm combining with a variable number ofcyclic exchanges Experiments show that the algorithm isefficient for 119875119898 119862max Akyol and Bayhan [19] proposeda dynamical neural network that employs parameters oftime varying penalty The simulation results showed thatthe proposed algorithm generated feasible solutions and itfound better makespan when compared to LPT Kashan andKarimi [20] presented discrete particle swarm optimization(DPSO) algorithm for makespan minimization Computa-tional results showed that hybridized DPSO (HDPSO) algo-rithmoutperforms both SA andDPSO algorithms Sevkli andUysal [21] proposed modified variable neighborhood search(MVNS) which is based on exchange and move neighbor-hood structures Computational results demonstrated thatthe proposed algorithm outperforms both GA and LPTalgorithms Min and Cheng [14] proposed a harmony search(HS) algorithmwith dynamic subpopulation (DSHS) Resultsshow that DSHS algorithm outperforms SA and HDPSOfor many instances Moreover the execution time is lessthan 1 sec for all computations Chen et al [22] proposeddiscrete harmony search (DHS) algorithm that uses discreteencoding scheme to initialize the harmony memory (HM)then the improvisation scheme for generating new harmonyis redefined for suitability for solving the combinational opti-mization problem In addition the study made hybridizing a

local search method with DHS to increase the speed of localsearch Computational results show that theDHS algorithm isvery competitive when compared with other heuristics in theliterature Jing and Jun-qing [23] proposed efficient variableneighborhood search that uses four neighborhood structuresand has two versions One version uses LPT sequence as aninitial solution The other version uses random sequence asan initial solution A computational result demonstrates thatEVNS is efficient in searching global or near global optimumM Sevkli and A Z Sevkli [24] proposed stochasticallyperturbed particle swarm optimization algorithm (SPPSO)The algorithm compared two recent PSO algorithms Itis concluded that SPPSO algorithm has produced betterresults than DPSO and PSOspv in terms of the optimalsolutions number Laha [25] proposed an improved simulatedannealing (SA) heuristic Computational results show thatthe proposed heuristic is better than that produced by thebest-known heuristic in the literature Other advantages ofit are the ease of implementation In this paper the proposedalgorithm of Jing and Jun-qing [23] in their paper ldquoefficientvariable neighborhood search for identical parallel machinesschedulingrdquo is used with some changes on it One of thechanges is changing in the order of the neighborhood struc-tures and the other change is adding another neighborhoodstructure to get five neighborhood structures in our proposedalgorithm

The remaining sections of this paper are organized asfollows In Section 2 a brief description of IPMS problemis mentioned In Section 3 the steps of proposed algorithmare described in detail and the neighborhood structuresof this proposed algorithm are explained In Section 4computational results are discussed Conclusion is made inSection 5

2 Problem Description

The identical parallel machine scheduling (IPMS) problemcan be described as follows

A set 119899 of an independent jobs 119869 = 1198691 1198692 119869119899 tobe processed on 119898 identical parallel machines 119872 = 11987211198722 119872119898with the processing time of job 119894on any identicalmachine is given by 119901119894

A job can only be processed on one machine simultane-ously and a machine cannot process more than one job ata time Priority and precedence constraints are not allowedThere is no job cancellation and a job completes its processingon a machine without interruption

The objective is to minimize the total completion timeldquothe makespanrdquo of scheduling jobs on the machines

This scheduling problem can be described by a triple 120572 |120573 | 120574 as follows

119875119898 119862max (1)

where 119875 indicates parallel machine environment119898 indicatesnumber of machines 120573 indicates no constraints in thisproblem and 119862max indicates that the objective is to minimizethe makespan

Mathematical Problems in Engineering 3

1st version

LPT initial solutionHIVNS1

Random initial solutionRIVNS1

LPT initial solutionHIVNS2

Random initial solutionRIVNS2

2nd version

Figure 1 The two versions of the proposed algorithms

This problem is interesting because minimizing themakespan has the effect of balancing the load over the variousmachines which is an important goal in practice

3 Development of the Proposed(IVNS) Algorithm

31 Basic VNS Variable neighborhood search (VNS) is ametaheuristic proposed by Mladenovic and Hansen [26] toenhance the solution quality by systematic neighborhoodschanges The main VNS algorithm steps can be summarizedas follows Initialization choose the neighborhood structuresset (NK1015840) 119896 = 1 2 1198961015840max obtain an initial solution andselect a stopping condition Repeat the next steps until thestopping condition is satisfied

(1) Set 119896 larr 1(2) Repeat the following steps until 119896 = 119896max

(a) Shaking generate a point 1199091015840 at random from the119896th neighborhood of 119909 (1199091015840 isin 119873119896(119909))

(b) Local search apply some local search methodwith 1199091015840 as initial solution denote with 11990910158401015840 the soobtained local optimum

(c) Move or not if the local optimum 11990910158401015840 is betterthan the incumbent 119909 move there (119909 larr 11990910158401015840)and continue the search with 1198731 (119896 larr 1)otherwise set 119896 larr 119896 + 1 improved variableneighborhood search (IVNS) algorithm

As we mentioned earlier the proposed algorithm is anaddition of the proposed algorithm of Jing and Jun-qing [23]

The proposed algorithms have two versions and twotypes for each version as shown in Figure 1 In the firstversion a new neighborhood structure was added to thefour neighborhood structures which are proposed by Jingand Jun-qing [23] while in the second version the order ofthese neighborhood structures was changed Both versionsuse LPT [20] and random initial solutions and are referredto as ldquoHIVNSrdquo and ldquoRIVNSrdquo respectively All these versionsof the proposed algorithm use the same five neighborhoodstructures These neighborhood structures will be discussedin the following section

32 Neighborhood Structures Determining the neighbor-hood structures is critical in the VNS algorithm To enhance

the local searching abilities five different kinds of neigh-borhoods are utilized to find better solutions on a givenschedule in the proposed algorithm which are designedbased on such an idea that a given solution can be improvedby moving or swapping jobs between the problem machines(themachines with their finished time equal to themakespanof the solution) and any other nonproblem machines (themachines with their finished time less than the makespan ofthe solution)

The five neighborhood structures are illustrated as fol-lows

(1) Move move a job 119869119894 from 119872119901 to 119872119899119901 if condition(119862119872119901 minus 119862119872119899119901 gt 119901119894) is satisfied

(2) Exchange 1 exchange a job 119869119894 selected from119872119901 withanother job 119869119895 selected from119872119899119901 if (119901119894 minus 119901119895 gt 0) and(119862119872119901 minus 119862119872119899119901 gt 119901119894 minus 119901119895)

(3) Exchange 2 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with one job 119869119896 selected from119872119899119901 if (119901119894 +119901119895 minus 119901119896 gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 + 119901119895 minus 119901119896)

(4) Exchange 3 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with two jobs 119869119896 and 119869119905 selected from119872119899119901if (119901119894 + 119901119895 minus (119901119896 + 119901119905) gt 0) and (119862119872119901 minus 119862119872119899119901 gt119901119894 + 119901119895 minus (119901119896 + 119901119905))

(5) Exchange 4 exchange one job 119869119894 selected from 119872119901with two jobs 119869119894 and 119869119896 selected from 119872119899119901 if (119901119894 minus(119901119895 + 119901119896) gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 minus (119901119895 + 119901119896))

The orders assigned to the types proposed of the algo-rithm are as follows

(1) The order of ldquoHIVNS1rdquo and ldquoRIVNS1rdquo is ldquomoveexchange 1 exchange 2 exchange 3 and exchange 4rdquo

(2) The order of ldquoHIVNS2rdquo and ldquoRIVNS2rdquo is ldquoexchange3 exchange 1 move exchange 2 and exchange4rdquo Improved VNS (IVNS) flow chart is shown asFigure 2

33 Steps of IVNS The steps of IVNS for ldquoHIVNS1rdquo andldquoRIVNS1rdquo are shown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

4 Mathematical Problems in Engineering

No

Shake procedure get a new start

Yes

Yes

Neighborhood search

NoNo

Yes

Initialization set initial solution x0 max iteration

Output x0

point x

k = 1

Local search get a solution

k = 1

x = x

k = k + 1

if k gt kGR

x isin Nk(x)

i = i + 1

then x0 = x

num (MaxIterNum) i = 0

i lt MaxIterNum

if (f(x) lt f(x0))

If (f(x) lt f(x))

Figure 2 Flow chart of the basic VNS algorithm

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212Repeat

Switch (119870) 119870 = 1 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break

119870 = 3 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 4 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

The steps of IVNS for ldquoHIVNS2rdquo and ldquoRIVNS2rdquo areshown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Mathematical Problems in Engineering 5

Table 1 Number of machines and jobs

Number of machines 2 5 10 20Number of jobs 20 50 100 200 20 50 100 200 20 50 100 200 50 100 200

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212

Repeat

Switch (119870) 119870 = 1 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break119870 = 3 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 4 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

4 Computational Results and Comparison

In this section the results of two versions of the proposedalgorithm were compared with LPT [1] SA [17] GA [14]MVNS [21] and IVNS [23] algorithms from the literatureThe two versions of the improved variable neighborhoodsearch algorithms ldquoHIVNS1 and RIVNS1rdquo and ldquoHIVNS2 andRIVNS2rdquo were coded in MATLAB R2012a and executed oni5 CPU 5GHz with 6GB of RAM All of them were stoppedafter getting the lower bound or running for 100 iterations forRIVNS1 and RIVNS2 The number of machines and numberof jobs are shown in Table 1

The processing time of the jobs is the same as Jingand Jun-qing [23] As a result 15 instances were conductedand each instance was conducted with 10 generations ofdifferent processing times The total is 150 instances The

094098102106

11

AlgorithmsAverageLower boundMaximum average

mak

espa

n m

eans

Aver

age o

f

LPT SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

Figure 3 Histogram of averages of makespan means

performance of the algorithms is measured with respectto the average makespan (mean) and average CPU time(Avg time in second) The ldquomeanrdquo performance is a relativequality measure of solutions computed by 119862LB where 119862is the average makespan obtained for each instance with 10generations given by the algorithm and LB is the lower boundof the instance calculated in equation that mentioned inSection 33 Step 2 The ldquoAvg timerdquo refers to the total time ittakes for the algorithm to finish the solution Table 2 presentsthe results of the previous algorithms from the literaturewhileTable 3 presents the results from the proposed algorithms Bycomparing the results of the average means of the makespanIt is obvious that the proposed algorithms outperform all thealgorithms in Table 2 from the literature It is worth notingthat for each instance the proposed algorithms obtain noworse than algorithms inTable 2 except at 10machines and 20jobs instance in only HIVNS algorithm and that is due to thedifficulty facing the proposed algorithms to get a lower boundwhen the difference between the number ofmachines and thenumber of jobs is relatively small In addition by comparingthe two proposed algorithms we can see that the versionshave the same average means of the makespan in case ofrandom initial solution while as the 2nd version outperformsthe 1st version in case of LPT initial solution in both averagemeans of makespan and average CPU time

Figure 3 shows the averages of makespan mean themaximum averages for each algorithm and the lower boundIt can be observed that the two proposed versions algorithmshave makespan means averages which closed to the lowerbound especially in RIVNS1 and RIVNS2 that have the sameaverage (10054)

Figure 4 shows the averages of Avg time means for allalgorithms We can see that RIVNS1 and RIVNS2 have thesmallest Avg time which are 0008 and 0007 respectivelyMoreover it can be observed that the Avg time of HIVNS1and HIVNS2 is much higher than HIVNS because in thispaper Matlab is used to construct the code of HIVNS1 and

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

2 Mathematical Problems in Engineering

Wong [10] introduce Two-Machine Optimal Schedulingwhich uses lexicographic searchTheir method performs bet-ter that LPTMULTIFIT andMULTIFIT extension algorithmand it takes less amount of CPU time than MULTIFIT andMULTIFIT extension algorithms

Riera et al [11] proposed two approximate algorithmsthat use LPT as an initial solution and compare themwith dynamic programming and MULTIFIT algorithmsAlgorithm 1 uses exchange between two jobs to improvethe makespan Algorithm 2 schedules a job such that thecompletion time and process time of the selected job arenear the bound Their second algorithm is compared withMULTIFIT algorithms and results showed similarity to theMULTIFIT algorithm but their algorithm reduces CPUtime with respect to MULTIFIT heuristic Cheng and Gen[12] applied memetic algorithm to minimize the maximumweighted absolute lateness on PMS and showed that it out-performs genetic algorithm and the conventional heuristicsGhomi and Ghazvini [13] proposed a pairwise interchangealgorithm and it gave near optimal solution in a short periodof timeMin andCheng [14] proposed a genetic algorithmGAusing machine code They showed that GA outperforms LPTand SA and is suitable for large scale IPMS problems Guptaand Ruiz-Torres [15] proposed a LISTFIT heuristic basedon bin-backing and list scheduling The LISTFIT generatean optimal or near optimal solution and outperforms LPTMULTIFIT and COMBINE heuristics Costa et al [16]proposed algorithm inspired by the immune systems ofvertebrate animals Lee et al [17] proposed a simulatedannealing (SA) approach for makespan minimization onIPMS It chooses LPT as an initial solution Computationalresults showed that the SAheuristic outperforms the LISTFITand pairwise interchange (PI) algorithms Moreover it isefficient for large scale problems Tang and Luo [18] proposea new ILS algorithm combining with a variable number ofcyclic exchanges Experiments show that the algorithm isefficient for 119875119898 119862max Akyol and Bayhan [19] proposeda dynamical neural network that employs parameters oftime varying penalty The simulation results showed thatthe proposed algorithm generated feasible solutions and itfound better makespan when compared to LPT Kashan andKarimi [20] presented discrete particle swarm optimization(DPSO) algorithm for makespan minimization Computa-tional results showed that hybridized DPSO (HDPSO) algo-rithmoutperforms both SA andDPSO algorithms Sevkli andUysal [21] proposed modified variable neighborhood search(MVNS) which is based on exchange and move neighbor-hood structures Computational results demonstrated thatthe proposed algorithm outperforms both GA and LPTalgorithms Min and Cheng [14] proposed a harmony search(HS) algorithmwith dynamic subpopulation (DSHS) Resultsshow that DSHS algorithm outperforms SA and HDPSOfor many instances Moreover the execution time is lessthan 1 sec for all computations Chen et al [22] proposeddiscrete harmony search (DHS) algorithm that uses discreteencoding scheme to initialize the harmony memory (HM)then the improvisation scheme for generating new harmonyis redefined for suitability for solving the combinational opti-mization problem In addition the study made hybridizing a

local search method with DHS to increase the speed of localsearch Computational results show that theDHS algorithm isvery competitive when compared with other heuristics in theliterature Jing and Jun-qing [23] proposed efficient variableneighborhood search that uses four neighborhood structuresand has two versions One version uses LPT sequence as aninitial solution The other version uses random sequence asan initial solution A computational result demonstrates thatEVNS is efficient in searching global or near global optimumM Sevkli and A Z Sevkli [24] proposed stochasticallyperturbed particle swarm optimization algorithm (SPPSO)The algorithm compared two recent PSO algorithms Itis concluded that SPPSO algorithm has produced betterresults than DPSO and PSOspv in terms of the optimalsolutions number Laha [25] proposed an improved simulatedannealing (SA) heuristic Computational results show thatthe proposed heuristic is better than that produced by thebest-known heuristic in the literature Other advantages ofit are the ease of implementation In this paper the proposedalgorithm of Jing and Jun-qing [23] in their paper ldquoefficientvariable neighborhood search for identical parallel machinesschedulingrdquo is used with some changes on it One of thechanges is changing in the order of the neighborhood struc-tures and the other change is adding another neighborhoodstructure to get five neighborhood structures in our proposedalgorithm

The remaining sections of this paper are organized asfollows In Section 2 a brief description of IPMS problemis mentioned In Section 3 the steps of proposed algorithmare described in detail and the neighborhood structuresof this proposed algorithm are explained In Section 4computational results are discussed Conclusion is made inSection 5

2 Problem Description

The identical parallel machine scheduling (IPMS) problemcan be described as follows

A set 119899 of an independent jobs 119869 = 1198691 1198692 119869119899 tobe processed on 119898 identical parallel machines 119872 = 11987211198722 119872119898with the processing time of job 119894on any identicalmachine is given by 119901119894

A job can only be processed on one machine simultane-ously and a machine cannot process more than one job ata time Priority and precedence constraints are not allowedThere is no job cancellation and a job completes its processingon a machine without interruption

The objective is to minimize the total completion timeldquothe makespanrdquo of scheduling jobs on the machines

This scheduling problem can be described by a triple 120572 |120573 | 120574 as follows

119875119898 119862max (1)

where 119875 indicates parallel machine environment119898 indicatesnumber of machines 120573 indicates no constraints in thisproblem and 119862max indicates that the objective is to minimizethe makespan

Mathematical Problems in Engineering 3

1st version

LPT initial solutionHIVNS1

Random initial solutionRIVNS1

LPT initial solutionHIVNS2

Random initial solutionRIVNS2

2nd version

Figure 1 The two versions of the proposed algorithms

This problem is interesting because minimizing themakespan has the effect of balancing the load over the variousmachines which is an important goal in practice

3 Development of the Proposed(IVNS) Algorithm

31 Basic VNS Variable neighborhood search (VNS) is ametaheuristic proposed by Mladenovic and Hansen [26] toenhance the solution quality by systematic neighborhoodschanges The main VNS algorithm steps can be summarizedas follows Initialization choose the neighborhood structuresset (NK1015840) 119896 = 1 2 1198961015840max obtain an initial solution andselect a stopping condition Repeat the next steps until thestopping condition is satisfied

(1) Set 119896 larr 1(2) Repeat the following steps until 119896 = 119896max

(a) Shaking generate a point 1199091015840 at random from the119896th neighborhood of 119909 (1199091015840 isin 119873119896(119909))

(b) Local search apply some local search methodwith 1199091015840 as initial solution denote with 11990910158401015840 the soobtained local optimum

(c) Move or not if the local optimum 11990910158401015840 is betterthan the incumbent 119909 move there (119909 larr 11990910158401015840)and continue the search with 1198731 (119896 larr 1)otherwise set 119896 larr 119896 + 1 improved variableneighborhood search (IVNS) algorithm

As we mentioned earlier the proposed algorithm is anaddition of the proposed algorithm of Jing and Jun-qing [23]

The proposed algorithms have two versions and twotypes for each version as shown in Figure 1 In the firstversion a new neighborhood structure was added to thefour neighborhood structures which are proposed by Jingand Jun-qing [23] while in the second version the order ofthese neighborhood structures was changed Both versionsuse LPT [20] and random initial solutions and are referredto as ldquoHIVNSrdquo and ldquoRIVNSrdquo respectively All these versionsof the proposed algorithm use the same five neighborhoodstructures These neighborhood structures will be discussedin the following section

32 Neighborhood Structures Determining the neighbor-hood structures is critical in the VNS algorithm To enhance

the local searching abilities five different kinds of neigh-borhoods are utilized to find better solutions on a givenschedule in the proposed algorithm which are designedbased on such an idea that a given solution can be improvedby moving or swapping jobs between the problem machines(themachines with their finished time equal to themakespanof the solution) and any other nonproblem machines (themachines with their finished time less than the makespan ofthe solution)

The five neighborhood structures are illustrated as fol-lows

(1) Move move a job 119869119894 from 119872119901 to 119872119899119901 if condition(119862119872119901 minus 119862119872119899119901 gt 119901119894) is satisfied

(2) Exchange 1 exchange a job 119869119894 selected from119872119901 withanother job 119869119895 selected from119872119899119901 if (119901119894 minus 119901119895 gt 0) and(119862119872119901 minus 119862119872119899119901 gt 119901119894 minus 119901119895)

(3) Exchange 2 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with one job 119869119896 selected from119872119899119901 if (119901119894 +119901119895 minus 119901119896 gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 + 119901119895 minus 119901119896)

(4) Exchange 3 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with two jobs 119869119896 and 119869119905 selected from119872119899119901if (119901119894 + 119901119895 minus (119901119896 + 119901119905) gt 0) and (119862119872119901 minus 119862119872119899119901 gt119901119894 + 119901119895 minus (119901119896 + 119901119905))

(5) Exchange 4 exchange one job 119869119894 selected from 119872119901with two jobs 119869119894 and 119869119896 selected from 119872119899119901 if (119901119894 minus(119901119895 + 119901119896) gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 minus (119901119895 + 119901119896))

The orders assigned to the types proposed of the algo-rithm are as follows

(1) The order of ldquoHIVNS1rdquo and ldquoRIVNS1rdquo is ldquomoveexchange 1 exchange 2 exchange 3 and exchange 4rdquo

(2) The order of ldquoHIVNS2rdquo and ldquoRIVNS2rdquo is ldquoexchange3 exchange 1 move exchange 2 and exchange4rdquo Improved VNS (IVNS) flow chart is shown asFigure 2

33 Steps of IVNS The steps of IVNS for ldquoHIVNS1rdquo andldquoRIVNS1rdquo are shown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

4 Mathematical Problems in Engineering

No

Shake procedure get a new start

Yes

Yes

Neighborhood search

NoNo

Yes

Initialization set initial solution x0 max iteration

Output x0

point x

k = 1

Local search get a solution

k = 1

x = x

k = k + 1

if k gt kGR

x isin Nk(x)

i = i + 1

then x0 = x

num (MaxIterNum) i = 0

i lt MaxIterNum

if (f(x) lt f(x0))

If (f(x) lt f(x))

Figure 2 Flow chart of the basic VNS algorithm

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212Repeat

Switch (119870) 119870 = 1 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break

119870 = 3 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 4 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

The steps of IVNS for ldquoHIVNS2rdquo and ldquoRIVNS2rdquo areshown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Mathematical Problems in Engineering 5

Table 1 Number of machines and jobs

Number of machines 2 5 10 20Number of jobs 20 50 100 200 20 50 100 200 20 50 100 200 50 100 200

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212

Repeat

Switch (119870) 119870 = 1 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break119870 = 3 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 4 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

4 Computational Results and Comparison

In this section the results of two versions of the proposedalgorithm were compared with LPT [1] SA [17] GA [14]MVNS [21] and IVNS [23] algorithms from the literatureThe two versions of the improved variable neighborhoodsearch algorithms ldquoHIVNS1 and RIVNS1rdquo and ldquoHIVNS2 andRIVNS2rdquo were coded in MATLAB R2012a and executed oni5 CPU 5GHz with 6GB of RAM All of them were stoppedafter getting the lower bound or running for 100 iterations forRIVNS1 and RIVNS2 The number of machines and numberof jobs are shown in Table 1

The processing time of the jobs is the same as Jingand Jun-qing [23] As a result 15 instances were conductedand each instance was conducted with 10 generations ofdifferent processing times The total is 150 instances The

094098102106

11

AlgorithmsAverageLower boundMaximum average

mak

espa

n m

eans

Aver

age o

f

LPT SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

Figure 3 Histogram of averages of makespan means

performance of the algorithms is measured with respectto the average makespan (mean) and average CPU time(Avg time in second) The ldquomeanrdquo performance is a relativequality measure of solutions computed by 119862LB where 119862is the average makespan obtained for each instance with 10generations given by the algorithm and LB is the lower boundof the instance calculated in equation that mentioned inSection 33 Step 2 The ldquoAvg timerdquo refers to the total time ittakes for the algorithm to finish the solution Table 2 presentsthe results of the previous algorithms from the literaturewhileTable 3 presents the results from the proposed algorithms Bycomparing the results of the average means of the makespanIt is obvious that the proposed algorithms outperform all thealgorithms in Table 2 from the literature It is worth notingthat for each instance the proposed algorithms obtain noworse than algorithms inTable 2 except at 10machines and 20jobs instance in only HIVNS algorithm and that is due to thedifficulty facing the proposed algorithms to get a lower boundwhen the difference between the number ofmachines and thenumber of jobs is relatively small In addition by comparingthe two proposed algorithms we can see that the versionshave the same average means of the makespan in case ofrandom initial solution while as the 2nd version outperformsthe 1st version in case of LPT initial solution in both averagemeans of makespan and average CPU time

Figure 3 shows the averages of makespan mean themaximum averages for each algorithm and the lower boundIt can be observed that the two proposed versions algorithmshave makespan means averages which closed to the lowerbound especially in RIVNS1 and RIVNS2 that have the sameaverage (10054)

Figure 4 shows the averages of Avg time means for allalgorithms We can see that RIVNS1 and RIVNS2 have thesmallest Avg time which are 0008 and 0007 respectivelyMoreover it can be observed that the Avg time of HIVNS1and HIVNS2 is much higher than HIVNS because in thispaper Matlab is used to construct the code of HIVNS1 and

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

Mathematical Problems in Engineering 3

1st version

LPT initial solutionHIVNS1

Random initial solutionRIVNS1

LPT initial solutionHIVNS2

Random initial solutionRIVNS2

2nd version

Figure 1 The two versions of the proposed algorithms

This problem is interesting because minimizing themakespan has the effect of balancing the load over the variousmachines which is an important goal in practice

3 Development of the Proposed(IVNS) Algorithm

31 Basic VNS Variable neighborhood search (VNS) is ametaheuristic proposed by Mladenovic and Hansen [26] toenhance the solution quality by systematic neighborhoodschanges The main VNS algorithm steps can be summarizedas follows Initialization choose the neighborhood structuresset (NK1015840) 119896 = 1 2 1198961015840max obtain an initial solution andselect a stopping condition Repeat the next steps until thestopping condition is satisfied

(1) Set 119896 larr 1(2) Repeat the following steps until 119896 = 119896max

(a) Shaking generate a point 1199091015840 at random from the119896th neighborhood of 119909 (1199091015840 isin 119873119896(119909))

(b) Local search apply some local search methodwith 1199091015840 as initial solution denote with 11990910158401015840 the soobtained local optimum

(c) Move or not if the local optimum 11990910158401015840 is betterthan the incumbent 119909 move there (119909 larr 11990910158401015840)and continue the search with 1198731 (119896 larr 1)otherwise set 119896 larr 119896 + 1 improved variableneighborhood search (IVNS) algorithm

As we mentioned earlier the proposed algorithm is anaddition of the proposed algorithm of Jing and Jun-qing [23]

The proposed algorithms have two versions and twotypes for each version as shown in Figure 1 In the firstversion a new neighborhood structure was added to thefour neighborhood structures which are proposed by Jingand Jun-qing [23] while in the second version the order ofthese neighborhood structures was changed Both versionsuse LPT [20] and random initial solutions and are referredto as ldquoHIVNSrdquo and ldquoRIVNSrdquo respectively All these versionsof the proposed algorithm use the same five neighborhoodstructures These neighborhood structures will be discussedin the following section

32 Neighborhood Structures Determining the neighbor-hood structures is critical in the VNS algorithm To enhance

the local searching abilities five different kinds of neigh-borhoods are utilized to find better solutions on a givenschedule in the proposed algorithm which are designedbased on such an idea that a given solution can be improvedby moving or swapping jobs between the problem machines(themachines with their finished time equal to themakespanof the solution) and any other nonproblem machines (themachines with their finished time less than the makespan ofthe solution)

The five neighborhood structures are illustrated as fol-lows

(1) Move move a job 119869119894 from 119872119901 to 119872119899119901 if condition(119862119872119901 minus 119862119872119899119901 gt 119901119894) is satisfied

(2) Exchange 1 exchange a job 119869119894 selected from119872119901 withanother job 119869119895 selected from119872119899119901 if (119901119894 minus 119901119895 gt 0) and(119862119872119901 minus 119862119872119899119901 gt 119901119894 minus 119901119895)

(3) Exchange 2 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with one job 119869119896 selected from119872119899119901 if (119901119894 +119901119895 minus 119901119896 gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 + 119901119895 minus 119901119896)

(4) Exchange 3 exchange two jobs 119869119894 and 119869119895 selectedfrom119872119901 with two jobs 119869119896 and 119869119905 selected from119872119899119901if (119901119894 + 119901119895 minus (119901119896 + 119901119905) gt 0) and (119862119872119901 minus 119862119872119899119901 gt119901119894 + 119901119895 minus (119901119896 + 119901119905))

(5) Exchange 4 exchange one job 119869119894 selected from 119872119901with two jobs 119869119894 and 119869119896 selected from 119872119899119901 if (119901119894 minus(119901119895 + 119901119896) gt 0) and (119862119872119901 minus 119862119872119899119901 gt 119901119894 minus (119901119895 + 119901119896))

The orders assigned to the types proposed of the algo-rithm are as follows

(1) The order of ldquoHIVNS1rdquo and ldquoRIVNS1rdquo is ldquomoveexchange 1 exchange 2 exchange 3 and exchange 4rdquo

(2) The order of ldquoHIVNS2rdquo and ldquoRIVNS2rdquo is ldquoexchange3 exchange 1 move exchange 2 and exchange4rdquo Improved VNS (IVNS) flow chart is shown asFigure 2

33 Steps of IVNS The steps of IVNS for ldquoHIVNS1rdquo andldquoRIVNS1rdquo are shown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

4 Mathematical Problems in Engineering

No

Shake procedure get a new start

Yes

Yes

Neighborhood search

NoNo

Yes

Initialization set initial solution x0 max iteration

Output x0

point x

k = 1

Local search get a solution

k = 1

x = x

k = k + 1

if k gt kGR

x isin Nk(x)

i = i + 1

then x0 = x

num (MaxIterNum) i = 0

i lt MaxIterNum

if (f(x) lt f(x0))

If (f(x) lt f(x))

Figure 2 Flow chart of the basic VNS algorithm

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212Repeat

Switch (119870) 119870 = 1 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break

119870 = 3 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 4 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

The steps of IVNS for ldquoHIVNS2rdquo and ldquoRIVNS2rdquo areshown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Mathematical Problems in Engineering 5

Table 1 Number of machines and jobs

Number of machines 2 5 10 20Number of jobs 20 50 100 200 20 50 100 200 20 50 100 200 50 100 200

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212

Repeat

Switch (119870) 119870 = 1 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break119870 = 3 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 4 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

4 Computational Results and Comparison

In this section the results of two versions of the proposedalgorithm were compared with LPT [1] SA [17] GA [14]MVNS [21] and IVNS [23] algorithms from the literatureThe two versions of the improved variable neighborhoodsearch algorithms ldquoHIVNS1 and RIVNS1rdquo and ldquoHIVNS2 andRIVNS2rdquo were coded in MATLAB R2012a and executed oni5 CPU 5GHz with 6GB of RAM All of them were stoppedafter getting the lower bound or running for 100 iterations forRIVNS1 and RIVNS2 The number of machines and numberof jobs are shown in Table 1

The processing time of the jobs is the same as Jingand Jun-qing [23] As a result 15 instances were conductedand each instance was conducted with 10 generations ofdifferent processing times The total is 150 instances The

094098102106

11

AlgorithmsAverageLower boundMaximum average

mak

espa

n m

eans

Aver

age o

f

LPT SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

Figure 3 Histogram of averages of makespan means

performance of the algorithms is measured with respectto the average makespan (mean) and average CPU time(Avg time in second) The ldquomeanrdquo performance is a relativequality measure of solutions computed by 119862LB where 119862is the average makespan obtained for each instance with 10generations given by the algorithm and LB is the lower boundof the instance calculated in equation that mentioned inSection 33 Step 2 The ldquoAvg timerdquo refers to the total time ittakes for the algorithm to finish the solution Table 2 presentsthe results of the previous algorithms from the literaturewhileTable 3 presents the results from the proposed algorithms Bycomparing the results of the average means of the makespanIt is obvious that the proposed algorithms outperform all thealgorithms in Table 2 from the literature It is worth notingthat for each instance the proposed algorithms obtain noworse than algorithms inTable 2 except at 10machines and 20jobs instance in only HIVNS algorithm and that is due to thedifficulty facing the proposed algorithms to get a lower boundwhen the difference between the number ofmachines and thenumber of jobs is relatively small In addition by comparingthe two proposed algorithms we can see that the versionshave the same average means of the makespan in case ofrandom initial solution while as the 2nd version outperformsthe 1st version in case of LPT initial solution in both averagemeans of makespan and average CPU time

Figure 3 shows the averages of makespan mean themaximum averages for each algorithm and the lower boundIt can be observed that the two proposed versions algorithmshave makespan means averages which closed to the lowerbound especially in RIVNS1 and RIVNS2 that have the sameaverage (10054)

Figure 4 shows the averages of Avg time means for allalgorithms We can see that RIVNS1 and RIVNS2 have thesmallest Avg time which are 0008 and 0007 respectivelyMoreover it can be observed that the Avg time of HIVNS1and HIVNS2 is much higher than HIVNS because in thispaper Matlab is used to construct the code of HIVNS1 and

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

4 Mathematical Problems in Engineering

No

Shake procedure get a new start

Yes

Yes

Neighborhood search

NoNo

Yes

Initialization set initial solution x0 max iteration

Output x0

point x

k = 1

Local search get a solution

k = 1

x = x

k = k + 1

if k gt kGR

x isin Nk(x)

i = i + 1

then x0 = x

num (MaxIterNum) i = 0

i lt MaxIterNum

if (f(x) lt f(x0))

If (f(x) lt f(x))

Figure 2 Flow chart of the basic VNS algorithm

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212Repeat

Switch (119870) 119870 = 1 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break

119870 = 3 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 4 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

The steps of IVNS for ldquoHIVNS2rdquo and ldquoRIVNS2rdquo areshown as follows

Step 1 Generate initial schedule 119883 (generated randomly orobtain from the LPT rule) MaxIterNum = 100 119894 = 0 and119896max = 5

Mathematical Problems in Engineering 5

Table 1 Number of machines and jobs

Number of machines 2 5 10 20Number of jobs 20 50 100 200 20 50 100 200 20 50 100 200 50 100 200

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212

Repeat

Switch (119870) 119870 = 1 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break119870 = 3 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 4 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

4 Computational Results and Comparison

In this section the results of two versions of the proposedalgorithm were compared with LPT [1] SA [17] GA [14]MVNS [21] and IVNS [23] algorithms from the literatureThe two versions of the improved variable neighborhoodsearch algorithms ldquoHIVNS1 and RIVNS1rdquo and ldquoHIVNS2 andRIVNS2rdquo were coded in MATLAB R2012a and executed oni5 CPU 5GHz with 6GB of RAM All of them were stoppedafter getting the lower bound or running for 100 iterations forRIVNS1 and RIVNS2 The number of machines and numberof jobs are shown in Table 1

The processing time of the jobs is the same as Jingand Jun-qing [23] As a result 15 instances were conductedand each instance was conducted with 10 generations ofdifferent processing times The total is 150 instances The

094098102106

11

AlgorithmsAverageLower boundMaximum average

mak

espa

n m

eans

Aver

age o

f

LPT SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

Figure 3 Histogram of averages of makespan means

performance of the algorithms is measured with respectto the average makespan (mean) and average CPU time(Avg time in second) The ldquomeanrdquo performance is a relativequality measure of solutions computed by 119862LB where 119862is the average makespan obtained for each instance with 10generations given by the algorithm and LB is the lower boundof the instance calculated in equation that mentioned inSection 33 Step 2 The ldquoAvg timerdquo refers to the total time ittakes for the algorithm to finish the solution Table 2 presentsthe results of the previous algorithms from the literaturewhileTable 3 presents the results from the proposed algorithms Bycomparing the results of the average means of the makespanIt is obvious that the proposed algorithms outperform all thealgorithms in Table 2 from the literature It is worth notingthat for each instance the proposed algorithms obtain noworse than algorithms inTable 2 except at 10machines and 20jobs instance in only HIVNS algorithm and that is due to thedifficulty facing the proposed algorithms to get a lower boundwhen the difference between the number ofmachines and thenumber of jobs is relatively small In addition by comparingthe two proposed algorithms we can see that the versionshave the same average means of the makespan in case ofrandom initial solution while as the 2nd version outperformsthe 1st version in case of LPT initial solution in both averagemeans of makespan and average CPU time

Figure 3 shows the averages of makespan mean themaximum averages for each algorithm and the lower boundIt can be observed that the two proposed versions algorithmshave makespan means averages which closed to the lowerbound especially in RIVNS1 and RIVNS2 that have the sameaverage (10054)

Figure 4 shows the averages of Avg time means for allalgorithms We can see that RIVNS1 and RIVNS2 have thesmallest Avg time which are 0008 and 0007 respectivelyMoreover it can be observed that the Avg time of HIVNS1and HIVNS2 is much higher than HIVNS because in thispaper Matlab is used to construct the code of HIVNS1 and

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

Mathematical Problems in Engineering 5

Table 1 Number of machines and jobs

Number of machines 2 5 10 20Number of jobs 20 50 100 200 20 50 100 200 20 50 100 200 50 100 200

Step 2 Compute lower bound LB(119862max) = max[(1119898)sum119899119894=1 119901119894]max119894119901119894

Step 3 Repeat until 119891(119883) (makespan of 119883) is equal to LB or119894 gtMaxIterNum

Step 31 For schedule119883 distinguish the problemmachine set(119878pm) and the nonproblem machine set (119878npm)

Step 32 For each machine119872119901 in 119878pm do

Step 321 For each machine119872119899119901 in 119878npm do

Step 3211 Set 119870 = 1 finish = false

Step 3212

Repeat

Switch (119870) 119870 = 1 1198831015840 = Exchange 3(119872119901119872119899119901 119883)break119870 = 2 1198831015840 = Exchange 1(119872119901119872119899119901 119883)break119870 = 3 1198831015840 = Move(119872119901119872119899119901 119883) break119870 = 4 1198831015840 = Exchange 2(119872119901119872119899119901 119883)beak119870 = 5 1198831015840 = Exchange 4(119872119901119872119899119901 119883)if (119891(1198831015840) lt 119891(119883)) then set119883 = 1198831015840 finish =true go to Step 3 otherwise set 119896 = 119896 + 1

until 119896 = 119896max

Step 33 If finish = false then 119894 = 119894 + 1

Step 4 Output the best solution119883 found so far

4 Computational Results and Comparison

In this section the results of two versions of the proposedalgorithm were compared with LPT [1] SA [17] GA [14]MVNS [21] and IVNS [23] algorithms from the literatureThe two versions of the improved variable neighborhoodsearch algorithms ldquoHIVNS1 and RIVNS1rdquo and ldquoHIVNS2 andRIVNS2rdquo were coded in MATLAB R2012a and executed oni5 CPU 5GHz with 6GB of RAM All of them were stoppedafter getting the lower bound or running for 100 iterations forRIVNS1 and RIVNS2 The number of machines and numberof jobs are shown in Table 1

The processing time of the jobs is the same as Jingand Jun-qing [23] As a result 15 instances were conductedand each instance was conducted with 10 generations ofdifferent processing times The total is 150 instances The

094098102106

11

AlgorithmsAverageLower boundMaximum average

mak

espa

n m

eans

Aver

age o

f

LPT SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

Figure 3 Histogram of averages of makespan means

performance of the algorithms is measured with respectto the average makespan (mean) and average CPU time(Avg time in second) The ldquomeanrdquo performance is a relativequality measure of solutions computed by 119862LB where 119862is the average makespan obtained for each instance with 10generations given by the algorithm and LB is the lower boundof the instance calculated in equation that mentioned inSection 33 Step 2 The ldquoAvg timerdquo refers to the total time ittakes for the algorithm to finish the solution Table 2 presentsthe results of the previous algorithms from the literaturewhileTable 3 presents the results from the proposed algorithms Bycomparing the results of the average means of the makespanIt is obvious that the proposed algorithms outperform all thealgorithms in Table 2 from the literature It is worth notingthat for each instance the proposed algorithms obtain noworse than algorithms inTable 2 except at 10machines and 20jobs instance in only HIVNS algorithm and that is due to thedifficulty facing the proposed algorithms to get a lower boundwhen the difference between the number ofmachines and thenumber of jobs is relatively small In addition by comparingthe two proposed algorithms we can see that the versionshave the same average means of the makespan in case ofrandom initial solution while as the 2nd version outperformsthe 1st version in case of LPT initial solution in both averagemeans of makespan and average CPU time

Figure 3 shows the averages of makespan mean themaximum averages for each algorithm and the lower boundIt can be observed that the two proposed versions algorithmshave makespan means averages which closed to the lowerbound especially in RIVNS1 and RIVNS2 that have the sameaverage (10054)

Figure 4 shows the averages of Avg time means for allalgorithms We can see that RIVNS1 and RIVNS2 have thesmallest Avg time which are 0008 and 0007 respectivelyMoreover it can be observed that the Avg time of HIVNS1and HIVNS2 is much higher than HIVNS because in thispaper Matlab is used to construct the code of HIVNS1 and

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

6 Mathematical Problems in Engineering

Table2Th

emakespanresults

before

thec

hang

ingordero

fthe

neighb

orho

odstructures

[23]

Instance

number

119898119899

LPT

SAGA

MVNS

RIVNS

HIV

NS

Mean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

eMean

Avgtim

e1

2

2010

033

1000

610

149

1000

002215

1000

000106

1000

0000

0910

000

000

062

5010

001

1000

010

738

1000

005131

1000

000226

1000

000016

1000

0000

033

100

1000

010

000

19382

1000

027948

1000

0004

4510

000

00031

1000

0000

054

200

1000

010

000

09218

1000

023145

1000

000905

1000

000059

1000

0000

035

5

2010

315

10264

13614

10356

03412

10127

00124

10045

00083

10043

00073

650

10053

10045

25502

10312

06935

10050

00283

10012

00027

1000

900022

7100

10005

10005

42271

10242

13862

10052

00534

10002

000

4210

000

00015

8200

10003

10003

144302

10168

26560

10084

01064

10001

00075

1000

000035

9

10

2010

794

10792

19328

11336

03575

10987

00131

10799

004

0210

590

00384

1050

10207

10207

41428

11169

16303

10315

00315

10078

00075

10036

000

6711

100

10110

10110

2013

0111421

14547

10181

006

4710

031

000

6110

011

000

4312

200

10007

10007

10165

10611

33761

10127

01561

1000

400117

10002

00034

1320

5010

510

10510

24803

13167

09728

11671

004

4110

312

01389

10304

01392

14100

10123

10123

42799

12270

18111

10857

00905

10087

00179

1004

600087

15200

10063

10063

146684

11353

344

9210

415

01840

10039

00163

10024

00113

Average

10148

10142

507

7910

827

1598

210

319

006

3510

095

00182

1007

100152

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

Mathematical Problems in Engineering 7

Table 3 The makespan results after the effect of the changing order of the neighborhood structures

Instance number 119898 119899RIVNS1 HIVNS1 RIVNS2 HIVNS2

Mean Avg time Mean Avg time Mean Avg time Mean Avg time1

2

20 10000 00005 10000 00004 10000 00008 10000 000052 50 10000 00002 10000 00002 10000 00002 10000 000023 100 10000 00004 10000 00005 10000 00004 10000 000834 200 10000 00002 10000 00004 10000 00002 10000 070305

5

20 10005 00014 10000 00025 10005 00014 10000 000166 50 10000 00005 10000 00028 10000 00008 10000 001757 100 10000 00235 10000 01104 10000 00258 10000 017058 200 10000 00002 10000 25874 10000 00002 10000 250579

10

20 10625 00016 10682 00040 10625 00015 10602 0004010 50 10004 00103 10000 00156 10000 00064 10000 0058011 100 10000 00052 10000 01122 10000 00054 10000 0071312 200 10000 00238 10000 18051 10000 00207 10000 155613

2050 10175 00158 10256 00542 10175 00153 10263 00672

14 100 10003 00288 10000 02312 10003 00270 10000 0022215 200 10000 00004 10000 10538 10000 00004 10000 07414

Average 10054 00075 10063 03987 10054 00071 10058 03951

SA GA

MV

NS

RIV

NS

HIV

NS

RIV

NS1

HIV

NS1

RIV

NS2

HIV

NS2

CPU time 5078 1598 0064 0018 0015 0008 0399 0007 0395

0000

1000

2000

3000

4000

5000

6000

Tim

e (se

cond

)

Figure 4 Histogram of averages of Avg time means

HIVNS2 and the authors of [23] used C++ to constructHIVNS Thus the benefits of Avg time offered by C++ faroutweigh the simplicity of Matlab Avg time is especiallyimportant when dealing with algorithms sincemany calcula-tions involve immense optimization with complex equationsand algorithms or calculations with a large number ofiterations As the amount of data increases the computationtime (Avg time) for Matlab code increases significantlythereforeMatlab code takesmore time for those calculationsfor example in Table 3 the cases of 119898 = 2 5 10 20 119899 = 200are need to large number of iterations C++ is theway to go foralgorithms calculations because of its speed and versatility

5 Conclusion

In this paper two versions of improved variable neighbor-hood search (IVNS) algorithms are proposed for scheduling

identical parallelmachines IPMwith the objective of studyingthe effect of adding a newneighborhood structure and chang-ing the order of neighborhood structures on minimizing themakespan119862max In the proposed algorithms amachine basedencoding method and five neighborhood structures are usedto enhance the quality of the final solution Computationalresults showed that the proposed algorithms outperform allthe algorithms in the literature and obtain no worse thanalgorithms except when the number of machines and thenumber of jobs are relatively small which is due to thedifficulty facing the proposed algorithms to get a lower boundin that case In addition we concluded that the second versionoutperforms the first version in case of LPT initial solutionand therefore changing the order of neighborhood structureshas an effect on minimizing the makespan Further researchis to implement the proposed algorithms in scheduling ofunrelated parallel machines

Conflicts of Interest

The authors declare that there are no conflicts of interestregarding the publication of this paper

Acknowledgments

The authors extend their appreciation to the Deanship ofScientific Research at King Saud University for funding thiswork through Research Group no RG-1439-009

References

[1] M R Gary and D S Johnson Computers and IntractabilityA Guide to the Theory of NP-completeness WH Freeman andCompany New York NY USA 1979

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

8 Mathematical Problems in Engineering

[2] R L Graham ldquoBounds on multiprocessing timing anomaliesrdquoSIAM Journal onAppliedMathematics vol 17 pp 416ndash429 1969

[3] S L van de Velde ldquoDuality-based algorithms for schedulingunrelated parallel machinesrdquo ORSA Journal on Computing vol5 no 2 pp 192ndash205 1993

[4] E Mokotoff ldquoAn exact algorithm for the identical parallelmachine scheduling problemrdquo European Journal of OperationalResearch vol 152 no 3 pp 758ndash769 2004

[5] J Coffman M R Garey and D S Johnson ldquoAn application ofbin-packing to multiprocessor schedulingrdquo SIAM Journal onComputing vol 7 no 1 pp 1ndash17 1978

[6] J H Blackstone Jr andD T Phillips ldquoAn improved heuristic forminimizing makespan among m identical parallel processorsrdquoComputers amp Industrial Engineering vol 5 no 4 pp 279ndash2871981

[7] C-Y Lee and J D Massey ldquoMultiprocessor scheduling com-bining LPT andMULTIFITrdquoDiscrete Applied Mathematics vol20 no 3 pp 233ndash242 1988

[8] M Y Yue ldquoOn the exact upper bound for the multifit processorscheduling algorithmrdquo Annals of Operations Research vol 24no 1ndash4 pp 233ndash259 1990

[9] C-Y Lee and J D Massey ldquoMultiprocessor scheduling Anextension of the MULTIFIT algorithmrdquo Journal of Manufactur-ing Systems vol 7 no 1 pp 25ndash32 1988

[10] J C Ho and J S Wong ldquoMakespanminimization for m parallelidentical processorsrdquo Naval Research Logistics (NRL) vol 42no 6 pp 935ndash948 1995

[11] J Riera D Alcaide and J Sicilia ldquoApproximate algorithms forthe 119875 119862max problemrdquo TOP vol 4 no 2 pp 345ndash359 1996

[12] R Cheng and M Gen ldquoParallel machine scheduling problemsusing memetic algorithmsrdquo Computers amp Industrial Engineer-ing vol 33 no 3-4 pp 761ndash764 1997

[13] S F Ghomi and F J Ghazvini ldquoA pairwise interchange algo-rithm for parallel machine schedulingrdquo Production Planning ampControl vol 9 pp 685ndash689 1998

[14] L Min and W Cheng ldquoA genetic algorithm for minimizing themakespan in the case of scheduling identical parallel machinesrdquoArtificial Intelligence in Engineering vol 13 no 4 pp 399ndash4031999

[15] J N D Gupta and J Ruiz-Torres ldquoA LISTFIT heuristic for min-imizing makespan on identical parallel machinesrdquo ProductionPlanning and Control vol 12 no 1 pp 28ndash36 2001

[16] A M Costa P A Vargas F J Von Zuben and P M FrancaldquoMakespan minimization on parallel processors An immune-based approachrdquo in Proceedings of the 2002 Congress on Evolu-tionary Computation (CEC rsquo02) pp 920ndash925 IEEE May 2002

[17] W-C Lee C-C Wu and P Chen ldquoA simulated annealingapproach to makespan minimization on identical parallelmachinesrdquo The International Journal of Advanced Manufactur-ing Technology vol 31 no 3-4 pp 328ndash334 2006

[18] L Tang and J Luo ldquoA new ILS algorithm for parallel machinescheduling problemsrdquo Journal of Intelligent Manufacturing vol17 no 5 pp 609ndash619 2006

[19] D E Akyol andGM Bayhan ldquoMinimizingmakespan on iden-tical parallel machines using neural networksrdquo in Proceedings ofthe International Conference on Neural Information Processingvol 4234 of Lecture Notes in Computer Science pp 553ndash562Springer 2006

[20] A H Kashan and B Karimi ldquoA discrete particle swarmoptimization algorithm for scheduling parallelmachinesrdquoCom-puters amp Industrial Engineering vol 56 no 1 pp 216ndash223 2009

[21] M Sevkli and H Uysal ldquoA modified variable neighbor-hood search for minimizing the makespan onidentical parallelmachinesrdquo in Proceedings of the 2009 International Conferenceon Computers and Industrial Engineering (CIE rsquo09) pp 108ndash111IEEE July 2009

[22] J Chen Q-K Pan and H Li ldquoHarmony search algorithmwith dynamic subpopulations for scheduling identical parallelmachinesrdquo in Proceedings of the 2010 6th International Confer-ence on Natural Computation (ICNC rsquo10) pp 2369ndash2373 IEEEAugust 2010

[23] C Jing and L Jun-qing ldquoEfficient variable neighborhood searchfor identical parallel machines schedulingrdquo in Proceedings of theControl Conference (CCC rsquo12) pp 7228ndash7232 IEEE 2012

[24] M Sevkli and A Z Sevkli A Stochastically Perturbed ParticleSwarm Optimization for Identical Parallel Machine SchedulingProblems INTECH Open Access Publisher 2012

[25] D Laha ldquoA simulated annealing heuristic for minimizingmakespan in parallel machine schedulingrdquo in Proceedings of theInternational Conference on Swarm Evolutionary and MemeticComputing pp 198ndash205 Springer 2012

[26] N Mladenovic and P Hansen ldquoVariable neighborhood searchrdquoComputers amp Operations Research vol 24 no 11 pp 1097ndash11001997

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom

Hindawiwwwhindawicom Volume 2018

MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Mathematical Problems in Engineering

Applied MathematicsJournal of

Hindawiwwwhindawicom Volume 2018

Probability and StatisticsHindawiwwwhindawicom Volume 2018

Journal of

Hindawiwwwhindawicom Volume 2018

Mathematical PhysicsAdvances in

Complex AnalysisJournal of

Hindawiwwwhindawicom Volume 2018

OptimizationJournal of

Hindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom Volume 2018

Engineering Mathematics

International Journal of

Hindawiwwwhindawicom Volume 2018

Operations ResearchAdvances in

Journal of

Hindawiwwwhindawicom Volume 2018

Function SpacesAbstract and Applied AnalysisHindawiwwwhindawicom Volume 2018

International Journal of Mathematics and Mathematical Sciences

Hindawiwwwhindawicom Volume 2018

Hindawi Publishing Corporation httpwwwhindawicom Volume 2013Hindawiwwwhindawicom

The Scientific World Journal

Volume 2018

Hindawiwwwhindawicom Volume 2018Volume 2018

Numerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisNumerical AnalysisAdvances inAdvances in Discrete Dynamics in

Nature and SocietyHindawiwwwhindawicom Volume 2018

Hindawiwwwhindawicom

Dierential EquationsInternational Journal of

Volume 2018

Hindawiwwwhindawicom Volume 2018

Decision SciencesAdvances in

Hindawiwwwhindawicom Volume 2018

AnalysisInternational Journal of

Hindawiwwwhindawicom Volume 2018

Stochastic AnalysisInternational Journal of

Submit your manuscripts atwwwhindawicom