Robust Optimization of Performance Scheduling Problem under Accepting Strategy

ABSTRACT

In this paper, the problem of program performance scheduling with accepting strategy is studied. Considering the uncertainty of actual situation, the duration of a program is expressed as a bounded interval. Firstly, we decide which programs are accepted. Secondly, the risk preference coefficient of the decision maker is introduced. Thirdly, the min-max robust optimization model of the uncertain program show scheduling is built to minimize the performance cost and determine the sequence of these programs. Based on the above model, an effective algorithm for the original problem is proposed. The computational experiment shows that the performance’s cost (revenue) will increase (decrease) with decision maker’s risk aversion.

In this paper, the problem of program performance scheduling with accepting strategy is studied. Considering the uncertainty of actual situation, the duration of a program is expressed as a bounded interval. Firstly, we decide which programs are accepted. Secondly, the risk preference coefficient of the decision maker is introduced. Thirdly, the min-max robust optimization model of the uncertain program show scheduling is built to minimize the performance cost and determine the sequence of these programs. Based on the above model, an effective algorithm for the original problem is proposed. The computational experiment shows that the performance’s cost (revenue) will increase (decrease) with decision maker’s risk aversion.

KEYWORDS

Performance Scheduling, Robust Optimization, Duality Theory, 0 - 1 Mixed Linear Programming

Performance Scheduling, Robust Optimization, Duality Theory, 0 - 1 Mixed Linear Programming

1. Introduction

In the planning stage of a performance, many programs will sign up for the performance, and each program bring different revenue and different cost. The program group needs to comprehensively consider the performance revenue and cost, and decides which programs to be accepted. Each program requires multiple performers, and a performer can also participate in multiple programs. If the programs which the same performer participates in are not consecutively arranged, the waiting cost will occur. Therefore, after determining the set of the accepted programs, the reasonable scheduling of these accepted programs to increase the performance revenue is the main concern of the decision maker.

Scholars have done extensive researches on program performance scheduling and the similar film production scheduling problem. Cheng et al. [1] first mentioned scheduling problems in the film production process. They assumed that the shooting duration of each scene is a certain value, considering the actor waiting cost problem during film shooting. They proved that the problem is strong NP-hard, and used the branch and bound algorithm to solve the minimum cost problem. Nordström and Tufekci [2] proposed several hybrid algorithms which use limited pairwise interchange procedure within the simple genetic algorithm framework to solve the problem proposed by Cheng et al. [1] . Their algorithms outperformed in terms of quality of solution and computational time. Bomsdorf and Derigs [3] presented the movie shooting scheduling problem and formulated a conceptual model. And they proposed a meta-heuristic algorithm to generate a timeline for film. Stuckey et al. [4] applied a dynamic programming to the scheduling problem to minimize the cost of the talent. They assumed that the actor’s appearance fee is different, and the shooting time of each scene is a certain value, showed a number of ways to improve the dynamic programming solution by preprocessing and restricting the search. Wang et al. [5] generalized a scheduling model by incorporating the performers waiting cost and operating cost in film shooting. They used the next fit (NF) algorithm and the first fit decreasing (FFD) algorithm to allocate scenes to work days so as to provide initial solutions for further improvements. Dynamic programming, iterated local search, and tabu search are adopted to constitute the second-phase improvement procedures. Qin et al. [6] formulated the talent scheduling problem as an integer linear programming model and designed an improved branch and bound method to deal with it. Sakulsom and Tharmmaphornphilas [7] studied a music performance scheduling problem. The objective is to minimize the total number of days that all performers have to show up, and sequence the music pieces within each day to minimize the total waiting time of the performers. They proposed a 2-stage methodology to schedule music pieces, which is a combination of a cell formation technique and an integer-programming model.

Based on the above research results, it can be concluded that the film shooting or program rehearsal duration of the predecessors is assumed to be a certain value. In the actual program performances, due to factors such as staff absence, equipment failure, performance effects, etc., the duration of each program is uncertain. Therefore, the results obtained by the deterministic research method are greatly deviated from the actual situation. Zhen et al. [8] proposed the performance scheduling problem of mixed duration, and expressed the performance time as the sum of the certain performance time and the uncertain adjustment time. They assumed that the uncertain adjustment time obeyed the normal distribution. The objective is to minimize the total waiting cost of the performers. However, in the actual situations, a large amount of data is often required to obtain a distribution function of a random parameter. This paper breaks through this limitation and uses the robust optimization method to express the uncertain duration as a continuous bounded interval. We only need to know the upper and lower bound of the program performance duration.

Robust optimization is an effective method to solve the uncertain problem and has been widely used. Wang and Tang [9] proposed a two-stage robust optimization method for interval-type surgical scheduling problem, effectively reducing the adverse effects of service time uncertainty on hospital revenue. Xu et al. [10] built a robust scheduling model for homogeneous parallel machines based on the min-max regret criterion under the condition that only knew the processing time interval. Qiu et al. [11] used the robust optimization method to solve the order policy of the integrated supply chain and the contract coordination policy of the distributed supply chain under the min-max regret value criterion. Zhang et al. [12] built an emergency rescue network based on the scenario of min-max regret value criteria, constructed a robust optimization model, and transformed it into a mixed integer programming model, and they used the scenario relaxation algorithm to solve this model.

This paper considers the uncertain duration program performance scheduling problem under accepting strategy (recorded as P0), the remainder is organized as follows. In Section 2, we use a simple example to describe the program performance scheduling problem and introduce the application of the min-max robust optimization method in this paper. In the case of determining the set of the accepted performance programs, the decision maker’s risk preference coefficient [13] [14] is introduced and a robust performance scheduling model (RPSM) is built for these accepted programs in Section 3. Then we transform the RPSM into a 0 - 1 mixed linear programming model to minimize the performance cost, and based on the algorithm for solving RPSM, the algorithm H of P0 is constructed to determine the accepted or rejected programs and the performance sequence of these accepted programs in Section 4. Finally we use Matlab software to carry out numerical experiments, verify the actual performance of algorithm H, and compare the influence of decision maker’s risk preference on performance cost.

2. Problem Description

In the actual performance, due to limitations of time or layout and so on, decision maker needs to consider program revenue and cost in a comprehensive manner, select a part of the programs from many registration programs to perform, and arrange the performance sequence of these accepted programs. Accepting a program will generate revenue, and rejecting a program will occur penalty cost. After determining the set of the accepted performance programs, multiple performers will participate in the performance. If there is no comprehensive arrangement for a performer who participates in different programs, the performer will generate waiting time and waiting cost. The objective of this paper is to find a solution that maximizes the overall revenue of the program group, including determining the set of the accepted performance programs and scheduling these accepted performance programs. The overall revenue of the program group consists of the following two parts: 1) the revenue value of the accepted programs and the penalty cost of the rejected programs. 2) the performance cost of the performers at the performance scene, including the appearance fees and the waiting cost.

Let’s look at the following example. A program group receives 6 programs $\left\{{s}_{1},{s}_{2},{s}_{3},{s}_{4},{s}_{5},{s}_{6}\right\}$ . These programs are performed by 3 performers $\left\{{a}_{1},{a}_{2},{a}_{3}\right\}$, and the performers participation list is shown in Table 1 (In this paper, √ indicates that the performer participates in the program, and × indicates that the performer does not participate in the program). The performance duration of the programs $\left\{{s}_{1},{s}_{2},{s}_{3},{s}_{4},{s}_{5},{s}_{6}\right\}$ is $\left\{\left[3,5\right],\left[5,8\right],\left[4,6\right],\left[3,6\right],\left[4,5\right],\left[5,7\right]\right\}$, respectively, the performance revenue from the accepted programs $\left\{{s}_{1},{s}_{2},{s}_{3},{s}_{4},{s}_{5},{s}_{6}\right\}$ is $\left\{175,100,145,155,110,160\right\}$, respectively, and the penalty cost of rejecting them is $\left\{30,20,33,40,25,35\right\}$, respectively.

Assume that the program group decides to accept the programs
$\left\{{s}_{1},{s}_{3},{s}_{4},{s}_{6}\right\}$ to perform, the performance revenue from these accepted programs
$\left\{{s}_{1},{s}_{3},{s}_{4},{s}_{6}\right\}$ is 635, and the penalty cost for these rejected programs
$\left\{{s}_{2},{s}_{5}\right\}$ is 45. After determining the set of the accepted programs, in order to improve efficiency and reduce cost, we assume that performers show up on time before the first program they play starts, leave immediately after the last program they play finishes. In this paper, the robust optimization method is used to arrange the sequence of these accepted programs
$\left\{{s}_{1},{s}_{3},{s}_{4},{s}_{6}\right\}$, which is based on the principle of “min-max” [15] . “Min-max” means that after formulating a program performance scheduling scheme, due to the uncertainty of the performance duration, the overall revenue of the program group may have multiple possibilities. We calculate the maximum performance cost for each possible scheduling scheme, and then select the least one in these maximum costs. Table 2 lists the two scheduling schemes (the blank space in the table indicates that the performer has finished performance and left the scene).
${t}_{1},{t}_{2},{t}_{3}$ and
${t}_{4}$ indicate the time series. In the scheme 1, a_{1} is idle in t_{3}, waiting time is the performance duration of s_{4}. a_{2} and a_{3} are idle in t_{2}, waiting time is the performance duration of s_{3}. The

Table 1. Performer participation list.

Table 2. Program performance schemes.

maximum performance cost under scheme 1 is calculated as c_{1}. In the scheme 2, only a_{1} is idle at t_{2}, and the waiting time is the performance duration of s_{4}, and the maximum performance cost is calculated as c_{2}. Our robust optimization method is to calculate the maximum performance cost of all feasible solutions, then pick the solution with the least cost in the maximum performance cost to determine the performance sequence of these accepted programs. Based on the above, the overall revenue of the program group is obtained by comprehensively considering the revenue of the accepted programs and the penalty cost of the rejected programs.

3. Performance Scheduling Model

3.1. Model Hypothesis

1) n performers $A=\left\{{a}_{1},{a}_{2},\cdots ,{a}_{n}\right\}$ participate in k programs. If the performer ${a}_{i}$ participates in the program ${s}_{j}$, defining ${w}_{ij}$ to be 1, 0 otherwise, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,k\right\}$ . The program group needs to select several programs from k registration programs to perform. Suppose that the number of accepted programs is m (m is part of our decision). $T=\left\{1,2,\cdots ,m\right\}$ represents the set of time series, arranged from small to large, each program occupies a time series.

2) The duration of program ${q}_{j}$ is an uncertain number, ${q}_{j}\in \left[l{q}_{j},u{q}_{j}\right]$, $l{q}_{j}$, $u{q}_{j}$ value is known.

3) Unit time waiting cost of performer ${a}_{i}$ is $l{c}_{i}$, unit time appearance fee of performer ${a}_{i}$ is $u{c}_{i}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$ .

4) Performer ${a}_{i}$ shows up on time before the first program he play starts, leave immediately after the last program he play finishes.

5) The revenue of the accepted program is ${p}_{j}$, the penalty cost of the rejected program is ${f}_{j}$, $\forall j\in \left\{1,2,\cdots ,k\right\}$ . The objective function of this paper is the revenue of accepted programs minus the penalty cost of rejected programs, and then minus the performance cost of these accepted programs. The performance cost of the accepted programs is the appearance fees and waiting cost of the performers at the scene.

3.2. Robust Performance Scheduling Model of the Accepted Program (RPSM)

After determining the set of the accepted performance program, it is important to properly schedule these accepted programs and reduce the performance cost so that the program group can achieve a better overall revenue. Without loss of generality, let’s assume that the accepted programs are $\left\{{s}_{1},{s}_{2},\cdots ,{s}_{m}\right\}$, the decision variables are as follows:

${x}_{jt}=1$ if program ${s}_{j}$ performs in time series t, 0 otherwise.

${y}_{it}=1$ if performer ${a}_{i}$ performs in time series t, 0 otherwise.

${a}_{it}=1$ if performer ${a}_{i}$ arrives at the scene before time series t (including t), 0 otherwise.

${l}_{it}=1$ if performer ${a}_{i}$ leaves the scene after time series t (including t), 0 otherwise.

${d}_{it}=1$ if performer ${a}_{i}$ waits at the scene in time series t, 0 otherwise.

Let ${\theta}_{j}=\frac{{q}_{j}-l{q}_{j}}{u{q}_{j}-l{q}_{j}}$ be the degree which the performance duration ${q}_{j}$ of

program ${s}_{j}$ deviates from the lower bound $l{q}_{j}$, ${\theta}_{j}\in \left[0,1\right]$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ . In this paper, the idea of Bertsimas and Sim [13] [14] is used to control the conservative degree of robust optimal scheduling model, and the risk preference

coefficient $\mu $ of decision maker is introduced, $\underset{j=1}{\overset{m}{\sum}}{\theta}_{j}}={\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{{q}_{j}-l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\le \mu $,

$0\le \mu \le m$, indicating that up to $\mu $ programs which duration reaches the upper bound at the same time. The $\mu $ value is given by the decision maker in advance, and the more conservative the decision maker is, the larger the $\mu $ value is. Therefore, the uncertain set of the program duration ${q}_{j}$ is expressed as:

$\left\{{q}_{j}|l{q}_{j}\le {q}_{j}\le u{q}_{j},{\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{{q}_{j}-l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\le \mu ,\forall j\in \left\{1,2,\cdots ,m\right\}\right\}$ (1)

We build the following min-max robust performance scheduling model for the accepted programs (RPSM):

$\mathrm{min}\left\{\mathrm{max}\left\{{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {\displaystyle \underset{j=1}{\overset{m}{\sum}}{q}_{j}\cdot {x}_{jt}+{\displaystyle \underset{i=1}{\overset{n}{\sum}}u{c}_{i}{\displaystyle \underset{j=1}{\overset{m}{\sum}}{q}_{j}\cdot {w}_{ij}}}}}}\right\}\right\}$ (2)

s.t.

$\underset{t=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (3)

$\underset{j=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall t\in T$ (4)

${y}_{it}-{\displaystyle \underset{j=1}{\overset{m}{\sum}}{x}_{jt}\cdot {w}_{ij}=0}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (5)

${a}_{it}+{l}_{it}-{h}_{it}\le 1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (6)

${a}_{it}-{a}_{i,t+1}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (7)

${l}_{i,t+1}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (8)

$-{a}_{it}-{l}_{it}\le -1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (9)

${y}_{it}-{a}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (10)

${y}_{it}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (11)

${h}_{it}-{d}_{it}-{y}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (12)

$l{q}_{j}\le {q}_{j}\le u{q}_{j}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (13)

$\underset{j=1}{\overset{m}{\sum}}\frac{{q}_{j}-l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\le \mu $ (14)

${x}_{jt},{y}_{it},{h}_{it},{a}_{it},{l}_{it},{d}_{it}\in \left\{0,1\right\}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$, $\forall t\in T$ (15)

Equation (2) is the minimum performance cost in the worst case of the uncertain set, including the waiting cost and appearance fees of the performers at the scene. Equation (3) and Equation (4) force to schedule only one program at one time series and assign each program to only one time series. Equations (5)-(9) indicate the presence of performers in each time series. Equation (10) and Equation (11) indicate that each performer shows up on time before the first program he play starts, leave immediately after the last program he play finishes. Equation (12) determines if the performer is in a wait state at a time series. Equation (13) and Equation (14) constitute an uncertain set of program performance duration. Equation (15) indicates that the decision variables are a 0 - 1 variable.

4. Solve Model

This paper studies the uncertain duration program performance scheduling problem under accepting strategy (P0). Firstly, we decide to accept some programs from the registration programs; secondly, we build a robust performance scheduling model for the accepted programs (RPSM) to get the program performance sequence and performance cost; Finally, based on the above, we consider the revenue of the accepted programs and the penalty cost of the rejected programs comprehensively, and determine a feasible scheduling scheme for P0 to maximize the performance revenue.

4.1. Solve the Robust Performance Scheduling Model of the Accepted Programs (RPSM)

We observe the Equations (2)-(15) of RPSM, finding that only the objective function (2), Equation (13) and Equation (14) are affected by the uncertainty of duration ${q}_{j}$ . Therefore, RPSM can be expressed as a two-stage robust optimization model. The decision variables of the first stage are ${x}_{jt}$, ${y}_{it}$, ${h}_{it}$, ${a}_{it}$, ${l}_{it}$, ${d}_{it}$, the decision variable of the second stage is ${q}_{j}$ . Then the two-stage robust optimization model can be expressed as:

$\mathrm{min}\text{F}\left({\left\{{q}_{j}\right\}}_{j\in \left\{1,2,\cdots ,m\right\}}\right)$ (16)

The constraints are Equations (3)-(12) and Equation (15).

And:

$\text{F}\left({\left\{{q}_{j}\right\}}_{j\in \left\{1,2,\cdots ,m\right\}}\right)=\mathrm{max}\left\{{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {\displaystyle \underset{j=1}{\overset{m}{\sum}}{q}_{j}\cdot {x}_{jt}+{\displaystyle \underset{i=1}{\overset{n}{\sum}}u{c}_{i}{\displaystyle \underset{j=1}{\overset{m}{\sum}}{q}_{j}\cdot {w}_{ij}}}}}}\right\}$ (17)

s.t.

$l{q}_{j}\le {q}_{j}\le u{q}_{j}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (18)

$\underset{j=1}{\overset{m}{\sum}}\frac{{q}_{j}-l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\le \mu $ (19)

Noticing that (17)-(19) is a linear programming problem for ${q}_{j}$, which is equivalent to the following form:

$\mathrm{max}\left\{{\displaystyle \underset{j=1}{\overset{m}{\sum}}{q}_{j}\cdot \left({\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {x}_{jt}+{\displaystyle \underset{i=1}{\overset{n}{\sum}}u{c}_{i}\cdot {w}_{ij}}}}\right)}\right\}$ (20)

s.t.

$-{q}_{j}\le -l{q}_{j}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (21)

${q}_{j}\le u{q}_{j}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (22)

$\underset{j=1}{\overset{m}{\sum}}\frac{{q}_{j}}{u{q}_{j}-l{q}_{j}}\le \mu +{\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{l{q}_{j}}{u{q}_{j}-l{q}_{j}}}$ (23)

Since the feasible domain of the linear programming problem is bounded, according to the strong dual theory [16] , it can be known that the original maximization problem can be equivalent to the minimization problem. We define ${\rho}_{j}$ as the dual variable of ${q}_{j}$ in Equations (21)-(23), $\forall j\in \left\{1,2,\cdots ,2m+1\right\}$ . Then the dual programming problem of (20)-(23) is as the following:

$\mathrm{min}{\displaystyle \underset{j=1}{\overset{m}{\sum}}\left(-l{q}_{j}\right)\cdot {\rho}_{j}}+{\displaystyle \underset{j=m+1}{\overset{2m}{\sum}}u{q}_{j-m}\cdot {\rho}_{j}+{\rho}_{2m+1}\cdot \left(\mu +{\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\right)}$ (24)

s.t.

$-{\rho}_{j}+{\rho}_{m+j}+\frac{1}{u{q}_{j}-l{q}_{j}}\cdot {\rho}_{2m+1}={\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {x}_{jt}}}+{\displaystyle \underset{i=1}{\overset{n}{\sum}}{w}_{ij}\cdot u{c}_{i}}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (25)

${\rho}_{j}\ge 0$, $\forall j\in \left\{1,2,\cdots ,2m+1\right\}$ (26)

Bringing the Equations (24)-(26) to the RPSM to get the equivalent model RPSM1:

$\mathrm{min}{\displaystyle \underset{j=1}{\overset{m}{\sum}}\left(-l{q}_{j}\right)\cdot {\rho}_{j}}+{\displaystyle \underset{j=m+1}{\overset{2m}{\sum}}u{q}_{j-m}\cdot {\rho}_{j}+{\rho}_{2m+1}\cdot \left(\mu +{\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\right)}$ (27)

s.t.

$\underset{t=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (28)

$\underset{j=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall t\in T$ (29)

${y}_{it}-{\displaystyle \underset{j=1}{\overset{m}{\sum}}{x}_{jt}\cdot {w}_{ij}=0}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (30)

${a}_{it}+{l}_{it}-{h}_{it}\le 1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (31)

${a}_{it}-{a}_{i,t+1}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (32)

${l}_{i,t+1}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (33)

$-{a}_{it}-{l}_{it}\le -1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (34)

${y}_{it}-{a}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (35)

${y}_{it}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (36)

${h}_{it}-{d}_{it}-{y}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (37)

$-{\rho}_{j}+{\rho}_{m+j}+\frac{1}{u{q}_{j}-l{q}_{j}}\cdot {\rho}_{2m+1}={\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {x}_{jt}+{\displaystyle \underset{i=1}{\overset{n}{\sum}}{w}_{ij}\cdot u{c}_{i}}}}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (38)

${\rho}_{j}\ge 0$, $\forall j\in \left\{1,2,\cdots ,2m+1\right\}$ (39)

${x}_{jt},{y}_{it},{h}_{it},{a}_{it},{l}_{it},{d}_{it}\in \left\{0,1\right\}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$, $\forall t\in T$ (40)

It is observed that the Equation (38) contains nonlinear part $\underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {d}_{it}\cdot {x}_{jt}}$ .

In order to convert the nonlinear constraints into linear constraints, we introduce variable ${\Delta}_{it}^{j}$, let ${\Delta}_{it}^{j}\ge {d}_{it}+{x}_{jt}-1$, ${\Delta}_{it}^{j}\in \left\{0,1\right\}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$, $\forall t\in T$ . Then we can obtain the following 0-1 mixed linear programming model RPSM2:

$\mathrm{min}{z}_{1}={\displaystyle \underset{j=1}{\overset{m}{\sum}}\left(-l{q}_{j}\right)\cdot {\rho}_{j}}+{\displaystyle \underset{j=m+1}{\overset{2m}{\sum}}u{q}_{j-m}\cdot {\rho}_{j}+{\rho}_{2m+1}\cdot \left(\mu +{\displaystyle \underset{j=1}{\overset{m}{\sum}}\frac{l{q}_{j}}{u{q}_{j}-l{q}_{j}}}\right)}$ (41)

s.t.

$\underset{t=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (42)

$\underset{j=1}{\overset{m}{\sum}}{x}_{jt}=1$, $\forall t\in T$ (43)

${y}_{it}-{\displaystyle \underset{j=1}{\overset{m}{\sum}}{x}_{jt}\cdot {w}_{ij}=0}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (44)

${a}_{it}+{l}_{it}-{h}_{it}\le 1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (45)

${a}_{it}-{a}_{i,t+1}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (46)

${l}_{i,t+1}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in \left\{1,2,\cdots ,m-1\right\}$ (47)

$-{a}_{it}-{l}_{it}\le -1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (48)

${y}_{it}-{a}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (49)

${y}_{it}-{l}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (50)

${h}_{it}-{d}_{it}-{y}_{it}\le 0$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall t\in T$ (51)

$-{\rho}_{j}+{\rho}_{m+j}+\frac{1}{u{q}_{j}-l{q}_{j}}\cdot {\rho}_{2m+1}-{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\displaystyle \underset{t=1}{\overset{m}{\sum}}l{c}_{i}\cdot {\Delta}_{it}^{j}}}={\displaystyle \underset{i=1}{\overset{n}{\sum}}{w}_{ij}\cdot u{c}_{i}}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$ (52)

${x}_{jt}+{d}_{it}-{\Delta}_{it}^{j}\le 1$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$, $\forall t\in T$ (53)

${\rho}_{j}\ge 0$, $\forall j\in \left\{1,2,\cdots ,2m+1\right\}$ (54)

${x}_{jt},{y}_{it},{h}_{it},{a}_{it},{l}_{it},{d}_{it},{\Delta}_{it}^{j}\in \left\{0,1\right\}$, $\forall i\in \left\{1,2,\cdots ,n\right\}$, $\forall j\in \left\{1,2,\cdots ,m\right\}$, $\forall t\in T$ (55)

It is observed that the difference between the above two models is the Equation (38) in RPSM1 and the Equation (52), Equation (53) in RPSM2. However, RPSM1 and RPSM2 have the same optimal solution and objective function value. The reason is the coefficients symbol of ${\rho}_{j}$, ${\rho}_{m+j}$, ${\rho}_{2m+1}$ in Equation (41) and Equation (52) are the same (positive or negative simultaneously). So the change of the variable value in Equation (52) will have the same effect in Equation (41). From Equation (53), it can be known that when at least one of ${d}_{it}$ and ${x}_{jt}$ is 0, ${\Delta}_{it}^{j}$ can be 1 or 0. However, in the process of finding the optimal solution, ${\Delta}_{it}^{j}$ should be 0 as much as possible, because if ${\Delta}_{it}^{j}=1$, then in the Equation (52), ${\rho}_{j}$ decreases or ${\rho}_{m+j}$, ${\rho}_{2m+1}$ increases, in the corresponding Equation (41), ${\rho}_{j}$ decreases or ${\rho}_{m+j}$, ${\rho}_{2m+1}$ increases, the objective value increases(When other variables remain unchanged). The performance of this solution is worse than the solution corresponding to ${\Delta}_{it}^{j}=0$ . So when at least one of ${d}_{it}$ and ${x}_{jt}$ is 0, ${\Delta}_{it}^{j}=0$, only when ${d}_{it}={x}_{jt}=1$, ${\Delta}_{it}^{j}=1$ . In summary, variable ${\Delta}_{it}^{j}$ is introduced to transform the nonlinear programming RPSM1 into linear programming RPSM2 successfully.

4.2. Performance Scheduling Algorithm H

We denote $\stackrel{\xaf}{{p}_{j}}=\frac{{p}_{j}}{\frac{l{q}_{j}+u{q}_{j}}{2}}$, and sort the registration programs from large to

small according to $\stackrel{\xaf}{{p}_{j}}$ value. Without loss of generality, we assume that the program sequence is $\left\{{s}_{1},{s}_{2},\cdots ,{s}_{k}\right\}$ . The algorithm H for solving P0 is described as follows:

Algorithm H:

Step 1. Select a set of program number $\left\{{m}_{1},{m}_{2},\cdots ,{m}_{x}\right\}$ that may be accepted, $1\le {m}_{1},{m}_{2},\cdots ,{m}_{x}\le k$ . $\tau =1$ .

Step 2. The number of the accepted programs is ${m}_{\tau}$, the set of performance program is $\left\{{s}_{1},{s}_{2},\cdots ,{s}_{{m}_{\tau}}\right\}$ .

Step 2.1. Build a robust performance scheduling model RPSM of these accepted programs. Through the method in Section 4.1, we can get the performance sequence of $\left\{{s}_{1},{s}_{2},\cdots ,{s}_{{m}_{\tau}}\right\}$ and the performance cost ${z}_{1}$ .

Step 2.2. Calculate the difference between the revenue of ${m}_{\tau}$ accepted programs

and the penalty cost of $\left(k-{m}_{\tau}\right)$ rejected programs, namely ${z}_{2}={\displaystyle \underset{j=1}{\overset{{m}_{\tau}}{\sum}}{p}_{j}}-{\displaystyle \underset{j={m}_{\tau}+1}{\overset{k}{\sum}}{f}_{j}}$,

and get the performance revenue $z\left({m}_{\tau}\right)={z}_{2}-{z}_{1}$ .

Step 3. $\tau =\tau +1$ . If $\tau \le x$, go back step 2, otherwise go step 4.

Step 4. Denote $\stackrel{\xaf}{m}=\mathrm{arg}\mathrm{max}\left\{z\left({m}_{1}\right),z\left({m}_{2}\right),\cdots ,z\left({m}_{x}\right)\right\}$, select program set $\left\{{s}_{1},{s}_{2},\cdots ,{s}_{\stackrel{\xaf}{m}}\right\}$ to perform, get the maximum performance revenue $z\left(\stackrel{\xaf}{m}\right)$ and the performance sequence.

5. Numerical Experiment

We use Matlab 2017b software for numerical experiments. The experiment was carried out under the Windows 10 Professional 64-bit i5-3230M 8GB RAM operation environment. RPSM2 is a 0 - 1 mixed linear programming model, so it can be solved by the Intlinprog function in Matlab 2017b.

Although RPSM2 contains $\left(n+1\right){m}^{2}+5mn+2m+1$ variables, $n{m}^{2}+8mn+3m-2n$ constraints, after observation, we find that the zero element in the coefficient matrix of the constraints is the majority, which is a sparse matrix [17] . Its density is the ratio of non-zero elements to total elements in the matrix, namely

$\frac{\left(2n+4\right){m}^{2}+\left(18n+3\right)m-4n}{\left[\left(n+1\right){m}^{2}+5mn+2m+1\right]\left[n{m}^{2}+8mn+3m-2n\right]}$ .

And as the values of m and n increase, the density of the RPSM2 coefficient matrix decreases sharply. For sparse matrices, Matlab only stores non-zero element values and their positions. Therefore, we use the sparse feature of RPSM2 coefficient matrix to reduce the variable storage space of the computer and improve the running speed of the procedure.

The parameter setting of the problem P0 is as the following:

1) There are 15 registration programs. These programs require 30 performers to participate in. The relationship between the performers and the programs is a

0 - 1 matrix, defined as $W={\left({w}_{ij}\right)}_{30\times 15}$ . The value of each element is generated by

Matlab according to the random uniform probability.

2) The duration ${q}_{j}$ of program ${s}_{j}$ is a bounded interval value, $\forall j=\left\{1,2,\cdots ,15\right\}$ . The lower bound $l{q}_{j}$ of ${q}_{j}$ obeys the random uniform distribution between the interval [3, 5], and the upper bound $u{q}_{j}$ obeys the random uniform distribution between the interval [6, 10].

3) The unit time waiting cost $l{c}_{i}$ obeys the random uniform distribution between the interval [10, 20], and the unit time appearance fee $u{c}_{i}$ obeys the random uniform distribution between the interval [50, 80].

4) The revenue ${p}_{j}$ of the accepted program ${s}_{j}$ obeys the random uniform distribution between the interval [5000, 20000], and the penalty cost ${f}_{j}$ of the rejected program ${s}_{j}$ obeys the random uniform distribution between the interval [800, 1000].

Due to time and layout restrictions, the decision maker decides to accept 8 to 13 programs, namely $m\in \left\{8,9,10,11,12,13\right\}$ . In order to explain the influence of

decision maker’ risk preference on performance cost, $\mu =0,\frac{m}{2},m$ is selected to

conduct experiments, which represent that decision maker is extremely preferences, moderate risk preferences and very conservative.

The numerical experiment results are shown in Table 3. The first column is the number of the accepted programs, and the second column is the overall performance revenue z, performance cost ${z}_{1}$ and procedure running time when the decision maker prefers the risk. Columns 3 and 4 and so on. (A negative value in the table indicates that the performance revenue is negative).

In Table 3, as the number of accepted programs increases, the variables and constraints in the model increase correspondingly. However, our procedure gets the scheduling scheme within the acceptable time, indicating that the algorithm H has a robust practicality in the actual performance scheduling. The horizontal axis of Figure 1 is the $\mu $ value, indicating risk preference degree of the decision

Table 3. Numerical experiment results.

Figure 1. Relationship between performance cost and risk preference of decision maker.

maker, and the vertical axis is the performance cost
${z}_{1}$ . In the figure,
$m=8$ means that 8 programs are accepted, others and so on. It can be seen from this figure that in the case of determining the set of the accepted programs, the performance cost increases with the increase of

$\mu $, namely if the decision maker avoids the risk, the performance cost will increase and the performance revenue will reduce. Therefore, decision maker can obtain an ideal performance scheduling solution based on their own risk preference.

6. Conclusions

This paper studies a problem of uncertain duration performance scheduling under accepting strategy, the accepted programs will bring in revenue, and the rejected programs will produce corresponding penalty cost. After determining the set of the accepted performance programs, the decision maker’s risk preference coefficient is introduced, and the min-max robust performance scheduling model of these accepted programs is built, and then it is transformed into a 0 - 1 mixed linear programming model to minimize the performance cost. Based on this, we combined with the revenue of the accepted programs and the penalty cost of the rejected programs, the algorithm H for solving the performance scheduling problem under accepting strategy is proposed, which provides a reference for decision maker to choose the ideal program scheduling scheme.

This article does not limit the sequence of program performance, but in the actual world, the decision maker will arrange a program at the opening or finale. Or depending on the type of program, some programs must not be adjacent. In addition, multi-objective functions can also be studied, such as maximizing overall revenue on the basis of ensuring that the performance cost does not exceed the budget.

Cite this paper

Ding, H. , Fan, Y. and Zhong, W. (2018) Robust Optimization of Performance Scheduling Problem under Accepting Strategy.*Open Journal of Optimization*, **7**, 65-78. doi: 10.4236/ojop.2018.74004.

Ding, H. , Fan, Y. and Zhong, W. (2018) Robust Optimization of Performance Scheduling Problem under Accepting Strategy.

References

[1] Cheng, T.C.E., Diamond, J.E. and Lin, B.M.T. (1993) Optimal Scheduling in Film Production to Minimize Talent Hold Cost. Journal of Optimization Theory and Applications, 79, 479-492.

https://doi.org/10.1007/BF00940554

[2] Nordström, A.L. and Tufekci, S. (1994) A Genetic Algorithm for the Talent Scheduling Problem. Computers & Operations Research, 21, 927-940.

https://doi.org/10.1016/0305-0548(94)90021-3

[3] Bomsdorf, F. and Derigs, U. (2008) A Model, Heuristic Procedure and Decision Support System for Solving the Movie Shoot Scheduling Problem. OR Spectrum, 30, 751-772.

https://doi.org/10.1007/s00291-007-0103-6

[4] Garcia De La Banda, M., Stuckey, P.J. and Chu, G. (2011) Solving Talent Scheduling with Dynamic Programming. INFORMS Journal on Computing, 23, 120-137.

https://doi.org/10.1287/ijoc.1090.0378

[5] Wang, S.Y., Chuang, Y.T. and Lin, B.M.T. (2016) Minimizing Talent Cost and Operating Cost in Film Production. Journal of Industrial and Production Engineering, 33, 17-31.

[6] Qin, H., Zhang, Z.Z., Lim, A., et al. (2016) An Enhanced Branch-and-Bound Algorithm for the Talent Scheduling Problem. European Journal of Operational Research, 250, 412-426.

https://doi.org/10.1016/j.ejor.2015.10.002

[7] Sakulsom, N. and Tharmmaphornphilas, W. (2014) Scheduling a Music Performance Problem with Unequal Music Piece Length. Computers & Industrial Engineering, 70, 20-30.

https://doi.org/10.1016/j.cie.2013.12.017

[8] Zhen, L., Liu, B. and Wang, W.C. (2017) Scheduling of Performance with Mixing Piece Length. Journal of Systems Management, 26, 850-856.

[9] Wang, Y. and Tang, J.F. (2016) A Two-Stage Robust Optimization Method for Solving Surgery Scheduling Problem. Journal of Systems Engineering, 31, 431-440.

[10] Xu, X.Q., Cui, W.T., Lin, J. and Qian, Y.J. (2013) Robust Identical Parallel Machines Scheduling Model Based on Min-Max Regret Criterion. Journal of Systems Engineering, 28, 729-737.

[11] Qiu, R.Z. and Huang X.Y. (2011) A Robust Supply Chain Coordination Model Based on Minimax Regret Criterion. Journal of Systems Management, 20, 296-302.

[12] Zhang, L., Chen, T. and Huang, J. (2014) Emergency Network Model and Algorithm Based on Min-Max Regret Robust Optimization. Chinese Journal of Management Science, 22, 131-139.

[13] Bertsimas, D. and Sim, M. (2004) The Price of Robustness. Operations Research, 52, 35-53.

https://doi.org/10.1287/opre.1030.0065

[14] Bertsimas, D. and Sim, M. (2003) Robust Discrete Optimization and Network Flows. Mathematical Programming, 98, 49-71.

https://doi.org/10.1007/s10107-003-0396-4

[15] Charalambous, C. and Conn, A.R. (1978) An Efficient Method to Solve the Minimax Problem Directly. SIAM Journal on Numerical Analysis, 15, 162-187.

https://doi.org/10.1137/0715011

[16] Yao, E.Y., He, Y. and Cheng, S.P. (2001) Mathematical Planning and Combination Optimization. Zhejiang University Press, Hangzhou, 31-35.

[17] Liu, H. (2013) MATLAB R2012a Completely Self-Study. Electronic Industry Press, Beijing.

[1] Cheng, T.C.E., Diamond, J.E. and Lin, B.M.T. (1993) Optimal Scheduling in Film Production to Minimize Talent Hold Cost. Journal of Optimization Theory and Applications, 79, 479-492.

https://doi.org/10.1007/BF00940554

[2] Nordström, A.L. and Tufekci, S. (1994) A Genetic Algorithm for the Talent Scheduling Problem. Computers & Operations Research, 21, 927-940.

https://doi.org/10.1016/0305-0548(94)90021-3

[3] Bomsdorf, F. and Derigs, U. (2008) A Model, Heuristic Procedure and Decision Support System for Solving the Movie Shoot Scheduling Problem. OR Spectrum, 30, 751-772.

https://doi.org/10.1007/s00291-007-0103-6

[4] Garcia De La Banda, M., Stuckey, P.J. and Chu, G. (2011) Solving Talent Scheduling with Dynamic Programming. INFORMS Journal on Computing, 23, 120-137.

https://doi.org/10.1287/ijoc.1090.0378

[5] Wang, S.Y., Chuang, Y.T. and Lin, B.M.T. (2016) Minimizing Talent Cost and Operating Cost in Film Production. Journal of Industrial and Production Engineering, 33, 17-31.

[6] Qin, H., Zhang, Z.Z., Lim, A., et al. (2016) An Enhanced Branch-and-Bound Algorithm for the Talent Scheduling Problem. European Journal of Operational Research, 250, 412-426.

https://doi.org/10.1016/j.ejor.2015.10.002

[7] Sakulsom, N. and Tharmmaphornphilas, W. (2014) Scheduling a Music Performance Problem with Unequal Music Piece Length. Computers & Industrial Engineering, 70, 20-30.

https://doi.org/10.1016/j.cie.2013.12.017

[8] Zhen, L., Liu, B. and Wang, W.C. (2017) Scheduling of Performance with Mixing Piece Length. Journal of Systems Management, 26, 850-856.

[9] Wang, Y. and Tang, J.F. (2016) A Two-Stage Robust Optimization Method for Solving Surgery Scheduling Problem. Journal of Systems Engineering, 31, 431-440.

[10] Xu, X.Q., Cui, W.T., Lin, J. and Qian, Y.J. (2013) Robust Identical Parallel Machines Scheduling Model Based on Min-Max Regret Criterion. Journal of Systems Engineering, 28, 729-737.

[11] Qiu, R.Z. and Huang X.Y. (2011) A Robust Supply Chain Coordination Model Based on Minimax Regret Criterion. Journal of Systems Management, 20, 296-302.

[12] Zhang, L., Chen, T. and Huang, J. (2014) Emergency Network Model and Algorithm Based on Min-Max Regret Robust Optimization. Chinese Journal of Management Science, 22, 131-139.

[13] Bertsimas, D. and Sim, M. (2004) The Price of Robustness. Operations Research, 52, 35-53.

https://doi.org/10.1287/opre.1030.0065

[14] Bertsimas, D. and Sim, M. (2003) Robust Discrete Optimization and Network Flows. Mathematical Programming, 98, 49-71.

https://doi.org/10.1007/s10107-003-0396-4

[15] Charalambous, C. and Conn, A.R. (1978) An Efficient Method to Solve the Minimax Problem Directly. SIAM Journal on Numerical Analysis, 15, 162-187.

https://doi.org/10.1137/0715011

[16] Yao, E.Y., He, Y. and Cheng, S.P. (2001) Mathematical Planning and Combination Optimization. Zhejiang University Press, Hangzhou, 31-35.

[17] Liu, H. (2013) MATLAB R2012a Completely Self-Study. Electronic Industry Press, Beijing.