Top Banner
Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC Algorithm with Genetic and Local Search for solving QAP Suman Samanta * , Department of Industrial and Management Engineering Indian Institute of Technology, Kanpur Kanpur, India E-mail (*corresponding author): [email protected] Deepu Philip Department of Industrial and Management Engineering Indian Institute of Technology Kanpur-208016, India. J. Ramkumar Department of Mechanical Engineering and Design Program Indian Institute of Technology Kanpur-208016, India Abstract—The two most commonly used local search operators in heuristic or metaheuristic approach to solve combinatorial optimization problems are genetic search operators and neighborhood search operators. In this work both the operators are combined to develop two modified artificial bee colony algorithms (ABC) for solving combinatorial optimization problems. We used those modified ABCs to solve some standard quadratic assignment problems (QAP) available in QAP library. The performance of the algorithms are compared with some state of the art algorithms by means of percentage variation from the known optima and minimum time requirement to reach the optima. Keywords— Combinatorial optimization problem, quadratic assignment problem, artificial bee colony algorithm. I. INTRODUCTION QAP is to allocate n facilities in n locations, taken from n available locations in such way that minimizes the interdepartmental traveling distance, cost or time for any facility. It is one of the fundamental combinatorial optimization problems from the category of facility layout optimization. This problem was first formalized by Koopmans and Beckmann [1]. The problem is formulated as shown in Equation 1. () ( ) ∑∑ = = n i n k k i ik n b a S 1 1 min ϕ ϕ ϕ (1) Where n S is the set of all n permutations of the integers 1, 2, …..,n. Each individual value of the product () ( ) k i ik b a ϕ ϕ denotes the transportation distance caused by assigning facility i to location () i ϕ and facility k to location () k ϕ . QAP is an NP-hard problem by nature [2 – 4]. Because of this, heuristic and meta-heuristic approaches are widely used to solve this kind of problem. Two most commonly used local search operators these heuristics or metaheuristics commonly adopt are neighborhood search operator and genetic search operator. The neighborhood search operators try to improve the solution by recombination of locations within itself. Which means the child (new solution) is generated by exploiting only one parent. Some example of such operators are swap or exchange, insert, inverse, translocation, scramble etc. [5]. Among these search operators, ‘swap’ and ‘insert’ are most commonly used. Some of the algorithms that uses neighborhood search operator for local search are: tabu search (TS) [6, 7], robust tabu search (ROTS) [8], simulated annealing (SA) [9, 10], go with the winner (GWW) [11], variable neighborhood search (VNS) [12, 13] etc. The genetic search operator uses crossover between more than one (mostly two) initial solutions (parents) to form the offspring or child. Some example of such operators are: position based crossover (PBX), order based crossover (OBX), one point crossover (1PX), Cycle crossover operator (CX), Order crossover (OX), Linear order crossover operator (LOX), Partially mapped crossover operator (PMX), Edge recombination crossover operator (ER), Two-point crossover operator (2PX_V1, V2, and V3). Kellego¨z et. al [14] proved that PBX and OBX are better than other genetic operators by using fifteen one machine total weighted tardiness problem. Some 1207 © IEOM Society International
10

Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Mar 15, 2018

Download

Documents

dinhhanh
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

A Modified ABC Algorithm with Genetic and Local Search for solving QAP

Suman Samanta*, Department of Industrial and Management Engineering

Indian Institute of Technology, Kanpur Kanpur, India

E-mail (*corresponding author): [email protected]

Deepu Philip Department of Industrial and Management Engineering

Indian Institute of Technology Kanpur-208016, India.

J. Ramkumar Department of Mechanical Engineering and Design Program

Indian Institute of Technology Kanpur-208016, India

Abstract—The two most commonly used local search operators in heuristic or metaheuristic approach to solve combinatorial optimization problems are genetic search operators and neighborhood search operators. In this work both the operators are combined to develop two modified artificial bee colony algorithms (ABC) for solving combinatorial optimization problems. We used those modified ABCs to solve some standard quadratic assignment problems (QAP) available in QAP library. The performance of the algorithms are compared with some state of the art algorithms by means of percentage variation from the known optima and minimum time requirement to reach the optima.

Keywords— Combinatorial optimization problem, quadratic assignment problem, artificial bee colony algorithm.

I. INTRODUCTION

QAP is to allocate n facilities in n locations, taken from n available locations in such way that minimizes the interdepartmental traveling distance, cost or time for any facility. It is one of the fundamental combinatorial optimization problems from the category of facility layout optimization. This problem was first formalized by Koopmans and Beckmann [1]. The problem is formulated as shown in Equation 1.

( ) ( )∑∑= =∈

n

i

n

kkiik

n

baS 1 1

minϕϕϕ

(1)

Where nS is the set of all n permutations of the integers 1, 2, …..,n. Each individual value of the product ( ) ( )kiikba ϕϕ denotesthe transportation distance caused by assigning facility i to location ( )iϕ and facility k to location ( )kϕ .

QAP is an NP-hard problem by nature [2 – 4]. Because of this, heuristic and meta-heuristic approaches are widely used to solve this kind of problem. Two most commonly used local search operators these heuristics or metaheuristics commonly adopt are neighborhood search operator and genetic search operator. The neighborhood search operators try to improve the solution by recombination of locations within itself. Which means the child (new solution) is generated by exploiting only one parent. Some example of such operators are swap or exchange, insert, inverse, translocation, scramble etc. [5]. Among these search operators, ‘swap’ and ‘insert’ are most commonly used. Some of the algorithms that uses neighborhood search operator for local search are: tabu search (TS) [6, 7], robust tabu search (ROTS) [8], simulated annealing (SA) [9, 10], go with the winner (GWW) [11], variable neighborhood search (VNS) [12, 13] etc. The genetic search operator uses crossover between more than one (mostly two) initial solutions (parents) to form the offspring or child. Some example of such operators are: position based crossover (PBX), order based crossover (OBX), one point crossover (1PX), Cycle crossover operator (CX), Order crossover (OX), Linear order crossover operator (LOX), Partially mapped crossover operator (PMX), Edge recombination crossover operator (ER), Two-point crossover operator (2PX_V1, V2, and V3). Kellego¨z et. al [14] proved that PBX and OBX are better than other genetic operators by using fifteen one machine total weighted tardiness problem. Some

1207© IEOM Society International

Page 2: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

application of genetic searches for combinatorial optimization is: genetic algorithm (GA) [15], genetic search based artificial bee colony algorithm (g-ABC) [16], best-so-far ABC [17] etc. However, to solve medium to high difficulty level problems, there has an increasing trend of combining both the genetic and local search operators; eg. Genetic hybrid algorithms [18-24].

The two strategies that various neighborhood search algorithms adopt while neighborhood search are any neighbor search and best neighbor search. Any neighbor search does not compare among the neighbors and randomly select any neighbor. Some algorithms that use the any neighbor search are VNS, GWW, genetic local search (GLS) [21] etc. On the other hand, the best neighbor strategy considers the best possible neighbor under the given conditions. TS, ROTS etc. are the examples of such algorithms which explore the best possible neighbor.

In this work, the genetic and neighborhood search operators are used in the two basic local search steps (viz. employed step and onlooker step) correspondingly of artificial bee colony algorithm to solve QAP. Two kinds of neighborhood search strategies, viz. any and best, are utilized. The two resulted algorithms are named as m-ABC1 and m-ABC2 here after. The algorithms are used to solve various standard QAPs available in the QAP library [4]. The performance of the algorithms are compared with other state of the algorithms viz. ROTS [8], GWW [10], VNS –II [13], GA [15], g-ABC [16], best-so-far ABC (b-s-f ABC) [17], Genetic local search (GLS) [21] by two criteria: the optimum value found by the algorithm under similar parametric settings and minimum time to reach the optimum value. The result shows the developed algorithms performs significantly better than most of the other state of the art algorithms.

II. MODIFIED ABC WITH GENETIC AND NEIGHBORHOOD SEARCH OPERATOR FOR QAP Artificial bee colony algorithm (ABC) is a non-conventional optimization algorithm developed by Karaboga [25]. With its

two-step local search (employed and onlooker) and one step global search (scout) methodology, it became a very popular choice in solving various optimization problems; viz. power flow optimization and optimal sizing of photovoltaic systems [26, 27], control engineering [28], image processing [29], design optimization [30], parametric optimization non-traditional manufacturing processes[31] etc. In this work, we used the main essence of ABC algorithm while using a genetic operator in the employed step and neighborhood search operator in the onlooker step for solving QAP. We also have used two different search strategy in the neighborhood search (onlooker step) two form two variations of the modified ABC – m-ABC1, and m-ABC2. The various steps involved in both the variation are discussed below. The pseudocodes of the algorithm is described in Fig. 1-5.

A. Initialization This step selects random food source (feasible solution) within the search space for the formation of the initial population.

Each food source is assigned to an employed bee for further improvement.

B. Employed step In this step, each employed bee selects a neighborhood position randomly and performs a local search by using genetic

search operator PBX [13]. Greedy selection approach is applied to select better food source among the newly explored food sources and old food sources.

C. Onlooker step The onlooker selects the food sources based on their fitness. The formulae for calculation of the fitness and finding out the

probability of selection of a food source is similar to as described by Karaboga [24]. For m-ABC1, we use neighborhood search operator SWAP or INSERT [5], each with a probability of 0.5 for each selected food source to do the required local search. While for m-ABC2, we use the best neighbor selection approach instead of any neighbor. The nectar amount of various neighbors is identified by the way described by Taillard [8].

D. Scoutstep Similar to the ABC described by Karaboga [24], here also a food source is abandoned for further search when it does not

show any improvement for a predefined number of trials (limit). This kind of food source is replaced by new randomly selected food source. However, the insertion of a new food source to the population is restricted to maximum one in a cycle.

1208© IEOM Society International

Page 3: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

Fig. 1. Flow for m-ABC algorithms

Fig. 2.Pseudocode for send employed bees

Send Employed Bees (colony) For i= 1 to CS

For each food source xi Select a neighbor xk randomly from the colony (neighbor position ≠ i) do PBX between xi and xk to identify new neighboring solution x2i

end For do CalculateFitness (xi and x2i) to calculate fitness value of each food source do GreedySelection between xi and x2i

Algorithm: m-ABC for QAP o Input

Colony size (CS): the size of the initial population limit: the abandonment criterion MaxCycle: maximum cycle number runtime: number of runs

o Output: Best solution found

1. Initialization: Initialize the population

2. For r=1 to runtime while t <= MaxCycle send employed bees (colony) send onlooker bees (colony, fitness) send scout bee (colony, limit)

3. Return Best solution found

1209© IEOM Society International

Page 4: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

Fig. 3.Pseudocode for send onlooker bees for m-ABC1

III. SIMULATION RESULTS AND ANALYSIS

A. Benchmarks In this work, we have used 20 benchmark problems starting from n=12 to n=64 having difficulty level low to medium from

the QAP Library [4] for accessing the performance of each algorithm. The problems are: Chr12c, Had14, Nug15, Had16, Nug17, Had18, Els19, Scr20, Nug21, Nug22, Nug24, Nug25, Bur26h, Nug27, Nug28, Lipa30a, Esc32h, Lipa40b, Lipa50b, and Esc64a.

B. Performance criteria Two performance criteria are used to compare among the algorithms:

1) Percentage difference to the optima (dif): The percentage difference between the optima found by the algorithm and known optima [4].

2) The minimum time is taken by the algorithm to reach the optima (Mintime).

Send Onlooker Bees (colony, fitness)(m-ABC1) Selection probability of each food source i, probi is

∑=

= CS

ii

ii

fit

fitprob

1

i = 1, t =0; while t<CS

if rand<probi t=t+1; xt = xi; if rand < 0.5

do SWAP on xt to generate new food source xt'

else, do INSERT on xt to generate new food source xt

' end if

end if i=i+1 if i==CS+1 i=1 end if

do GreedySelection between xi and xt'

end while do CalculateFitness (xi)

1210© IEOM Society International

Page 5: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

Fig. 4. Pseudocode send onlooker bees for m-ABC2

Fig. 5.Pseudocode for send scout bees

C. Parameter settings Since we are calculating the minimum time requirement to reach the optima as a comparison criterion, it is important to

maintain equivalent settings for each algorithms and also to run on the same platform. To do that, the initial population size for all the population based algorithms are made similar. All the algorithm codes are reproduced in matlab, and run them on Matlab R2014a interface in an Intel i7 -4770s CPU with cloak speed 3.10 GHz and 8GB DDR3 ram system. All the parametric settings for different algorithms are summarized in TABLE I.1) Fig. 6 and 7 shows the convergence of various algorithms for Nug15 and Nug28 respectively.

Send scout bee (colony, limit) Find limitk==max(limt) If max(limit)>MaxLimit

limitk=0; xk= randomly selected new food source evaluate and include xk in the colony

end if

Send Onlooker Bees (colony, fitness)(m-ABC2) Selection probability of each food source i, probi is

∑=

= CS

ii

ii

fit

fitprob

1

i = 1, t =0, pi=0; while t<CS

if rand<probi t=t+1; pi= pi+1; xt = xi; xt

' = best neighbor of xt end if i=i+1 if i==CS+1 i=1 end if

do GreedySelection between xi and xt'

end while do CalculateFitness (xi)

1211© IEOM Society International

Page 6: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

TABLE I. PARAMETER SETTING FOR VARIOUS ALGORITHMS

Algorithm Colony size/population size

Limit/mutation rate/random walk (L)

Maximum cycle number

m-ABC1 2n for n<28 4n for n>= 28 Limit: 100 100n

g-ABC -do- -do- -do- b-s-f ABC -do- -do- -do-

GLS -do- -NA- -do- m-ABC2 2n Limit:100 -do-

GWW n for n <=20 2n for n>20

L=50 for n<=16 L=100 for n>16

Till the algorithm reaches to the local optima

ROTS -NA- -NA- 100n

GA 2n for n<30 4n for n>= 30 Mutation rate: 5% 100n

VNS-II -do- Shaking limit:100 -do-

D. Simulation Results for the Benchmarked Problems All the algorithms were run for 30 times with the above mentioned parametric settings for each problem instances.

TABLEII summarizes the run results. The red colored timings show the best found timings and green colored timings show the second best timing to reach the optima of a specific problem instance.

E. Analysis of the results Simulation results of the various algorithms under the stated parametric settings reveals various interesting insights.

1) With the lower value of n (<17) all the algorithms have converged to the optima, whereas, with the increase of n, all the algorithms could not reach to the optima.

2) Four algorithms have provided optimum results in all the benchmarked problems considered, viz. m-ABC1 and 2, GWW, and ROTS.

3) ROTS scored minimum time in total 16 problems, however, m-ABC2 scored the minimum time in four problems (Chrl12c, Nug15, Bur26h, Nug28).m-ABC1 could not score the fastest time in any of the instances.

4) The developed m-ABC2 gives 2nd fastest to reach optima after ROTS in all the other instances

5) The convergence curves in Fig. 6 and 7 depict that m-ABC2 and ROTS have very similar rate of convergence at the beginning, however, the m-ABC2 has stiffer rate of convergence afterwards.

6) m-ABC1 has less rate of convergence than ROTS and m-ABC2. However considering both the figure, it can easily be said that this algorithm has third best convergence rate.

7) GWW has very less convergence rate at the beginning, however, the rate changes drastically with increase in cycle number, which results in very stiff convergence curve.

8) GA, GLS, and VNS-II also provide optimum value in most of the cases and reaches very near to the optima in almost all the cases under the tested parameters (TABLE II). The convergence of these algorithms are also seen to be good (especially GA and VNS-II) and lies very near to m-ABC1 (Fig.6 and 7). However, the least time taken by these algorithms to reach the optima is higher in almost all the cases in comparison with m-ABC1 (exception: Had18).

9) TABLE II and Fig.6 and 7 also reveals that g-ABC and b-s-f ABC performed worse under the given conditions and problem instances.

1212© IEOM Society International

Page 7: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proc

eedi

ngs o

f the

201

6 In

tern

atio

nal C

onfe

renc

e on

Indu

stria

l Eng

inee

ring

and

Ope

ratio

ns M

anag

emen

t Ku

ala

Lum

pur,

Mal

aysia

, Mar

ch 8

-10,

201

6 TA

BLE

II.

SIM

ULA

TIO

N R

ESU

LTS

FOR

TH

E B

ENC

HM

AR

KED

PR

OB

LEM

S

Prob

lem

n

Opt

ima

m-A

BC

1 m

AB

C2

g-A

BC

b-

s-f A

BC

G

LS

GW

W

RO

TS

GA

V

NS

dif

Min

time

dif

Min

time

dif

Min

time

dif

Min

time

dif

Min

time

dif

Min

time

di)

Min

time

dif

Min

time

dif

Min

time

Chr

12c

12

11

156

0 0.

06

0 0.

05

0.00

0.

84

0.00

2.

20

0.00

0.

29

0 0.

36

0 0.

07

0.00

0.

23

0.00

0.

09

Had

14

14

2724

0

0.12

0

0.09

0.

00

0.64

0.

00

1.45

0.

00

0.65

0

0.61

0

0.03

0.

00

0.35

0.

00

0.15

Nug

15

15

1150

0

0.25

0

0.12

0.

00

2.65

0.

00

2.65

0.

00

1.21

0

1.25

0

0.13

0.

00

0.73

0.

00

0.60

Had

16

16

3720

0

0.27

0

0.16

0.

00

2.58

0.

00

2.44

0.

00

0.90

0

0.78

0

0.03

0.

00

0.50

0.

00

0.33

Nug

17

17

1732

0

0.45

0

0.40

1.

50

7.05

1.

27

3.94

0.

35

4.41

0

1.82

0

0.11

0.

00

3.14

0.

00

2.28

Had

18

18

5358

0

0.72

0

0.29

0.

04

6.11

0.

00

5.30

0.

00

2.97

0

1.96

0

0.16

0.

00

0.75

0.

00

0.48

Els

19

19

1721

2548

0

0.44

0

0.39

0.

00

1.76

0.

00

3.20

0.

00

1.36

0

2.82

0

0.16

0.

00

1.00

0.

00

1.81

Scr2

0 20

11

0030

0

1.04

0

0.33

0.

00

10.0

4 0.

00

2.86

0.

03

3.02

0

3.33

0

0.14

0.

76

3.16

0.

07

2.10

Nug

21

21

2438

0

1.55

0

0.70

1.

39

12.6

0 1.

72

11.3

0 0.

33

2.31

0

7.76

0

0.40

0.

41

6.67

0.

49

5.55

Nug

22

22

3596

0

1.53

0

0.66

2.

11

14.2

9 1.

95

8.21

0.

72

9.17

0

8.78

0

0.46

0.

00

4.00

0.

39

1.68

Nug

24

24

3488

0

1.18

0

1.33

3.

38

9.47

3.

44

16.6

7 0.

86

9.21

0

9.21

0

0.86

0.

63

7.09

1.

15

4.40

Nug

25

25

3744

0

1.37

0

1.33

8.

23

10.1

2 10

.04

15.4

2 0.

21

7.74

0

9.89

0

0.98

0.

05

7.61

0.

80

2.28

Bur

26h

26

7098

658

0 2.

78

0 1.

41

0.42

30

.25

0.40

24

.88

0.00

11

.11

0 11

.56

0 1.

99

0.01

17

.61

0.01

6.

55

Nug

27

27

5234

0

2.23

0

2.47

3.

32

17.7

5 3.

44

27.0

7 1.

72

12.1

9 0

11.7

8 0

1.48

0.

84

12.0

7 1.

83

7.07

Nug

28

28

5166

0

12.0

5 0

3.78

9.

52

74.6

7 7.

05

75.8

6 2.

71

7.97

0

13.1

1 0

7.98

0.

93

30.6

7 2.

28

12.0

6

Lip

a30a

30

13

178

0 8.

17

0 3.

37

2.77

86

.39

2.77

35

.26

1.62

24

.76

0 17

.46

0 2.

71

1.52

37

.11

1.68

16

.20

Esc

32h

32

438

0 2.

73

0 1.

41

8.22

63

.03

5.94

49

.98

0.00

25

.71

0 6.

72

0 1.

02

0.46

18

.66

0.00

23

.61

Lip

a40b

40

47

6581

0

24.4

8 0

8.33

23

.35

77.5

1 23

.15

194.

65

0.00

75

.82

0 44

.69

0 2.

53

0.00

35

.61

8.70

34

.82

Lip

a50b

50

12

1024

4 0

48.4

7 0

26.3

2 23

.69

196.

12

23.8

7 17

6.74

19

.52

136.

21

0 83

.40

0 3.

08

17.3

2 13

4.04

19

.0

29.9

4

Esc

64a

64

116

0 10

.01

0 7.

83

13.7

9 27

9.16

13

.79

130.

59

0.00

40

.21

0 17

.68

0 0.

54

0.00

20

.39

0.00

16

.57

1213© IEOM Society International

Page 8: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

Fig. 6. Convergence of various algorithms for Nug15

Fig. 7. Convergence of various algorithms for Nug28

IV. CONCLUSIONS This paper suggests two new variations of artificial bee colony algorithm named as m-ABC1 and m-ABC2 for solving

quadratic assignment problems with low to medium difficulty level. In this developed algorithm, two well-known local search approaches namely, genetic search and neighborhood search were used in two local search steps of ABC viz. employed step and onlooker step. The variation is made in the onlooker step by means of any neighborhood search and best neighborhood search, to form two variety of ABC.

The m-ABC1 and m-ABC2 were evaluated on 20 well-known test problems available in QAP library and also compared with other seven known algorithms for combinatorial optimization. The test result proves the superiority of the developed algorithms over other six algorithms than ROTS in terms of finding the optima in minimum time. However the convergence rate of m-ABC2 is found to be very similar and in some cases more stiff than ROTS.

In conclusion, it can be said that the results are very much encouraging for a further improvement of this algorithm for difficult QAPs. Although the algorithm is used here to solve QAP, it can easily be used for solving other combinatorial optimization problems with very less modification. How this algorithm works for solving other combinatorial optimization problems remain a future study of this work.

1100

1150

1200

1250

1300

1350

1400

1450

1500

0 100 200 300 400 500 600 700

m-ABC1

g-ABC

bsf-ABC

GLS

GWW

ROTS

GA

VNS-II

m-ABC2

5000

5200

5400

5600

5800

6000

6200

6400

6600

6800

0 200 400 600 800 1000 1200 1400

mABC1GAVNS-IIg-ABCbsf-ABCGLSGWWROTSmABC2

1214© IEOM Society International

Page 9: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

REFERENCES [1] T. C. Koopmans and M. Beckman, “Assignment problems and the location of economic activities,” Econometrica, vol. 25, pp. 53–

76, 1957. [2] Sahni, S. and Gonzalez, P., “P-Complete approximation problems”, J. ACM, vol. 23, pp. 555–565, 1976. [3] Burkard, R.E. and Stratmann, K.H., “Numerical investigations on quadratic assignment problems” Nav. Res. Logist. Q., vol. 25,

pp. 129–148, 1978. [4] Burkard, R.E., Karisch, S. and Rendl, F., “QAPLIB-a quadratic assignment problem library”. J. Global. Optim. , 10, 391–403,

1997. [5] Pitzer, E., Beham, A., Affenzeller, M., “Algorithm selection for QAP using FLA”, In Proc of 13th European Conference, EvoCOP

2013, 2013. [6] Glover, F. “Tabu Search-part I”, ORSA Journal on Computing vol.1, pp. 190-206, 1989. [7] Glover, F. “Tabu Search-Part II” ORSA Journal on Computing vol. 2, pp. 4-32, 1990. [8] Taillard, E., “Robust taboo search for quadratic assignment problem”, parallel computing, vol. 17, 443-455, 1991. [9] Kirkpatrick, S.; Gelatt Jr, C. D.; Vecchi, M. P., "Optimization by Simulated Annealing". Science vol. 220 (4598), pp 671–680,

1983. [10] Černý, V., "Thermodynamical approach to the traveling salesman problem: An efficient simulation algorithm". Journal of

Optimization Theory and Applications vol 45, pp 41–51, 1985. [11] Aldous, D., Vazirani, U., “Go with the winners algorithms”. In Proceedings of the 35th IEEE Symposium on Foundations of

Computer Science, pp. 492-501, 1994. [12] Mladenovic, N., and Hansen, P., "Variable neighborhood search", Computers and Operations Research vol. 24 (11), pp. 1097–

1100, 1997. [13] Sevkli, M., and Aydin, M.E., ‘‘Variable Neighborhood Search for Job Shop Scheduling Problems’’, Journal of software, vol. 1(2),

pp. 34-39, 2006. [14] Kellego¨z, T., Toklu, B., and Wilson, J., “Comparing efficiencies of genetic crossover operators for one machine total weighted

tardiness problem”, Applied Mathematics and Computation, vol. 199, pp. 590-598, (2008). [15] Anderson, E. J., Ferris, M. C., “Genetic Algorithms for Combinatorial Optimization: The Assemble Line Balancing Problem”,

ORSA Journal on Computing vol. 6(2), pp161-173, 1994. [16] Samanta, S., Onkar, Philip, D., “Facility Layout Optimization of an Indian SME Using Artificial Bee Colony Algorithm”, In proc

of Int'l Conference on Artificial Intelligence, Energy and Manufacturing Engineering (ICAEME’2014), pp 16- 22, 2014. [17] Banharnsakun, A., Sirinaovakul, B., and Achalakul, T., “Job Shop Scheduling with the Best-so-far ABC”, Engineering

applications of artificial intelligence, vol. 25, pp583-593, 2012. [18] Fleurent, C., and Ferland, J. A., “Genetic hybrids for the quadratic assignment problem,” In p. Pardalos and h. Wolkowicz,

editors, quadratic assignment and related problems, Dimacs series in discrete mathematics and theoretical computer science, vol. 16, pp.173-187. 1994.

[19] Ostrowski, T., and Ruoppila, V.T., “Genetic annealing search for index assignment in vector quantization”, Technical report, pattern recognition letters, vol.18 (4), pp.311-318, 1997.

[20] Taillard, E.D., and Gambardella, L.M., “Adaptive memories for the quadratic assignment problem”, Research report, Idsia Lugano, Switzerland, 1997.

[21] Jaszkiewicz, A., “Genetic local search for multi objective combinatorial optimization”, European journal of operations research, vol. 137, pp. 50-71, 2002.

[22] Tseng, L., Liang, S., “A hybrid metaheuristic for the quadratic assignment problem”, Comput. Optim. Appl. vol.34 (1), pp85–113, 2005.

[23] Tosun, U., Dokeroglu, T., Cosar, A., “A new robust island parallel genetic algorithm for the quadratic assignment problem”. Int. J. Prod. Res. vol. 51 (14), pp4117–4133, 2013.

[24] Tosun, U., “On the performance of parallel hybrid algorithms for the solution of the quadratic assignment problem”, Engineering applications of artificial intelligence, vol. 39, pp. 267-278, 2015.

[25] D Karaboga., “An Idea Based On Honey Bee Swarm For Numerical Optimization”, Technical Report-TR06, Erciyes University, Computer Engineering Department, 2005.

[26] Kılıç, U., and Ayan, K., “Optimizing power flow of AC–DC power systems using artificial bee colony algorithm”, Int. J. Electr. Power Energy Syst., vol. 53, pp. 592–602, 2013

[27] Mohamed, A.F., Elarini, M.M., Othman, A.M., “A new technique based on artificial bee colony algorithm for optimal sizing of stand-alone photovoltaic system”, J. Adv. Res. vol. 5, pp. 397–408, 2014.

[28] ] Sahed, O.A., Kara, K., Benyoucef, A., “Artificial bee colony-based predictive control for non-linear systems”, Trans. Instit. Meas. Control, 2014.

1215© IEOM Society International

Page 10: Kuala Lumpur, Malaysia, March 8-10, 2016 A Modified ABC ...ieomsociety.org/ieom_2016/pdfs/333.pdf · combined to develop two modified artificial bee colony algorithms ... [6, 7],

Proceedings of the 2016 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, March 8-10, 2016

[29] Cuevas, E., Zaldívar, D., Pérez-Cisneros, M., Sossa, H., Osuna, V., “Block matching algorithm for motion estimation based onartificial bee colony (ABC)”, Appl. Soft Comput., vol. 13, pp. 3047–3059, 2013.

[30] Encan, A. S., ahin, S., Kılıc, B., Kılıc, U., “Design and economic optimization of shell and tube heat exchangers using artificialbee colony (ABC) algorithm”, Energy Convers. Manag. vol. 52, pp. 3356–3362, 2011.

[31] Samanta, S., and Chakraborty, S., “Parametric optimization of some non-traditional machining processes using artificial beecolony algorithm”, Eng. Appl. Artif. Intell. Vol.24 pp.946–957, 2011.

BIOGRAPHY

Suman Samantais Doctoral student in the department of Industrial and Management Engineering, Indian Institute of Technology, Kanpur. He completed is masters of production engineering (with university medal) in the year 2010 from the department of Production Engineering, Jadavpur University, Kolkata. He completed his graduation B.Tech (Hons) in Production Engineering and Management in the year 2007 from the National Institute of Technology, Jamshedpur. His research interests includes process simulation modeling, optimization, operations management etc. He is presently working in the areas of facility layout optimization, process simulation modeling and optimization, productivity improvement of various SMEs nearer to Kanpur, India.

Deepu Philip is currently an Associate Professor and teaches Production and Operations Management, Systems Simulation, and Automated Analytics. He has a bachelors in Industrial Engineering with his doctorate in Industrial & Management Engineering from MSU Bozeman. He has an outstanding combination of academic and industrial experience with leading organizations of the world. His research and teaching focus is on innovation solutions using cutting-edge digital technologies that employ the highest standards of application development. He has experience in designing and implementing complex system of systems in different fields including defense, aviation, transportation, banking, automation, health care, energy, communication, and decision making. In addition, he works on planning and scheduling, local search algorithms, systems simulation, decision support systems, distributed systems, and Fuzzy databases.

J.Ramkumar is currently a Professor and teaches manufacturing science, micro/Nano technology, and new product development. He has a bachelors inProduction Engineering with his doctorate in Defect quantification in drilling of composites from IIT Madras, India with a best thesis award. Over the yearshis contribution in teaching and research is remarkable. He has worked for BOSCH group and improved the productivity of the company. His research andteaching focus is on Nano technology and inclusive design. He has several international and national patents in his credit and has published more than 100journal papers.

1216© IEOM Society International