MathLabs
定理証明済み

カントロビッチ双対定理

内容

XX 上の確率測度 μ\mu、YY(ともにポーランド空間)上の確率測度 ν\nu、および連続で下に有界な費用 cc に対して、カントロビッチ問題は、各点で φ(x)+ψ(y)≤c(x,y)\varphi(x) + \psi(y) \le c(x,y) を満たすカントロビッチ・ポテンシャルの組 φ:X→R\varphi:X\to\mathbb R、ψ:Y→R\psi:Y\to\mathbb R にわたる双対最大化問題に等しい:min⁡γ∈Π(μ,ν)∫X×Yc dγ  =  sup⁡φ(x)+ψ(y)≤c(x,y)(∫Xφ dμ+∫Yψ dν),\min_{\gamma\in\Pi(\mu,\nu)}\int_{X\times Y} c\,d\gamma \;=\; \sup_{\varphi(x)+\psi(y)\le c(x,y)}\left(\int_X\varphi\,d\mu+\int_Y\psi\,d\nu\right), 最小値と上限はともに達成される。

なぜ正しいのか?

φ(x)\varphi(x) を xx で質量を引き取る際に課される価格、ψ(y)\psi(y) を yy で引き渡す際に支払われる価格とみなし、独立した運送業者たちによる分散市場を考えてみよう。どの業者も直接輸送のルートを出し抜いて利益を得ることはできない:φ(x)\varphi(x) で買い ψ(y)\psi(y) で売る取引は、実際に質量を運んで c(x,y)c(x,y) を支払うことに決して勝てない——これがまさに制約 φ(x)+ψ(y)≤c(x,y)\varphi(x)+\psi(y)\le c(x,y) である。市場が最適に清算されるとき、等式 φ(x)+ψ(y)=c(x,y)\varphi(x)+\psi(y)=c(x,y) は最適な計画が実際に使う経路上でちょうど成り立つ。

証明の概略

有限離散測度 μ=∑ipiδxi\mu=\sum_i p_i\delta_{x_i}、ν=∑jqjδyj\nu=\sum_j q_j\delta_{y_j} に対して、カントロビッチ問題はまさに有限線形計画 min⁡∑i,jcijγij\min \sum_{i,j} c_{ij}\gamma_{ij}(制約 γij≥0\gamma_{ij}\ge0、∑jγij=pi\sum_j\gamma_{ij}=p_i、∑iγij=qj\sum_i\gamma_{ij}=q_j)である。その線形計画双対は、すべての i,ji,j について φi+ψj≤cij\varphi_i+\psi_j\le c_{ij} を満たす下で max⁡∑ipiφi+∑jqjψj\max \sum_i p_i\varphi_i+\sum_j q_j\psi_j を求める問題に正確に一致する——原問題の各周辺制約につき1つの双対変数が対応する。原問題は実行可能であり(直積カップリング piqjp_iq_j が常に条件を満たす)、00 で下に有界であり、実行可能多面体 {γij≥0}∩Π(μ,ν)\{\gamma_{ij}\ge0\}\cap\Pi(\mu,\nu) はコンパクトであるから、線形計画の強双対性により2つの最適値は等しく、両方とも達成される。これにより離散の場合に定理が厳密に証明される。ポーランド空間上の一般的な主張は、同じ主・双対の枠組みを Cb(X×Y)C_b(X\times Y) 上の凸汎関数に対するフェンシェル・ロックアフェラー双対として適用することで従う(カントロビッチ、1942年;完全な議論は Villani 2009, 定理5.10 を参照)。

この定理を使うトピック

ステップごとの証明

この定理のステップごとの証明はまだありません。

参考文献

  1. Cédric Villani (2009). Optimal Transport: Old and New
  2. Marco Cuturi (2013). Sinkhorn Distances: Lightspeed Computation of Optimal Transport · arXiv:1306.0895
  3. Yann Brenier (1991). Polar Factorization and Monotone Rearrangement of Vector-Valued Functions · DOI:10.1002/cpa.3160440402
  4. John Lott, Cédric Villani (2009). Ricci Curvature for Metric-Measure Spaces via Optimal Transport