Integrated pollster and vehicle routing S. Gutirrez, A. Miniguano, - - PowerPoint PPT Presentation

integrated pollster and vehicle routing
SMART_READER_LITE
LIVE PREVIEW

Integrated pollster and vehicle routing S. Gutirrez, A. Miniguano, - - PowerPoint PPT Presentation

Integrated pollster and vehicle routing S. Gutirrez, A. Miniguano, D. Recalde, L. M. Torres, R. Torres Centro de Modelizacin Matemtica - ModeMat Escuela Politcnica Nacional - Quito CO@Work 2020 TU Berlin - Berlin Mathematical School,


slide-1
SLIDE 1

Integrated pollster and vehicle routing

  • S. Gutiérrez, A. Miniguano, D. Recalde, L. M. Torres, R. Torres

Centro de Modelización Matemática - ModeMat Escuela Politécnica Nacional - Quito

CO@Work 2020 TU Berlin - Berlin Mathematical School, September 14th - 25th, 2020

slide-2
SLIDE 2

CO@Work 2020 Integrated pollster and vehicle routing

Outline

2

Motivation and problem definition Modeling via mixed integer programming Computational results Conclusions

slide-3
SLIDE 3

CO@Work 2020 Integrated pollster and vehicle routing

The problem

3

The National Statistics Bureau of Ecuador (INEC) is responsible for constructing the Consumer Price Index. To collect information a sample set of stores must be visited monthly. Stores are visited by a set of pollsters, transported by a fleet of hired vehicles. Pollsters also walk between stores. Tasks: schedule visits to stores within time horizon schedule daily service duties for pollsters define daily routes for vehicles

slide-4
SLIDE 4

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

4

1 2 1 2

slide-5
SLIDE 5

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

5

1 2 1 2

t = 0

slide-6
SLIDE 6

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

6

1 2 1 2

t = 2

slide-7
SLIDE 7

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

7

1 2 1 2

t = 5

slide-8
SLIDE 8

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

8

1 2 2

t = 6

1

slide-9
SLIDE 9

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

9

1 2 1 2

t = 15

slide-10
SLIDE 10

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

10

1 2 1 2

t = 17

slide-11
SLIDE 11

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

11

2 2

t = 19

slide-12
SLIDE 12

CO@Work 2020 Integrated pollster and vehicle routing

A Mixed Linear and Integer Programming Problem

12

2 2

t = 20

slide-13
SLIDE 13

CO@Work 2020 Integrated pollster and vehicle routing

E: pollsters K: vehicles set of stores S: days within time horizon for planning

n

Definitions: sets

13

slide-14
SLIDE 14

CO@Work 2020 Integrated pollster and vehicle routing

Definitions: network

14

two nodes for each store (customer)

C− := {1,…, n} C+ := {n + 1,…,2n}

(arrival at store) (departure from store) two nodes for depot

0, 2n + 1

three sets of weighted arcs

AS := {(i, i + n) : i ∈ C−} AW := {(i, j) : i ∈ C+, j ∈ C−, j ≠ i + n} AV := {(i, j) : i, j ∈ C} ∪ {(0,i) : i ∈ C} ∪ {(i,2n + 1) : i ∈ C}

service arcs; service times ti,i+n walking arcs; walk times ti,j vehicle transportation arcs; travel times τi,j

C := C− ∪ C+, V := C ∪ {0,2n + 1}

slide-15
SLIDE 15

CO@Work 2020 Integrated pollster and vehicle routing

Definitions: paths and routes

15

A walking path for a pollster is a simple path from with alternating arcs from the sets A vehicle path is a directed path between two nodes in using

  • nly arcs from

A service route for a pollster is a dipath from to consisting of an alternating sequence of vehicle and walking subpaths. A vehicle route is a vehicle path from to The duration of a route is the sum of its arc weights. A service route is feasible if it does not exceed a maximum allowed duration, including time for a lunch break.

i ∈ C− to j ∈ C+ AS and AW . V AV . 2n + 1 2n + 1.

slide-16
SLIDE 16

CO@Work 2020 Integrated pollster and vehicle routing

Definitions: schedules

16

Vehicles pick-up and deliver pollsters to certain stores; pollsters can share vehicles. Vehicle fleet is homogeneous, with vehicle capacity Pollster “fleet” is homogeneous: any pollster can visit any store. A daily schedule consists of a set of feasible pollster routes and “compatible” vehicle routes, i.e., for any

if is contained in some service route, then it is contained in a vehicle route is not contained in more than service routes

Q . (i, j) ∈ AV :

(i, j) (i, j) Q

slide-17
SLIDE 17

CO@Work 2020 Integrated pollster and vehicle routing

The IPVRP

17

Each store is visited once. Number of working days is minimized. Number of service routes is minimized. Number of vehicle routes is minimized.

Task:

Find a set of daily schedules, at most one for each day in a given time horizon, such that:

slide-18
SLIDE 18

CO@Work 2020 Integrated pollster and vehicle routing

The INEC instance

18

slide-19
SLIDE 19

CO@Work 2020 Integrated pollster and vehicle routing

Outline

19

Motivation and problem definition Modeling via mixed integer programming Computational results Conclusions

slide-20
SLIDE 20

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: variables

20

Binary variables:

xe,s

i,j , ∀(i, j) ∈ AW, e ∈ E, s ∈ S :

xe,s

i,j = 1 ⇔ e walks from i to j on day s

ze,s

i,j , ∀(i, j) ∈ AV, e ∈ E, s ∈ S :

ze,s

i,j = 1 ⇔ e is transported from i to j on day s

be,s

i , ∀i ∈ C−, e ∈ E, s ∈ S :

be,s

i

= 1 ⇔ i is start of a walking path for e on day s f e,s

i , ∀i ∈ C+, e ∈ E, s ∈ S :

f e,s

i

= 1 ⇔ i is end of a walking path of e on day s xe,s

i,i+n, ∀(i, i + n) ∈ AS, e ∈ E, s ∈ S :

xe,s

i,i+n = 1 ⇔ e visits store i on day s

us, s ∈ S : us = 1 ⇔ s has a daily schedule assigned to it

slide-21
SLIDE 21

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: constraints

21

( j,i)∈AV

ze,s

j,i ≤ 1 − xe,s i,i+n + be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

i∈C

z e,s

0,i ≤ us,

∀e ∈ E, s ∈ S . ∑

(i,j)∈AV

ze,s

i,j ≤ 1 − xe,s i−n,i + f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = − f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AW

xe,s

j,i − xe,s i,i+n = − be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, xe,s

i−n,i −

(i,j)∈AW

xe,s

i,j = f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

s∈S ∑ e∈E

xe,s

i,i+n = 1,

∀i ∈ C−,

each store is visited by one pollster on one day

slide-22
SLIDE 22

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: constraints

22

“multicommodity flow demand” constraints for walking paths: sync x variables with b and f variables

( j,i)∈AW

xe,s

j,i − xe,s i,i+n = − be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, xe,s

i−n,i −

(i,j)∈AW

xe,s

i,j = f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i ≤ 1 − xe,s i,i+n + be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

i∈C

z e,s

0,i ≤ us,

∀e ∈ E, s ∈ S . ∑

(i,j)∈AV

ze,s

i,j ≤ 1 − xe,s i−n,i + f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = − f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

s∈S ∑ e∈E

xe,s

i,i+n = 1,

∀i ∈ C−,

slide-23
SLIDE 23

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: constraints

23

( j,i)∈AV

ze,s

j,i ≤ 1 − xe,s i,i+n + be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

(i,j)∈AV

ze,s

i,j ≤ 1 − xe,s i−n,i + f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

i∈C

z e,s

0,i ≤ us,

∀e ∈ E, s ∈ S . ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = − f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AW

xe,s

j,i − xe,s i,i+n = − be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, xe,s

i−n,i −

(i,j)∈AW

xe,s

i,j = f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

s∈S ∑ e∈E

xe,s

i,i+n = 1,

∀i ∈ C−,

degree constraints for vehicle transportation: forbid arcs that cannot connect properly to walking paths

slide-24
SLIDE 24

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: constraints

24

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = − f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i ≤ 1 − xe,s i,i+n + be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

(i,j)∈AV

ze,s

i,j ≤ 1 − xe,s i−n,i + f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

i∈C

z e,s

0,i ≤ us,

∀e ∈ E, s ∈ S . ∑

( j,i)∈AW

xe,s

j,i − xe,s i,i+n = − be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, xe,s

i−n,i −

(i,j)∈AW

xe,s

i,j = f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

s∈S ∑ e∈E

xe,s

i,i+n = 1,

∀i ∈ C−,

“multicommodity flow demand” constraints for vehicle transportation (pickup and delivery of pollsters): sync z variables with b and f variables

slide-25
SLIDE 25

CO@Work 2020 Integrated pollster and vehicle routing

Pollster routing: constraints

25

i∈C

z e,s

0,i ≤ us,

∀e ∈ E, s ∈ S . ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i − ∑ (i,j)∈AV

ze,s

i,j = − f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AV

ze,s

j,i ≤ 1 − xe,s i,i+n + be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

(i,j)∈AV

ze,s

i,j ≤ 1 − xe,s i−n,i + f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

( j,i)∈AW

xe,s

j,i − xe,s i,i+n = − be,s i ,

∀i ∈ C−, e ∈ E, s ∈ S, xe,s

i−n,i −

(i,j)∈AW

xe,s

i,j = f e,s i ,

∀i ∈ C+, e ∈ E, s ∈ S, ∑

s∈S ∑ e∈E

xe,s

i,i+n = 1,

∀i ∈ C−,

degree constraint at depot (at most one route per pollster and day)

slide-26
SLIDE 26

CO@Work 2020 Integrated pollster and vehicle routing

Vehicle routing: variables

26

Binary variables:

yk,s

i,j , ∀(i, j) ∈ AV, k ∈ K, s ∈ S :

yk,s

i,j = 1 ⇔ k travels through (i, j) on day s

Remind that…

ze,s

i,j , ∀(i, j) ∈ AV, e ∈ E, s ∈ S :

ze,s

i,j = 1 ⇔ e is transported from i to j on day s

be,s

i , ∀i ∈ C−, e ∈ E, s ∈ S :

be,s

i

= 1 ⇔ i is start of a walking path for e on day s f e,s

i , ∀i ∈ C+, e ∈ E, s ∈ S :

f e,s

i

= 1 ⇔ i is end of a walking path of e on day s us, s ∈ S : us = 1 ⇔ i if day s has a daily schedule assigned to it

slide-27
SLIDE 27

CO@Work 2020 Integrated pollster and vehicle routing

Vehicle routing: constraints

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

be,s

i ,

∀i ∈ C−, ∀s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

f e,s

i ,

∀i ∈ C+, ∀s ∈ S, ∑

( j,i)∈AV

yk,s

j,i − ∑ (i,j)∈AV

yk,s

i,j = 0,

∀i ∈ C, k ∈ K, s ∈ S, ∑

i∈C

yk,s

0,i ≤ us,

∀k ∈ K, s ∈ S, ∑

e∈E

ze,s

i,j ≤ Q ∑ k∈K

yk,s

i,j ,

∀(i, j) ∈ AV, s ∈ S,

degree constraints for vehicle routes:

  • ut-degree is required to be 1 if node is starting or ending node of some

walking path, 0 otherwise

slide-28
SLIDE 28

CO@Work 2020 Integrated pollster and vehicle routing

Vehicle routing: constraints

( j,i)∈AV

yk,s

j,i − ∑ (i,j)∈AV

yk,s

i,j = 0,

∀i ∈ C, k ∈ K, s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

be,s

i ,

∀i ∈ C−, ∀s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

f e,s

i ,

∀i ∈ C+, ∀s ∈ S, ∑

i∈C

yk,s

0,i ≤ us,

∀k ∈ K, s ∈ S, ∑

e∈E

ze,s

i,j ≤ Q ∑ k∈K

yk,s

i,j ,

∀(i, j) ∈ AV, s ∈ S,

multicommodity flow conservation constraints at stores

slide-29
SLIDE 29

CO@Work 2020 Integrated pollster and vehicle routing

Vehicle routing: constraints

i∈C

yk,s

0,i ≤ us,

∀k ∈ K, s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

be,s

i ,

∀i ∈ C−, ∀s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

f e,s

i ,

∀i ∈ C+, ∀s ∈ S, ∑

( j,i)∈AV

yk,s

j,i − ∑ (i,j)∈AV

yk,s

i,j = 0,

∀i ∈ C, k ∈ K, s ∈ S, ∑

e∈E

ze,s

i,j ≤ Q ∑ k∈K

yk,s

i,j ,

∀(i, j) ∈ AV, s ∈ S,

degree constraint at depot: at most one route allowed per vehicle, and only on days used in the schedule

slide-30
SLIDE 30

CO@Work 2020 Integrated pollster and vehicle routing

Vehicle routing: constraints

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

be,s

i ,

∀i ∈ C−, ∀s ∈ S, ∑

k∈K ∑ (i,j)∈AV

yk,s

i,j = ∑ e∈E

f e,s

i ,

∀i ∈ C+, ∀s ∈ S, ∑

( j,i)∈AV

yk,s

j,i − ∑ (i,j)∈AV

yk,s

i,j = 0,

∀i ∈ C, k ∈ K, s ∈ S, ∑

i∈C

yk,s

0,i ≤ us,

∀k ∈ K, s ∈ S, ∑

e∈E

ze,s

i,j ≤ Q ∑ k∈K

yk,s

i,j ,

∀(i, j) ∈ AV, s ∈ S,

each arc of the vehicle network with positive transportation demand must be covered by a vehicle route; account for vehicle capacity Q

slide-31
SLIDE 31

CO@Work 2020 Integrated pollster and vehicle routing

Shift-length: Variables

31

Variables

Bi, ∀i ∈ V : arrival Tme at store i, if i ∈ C− departure Tme from store i − n, if i ∈ C+ duraTon of longest service route, if i = 2n + 1 equals to 0 for consistency, if i = 0 we,s

i , ∀i ∈ C−, e ∈ E, s ∈ S :

we,s

i

= 1 ⇔ e takes break aUer visiTng i on day s (in this case, Bi+n is departure Tme aUer break)

Remind that…

yk,s

i,j , ∀(i, j) ∈ AV, k ∈ K, s ∈ S :

yk,s

i,j = 1 ⇔ k travels through (i, j) on day s

xe,s

i,j , ∀(i, j) ∈ AW, e ∈ E, s ∈ S :

xe,s

i,j = 1 ⇔ e walks from i to j on day s

slide-32
SLIDE 32

CO@Work 2020 Integrated pollster and vehicle routing

Shift-length: constraints

32

Bi+n ≥ Bi + ti,i+n + P∑

e∈E ∑ s∈S

we,s

i ,

∀(i, i + n) ∈ AS, Bj ≥ Bi + ti,j − M(1 − ∑

e∈E ∑ s∈S

xe,s

i,j ),

∀(i, j) ∈ AW, Bj ≥ Bi + τi,j − M(1 − ∑

k∈K ∑ s∈S

yk,s

i,j ),

∀(i, j) ∈ AV, B0 = 0, B2n+1 ≤ Bmax,

account for service times and duration of lunch breaks (when present) Parameters: ti,i+n : service Tme at store i

P : duraTon of lunch break

slide-33
SLIDE 33

CO@Work 2020 Integrated pollster and vehicle routing

account for walking times between consecutive stores Parameters:

ti,j : walking Tme from store i to store j

M : sufficiently large constant

Shift-length: constraints

33

Bj ≥ Bi + ti,j − M(1 − ∑

e∈E ∑ s∈S

xe,s

i,j ),

∀(i, j) ∈ AW, Bi+n ≥ Bi + ti,i+n + P∑

e∈E ∑ s∈S

we,s

i ,

∀(i, i + n) ∈ AS, Bj ≥ Bi + τi,j − M(1 − ∑

k∈K ∑ s∈S

yk,s

i,j ),

∀(i, j) ∈ AV, B0 = 0, B2n+1 ≤ Bmax,

slide-34
SLIDE 34

CO@Work 2020 Integrated pollster and vehicle routing

Shift-length: constraints

34

Bj ≥ Bi + τi,j − M(1 − ∑

k∈K ∑ s∈S

yk,s

i,j ),

∀(i, j) ∈ AV, Bj ≥ Bi + ti,j − M(1 − ∑

e∈E ∑ s∈S

xe,s

i,j ),

∀(i, j) ∈ AW, Bi+n ≥ Bi + ti,i+n + P∑

e∈E ∑ s∈S

we,s

i ,

∀(i, i + n) ∈ AS, B0 = 0, B2n+1 ≤ Bmax,

account for vehicle transportation times when j=2n+1, maximum duration of a service route is computed Parameters:

τi,j : driving Tme from node i to node j

M : sufficiently large constant

slide-35
SLIDE 35

CO@Work 2020 Integrated pollster and vehicle routing

Shift-length: constraints

35

Bi+n ≥ Bi + ti,i+n + P∑

e∈E ∑ s∈S

we,s

i ,

∀(i, i + n) ∈ AS, Bj ≥ Bi + ti,j − M(1 − ∑

e∈E ∑ s∈S

xe,s

i,j ),

∀(i, j) ∈ AW, Bj ≥ Bi + τi,j − M(1 − ∑

k∈K ∑ s∈S

yk,s

i,j ),

∀(i, j) ∈ AV, B0 = 0, B2n+1 ≤ Bmax,

specify upper bound for route length Parameters:

Bmax : maximum allowed duraTon for a route

slide-36
SLIDE 36

CO@Work 2020 Integrated pollster and vehicle routing

each break must start within the prescribed time window Parameters: [T0, T1] : Tme window for starTng of break

ti,i+n : service Tme at store i

M : sufficiently large constant

Pollster breaks

36

T0∑

e∈E ∑ s∈S

we,s

i

≤ Bi + ti,i+n ≤ T1 + M(1 − ∑

e∈E ∑ s∈S

we,s

i ),

∀i ∈ C−, we,s

i

≤ xe,s

i,i+n,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

i∈C−

we,s

i

= ∑

j∈C

ze,s

0,j,

∀e ∈ E, s ∈ S .

slide-37
SLIDE 37

CO@Work 2020 Integrated pollster and vehicle routing

Pollster breaks

37

we,s

i

≤ xe,s

i,i+n,

∀i ∈ C−, e ∈ E, s ∈ S, T0∑

e∈E ∑ s∈S

we,s

i

≤ Bi + ti,i+n ≤ T1 + M(1 − ∑

e∈E ∑ s∈S

we,s

i ),

∀i ∈ C−, ∑

i∈C−

we,s

i

= ∑

j∈C

ze,s

0,j,

∀e ∈ E, s ∈ S .

pollster e can have a break after visiting store i on day s only if e serves i on that day sync variables we,s

i and xe,s i,i+n :

slide-38
SLIDE 38

CO@Work 2020 Integrated pollster and vehicle routing

Pollster breaks

38

T0∑

e∈E ∑ s∈S

we,s

i

≤ Bi + ti,i+n ≤ T1 + M(1 − ∑

e∈E ∑ s∈S

we,s

i ),

∀i ∈ C−, we,s

i

≤ xe,s

i,i+n,

∀i ∈ C−, e ∈ E, s ∈ S, ∑

i∈C−

we,s

i

= ∑

j∈C

ze,s

0,j,

∀e ∈ E, s ∈ S .

each pollster must have exactly one break on each service route

slide-39
SLIDE 39

CO@Work 2020 Integrated pollster and vehicle routing

Objective function

39

min κ0∑

s∈S

us + κ1∑

s∈S ∑ k∈K ∑ i∈C

yk,s

0,i + κ2∑ s∈S ∑ k∈K ∑ i∈C

ze,s

0,i

number of working days total number of required vehicle routes total number of required service routes Weighted sum of three components:

slide-40
SLIDE 40

CO@Work 2020 Integrated pollster and vehicle routing

Outline

40

Motivation and problem definition Modeling via mixed integer programming Computational results Conclusions

slide-41
SLIDE 41

CO@Work 2020 Integrated pollster and vehicle routing

A “toy” example

41

slide-42
SLIDE 42

CO@Work 2020 Integrated pollster and vehicle routing

A “toy” example

42

8 stores, 2 pollsters, 2 vehicles, 3 days in time horizon Service times in [2; 34] Walking times in [2; 28] Vehicle transportation times in [0.5; 7] Shift-duration Pollster break duration , and must be taken in time window [50; 90] Objective coefficients: .

Bmax = 120. P = 20 κ0 = 200, κ1 = 100, κ2 = 40

slide-43
SLIDE 43

CO@Work 2020 Integrated pollster and vehicle routing

A “toy” example

43

MIP model with 2923 variables and 1683 rows Solved on a MacBook Pro i7 2.6GHz with 16 GB RAM, OSX Catalina, using Gurobi 9.0.2. as MIP solver, TimeLimit= 3600 (best solution and LB found within 2 minutes) Feasible solution with Gap= 52.8%, uses 1 vehicle, 2 pollsters, 2 days. Day 1 Day 2

slide-44
SLIDE 44

CO@Work 2020 Integrated pollster and vehicle routing

The INEC complete instance

44

C = 820 S = 17

E = 5 K = 3 Bmax = 84

P = 8

Q = 3 A ≈ 3.4 mio. > 420 mio. binary variables!!! > 49 mio. constraints!!!

slide-45
SLIDE 45

CO@Work 2020 Integrated pollster and vehicle routing

Solution approach

45

  • 1. Partition of the instance

into “half-daily snapshots” (Graph partitioning techniques balancing total node load).

  • 2. Route pollsters and vehicles

in each partition without considering lunch break.

  • 3. Link partial solutions into a

complete schedule (matching problem).

slide-46
SLIDE 46

CO@Work 2020 Integrated pollster and vehicle routing

Results

46

30 partitions were constructed, with aggregated service times varying in [90; 210] Each partition could be solved using 2 pollsters and 1 vehicle. Global schedule with 15 days, 2 pollsters, 1 vehicle (previous: 17d, 5p, 3v)

slide-47
SLIDE 47

CO@Work 2020 Integrated pollster and vehicle routing

Outline

47

Motivation and problem definition Modeling via mixed integer programming Computational results Conclusions

slide-48
SLIDE 48

CO@Work 2020 Integrated pollster and vehicle routing

Conclusions & Outlook

48

Presented a mixed integer programming model for scheduling the visits of stores and for the integrated routing of pollsters and vehicles. Model is hard to solve even for small “toy” instances. Three-phase solution approach: partitioning of stores into half-day instances, partial (vehicle + pollster) routing, matching of routes into daily schedules. Tighten linear relaxation (cutting-planes) Alternative formulations? Column generation?

Future work

slide-49
SLIDE 49

CO@Work 2020 Integrated pollster and vehicle routing 49

Thank you for your attention!!!