A Robust Lot Sizing Problem with Ill-known Demands Thanks: The second author of the paper was partially supported by Polish Committee for Scientific Research, grant N N206 492938.

Abstract

The paper deals with a lot sizing problem with ill-known demands modeled by fuzzy intervals whose membership functions are possibility distributions for the values of the uncertain demands. Optimization criteria, in the setting of possibility theory, that lead to choose robust production plans under fuzzy demands are given. Some algorithms for determining optimal robust production plans with respect to the proposed criteria, and for evaluating production plans are provided. Some computational experiments are presented.

Keywords: Fuzzy Optimization, Dynamic Lot Sizing, Uncertain Demand, Possibility Theory

1 Introduction

Nowadays, companies do not compete as independent entities but as a part of collaborative supply chains. Uncertainty in demands creates a risk in a supply chain as backordering, obsolete inventory due to the bullwhip effect [1]. To reduce this risk two different approaches exist that are considered here. The first approach consists in a collaboration between the customer and the supplier and the second one consists in an integration of uncertainty into a planning process.

The collaborative processes mainly aim to reduce a risk in a supply chain [2]. This is done by enforcing a coordination in a supply chain. Two approaches can be applied: vertical and horizontal. The vertical approach is a centralized decision making that synchronizes a supply chain (the most common way to coordinate within companies). The horizontal one refers to the collaborative planning, in which a supply chain can be seen as a chain, where actors are independent entities [3]. The industrial collaborative planning has been standardized for implementing a cooperation between retailers and manufactures. This process is called Collaborative Planning, Forecasting and Replenishment (CPFR) [4]. More precisely, the collaborative processes are usually characterized by a set of point-to-point (customer/supplier) relationships with a partial information sharing [2, 5]. In the collaborative supply chain, a procurement plan is built and propagated through a supply chain. Namely, the procurement plan is composed of three horizons: freezing, flexible and free ones [2]. Quantities in the freezing horizon are crisp and can not be modified, quantities in the flexible horizon are intervals and can be modified under constraints imposed by a previous procurement plan. In the free horizon quantities can be modified without constraints. Another way to reduce a risk in a supply chain is to integrate the uncertainty in a planning process. In the literature, three different sources of uncertainty are distinguished (see [6] for a review): demand, process and supply. These uncertainties are due to difficulties to access to available historical data allowing to determine a probability distribution.

In this paper, we focus on the collaborative supply chain (a supply chain, where actors are independents entities) under uncertain demands. In most companies today, especially in aeronautic companies, actors use the Manufacturing Resource Planning (MRPII) to plan theirs production. MRPII is a planning control process composed of three processes (the production process, the procurement process and the distribution process) and three levels [7]: the strategic level (computing commercial and industrial plans), the tactical level (the Master Production Scheduling (MPS) and the Material Requirement Planning (MRP)) and the operational level (a detailed scheduling and a shop floor control). MRPII have been also extended to take into account: the imprecision on quantities of demands (MPS) [8], the imprecision on quantities of demands and uncertain orders [9] (MRP) and the imprecision on quantities and on dates of demands with uncertain order dates [10] (MRP).

In this paper, we wish to investigate the part of the MRPII process. Namely, the procurement process in the tactical level in the collaborative context. Our purpose is to help the decision maker of a procurement service to evaluate a performance of a given procurement plan with ill-known gross requirements and to compute a procurement plan in a collaborative supply chain (with and without supplier capacity sharing due to a procurement contract) with ill-known gross requirements.

Several production planning problems have been adapted to the case of fuzzy demands: economic order quantity [11, 12], multi-period planning [8, 9, 10, 13, 14, 15, 16, 17], and the problem of supply chain planning (production distribution, centralized supply chain) [18, 19, 20, 21, 22, 23]. In the literature, there are two popular families of approaches for coping with fuzzy parameters. In the first family, a defuzzification is first performed and then deterministic optimization methods are used [20, 21]. In the second one, the objective is expressed in the setting of possibility theory [24] and credibility theory [25]. We can distinguish: the possibilistic programming (a fuzzy mathematical programming) in which a solution optimizing a criterion based on the possibility measure is built [16, 17], the credibility measure based programming in which the credibility measure is used to guaranty a service level (chance constraints on the inventory level) [26] or the goal is to choose a solution that optimizes a criterion based on the credibility measure [13] and a decision support based on the propagation of the uncertainty to the inventory level and backordering level [8, 9, 10]. Here, we restrict our attention to uncertainty propagation in MRP (the tactical level) [8, 9, 10] and we propose methods both for evaluating a procurement plan in terms of costs under uncertain demands and for computing a procurement plan which minimizes the impact of uncertainty on costs, since the approaches proposed in the literature are not able to do this.

Popular setting of problems for hedging against uncertainty of parameters is robust optimization [27]. In the robust optimization setting the uncertainty is modeled by specifying a set of all possible realizations of the parameters called scenarios. No probability distribution in the scenario set is given. The value of each parameter may fall within a given closed interval and the set of scenarios is the Cartesian product of these intervals. Then, in order to choose a solution, two optimization criteria, called the min-max and the min-max regret, can be adopted. Under the min-max criterion, we seek a solution that minimizes the largest cost over all scenarios. Under the min-max regret criterion we wish to find a solution, which minimizes the largest deviation from optimum over all scenarios.

In this paper, we are interested in computing a robust procurement plan (with and without delivering capacity of the supplier sharing). The delivering capacity are composed of two bounds: the lower one being the minimal accepted quantity that is sent to the customer and the upper bound which is due to a production capacity of the supplier. Moreover the customer accepts to have backordering but it is more penalized that inventory. This problem is equivalent to the problem of production planning with backordering, more precisely to a certain version of the lot sizing problem (see, e.g., [28, 29]), where: the procured quantities are production quantities, a production plan; delivering constraints are production constraints, capacity limits on production plans; and the gross requirements are demands. Thus, the problem consists in finding a production plan that fulfills capacity limits and minimizes the total cost of storage and backordering subject to the conditions of satisfying each demand. It is efficiently solvable when the demands are precisely known (see, e.g., [30, 31, 32]). However, the demands are seldom precisely known in advance and the uncertainty must be taken into account.

In this paper, we consider the above problem with uncertain demands modeled by fuzzy intervals. The membership function of a fuzzy interval is a possibility distribution describing, for each value of the demand, the extent to which it is a possible value. In other words, it means that the value of this demand belongs to a $\lambda$-cut of the fuzzy interval with the degree of necessity (confidence) $1-\lambda$. To evaluate a production plan, we assign to it, degrees of possibility and necessity that its cost does not exceed a given threshold and a degree of necessity that costs of the plan fall within a given fuzzy goal. In order to find “robust solutions” under fuzzy demands, we apply two criteria. The first one consists in choosing a production plan which maximizes the degree of necessity (certainty) that its cost does not exceed a given threshold. The second criterion is weaker than the first one and consists in choosing a plan with the maximum degree of necessity that costs of the plan fall within a given fuzzy goal. A similar criterion has been proposed in [33] for discrete optimization problems with fuzzy costs. We provide some methods for finding a robust production plan with respect to the proposed criteria as well as for evaluating a given production plan under fuzzy-valued demands which heavy rely on methods for finding a robust production plan, called optimal robust production plan, in the problem of production planning under interval-valued demands with the robust min-max criterion. Namely, it turns out that the considered fuzzy problems can be reduced to examining a family of the interval problems with the min-max criterion. Therefore, we generalize in this way the min-max criterion under the interval structure of uncertainty to the fuzzy case.

The paper is organized as follows. In Section 2, we recall some notions of possibility theory. In Section 3, we present a lot-size problem with backorders and precise demands. In Section 4, we present our results. Namely, we investigate the interval case, that is the lot-size problem with backorders in which uncertain demands are specified as closed intervals. We construct algorithms for finding an optimal robust production plan (a polynomial algorithm for the case without capacity limits and an iterative algorithm for the case with capacity limits) and for evaluating a given production plan (linear and mixed integer programing methods, a pseudopolynomial algorithm). An experimental evidence of the efficiency of the proposed algorithms is provided. In Section 5, we extend our results from the previous section to the fuzzy case. We study the lot-size problem with backorders with uncertain demands modeled by fuzzy intervals in a setting of possibility theory. We provide methods for seeking a robust production plan with respect to two proposed criteria as well as for evaluating a given production plan under fuzzy-valued demands (the methods heavily rely on the ones from the interval case). The efficiency of the methods is confirmed experimentally.

2 Selected Notions of Possibility Theory

A fuzzy interval $\widetilde{A}$ is a fuzzy set in $\mathbb{R}$ whose membership function $\mu_{\widetilde{A}}$ is normal, quasi concave and upper semicontinuous. Usually, it is assumed that the support of a fuzzy interval is bounded. The main property of a fuzzy interval is the fact that all its $\lambda$-cuts, that is the sets $\widetilde{A}^{[\lambda]}=\{x:\mu_{\widetilde{A}}(x)\geq\lambda\}$, $\lambda\in(0,1]$, are closed intervals. We will assume that $\widetilde{A}^{[0]}$ is the smallest closed set containing the support of $\widetilde{A}$. So, every fuzzy interval $\widetilde{A}$ can be represented as a family of closed intervals $\widetilde{A}^{[\lambda]}=[a^{-[\lambda]},a^{+[\lambda]}]$, parametrized by the value of $\lambda\in[0,1]$. In many applications, the class of trapezoidal fuzzy intervals is used. A trapezoidal fuzzy interval, denoted by a quadruple $\widetilde{A}=(a,b,c,d)$ and its membership function has the following form:

$$\mu_{\widetilde{A}}(z)=\begin{cases}0&\text{if $z\leq a$,}\\ \frac{z-a}{b-a}&\text{if $a<z<b$,}\\ 1&\text{if $b\leq z\leq c$,}\\ \frac{d-z}{d-c}&\text{if $c<z<d$,}\\ 0&\text{if $z\geq d$.}\end{cases}$$

Its $\lambda$-cuts are simply $[a+\lambda(b-a),d-\lambda(d-c)]$ for $\lambda\in[0,1]$. Notice that this representation contains triangular fuzzy intervals ($b=c$).

Let us now recall the possibilistic interpretation of fuzzy intervals. Possibility theory [24] is an approach to handle incomplete information and it relies on two dual measures: possibility and necessity, which express plausibility and certainty of events. Both measures are built from a possibility distribution. Let a fuzzy interval $\widetilde{A}$ be attached with a single-valued variable $a$ (uncertain real quantity). The membership function $\mu_{\widetilde{A}}$ is understood as a possibility distribution, $\pi_{a}=\mu_{\widetilde{A}}$, which describes the set of more or less plausible, mutually exclusive values of the variable $a$. It can encode a family of probability functions [34]. In particular, a degree of possibility can be viewed as the upper bound of a degree of probability [34]. The value of $\pi_{a}(v)$ represents the possibility degree of the assignment $a=v$, i.e. $\mathrm{\Pi}(a=v)=\pi_{a}(v)=\mu_{\widetilde{A}}(v)$, where $\mathrm{\Pi}(a=v)$ is the possibility of the event that $a$ will take the value of $v$. In particular, $\pi_{a}(v)=0$ means that $a=v$ is impossible and $\pi_{a}(v)>0$ means that $a=v$ is plausible. Equivalently, it means that the value of $a$ belongs to a $\lambda$-cut $\widetilde{A}^{[\lambda]}$ with confidence (or degree of necessity) $1-\lambda$. A detailed interpretation of the possibility distribution and some methods of obtaining it from the possessed knowledge are described in [24, 35]. Let $\widetilde{G}$ be a fuzzy interval. Then “$a\in\widetilde{G}$” is a fuzzy event. The possibility of “$a\in\widetilde{G}$”, denoted by $\Pi(a\in\widetilde{G})$, is as follows [36]:

$$\Pi(a\in\widetilde{G})=\sup_{v\in\mathbb{R}}\min\{\pi_{a}(v),\mu_{\widetilde{G}}(v)\}.$$ (1)

$\Pi(a\in\widetilde{G})$ evaluates the extent to which “$a\in\widetilde{G}$” is possibly true. The necessity of event “$a\in\widetilde{G}$”, denoted by $\mathrm{N}(a\in\widetilde{G})$, is as follows:

$$\mathrm{N}(a\in\widetilde{G}) =1-\Pi(a\not\in\widetilde{G})=1-\sup_{v\in\mathbb{R}}\min\{\pi_{a}(v),1-\mu_{\widetilde{G}}(v)\}$$ (2)

$$=\inf_{v\in\mathbb{R}}\max\{1-\pi_{a}(v),\mu_{\widetilde{G}}(v)\},$$

where $1-\mu_{\widetilde{G}}$ is the membership function of the complement of the fuzzy set $\widetilde{G}$. $\mathrm{N}(a\in\widetilde{G})$ evaluates the extent to which “$a\in\widetilde{G}$” is certainly true. Observe that if $G$ is a classical set, then $\Pi(a\in G)=\sup_{v\in G}\pi_{a}(v)$ and $\mathrm{N}(a\in G)=1-\sup_{v\notin G}\pi_{a}(v)$.

3 The Deterministic Problem

We are given $T$ periods. For period $t$, $t=1,\ldots,T$, let $d_{t}$ be the demand in period $t$, $d_{t}\geq 0$ (here we assume that the demands are precise), $x_{t}$ the production amount in period $t$, $x_{t}\geq 0$, $l_{t}$, $u_{t}$ the production capacity limits on $x_{t}$. Let $\mathbb{X}\subseteq\mathbb{R}^{T}_{+}$ be the set of feasible production amounts. Two cases are distinguished, the case with no capacity limits $\mathbb{X}=\{(x_{1},\ldots,x_{T})\;:\;x_{t}\geq 0,t=1,\ldots,T\}$ and the one with capacity limits $\mathbb{X}=\{(x_{1},\ldots,x_{T})\;:\;l_{t}\leq x_{t}\leq u_{t},t=1,\ldots,T\}$. Set $\mathbf{D}_{t}=\sum_{i=1}^{t}d_{i}$ and $\mathbf{X}_{t}=\sum_{i=1}^{t}x_{i}$, $\mathbf{D}_{t}$ and $\mathbf{X}_{t}$ stand for the cumulative demand up to period $t$ and the production level up to period $t$, respectively. Obviously, $\mathbf{X}_{t-1}\leq\mathbf{X}_{t}$ and $\mathbf{D}_{t-1}\leq\mathbf{D}_{t}$, $t=2,\ldots,T$. The costs of carrying one unit of inventory from period $t$ to period $t+1$ is given by $c^{I}_{t}\geq 0$ and the costs of backordering one unit from period $t+1$ to period $t$ is given by $c^{B}_{t}\geq 0$. The nonnegative real function $L_{t}(u,v)$ represents either the cost of storing inventory from period $t$ to period $t+1$ or the cost of backordering quantity from period $t+1$ to period $t$, namely $L_{t}(\mathbf{X}_{t},\mathbf{D}_{t})=c^{I}_{t}(\mathbf{X}_{t}-\mathbf{D}_{t})$ if $\mathbf{X}_{t}\geq\mathbf{D}_{t}$; $c^{B}_{t}(\mathbf{D}_{t}-\mathbf{X}_{t})$ otherwise. The function has the form $L_{t}(\mathbf{X}_{t},\mathbf{D}_{t})=\max\{c^{I}_{t}(\mathbf{X}_{t}-\mathbf{D}_{t}),c^{B}_{t}(\mathbf{D}_{t}-\mathbf{X}_{t})\}$.

Our production planning problem with the deterministic (precise) demands consists in finding a feasible production plan $\boldsymbol{x}=(x_{1},\ldots,x_{T})$, $\boldsymbol{x}\in\mathbb{X}$, that minimizes the total cost of storage and backordering subject to the conditions of satisfying each demand, namely

$$\min_{\boldsymbol{x}\in\mathbb{X}}F(\boldsymbol{x})=\min_{\boldsymbol{x}\in\mathbb{X}}\sum_{t=1}^{T}L_{t}(\mathbf{X}_{t},\mathbf{D}_{t}).$$ (3)

Obviously, the problem (3) is a version of the classical dynamic lot-size problem with backorders (see, e.g., [28, 29]). Without loss of generality, we can assume that an initial inventory $I_{0}$ and an initial backorder $B_{0}$ are equal to zero. Otherwise, one can append period $0$ and assign $x_{0}=I_{0}$ and $d_{0}=0$ with zero inventory cost if $I_{0}>0$ or assign $x_{0}=0$ and $d_{0}=B_{0}$ with zero backorder cost if $B_{0}>0$.

In the case with no capacity limits problem (3) has a trivial optimal solution equal to $(d_{1},\ldots,d_{T})$. In the case with capacity limits, (3) can be formulated as the minimum cost flow problem (see, e.g., [37]):

$$\begin{array}[]{lll}\min&\displaystyle\sum_{t=1}^{T}(c^{I}_{t}I_{t}+c^{B}_{t}B_{t})&\\ \text{s.t. }&B_{t}-I_{t}=\sum_{j=1}^{t}(d_{j}-x_{j}),&t=1,\ldots,T,\\ &l_{t}\leq x_{t}\leq u_{t},&t=1,\ldots,T,\\ &B_{t},I_{t}\geq 0,&t=1,\ldots,T.\end{array}$$ (4)

Problem (4) can efficiently solved, for instance, by an algorithm presented in [32] that takes into account a special structure of the underling network.

4 Robust Problem

Assume that demands $d_{t}$, $t=1,\ldots,T$, in problem (3), are only known to belong to intervals $D_{t}=[d^{-}_{t},d^{+}_{t}]$, $d^{-}_{t}\geq 0$. This means that we neither know the exact demands, nor can we set them precisely. We assume that the demands are unrelated to one another and there is no probability distribution in $D_{t}$, $t=1,\ldots,T$. A vector $S=(d_{1},\ldots,d_{T})$, $d_{t}\in D_{t}$, that represents an assignment of demands $d_{t}$ to periods $t$, $t=1,\ldots,T$, is called a scenario. Thus every scenario expresses a realization of the demands. We denote by $\Gamma$ the set of all the scenarios, i.e. $\Gamma=[d^{-}_{1},d^{+}_{1}]\times\cdots\times[d^{-}_{T},d^{+}_{T}]$. Among the scenarios of $\Gamma$ extreme scenarios can be distinguished, that is the ones, which belong to $\{d^{-}_{1},d^{+}_{1}\}\times\cdots\times\{d^{-}_{T},d^{+}_{T}\}$, the set of extreme scenarios is denoted by $\Gamma_{\text{ext}}$. We denote by $S^{+}$ (resp. $S^{-}$) the extreme scenario in which all the demands are set to their upper (resp. lower) bounds. The demand and the cumulative demand in period $t$ under scenario $S$ are denoted by $d_{t}(S)\in D_{t}$ and $\mathbf{D}_{t}(S)$, respectively, $\mathbf{D}_{t}(S)=\sum_{i=1}^{t}d_{i}(S)$. Clearly, for every $S\in\Gamma$ it holds $\mathbf{D}_{t-1}(S)\leq\mathbf{D}_{t}(S)$, $t=2,\ldots,T$, and $\mathbf{D}_{t}(S)\in[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$. The function $L_{t}(\mathbf{X}_{t},\mathbf{D}_{t}(S))=\max\{c^{I}_{t}(\mathbf{X}_{t}-\mathbf{D}_{t}(S)),c^{B}_{t}(\mathbf{D}_{t}(S)-\mathbf{X}_{t})\}$, represents either the cost of storing inventory from period $t$ to period $t+1$ or the cost of backordering quantity from period $t+1$ to period $t$ under scenario $S$. Now $F(\boldsymbol{x},S)$ denotes the total cost of a production plan $\boldsymbol{x}\in\mathbb{X}$ under scenario $S$, i.e. $F(\boldsymbol{x},S)=\sum_{t=1}^{T}L_{t}(\mathbf{X}_{t},\mathbf{D}_{t}(S))$.

In order to choose a robust production plan, one of robust criteria, called the min-max can be adopted (see, e.g. [27]). In the min-max version of problem (3), we seek a feasible production plan with the minimum the worst total cost over all scenarios, that is

$$\textsc{ROB}:\;\min_{\boldsymbol{x}\in\mathbb{X}}A(\boldsymbol{x})=\min_{\boldsymbol{x}\in\mathbb{X}}\max_{S\in\Gamma}F(\boldsymbol{x},S)=\min_{\boldsymbol{x}\in\mathbb{X}}\max_{S\in\Gamma}\sum_{t=1}^{T}L_{t}(\mathbf{X}_{t},\mathbf{D}_{t}(S)).$$

In other words, we wish to find among all production plans the one that minimizes the maximum production plan cost over all scenarios, that minimizes $A(\boldsymbol{x})$, $A(\boldsymbol{x})$ is the maximal cost of production plan $\boldsymbol{x}$. An optimal solution $\boldsymbol{x}^{r}$ to the problem ROB is called optimal robust production plan.

Let $\boldsymbol{x}\in\mathbb{X}$ be a given production plan. A scenario $S^{o}\in\Gamma$ that minimizes the total cost $F(\boldsymbol{x},S)$ of the production plan $\boldsymbol{x}$ is called optimistic scenario. A scenario $S^{w}\in\Gamma$ that maximizes the total cost $F(\boldsymbol{x},S)$ of the production plan $\boldsymbol{x}$ is called the worst case scenario.

4.1 Evaluating Production Plan

In this section, we show how to evaluate a given production plan $\boldsymbol{x}^{*}\in\mathbb{X}$. We first consider the problem of determining the optimal interval, $F_{\boldsymbol{x}^{*}}=[f_{\boldsymbol{x}^{*}}^{-},f_{\boldsymbol{x}^{*}}^{+}]$, containing possible values of costs of the production plan $\boldsymbol{x}^{*}$ which can be rigorously defined as the following optimization problems:

$$f_{\boldsymbol{x}^{*}}^{-}= \min_{S\in\Gamma}F(\boldsymbol{x}^{*},S),$$ (5)

$$f_{\boldsymbol{x}^{*}}^{+}= \max_{S\in\Gamma}F(\boldsymbol{x}^{*},S).$$ (6)

It is easily seen that the problem of computing the optimal lower bound on costs of $\boldsymbol{x}^{*}$ (5) is equivalent to the one of determining an optimistic scenario $S^{o}$ for $\boldsymbol{x}^{*}$, namely $f_{\boldsymbol{x}^{*}}^{-}=F(\boldsymbol{x}^{*},S^{o})=\min_{S\in\Gamma}F(\boldsymbol{x}^{*},S)$. Similarly, the problem of computing the optimal upper bound on costs of $\boldsymbol{x}^{*}$ (6) is equivalent to the problem of determining a worst case scenario $S^{w}$ for $\boldsymbol{x}^{*}$, i.e. $f_{\boldsymbol{x}^{*}}^{+}=A(\boldsymbol{x}^{*})=F(\boldsymbol{x}^{*},S^{w})=\max_{S\in\Gamma}F(\boldsymbol{x}^{*},S)$. Thus

$$F_{\boldsymbol{x}^{*}}=[f_{\boldsymbol{x}^{*}}^{-},f_{\boldsymbol{x}^{*}}^{+}]=[F(\boldsymbol{x}^{*},S^{o}),F(\boldsymbol{x}^{*},S^{w})].$$ (7)

Using the optimal interval $F_{\boldsymbol{x}^{*}}$ of possible values of costs of production plan $\boldsymbol{x}^{*}$ allows us to evaluate possibility and necessity that the cost of the plan does not exceed a given threshold under uncertain demands modeled by intervals. Hence, in order to assert possibility that the cost of the plan does not exceed a given threshold $g$, i.e. to assert whether there exits a scenario $S\in\Gamma$ for which $F(\boldsymbol{x}^{*},S)\leq g$, it suffices to determine an optimistic scenario $S^{o}$, the optimal lower bound $f_{\boldsymbol{x}^{*}}^{-}=F(\boldsymbol{x}^{*},S^{o})$ and evaluate $f_{\boldsymbol{x}^{*}}^{-}\leq g$. If the inequality holds then there exists a scenario; otherwise not. Similarly, evaluating necessity that the cost of the plan does not exceed a given threshold $g$, i.e. asserting whether $F(\boldsymbol{x}^{*},S)\leq g$ for every scenario $S\in\Gamma$, we only need to determine worst case scenario $S^{w}$, the optimal upper bound $f_{\boldsymbol{x}^{*}}^{+}=F(\boldsymbol{x}^{*},S^{w})$ and evaluate $,f_{\boldsymbol{x}^{*}}^{+}\leq g$. Thus, evaluating a production plan boils down to computing its optimistic and worst case scenarios.

Let us consider the problem of computing an optimistic scenario for a given production plan $\boldsymbol{x}^{*}\in\mathbb{X}$, that is the problem (5). Its minimum is attained for some $S\in\Gamma$, since $F(\boldsymbol{x},S)$ is a continuous function on the bounded closed set $\Gamma$. Problem (5) can be formulated as a linear programming problem:

$$\begin{array}[]{llll}f_{\boldsymbol{x}^{*}}^{-}=&\min&\sum_{t=1}^{T}(c^{I}_{t}I_{t}+c^{B}_{t}B_{t})&\\ &\text{s.t. }&B_{t}-I_{t}=\sum_{j=1}^{t}(s_{j}-x^{*}_{j}),&t=1,\ldots,T,\\ &&s_{t}\in[d^{-}_{t},d^{+}_{t}],&t=1,\ldots,T,\\ &&B_{t},I_{t}\geq 0,&t=1,\ldots,T.\end{array}$$ (8)

If $s^{o}_{t}$, $B^{o}_{t}$ and $I^{o}_{t}$, $t=1,\ldots,T$, is an optimal solution to problem (8), then $S^{o}=(s^{o}_{1},\ldots,s^{o}_{T})$ is an optimistic scenario for $\boldsymbol{x}^{*}$ (an optimistic realization of uncertain demands) and $I^{o}_{t}$ is storing inventory amount from period $t$ to period $t+1$ and $B^{o}_{t}$ represents backordering amount from period $t+1$ to period $t$ under the optimistic scenario $S^{o}$. The problem (8) can be reduced to the classical minimum cost flow problem and effectively solved by algorithms that take into a special structure of the underling network (see, e.g., [32]). Hence, and fact that $c^{I}_{t},c^{B}_{t}\geq 0$ it follows that for $t=1,\ldots,T$ one of $I^{o}_{t}$ and $B^{o}_{t}$ is zero.

Let us study the problem of computing a worst case scenario for a given production plan $\boldsymbol{x}^{*}\in\mathbb{X}$, that is the problem (6). Since $F(\boldsymbol{x}^{*},S)$ is a continuous function on the bounded closed set $\Gamma$, it attains maximum for some $S\in\Gamma$. The problem (6) can be formulated as a mixed integer programming problem (MIP):

$$\begin{array}[]{llll}f_{\boldsymbol{x}^{*}}^{+}=&\max&\sum_{t=1}^{T}(c^{I}_{t}I_{t}+c^{B}_{t}B_{t})&\\ &\text{s.t. }&B_{t}-I_{t}=\sum_{j=1}^{t}(s_{j}-x^{*}_{j}),&t=1,\ldots,T,\\ &&s_{t}\in[d^{-}_{t},d^{+}_{t}],&t=1,\ldots,T,\\ &&B_{t},I_{t}\geq 0,&t=1,\ldots,T,\\ &&I_{t}\leq(1-\delta_{t})\sum_{j=1}^{t}(x^{*}_{j}-d^{-}_{j}),&t=1,\ldots,T,\\ &&B_{t}\leq\delta_{t}\sum_{j=1}^{t}(d^{+}_{j}-x^{*}_{j}),&t=1,\ldots,T,\\ &&\delta_{t}\in\{0,1\},&t=1,\ldots,T.\end{array}$$ (9)

Let $s^{w}_{t}$, $B^{w}_{t}$, $I^{w}_{t}$ and $\delta_{t}$, $t=1,\ldots,T$, be an optimal solution to problem (9). Then $S^{w}=(s^{w}_{1},\ldots,s^{w}_{T})$ is a worst case scenario for $\boldsymbol{x}^{*}$ (a pessimistic realization of uncertain demands) and $I^{w}_{t}$ is storing inventory amount from period $t$ to period $t+1$ and $B^{w}_{t}$ is backordering amount from period $t+1$ to period $t$ under the worst case scenario $S^{w}$. The last two constraints model (9) and the binary variables $\delta_{t}$ ensure that storing inventory from period $t$ to period $t+1$ and backordering from period $t+1$ to period $t$ is not performed simultaneously (either $I_{t}>0$ or $B_{t}>0$). If $\delta_{t}=1$ then backordering is performed $B_{t}>0$; otherwise storing inventory is performed $I_{t}>0$. Thus, the problem (6) turns out to be much harder than (5).

We now solve the problem (6) by means of dynamic programming. Let us present a result which shows that determining a worst case scenario $S^{w}$ can be restricted to the vertices of $\Gamma$, that is to the set of extreme scenarios $\Gamma_{\text{ext}}$. We prove the convexity of the cost function on $\Gamma$.

Proposition 1.

Function $F(\boldsymbol{x}^{*},S)$ is convex on $\Gamma$ for any fixed production plan $\boldsymbol{x}^{*}\in\mathbb{X}$.

Proof.

Function $c^{I}_{t}(\mathbf{X}^{*}_{t}-\mathbf{D}_{t}(S))$ and $c^{B}_{t}(\mathbf{D}_{t}(S)-\mathbf{X}^{*}_{t})$ are convex on $\Gamma$. Then so are $\max\{c^{I}_{t}(\mathbf{X}^{*}_{t}-\mathbf{D}_{t}(S)),c^{B}_{t}(\mathbf{D}_{t}(S)-\mathbf{X}^{*}_{t})\}$ and $\sum_{t=1}^{T}\max\{c^{I}_{t}(\mathbf{X}^{*}_{t}-\mathbf{D}_{t}(S)),c^{B}_{t}(\mathbf{D}_{t}(S)-\mathbf{X}^{*}_{t})\}$. ∎

The following result allows us to reduce the set of scenarios $\Gamma$ to the set of extreme scenarios $\Gamma_{\text{ext}}$.

Proposition 2.

An optimal scenario for problem (6) (a worst case scenario) is an extreme one.

Proof.

Function $F(\boldsymbol{x}^{*},S)$ attains its maximum in $\Gamma$. Since $F(\boldsymbol{x}^{*},S)$ is convex (Proposition 1) and $\Gamma$ is the hyper-rectangle, an optimal scenario for problem (6) is an extreme one (see, e.g., [38]). ∎

Applying Proposition 2, we can rewrite problem (6) as:

$$f_{\boldsymbol{x}^{*}}^{+}=A(\boldsymbol{x}^{*})=F(\boldsymbol{x}^{*},S^{w})=\max_{S\in\Gamma_{\text{ext}}}F(\boldsymbol{x}^{*},S).$$ (10)

We are now ready to give a dynamic programming based algorithm for solving problem (10). Let $\mathbb{D}_{t}$ be the set of feasible cumulative demand levels in period $t$, $t=1,\ldots,T$, i.e. $\mathbb{D}_{t}=\{\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{-})+1,\ldots,\mathbf{D}_{t}(S^{+})\}$, let $\mathcal{L}_{t-1}(\mathbf{D})$ be the maximal cost of a given production plan $\boldsymbol{x}^{*}$ over periods $t,\ldots,T$, when the cumulative demand level up to period $t-1$ is equal to $\mathbf{D}$, $\mathbf{D}\in\mathbb{D}_{t-1}$, $\mathcal{L}_{t-1}:\mathbb{D}_{t-1}\rightarrow\mathbb{R}_{+}$. Set $\mathbb{D}_{0}=\{0\}$. It is evident that:

$$\mathcal{L}_{T}(\mathbf{D}) =0 \mathbf{D}\in\mathbb{D}_{T},$$ (11)

$$\mathcal{L}_{t-1}(\mathbf{D}) =\max\left\{\begin{array}[]{l}L_{t}(\mathbf{X}^{*}_{t},\mathbf{D}+d^{-}_{t})+\mathcal{L}_{t}(\mathbf{D}+d^{-}_{t})\\ L_{t}(\mathbf{X}^{*}_{t},\mathbf{D}+d^{+}_{t})+\mathcal{L}_{t}(\mathbf{D}+d^{+}_{t})\end{array}\right\} \mathbf{D}\in\mathbb{D}_{t-1},$$

$$t=T,\ldots,1.$$

The maximal cost of production plan $\boldsymbol{x}^{*}$ over period $1,\ldots,T$ is equal to $\mathcal{L}_{0}(0)$, $\mathcal{L}_{0}(0)=f_{\boldsymbol{x}^{*}}^{+}$, which is computed according to the backward recursion (11) and (4.1). The corresponding to $\boldsymbol{x}^{*}$ worst case scenario $S^{w}$ can be determined by a forward recursion technique. It is sufficient to store for each $\mathbf{D}\in\mathbb{D}_{t-1}$ the value for which the maximum in (4.1) is attained, that is either $\mathbf{D}+d^{-}_{t}$ or $\mathbf{D}+d^{+}_{t}$. The running time of the dynamic programming based algorithm is $O(T\cdot\mathbf{D}_{T})$, which is pseudo-polynomial. We have thus proved the following theorem.

Theorem 1.

There is an algorithm for computing the maximal cost of a given production plan $\boldsymbol{x}^{*}$ and its a worst case scenario $S^{w}$, which runs in $O(T\cdot\mathbf{D}_{T})$.

Since finding a worst case scenario requires taking into account only extreme demand scenarios (see Proposition 2), the running time of the above algorithm may be additionally refined by reducing the cardinalities of sets $\mathbb{D}_{t}$, $t=1,\ldots,T$, in (11) and (4.1). Note that we need only consider cumulative demand levels $\mathbf{D}$ which can be obtained by summing instantiated demands, at their lower or upper bounds, in the periods up to $t$. Namely, $\mathbf{D}\in\mathbb{D}_{t}$ if and only if $\mathbf{D}=\sum_{k=1}^{t}d_{k}$, where $d_{k}\in\{d^{-}_{k},d^{+}_{k}\}$. Hence, each reduced set of possible cumulative demand levels $\mathbb{D}_{t}$ has form $\{\mathbf{D}^{1}_{t},\ldots,\mathbf{D}^{l}_{t}\}$. Of course, $l$ is bound by $\mathbf{D}_{t}(S^{+})$. An idea of the improved dynamic programming based algorithm can be outlined as follows:

Step 1 (reducing sets $\mathbb{D}_{t}$): built a directed acyclic network composed of node that is associated with $\mathbb{D}_{0}=\{0\}$ and $T$ layers that are associated with the reduced sets of possible cumulative demand levels $\mathbb{D}_{t}=\{\mathbf{D}^{1}_{t},\ldots,\mathbf{D}^{l}_{t}\}$, $t=1,\ldots,T$. The $t$-th layer is composed of nodes corresponding to cumulative demand levels $\mathbf{D}$, $\mathbf{D}\in\mathbb{D}_{t}$. The arc between $\mathbf{D}^{u}_{t-1}$ and $\mathbf{D}^{v}_{t}$ exists if and only if either $\mathbf{D}^{v}_{t}=\mathbf{D}^{u}_{t-1}+d^{-}_{t}$ or $\mathbf{D}^{v}_{t}=\mathbf{D}^{u}_{t-1}+d^{+}_{t}$. Observe that, except for nodes from the last layer, each node in the network has exactly two outgoing arcs. An example of the constructed network is presented in Figure 1.

Step 2: compute the maximal cost of production plan $\boldsymbol{x}^{*}$ over period $1,\ldots,T$ according to the backward recursion (11) and (4.1) in the constructed network with the layers corresponding to reduced sets $\mathbb{D}_{t}$ and store for each $\mathbf{D}\in\mathbb{D}_{t-1}$ the value $\mathbf{D}\in\mathbb{D}_{t}$ for which the maximum in (4.1) is attained.

Step 3: determine a worst case scenario for production plan $\boldsymbol{x}^{*}$ by performing a simple forward recursion in the constructed network using the stored (in Step 2) values for which the maxima in (4.1) are attained.

Figure 1: An example of the constructed network in Step 1.

The network in Step 1 can be built in $O(T\cdot\max_{t=1,\ldots,T}|\mathbb{D}_{t}|)$ time. The running time of Step 2 is the same time as Step 1. Step 3 can be done in $O(T)$ time. Hence, the overall running time of the improved algorithm is $O(T\cdot\max_{t=1,\ldots,T}|\mathbb{D}_{t}|)$. It is easily seen that $\max_{t=1,\ldots,T}|\mathbb{D}_{t}|$ is upper bounded by $\mathbf{D}_{T}(S^{+})$ and at the worst case $\max_{t=1,\ldots,T}|\mathbb{D}_{t}|=\mathbf{D}_{T}(S^{+})$.

Furthermore, the running time can be reduced if $d^{+}_{1}-d^{-}_{1}=\cdots=d^{+}_{T}-d^{-}_{T}=h$. Then we find $\mathbb{D}_{t}=\{\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{-})+h,\ldots,\mathbf{D}_{t}(S^{-})+th\}$, $t=1,\ldots,T$. Now the running time is $O(T^{2})$, which is polynomial.

4.2 Solving the Robust Problem

Let us consider the problem ROB with no capacity limits, i.e. the problem with the set $\mathbb{X}=\{(x_{1},\ldots,x_{T})\;:\;x_{t}\geq 0,t=1,\ldots,T\}$. In this case, we make the assumption: the costs of carrying one unit of inventory from period $t$ to period $t+1$ for every $t=1,\ldots,T$ are equal, we denote it by $c^{I}$ and the costs of backordering one unit from period $t+1$ to period $t$ for every $t=1,\ldots,T$ are equal, we denote it by $c^{B}$. Note that function $F(\boldsymbol{x},S)$ is continuous on $\mathbb{X}$ and $\Gamma$, $\Gamma$ is a closed bounded set, and so $A(\boldsymbol{x})$ is well defined continuous function on $\mathbb{X}$ (see, e.g., [39, Theorem 1.4]). We show that an optimal robust production plan $\hat{\boldsymbol{x}}=(\hat{x}_{1},\ldots,\hat{x}_{T})$, i.e. $\hat{\boldsymbol{x}}=\text{arg }\min_{\boldsymbol{x}\in\mathbb{X}}A(\boldsymbol{x})$, exists and can be computed by the following formulae:

$$\begin{array}[]{lll}\hat{\mathbf{X}}_{1}:=\frac{c^{B}\mathbf{D}_{1}(S^{+})+c^{I}\mathbf{D}_{1}(S^{-})}{c^{B}+c^{I}},&\hat{x}_{1}:=\hat{\mathbf{X}}_{1}&\\ \hat{\mathbf{X}}_{t}:=\frac{c^{B}\mathbf{D}_{t}(S^{+})+c^{I}\mathbf{D}_{t}(S^{-})}{c^{B}+c^{I}},&\hat{x}_{t}:=\hat{\mathbf{X}}_{t}-\hat{\mathbf{X}}_{t-1},&t\geq 2.\end{array}$$ (14)

An algorithm for determining a production plan $\hat{\boldsymbol{x}}$ according to (14) can be implemented in $O(T)$ time. Before we show that $\hat{\boldsymbol{x}}$ is optimal to problem ROB (an optimal robust production plan) with no capacity limits we prove the following proposition.

Proposition 3.

Let $\hat{\boldsymbol{x}}$ be a plan computed according to (14). Then $\hat{\boldsymbol{x}}$ is feasible and $S^{-}$ and $S^{+}$ are the worst case scenarios for plan $\hat{\boldsymbol{x}}$, i.e. $A(\hat{\boldsymbol{x}})=F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})$.

Proof.

See Appendix A. ∎

We are now ready to prove that $\hat{\boldsymbol{x}}$ is an optimal robust production plan.

Theorem 2.

A production plan determined by formulae (14) is an optimal one for problem ROB with no capacity limits.

Proof.

See Appendix A. ∎

Note that if an initial backorder $B_{0}$ or an initial inventory $I_{0}$ are not equal to zero then one can modify the interval demand $D_{1}$ as follows: $D_{1}:=[d^{-}_{1}+B_{0},d^{+}_{1}+B_{0}]$ if $B_{0}>0$ or $D_{1}:=[d^{-}_{1}-I_{0},d^{+}_{1}-I_{0}]$ if $0<I_{0}\leq d^{-}_{1}$, and apply formulae (14) to determine an optimal robust production plan for problem ROB with the modified demand. If $I_{0}>d^{-}_{1}$ then one appends period $0$, as it has been described in Section 3, and applies an algorithm (Algorithm 1) for the case with capacity limits, $l_{t}=0$, $u_{t}=M$, $t=1,\ldots,T$, where $M$ is a large number.

Let us turn to the problem ROB with capacity limits, i.e. the problem with the set $\mathbb{X}=\{(x_{1},\ldots,x_{T})\;:\;l_{t}\leq x_{t}\leq u_{t},t=1,\ldots,T\}$. Notice $\Gamma$ is a bounded closed set. Function $F(\boldsymbol{x},S)$ is continuous on $\mathbb{X}$ and $\Gamma$ and hence $A(\boldsymbol{x})$ is continuous function on $\mathbb{X}$ (see, e.g., [39, Theorem 1.4]). From this and the fact $\mathbb{X}$ is a bounded closed set it follows that $A(\boldsymbol{x})$ attains its minimum on $\mathbb{X}$.

We now construct an iterative algorithm for solving problem ROB based on on iterative relaxation scheme for min-max problems proposed in [40]. Similar methods were developed for min-max regret linear programming problems with an interval objective function [41, 42]. Let us consider the problem (RX-ROB) being a relaxation of problem ROB that consists in replacing a given scenario set $\Gamma$ with a discrete scenario set $\Gamma_{\text{dis}}=\{S^{1},\ldots,S^{K}\}$, $\Gamma_{\text{dis}}\subseteq\Gamma$:

$$\begin{array}[]{rll}\textsc{RX-ROB:}&a^{*}=\min a&\\ \text{s.t. }&a\geq F(\boldsymbol{x},S^{k})&\forall S^{k}\in\Gamma_{\text{dis}},\\ &\boldsymbol{x}\in\mathbb{X},&\end{array}$$ (15)

where $S^{k}=(s^{k}_{t})_{t=1}^{T}$. The constraint $a\geq F(\boldsymbol{x},S^{k})$, called scenario cut, is associated with exactly one scenario $S^{k}\in\Gamma_{\text{dis}}$. Since $\Gamma_{\text{dis}}\subseteq\Gamma$, the maximal cost $a^{*}$ of an optimal solution $\boldsymbol{x}^{*}$ of problem RX-ROB over discrete scenario set $\Gamma_{\text{dis}}$ is a lower bound on the maximal cost of an optimal robust production plan $\boldsymbol{x}^{r}$ for problem ROB, i.e. $a^{*}\leq A(\boldsymbol{x}^{r})$. Note that the scenario cut, $a\geq F(\boldsymbol{x},S^{k})$, associated with $S^{k}$ is not a linear constraint. One can linearize the cut by replacing it in RX-ROB with the following $T+1$ constraints and $2T$ new decision variables:

$$\begin{array}[]{ll}a\geq\sum_{t=1}^{T}(c^{I}_{t}I^{S^{k}}_{t}+c^{B}_{t}B^{S^{k}}_{t}),&\\ B^{S^{k}}_{t}-I^{S^{k}}_{t}=\sum_{j=1}^{t}(s^{k}_{j}-x_{j}),&t=1,\ldots,T,\\ B^{S^{k}}_{t},I^{S^{k}}_{t}\geq 0,&t=1,\ldots,T.\end{array}$$

Our algorithm (Algorithm 1) starts with zero lower bound on the maximal cost of an optimal robust production plan $\boldsymbol{x}^{r}$, $LB=0$, a candidate $\boldsymbol{x}^{*}\in\mathbb{X}$ for an optimal solution for ROB and empty discrete scenario set, $\Gamma_{\text{dis}}=\emptyset$. At each iteration, a worst case scenario $S^{w}$ for $\boldsymbol{x}^{*}$ is computed by applying the method (9) or the dynamic programming based algorithm presented in Section 4.1. Clearly, $A(\boldsymbol{x}^{*})=F(\boldsymbol{x}^{*},S^{w})$ is an upper bound on $A(\boldsymbol{x}^{r})$, $A(\boldsymbol{x}^{r})\leq A(\boldsymbol{x}^{*})$. If a termination criterion is fulfilled (usually $(A(\boldsymbol{x}^{*})-LB)/LB\leq\epsilon$ if $LB>1$; $A(\boldsymbol{x}^{*})-LB\leq\epsilon$ otherwise, $\epsilon>0$ is a given tolerance) then algorithm stops with production plan $\boldsymbol{x}^{*}$, which is an approximation of an optimal robust production plan. Otherwise the worst case scenario $S^{w}=(s^{w}_{t})_{t=1}^{T}$ is added to $\Gamma_{\text{dis}}$, the corresponding to $S^{w}$ scenario cut is appended to problem RX-ROB. Next the updated linear programming problem RX-ROB is solved to obtain a better candidate $\boldsymbol{x}^{*}$ for an optimal solution for ROB and new lower bound $LB=a^{*}$. Since set $\Gamma_{\text{dis}}$ is updated during the course of the algorithm, the computed values of lower bounds are nondecreasing sequence of their values. Then new iteration is started.

Input: Interval demands $D_{t}=[d^{-}_{t},d^{+}_{t}]$, costs $c^{I}_{t}$, $c^{B}_{t}$, $t=1,\ldots T$, initial production plan $\boldsymbol{x}^{*}\in\mathbb{X}$, a convergence tolerance parameter $\epsilon>0$.

Output: A production plan $\hat{\boldsymbol{x}}^{r}$, an approximation of an optimal robust production plan, and its worst case scenario $S^{w}$.

Step 0. $k:=0$, $LB:=0$, $\Gamma_{\text{dis}}:=\emptyset$.
Step 1. $\boldsymbol{x}^{k}:=\boldsymbol{x}^{*}$.
Step 2. Compute a worst case scenario $S^{w}$ for $\boldsymbol{x}^{k}$ by applying the method (9) or the dynamic programming based algorithm presented in Section 4.1.
Step 3. $\Delta:=F(\boldsymbol{x}^{k},S^{w})-LB$. If $LB>1$ then $\Delta:=\Delta/LB$.

If $\Delta\leq\epsilon$ then output $\boldsymbol{x}^{k},S^{w}$ and STOP.
Step 4. $k:=k+1$.
Step 5. $S^{k}:=S^{w}$, $\Gamma_{\text{dis}}:=\Gamma_{\text{dis}}\cup\{S^{k}\}$ and append scenario cut $a\geq F(\boldsymbol{x},S^{k})$ to problem RX-ROB.
Step 6. Compute an optimal solution $(\boldsymbol{x}^{*},a^{*})$ for RX-ROB, $LB:=a^{*}$, and go to Step 1.

Algorithm 1 Solving problem ROB.

In order to choose a good initial production plan $\boldsymbol{x}^{*}\in\mathbb{X}$ in Algorithm 1, we suggest to solve the classical production planning problem (3) with capacity limits (the model (4)) under the midpoint demand scenario $S^{m}$, i.e. $d_{t}(S^{m})=(d^{-}_{t}+d^{+}_{t})/2$, $t=1,\ldots,T$ and take an optimal production plan under the midpoint scenario as an initial production plan.

Theorem 3.

Algorithm 1 terminates in a finite number of steps for any given $\epsilon>0$.

Proof.

See Appendix A. ∎

Note that if an initial inventory $I_{0}$ or an initial backorder $B_{0}$ are not equal to zero then one appends period $0$, as it has been described in Section 3, and applies Algorithm 1 for $T+1$ periods.

Let us illustrate, by the following example, that solving problem ROB leads to a robust production plan. We are given $5$ periods with the production capacity limits on a production plan: $l_{1}=40$, $u_{1}=50$, $l_{2}=30$, $u_{2}=40$, $l_{3}=30$, $u_{3}=40$, $l_{4}=10$, $u_{4}=35$ and $l_{5}=10$, $u_{5}=35$. The costs of carrying one unit of inventory from period $t$ to period $t+1$, $c^{I}_{t}$, for every $t=1,\ldots,5$ equal $1$ and the costs of backordering one unit from period $t+1$ to period $t$, $c^{B}_{t}$, for every $t=1,\ldots,5$ equal $5$. The knowledge about demands in each period is represented by the intervals: $D_{1}=[30,45]$, $D_{2}=[5,15]$, $D_{3}=[10,30]$, $D_{4}=[20,40]$ and $D_{5}=[20,40]$. The scenario set $\Gamma$ (states of the world) is $\Gamma=[30,45]\times[5,15]\times[10,30]\times[20,40]\times[20,40]$ (see Figure 2). The execution of Algorithm 1 gives a production plan: $x^{\text{opt}}_{1}=40$, $x^{\text{opt}}_{2}=30$, $x^{\text{opt}}_{3}=30$, $x^{\text{opt}}_{4}=27.9167$, $x^{\text{opt}}_{5}=10$ with the total cost of $215.833$ ($\boldsymbol{x}^{\text{opt}}$ is an approximation of an optimal robust production plan with convergence tolerance parameter $\epsilon=0.0001$, the maximal cost of $\boldsymbol{x}^{\text{opt}}$ is no more than $0.01\%$ from optimality). The worst case scenario $S^{w}\in\Gamma$ is $d_{1}(S^{w})=30$, $d_{2}(S^{w})=5$, $d_{3}(S^{w})=10$, $d_{4}(S^{w})=20$, $d_{5}(S^{w})=20$ and $\max_{S\in\Gamma}F(\boldsymbol{x}^{\text{opt}},S)=F(\boldsymbol{x}^{\text{opt}},S^{w})=215.833$ (see Figure 2). This means that total costs of production plan $\boldsymbol{x}^{\text{opt}}$ do not exceed the value of $215.833$ over the set of scenarios. Moreover, the plan $\boldsymbol{x}^{\text{opt}}$ has the best worst performance, i.e. it minimizes the total cost over the all scenarios. Additionally, making use of the methods presented in Section 4.1, one gets complete information about all possible values of costs of the production plan $\boldsymbol{x}^{\text{opt}}$ over the set of scenarios $\Gamma$, by determining the optimal interval $F_{\boldsymbol{x}^{\text{opt}}}=[f_{\boldsymbol{x}^{\text{opt}}}^{-},f_{\boldsymbol{x}^{\text{opt}}}^{+}]$ that contains these values (see (7)). This interval equals $[40,215.833]$. A popular approach for solving a problem with uncertain parameters modeled by the classical intervals is taking the midpoints of the intervals (the average values of the possible parameter values) and solving the problem with these deterministic parameters. In our example, the midpoint scenario has the form $d_{1}(S^{\text{mid}})=37.5$, $d_{2}(S^{\text{mid}})=10$, $d_{3}(S^{\text{mid}})=20$, $d_{4}(S^{\text{mid}})=30$, $d_{5}(S^{\text{mid}})=30$. An algorithm for the problem with the midpoint demands (see Section 3) returns an optimal production plan: $x^{\text{mid}}_{1}=40$, $x^{\text{mid}}_{2}=30$, $x^{\text{opt}}_{3}=30$, $x^{\text{mid}}_{4}=10$, $x^{\text{mid}}_{5}=17.5$ with the total cost of $70$. But, if scenario $d_{1}(S^{w})=45$, $d_{2}(S^{w})=15$, $d_{3}(S^{w})=30$, $d_{4}(S^{w})=40$, $d_{5}(S^{w})=40$ (a worst case scenario for $\boldsymbol{x}^{\text{mid}}$) occurs then the cost will be equal to $357.5$ ($\max_{S\in\Gamma}F(\boldsymbol{x}^{\text{mid}},S)=357.5$). Note that $357.5\gg 215.833$. Similar situation is for two extreme scenarios $S^{+}$ and $S^{-}$, the scenarios in which all the demands are set to their upper bounds and the lower bounds, respectively. Again running an algorithm for the crisp problem with the demands under scenario $S^{+}$ and $S^{-}$, we obtain optimal solutions: $x^{+}_{1}=45$, $x^{+}_{2}=30$, $x^{+}_{3}=30$, $x^{+}_{4}=30$, $x^{+}_{5}=35$ with the total cost of $35$ under $S^{+}$ and $x^{-}_{1}=40$, $x^{-}_{2}=30$, $x^{-}_{3}=30$, $x^{-}_{4}=10$, $x^{-}_{5}=10$, under $S^{-}$, with the total cost of $180$. It turns out that if scenario $d_{1}(S^{w})=30$, $d_{2}(S^{w})=5$, $d_{3}(S^{w})=10$, $d_{4}(S^{w})=20$, $d_{5}(S^{w})=20$, $S^{w}\in\Gamma$, (a worst case scenario for $\boldsymbol{x}^{+}$) occurs then the cost of  $\boldsymbol{x}^{+}$ will be equal to $270$. Similarly, if $d_{1}(S^{w})=45$, $d_{2}(S^{w})=15$, $d_{3}(S^{w})=30$, $d_{4}(S^{w})=40$, $d_{5}(S^{w})=40$, $S^{w}\in\Gamma$, (a worst case scenario for $\boldsymbol{x}^{-}$) occurs then the cost of  $\boldsymbol{x}^{-}$ will be equal to $395$. The summary of the input and output data of the above illustrative example is shown in Table 1. Accordingly, we have no doubts that the computed plan $\boldsymbol{x}^{\text{opt}}$ with respect to the min-max criterion (problem ROB) is a robust one.

Figure 2: The computed production plan $\boldsymbol{x}^{\text{opt}}$ in the illustrative example.

Table 1: Summary of the data and results of the illustrative example

	data†			results
	capacity		interval	robust		plan		plan		plan
	limits		demands	plan		under $S^{\text{mid}}$		under $S^{+}$		under $S^{-}$
$t$	$l_{t}$	$u_{t}$	$D_{t}$	$x^{\text{opt}}_{t}$	$d_{t}(S^{w})$	$x^{\text{mid}}_{t}$	$d_{t}(S^{w})$	$x^{+}_{t}$	$d_{t}(S^{w})$	$x^{-}_{t}$	$d_{t}(S^{w})$
1	40	50	[30,45]	40	30	40	45	45	30	40	45
2	30	40	[5,15]	30	5	30	15	30	5	30	15
3	30	40	[10,30]	30	10	30	30	30	10	30	30
4	10	35	[20,40]	27.9167	20	10	40	30	20	10	40
5	10	35	[20,40]	10	20	17.5	40	35	20	10	40
${}^{{\dagger}}c^{I}_{t}=1,c^{B}_{t}=5,t=1,\ldots,5$				the worst costs for plans: $\boldsymbol{x}^{\text{opt}}$, $\boldsymbol{x}^{\text{mid}}$, $\boldsymbol{x}^{+}$, $\boldsymbol{x}^{-}$
				$F(\boldsymbol{x}^{\text{opt}},S^{w})$		$F(\boldsymbol{x}^{\text{mid}},S^{w})$		$F(\boldsymbol{x}^{+},S^{w})$		$F(\boldsymbol{x}^{-},S^{w})$
				215.833		357.5		270		395

In order to check the efficiency of Algorithm 1, we performed some computational tests. For every $T=100,200,\ldots,1000$, ten instances of the problem ROB with capacity limits were generated. In every instance, inventory costs were randomly chosen from the set $\{1,2,\ldots,10\}$, backorder costs were randomly chosen from the set $\{20,21,\ldots,50\}$, the demands and the production capacities were randomly generated intervals $[X,Y]$, where $X$ is an integer-valued random variable uniformly distributed in $\{0,1,\ldots,99\}$ and $Y$ is an integer-valued random variable uniformly distributed in $\{100,101,\ldots,199\}$. To solve the generated instances, we used IBM ILOG CPLEX 12.2 library (parallel using up 2 threads) [43] and a computer equipped with Intel Core 2 Duo 2.5 GHz. In Table 2, minimal, average and maximal computation times in seconds, required to find approximations of optimal robust production plans with convergence tolerance parameter $\epsilon=0.0001$, are presented. Thus, the maximal costs of the computed production plans are no more than $0.01\%$ from optimality. All computations finished in a few iterations and about $98\%$ of the total running time was spent on computing worst case scenarios by MIP model (9). We also implemented and ran the improved dynamic programming based algorithm for computing worst case scenarios presented in Section 4.1, but it turned out that solving MIP model for determining worst case scenarios was much faster than computing them by the dynamic programming algorithm. As we can see from the obtained results, Algorithm 1 allows us to solve quite large problems having up to 1000 periods in reasonable time.

Table 2: Minimal, average and maximal computation times in seconds

T	100	200	300	400	500	600	700	800	900	1000
min	0.25	0.63	0.96	2.04	2.64	23.55	9.64	24.54	61.26	41.44
avg	0.39	0.84	5.68	9.18	11.47	61.44	72.20	180.78	282.20	372.58
max	0.65	1.05	13.40	20.47	21.78	147.32	188.36	328.61	573.26	893.01

5 Fuzzy Problem

In this section, we apply a more elaborate approach to model uncertain demands. Namely, the uncertain demands, in problem (3), are modeled by fuzzy intervals $\widetilde{D}_{t}$, $t=1,\ldots,T$. Here, a membership function of $\widetilde{D}_{t}$ is regarded as a possibility distribution for the values of the unknown demand $d_{t}$ (see Section 2). The possibility degree of the assignment $d_{t}=s$ is $\Pi(d_{t}=s)=\pi_{d_{t}}(s)=\mu_{\widetilde{D}_{t}}(s)$. Let $S=(s_{t})_{t=1}^{T}$ be a scenario that represents a state of the world where $d_{t}=s_{t}$, for $t=1,\ldots,T$. It is assumed that the demands are unrelated one to each other. Hence, the possibility distributions associated with the demands induce the following possibility distribution over all scenarios in $S\in\mathbb{R}^{T}$ (see [44]):

$$\pi(S)=\Pi((d_{1}=s_{1})\wedge\cdots\wedge(d_{T}=s_{T}))=\min_{t=1,\ldots,T}\Pi(d_{t}=s_{t})=\min_{t=1,\ldots,T}\mu_{\widetilde{D}_{t}}(s_{t}).$$ (16)

The value of $\pi(S)$ stands for the possibility of the event that scenario $S\in\mathbb{R}^{T}$ will occur. We have thus extended scenario set $\Gamma$ given by the intervals (see Section 4) to the fuzzy case and now $\widetilde{\Gamma}$ is a fuzzy set of scenarios with membership function $\mu_{\widetilde{\Gamma}}(S)=\pi(S)$, $S\in\mathbb{R}^{T}$. We see at once that the $\lambda$-cuts of $\widetilde{\Gamma}$ for every $\lambda\in(0,1]$ fulfill the following equality:

$$\widetilde{\Gamma}^{[\lambda]}=\{S:\;\pi(S)\geq\lambda\}=[d^{-[\lambda]}_{1},d^{+[\lambda]}_{1}]\times\cdots\times[d^{-[\lambda]}_{T},d^{+[\lambda]}_{T}],$$

which is from (16) and the definition of $\lambda$-cut. We also define $\widetilde{\Gamma}^{[0]}=[d^{-[0]}_{1},d^{+[0]}_{1}]\times\cdots\times[d^{-[0]}_{T},d^{+[0]}_{T}]$. Notice that $\widetilde{\Gamma}^{\lambda}$, $\lambda\in[0,1]$, is the classical scenario set containing all scenarios whose possibility of occurrence is not less than $\lambda$.

5.1 Evaluating Production Plan

In order to choose a reasonable production plan under fuzziness, we first show how to evaluate a given production plan $\boldsymbol{x}\in\mathbb{X}$. Notice that a cost of production plan $\boldsymbol{x}$ is unknown quantity, denoted by $f_{\boldsymbol{x}}$, since demands are unknown and modeled by fuzzy intervals in the setting of possibility theory. Thus, the unknown cost $f_{\boldsymbol{x}}$ falls within fuzzy interval $\widetilde{F}_{\boldsymbol{x}}$, called fuzzy cost of plan $\boldsymbol{x}$, whose membership function $\mu_{\widetilde{F}_{\boldsymbol{x}}}$ is a possibility distribution for the values of fuzzy variable $f_{\boldsymbol{x}}$ (unknown cost of $\boldsymbol{x}$), $\pi_{\widetilde{F}_{\boldsymbol{x}}}=\mu_{\widetilde{F}_{\boldsymbol{x}}}$, defined as follows:

$$\mu_{\widetilde{F}_{\boldsymbol{x}}}(v)=\Pi(f_{\boldsymbol{x}}=v)=\sup_{\{S:\;F(\boldsymbol{x},S)=v\}}\pi(S),\;\;v\in\mathbb{R}.$$ (17)

Making use of (17), we can define degrees of possibility and necessity that a cost of a given plan $\boldsymbol{x}\in\mathbb{X}$ does not exceed a given threshold $g$:

$$\Pi(f_{\boldsymbol{x}}\leq g) =\sup_{v\leq g}\pi_{\widetilde{F}_{\boldsymbol{x}}}(v)=\sup_{\{S:\;F(\boldsymbol{x},S)\leq g\}}\pi(S),$$ (18)

$$\mathrm{N}(f_{\boldsymbol{x}}\leq g) =1-\Pi(f_{\boldsymbol{x}}>g)=1-\sup_{v>g}\pi_{\widetilde{F}_{\boldsymbol{x}}}(v)=1-\sup_{\{S:\;F(\boldsymbol{x},S)>g\}}\pi(S).$$ (19)

It is easily seen that $\Pi(f_{\boldsymbol{x}}\leq g)=\lambda$ means that there exists a scenario $S$ such that $\pi(S)=\lambda$ in which the cost of plan $\boldsymbol{x}$ does not exceed threshold $g$, $F(\boldsymbol{x},S)\leq g$. $\mathrm{N}(f_{\boldsymbol{x}}\leq g)=1-\lambda$ means that for all scenarios $S$ such that $\pi(S)>\lambda$, costs of plan $\boldsymbol{x}$ under these scenarios do not exceed $g$.

We now consider the problem of computing the degrees (18) and (19) of a given plan $\boldsymbol{x}$. Write $\widetilde{F}_{\boldsymbol{x}}^{[\lambda]}=[f^{-[\lambda]}_{\boldsymbol{x}},f^{+[\lambda]}_{\boldsymbol{x}}]$. Note that the interval $\widetilde{F}_{\boldsymbol{x}}^{[\lambda]}$ ($\lambda$-cut of the fuzzy cost) is the optimal interval of possible costs of $\boldsymbol{x}$ (see (7)) in problem ROB under interval scenario set $\widetilde{\Gamma}^{[\lambda]}$. Hence there exists a link between the interval and the fuzzy cases:

$$\Pi(f_{\boldsymbol{x}}\leq g) = \sup\{\lambda\in[0,1]:f^{-[\lambda]}_{\boldsymbol{x}}\leq g\},$$ (20)

$$\mathrm{N}(f_{\boldsymbol{x}}\leq g) = 1-\inf\{\lambda\in[0,1]:f^{+[\lambda]}_{\boldsymbol{x}}\leq g\}.$$ (21)

From equations (20) and (21), we obtain methods for computing the degrees. So, in order to compute $\Pi(f_{\boldsymbol{x}}\leq g)$ (resp. $\mathrm{N}(f_{\boldsymbol{x}}\leq g)$) we need to find the largest (resp. smallest) value of $\lambda$ such that there exits a scenario $S\in\widetilde{\Gamma}^{[\lambda]}$ for which $F(\boldsymbol{x},S)\leq g$ (resp. for every scenario $S\in\widetilde{\Gamma}^{[\lambda]}$ inequality $F(\boldsymbol{x},S)\leq g$ holds), which is equivalent to determine an optimistic scenario $S^{o}\in\widetilde{\Gamma}^{[\lambda]}$ (resp. a worst case scenario $S^{w}\in\widetilde{\Gamma}^{[\lambda]}$) for $\boldsymbol{x}$ by solving (8) (resp. (9)) and evaluating $F(\boldsymbol{x},S^{o})\leq g$ (resp. $F(\boldsymbol{x},S^{w})\leq g$). Notice $f^{-[\lambda]}_{\boldsymbol{x}}=F(\boldsymbol{x},S^{o})$ (resp. $f^{+[\lambda]}_{\boldsymbol{x}}=F(\boldsymbol{x},S^{w})$). Since $f^{-[\lambda]}_{\boldsymbol{x}}$ (resp. $f^{+[\lambda]}_{\boldsymbol{x}}$) is nondecreasing (resp. nonincreasing) function of $\lambda$, we can apply a binary search technique on $\lambda\in[0,1]$.

The fuzzy cost of production plan $\boldsymbol{x}$ (the possibility distribution for costs of $\boldsymbol{x}$), $\widetilde{F}_{\boldsymbol{x}}$ can be determined approximately, if necessary, via the use of $\lambda$-cuts. Namely, the optimal intervals of possible costs of $\widetilde{F}^{[\lambda]}_{\boldsymbol{x}}=[f^{-[\lambda]}_{\boldsymbol{x}},f^{+[\lambda]}_{\boldsymbol{x}}]$ under $\widetilde{\Gamma}^{[\lambda]}$ are computed for suitably chosen $\lambda$-cuts. Then fuzzy cost $\widetilde{F}_{\boldsymbol{x}}$ is reconstructed from their $\lambda$-cuts. This approach makes sense since intervals $\widetilde{F}^{[\lambda]}_{\boldsymbol{x}}$ are nested.

5.2 Fuzzy Robust Problem

We now propose two criteria of choosing a robust solution in the fuzzy-valued problem (3).

We are given a threshold $g$, and we would like to find a production plan which maximizes the degree of certainty (necessity) that its cost does not exceed threshold $g$. Thus, we would like to solve the following problem:

$$\max_{\boldsymbol{x}\in\mathbb{X}}\mathrm{N}(f_{\boldsymbol{x}}\leq g).$$ (22)

There are no doubts that an optimal production plan computed according to (22) is a robust one, since with the highest degree of certainty costs of the plan over scenarios will not exceed threshold $g$. By (21), it is easy to check that problem (22) is equivalent to the following mathematical programming problem:

$$\begin{array}[]{ll}\min&\lambda\\ \text{s.t.}&f^{+[\lambda]}_{\boldsymbol{x}}\leq g,\\ &\lambda\in[0,1],\\ &\boldsymbol{x}\in\mathbb{X}.\end{array}$$ (23)

If $\lambda^{*}$ is the optimal objective and $\boldsymbol{x}^{*}$ is an optimal solution for problem (23) then $\mathrm{N}(f_{\boldsymbol{x}^{*}}\leq g)=1-\lambda^{*}$ and $\boldsymbol{x}^{*}$ is an optimal production plan for (22). If (23) is infeasible then $\mathrm{N}(f_{\boldsymbol{x}}\leq g)=0$ for all $\boldsymbol{x}\in\mathbb{X}$.

We now present a more general criterion of choosing a robust production plan than (22). Namely, suppose that a decision maker knows her/his preferences about a cost of a production plan $f_{\boldsymbol{x}}$ and expresses it by a fuzzy goal $\widetilde{G}$, which is a fuzzy interval with a bounded support and a nonincreasing upper semicontinuous membership function $\mu_{\widetilde{G}}:\mathbb{R}\rightarrow[0,1]$ such that $\mu_{\widetilde{G}}(v)=1$ for $v\in[0,g]$. The value of $\mu_{\widetilde{G}}(f_{\boldsymbol{x}})$ is the extent to which cost $f_{\boldsymbol{x}}$ of $\boldsymbol{x}$ satisfies the decision maker. Now the requirement “$f_{\boldsymbol{x}}\leq g$” is replaced softer one, i.e. “$f_{\boldsymbol{x}}\in\widetilde{G}$”. So, by (2) and (17) the necessity that event “$f_{\boldsymbol{x}}\in\widetilde{G}$” holds can be expressed as follows:

$$\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G}) = 1-\Pi(f_{\boldsymbol{x}}\not\in\widetilde{G})$$ (24)

$$= 1-\sup_{v\in\mathbb{R}}\min\{\pi_{f_{\boldsymbol{x}}}(v),1-\mu_{\widetilde{G}}(v)\}$$

$$= 1-\sup_{S}\min\{\pi(S),1-\mu_{\widetilde{G}}(F(\boldsymbol{x},S))\}.$$

Thus, if $\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G})=1-\lambda$ means that for all scenarios $S$ such that $\pi(S)>\lambda$, the degree that costs of plan $\boldsymbol{x}$ fall within fuzzy goal $\widetilde{G}$, is not less than $1-\lambda$. Note that $\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G})$ is more general and weaker than $\mathrm{N}(f_{\boldsymbol{x}}\leq g)$. If $\mu_{\widetilde{G}}(v)=0$ for $v>g$ then they are the same. Moreover, $\mathrm{N}(f_{\boldsymbol{x}}\leq g)\leq\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G})$.

Let us give the second criterion of choosing a robust plan. We are given a fuzzy goal $\widetilde{G}$, and we wish to find a production plan which maximizes the necessity degree that costs of the plan fall within fuzzy goal $\widetilde{G}$. Thus we need to solve the following optimization problem,

$$\max_{\boldsymbol{x}\in\mathbb{X}}\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G}).$$ (25)

We check at once that problem (25) is equivalent to the following mathematical programming problem, which is from (21) and (24):

$$\begin{array}[]{ll}\min&\lambda\\ \text{s.t.}&f^{+[\lambda]}_{\boldsymbol{x}}\leq g^{+[1-\lambda]},\\ &\lambda\in[0,1],\\ &\boldsymbol{x}\in\mathbb{X}.\end{array}$$ (26)

If $(\boldsymbol{x}^{*},\lambda^{*})$ is an optimal solution for problem (26), then $\mathrm{N}(f_{\boldsymbol{x}^{*}}\in\widetilde{G})=1-\lambda^{*}$. If (26) is infeasible then $\mathrm{N}(f_{\boldsymbol{x}}\in\widetilde{G})=0$ for all $\boldsymbol{x}\in\mathbb{X}$.

An algorithm for solving problem (26) (resp. (23)) is based on the standard binary search technique in $[0,1]$ (the interval of possible values of $\lambda$) which follows from the fact that $f^{+[\lambda]}_{\boldsymbol{x}}$ is nonincreasing and $g^{+[1-\lambda]}$ (resp. $g$) is nondecreasing function of $\lambda$. We call the algorithm the binary search based algorithm. To find an optimal $(x^{*},\lambda^{*})$, $x^{*}\in\mathbb{X}$, $\lambda^{*}\in[0,1]$, with a given error tolerance $\xi>0$, we seek at each iteration, for a fixed $\lambda$, a plan $\boldsymbol{x}\in\mathbb{X}$ satisfying $f^{+[\lambda]}_{\boldsymbol{x}}\leq g^{+[1-\lambda]}$ (resp. $f^{+[\lambda]}_{\boldsymbol{x}}\leq g$), which boils down to seeking an optimal robust production plan $\boldsymbol{x}^{r}$ and its worst case scenario under scenario set $\widetilde{\Gamma}^{[\lambda]}$, i.e. to solving problem ROB (see Section 4). Note that $f^{+[\lambda]}_{\boldsymbol{x}}=A(\boldsymbol{x})$ and also that $f^{+[\lambda]}_{\boldsymbol{x}}\leq g^{+[1-\lambda]}$ (resp. $f^{+[\lambda]}_{\boldsymbol{x}}\leq g$) holds for some $\boldsymbol{x}\in\mathbb{X}$ if and only if it holds for an optimal robust production plan under $\widetilde{\Gamma}^{[\lambda]}$. Thus, at each iteration, we use either formulae (14) for the case without capacity limits or Algorithm 1 for the case with capacity limits. If the length of determined interval of possible values of $\lambda$ is less than $\xi$, then an optimal robust production plan $\boldsymbol{x}^{r}$ for a fixed $\lambda$, $(\boldsymbol{x}^{r},\lambda)$, is an approximation of an optimal solution for problem (26) (resp. (23)) with precision $\xi$. The running time of the above algorithm is $O(I(T)\log\xi^{-1})$ time, where $\xi>0$ is a given error tolerance and $I(|T|)$ is time required for finding an optimal robust production plan and its worst case scenario under interval scenario set $\widetilde{\Gamma}^{[\lambda]}$ (the running time of either (14) or Algorithm 1).

We now show, by the following illustrative example, that determining a production plan maximizing the necessity degree that costs of the plan fall within fuzzy goal $\widetilde{G}$ (problem (25)) is a robust one under uncertain demands modeled by fuzzy intervals. We are given $5$ periods with the production capacity limits on a production plan: $l_{1}=40$, $u_{1}=50$, $l_{2}=30$, $u_{2}=40$, $l_{3}=30$, $u_{3}=40$, $l_{4}=10$, $u_{4}=35$ and $l_{5}=10$, $u_{5}=35$ and the same the costs of carrying one unit of inventory from period $t$ to period $t+1$, $c^{I}_{t}=1$, $t=1,\ldots,5$, and the same the costs of backordering one unit from period $t+1$ to period $t$ ,$c^{B}_{t}=5$, $t=1,\ldots,5$. The demand uncertainty in each period is represented by the triangular fuzzy intervals: $\widetilde{D}_{1}=(30,37.5,45)$, $\widetilde{D}_{2}=(5,10,15)$, $\widetilde{D}_{3}=(10,20,30)$, $\widetilde{D}_{4}=(20,30,40)$ and $\widetilde{D}_{5}=(20,30,40)$, regarded as possibility distributions for the values of the unknown demands. The fuzzy set of scenarios has the membership function: $\mu_{\widetilde{\Gamma}}(S)=\pi(S)=\min_{t=1,\ldots,5}\mu_{\widetilde{D}_{t}}(s_{t})$, $S\in\mathbb{R}^{5}$. The fuzzy goal $\widetilde{G}$ is trapezoidal fuzzy interval $\widetilde{G}=(0,0,195.83,215.42)$, where $195.833$ is the maximal cost of an optimal robust production plan for the problem ROB without capacity limits (an ideal supplier) and under the supports (the interval demands) of the fuzzy demands, i.e. $\widetilde{D}^{[0]}_{t}$, $t=1,\ldots 5$. Thus, a production plan with the cost less than $195.833$ is totally accepted and with the cost greater than $215.42$ is not at all accepted. The binary search based algorithm outputs a production plan (with  $\xi=0.01$ and $\epsilon=0.0001$ for Algorithm 1): $x^{\text{opt}}_{1}=40$, $x^{\text{opt}}_{2}=30$, $x^{\text{opt}}_{3}=30$, $x^{\text{opt}}_{4}=25.3776$, $x^{\text{opt}}_{5}=10$, that maximizes the necessity degree that costs of the plan fall within fuzzy goal $\widetilde{G}$, $\mathrm{N}(f_{\boldsymbol{x}^{\text{opt}}}\in\widetilde{G})=1-\lambda=0.883$. This means that for all scenarios $S$ whose possibility of occurrence is greater than $0.117$, $\pi(S)>\lambda=0.117$, the degree of necessity (the degree of certainty) that total costs of plan $\boldsymbol{x}^{\text{opt}}$ fall within fuzzy goal $\widetilde{G}$, is not less than $0.883$. Furthermore, we are sure that the total costs of the plan do not exceed $196.986$ (the total cost at $\lambda$-cut equal to $0.117$) for every scenario $S$ such that $\pi(S)>0.117$. We now apply existing approaches to our example. We first consider methods based on a defuzzification which take into account only one scenario resulting from a defuzzification of fuzzy parameters (demands) – see, e.g., [20, 21]. Applying, for instance, the index of Yager [45], we get crisp demands: $d_{1}(S^{\text{Y}})=37.5$, $d_{2}(S^{\text{Y}})=10$, $d_{3}(S^{\text{Y}})=20$, $d_{4}(S^{\text{Y}})=30$, $d_{5}(S^{\text{Y}})=30$. An algorithm for the crisp dynamic lot-size problem with these demands (see Section 3) returns an optimal production plan: $x^{\text{Y}}_{1}=40$, $x^{\text{Y}}_{2}=30$, $x^{\text{Y}}_{3}=30$, $x^{\text{Y}}_{4}=10$, $x^{\text{Y}}_{5}=17.5$ with the total cost of $70$. However, the cost of $\boldsymbol{x}^{\text{Y}}$ may be even $313.262$ for scenarios $S$ such that $\pi(S)>0.117$, and so $313.262\gg 196.986$. The necessity degree that costs of $\boldsymbol{x}^{\text{Y}}$ fall within $\widetilde{G}$ equals 0.593, $\mathrm{N}(f_{\boldsymbol{x}^{\text{Y}}}\in\widetilde{G})=0.593$, which gives $\mathrm{N}(f_{\boldsymbol{x}^{\text{Y}}}\in\widetilde{G})<\mathrm{N}(f_{\boldsymbol{x}^{\text{opt}}}\in\widetilde{G})$. Let us examine a possibilistic programming (a mathematical programming with fuzzy parameters), where solution concepts are based on the Bellman-Zadeh approach [46] - see, e.g., [16, 17]. In this way, the assertion of the form “$f_{\boldsymbol{x}}\in\widetilde{G}$”, where $f_{\boldsymbol{x}}$ is a cost of production plan $\boldsymbol{x}$, is treated as a fuzzy constraint and values of the membership function $\mu_{\widetilde{G}}$ stand for degrees of satisfaction of this constraint - the fuzzy goal. In other words, the assertion induces $\mathrm{\Pi}(f_{\boldsymbol{x}}\in\widetilde{G})$. The joint possibility distribution generated by the fuzzy goal as well as constraints with fuzzy demands is the minimum of the possibility of satisfaction of the fuzzy goal and the possibility of feasibility of the constraints and thus an optimal production plan is a plan, denoted by $\boldsymbol{x}^{\text{BZ}}$, that maximizes the possibility degree of satisfaction both goal and the constraints. A trivial verification shows that for production plan $\boldsymbol{x}^{\text{Y}}$ with the total cost of $70$, the possibility degree of satisfaction the goal as well as the constraints is equal to 1 in our fuzzy problem under consideration - in problem (3) with triangular fuzzy demands $\widetilde{D}_{t}$, $t=1,\ldots,5$, and so $\boldsymbol{x}^{\text{BZ}}=\boldsymbol{x}^{\text{Y}}$. However, as we have seen above, this plan is not a robust one. The summary of the data and results of the example is given in Table 3.

Table 3: Summary of the data and results of the illustrative example

	data†			results
	capacity		interval	robust	plan	plan
	limits		demands	plan	(index of Yager)	(Bellman-Zadeh)
$t$	$l_{t}$	$u_{t}$	$\widetilde{D}_{t}$	$x^{\text{opt}}_{t}$	$x^{\text{Y}}_{t}$	$x^{\text{BZ}}_{t}$
1	40	50	(30,37.5,45)	40	40	40
2	30	40	(5,10,15)	30	30	30
3	30	40	(10,20,30)	30	30	30
4	10	35	(20,30,40)	25.3776	10	10
5	10	35	(20,30,40)	10	17.5	17.5
${}^{{\dagger}}c^{I}_{t}=1,c^{B}_{t}=5,t=1,\ldots,5$				the worst costs under $S$ such that $\pi(S)>0.117$ for
$\widetilde{G}=(0,0,195.83,215.42)$				$\boldsymbol{x}^{\text{opt}}$	$\boldsymbol{x}^{\text{Y}}$	$\boldsymbol{x}^{\text{BZ}}$
				196.986	313.262	313.262
				$\mathrm{N}(f_{\boldsymbol{x}^{\text{opt}}}\in\widetilde{G})$	$\mathrm{N}(f_{\boldsymbol{x}^{\text{Y}}}\in\widetilde{G})$	$\mathrm{N}(f_{\boldsymbol{x}^{\text{BZ}}}\in\widetilde{G})$
				0.883	0.593	0.593

In order to evaluate the efficiency of the binary search based algorithm, we show some results of computational experiments. For every number of periods $T=100,200,\ldots,1000$, ten instances of the problem (25) with capacity limits were generated. In every instance, inventory costs were randomly chosen from the set $\{1,2,\ldots,10\}$, backorder costs were randomly chosen from the set $\{20,21,\ldots,50\}$, the production capacities were randomly generated intervals $[X,Y]$, where $X$ is an integer-valued random variable uniformly distributed in $\{0,1,\ldots,99\}$ and $Y$ is an integer-valued random variable uniformly distributed in $\{100,101,\ldots,199\}$, the demands are triangular fuzzy intervals with the supports equal to $[0,199]$ and the modal values equal to $Z$, where $Z$ is an integer-valued random variable uniformly distributed in $\{0,1,\ldots,199\}$, the fuzzy goal $\widetilde{G}$ was modeled as a trapezoidal fuzzy interval $\widetilde{G}=(0,0,c,d)$, where $c$ was chosen as the maximal cost of an optimal robust production plan for the problem ROB without capacity limits and under the interval demands being the supports of the generated triangular fuzzy demands. The values of $d$ were equal to $c+\beta$, for $\beta\in\{0.00\cdot c,0.25\cdot c,0.50\cdot c,0.75\cdot c,1.00\cdot c\}$ and the error tolerance $\xi=0.01$. We used IBM ILOG CPLEX 12.2 library (parallel using up 2 threads) [43] and a computer equipped with Intel Core 2 Duo 2.5 GHz to solve the generated instances. In Table 4 average computation times in seconds are presented. As we can see from the obtained results, the binary search based algorithm, which calls Algorithm 1 at each iteration with convergence tolerance parameter $\epsilon=0.0001$, can solve efficiently the problem (25), with capacity limits, having up to 1000 periods.

Table 4: Average computation times in seconds

	$\beta$
T	$0.00\cdot c$	$0.25\cdot c$	$0.50\cdot c$	$0.75\cdot c$	$1.00\cdot c$
100	0.96	0.92	0.88	0.83	0.92
200	2.36	2.29	2.26	2.11	2.10
300	7.22	6.93	6.81	6.76	6.86
400	14.68	14.67	14.55	14.57	14.55
500	27.51	27.06	26.60	26.69	26.64
600	46.73	46.85	46.72	46.38	46.49
700	111.72	109.70	109.41	109.71	109.43
800	223.22	217.24	217.89	217.65	217.85
900	243.51	243.28	243.87	244.75	245.21
1000	417.78	416.86	416.14	414.14	415.80

6 Conclusion

In this paper, we have proposed methods to compute a robust procurement plan in the collaborative supply chain, where the customer uses a version of MRP with ill-known demands to plan a production. This problem is a certain version of the lot sizing problem with ill-known demands modeled by fuzzy intervals, whose membership functions are regarded as possibility distributions for the values of the unknown demands. We have introduced, in this setting, the degrees of possibility and necessity that the cost of a plan does not exceed a given threshold and a degree of necessity that costs of a plan fall within a given fuzzy goal, which allows us to evaluate a given production plan. Moreover, we have provided methods for computing these degrees. For finding robust production plans under fuzzy demands, we have proposed two criteria: the first one consists in choosing a production plan which maximizes the degree of necessity that its cost does not exceed a given threshold, the second criterion is softer than the first one and consists in choosing a plan with the maximum degree of necessity that costs of the plan fall within a given fuzzy goal. We have constructed the algorithms for determining optimal robust production plans with respect to the criteria and confirmed their efficiency experimentally. The criteria are a generalization, to the fuzzy case, of the known from literature the min-max criterion. Consequently, we have shown in the paper that there exists a link between interval uncertainty with the min-max criterion and possibilistic uncertainty with the necessity based criteria. It turns out that the evaluation of a production plan and choosing a plan in the fuzzy-valued problem are not harder than in the interval-valued case. The difficulty of solving the fuzzy problems lies in the interval case, since it is reduced to solving a small number of interval problems. Therefore, we have discussed first the interval-valued case. In this case, we have considered the problem of determining the optimal interval of possible costs of a production plan, which allowed us to evaluate the plan. Determining the optimal bounds of the interval boils down to computing optimistic and worst case scenarios. We have proposed linear programing based method for computing an optimistic scenario and mixed integer programing and dynamic programming methods for computing a worst case scenario. We have also identified a polynomial solvable case. For computing an optimal robust production plan, we have provided a polynomial algorithm and iterative one for the cases: with no capacity limits and with capacity limits, respectively. Then we have extended the methods introduced for the interval-valued problem to the fuzzy-valued one.

There is still an open question concerning the complexity status of computing a worst case scenario of a given production plan. The problem is pseudopolynomially solvable and polynomially solvable under certain assumptions and seems to be a core of most of the problems considered in the paper. These assumptions are nearly realistic and make possible extension of our approach to the case where a procurement plan is given for a family of product. In other words, when the sum of quantities procured has to respect supplier capacity constraints which are computed from a previous procurement plan. This problem is equivalent to the multi-item capacitated lot sizing problem. The fact that the complexity status is still open creates the possibility to find a polynomial algorithm and to extend our approach to the multi-item, multi-level capacitated lot sizing problem without the assumptions. So, it is an interesting topic of further research.

Appendix A Appendix

Proposition 3.

It is easily seen that $\hat{\mathbf{X}}_{t}\in[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$, $c^{I}(\hat{\mathbf{X}}_{t}-\mathbf{D}_{t}(S^{-}))=c^{B}(\mathbf{D}_{t}(S^{+})-\hat{\mathbf{X}}_{t})$, $t=1,\ldots,T$, and $\hat{x}_{1}\geq 0$. Since $\mathbf{D}_{t-1}(S^{-})\leq\mathbf{D}_{t}(S^{-})$ and $\mathbf{D}_{t-1}(S^{+})\leq\mathbf{D}_{t}(S^{+})$, $t\geq 2$, (14) shows that $\hat{x}_{t}\geq 0$, $t\geq 2$. Hence $L_{t}(\hat{\mathbf{X}}_{t},\mathbf{D}_{t}(S^{-}))=c^{I}(\hat{\mathbf{X}}_{t}-\mathbf{D}_{t}(S^{-}))$ and $L_{t}(\hat{\mathbf{X}}_{t},\mathbf{D}_{t}(S^{+}))=c^{B}(\mathbf{D}_{t}(S^{+})-\hat{\mathbf{X}}_{t})$, $t=1,\ldots,T$, and so $F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})$. Let $S$ be any scenario, $S\in\Gamma$. Therefore, $F(\hat{\boldsymbol{x}},S)=\sum_{t=1}^{T}\max\{c^{I}(\hat{\mathbf{X}}_{t}-\mathbf{D}_{t}(S)),c^{B}(\mathbf{D}_{t}(S)-\hat{\mathbf{X}}_{t})\}\leq\sum_{t=1}^{T}\max\{c^{I}(\hat{\mathbf{X}}_{t}-\mathbf{D}_{t}(S^{-})),c^{B}(\mathbf{D}_{t}(S^{+})-\hat{\mathbf{X}}_{t})\}=\sum_{t=1}^{T}L_{t}(\hat{\mathbf{X}}_{t},\mathbf{D}_{t}(S^{-}))=\sum_{t=1}^{T}L_{t}(\hat{\mathbf{X}}_{t},\mathbf{D}_{t}(S^{+}))\\ =F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})$. ∎

Theorem 2.

We show that $A(\boldsymbol{x})\geq A(\hat{\boldsymbol{x}})$ for every $\boldsymbol{x}\in\mathbb{X}$. Consider any $\boldsymbol{x}^{{}^{\prime}}\in\mathbb{X}$. Let us modify $\boldsymbol{x}^{{}^{\prime}}$ in the following way:

$$\mathbf{X}^{{}^{\prime\prime}}_{t}=\begin{cases}\mathbf{D}_{t}(S^{-})&\text{if $\mathbf{X}^{{}^{\prime}}_{t}<\mathbf{D}_{t}(S^{-})$,}\\ \mathbf{X}^{{}^{\prime}}_{t}&\text{if $\mathbf{D}_{t}(S^{-})\leq\mathbf{X}^{{}^{\prime}}_{t}\leq\mathbf{D}_{t}(S^{+})$, $t=1,\ldots,T$,}\\ \mathbf{D}_{t}(S^{+})&\text{if $\mathbf{X}^{{}^{\prime}}_{t}>\mathbf{D}_{t}(S^{+})$.}\end{cases}$$

Now $\mathbf{X}^{{}^{\prime\prime}}_{t}\in[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$. From the feasibility of $\boldsymbol{x}^{{}^{\prime}}$, it follows that $\mathbf{X}^{{}^{\prime}}_{t-1}\leq\mathbf{X}^{{}^{\prime}}_{t}$, $t\geq 2$. Hence and $\mathbf{D}_{t-1}(S^{-})\leq\mathbf{D}_{t}(S^{-})$ and $\mathbf{D}_{t-1}(S^{+})\leq\mathbf{D}_{t}(S^{+})$, $t\geq 2$, we obtain $\mathbf{X}^{{}^{\prime\prime}}_{t-1}\leq\mathbf{X}^{{}^{\prime\prime}}_{t}$, $t\geq 2$ and, in consequence $\boldsymbol{x}^{{}^{\prime\prime}}\in\mathbb{X}$. Furthermore, it is easy to check that $A(\boldsymbol{x}^{{}^{\prime}})\geq A(\boldsymbol{x}^{{}^{\prime\prime}})$. By the definition of $A(\boldsymbol{x}^{{}^{\prime\prime}})$, we have $A(\boldsymbol{x}^{{}^{\prime\prime}})\geq\max\{F(\boldsymbol{x}^{{}^{\prime\prime}},S^{-}),F(\boldsymbol{x}^{{}^{\prime\prime}},S^{+})\}$. We now only need to show that $\max\{F(\boldsymbol{x}^{{}^{\prime\prime}},S^{-}),F(\boldsymbol{x}^{{}^{\prime\prime}},S^{+})\}\geq F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})$.

Let us focus on a production plan $\hat{\boldsymbol{x}}$ (determined by formulae (14)). Notice that function $F(\boldsymbol{x},S^{-})=\sum_{t=1}^{T}c^{I}(\mathbf{X}_{t}-\mathbf{D}_{t}(S^{-}))$ (resp. $F(\boldsymbol{x},S^{+})=\sum_{t=1}^{T}c^{B}(\mathbf{D}_{t}(S^{+})-\mathbf{X}_{t})$) are the sum of linear and increasing (resp. decreasing) functions with respect to $\mathbf{X}_{t}$, $c^{I}(\mathbf{X}_{t}-\mathbf{D}_{t}(S^{-}))\geq 0$ (resp. $c^{B}(\mathbf{D}_{t}(S^{+})-\mathbf{X}_{t})\geq 0$) for $\mathbf{X}_{t}\in[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$, $t=1,\ldots,T$. Hence, for each $t=1,\ldots,T$ there exists the intersection point in interval $[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$. It is easy to check that $\hat{\mathbf{X}}_{t}\in[\mathbf{D}_{t}(S^{-}),\mathbf{D}_{t}(S^{+})]$ and $c^{I}(\hat{\mathbf{X}}_{t}-\mathbf{D}_{t}(S^{-}))=c^{B}(\mathbf{D}_{t}(S^{+})-\hat{\mathbf{X}}_{t})$, $t=1,\ldots,T$. Therefore, the points $\hat{\mathbf{X}}_{t}$, $t=1,\ldots,T$, are the intersection ones and $\max\{F(\boldsymbol{x}^{{}^{\prime\prime}},S^{-}),F(\boldsymbol{x}^{{}^{\prime\prime}},S^{+})\}\geq F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})$. Proposition 3 shows that $F(\hat{\boldsymbol{x}},S^{-})=F(\hat{\boldsymbol{x}},S^{+})=A(\hat{\boldsymbol{x}})$ and we thus get $A(\boldsymbol{x}^{{}^{\prime}})\geq A(\boldsymbol{x}^{{}^{\prime\prime}})\geq A(\hat{\boldsymbol{x}})$. ∎

Theorem 3.

The proof is almost the same as those given in [47, Theorem 2.5] and [40, Theorem 3]. We will denote by $\{(\boldsymbol{x}^{k},a^{k})\}$ the sequence of optimal solution $(\boldsymbol{x}^{*},a^{*})$ computed in consecutive iterations in Step 6 of Algorithm 1, $k$ stands for the $k$-th iteration. By picking a subsequence, if necessary, the sequence $\{(\boldsymbol{x}^{k},a^{k})\}$ converges to point $(\hat{\boldsymbol{x}},\hat{a})$, $\hat{\boldsymbol{x}}\in\mathbb{X}$, which follows from the fact that sequence $\{\boldsymbol{x}^{k}\}$ belong to the bounded and closet set $\mathbb{X}\subseteq\mathbb{R}^{T}_{+}$ ($\mathbb{X}$ is a compact set) and $\{a^{k}\}$ is a nondecreasing sequence bounded above. Similar considerations apply to the sequence $\{S^{k}\}$ of worst case scenarios determined in Step 2 of Algorithm 1. The set $\Gamma$ is a closed and bounded (compact) and hence $\{S^{k}\}$ converges to $\hat{S}\in\Gamma$. Since scenario cuts are appended to problem RX-ROB, the inequality $a^{k+1}\geq F(\boldsymbol{x}^{k+1},S^{k})$ holds. By continuity of $F$, we have

$$\hat{a}\geq F(\hat{\boldsymbol{x}},\hat{S}).$$ (27)

Let us define the set $\mathcal{S}^{w}(\boldsymbol{x})$ of worst case scenarios for $\boldsymbol{x}\in X$, i.e. $\mathcal{S}^{w}(\boldsymbol{x})=\{S^{w}\,|\,S^{w}=\text{arg }\max_{S\in\Gamma}F(\boldsymbol{x},S)\}$, which is the point-to-set mapping. The set $\mathcal{S}^{w}(\boldsymbol{x})$ is nonempty for every for $\boldsymbol{x}\in X$. By [39, Theorem 1.5], $\mathcal{S}^{w}$ is upper semicontinuous at $\hat{\boldsymbol{x}}$ and so $\hat{S}\in\mathcal{S}^{w}(\hat{\boldsymbol{x}})$. Therefore

$$A(\hat{\boldsymbol{x}})=\max_{S\in\Gamma}F(\hat{\boldsymbol{x}},S)=F(\hat{\boldsymbol{x}},\hat{S}).$$ (28)

Combining (27) and (28), we obtain $\hat{a}\geq A(\hat{\boldsymbol{x}})$. By [39, Lemma 1.2] $A$ is upper semicontinuous, which yields

$$A(\boldsymbol{x}^{k})=F(\boldsymbol{x}^{k},S^{k})\leq a^{k}+\epsilon,\text{ for some sufficiently large }k.$$

This implies that the termination criterion in Step 3 of Algorithm 1 will be fulfilled in a finite number of iterations. ∎