Alg lgebraic Crossover Operators for Permutations Valentino - - PowerPoint PPT Presentation

β–Ά
alg lgebraic crossover operators
SMART_READER_LITE
LIVE PREVIEW

Alg lgebraic Crossover Operators for Permutations Valentino - - PowerPoint PPT Presentation

Alg lgebraic Crossover Operators for Permutations Valentino Santucci 1,2 , Marco Baioletti 2 , Alfredo Milani 2 1 University for Foreigners of Perugia, Italy 2 University of Perugia, Italy Finitely Generated Group A group is a set


slide-1
SLIDE 1

Alg lgebraic Crossover Operators for Permutations

Valentino Santucci 1,2, Marco Baioletti 2, Alfredo Milani 2

1 University for Foreigners of Perugia, Italy 2 University of Perugia, Italy

slide-2
SLIDE 2

Finitely Generated Group

  • A group is a set π‘Œ endowed with an operation ∘: π‘Œ Γ— π‘Œ β†’ π‘Œ that:
  • is associative

𝑦 ∘ 𝑧 ∘ 𝑨 = 𝑦 ∘ 𝑧 ∘ 𝑨 βˆ€π‘¦, 𝑧, 𝑨 ∈ π‘Œ

  • has a neutral element

𝑓 ∈ π‘Œ s.t. 𝑦 ∘ 𝑓 = 𝑓 ∘ 𝑦 = 𝑦 βˆ€π‘¦ ∈ π‘Œ

  • has inverse elements

βˆƒ π‘¦βˆ’1 s.t. 𝑦 ∘ π‘¦βˆ’1 = π‘¦βˆ’1 ∘ 𝑦 = 𝑓 βˆ€π‘¦ ∈ π‘Œ

  • A group π‘Œ is finitely generated if there exists a finite generating set 𝐼 βŠ† π‘Œ such that

every 𝑦 ∈ π‘Œ can be expressed as a finite composition of the generators in 𝐼, i.e.,

𝑦 = β„Žπ‘—1 ∘ β„Žπ‘—2 ∘ β‹― ∘ β„Žπ‘—π‘™ with β„Žβˆ— ∈ 𝐼

  • The length of a minimal decomposition of 𝑦 in terms of H is the weight of 𝑦, we

denote it with |𝑦|

  • 𝑦 βŠ‘ 𝑧 iff there exists (at least) a minimal decomposition of x that is a prefix of a

minimal decomposition of y

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 2

slide-3
SLIDE 3

Permutations form a group

  • A permutation of [π‘œ] = {1,2, … , π‘œ} is a bijective discrete function from [π‘œ]

to [π‘œ], thus it is possible to compose permutations: 𝑨 = 𝑦 ∘ 𝑧 iff 𝑨 𝑗 = 𝑦(𝑧 𝑗 ) for 1 ≀ 𝑗 ≀ π‘œ

  • The composition:
  • is associative

𝑦 ∘ 𝑧 ∘ 𝑨 = 𝑦 ∘ (𝑧 ∘ 𝑨)

  • has neutral element

𝑓 = 1,2, … , π‘œ

  • has inverse elements

π‘¦βˆ’1 𝑗 = π‘˜ iff 𝑦 π‘˜ = 𝑗

  • The permutations of [π‘œ], together with the ∘ operation, form a group

structure called the symmetric group 𝒯(π‘œ)

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 3

slide-4
SLIDE 4

Permutations form a F.G. group

  • Adjacent swap moves as generating set

𝐡𝑇𝑋 = 𝜏1, 𝜏2, … , πœπ‘œβˆ’1 where, for any 1 ≀ 𝑗 < π‘œ: πœπ‘— π‘˜ = ࡞ 𝑗 + 1 if π‘˜ = 𝑗 𝑗 if π‘˜ = 𝑗 + 1 π‘˜

  • therwise
  • (𝑦 ∘ πœπ‘—) is the permutation 𝑦 where the items at positions 𝑗 and 𝑗 + 1

have been swapped

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 4

slide-5
SLIDE 5

Cayley graph

  • A finitely generated group induces a colored digraph (namely, a Cayley

graph) where:

  • π‘Œ are the vertices
  • there exists an arc 𝑦 β†’ (𝑦 ∘ β„Ž) colored by β„Ž for every 𝑦 ∈ π‘Œ and β„Ž ∈ 𝐼 βŠ† π‘Œ
  • Connection between a Cayley graph and a combinatorial search space:
  • π‘Œ is the set of solutions
  • 𝐼 βŠ† π‘Œ represents simple search moves in the space of solutions
  • The Cayley graph induces:
  • neighborhood relationships among solutions
  • a distance between solutions (shortest path distance)
  • a single representation for both solutions and displacements

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 5

slide-6
SLIDE 6

Cayley graph

Generators: <2134> <1324> <1243>

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 6

slide-7
SLIDE 7

Vector-like operations

  • Paths on the Cayley graph can be encoded by composing the labels on

the arcs => paths can be encoded using permutations!

  • Permutations encode both Β«pointsΒ» and Β«vectorsΒ» in the search

space

  • Let’s define βŠ•, βŠ–, βŠ™ in such a way that they work consistently w.r.t.

their usual numerical counterparts!

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 7

slide-8
SLIDE 8

Sum and Difference of Permutations

  • Given two n-length permutations x and y:
  • 𝑦 βŠ• 𝑧 ≔ 𝑦 ∘ 𝑧
  • 𝑦 βŠ– 𝑧 ≔ π‘§βˆ’1 ∘ 𝑦
  • They are geometrically meaningful:
  • 𝑦 βŠ• 𝑧 is the point reachable starting from point x and following the path y
  • 𝑦 βŠ– 𝑧 represents a path connecting the point y to the point x
  • They are algebraically consistent:

𝑦 = 𝑧 βŠ• 𝑦 βŠ– 𝑧 = 𝑧 ∘ π‘§βˆ’1 ∘ 𝑦 = 𝑦

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 8

slide-9
SLIDE 9

Multiply a permutation by a scalar in [0,1]

  • Given a n-length permutation z and a scalar 𝑏 ∈ 0,1
  • Let’s assume:
  • z represents a path
  • 𝑨 = πœπ‘—1 ∘ πœπ‘—2 ∘ β‹― ∘ πœπ‘—π‘€ is a minimal decomposition of 𝑨 with length 𝑀
  • 𝑏 βŠ™ 𝑨 ≔ πœπ‘—1 ∘ πœπ‘—2 ∘ β‹― ∘ πœπ‘—π‘™ where 𝑙 = 𝑏 βˆ™ 𝑀
  • Geometrically, 𝑏 βŠ™ 𝑨 is a Β«truncationΒ» of the path z

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 9

slide-10
SLIDE 10

How to compute them?

  • βŠ• and βŠ– simply require permutations composition and inversion
  • βŠ™ requires an algorithm for computing minimal decomposition(s)
  • There can be multiple minimal decompositions
  • They can be computed using a Β«bubble sortΒ»-like algorithm:

Iteratively choose (and apply) an adjacent swap moving the permutation closer to the identity permutation (the only sorted permutation)

  • Two different strategies:
  • RandBS: randomly choose one suitable adjacent swap
  • GreedyBS: choose the best (and suitable) adjacent swap basing on the fitness

function

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 10

slide-11
SLIDE 11

RandBS and GreedyBS

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 11

slide-12
SLIDE 12

What can we do with them?

  • Algebraic Differential Evolution
  • Algebraic Particle Swarm Optimization
  • Discretize numerical EAs whose Β«move rulesΒ» are linear combinations
  • f solutions
  • … design an algebraic crossover

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 12

slide-13
SLIDE 13

Algebraic Crossovers

  • We propose 3 classes of algebraic crossovers:
  • Group-based algebraic crossovers
  • Lattice-based algebraic crossovers
  • Hybrid algebraic crossovers

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 13

slide-14
SLIDE 14

Group-based Algebraic Crossovers (AXG)

  • Reasonably, a crossover between two permutations x and y should return a

permutation z that is Β«in the middleΒ» between x and y

  • The interval [x,y] can be formally defined in many equivalent ways:
  • A group-based algebraic crossover operator AXG can be abstractly defined

as any operator which, given two permutations x and y, returns a permutation z = AXG(x,y) such that 𝑨 ∈ [𝑦, 𝑧]

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 14

slide-15
SLIDE 15

A property on pairwise precedences of items

  • z = AXG(x,y)
  • AXG is precedence-respectful:

z contains all the common precedences between x and y

  • AXG transmits precedences:

all the common precedences of z come from x or y (no new precedence is generated)

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 15

slide-16
SLIDE 16

How to compute AXG?

  • Ideally:
  • 𝑨 = 𝑦 βŠ• 𝑏 βŠ™ (𝑧 βŠ– 𝑦), with 𝑏 ∈ [0,1]
  • Different strategies basing on the value of a and the βŠ™ computation strategy
  • Practically:
  • We enumerate all the permutations in a shortest path from x to y by running

RandBS (R) or GreedyBS (G) on 𝑧 βŠ– 𝑦

  • We select one permutation in the path: a random one (R), the middle one (T),

the best one (B)

  • 6 different implementations:
  • AXG-RR, AXG-RT, AXG-RB
  • AXG-GR, AXG-GT, AXG-GB

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 16

slide-17
SLIDE 17
  • Meet: 𝑨 = 𝑦 ∧ 𝑧 iff
  • 1. 𝑨 βŠ‘ 𝑦 and 𝑨 βŠ‘ 𝑧
  • 2. z is the Β«longestΒ» permutation

with prop.1

  • Join: 𝑨 = 𝑦 ∨ 𝑧 iff
  • 1. 𝑦 βŠ‘ 𝑨 and 𝑧 βŠ‘ 𝑨
  • 2. z is the Β«shortestΒ» permutation

with prop.1

The partial order βŠ‘ is a Lattice

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 17

slide-18
SLIDE 18

Meet and Join as Crossover Operators

AXL-Join exploits the Β«De MorganΒ»-like property:

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 18

slide-19
SLIDE 19

Hybrid Algebraic Crossovers

  • Let AXG be any group-based algebraic crossover
  • Let m = AXL-Meet(x,y)
  • Let j = AXL-Join(x,y)
  • An hybrid algebraic crossover AXH is defined as

AXH(x,y) := AXG(m,j)

  • 6 hybrid alg. crossovers: one for each AXG crossover

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 19

slide-20
SLIDE 20

AXH-* are more explorative than AXG-*

  • AXH crossovers produce an offspring 𝑨 ∈ [𝑦 ∨ 𝑧, 𝑦 ∧ 𝑧]
  • [𝑦, 𝑧] βŠ† [𝑦 ∨ 𝑧, 𝑦 ∧ 𝑧]
  • It is possible to introduce precedences not appearing in the parents

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 20

slide-21
SLIDE 21

Experimental Setting

  • Benchmark instances from LOP, PFSP, QAP, TSP
  • Comparison against 7 popular permutation crossover from literature:

PMX, OX1, OX2, CX, AP, POS, ER

  • Three scenarios:
  • Randomly generated permutations
  • Local optima permutations
  • Crossovers embedded in standard steady-state GA

(population size = 50, random selection, crossover prob = 1, mutation prob = 0.05, +1 replacement)

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 21

slide-22
SLIDE 22

Experimental Results

AXG crossovers loss diversity very quickly when inside a GA AXH crossovers slow down the diversity loss but only slightly

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 22

slide-23
SLIDE 23

Conclusion and Future Work

  • What’s new:

Group and lattice structures of the search space exploited to build crossovers and derive some properties of them

  • Issues to address:

Need to better balance intensification and diversification

  • Other future work:
  • AXG operators using different generating sets (EXC and INS)
  • AXL operators using semi-lattice
  • Build a clever GA using AX operators
  • Precedence properties may help to design a good recombination for LOP

11 Jul 2018 - WCCI 2018 - Rio de Janeiro (Brasil) Algebraic Crossover Operators for Permutations - V. Santucci, M. Baioletti, A. Milani 23

slide-24
SLIDE 24

Thanks!!! !!!