Oscar Cubo M edina fi.upm.es)

Similar documents
Metaheuristics and Local Search

Metaheuristics and Local Search. Discrete optimization problems. Solution approaches

Methods for finding optimal configurations

Nov Julien Michel

Algorithms and Complexity theory

Methods for finding optimal configurations

Hill climbing: Simulated annealing and Tabu search

Artificial Intelligence Heuristic Search Methods

3.4 Relaxations and bounds

GENETIC ALGORlIHMS WITH AN APPLICATION TO NONLINEAR TRANSPORTATION PROBLEMS

Zebo Peng Embedded Systems Laboratory IDA, Linköping University

Finding optimal configurations ( combinatorial optimization)

Chapter - 3. ANN Approach for Efficient Computation of Logarithm and Antilogarithm of Decimal Numbers

Algorithms. NP -Complete Problems. Dong Kyue Kim Hanyang University

Algorithm Design Strategies V

Totally unimodular matrices. Introduction to integer programming III: Network Flow, Interval Scheduling, and Vehicle Routing Problems

CS 331: Artificial Intelligence Local Search 1. Tough real-world problems

Basic Semantics of Prolog. CAS 706 Program m ing Language. Xinjun Wu mcmaster.ca

Local Search. Shin Yoo CS492D, Fall 2015, School of Computing, KAIST

Single Solution-based Metaheuristics

Random Search. Shin Yoo CS454, Autumn 2017, School of Computing, KAIST

C o r p o r a t e l i f e i n A n c i e n t I n d i a e x p r e s s e d i t s e l f

CS325: Analysis of Algorithms, Fall Final Exam

CS 6901 (Applied Algorithms) Lecture 3

More Approximation Algorithms

PROBLEM SOLVING AND SEARCH IN ARTIFICIAL INTELLIGENCE

Introduction to integer programming III:

APPLICATION OF AUTOMATION IN THE STUDY AND PREDICTION OF TIDES AT THE FRENCH NAVAL HYDROGRAPHIC SERVICE

Swarm Intelligence Traveling Salesman Problem and Ant System

CS/COE

Lecture 4: Simulated Annealing. An Introduction to Meta-Heuristics, Produced by Qiangfu Zhao (Since 2012), All rights reserved

Chapter 5 Workshop on Fitting of Linear Data

ECS122A Handout on NP-Completeness March 12, 2018

Lecture H2. Heuristic Methods: Iterated Local Search, Simulated Annealing and Tabu Search. Saeed Bastani

Large chunks. voids. Use of Shale in Highway Embankments

CS 6783 (Applied Algorithms) Lecture 3

Dynamic Programming. Weighted Interval Scheduling. Algorithmic Paradigms. Dynamic Programming

Vehicle Routing and Scheduling. Martin Savelsbergh The Logistics Institute Georgia Institute of Technology

P C max. NP-complete from partition. Example j p j What is the makespan on 2 machines? 3 machines? 4 machines?

Multiple Alignment. Slides revised and adapted to Bioinformática IST Ana Teresa Freitas

P, NP, NP-Complete. Ruth Anderson

Scheduling with AND/OR Precedence Constraints

There are self-avoiding walks of steps on Z 3

TIDAL PREDICTION WITH A SMALL PERSONAL COMPUTER

Optimization. Announcements. Last Time. Searching: So Far. Complete Searching. Optimization

Approximation Algorithms (Load Balancing)

Designing the Human Machine Interface of Innovative Emergency Handling Systems in Cars

CSE 421 Introduction to Algorithms Final Exam Winter 2005

Principles of AI Planning

Bio nformatics. Lecture 3. Saad Mneimneh

Embedded Systems 14. Overview of embedded systems design

A tabu search algorithm for the minmax regret minimum spanning tree problem with interval data

Recoverable Robustness in Scheduling Problems

NOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or

Spin Glas Dynamics and Stochastic Optimization Schemes. Karl Heinz Hoffmann TU Chemnitz

Principles of AI Planning

Artificial Intelligence Methods (G5BAIM) - Examination

Computational Complexity

Bio nformatics. Lecture 23. Saad Mneimneh

Local Search & Optimization

Dynamic Scheduling with Genetic Programming

Fundamentals of Metaheuristics

VIII. NP-completeness

Check off these skills when you feel that you have mastered them. Write in your own words the definition of a Hamiltonian circuit.

Improved Priority Assignment for the Abort-and-Restart (AR) Model

An Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse

Una descrizione della Teoria della Complessità, e più in generale delle classi NP-complete, possono essere trovate in:

Early-Tardy Minimization for Joint Scheduling of Jobs and Maintenance Operations on a Single Machine

CS Analysis of Recursive Algorithms and Brute Force

Protein folding. α-helix. Lecture 21. An α-helix is a simple helix having on average 10 residues (3 turns of the helix)

Multi-scale MetUM simulations ofa M esoscaleconvective System Sam antha Clarke

Introduction to Computer Science and Programming for Astronomers

NP-Complete Problems and Approximation Algorithms

Part III: Traveling salesman problems

Logic, Optimization and Data Analytics

A pruning pattern list approach to the permutation flowshop scheduling problem

CSE 431/531: Analysis of Algorithms. Dynamic Programming. Lecturer: Shi Li. Department of Computer Science and Engineering University at Buffalo

CS 580: Algorithm Design and Analysis. Jeremiah Blocki Purdue University Spring 2018

Constraint satisfaction search. Combinatorial optimization search.

Doctoral Course in Speech Recognition. May 2007 Kjell Elenius

Analysis of Uninformed Search Methods

Local Search & Optimization

Dynamic Programming. Reading: CLRS Chapter 15 & Section CSE 6331: Algorithms Steve Lai

A SEQUENTIAL ELIMINATION ALGORITHM FOR COMPUTING BOUNDS ON THE CLIQUE NUMBER OF A GRAPH

Part III: Traveling salesman problems

Scheduling with Common Due Date, Earliness and Tardiness Penalties for Multi-Machine Problems: A Survey

(tree searching technique) (Boolean formulas) satisfying assignment: (X 1, X 2 )

GRASP A speedy introduction

Unit 1A: Computational Complexity

University of Toronto Department of Electrical and Computer Engineering. Final Examination. ECE 345 Algorithms and Data Structures Fall 2016

NP-Completeness. Andreas Klappenecker. [based on slides by Prof. Welch]

Artificial Intelligence Methods (G5BAIM) - Examination

0 1 d 010. h 0111 i g 01100

This means that we can assume each list ) is

Introduction into Vehicle Routing Problems and other basic mixed-integer problems

Greedy Homework Problems

ARTIFICIAL INTELLIGENCE

A Polynomial-Time Algorithm to Find Shortest Paths with Recourse

Bi-criteria Scheduling Problems on Parallel Machines

A parallel metaheuristics for the single machine total weighted tardiness problem with sequence-dependent setup times

Transcription:

. euristics. reedy. LocalSearch. illimbing. Meta-heuristics. Tabu Search. Sim ulated nnealing Oscar ubo M edina (ocubo@ fi.upm.es) ased on Neighbourhood Search euristic derives from the verb heuriskein (ενρισκειν)which m eans to find aster than mathem aticaloptimization (branch & bound,simplex,etc) M eta m eans beyond,in an upper lever Meta-heuristics are strategies that guide the search process The goalis to explore the search space in order to find (near-)optimalsolutions Meta-heuristics are notproblem -specific The basic concepts ofm eta-heuristics perm itan abstract leveldescription They may incorporate mechanism s to avoid getting trapped in confined areas ofthe search space

euristic onstructive algorithm s (greedy) Localsearch algorithm s (hil-clim bing ) Meta-heuristic Trajectory methods: escribe a trajectory in the search space during the search process Variable Neighbourhood Search Iterated LocalSearch Sim ulated nnealing Tabu Search Population-based: Perform search processes which describe the evolution ofa setofpoints in the search space volutionary com putation enerate solutions from scratch by adding (to an initia ly em pty partialsolution) com ponents,untilthe solution is com plete greedy algorithm works w in phases.teach phase: You take the bestyou can getrightnow,w ithout regard for future consequences You hope thatby choosing a localoptimum at each step,you wilend up ata globaloptim um Scheduling: processors 9 jobs (,,,0,,,, and minutes) reedy assign to each free processor the longest-running job Scheduling: processors 9 jobs (,,,0,,,, and minutes) reedy assign to each free processor the shortest-running job 0 Total time: m 0 Total time: 0 m

Longest-running job: minutes Shortest-running running job: 0 minutes 0 Optimal: minutes 0 0 Iterative algorithm : Startfrom som e initialsolution xplore the neighbourhood ofthe current solution Replace the currentsolution by a better solution Neighbourhood: being X the search space,we define the neighbourhood system N in X as N : x N (x ) ifferentprocedures depending on choice criteria and term ination criteria Stochastic: choose atrandom ilclimbing: only perm itmoves to neighbours thatimprove the current reedy the bestneighbour nxious the first neighbour improving Sidew ays moves alow s moves with sam e fitness s= enerateinitialsolution() while Term ination criteria notm et s =PicktRandom ( N(x) ) If ( f(s ) < f(s) ) then s= s ndif end while

itness itness: : itness itness: : itness itness: :

itness: : asic principles: Keep only a single state (solution) in m em ory enerate only the neighbours ofthatstate Keep one ofthe neighbours and discard others Key features: No system atic No increm ental Key advantages: Use very little memory (constantam ount) ind solutions in search spaces too large for system atic algorithm s Success ofhil-clim bing depends on shape of landscape Shape oflandscape depends on problem form ulation and fitness function Landscapes for realistic problem s often look like a worst-case scenario NP-hard problem s typica ly have exponential num ber oflocal-m inima ailing on neighbourhood search Propensity to deliver solutions which are only localoptima Solutions depend on the initialsolution pplies localsearch to an initialsolution until itfinds the localoptimum; then itperturbs the solution and itrestartslocalsearch

Strategy based on dynam ica ly changing neighborhood structure Steepestdescentw ith m em ory Moves through solution space Uses memory techniques to avoid cycling The overa lapproach is to avoid entrainm entin cycles by forbidding or penalizing moves which take the solution,in the nextiteration,to points in the solution space previously visited (red lover,om puters and Operations Research,98) To avoid revelsalm oves the lastm oves are makedas Tabu hange value ofx from false totrue: [x = false] istabu Sw ap elem entsiand j: [hange j,i] istabu rop iand add j: [add i] and [drop j] are tabu MinimalSpanning Tree: Totalcostofthe links used is a minimum lthe points are connected together onstraints: x + x + x x x Penalty: 0 x 9 x 0 x x x x x 7 8

The move considered is the standard edge sw ap move m n edge is listed as tabu ifitwas added within the lasttwo iterations ifitwas deleted within the lasttwo iterations Iteration x + x + x x x Penalty: 0 rop x 9 x 0 x x x x x 7 8 (X) = +00 dd Iteration x + x + x x x Penalty: 0 Tabu x 9 x 0 x (X) = 8 Tabu Iteration x + x + x x x Penalty: 0 Tabu dd x 9 x x 0 (X) = Tabu rop x x x x rop x x dd 7 8 x x 7 Tabu 8 Tabu

Iteration x + x + x x x Penalty: 0 Tabu dd x 9 x x 0 (X) = Tabu rop Iteration x + x + x x x Penalty: 0 Tabu x 9 x x 0 (X) = Tabu x x x x x x 7 Tabu 8 Tabu x x 7 Tabu 8 Tabu Pros: Tabu Search yields relatively good solutions to previously intractable problem s Tabu Search provides com parable or superior solutions to other optim ization techniques ons: Tabu Search does notguarantee optim ality Tabu Search is aw kw ard for problem s with continuous variables Tabu Search assum es fastperform ance evaluation The construction oftabu listis heuristic ased on the cooling ofmaterialin a heatbath LocalSearch butwe a low moves resulting in solutions ofworse quality than the currentsolution void from localsolutions

s= enerateinitialsolution() t=t(0) while Term ination criteria notmet s =PicktRandom ( N(x) ) If ( f(s ) < f(s) ) then else s= s ccepts as new solution with probability ndif Update(T) end while e f ( s') f ( s) T Initialtem perature Should be suitable high.mostofthe initial moves mustbe accepted (> 0% ) ooling schedule Tem perature is reduced after every m ove Tw o main methods: T T = + β T