Research Article  Open Access
Qi Zeng, XiuZhen Xu, "Assessing the Reliability of a Multistate Logistics Network under the Transportation Cost Constraint", Discrete Dynamics in Nature and Society, vol. 2016, Article ID 2628950, 8 pages, 2016. https://doi.org/10.1155/2016/2628950
Assessing the Reliability of a Multistate Logistics Network under the Transportation Cost Constraint
Abstract
Reliability and transportation cost are two important indicators to measure the performance of logistics network. As a combination of reliability and transportation cost, the performance index of logistics network is defined as the probability that at least units of flow demand can be successfully transmitted from the source to the destination with the total transportation cost less than or equal to . In this paper, an algorithm is developed to calculate in terms of minimal paths (MPs for short). The proposed algorithm employs a decomposition technique to divide the search space of MPs such that the search space of MPs can be dramatically reduced, and thereby MPs can be efficiently obtained. An example is provided to illustrate the proposed algorithm. Finally, computational experiments conducted on one benchmark network indicate that the proposed method has an advantage over the existing methods.
1. Introduction
With the swift development of Ecommerce, express logistics in China is presenting a vigorous development trend. Meanwhile, most of enterprises begin to pay more attention to logistics performance for improving business efficiency. Logistics networks which provide the infrastructure for the storage and distribution of products undertake the mission of promoting the efficient and safe movement of goods over time and space and thus play an increasingly important part in sustaining the national economic and social development. It is widely accepted that an efficient, reliable, and costeffective logistics network not only affects a firm’s shortterm benefits but also influences a firm’s longterm development.
The performance assessment of logistics networks is a popular issue in the field of logistics and supply chain management. In a realworld environment, the performance of logistics networks is always affected by various unexpected events, such that it is subject to degradation. Hence, it is of importance from the perspective of logistics management to assess the capability of logistics networks to ensure the delivery of required quantity of goods to the right demand point. In addition, the transportation cost is also a major concern for logistics providers. Therefore, combining capacity reliability and transportation cost for performance measure of logistics networks, in a sense, is of significance to possible improvement of logistics performance.
Theoretically, logistics network can be represented as sets of both nodes and arcs, where each node stands for a supplier, a transfer center, or a market, and each arc connecting a pair of nodes stands for transportation medium (traffic tools, traffic routes, or both) [1]. In fact, the transportation medium, which is multistate due to the nature of the traffic tools (truck, railway, cargo ship, etc.), may be in a failure state, partial failure state, or maintenance state [1, 2]. That is, the number of available traffic tools is not fixed in some sense, and thus each arc has several possible capacities. Therefore, logistics network can be modeled as a multistate flow network in which arcs are associated with multiple integer capacities, operational reliability, and unit transportation cost. And, the goods transported through such logistics network are looked upon as a flow. The reliability index of logistics network is defined as the probability that at least units of flow demand can be successfully transmitted from the source to the destination with the total transportation cost less than or equal to c [3–9]. One of the general algorithms for computing is using minimal paths (MPs) [4–9]. A MP, x, is a minimal state vector meeting the demand and the cost constraint c, which means that, for any , does not meet the demand or the cost constraint [6]. If all MPs are known, the wellknown InclusionExclusion rule is available to calculate [7–9]. So, the most important work is how to efficiently determine all MPs.
A minimal path (MP) is a subset of arcs, such that if any arc is removed from this set, the remaining set is no longer a path [10]. Based on MPs, Lin [5] proposed a simple algorithm to search for all MPs. Lin’s algorithm first uses the enumeration algorithm to search for MP candidates. Then, all MP candidates are verified whether they are MPs via the comparison method. It is time consuming to check MP candidates by a comparison method due to the exponentially growing number of MP candidates [7]. Lin [6] extended the work to the case with unreliable nodes. The method of Lin uses MPs to assign the flow to each component (arc or node). All MPs can be obtained by the comparison algorithm. Instead of the comparison method, Yeh [7] proposed a cyclechecking method to verify MP candidates and demonstrated that his method is more efficient. Without requiring MPs information, Yeh [8] proposed an algorithm to search for MPs by solving a simple model and proved that his algorithm is more efficient than the MPs based algorithms [5–7]. And yet, the algorithm by Yeh [8] is an exhaustive enumeration method in itself, so its searching efficiency is low. By adding some constraints to the minimal capacities of arcs, Niu and Xu [9] proposed an algorithm to improve the searching efficiency of solving MPs.
The above discussions indicate that the methods in [5–7] for solving MPs require all MPs information. However, it is a cumbersome task to find all MPs, because the MP problem is proven to be NPhard [11]. Furthermore, the time complexity of the MPs based algorithms [5–7] is directly proportional to the number of MPs which grows exponentially with the scale of the network; thus, the computational efficiency of these algorithms in [5–7] is far from satisfactory. The methods in [8, 9] require no MPs information and thus are an improvement to the ones in [5–7]. But, one remarkable limitation of Yeh’s method [8] is that a huge number of state vectors need to be enumerated, which results in the lower searching efficiency. While Niu and Xu [9] have made some efforts to advance the searching efficiency, the method by Niu and Xu is even inferior to the one by Yeh when the demand level is low, and thus its searching efficiency is also not satisfactory. Hence, the purpose of this paper is to propose a new efficient algorithm for solving all MPs. The major contribution of this paper is that we develop a new approach to shorten the search space of MP. In particular, to decrease the number of enumerated state vectors, an effective decomposition technique is employed to divide the search space of MPs, such that the search subspace of MPs can be dramatically reduced; thereby MPs can be efficiently obtained. A real example is provided to illustrate the proposed algorithm. Moreover, the experimental results clearly indicate that the proposed algorithm is superior to the abovementioned methods [5–9] in solving the MP problem.
The rest of this paper is organized as follows. Section 2 introduces the network model and basic results. In Section 3, the proposed algorithm is described in detail, and its time complexity is also analyzed. In Section 4, an illustrative example is provided to demonstrate the proposed algorithm. Computational experiments are conducted in Section 5 to compare the proposed algorithm with the existing methods. The final section presents the concluding remarks.
2. Preliminaries
2.1. Multistate Network Model
A multistate network is represented by with the source node and the destination node t, where is the set of nodes with denoting the number of nodes except and , is the set of arcs with denoting the number of arcs, is the largest capacity vector, and is the cost vector, where is the unit transportation cost of . The state of an arc is denoted by which takes integer values from 0 to , where denotes the largest capacity of arc . A state vector indicates the current state of each arc. Let denote a set of state vectors, and the smallest and largest state vectors in set are denoted by and , respectively. Let denote the maxflow of the network under and denote the flows through for ; then the total transportation cost of the network under is . Specifically, if is a feasible state vector (i.e., it satisfies the flow conservation law), then . The network model satisfies the following assumptions [7–9]: the state of each arc is a random variable which takes integer values from 0 to according to a given distribution; the states of different arcs are statistically independent; all flows in the network obey the conservation law: that is, total flows into and from a node (other than the source and destination nodes) are all equal.
To facilitate the understanding, Figure 1 and Table 1 are used to illustrate several notations. Figure 1 shows , , , and . Table 1 indicates and . Given a state vector , which indicates that the current states of , , , , , and are 2, 1, 1, 0, 1, and 2, respectively, the maxflow of the network under is . Obviously, is a feasible state vector; then the total transportation cost of the network under is 2 3 + 1 1 + 1 1 + 0 1 + 1 1 + 2 3 . Let , , , , , ; then the smallest state vector in is , and the largest state vector in is . The maxflow of the network under and is = 0 and = 2.

2.2. Evaluation of in terms of MPs
The InclusionExclusion method is introduced to show how to calculate in terms of MPs. Assume are all MPs, and let , , where , ), and means that ≥ for [6]; then can be evaluated via the InclusionExclusion method as follows: where = , = , = = .
2.3. Fundamental Results on the MP Problem
A state vector x = (, ) is a MP if and only if ; for each > 0, where : that is, state is 1 for and 0 for other arcs; and [8, 9]. The following lemma which is originated from the work of Niu and Xu [9] is a basis for searching for MPs.
Lemma 1. A state vector is a MP if and only if satisfies the following conditions:where is the set of arcs emanating from node and is the set of arcs pointing to node . It should be pointed out that (2)–(4) are established on the basis of the wellknown flow conservation law [9]. By lemma, it is easy to obtain the following conclusion.
Corollary 2. For a state vector with , if satisfies conditions (3), (4), (5), and (6), then is a MP.
By Lemma 1, the enumeration method can be used to search for MPs. Given a multistate network , (4) shows that the total number of state vectors contained in search space of MPs is which is a pretty huge number. Thus, it is inefficient to search for MPs by Lemma 1. Below, we will introduce how to divide the search space into disjoint subspaces from which all MPs can be efficiently obtained.
3. The Proposed Algorithm
Given a set of state vectors , for , where is the set of states of with respect to , that is, , if we limit the search space of MP in set , it is apparent that Lemma 1 can be transformed into the following conclusion.
Corollary 3. A state vector is a MP in set if and only if satisfies the following conditions:
It should be noted that the only difference between Lemma 1 and Corollary 3 is the state range of . Next, we will consider a special state vector in set X: the smallest state vector . It is clear that there are three cases with respect to Case : ; Case : ; and Case : .
Case 1 (). If > d, we have the following statement.
Theorem 4. For set X, if , there is no MP in .
Proof. Since is the smallest state vector in , we have for any vector . If , then for any . That is, for any ; then, by the definition of MP, there is no MP in .
According to Theorem 4, once we have verified that holds, set will be discarded (no MP exists in ). Thus, the search space of MPs will be reduced, which is extremely beneficial in solving MPs.
Case 2 (). In such case, we can claim that the following theorem holds.
Theorem 5. For set , if , there exists no more than one MP in set . Moreover, if there exists one MP, is the unique MP in set X.
Proof. Suppose vector is a MP in set X; then M(x) = d. Note that , and , which is contrary to the definition of MP. Therefore, only has the chance to be a MP in X.
Theorem 5 provides an effective method to search for MPs in set . If = d holds, we only need to check whether is a MP by Corollary 3. If is a MP, it is the unique MP in set . Otherwise, there is no MP in set X.
Case 3 (). In such case, we need to analyze the largest state vector in set . If , we have the following statement.
Theorem 6. For set X, if , there is no MP in set X.
Proof. Since is the largest state vector in set X, we have for any . If , it is clear to derive for any . That is, for any ; then, by the definition of MP, there is no dMP in set X.
It is noteworthy that all of the three theorems contribute to reducing the search space of MPs. If , then we have . In such case, the enumeration algorithm can be used to search for MPs in set by Corollary 3. Now we consider a special case wherein set satisfies and simultaneously. The following theorem shows that a dflow (,) can be derived from using the maxflow algorithm, where dflow (,) consists of flows through each arc for when the maxflow of the network from to is d.
Theorem 7. For set X, if and , then there exists at least one dflow , such that ,.
Proof. Since , at least units of flow can be transmitted from the source to the destination when the network is under . Also, such a flow can be obtained by the following steps [10, 12, 13]:(a)Add a fictitious node , and add an arc with a fixed state from to .(b)The state of each arc E is defined by .(c)Determine the maxflow from to . Since , and the state of the arc from to is d, the maxflow from to is d.(d)For each E, the flows through constitute a dflow whenever the maxflow from to is determined.Since the flows through are upper bounded by , we have for all . Meanwhile, it is apparent that for all E. Thus, we have for all ; that is, for all . Consequently, we have dflow .
The definition of flow indicates that a dflow satisfies the flow conservation law; thus, it satisfies (7)–(9). Accordingly, for such a dflow derived from by Theorem 7, if there is no directed cycle in it and , it is a MP. Otherwise, it is not a MP. Grounded on the obtained dflow , set can be decomposed into disjoint subsets from which other MPs can be derived. The decomposition technique which is proposed by Jane and Laih [13] will be described below.
Let ; then it is obvious that is nonempty. Suppose , by pivoting on arcs , , one by one; is decomposed into nonempty disjoint subsets , , , and as follows:(1)Pivot on arc (2)Pivot on arc (3)(q)Pivot on arc , , (q + 1) = , for .
Theorem 8. (i) ; (ii) subsets , , , and are nonempty and disjoint; (iii) dflow is the smallest state vector in .
Readers may refer to the work of Jane and Laih [13] for the proofs. Now, we first consider subset . Since the dflow derived from Theorem 7 is the smallest state vector in set , if it is a MP, it is the unique MP in set ; otherwise, there is no MP in set . As a result, there is no need to check other state vectors in ; that is, can be discarded as soon as dflow () has been checked, which will reduce the number of enumerated state vectors to a large extent.
Then, we discuss subsets . For each , if or for subset , there is no MP in by Theorem 6 or Theorem 4. Then is discarded, which will also cut down the search space of MPs. If , we only need to check whether is a MP by Theorem 5. Otherwise, we have , and thus the enumeration algorithm can be employed to search for MPs in by Corollary 3.
Finally, we discuss subset . It is noted that . If , is also discarded. Otherwise, subset satisfies and simultaneously. Consequently, Theorem 7 can be applied to obtain a dflow from . Based on the derived dflow, is decomposed into nonempty disjoint subsets in the same manner. As discussed earlier, all of these subsets can be handled similarly. The decomposition procedure will be terminated when no subset, , satisfies and simultaneously.
Based on the above discussions, all MPs can be found using the following steps.
Step . Let , for .
Step . Compute . If , there is no MP in (by Theorem 6), and halt.
Step . Seek a dflow by Theorem 7 and compute .
Step . For , is decomposed into , , , and by Theorem 8.
Step . For , if there is no directed cycle in it and , it is a MP; otherwise, it is not a MP.
Step . For i = 2 to q,(5.1)compute ; if , there is no MP in (by Theorem 6), and go to (5.4);(5.2)compute ; if , there is no MP in (by Theorem 4); if , check if is a MP by Corollary 3, and go to (5.4);(5.3)use the enumeration algorithm to search for MPs in by Corollary 3;(5.4).
Step . Let , and go to Step .
The time complexity of the proposed method is discussed below. Either Step or Step takes time. Step takes to evaluate [14]. A dflow can be found in time by the maxflow algorithm. It requires time to obtain [13]. Then, Step requires time in total. Step takes time to decompose . Step requires time. Therefore, the algorithm totally takes time from Steps to . It is apparent that the time complexity of Step is determined by the number of enumerated state vectors. The number of state vectors in Step is less than ; thus, the time complexity of Step is less than ). Therefore, the time complexity of the proposed algorithm is less than ).
4. An Illustrative Example
In this section, we use a simple logistics network to demonstrate the proposed algorithm step by step. The network, as shown in Figure 1, is cited from the literature in [6–9]. The data of each arc is given in Table 1. If the manager is assigned a budget of 14 (c = 14) to transport 3 units of commodity (demand level d = 3) from to t, the reliability index R_{(3,14)} can be obtained as follows (for convenience, the set of state vectors is denoted by its smallest and largest state vectors as [, ]).
Step . Let .
Step . .
Step . Seek a 3flow (3, 2, 1, 0, 0, 1), and compute .
Step . X is decomposed into , , , , and , where , , , ,, and .
Step . Since , is not a (3, 14)MP.
Step . .(5.1), and .(5.2).(5.3)Use the enumeration algorithm to search for (3, 14)MPs in , but there is no (3, 14)MP.(5.4).(5.1), and .(5.2).(5.3)Use the enumeration algorithm to search for (3, 14)MPs in , but there is no (3, 14)MP.(5.4).(5.1), and ; then there is no (3, 14)MP in .
Step . Let , and go to Step .
Step . .
Step . Seek a 3flow (2, 2, 0, 0, 1, 1), .
Step . X is decomposed into , , , , and , where = [(0, 0, 0, 0, 0, 0), (1, 2, 1, 1, 2, 2)], = [(2, 0, 0, 0, 0, 0), (2, 1, 1, 1, 2, 2)], = [(2, 2, 0, 0, 0, 0), (2, 2, 1, 1, 0, 2)], = [(2, 2, 0, 0, 1, 0), (2, 2, 1, 1, 2, 0)], and = [(2, 2, 0, 0, 1, 1), (2, 2, 1, 1, 2, 2)].
⫽ the disjoint sets of state vectors are denoted here as , , , , and ⫽
Step . Since there is no directed cycle in (2, 2, 0, 0, 1, 1), and C() = 12 < 14, = (2, 2, 0, 0, 1, 1) is a (3, 14)MP.
Step . .(5.1), and M() = 3.(5.2).(5.3)Use the enumeration algorithm to search for (3, 14)MPs in , but there is no (3, 14)MP.(5.4).(5.1) = [(2, 2, 0, 0, 0, 0), (2, 2, 1, 1, 0, 2)], and , then there is no (3, 14)MP in , and go to (5.4).(5.4).(5.1), and ; then there is no (3, 14)MP in .
Step . Let X = = [(0, 0, 0, 0, 0, 0), (1, 2, 1, 1, 2, 2)], and go to Step .
Step . .
Step . Seek a 3flow (1, 1, 0, 0, 2, 2), .
Step . is decomposed into , , , , and , where , , , , and .
⫽the disjoint sets of state vectors are denoted here as , , , , and ⫽
Step . Since there is no directed cycle in (1, 1, 0, 0, 2, 2), and C() = 12 < 14, is a (3, 14)MP.
Step . .(5.1), and , then there is no (3, 14)MP in , and go to (5.4).(5.4) i = 3.(5.1), and , then there is no (3, 14)MP in , and go to (5.4).(5.4).(5.1), and .(5.2) M() = 1.(5.3)Use the enumeration algorithm to search for (3, 14)MPs in , and (1, 2, 0, 1, 2, 1) is a (3, 14)MP.
Step . Let X = , and go to Step .
Step . , and halt.
Consequently, (2, 2, 0, 0, 1, 1), (1, 1, 0, 0, 2, 2), and (1, 2, 0, 1, 2, 1) are all (3, 14)MPs. Let , , and ; then, according to (1), = − − = + − +
That is, the probability of transmitting 3 units of commodity demand from to with the total transportation cost less than or equal to 14 is 0.64005 (note that the final result in [7, 9] is incorrect). The value, 0.64005, gives the manager information on the capability of the network to ensure the shipment of 3 units of commodity within the budget 14 from to . If this value is below the threshold set by the manager, more investments are needed to improve the logistics network. Otherwise, the performance of the logistics network is desirable.
5. Computational Experiments
As pointed out in Section 2, the algorithms of Yeh [8] and Niu and Xu [9] solve the MP problem without knowing all MPs in advance. Therefore, we conduct numerical experiments to check the performance of the proposed algorithm by comparing it with the algorithms of Yeh and Niu and Xu. A multistate network, as shown in Figure 2, is cited from the literature in [9]. The largest state vector and the cost vector are and , respectively. All of the three algorithms are implemented in a MATLAB program with a PC (AMD Athlon (tm) 64 X2 1.71 GHz C PU). To comprehensively embody the performance of the three algorithms, six different combinations of demand level and cost constraint are randomly chosen. The computational results are listed in Table 2.

It is observed from Table 2 that the proposed algorithm spends less time than Yeh’s [8] in searching for different MPs. In particular, the proposed algorithm has a better searching efficiency as the demand level grows. This is as expected, because Yeh’s algorithm is an exhaustive enumeration method in nature, whereas the proposed algorithm utilizes a decomposition technique to divide the search space, such that the search space of MPs can be dramatically reduced, and thereby MPs can be efficiently obtained. Furthermore, Table 2 shows that the proposed algorithm displays an advantage over the algorithm of Niu and Xu [9]. Both of them are partialenumeration methods, but the proposed technique is more effective in shortening the search space of MPs. Therefore, the proposed algorithm outperforms the two methods in solving MPs.
6. Conclusions
Reliability analysis is a fundamental tool for understanding the operational level of logistics networks. Logistics network can be modeled as a typical multistate flow network in which arcs are associated with multiple integer capacities, operational reliability, and unit transportation cost. The reliability index which takes into account both reliability and transportation cost can serve as a comprehensive performance assessment of logistics networks. This paper is devoted to proposing a new efficient algorithm for solving all MPs. A major contribution of this paper is that we develop a new approach to shorten the search space of MP. Specifically, the proposed algorithm employs an effective decomposition technique to divide the search space of MPs such that the search space of MPs can be dramatically reduced. A real example is used to illustrate the applicability of the proposed algorithm. Also, to check the performance of the proposed algorithm, computational experiments on a mediumsized logistics network are conducted, and the results show that the proposed algorithm has an advantage over the existing methods [5–9].
For future research, there is still a room for developing efficient methods for . For example, as with the existing algorithms [4–9], the proposed algorithm is also an enumeration method in nature, and thus it still needs to perform the enumeration of state vectors. Therefore, it is more desirable to develop new methods which are not grounded on the enumeration of state vectors. Furthermore, we only consider the transmission of single type of commodity in the network. Thus, it is more practical to extend the proposed algorithm to the reliability evaluation of multicommodity multistate networks under cost constraint in which multiple types of commodity are transmitting from the source to the destination.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work is partly supported by the National Natural Science Foundation of China (61300124 and 61403128), the Key Project of Science and Technology Research from the Education Department of Henan Province (13B120022 and 13B630034), and the Young Foundation of Henan Polytechnic University (Q201409).
References
 Y.K. Lin, “System reliability evaluation for a multistate supply chain network with failure nodes using minimal paths,” IEEE Transactions on Reliability, vol. 58, no. 1, pp. 34–40, 2009. View at: Publisher Site  Google Scholar
 Y.K. Lin, “Performance evaluation for the logistics system in case that capacity weight varies from arcs and types of commodity,” International Journal of Production Economics, vol. 107, no. 2, pp. 572–580, 2007. View at: Publisher Site  Google Scholar
 C.C. Jane, “Performance evaluation of logistics systems under cost and reliability considerations,” Transportation Research Part E: Logistics and Transportation Review, vol. 47, no. 2, pp. 130–137, 2011. View at: Publisher Site  Google Scholar
 C.C. Jane and Y.W. Laih, “Evaluating cost and reliability integrated performance of stochastic logistics systems,” Naval Research Logistics, vol. 59, no. 7, pp. 577–586, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 J.S. Lin, “Reliability evaluation of capacitatedflow networks with budget constraints,” IIE Transactions, vol. 30, no. 12, pp. 1175–1180, 1998. View at: Publisher Site  Google Scholar
 Y.K. Lin, “Reliability of a stochasticflow network with unreliable branches & nodes, under budget constraints,” IEEE Transactions on Reliability, vol. 53, no. 3, pp. 381–387, 2004. View at: Publisher Site  Google Scholar
 W.C. Yeh, “A new approach to evaluate reliability of multistate networks under the cost constraint,” Omega, vol. 33, no. 3, pp. 203–209, 2005. View at: Publisher Site  Google Scholar
 W. C. Yeh, “An improved method for multistate flow network reliability with unreliable nodes and a budget constraint based on path set,” IEEE Transactions on Systems, Man, and Cybernetics—Part A: Systems and Humans, vol. 41, no. 2, pp. 350–355, 2011. View at: Google Scholar
 Y.F. Niu and X.Z. Xu, “Reliability evaluation of multistate systems under cost consideration,” Applied Mathematical Modelling, vol. 36, no. 9, pp. 4261–4270, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 L. R. Ford and D. R. Fulkerson, Flows in Networks, Princeton University Press, Princeton, NJ, USA, 1962. View at: MathSciNet
 W.C. Yeh, “A novel method for the network reliability in terms of capacitatedminimumpaths without knowing minimumpaths in advance,” Journal of the Operational Research Society, vol. 56, no. 10, pp. 1235–1240, 2005. View at: Publisher Site  Google Scholar
 R. K. Ahuja, T. L. Magnanti, and J. B. Orlin, Network Flows: Theory, Algorithms, and Applications, Prentice Hall, Englewood Cliffs, NJ, USA, 1993. View at: MathSciNet
 C.C. Jane and Y.W. Laih, “A practical algorithm for computing multistate twoterminal reliability,” IEEE Transactions on Reliability, vol. 57, no. 2, pp. 295–302, 2008. View at: Publisher Site  Google Scholar
 R. K. Ahuja, M. Kodialam, A. K. Mishra, and J. B. Orlin, “Computational investigations of maximum flow algorithms,” European Journal of Operational Research, vol. 97, no. 3, pp. 509–542, 1997. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Qi Zeng and XiuZhen Xu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.