[論文レビュー] 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

ひとこと要約

本稿では、流量需要をパrameterとするネットワークフロー・モデルにおけるすべての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 区分線形で、場合によっては不連続なエッジコストを有するネットワークにおいて、流量需要をパrameterとするすべてのWardrop均衡をどのように体系的に計算できるか？
• RQ2 ホモトピーに基づくアプローチは、ネットワーク均衡におけるフロー超過と頂点ポテンシャルの構造をどのように保持できるか？
• RQ3 複数のエッジで非微分性点が同時に発生しない解曲線において、どのような計算複雑度の保証が達成可能か？
• RQ4 アルゴリズムは、多商品および退化解曲線を扱うためにどのように拡張できるか？

主な発見

• アルゴリズムは、選択された超過ベクトルの経路に沿って、無向単商品ネットワークにおけるすべてのWardrop均衡を正しく計算する。
• 不連続コスト関数をジャンプを伴う区分線形関数としてモデル化することにより、有向ネットワークおよびエッジ容量に対しても拡張可能である。
• 非退化インスタンス（同時に非微分性を持つエッジが1つまで）では、アルゴリズムは出力多項式時間で実行される。
• 退化ケースでは、双対写像の線形セグメントを凸計画法で計算することで、アルゴリズムは依然として出力多項式時間のままである。
• 凸計画法のアプローチにより、アルゴリズムは多商品ネットワークへの拡張が可能になり、計算効率が保たれる。