Scheduling and Timetabling Integer linear programming models Marjan - - PowerPoint PPT Presentation

scheduling and timetabling
SMART_READER_LITE
LIVE PREVIEW

Scheduling and Timetabling Integer linear programming models Marjan - - PowerPoint PPT Presentation

Scheduling and Timetabling Integer linear programming models Marjan van den Akker 1 I ntro. Marjan van den Akker Lecturer/research Algorithms and Complexity group: Coordination Software- and Gameproject Master courses :


slide-1
SLIDE 1

1

Scheduling and Timetabling

Integer linear programming models Marjan van den Akker

slide-2
SLIDE 2

2

I ntro…….

Marjan van den Akker Lecturer/research Algorithms and Complexity group:

 Coordination Software- and Gameproject  Master courses :

 Algorithms for decision support (COSC),  Advanced Linear Programming (Mastermath)

 Research on planning/scheduling algorithms:

 Public transportation  Sustainable energy systems  integer linear programming and local search  simulation

4/24/2018

slide-3
SLIDE 3

3

Contents

 Integer linear programming is well-known modelling and

solution technique

 First, linear programming  Modelling integer linear programming problems  This is material from the course Algorithms for Decision

Support

 NB: we will work a lot on the blackboard.

You are strongly advised to take notes.

 A beautiful reader is available at

http://www.cs.uu.nl/docs/vakken/mads/LectureNotesILP.pdf

ADS, ILP 1

slide-4
SLIDE 4

4 SCM, chapter 3

4 Exam ple: Ajax

 Three types of computers: Alpha, Beta, and Gamma.  Net profit: $350,- per Alpha, $470,- per Beta, and

$610,- per Gamma.

 Every computer can be sold at the given profit.  Testing: Alpha and Beta computers on the A-line,

Gamma computers on the C-line.

 Testing takes 1 hour per computer.  Capacity A-line: 120 hours; capacity C-line: 80

hours.

 Required labor: 10 hours per Alpha, 15 hours per

Beta, and 20 hours per Gamma.

 Total amount of labor available: 2000 hours.

slide-5
SLIDE 5

5 SCM, chapter 3

5 Exam ple: Ajax

) labor ( 2000 ) line C ( 48 line) (A 120 MC , MB , MA MC 20 MB 15 MA 10 MC MB MA subject to MC 610 MB 470 MA 350 Z max            

Decision variables: MA number of alpha’s produced, etc Objective function Constraints

slide-6
SLIDE 6

6

Linear program m ing

Min cTx s.t. Ax ≤ b x ≥ 0 With ∈ , ∈ ℚ, A ∈ ℚ, and ∈ ℚ

ADS, ILP 1

slide-7
SLIDE 7

7 ADS, ILP 1

LP:Geom etric

5 x2 

1

2 1

  x x

x1 

x2 

2 1

max x x  2

(4,4)

3

2 1

  x x

2 1 3 1 2 3

slide-8
SLIDE 8

8

The sim plex m ethod

 Example dictionaries

ADS, ILP 1

slide-9
SLIDE 9

9 ADS, ILP 1

slide-10
SLIDE 10

10 ADS, ILP 1

slide-11
SLIDE 11

11 ADS, ILP 1

slide-12
SLIDE 12

12 ADS, ILP 1

slide-13
SLIDE 13

13 ADS, ILP 1

slide-14
SLIDE 14

14 ADS, ILP 1

slide-15
SLIDE 15

15 ADS, ILP 1

slide-16
SLIDE 16

16

Solution options

A linear programming problem can

 be infeasible

 Example:

max 6 4 3, 2 2 8, , 0

 be unbounded

 Example:

max 6 4 3, 2 2 8, , 0 for any 0 we have that , is feasible

 have a bounded optimum

ADS, ILP 1

slide-17
SLIDE 17

17

Solution m ethod for linear program m ing

 Simplex method

 Slower than polynomial  Practical

 Ellipsoid method

 Polynomial (Khachian, 1979)  Not practical

 Interior points methods

 Polynomial (Karmakar, 1984)  Outperforms Simplex for very large instances

P LP

ADS, ILP 1

Seems not too hard to

  • implement. But, for larger

problems you run into numerical problems. Use a standard solver (Gurobi, CPLEX, GLPK)

slide-18
SLIDE 18

18

Knapsack problem Item 1:paper 2:book 3:bread 4:smart

  • phone

5:water Utility 8 12 7 15 12 Volume 4 8 5 2 6

Knapsack with volume 15 What should you take with you to maximize utility?

ADS, ILP 1

18

slide-19
SLIDE 19

19

Knapsack problem ( 2 )

max z= 8 x1 + 12 x2 + 7 x3 + 15 x4 + 12 x5 subject to 4 x1 + 8 x2 + 5 x3 + 2 x4 + 6 x5 ≤ 15 x1, x2, x3 , x4 , x5 Є {0,1}

x1 = 1 if item 1 is selected, 0 otherwise, x2, ……

ADS, ILP 1

slide-20
SLIDE 20

20

( Mixed) I nteger linear program m ing

Min cTx + dT y s.t. Ax + By ≤ b x,y ≥ 0 x integral (or binary) Extension of LP:

 Good news: more possibilities for modelling  Bad news: larger solution times

ADS, ILP 1

slide-21
SLIDE 21

21

( Mixed) I nteger linear program Min cTx s.t. Ax + By ≤ b x,y ≥ 0 x integral (or binary) LP-relaxation Min cTx s.t. Ax + By ≤ b x,y ≥ 0 Lower bound (or upper bound in case of maximization)

ADS, ILP 1

slide-22
SLIDE 22

22 ADS, ILP 1

slide-23
SLIDE 23

23

Modeling

Decision variables Objective function Constraints

ADS, ILP 1

slide-24
SLIDE 24

24

Assignm ent problem

 n persons, n jobs.  Each person can do at most one job  Each job has to be executed  Cij cost if person i performs job j  We want to minimize cost

ADS, ILP 1

slide-25
SLIDE 25

25

Maxim um independent set

 Given a graph G=(V,E)  V: nodes  E: edges  An independent set I is a set of nodes, such that every

edge has at most one nodes in I, i.e., no pair of nodes in I is connected.

 What is the maximum number of nodes in an independent

set?

ADS, ILP 1

slide-26
SLIDE 26

26

Facility location Possible locations: n Customers:m

ADS, ILP 1

slide-27
SLIDE 27

27

Capacitated facility location

 Data:

 m customers,  Customer demand: Di  n possible locations of depots (facilities)  cij unit cost of serving customer i by depot j  Capacity depot: Cj  Fixed cost for opening depot DC: Fj

 Which depots are opened and which customer is served by which

depot?

ADS, ILP 1

slide-28
SLIDE 28

28

Uncapacitated facility location

 Data:

 m customers, n possible locations of depot  Each customer is assigned to one depot  dij cost of serving customer i by depot j  Fixed cost for opening depot DC: Fj

 Which depots are opened and which customer is served by which

depot?

ADS, ILP 1

slide-29
SLIDE 29

29

Uncapacitated facility location

 Two models FL and AFL  Same optimal value ZIP for ILP  Set of feasible solutions for LP-relaxation satisfy:

  •  Set of optimal values satisfy:

 The LP-relaxation of FL gives a stronger bound

ADS, ILP 1

slide-30
SLIDE 30

30

Single m achine scheduling exam ple

 n jobs, 1 machine  Job j requires uninterrupted processing time Pj  Job j has release date rj  Machine can process at most one job at a time  We want to minimize the weighted sum of the completion

times

 This is going to be denoted by

ADS, ILP 1