# Allocating Cost to Shippers in Pickup and Delivery Service

Authors: Ning Liu; Yaorong Cheng
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

This paper is going to study the shipper collaboration cost allocation problem in pickup and delivery (SCAPPD), which includes two crucial aspects, namely, the optimal redistribution of transport requests for minimization of the total cost of the alliance and equitable allocation of the shipper’s cost for keeping the stability of the alliance. We focus on the second issue and use core concept from cooperative game theory to develop an allocation methods, called Location Cost Allocation Model (LCAM), which take geographical location factor into account. As the goal of this study is to design an allocation tool that will be accepted in practice, we use the criteria of stability, consistency, and robustness to evaluate the LCAM, Shapley, Nucleolus, and the equal profit method (EPM). All four allocation methods are applied to both a case study and randomly generated instances. Our computational results show that the stability of the LCAM is as good as the nucleolus and EPM. In terms of consistency the LCAM performs best. LCAM also performs well in terms of robustness.


---

## Body

## 1. Introduction

In order to lower the total cost and ensure the quality of logistics, horizontal partners or even competitors usually apply the collaborative transportation strategies [1]. Products that come from various corporations are transported on the same route by a logistics company for the purpose of shipping cost reduction. Since the article published by Dantzig and Ramser [2], literature on the vehicle routing problem (VRP) has enjoyed a significant growth. Most of the literature attempts to find the overall optimal solution for VRP [3]. However, current VRP models aim to incorporate real-life complexities and different kinds of models have different application scenarios and conditions. This means that different VRP models have different cost allocation rules that cannot be generalized. Pickup and delivery problem (PDP) can be recognized as important vehicle routing problem where objects or people must be transported between origin and destination. In this article, we use the PDP definition proposed by Dumas et al. [4]. Most of the literature attempts to find the overall optimal solution for PDP (e.g., [5, 6]). Although significant discounts are offered to the group as a whole through collaboration transportation, each entity engaged cares about its own interest, which leads to the problem how the combined transportation costs on the same route should be distributed among partners.Mainstream approach that may be used to solve cost allocation problems for collaborative logistics is cooperative game approach. This problem has drawn widespread attention in recent years. Game theory provides several famous concepts of allocation, such as the core generated by Gilles [7], the Shapley value suggested by Shapley [8], and the nucleolus studied by Schmeidler [9]. The core theory requires individual, group, and coalition rationality, and these properties are the basis of cost allocation in the vast majority of the literature. But two drawbacks can be seen for core is not unique and may be empty. When it comes to certain calculation, Shapley values and nucleolus were compared with those suggested by the authors. Krajewska et al. [10] distributed the profit margins from the collaboration among partners by applying the Shapley value to cooperative game theory. There are many literatures using the Shapley value method for cost allocation or profit-sharing. (e.g., [11–14]). In order to discuss the cost allocation of gas and oil companies, Engevall et al. [15] developed the concept of demand nucleolus in travel salesman games. Yin et al. [16] studied the nucleolus cost allocation method for the vehicle routing problem with split delivery. Some new methods are discussed by scholars based on classic solution concept. Frisk et al. [17] conducted a research on a cooperative forest transportation planning problem as well as some typical cost allocation methods (including Shapley values and nucleolus), which put forward an equal profit method (EPM). Dai and Chen [18] paid close attention to the problem of carrier collaboration in pickup and delivery services and suggested three profit distribution mechanisms based on the Shapley value, the concept of proportional allocation and the contribution of each bearer in providing and serving requests.In terms of logistics cost analysis, the transportation cost allocation problem is basic for it is solely based on the route location. However, it causes some concerns about the method of distribution. A route cost allocation approach as well as a number of fairness standards was put forward by Fishburn and Pollak [19], where only the role of a single stop delivery cost was considered while other important factors were neglected, such as the relative position of clients. An approximate core was recognized and a geometric cost allocation approach with the idea of a moat was proposed by Faigle et al. [20]. According to the studies of Faigle et al. [20], Sun et al. [21] developed five fairness standards and contribution constrained packing model, which took the multiple fairness standards into consideration in cost allocation.This paper has the following structure. In the next section we describe the problem. The cost allocation model will be discussed in Section3, while, in Section 4, we conduct the comparison between our approach and others’ in terms of performance. And in Section 5, we will come to conclusion.


## 2. Problem Description

In terms of the shipper cooperation in the pickup and delivery services, a centralized planning framework is given to achieve the optimization of redistribution of transportation demands among the shippers to achieve the least amount of all-in costs after their alliance. In this case, it leads to a lower cost for the alliance as a whole comparing with the decentralized approach [22]. Previously, each shipper was only responsible for the fulfillment of his own requirement, which consists of the pickup node, delivery node, time window, and quantity. Before the alliance was set up, each shipper formulated his separate transportation planning to serve his own requirement with his own available transportation resources (vehicles). The calculation of the costs of shipper before the mutual partnership involved the individual cost of each shipper, respectively. Nonetheless, with the partnership relationship, the integration of the requests of all the shippers is made by a carrier, accompanied by the overall centralized transportation planning that is jointly worked out to reduce the total expenses to the smallest amount for all shippers. The total expenses of all the shippers after cooperating with one another can thus be derived from it. According to the collaborative transportation plan of the parties, the total costs shall be shared among the shippers with the aid of the cost allocation mechanism.As was illustrated earlier, two major problems can be summarized to be solved in the SCAPPD. One problem is the collaborative planning. In order to achieve the lowest overall cost for the shippers, vehicle tours are supposed to meet the requirements within the capacity limit of each vehicle and the time window limit of each request and a series of demands and requirements (pickup/delivery) should be guaranteed. Each vehicle tour consisting of multiple arcs in a transportation network is required to set off and get back to the depot node of the corresponding carrier. To date, plenty of research has been carried out on the algorithm to solve PDPTW, for instance, Ropke and Pisinger in 2006 and Bent and Van in 2006. Another problem lies in the equitable allocation of the total cost amount to each shipper concerned, leading to a lower cost for each shipper as a result of the cooperation. The paper emphatically probes into the second problem.To further explain the aforesaid SCAPPD, an instance analysis is conducted below. Given 8 nodes and one depot in the transportation network, the vehicle can hold 60 units. Assume that the point-to-point distance meets the Euclidean distance (see the transportation network as shown in Figures1 and 2). A total of 4 transportation requests have been received for the shippers, which are signified with r1 to r4, respectively. Time windows, quantities, and node of the requests can be seen in Table 1. For example, 20 units of freight at node 7 is extracted by request r1 and the same quantity is delivered to node 8.Table 1
Time windows, quantities, and node of the requests.
Request number 1 2 3 4 Pickup node 7 1 5 3 Delivery node 8 6 2 4 Pickup time window [0,30] [90,100] [160,180] [260,290] Delivery time window [40,80] [110,150] [170,250] [280,340] Quantity 20 40 30 50 Pickup node coordinate (42,58) (67,85) (80,65) (85,35) Delivery node coordinate (55,85) (58,70) (95,35) (64,46)Figure 1
Precollaboration transportation network.Figure 2
Postcollaboration transportation network.


## 3. Allocation Models and Evaluation Criteria

### 3.1. Basic Definitions and Properties

It is said that each solution and concept for the total cost allocation can meet some fairness standards. However, none of them is able to meet all the standards according to the literature till now. Part of the properties that is most frequently used is enumerated as follows. A coalitionS is equal to a subset of participants, while the grand coalition N refers to all players. We suppose that the opportunity to form and work together with others is provided for each participant in the coalitions. When the collaboration is established in the coalition S, the total (normal) cost CS shall be produced. In cooperative game theory, the cost function is seen as the Characteristic Cost Function and each participant is identified as a participant. It is possible to make a comparison between the cost allocation issue and the cooperative game.A cost allocation method is applied to the separation of total cost, i.e.,CN. Among all the participants, j∈N is said to be efficient, which can be expressed as ∑j∈Nyj=C(N), where yj is the share of cost allocated to participant j. The cost allocation can be recognized separate and reasonable, where no participant pays more than the independent cost, which is equal to the previous cost alone before the establishment of the alliance. Property can be presented as yj≤C(j). The core of the game is identified as those who meet the requirements of the cost allocation.(1)yj≤CS∑j∈Nyj=CN


### 3.2. The Shapley Value

The Shapley value [8] concept gives us a unique solution to deal with the cost allocation problem. As can be seen from the computation formula below, it is assumed that only one participant j is accessible to the grand coalition each time. The participant in the coalition is responsible for the marginal cost which means his accession will increase the total cost of the entire coalition. The amount received by participants in this program depends on the order in which participants joined. The cost assigned to participant j is equal(2)yj=∑S⊂N,j∈SS-1!N-S!N!cS-cS-j· shows the number of participants in the coalition considered. The sum of this formula is equal to that over all coalitions S that include the participant j. The value of cS-c(S-{j}) refers to the increase of the cost of the alliance with the increased the participant j, which is expressed by the marginal cost of participant j regarding the coalition S-{j}.


### 3.3. The Nucleolus

The cost allocation is identified in the calculation of nucleolus of a game to reduce the worst inequity to the lowest degree, so that the individual rationality can be fulfilled. That is to say, aiming at reducing the maximum dissatisfaction of any coalition as much as possible, we investigated how dissatisfied it is with the proposed allocationy among the coalition S. We delivered the dissatisfaction of cost allocation y about the coalition S, as a measurement of the amount by which coalition S falls short of its potential cS in the allocation y. The nucleolus refers to the cost allocation y with the largest correlated excess vector in dictionary order. The concept was further defined by Schmeidler [9]; i.e., the nucleolus indicates the unique vector which has the nonincreasing sorted vector of excess lexicographically minimized. The vector of excesses of an allocation y∈Rn is the vector θy∈R2n-1 defined as(3)θy=eS,ywitheS,y=cS-∑i∈Syi,∀S⊂N,S≠∅


### 3.4. Equal Profit Method

In 2010, for the purpose of minimization of maximum difference between pairs of related savings, Frisk et. al 2010 came up with a method inspired by finding a stable allocation, which was identifies as equal profit method (EPM). The relative saving of participants is expressed as follows:ci-yi/c({i})=1-yi/c({i}). Assuming the cost allocation is stable, we have that c({i})≥yi. Thus, yi/c({i})-yj/c({j}) refers to the difference in relative savings between two participants of i andj. LP problem is supposed to be addressed so as to find the allocation.(4)minfs.t.f≥yici-yjcj,∀i,j,∑j∈Syj≤cS,S⊂N,∑j∈Syj=cN,-yj≥0,∀iThe first constraint set acts as a measurement of the pairwise difference between the revenues of participants. Variablef in the goal minimize the biggest differences. All the stable allocations are defined with the other two constraint sets.


### 3.5. Location Cost Allocation Model (LCAM)

The formulation of LCAM is in fact a two-stage model. The first stage, we build the following linear programming model to achieve the minimum excess of every coalition under the condition of efficiency. The model (M1) is as follows:(5)minε(6)s.t.∑i∈Nπi=CN(7)∑i,n+i∈Sπi≤CS+ε,S⊂N(8)ε≥0,πi≥0Objective function (5) requires minimum ε. Constraint (6) corresponds to the efficiency condition, which states total route cost should be completely allocated. Constraint (7) corresponds to the coalition rationality conditions. By introducing ε, we add minimal penalties to the constraints that define the core to maintain the grand coalition “stability” and to obtain suboptimal solutions. Constraint (8) demonstrates the value scope of ε and πi. By solving this model, we can get the value of ε, denoted as ε∗.The second stage, as the geometric characteristics of the route, is the main basis for cost allocation; the objective of the model (M2) is to achieve a stable allocation set, considering the marginal cost and the dispersion degree of the nodes. The model (M2) is as follows:(9)minf(10)s.t.Ωi=πi+πn+i,i∈P+(11)Ωi′=c0i+ci,n+i+cn+i,0,i∈P+(12)∑i∈Nπi=CN(13)∑i,n+i∈Sπi≤CS+ε∗,S⊂N(14)Ωi≤Ωi′,i∈P+(15)Ωi-CN×ri+rn+i2≤f,i∈P+(16)ri=ri1+ri2,i∈N(17)ri1=mi∑i∈Nmi,i∈N(18)ri2=∑j∈Ncij∑i∈N∑j∈Ncij,i∈N,j∈N(19)mi=CS-CS-i,i∈S,S∈N(20)Ωk≤Ωk+1,∀k=1,2,…,n-1(21)f≥0,πi≥0Objective function (9) requires minimum f. Constraint (10) shows that the cost to the shipper allocated by alliance equals the cost combination of relative pickup and delivery nodes. Constraint (11) shows shipper i single stop delivery cost from depot. Constraint (12) corresponds to the efficiency condition, which states total route cost should be completely allocated. Constraint (13) shows that, when ε∗=0, the core exists; we select the solution meeting constraints (10) through (21) in the core. When ε∗>0, the core is empty or the solution meeting constraints (10) through (21) is not in the core. Constraint (14) corresponds to the individual condition that each shipper’s cost should not exceed its single stop delivery cost. Constraint (15) shows that the cost allocation should be positively correlated with the cost contribution factors and the correlation should be as high as possible. Constraint (16) represents the cost contribution factors, including ri1 and ri2. In Constraint (17), ri1 demonstrates the marginal cost contribution rate of node i. In Constraint (18), ri2 represents the dispersion degree of node i with other nodes. In Constraint (19), the marginal cost mi is the reduction of the total route length after the node i is removed, while preserving the original access sequence. Constraint (20) reflects the consistent relation between cost allocation and cost contribution; namely, the larger the cost contribution, the higher the cost allocation to the shipper. Constraint (21) demonstrates the value scope of f,πi. Relative parameter and variable are shown in Table 2.Table 2
Parameter and Variable.
Parameter C ( N ) Total cost is calculated by PDP algorithm (e.g. [5, 6]). Ω i Cost of shipperi allocated by alliance. Ω i ′ Shipperi single stop delivery cost from depot. c i j Cost between nodesi and j, without loss of generality, we set cij=dij. dij is Euclidean distance between nodes i and j. r i Cost contribution factors, includingri1 and ri2. r i 1 Marginal cost contribution rate of nodei. r i 2 Dispersion degree of nodei with other nodes. m i Marginal cost, which is the decreased value of the total route length after removing nodei while keeping others with the original visiting sequence. ( k ) Let(k) be the index of the kth shipper in the sequence of ordered cost contribution factors. S Set of node delivered by same vehicle. N The set of nodes. P + The set of pickup nodes.  Variable  π i Cost of nodei allocated by alliance. ε A minimum penalized slack in the constraints defining the core to keep the grand coalition “stable”. f Minimum difference between the cost of shipper and its cost based on geographical location.In this allocation model, the two factors ofri1 and ri2 are considered. The reason for considering both factors is because they represent different aspects. ri1 is calculated from mi, and mi is marginal cost of the node, which considers decreased value of the total route length after removing the node while keeping others with the original visiting sequence. ri2 is calculated from ∑j∈Ncij and ∑j∈Ncij can be understood as a geographical isolation value, which is a distance between node i and every other node. Just considering only one of these factors has drawbacks. We will use two examples to illustrate that considering two factors at the same time will make the distribution result more fair and reasonable.In the first example, we analyze the situation mentioned in Figure3. The red line represents request 1 of Shipper A (r1, from node 1 to node 4) and the blue line represents request 2 of Shipper B (r2, from node 2 to node 3). Request 2 is just right on the straight path of request 1, and its marginal cost is 0. If we only consider the factor ri1 (marginal cost contribution rate of node i), the result of cost allocation for request 2 is 0. Normally, this kind of allocation would, in practice, be unacceptable to the affected shippers. In this example, Shipper A will not be willing to bear all the costs.Figure 3
The situation in the first example.In the second example, we analyze the situation mentioned in Figure4. The relevant data comes from Table 1. The only difference is that the position of the nodes has changed.Figure 4
The situation in the second example.The calculation results show that the dispersion degree of request 2 (from node 1 to node 6) is approximately equal to 0.256, which is equal to the dispersion degree of the request 4 (from node 3 to node 4). If we only consider the factorri2 (dispersion degree of node i with other nodes), request 2 and request 4 should be allocated the same cost. But the marginal cost contribution rate of request 2 is 0.41, while the marginal cost contribution rate of request 4 is 0.15, indicating that the impact of request 4 to total cost is much smaller, and it should be allocated with smaller costs.Model LCAM is a linear programming model, which is implemented in LINGO 11 using linear solver to solve the linear programming problems. There are many classical algorithms for linear programming problems, including simplex algorithm [23], ellipsoid method [24], and Karmarkar’s algorithm [25].


### 3.6. Evaluation Criteria

To evaluate the allocation methods we use several evaluation criteria. As the goal of this study is to design an allocation tool that will be accepted in practice, we have selected evaluation criteria to assess which tool is most likely to be accepted in practice. We focus mainly on important criteria from the perspective of the shippers. First of all, the method is preferred to generatestable allocations. Stable allocations are in the core of the game, a well-known concept from cooperative game theory. For such allocations, no subset of shippers has an incentive to withdraw from the collaborative planning on the basis of the allocation. Secondly, the allocation should beconsistent with the geographical location of the shippers. That is, if geographical location varies, the allocation should change accordingly in a consistent way. We evaluate the consistency by using an ordinary least squares regression. Thirdly, an allocation method is preferablyrobust. If a shipper propose a similar transport request periodically, they will resent being allocated significantly different allocation at each shipment. To quantify robustness we use the coefficient of variation of the allocation.


## 3.1. Basic Definitions and Properties

It is said that each solution and concept for the total cost allocation can meet some fairness standards. However, none of them is able to meet all the standards according to the literature till now. Part of the properties that is most frequently used is enumerated as follows. A coalitionS is equal to a subset of participants, while the grand coalition N refers to all players. We suppose that the opportunity to form and work together with others is provided for each participant in the coalitions. When the collaboration is established in the coalition S, the total (normal) cost CS shall be produced. In cooperative game theory, the cost function is seen as the Characteristic Cost Function and each participant is identified as a participant. It is possible to make a comparison between the cost allocation issue and the cooperative game.A cost allocation method is applied to the separation of total cost, i.e.,CN. Among all the participants, j∈N is said to be efficient, which can be expressed as ∑j∈Nyj=C(N), where yj is the share of cost allocated to participant j. The cost allocation can be recognized separate and reasonable, where no participant pays more than the independent cost, which is equal to the previous cost alone before the establishment of the alliance. Property can be presented as yj≤C(j). The core of the game is identified as those who meet the requirements of the cost allocation.(1)yj≤CS∑j∈Nyj=CN


## 3.2. The Shapley Value

The Shapley value [8] concept gives us a unique solution to deal with the cost allocation problem. As can be seen from the computation formula below, it is assumed that only one participant j is accessible to the grand coalition each time. The participant in the coalition is responsible for the marginal cost which means his accession will increase the total cost of the entire coalition. The amount received by participants in this program depends on the order in which participants joined. The cost assigned to participant j is equal(2)yj=∑S⊂N,j∈SS-1!N-S!N!cS-cS-j· shows the number of participants in the coalition considered. The sum of this formula is equal to that over all coalitions S that include the participant j. The value of cS-c(S-{j}) refers to the increase of the cost of the alliance with the increased the participant j, which is expressed by the marginal cost of participant j regarding the coalition S-{j}.


## 3.3. The Nucleolus

The cost allocation is identified in the calculation of nucleolus of a game to reduce the worst inequity to the lowest degree, so that the individual rationality can be fulfilled. That is to say, aiming at reducing the maximum dissatisfaction of any coalition as much as possible, we investigated how dissatisfied it is with the proposed allocationy among the coalition S. We delivered the dissatisfaction of cost allocation y about the coalition S, as a measurement of the amount by which coalition S falls short of its potential cS in the allocation y. The nucleolus refers to the cost allocation y with the largest correlated excess vector in dictionary order. The concept was further defined by Schmeidler [9]; i.e., the nucleolus indicates the unique vector which has the nonincreasing sorted vector of excess lexicographically minimized. The vector of excesses of an allocation y∈Rn is the vector θy∈R2n-1 defined as(3)θy=eS,ywitheS,y=cS-∑i∈Syi,∀S⊂N,S≠∅


## 3.4. Equal Profit Method

In 2010, for the purpose of minimization of maximum difference between pairs of related savings, Frisk et. al 2010 came up with a method inspired by finding a stable allocation, which was identifies as equal profit method (EPM). The relative saving of participants is expressed as follows:ci-yi/c({i})=1-yi/c({i}). Assuming the cost allocation is stable, we have that c({i})≥yi. Thus, yi/c({i})-yj/c({j}) refers to the difference in relative savings between two participants of i andj. LP problem is supposed to be addressed so as to find the allocation.(4)minfs.t.f≥yici-yjcj,∀i,j,∑j∈Syj≤cS,S⊂N,∑j∈Syj=cN,-yj≥0,∀iThe first constraint set acts as a measurement of the pairwise difference between the revenues of participants. Variablef in the goal minimize the biggest differences. All the stable allocations are defined with the other two constraint sets.


## 3.5. Location Cost Allocation Model (LCAM)

The formulation of LCAM is in fact a two-stage model. The first stage, we build the following linear programming model to achieve the minimum excess of every coalition under the condition of efficiency. The model (M1) is as follows:(5)minε(6)s.t.∑i∈Nπi=CN(7)∑i,n+i∈Sπi≤CS+ε,S⊂N(8)ε≥0,πi≥0Objective function (5) requires minimum ε. Constraint (6) corresponds to the efficiency condition, which states total route cost should be completely allocated. Constraint (7) corresponds to the coalition rationality conditions. By introducing ε, we add minimal penalties to the constraints that define the core to maintain the grand coalition “stability” and to obtain suboptimal solutions. Constraint (8) demonstrates the value scope of ε and πi. By solving this model, we can get the value of ε, denoted as ε∗.The second stage, as the geometric characteristics of the route, is the main basis for cost allocation; the objective of the model (M2) is to achieve a stable allocation set, considering the marginal cost and the dispersion degree of the nodes. The model (M2) is as follows:(9)minf(10)s.t.Ωi=πi+πn+i,i∈P+(11)Ωi′=c0i+ci,n+i+cn+i,0,i∈P+(12)∑i∈Nπi=CN(13)∑i,n+i∈Sπi≤CS+ε∗,S⊂N(14)Ωi≤Ωi′,i∈P+(15)Ωi-CN×ri+rn+i2≤f,i∈P+(16)ri=ri1+ri2,i∈N(17)ri1=mi∑i∈Nmi,i∈N(18)ri2=∑j∈Ncij∑i∈N∑j∈Ncij,i∈N,j∈N(19)mi=CS-CS-i,i∈S,S∈N(20)Ωk≤Ωk+1,∀k=1,2,…,n-1(21)f≥0,πi≥0Objective function (9) requires minimum f. Constraint (10) shows that the cost to the shipper allocated by alliance equals the cost combination of relative pickup and delivery nodes. Constraint (11) shows shipper i single stop delivery cost from depot. Constraint (12) corresponds to the efficiency condition, which states total route cost should be completely allocated. Constraint (13) shows that, when ε∗=0, the core exists; we select the solution meeting constraints (10) through (21) in the core. When ε∗>0, the core is empty or the solution meeting constraints (10) through (21) is not in the core. Constraint (14) corresponds to the individual condition that each shipper’s cost should not exceed its single stop delivery cost. Constraint (15) shows that the cost allocation should be positively correlated with the cost contribution factors and the correlation should be as high as possible. Constraint (16) represents the cost contribution factors, including ri1 and ri2. In Constraint (17), ri1 demonstrates the marginal cost contribution rate of node i. In Constraint (18), ri2 represents the dispersion degree of node i with other nodes. In Constraint (19), the marginal cost mi is the reduction of the total route length after the node i is removed, while preserving the original access sequence. Constraint (20) reflects the consistent relation between cost allocation and cost contribution; namely, the larger the cost contribution, the higher the cost allocation to the shipper. Constraint (21) demonstrates the value scope of f,πi. Relative parameter and variable are shown in Table 2.Table 2
Parameter and Variable.
Parameter C ( N ) Total cost is calculated by PDP algorithm (e.g. [5, 6]). Ω i Cost of shipperi allocated by alliance. Ω i ′ Shipperi single stop delivery cost from depot. c i j Cost between nodesi and j, without loss of generality, we set cij=dij. dij is Euclidean distance between nodes i and j. r i Cost contribution factors, includingri1 and ri2. r i 1 Marginal cost contribution rate of nodei. r i 2 Dispersion degree of nodei with other nodes. m i Marginal cost, which is the decreased value of the total route length after removing nodei while keeping others with the original visiting sequence. ( k ) Let(k) be the index of the kth shipper in the sequence of ordered cost contribution factors. S Set of node delivered by same vehicle. N The set of nodes. P + The set of pickup nodes.  Variable  π i Cost of nodei allocated by alliance. ε A minimum penalized slack in the constraints defining the core to keep the grand coalition “stable”. f Minimum difference between the cost of shipper and its cost based on geographical location.In this allocation model, the two factors ofri1 and ri2 are considered. The reason for considering both factors is because they represent different aspects. ri1 is calculated from mi, and mi is marginal cost of the node, which considers decreased value of the total route length after removing the node while keeping others with the original visiting sequence. ri2 is calculated from ∑j∈Ncij and ∑j∈Ncij can be understood as a geographical isolation value, which is a distance between node i and every other node. Just considering only one of these factors has drawbacks. We will use two examples to illustrate that considering two factors at the same time will make the distribution result more fair and reasonable.In the first example, we analyze the situation mentioned in Figure3. The red line represents request 1 of Shipper A (r1, from node 1 to node 4) and the blue line represents request 2 of Shipper B (r2, from node 2 to node 3). Request 2 is just right on the straight path of request 1, and its marginal cost is 0. If we only consider the factor ri1 (marginal cost contribution rate of node i), the result of cost allocation for request 2 is 0. Normally, this kind of allocation would, in practice, be unacceptable to the affected shippers. In this example, Shipper A will not be willing to bear all the costs.Figure 3
The situation in the first example.In the second example, we analyze the situation mentioned in Figure4. The relevant data comes from Table 1. The only difference is that the position of the nodes has changed.Figure 4
The situation in the second example.The calculation results show that the dispersion degree of request 2 (from node 1 to node 6) is approximately equal to 0.256, which is equal to the dispersion degree of the request 4 (from node 3 to node 4). If we only consider the factorri2 (dispersion degree of node i with other nodes), request 2 and request 4 should be allocated the same cost. But the marginal cost contribution rate of request 2 is 0.41, while the marginal cost contribution rate of request 4 is 0.15, indicating that the impact of request 4 to total cost is much smaller, and it should be allocated with smaller costs.Model LCAM is a linear programming model, which is implemented in LINGO 11 using linear solver to solve the linear programming problems. There are many classical algorithms for linear programming problems, including simplex algorithm [23], ellipsoid method [24], and Karmarkar’s algorithm [25].


## 3.6. Evaluation Criteria

To evaluate the allocation methods we use several evaluation criteria. As the goal of this study is to design an allocation tool that will be accepted in practice, we have selected evaluation criteria to assess which tool is most likely to be accepted in practice. We focus mainly on important criteria from the perspective of the shippers. First of all, the method is preferred to generatestable allocations. Stable allocations are in the core of the game, a well-known concept from cooperative game theory. For such allocations, no subset of shippers has an incentive to withdraw from the collaborative planning on the basis of the allocation. Secondly, the allocation should beconsistent with the geographical location of the shippers. That is, if geographical location varies, the allocation should change accordingly in a consistent way. We evaluate the consistency by using an ordinary least squares regression. Thirdly, an allocation method is preferablyrobust. If a shipper propose a similar transport request periodically, they will resent being allocated significantly different allocation at each shipment. To quantify robustness we use the coefficient of variation of the allocation.


## 4. Numerical Results

### 4.1. Case Study

We study the methods applied to a case study to gain insight into the performance in practice. Going back to Figures1 and 2 and Table 1, the specific description of the case has been described in Section 2. This instance allows us to evaluate the feasibility and effectiveness of LCAM method.In the first stage of LCAM, we can getε∗=0; the core consists of all the stable profit imputations; that is, Cv={x∈Iv∣∑i∈Sxi≥vS,∀S∈2N∖{∅}}. However, cost is negative profit. So we use the cost saving replacing cost obtained from the cooperation among the shippers to further study the nature of the game. The cost saving CSS of the coalition S is as follows: CSS=∑p∈STCp-TCS. As Table 3 indicates, each possible coalition S the minimum total cost TCS and related cost savings CSS can be calculated.Table 3
CoalitionS the minimum total cost TCS and corresponding cost savings CSS.
S { A } { B } { C } { D } { A B } { A C } { A D } { B C } T C S 76 89 133 95 95 161 145 159 C S S 0 0 0 0 70 48 26 63  S { B D } { C D } { A B C } { A B D } { A C D } { B C D } { A B C D } T C S 146 134 170 157 162 160 171 C S S 38 94 128 103 142 157 222Provided the cost savings for the individual shippers in the separate operation is zero, the cost savings generally increase with the expansion of the scale of the coalition. Ideally, largest cost savings would correspond with the grand coalition. For this instance, theCS function is monotonic, super additive, and balanced. As can be revealed in Figure 5, the core of this four-shipper game can be figured out by TUGlab [26]. The corresponding core vertices are shown in Table 4.Table 4
The corresponding core vertices.
Vertex y a y b y c y d 1 0 70 58 94 2 0 70 119 33 3 0 80 48 94 4 0 77 119 26 5 0 80 116 26 6 48 80 0 94 7 65 63 0 94 8 26 77 119 0 9 26 80 116 0 10 65 38 119 0 11 48 80 94 0 12 65 63 94 0Figure 5
The core of this four-shipper game.What is noteworthy is, as can be observed from this example, that the core-vertex solutions and core solutions are unfair when allocating small savings to one or more shippers. For example, Shipper A is allocated with nothing among five of the core vertices. In fact, this kind of distribution is unacceptable to related shippers as a rule. It justifies why it is far from enough to simply require the solution to be stable. It is more convenient to impose the fairness with the geographical location criteria as in the case of LCAM dose.As Table5 indicates, in this instance, the cost allocation shall meet the constraints of Shipper C > Shipper D > Shipper A > Shipper B according to Constraint (20). However, as shown in Table 5, the postcollaboration cost of each shipper given by Shapley or EPM is Shipper C > Shipper D > Shipper B > Shipper A, and they all do not satisfy condition (20). Nucleolus satisfy the condition (20) luckily. Because Shapley, nucleolus, and EPM do not consider the impact of geographical location of nodes on cost allocation by definition. So in terms of geographical location criterion the LCAM would be preferred over the other methods discussed in paper. It should be noted that, in some cases, there may be coincidences to meeting constraints (20) via Shapley value, nucleolus, or EPM.Table 5
Cost contribution factors.
Shippers r i 1 r i 2 r i A 17.95% 24.96% 42.91% B 18.76% 22.75% 41.51% C 39.34% 27.52% 66.86% D 23.95% 24.76% 48.71%


### 4.2. Randomly Generated Instances

In this section, we use three criteria presented in Section3.6 to evaluate the Shapley, nucleolus, EPM, and LCAM, through two randomly generated examples. All results are obtained by using an Intel® Core™ i5-5200 CPU @ 2.2GHz with 8 GB of RAM.Stability.First, the core turned out to be nonempty for each instance. Because the nucleolus, EPM, and LCAM guarantee a solution in the core whenever it is nonempty, their allocations are stable by definition. But the Shapley value is not constructed by implicitly using any of the core criteria. Therefore, from a stability point of view, the Shapley value performs worst. However, for each instance in Tables 6 and 8, all four mechanisms develop a cost allocation that stays at the core.Table 6
Results of the random instances from 1 to 20.
Instance 1 2 3 4 5 6 7 8 9 10 Pre-collaboration costs 76.3 143.9 355.2 282.3 224.8 168.9 86.4 593.5 632.5 638.6 88.6 147.7 125.1 311.8 271.1 359.0 141.0 743.1 808.1 614.8 133.3 147.3 118.4 266.3 203.5 252.9 136.0 114.0 129.5 83.0 95.5 173.9 118.8 256.4 262.8 315.6 119.1 394.3 422.3 401.1 Post-collaboration costs 171.0 306.0 439.3 524.7 347.8 502.7 176.0 1050.8 1097.3 1038.3 Cost collaboration of Shapley 32.9 62.9 298.1 155.6 81.6 50.3 26.0 366.9 363.7 450.1 36.2 57.0 50.4 179.2 108.4 183.1 55.3 446.3 516.8 390.1 61.9 68.1 41.9 92.5 65.6 96.9 55.5 54.8 53.0 43.7 40.0 118.0 48.9 97.4 92.2 172.4 39.2 182.8 163.8 154.4 Cost collaboration of Nucleolus 35.9 58.2 301.9 169.0 92.0 68.4 40.0 273.1 505.8 580.0 33.9 47.9 46.5 187.5 92.4 163.2 47.8 370.4 403.9 289.1 75.4 67.0 38.9 73.4 77.8 98.0 52.7 63.5 72.1 53.7 25.8 132.9 52.0 94.8 85.6 173.1 35.5 343.8 115.5 115.6 Cost collaboration of EPM 33.1 71.8 217.5 132.6 81.3 77.4 31.5 338.1 348.3 381.6 38.5 73.8 76.6 146.5 98.0 164.6 51.4 423.2 445.1 367.4 57.9 73.6 72.5 125.2 73.5 116.0 49.6 64.9 71.3 49.6 41.5 86.8 72.7 120.4 95.0 144.7 43.5 224.6 232.6 239.7 Cost collaboration of LCAM 36.7 64.3 261.4 157.1 81.7 60.3 23.8 342.3 338.5 415.4 35.5 56.6 63.7 178.2 113.1 176.0 57.2 394.8 459.2 373.8 57.2 67.3 54.9 92.2 64.3 95.5 58.1 102.7 125.7 62.2 41.6 117.8 59.3 97.2 88.7 170.9 36.9 211.0 173.9 186.9  Instance 11 12 13 14 15 16 17 18 19 20  Pre-collaboration costs 533.2 241.5 675.2 202.7 109.1 134.6 288.8 619.1 291.8 467.7 547.6 227.0 483.5 26.4 202.5 140.1 134.6 575.2 62.5 289.9 578.6 417.6 279.2 211.1 394.0 189.6 301.3 758.5 234.9 421.6 268.4 200.1 285.5 389.2 303.4 159.0 356.3 832.0 258.7 201.5 Post-collaboration costs 820.3 549.7 838.5 444.6 608.0 244.9 422.1 1002.1 338.0 758.2 Cost collaboration of Shapley 252.0 107.5 371.7 80.1 52.1 47.6 106.9 177.2 126.1 333.7 205.1 78.3 201.1 10.0 95.2 48.2 49.8 166.6 20.7 113.3 261.5 297.8 107.4 102.8 269.2 95.3 106.7 247.8 96.7 229.7 101.7 66.2 158.3 251.7 191.5 53.8 158.7 410.5 94.5 81.5 Cost collaboration of Nucleolus 279.8 105.1 392.3 178.3 73.9 55.2 163.8 192.5 207.9 367.7 142.9 69.9 110.8 14.3 82.4 47.5 75.3 193.5 34.2 94.5 226.8 306.6 135.7 79.6 265.3 96.1 67.0 200.3 62.1 182.1 170.8 68.1 199.7 172.4 186.4 46.1 116.0 415.8 33.8 113.9 Cost collaboration of EPM 226.9 122.2 328.5 108.6 65.8 52.9 112.7 222.8 116.3 256.8 233 114.9 235.3 14.1 122.0 55.0 52.6 207.0 24.9 159.2 246.2 211.3 135.8 113.2 237.4 74.5 117.7 273.0 93.6 231.5 114.2 101.3 138.9 208.7 182.8 62.5 139.1 299.3 103.2 110.7 Cost collaboration of LCAM 258.4 109.9 348.0 89.1 82.1 49.0 105.5 145.3 116.6 310.9 192.8 79.7 197.1 17.3 109.4 47.8 69.7 145.3 52.4 125.1 254.3 283.1 113.2 109.1 240.2 97.9 97.1 225.5 87.9 219.9 114.8 77 180.3 229.1 176.3 50.2 149.8 486.0 81.1 102.3Consistency.In order to assess the consistency of the different allocation methods, we need to design random instances. We evaluate consistency based on marginal distance and average distance. Marginal distance is the decreased value of the total route length after removing node i while keeping others with the original visiting sequence. Average distance is the average distance between node i and every other node. We evaluate the consistency by using an ordinary least squares regression.So far, without benchmark instance for SCAPPD, 20 randomly generated instances are applied to assess the performance of this mechanism. Each instance contains a vector. Set the number of nodes to 9 (1 node as depot node and 8 others as shipper node). The coordinates of each node are freely developed by the computer. The capacity of each vehicle is 60 units. Each request is related to 2 nodes, pickup node as well as delivery node. These requests are developed through random selection of one pickup nodei (i≠depot) and a delivery node (j≠iandj≠depot), each related to at most one request. Each request is related to a randomly introduced quantity of freight u units no bigger than a predefined number that is defined to 60. The time window for all request service is set to [0,480] (480 minutes = 8 hours).According to the experimental results, the cost reduction of each shipper resulting from the four mechanisms can be different, which is due to the fact that the core of each instance includes multiple possible allocations and that different definitions of fairness in each mechanism can result in different core allocations. The results of the LCAM pay more attention to the impact of geographical location on the allocation. We evaluate it based on the marginal distance and the average distance to other shippers. As these two factors simultaneously influence the allocate result, we have performed an ordinary least squares regression to evaluate each allocation method, where the allocate result is the dependent variable. The explanatory variables are a constant, the marginal distance, and the average distance to other shippers.A summary of the regression results is shown in Table7. For each allocation method we have reported the estimated coefficient and its significance level. We use R2 as a measure to compare the allocation methods. R2 is the proportion of variation in the dependent variable that is explained by the explanatory variables. We see that the LCAM has the highest R2 and hence performs best on the geographical location criterion, while Shapley and EPM perform relatively poor. Nucleolus performs worst.Table 7
Regression results to measure the geographical location consistency of the allocation methods.
Explanatory variable Shapley Nucleolus EPM LCAM Coeff. Sig. Coeff. Sig. Coeff. Sig. Coeff. Sig. Constant -15.3 0.153 -10.758 0.433 -2.509 0.819 -15.815 0.057 Marg. dist 0.107 0.000 0.104 0.000 0.065 0.000 0.086 0.000 Avg. dist. to other cust. 0.446 0.000 0.433 0.000 0.514 0.000 0.511 0.000  R 2 0.868 0.79 0.807 0.907Table 8
Results of the random instances similar to the instance 1.
Instance 1 2 3 4 5 6 7 8 9 10 Pre-collaboration costs 76.3 72.2 79.0 73.5 73.5 82.4 75.4 79.1 72.8 86.4 88.6 87.0 81.0 91.1 83.4 81.7 91.0 87.0 82.7 79.6 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 95.5 99.1 104.4 93.6 100.5 83.2 104.1 104.1 102.2 93.0 Post-collaboration costs 171.0 157.2 163.4 166.5 175.2 168.6 173.1 175.3 177.0 181.6 Cost collaboration of Shapley 32.8 25.0 31.4 26.2 33.9 37.3 35.1 35.9 32.6 44.0 36.2 32.2 29.6 38.2 34.2 29.5 35.4 35.0 34.0 30.3 62.0 60.5 58.7 63.8 62.8 68.3 57.5 60.4 62.5 64.4 40.0 39.5 43.6 38.3 44.4 33.5 45.1 44.0 47.9 43.0 Cost collaboration of Nucleolus 35.8 29.1 32.2 28.3 37.4 39.7 41.9 38.5 35.0 46.1 34.0 32.0 28.3 36.8 32.1 27.9 33.8 32.8 30.8 28.6 57.3 55.0 55.5 59.7 59.2 63.8 50.2 56.8 58.0 59.6 44.0 41.0 47.5 41.8 46.6 37.1 47.1 47.2 53.2 47.2 Cost collaboration of EPM 33.1 29.0 32.5 31.3 32.9 36.5 32.3 34.3 33.0 40.0 38.5 34.9 33.3 38.8 37.4 36.2 39.0 37.8 37.4 36.9 57.9 53.5 54.8 56.7 59.8 59.0 57.1 57.9 60.3 61.7 41.5 39.8 42.9 39.8 45.1 36.8 44.6 45.2 46.3 43.1 Cost collaboration of LCAM 36.7 28.0 34.0 30.0 38.6 39.8 39.6 39.8 36.6 45.7 35.5 31.9 30.1 37.8 34.0 30.3 35.2 35.0 34.6 30.7 57.2 57.0 54.5 59.3 57.1 62.8 52.4 56.1 57.6 59.6 41.6 40.2 44.8 39.4 45.5 35.8 46.0 44.4 48.2 45.5  Instance 11 12 13 14 15 16 17 18 19 20  Pre-collaboration costs 73.6 85.9 75.4 74.8 79.8 81.9 79.8 65.8 83.1 69.0 88.8 76.2 83.2 91.7 97.7 93.3 94.8 74.5 102.7 83.5 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 109.4 100.0 103.9 92.3 104.4 93.6 106.8 89.5 103.7 97.7 Post-collaboration costs 167.8 173.9 166.4 170.7 174.7 174.0 173.3 169.1 190.1 182.1 Cost collaboration of Shapley 29.6 39.7 31.4 31.3 30.7 34.9 33.8 31.3 40.7 36.0 35.1 27.6 31.6 37.6 42.0 37.8 38.3 30.4 45.4 36.4 55.6 61.5 58.9 64.5 58.1 63.1 56.7 66.8 60.7 64.7 47.4 45.2 44.5 37.3 44.0 38.2 44.6 40.7 43.3 45.0 Cost collaboration of Nucleolus 34.3 39.8 33.8 35.2 32.5 37.1 38.1 34.9 45.9 42.6 34.4 26.5 29.7 35.5 40.0 35.2 36.9 28.0 42.6 33.7 49.3 57.7 54.7 59.8 54.3 59.2 51.4 61.2 56.2 57.6 49.8 49.8 48.2 40.3 48.0 42.4 47.0 45.1 45.3 48.2 Cost collaboration of EPM 30.5 37.8 31.7 32.6 33.6 35.4 33.3 30.7 37.4 32.8 36.8 33.5 35.0 40.0 41.1 40.4 39.6 34.7 46.2 39.6 55.2 58.6 56.0 58.0 56.1 57.7 55.7 62.0 59.9 63.3 45.3 44.0 43.7 40.2 43.9 40.5 44.6 41.6 46.6 46.4 Cost collaboration of LCAM 33.4 41.8 34.7 35.4 34.2 38.2 38.1 35.6 45.3 41.2 34.0 29.0 32.1 37.1 40.9 37.4 37.2 31.1 44.0 35.8 52.4 56.5 54.4 59.5 54.5 58.3 53.0 60.2 56.9 59.5 48.1 46.6 45.3 38.7 45.1 40.1 45.0 42.2 43.9 45.6Robustness.In order to assess the robustness of the different allocation methods, we need to design another random example. We are interested in the variation of the target shipper’s allocation in a similar task. So we randomly generated 20 instances similar to the instance in Table 1, in which there is no change in request 3, and observe the variation of the allocation of Shipper C to different methods.We measure the robustness as the standard deviation relative to the average allocated emission, which is the coefficient of variation (CoV). To see which allocation method performs best in terms of robustness, we compare the CoV for the different methods. The coefficient of variation for the instances is shown in Table9. We see that EPM and LCAM are the best method in which CoV does not exceed 5%. LCAM performs worse than EPM in terms of robustness because LCAM considers more constraints; it is closer to the actual situation and seeks solutions in a smaller feasible domain.Table 9
Coefficient of variation (CoV) for the random instances similar to the instance in Table1.
Target shipper Allocation method Std. dev. of allocation Average allocation Cov Shipper C Shapley 3.4 61.6 0.0546 Nucleolus 3.7 56.8 0.0643 EPM 2.6 58.1 0.0443 LCAM 2.8 56.9 0.0492


## 4.1. Case Study

We study the methods applied to a case study to gain insight into the performance in practice. Going back to Figures1 and 2 and Table 1, the specific description of the case has been described in Section 2. This instance allows us to evaluate the feasibility and effectiveness of LCAM method.In the first stage of LCAM, we can getε∗=0; the core consists of all the stable profit imputations; that is, Cv={x∈Iv∣∑i∈Sxi≥vS,∀S∈2N∖{∅}}. However, cost is negative profit. So we use the cost saving replacing cost obtained from the cooperation among the shippers to further study the nature of the game. The cost saving CSS of the coalition S is as follows: CSS=∑p∈STCp-TCS. As Table 3 indicates, each possible coalition S the minimum total cost TCS and related cost savings CSS can be calculated.Table 3
CoalitionS the minimum total cost TCS and corresponding cost savings CSS.
S { A } { B } { C } { D } { A B } { A C } { A D } { B C } T C S 76 89 133 95 95 161 145 159 C S S 0 0 0 0 70 48 26 63  S { B D } { C D } { A B C } { A B D } { A C D } { B C D } { A B C D } T C S 146 134 170 157 162 160 171 C S S 38 94 128 103 142 157 222Provided the cost savings for the individual shippers in the separate operation is zero, the cost savings generally increase with the expansion of the scale of the coalition. Ideally, largest cost savings would correspond with the grand coalition. For this instance, theCS function is monotonic, super additive, and balanced. As can be revealed in Figure 5, the core of this four-shipper game can be figured out by TUGlab [26]. The corresponding core vertices are shown in Table 4.Table 4
The corresponding core vertices.
Vertex y a y b y c y d 1 0 70 58 94 2 0 70 119 33 3 0 80 48 94 4 0 77 119 26 5 0 80 116 26 6 48 80 0 94 7 65 63 0 94 8 26 77 119 0 9 26 80 116 0 10 65 38 119 0 11 48 80 94 0 12 65 63 94 0Figure 5
The core of this four-shipper game.What is noteworthy is, as can be observed from this example, that the core-vertex solutions and core solutions are unfair when allocating small savings to one or more shippers. For example, Shipper A is allocated with nothing among five of the core vertices. In fact, this kind of distribution is unacceptable to related shippers as a rule. It justifies why it is far from enough to simply require the solution to be stable. It is more convenient to impose the fairness with the geographical location criteria as in the case of LCAM dose.As Table5 indicates, in this instance, the cost allocation shall meet the constraints of Shipper C > Shipper D > Shipper A > Shipper B according to Constraint (20). However, as shown in Table 5, the postcollaboration cost of each shipper given by Shapley or EPM is Shipper C > Shipper D > Shipper B > Shipper A, and they all do not satisfy condition (20). Nucleolus satisfy the condition (20) luckily. Because Shapley, nucleolus, and EPM do not consider the impact of geographical location of nodes on cost allocation by definition. So in terms of geographical location criterion the LCAM would be preferred over the other methods discussed in paper. It should be noted that, in some cases, there may be coincidences to meeting constraints (20) via Shapley value, nucleolus, or EPM.Table 5
Cost contribution factors.
Shippers r i 1 r i 2 r i A 17.95% 24.96% 42.91% B 18.76% 22.75% 41.51% C 39.34% 27.52% 66.86% D 23.95% 24.76% 48.71%


## 4.2. Randomly Generated Instances

In this section, we use three criteria presented in Section3.6 to evaluate the Shapley, nucleolus, EPM, and LCAM, through two randomly generated examples. All results are obtained by using an Intel® Core™ i5-5200 CPU @ 2.2GHz with 8 GB of RAM.Stability.First, the core turned out to be nonempty for each instance. Because the nucleolus, EPM, and LCAM guarantee a solution in the core whenever it is nonempty, their allocations are stable by definition. But the Shapley value is not constructed by implicitly using any of the core criteria. Therefore, from a stability point of view, the Shapley value performs worst. However, for each instance in Tables 6 and 8, all four mechanisms develop a cost allocation that stays at the core.Table 6
Results of the random instances from 1 to 20.
Instance 1 2 3 4 5 6 7 8 9 10 Pre-collaboration costs 76.3 143.9 355.2 282.3 224.8 168.9 86.4 593.5 632.5 638.6 88.6 147.7 125.1 311.8 271.1 359.0 141.0 743.1 808.1 614.8 133.3 147.3 118.4 266.3 203.5 252.9 136.0 114.0 129.5 83.0 95.5 173.9 118.8 256.4 262.8 315.6 119.1 394.3 422.3 401.1 Post-collaboration costs 171.0 306.0 439.3 524.7 347.8 502.7 176.0 1050.8 1097.3 1038.3 Cost collaboration of Shapley 32.9 62.9 298.1 155.6 81.6 50.3 26.0 366.9 363.7 450.1 36.2 57.0 50.4 179.2 108.4 183.1 55.3 446.3 516.8 390.1 61.9 68.1 41.9 92.5 65.6 96.9 55.5 54.8 53.0 43.7 40.0 118.0 48.9 97.4 92.2 172.4 39.2 182.8 163.8 154.4 Cost collaboration of Nucleolus 35.9 58.2 301.9 169.0 92.0 68.4 40.0 273.1 505.8 580.0 33.9 47.9 46.5 187.5 92.4 163.2 47.8 370.4 403.9 289.1 75.4 67.0 38.9 73.4 77.8 98.0 52.7 63.5 72.1 53.7 25.8 132.9 52.0 94.8 85.6 173.1 35.5 343.8 115.5 115.6 Cost collaboration of EPM 33.1 71.8 217.5 132.6 81.3 77.4 31.5 338.1 348.3 381.6 38.5 73.8 76.6 146.5 98.0 164.6 51.4 423.2 445.1 367.4 57.9 73.6 72.5 125.2 73.5 116.0 49.6 64.9 71.3 49.6 41.5 86.8 72.7 120.4 95.0 144.7 43.5 224.6 232.6 239.7 Cost collaboration of LCAM 36.7 64.3 261.4 157.1 81.7 60.3 23.8 342.3 338.5 415.4 35.5 56.6 63.7 178.2 113.1 176.0 57.2 394.8 459.2 373.8 57.2 67.3 54.9 92.2 64.3 95.5 58.1 102.7 125.7 62.2 41.6 117.8 59.3 97.2 88.7 170.9 36.9 211.0 173.9 186.9  Instance 11 12 13 14 15 16 17 18 19 20  Pre-collaboration costs 533.2 241.5 675.2 202.7 109.1 134.6 288.8 619.1 291.8 467.7 547.6 227.0 483.5 26.4 202.5 140.1 134.6 575.2 62.5 289.9 578.6 417.6 279.2 211.1 394.0 189.6 301.3 758.5 234.9 421.6 268.4 200.1 285.5 389.2 303.4 159.0 356.3 832.0 258.7 201.5 Post-collaboration costs 820.3 549.7 838.5 444.6 608.0 244.9 422.1 1002.1 338.0 758.2 Cost collaboration of Shapley 252.0 107.5 371.7 80.1 52.1 47.6 106.9 177.2 126.1 333.7 205.1 78.3 201.1 10.0 95.2 48.2 49.8 166.6 20.7 113.3 261.5 297.8 107.4 102.8 269.2 95.3 106.7 247.8 96.7 229.7 101.7 66.2 158.3 251.7 191.5 53.8 158.7 410.5 94.5 81.5 Cost collaboration of Nucleolus 279.8 105.1 392.3 178.3 73.9 55.2 163.8 192.5 207.9 367.7 142.9 69.9 110.8 14.3 82.4 47.5 75.3 193.5 34.2 94.5 226.8 306.6 135.7 79.6 265.3 96.1 67.0 200.3 62.1 182.1 170.8 68.1 199.7 172.4 186.4 46.1 116.0 415.8 33.8 113.9 Cost collaboration of EPM 226.9 122.2 328.5 108.6 65.8 52.9 112.7 222.8 116.3 256.8 233 114.9 235.3 14.1 122.0 55.0 52.6 207.0 24.9 159.2 246.2 211.3 135.8 113.2 237.4 74.5 117.7 273.0 93.6 231.5 114.2 101.3 138.9 208.7 182.8 62.5 139.1 299.3 103.2 110.7 Cost collaboration of LCAM 258.4 109.9 348.0 89.1 82.1 49.0 105.5 145.3 116.6 310.9 192.8 79.7 197.1 17.3 109.4 47.8 69.7 145.3 52.4 125.1 254.3 283.1 113.2 109.1 240.2 97.9 97.1 225.5 87.9 219.9 114.8 77 180.3 229.1 176.3 50.2 149.8 486.0 81.1 102.3Consistency.In order to assess the consistency of the different allocation methods, we need to design random instances. We evaluate consistency based on marginal distance and average distance. Marginal distance is the decreased value of the total route length after removing node i while keeping others with the original visiting sequence. Average distance is the average distance between node i and every other node. We evaluate the consistency by using an ordinary least squares regression.So far, without benchmark instance for SCAPPD, 20 randomly generated instances are applied to assess the performance of this mechanism. Each instance contains a vector. Set the number of nodes to 9 (1 node as depot node and 8 others as shipper node). The coordinates of each node are freely developed by the computer. The capacity of each vehicle is 60 units. Each request is related to 2 nodes, pickup node as well as delivery node. These requests are developed through random selection of one pickup nodei (i≠depot) and a delivery node (j≠iandj≠depot), each related to at most one request. Each request is related to a randomly introduced quantity of freight u units no bigger than a predefined number that is defined to 60. The time window for all request service is set to [0,480] (480 minutes = 8 hours).According to the experimental results, the cost reduction of each shipper resulting from the four mechanisms can be different, which is due to the fact that the core of each instance includes multiple possible allocations and that different definitions of fairness in each mechanism can result in different core allocations. The results of the LCAM pay more attention to the impact of geographical location on the allocation. We evaluate it based on the marginal distance and the average distance to other shippers. As these two factors simultaneously influence the allocate result, we have performed an ordinary least squares regression to evaluate each allocation method, where the allocate result is the dependent variable. The explanatory variables are a constant, the marginal distance, and the average distance to other shippers.A summary of the regression results is shown in Table7. For each allocation method we have reported the estimated coefficient and its significance level. We use R2 as a measure to compare the allocation methods. R2 is the proportion of variation in the dependent variable that is explained by the explanatory variables. We see that the LCAM has the highest R2 and hence performs best on the geographical location criterion, while Shapley and EPM perform relatively poor. Nucleolus performs worst.Table 7
Regression results to measure the geographical location consistency of the allocation methods.
Explanatory variable Shapley Nucleolus EPM LCAM Coeff. Sig. Coeff. Sig. Coeff. Sig. Coeff. Sig. Constant -15.3 0.153 -10.758 0.433 -2.509 0.819 -15.815 0.057 Marg. dist 0.107 0.000 0.104 0.000 0.065 0.000 0.086 0.000 Avg. dist. to other cust. 0.446 0.000 0.433 0.000 0.514 0.000 0.511 0.000  R 2 0.868 0.79 0.807 0.907Table 8
Results of the random instances similar to the instance 1.
Instance 1 2 3 4 5 6 7 8 9 10 Pre-collaboration costs 76.3 72.2 79.0 73.5 73.5 82.4 75.4 79.1 72.8 86.4 88.6 87.0 81.0 91.1 83.4 81.7 91.0 87.0 82.7 79.6 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 95.5 99.1 104.4 93.6 100.5 83.2 104.1 104.1 102.2 93.0 Post-collaboration costs 171.0 157.2 163.4 166.5 175.2 168.6 173.1 175.3 177.0 181.6 Cost collaboration of Shapley 32.8 25.0 31.4 26.2 33.9 37.3 35.1 35.9 32.6 44.0 36.2 32.2 29.6 38.2 34.2 29.5 35.4 35.0 34.0 30.3 62.0 60.5 58.7 63.8 62.8 68.3 57.5 60.4 62.5 64.4 40.0 39.5 43.6 38.3 44.4 33.5 45.1 44.0 47.9 43.0 Cost collaboration of Nucleolus 35.8 29.1 32.2 28.3 37.4 39.7 41.9 38.5 35.0 46.1 34.0 32.0 28.3 36.8 32.1 27.9 33.8 32.8 30.8 28.6 57.3 55.0 55.5 59.7 59.2 63.8 50.2 56.8 58.0 59.6 44.0 41.0 47.5 41.8 46.6 37.1 47.1 47.2 53.2 47.2 Cost collaboration of EPM 33.1 29.0 32.5 31.3 32.9 36.5 32.3 34.3 33.0 40.0 38.5 34.9 33.3 38.8 37.4 36.2 39.0 37.8 37.4 36.9 57.9 53.5 54.8 56.7 59.8 59.0 57.1 57.9 60.3 61.7 41.5 39.8 42.9 39.8 45.1 36.8 44.6 45.2 46.3 43.1 Cost collaboration of LCAM 36.7 28.0 34.0 30.0 38.6 39.8 39.6 39.8 36.6 45.7 35.5 31.9 30.1 37.8 34.0 30.3 35.2 35.0 34.6 30.7 57.2 57.0 54.5 59.3 57.1 62.8 52.4 56.1 57.6 59.6 41.6 40.2 44.8 39.4 45.5 35.8 46.0 44.4 48.2 45.5  Instance 11 12 13 14 15 16 17 18 19 20  Pre-collaboration costs 73.6 85.9 75.4 74.8 79.8 81.9 79.8 65.8 83.1 69.0 88.8 76.2 83.2 91.7 97.7 93.3 94.8 74.5 102.7 83.5 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 133.3 109.4 100.0 103.9 92.3 104.4 93.6 106.8 89.5 103.7 97.7 Post-collaboration costs 167.8 173.9 166.4 170.7 174.7 174.0 173.3 169.1 190.1 182.1 Cost collaboration of Shapley 29.6 39.7 31.4 31.3 30.7 34.9 33.8 31.3 40.7 36.0 35.1 27.6 31.6 37.6 42.0 37.8 38.3 30.4 45.4 36.4 55.6 61.5 58.9 64.5 58.1 63.1 56.7 66.8 60.7 64.7 47.4 45.2 44.5 37.3 44.0 38.2 44.6 40.7 43.3 45.0 Cost collaboration of Nucleolus 34.3 39.8 33.8 35.2 32.5 37.1 38.1 34.9 45.9 42.6 34.4 26.5 29.7 35.5 40.0 35.2 36.9 28.0 42.6 33.7 49.3 57.7 54.7 59.8 54.3 59.2 51.4 61.2 56.2 57.6 49.8 49.8 48.2 40.3 48.0 42.4 47.0 45.1 45.3 48.2 Cost collaboration of EPM 30.5 37.8 31.7 32.6 33.6 35.4 33.3 30.7 37.4 32.8 36.8 33.5 35.0 40.0 41.1 40.4 39.6 34.7 46.2 39.6 55.2 58.6 56.0 58.0 56.1 57.7 55.7 62.0 59.9 63.3 45.3 44.0 43.7 40.2 43.9 40.5 44.6 41.6 46.6 46.4 Cost collaboration of LCAM 33.4 41.8 34.7 35.4 34.2 38.2 38.1 35.6 45.3 41.2 34.0 29.0 32.1 37.1 40.9 37.4 37.2 31.1 44.0 35.8 52.4 56.5 54.4 59.5 54.5 58.3 53.0 60.2 56.9 59.5 48.1 46.6 45.3 38.7 45.1 40.1 45.0 42.2 43.9 45.6Robustness.In order to assess the robustness of the different allocation methods, we need to design another random example. We are interested in the variation of the target shipper’s allocation in a similar task. So we randomly generated 20 instances similar to the instance in Table 1, in which there is no change in request 3, and observe the variation of the allocation of Shipper C to different methods.We measure the robustness as the standard deviation relative to the average allocated emission, which is the coefficient of variation (CoV). To see which allocation method performs best in terms of robustness, we compare the CoV for the different methods. The coefficient of variation for the instances is shown in Table9. We see that EPM and LCAM are the best method in which CoV does not exceed 5%. LCAM performs worse than EPM in terms of robustness because LCAM considers more constraints; it is closer to the actual situation and seeks solutions in a smaller feasible domain.Table 9
Coefficient of variation (CoV) for the random instances similar to the instance in Table1.
Target shipper Allocation method Std. dev. of allocation Average allocation Cov Shipper C Shapley 3.4 61.6 0.0546 Nucleolus 3.7 56.8 0.0643 EPM 2.6 58.1 0.0443 LCAM 2.8 56.9 0.0492


## 5. Conclusion

In this paper we study SCAPPD, in which we need to allocate the shipper’s cost fairly. We analyze three other allocation methods based on well-known solution concepts in cooperative game theory: the Shapley value, nucleolus, and EPM. Furthermore, LCAM are proposed which consider marginal cost and the dispersion degree of nodes in cost allocation on the basis of the core allocation method. Each allocation method is evaluated in terms of stability, consistency, and robustness. In our numerical experiments, the stability of the LCAM is as good as the nucleolus and EPM. In terms of consistency the LCAM performs best. LCAM also performs well in terms of robustness. LCAM offers more choices for the shippers to choose the cost allocation mechanism they want.In this study, each shipper is assumed to have only one pickup and delivery request in model LCAM. The model considering one-to-one (1-1) problems is the basis of many-to-many (M-M) problems. Hence, it is interesting to extend this research to consider a more general case where each shipper has multiple pickup and delivery requests.



---

*Source: 1283539-2019-01-29.xml*