Research Article  Open Access
A PathBased Gradient Projection Algorithm for the CostBased System Optimum Problem in Networks with Continuously Distributed Value of Time
Abstract
The costbased system optimum problem in networks with continuously distributed value of time is formulated as a pathbased form, which cannot be solved by the FrankWolfe algorithm. In light of magnitude improvement in the availability of computer memory in recent years, pathbased algorithms have been regarded as a viable approach for traffic assignment problems with reasonably large network sizes. We develop a pathbased gradient projection algorithm for solving the costbased system optimum model, based on GoldsteinLevitinPolyak method which has been successfully applied to solve standard user equilibrium and system optimum problems. The Sioux Falls network tested is used to verify the effectiveness of the algorithm.
1. Introduction
The traffic assignment problem consists in determining which routes to assign to the drivers who travel on a transportation network from some origins and some destinations. Wardrop [1] stated two principles for determining the assignment. The first principle leads to an equilibrium state in which no user can reduce his travel time by using an alternative route. The second principle induces a system optimal state in which the total travel time of all users (or the average journey time in a network) is minimum regarding the benefit of the society (i.e., whole traffic). It is well known that, in the standard traffic network equilibrium model, a socalled marginalcost toll (equivalent to the negative externality that an additional individual imposes on other users of the system) can drive a user equilibrium flow pattern to a system optimum [2]. This conventional traffic network equilibrium model typically assumes that users’ VOTs are identical, that is, homogeneous users. However, user heterogeneity is manifested in the fact that some travelers take slower paths to avoid tolls while others choose tolled roads to save time. Some studies pointed out that the VOT varies significantly across individuals because of different socioeconomic characteristics, trip purposes, attitudes, and inherent preferences [3–6].
The concept of value of time (VOT) plays a central role in road pricing analysis as it describes how users make tradeoffs between money and time in response to road toll charges. The network disutility can be measured in travel time or travel cost. It is obvious that different system optimal (SO) flow patterns will be obtained if we use different units (time or money) to measure the system disutility. Previous studies that address user heterogeneity can be classified into two categories [7]. The first is the multiclass approach in which the entire feasible VOT range is divided into several predetermined intervals according to a discrete VOT distribution or some socioeconomic characteristics [2, 8–10]. The second lets VOTs be continuously distributed across the population of trips, which is regarded to be more rational than others [11–15]. In the these studies, including those dealing with tworoute or highway/transit twomode problems [12, 16–19], the continuously distributed VOT between each OD pair was treated as a random variable following a probability density function [13, 15, 20–22].
For more realistically capturing the travelers’ path choice behavior in response to toll charges, Wu and Huang [23] assumed that the VOT varies significantly across individuals because of their different socioeconomic characteristics, trip purposes, attitudes, and inherent preferences. Furthermore, the VOT of each user is assumed to be deterministic and constant, because the factors influencing VOT keep unchanged within a certain time period. They extended the work of Yang and Huang [2] to the case with continuously distributed value of time across users, for finding anonymous tolls to realize target flow pattern in networks with continuously distributed value of time. To find anonymous tolls to realize system optimum in cost units, Wu and Huang [23] proposed a costbased system optimum model in general networks with continuously distributed value of time but did not study its traffic assignment algorithms. The calculation of the model is very difficult [24]. It is formulated in a pathbased form and cannot be solved by the FrankWolfe algorithm which is linkbased and cannot provide path flows.
In the past, pathbased algorithms, even those not enumerating all possible paths, were traditionally discarded by transportation researchers for solving largescale network problems because of intensive memory requirements and the difficulties in manipulating and storing paths [25]. However, pathbased algorithms automatically provide not only the linkflow solution but also the pathflow solution that may be required in certain applications. In light of magnitude improvement in the availability of computer memory in recent years, recent research on pathbased algorithms has demonstrated and established that it is a viable approach for traffic assignment problems with reasonably large network sizes [26–29]. Chen et al. [25] stated that much of the attention has been focused on two particular algorithms: the disaggregate simplicial decomposition (DSD) algorithm and the gradient projection (GP) algorithm. Application of GP to solve the traffic assignment problem is relatively new, but GP is as good as or better than DSD in direct comparisons.
In this paper, we develop a GP algorithm for the proposed costbased system optimum model, based on the GoldsteinLevitinPolyak (GLP) GP method formulated by Bertsekas [30] for general nonlinear multicommodity problems which is successfully used for solving traffic assignment problems [26]. The proposed GLP algorithm includes five key operations: (1) assigning the flow on each path to the links along to find the total link flows, (2) computing link travel times and path travel times and sorting path travel times in decreasing order, (3) computing the first derivative lengths (marginal social path travel costs) for all paths between each OD pair and finding the shortest first derivative lengths for each OD pair, (4) finding the second derivative lengths for each OD pair, and (5) updating the path flows using the second derivative lengths as scaling.
In the next section, a system optimum problem in cost units is formulated in fixed demand networks with continuously distributed value of time. In Section 3, we develop a GP algorithm for the costbased system optimum model. A numerical example is presented in Section 4 for testing the GP algorithm. Section 5 concludes the paper.
2. CostBased System Optimum Problem
In this section, we formulate the system optimum in cost units in a fixed demand network with heterogeneous users in terms of a continuously distributed VOT. Let be a directed network where is the set of nodes and the set of links. Each link has an associated flowdependent travel time which is assumed to be differentiable, convex, and monotonically increasing subject to the link flow . Let be the set of all OD pairs, the set of all paths connecting OD pair , and the number of paths between each OD pair. Let be the set of drivers (users) on path and the number of these drivers; . We have , , where equals 1 if path between OD pair contains link , and 0 otherwise. All users on link are charged by an exogenously given toll . Each path between each OD pair is associated with a travel time and a travel cost (dollar) , , and .
Let be the travel demand between each OD pair . Instead of assuming a unified VOT for the whole population, in this paper we set a unique, specific VOT for each tripmaker. Let the distribution of VOTs across individuals between each OD pair be characterized by a continuous function . Furthermore, let the population between each OD pair be ordered in decreasing order of their VOTs; that is, , where is the th user between each OD pair , as demonstrated in Figure 1. Note that here the VOT distribution is given. In fact, there exist many studies which address estimation of VOTs. By observing choices among alternative combinations of cost and travel time, based on stated preference (SP) analysis, due to the ability to control time and cost variables [31], information about the relative weighting of cost and time can be inferred, and from this the distribution of VOT can be derived. In general, literature on estimation of VOTs refers to two main models: multinomial logit model [32, 33] and mixed logit model [34, 35]. Recently, the mixed logit approach is popular because it does not have to assume the irrelevance of independent alternatives (IIA) property.
Let denote the paths between each OD pair , satisfying . The system optimum naturally requires that the users with higher VOTs should choose faster paths and those with lower VOTs choose slower paths. Otherwise, the total travel cost can be reduced by switching a higher VOT user on a slower path into a faster path. For the sake of notational consistency, we define and . The system optimum in cost units can be formulated as the following minimization problem: subject to where , . Recall that the users between each OD pair are arranged in decreasing order of their VOTs. This determines the integral formulation appearing in (1).
The firstorder optimality conditions of the above minimization problem are as follows: where is the minimal travel cost between each OD pair .
Note that the sum of the second and third terms on the left hand side of (4) is the total externality caused by the th user of OD pair , , or the th user of path . The second term states that this new tripmaker imposes additional costs on all users who may belong to other OD pairs and other paths and have their own VOTs but traverse the links of path . This is the traditional congestion externality reported in the literature. The third term states that the costs of paths 1 to are reduced since the th user’s path choice changes the VOTs of the users on these paths (this user chooses path rather than other longer paths). This is another kind of externality attributed to the VOT distribution of OD trips. Let
Also, for the remainder of the paper when we refer to the first derivative lengths we mean the first derivatives of the objective function, which can also be regarded as marginal social travel cost corresponding to path flow on path .
Equations (4) and (5) can then be rewritten as
Therefore, in a network where each user has a unique VOT, (7) state that, at optimality, the marginal social travel costs on all the used paths connecting a given OD pair are equal and less than or equal to those on all the unused paths. This is similar to a standard system optimum, in which, at optimality, the marginal social travel times on all the used paths between each given OD pair are equal and less than or equal to those on all the unused paths.
Clearly, the proposed model is more complicated than a standard system optimum. According to the extreme value theorem that a continuous function in the closed and bounded space attains its maximum and minimum, the above minimization program must attain its minimum value. The most common algorithm used to solve traffic assignment problems is linkbased FrankWolfe algorithm introduced by LeBlanc et al. [36]. However, the objective function is, in general, not convex since the path travel costs depending not only on its own path but also on other paths are inseparable and asymmetric in terms of the path flows. Fortunately, for the differentiable and bounded objective, the convexity of the constraints allows the development of a gradient projection algorithm.
3. PathBased Traffic Assignment Algorithm
In this section, we adopt the GoldsteinLevitinPolyak algorithm to the traffic assignment problem. In each iteration, the travel demand constraints (2) are eliminated by reformulating the pathflow variables in terms of nonshortest path flows in terms of the first derivative lengths to make projection operation simpler. This is implemented by partitioning into the shortest path flow and the nonshortest path flows and writing constraints (2) in terms of as follows:
It should be noted that is the shortest path in terms of the first derivative lengths but is the shortest path in terms of path travel times. Putting constraints (8) into the objective function, we obtain a new formulation with just the nonnegativity constraints on the nonshortest path flows as the decision variables. Consider where is the reformulated objective function only including the nonshortest path flows for all OD pairs. This reformulation is a program with only nonnegativity constraints. The first and second derivatives of the new objective function can be easily derived as follows: where
Note that the second derivative is, in general, not positive. This is different from the transformed objective function of the standard system optimum model [25]. Thus, in each iteration, the scaled GLP algorithm updates the path flows according to the following iteration equations: where is the iteration number, is the step size, is the shortest path in terms of the first derivative lengths between each OD pair , is a diagonal scaling, equivalent to the second derivative regarding path , that is, , and are the first derivative lengths along path and path between each OD pair , and denotes the projection operation.
Note that can be explained by (11), and (15) certainly satisfies the flow conservation constraints in (2) and nonnegativity conditions in (3).
With the above flow update equations, the complete algorithmic steps can be summarized as follows.
Step 1 (initialization). Set , for all , and perform allornothing assignments. This yields path flows , for all , and link flows , for all . Set iteration counter . Initialize the pathset with the shortest path for each OD pair .
Step 2 (update). Set , for all . Sort path travel time in decreasing order; that is, , where , .
Step 3. Update the first derivative lengths of all the paths in , for all .
Step 4 (direction finding). Find the shortest paths in terms of marginal social travel cost between each OD pair , based on . If , then add it to and record . Otherwise, tag the shortest among the paths in as .
Step 5 (move). Set the new path flows as follows:
where , for all , and is a scalar stepsize modifier.
Also, , for all , .
Update link flows .
Step 6 (convergence test). Determine the total deviation of marginal social travel costs between all OD pairs . If , then stop. Otherwise, set and go to Step 2.
For convenience, it is better to keep constant (i.e., ) since is used for scaling [26]. Given any starting set of path flows, there exists such that if the sequence generated by this algorithm converges to the objective function (1) [37].
4. Numerical Example
In the section, a numerical example is presented to illustrate the effectiveness of the proposed model and algorithm. The test network is the Sioux Falls network, which is a medium sized network with 24 nodes and 76 links, as shown in Figure 3. This paper considers only one origindestination pair from node 1 to node 20 for conveniently providing a complete picture of the travel pattern which demonstrates all used path flows, all used link flows, and path choice behaviors. We here use the standard Bureau of Public Road (BPR) link cost function for our numerical study. The functional form is given by where is the freeflow cost, is the flow, and is the link capacity.
Assume that the total demand is 50, and all network users which are differed by a continuously distributed VOT are given in Figure 2 as follows:
Note that, in the numerical example, we show that achieves very good convergence rate.
The GP algorithm provides a complete picture of the travel pattern and keeps track of the distribution of the OD flows among the different routes, as shown in Table 1 and Figure 3. In the Sioux Falls network, there are many paths between OD pair 1–20. However, the number of the used paths to define the system optimum keeps small; that is, there are only thirteen used paths which are given in a decreasing order in terms of path travel times. It can be seen that, at optimality, people with high VOTs would choose faster paths, whereas people with low VOTs would choose slower paths. This is consistent with the requirements for the system optimum in Section 2.

Figure 3 shows the optimal linkflow distribution. Note that the real lines refer to the used links and the dotted lines refer to the unused paths. The link flows equal the values beside the corresponding links and are also graphically demonstrated by the widths of the corresponding lines. It can be easily seen that in the network most links in the forward directions are used, but most links in the backward directions are unused.
Table 2 compares the total travel times and total travel costs under different types of traffic assignment, respectively. At the costbased system optimum, the total travel cost is the lowest whereas the total travel time is the highest. At the timebased system optimum, the total travel time is the lowest, and the total travel cost is more than the costbased system optimum but less than the user equilibrium.

Figure 4 depicts the pattern of convergence toward the minimum. This convergence pattern is demonstrated in terms of the reduction in the value of the objective function from iteration to iteration. After the 21st iteration, the marginal contribution of each successive iteration becomes smaller and smaller as the algorithm proceeds (this property is the basis for the convergence criterion used in the example). After the 40th iteration, the value of the objective function almost keeps unchanged and is approximately equal to the minimum value 13281, verifying the effectiveness of the proposed GP algorithm.
5. Conclusions
The VOT varies significantly across individuals because of the different socioeconomic characteristics, trip purposes, attitudes, and inherent preferences. Furthermore, the VOT of each user is assumed to be deterministic and constant because the factors influencing VOT keep unchanged within a certain time period. We provide a theoretical investigation of the system optimum problem in fixed demand networks with continuous VOT distribution. This system optimum problem is formulated as a minimization program in cost units, which is more complicated than standard system optimum. This is because its objective function is, in general, nonconvex since path travel costs depending not only on its own path but also on other paths are inseparable and asymmetric in terms of the path flows. Considering the convexity of the constraints, we have developed a pathbased GLP algorithm for solving this model. The test results in the Sioux Falls network show the effectiveness of the algorithm.
This study lays a solid foundation for road pricing to realize the costbased system optimum in fixed demand general networks with heterogeneous users. For future research, this framework can be extended to the case of elastic demand and further investigate how to determine anonymous tolls to realize the system optimum in cost units.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by grants from the National Basic Research Program of China (2012CB725401), the National Natural Science Foundation of China (71071004 and 71271004), and the National High Technology Research and Development Program of China (863 Program) (2012AA112401).
References
 J. G. Wardrop, “Some theoretical aspects of road traffic research,” Proceedings of the Institute of Civil Engineers II, no. 1, pp. 325–378, 1952. View at: Google Scholar
 H. Yang and H.J. Huang, “The multiclass, multicriteria traffic network equilibrium and systems optimum problem,” Transportation Research B, vol. 38, no. 1, pp. 1–15, 2004. View at: Publisher Site  Google Scholar
 K. A. Small and J. Yan, “The value of “value pricing” of roads: secondbest pricing and product differentiation,” Journal of Urban Economics, vol. 49, no. 2, pp. 310–336, 2001. View at: Publisher Site  Google Scholar
 D. Brownstone and K. A. Small, “Valuing time and reliability: assessing the evidence from road pricing demonstrations,” Transportation Research A, vol. 39, no. 4, pp. 279–293, 2005. View at: Publisher Site  Google Scholar
 K. A. Small, C. Winston, and J. Yan, “Uncovering the distribution of motorists' preferences for travel time and reliability,” Econometrica, vol. 73, no. 4, pp. 1367–1382, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Cirillo and K. W. Axhausen, “Evidence on the distribution of values of travel time savings from a sixweek diary,” Transportation Research A, vol. 40, no. 5, pp. 444–457, 2006. View at: Publisher Site  Google Scholar
 C.C. Lu, H. S. Mahmassani, and X. Zhou, “A bicriterion dynamic user equilibrium traffic assignment model and solution algorithm for evaluating dynamic road pricing strategies,” Transportation Research C, vol. 16, no. 4, pp. 371–389, 2008. View at: Publisher Site  Google Scholar
 R. Lindsey, “Existence, uniqueness, and trip cost function properties of user equilibrium in the bottleneck model with multiple user classes,” Transportation Science, vol. 38, no. 3, pp. 293–314, 2004. View at: Publisher Site  Google Scholar
 D. Han and H. Yang, “The multiclass, multicriterion traffic equilibrium and the efficiency of congestion pricing,” Transportation Research E, vol. 44, no. 5, pp. 753–773, 2008. View at: Publisher Site  Google Scholar
 A. Clark, A. Sumalee, S. Shepherd, and R. Connors, “On the existence and uniqueness of first best tolls in networks with multiple user classes and elastic demand,” Transportmetrica, vol. 5, no. 2, pp. 141–157, 2009. View at: Publisher Site  Google Scholar
 T. C. Lam and K. A. Small, “The value of time and reliability: measurement from a value pricing experiment,” Transportation Research E, vol. 37, no. 23, pp. 231–251, 2001. View at: Publisher Site  Google Scholar
 E. T. Verhoef and K. A. Small, “Product differentiation on roads: constrained congestion pricing with heterogeneous users,” Journal of Transport Economics and Policy, vol. 38, no. 1, pp. 127–156, 2004. View at: Google Scholar
 P. Marcotte and D. L. Zhu, “Existence and computation of optimal tolls in multiclass network equilibrium problems,” Operations Research Letters, vol. 37, no. 3, pp. 211–214, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. van den Berg and E. T. Verhoef, “Congestion tolling in the bottleneck model with heterogeneous values of time,” Transportation Research B, vol. 45, no. 1, pp. 60–78, 2011. View at: Publisher Site  Google Scholar
 D. Zhu, C. Li, and G. Chen, “Existence of strongly valid tolls for multiclass network equilibrium problems,” Acta Mathematica Scientia B, vol. 32, no. 3, pp. 1093–1101, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Mayet and M. Hansen, “Congestion pricing with continuously distributed values of time,” Journal of Transport Economics and Policy, vol. 34, no. 3, pp. 359–369, 2000. View at: Google Scholar
 F. Xiao and H. Yang, “Efficiency loss of private road with continuously distributed valueoftime,” Transportmetrica, vol. 4, no. 1, pp. 19–32, 2008. View at: Publisher Site  Google Scholar
 F. Xiao and H. M. Zhang, “Paretoimproving and selfsustainable pricing for the morning commute with nonidentical commuters,” Transportation Science, pp. 1–11, 2013. View at: Publisher Site  Google Scholar
 L. J. Tian, H. J. Huang, and H. Yang, “Tradable credit schemes for managing bottleneck congestion and modal split with heterogeneous users,” Transportation Research E, vol. 54, pp. 1–13, 2013. View at: Google Scholar
 R. Cole, Y. Dodis, and T. Roughgarden, “Pricing network edges for heterogeneous selfish users,” in Proceedings of the 35 Annual ACM Symposium on Theory of Computing, pp. 521–530, ACM, San Diego, Calif, USA, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 L. Fleischer, K. Jain, and M. Mahdian, “Tolls for heterogeneous selfish users in multicommodity networks and generalized congestion games,” in Proceedings of the 45th Annual IEEE Symposium on Foundations of Computer Science (FOCS '04), pp. 277–285, Rome, Italy, October 2004. View at: Google Scholar
 G. Karakostas and S. G. Kolliopoulos, “Edge pricing of multicommodity networks for heterogeneous selfish users,” in Proceedings of the 45th Annual IEEE Symposium on Foundations of Computer Science (FOCS '04), pp. 268–276, Rome, Italy, October 2004. View at: Google Scholar
 W. X. Wu and H. J. Huang, “Finding anonymous tolls to realize target flow pattern in networks with continuously distributed value of time,” submitted to Transportation Research B, 2013. View at: Google Scholar
 R. B. Dial, “Networkoptimized road pricing: part II: algorithms and examples,” Operations Research, vol. 47, no. 2, pp. 327–336, 1999. View at: Publisher Site  Google Scholar
 A. Chen, D.H. Lee, and R. Jayakrishnan, “Computational study of stateoftheart pathbased traffic assignment algorithms,” Mathematics and Computers in Simulation, vol. 59, no. 6, pp. 509–518, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Jayakrishnan, W. K. Tsai, J. N. Prashker, and S. Rajadhyaksha, “Faster pathbased algorithm for traffic assignment,” Transportation Research Record, no. 1443, pp. 75–83, 1994. View at: Google Scholar
 T. Larsson and M. Patriksson, “Simplicial decomposition with disaggregated representation for the traffic assignment problem,” Transportation Science, vol. 26, no. 1, pp. 4–17, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Sun, R. Jayakrishnan, and W. K. Tsai, “Computational study of a pathbased algorithm and its variants for static traffic assignment,” Transportation Research Record, no. 1537, pp. 106–115, 1996. View at: Google Scholar
 M. Tatineni, H. Edwards, and D. Boyce, “Comparison of disaggregate simplicial decomposition and FrankWolfe algorithms for useroptimal route choice,” Transportation Research Record, no. 1617, pp. 157–162, 1998. View at: Google Scholar
 D. P. Bertsekas, “On the GoldsteinLevitinPolyak gradient projection method,” IEEE Transactions on Automatic Control, vol. 21, no. 2, pp. 174–184, 1976. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 H. Gunn, “An introduction to the valuation of traveltime savings and losse,” in Handbook of Transport Modelling, Elsevier Science, 2000. View at: Google Scholar
 D. Brownstone and K. Train, “Forecasting new product penetration with flexible substitution patterns,” Journal of Econometrics, vol. 89, no. 12, pp. 109–129, 1998. View at: Publisher Site  Google Scholar
 K. E. Train, “Recreation demand models with taste differences over people,” Land Economics, vol. 74, no. 2, pp. 230–239, 1998. View at: Publisher Site  Google Scholar
 S. Algers, P. Bergstrom, M. Dahlberg, and J. L. Dillen, “Mixed logit estimation of the value of travel time,” UppsalaWorking Paper Series 1998:15, 1998. View at: Google Scholar
 S. Hess, M. Bierlaire, and J. W. Polak, “Estimation of value of traveltime savings using mixed logit models,” Transportation Research A, vol. 39, no. 23, pp. 221–236, 2005. View at: Publisher Site  Google Scholar
 L. J. LeBlanc, E. K. Morlok, and W. P. Pierskalla, “An efficient approach to solving the road network equilibrium traffic assignment problem,” Transportation Research, vol. 9, no. 5, pp. 309–318, 1975. View at: Publisher Site  Google Scholar
 D. Bertsekas and R. Gallager, Data Networks, Prentice Hall, Englewood Cliffs, NJ, USA, 2nd edition, 1992.
Copyright
Copyright © 2014 WenXiang Wu and HaiJun Huang. 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.