[论文解读] Computing all wardrop equilibria parametrized by the flow demand

Max Klimm, Philipp Warode|arXiv (Cornell University)|Jan 6, 2019

Advanced Differential Equations and Dynamical Systems被引用 4

一句话总结

该论文提出了一种基于同伦的算法，通过顶点势与流盈余向量之间的分段线性双射，将网络流模型中所有Wardrop均衡作为流量需求的函数进行计算。该算法可处理连续与不连续的成本函数，包括有向边和容量限制，并在非退化情况下实现输出多项式时间复杂度；通过凸规划方法，可进一步推广至多类别流量和退化情形。

ABSTRACT

We develop an algorithm that computes for a given undirected or directed network with flow-dependent piece-wise linear edge cost functions all Wardrop equilibria as a function of the flow demand. Our algorithm is based on Katzenelson's homotopy method for electrical networks. The algorithm uses a bijection between vertex potentials and flow excess vectors that is piecewise linear in the potential space and where each linear segment can be interpreted as an augmenting flow in a residual network. The algorithm iteratively increases the excess of one or more vertex pairs until the bijection reaches a point of non-differentiability. Then, the next linear region is chosen in a simplex-like pivot step and the algorithm proceeds. We first show that this algorithm correctly computes all Wardrop equilibria in undirected single-commodity networks along the chosen path of excess vectors. We then adapt our algorithm to also work for discontinuous cost functions which allows to model directed edges and/or edge capacities. Our algorithm is output-polynomial in non-degenerate instances where the solution curve never hits a point where the cost function of more than one edge becomes non-differentiable. For degenerate instances we still obtain an output-polynomial algorithm computing the linear segments of the bijection by a convex program. The latter technique also allows to handle multiple commodities.

研究动机与目标

• 开发一种系统化方法，用于计算网络流模型中随流量需求变化的所有Wardrop均衡。
• 将现有均衡计算技术扩展至处理不连续边成本函数，例如建模有向边或容量限制的情形。
• 通过在非退化情况下实现输出多项式时间复杂度，确保计算效率。
• 通过凸规划方法将该方法推广至多类别网络，以处理退化情形。

提出的方法

• 该算法采用Katzenberger的同伦方法并将其适配于网络流问题，建立了顶点势与流盈余向量之间的分段线性双射。
• 该双射的每个线性段对应于残差网络中的一个增广流，从而实现对解路径的逐步遍历。
• 算法通过迭代增加顶点对的盈余，直至到达不可微点，该点标志着线性区域的改变。
• 采用类似单纯形的枢轴步骤选择下一个线性区域，确保对解曲线的系统性探索。
• 对于不连续成本函数，该方法通过具有跳跃间断点的分段线性成本函数来建模边的行为。
• 在退化情形下，使用凸规划计算双射的线性段，从而实现向多类别流的扩展。

实验结果

研究问题

• RQ1 如何在具有分段线性、可能不连续的边成本的网络中，系统性地计算所有Wardrop均衡作为流量需求的函数？
• RQ2 能否将基于同伦的方法适配以保持网络均衡中流盈余与顶点势的结构？
• RQ3 当解曲线避开多个边同时出现不可微点的区域时，可实现何种计算复杂度保证？
• RQ4 如何将该算法扩展以处理多类别流量及退化解曲线？

主要发现

• 该算法正确计算了在选定盈余向量路径上，无向单类别网络中的所有Wardrop均衡。
• 通过使用具有跳跃的分段线性函数建模不连续成本函数，该方法可扩展至有向网络和边容量限制。
• 在非退化情形下（即任意时刻仅有一个边的成本函数不可微），该算法运行在输出多项式时间内。
• 在退化情形下，通过使用凸规划计算双射的线性段，该算法仍保持输出多项式时间复杂度。
• 凸规划方法使得该算法可扩展至多类别网络，同时保持计算效率。