[論文レビュー] General auction method for real-valued optimal transport

Joseph D. Walsh, Luca Dieci|arXiv (Cornell University)|May 18, 2017

Game Theory and Voting Systems参考文献 13被引用数 5

ひとこと要約

本稿では、整数値データに限らない一般化されたオークション法を、実数値最適輸送問題に適用するためのものである。この手法は有限回の反復で終了を保証し、事前誤差境界を提供する。また、実数値コストと重みを有する密行列かつ非対称な輸送問題において、従来のオークション手法と同等の性能を達成する。

ABSTRACT

Optimal transportation theory is an area of mathematics with real-world applications in fields ranging from economics to optimal control to machine learning. We propose a new algorithm for solving discrete transport (network flow) problems, based on classical auction methods. Auction methods were originally developed as an alternative to the Hungarian method for the assignment problem, so the classic auction-based algorithms solve integer-valued optimal transport by converting such problems into assignment problems. The general transport auction method we propose works directly on real-valued transport problems. Our results prove termination, bound the transport error, and relate our algorithm to the classic algorithms of Bertsekas and Castanon.

研究の動機と目的

• 従来のオークション法が有理数関係にある整数コストと重みを必要とするという制限を克服し、実数値最適輸送問題を直接解けるようにすること。
• 既存のオークションベースの手法が単位重みと整数コストに還元された割当問題に依存するというデータ制限を克服すること。
• 実数値輸送問題を扱える一般化されたオークションフレームワークを構築し、誤差の有界性や制御された近似といった主な利点を維持すること。
• 実数値最適輸送の文脈において収束性と誤差境界に関する理論的保証を確立すること。
• 密行列で非対称構造を示す実数値ネットワークフロー問題において、実用的な性能が従来のオークション手法と同等であることを実証すること。

提案手法

• 単位重みを有する割当問題に変換する必要がない、実数値輸送問題に直接作用する一般化オークション法を提案する。
• 双対変数と残差コストに基づく価格更新ルールを導入し、実行可能性を保ちつつ輸送計画の反復的改善を可能にする。
• 定理4.11を適用して、プライマルフローと双対価格ベクトルを用いた事後誤差境界を可能にし、制御された近似のためのオンラインεスケーリングを支援する。
• 解の精度と計算負荷のバランスを図るため、反復における観測誤差に応じてεを動的に調整する戦略を採用する。
• 部分輸送問題の解法としての拡張可能性を考慮し、逆一般オークションの変種を実装する。
• 半離散的最適輸送に一般的に見られる非対称問題の構造を活用し、オークション法が自然に効率的であるという特徴を活かす。

実験結果

リサーチクエスチョン

• RQ1 有理数関係にあるコストや整数重みを必要としないオークションベースの手法が、実数値最適輸送問題を解けるように拡張可能か。
• RQ2 実数値データ上で動作する一般化オークション法に対して、どのような理論的保証（例えば、終了性、誤差境界）を確立できるか。
• RQ3 密行列で実数値コストを有する輸送問題に対して、提案手法の計算複雑度と性能が、従来のオークション手法と比べてどの程度か。
• RQ4 調整可能なεスケーリングを通じて制御された近似を実現できるか、かつ解の品質を維持できるか。
• RQ5 一般化オークション法が並列計算や部分輸送問題への応用にどの程度適合できるか。

主な発見

• 一般化オークションアルゴリズムは有限時間で終了し、実数値最適輸送問題に対して事前誤差境界を提供する。
• 数値実験の結果、アルゴリズムはマシン精度に達する精度を達成しており、コストが10⁴を超える問題では、正確な解と11桁まで一致する最適コスト値を達成している。
• 時間計算量は約1.62×10⁻⁷·N².²⁸秒で、スケーリング指数b=2.28を示しており、同じ問題セットにおいてネットワークシンプレックス法のb=2.42と同等のスケーリング性能を示す。
• ストレージ計算量は、コスト値をオンデマンドで計算することで、大規模な事前計算された行列を避けるため、概ね線形に近い。
• アルゴリズムは、実数値コスト（例：pノルムから導出）と有理数関係にない重みを有する密行列かつ非対称な輸送問題を正常に解ける。従来のオークション手法では失敗する問題に対しても有効である。
• オープンソース実装が公開されており、オークション手法のバリエーションを直接比較可能であり、再現性の確保とさらなる開発を支援する。