MathLabs
定理証明済み

ブレニエの定理

内容

μ,ν\mu,\nu を Rn\mathbb R^n 上の確率測度で二次モーメントが有限であるとし、μ\mu がルベーグ測度に関して絶対連続であるとする。二次コスト c(x,y)=∣x−y∣2c(x,y) = |x-y|^2 に対して、加法定数を除いて一意な凸関数 φ:Rn→R\varphi:\mathbb R^n\to\mathbb R が存在し、T=∇φT = \nabla \varphi は μ\mu-ほとんど至るところで well-defined であり、T#μ=νT_{\#}\mu=\nu を満たし、カントロビッチ問題とモンジュの元の問題の両方の(本質的に)唯一の最小化元である。

なぜ正しいのか?

∣x−y∣2|x-y|^2 のような狭義凸のコストに対する任意の最適計画は、**cc-巡回単調**な集合の上に支えられている:有限個の対の組み替えでは総コストを下げられない。二次コストの場合、対の集合 (xi,yi)(x_i,y_i) の巡回単調性は、ある凸関数 φ\varphi が各 ii について xix_i における劣微分に yiy_i を含むという条件とまさに一致することが分かる——これは巡回単調な集合を凸関数の劣微分として特徴づけるロックアフェラーの定理である。μ\mu が絶対連続であるため、凸関数は μ\mu-ほとんど至るところで微分可能であり(アレクサンドロフの定理)、集合値の劣微分は真の勾配写像 T=∇φT=\nabla\varphi となる。

証明の概略

(概略、Brenier 1991 に従う。)ステップ1:c=∣x−y∣2c=|x-y|^2 に対するカントロビッチ双対性により、最適ポテンシャルを φ(x)=12∣x∣2−φ~(x)\varphi(x)=\tfrac12|x|^2-\tilde\varphi(x) を通じて書き直すと φ\varphi は凸となり、そのルジャンドル変換 φ∗(y)=sup⁡x(x⋅y−φ(x))\varphi^*(y)=\sup_x(x\cdot y-\varphi(x)) が第二のポテンシャルの役割を果たす;双対制約は最適な γ\gamma の台の上でちょうど y∈∂φ(x)y\in\partial\varphi(x) となる。ステップ2:ロックアフェラーの定理は、Rn×Rn\mathbb R^n\times\mathbb R^n の任意の巡回単調な部分集合が、ある凸で下半連続な φ\varphi の ∂φ\partial\varphi のグラフに含まれることを述べる;最適計画の台が巡回単調であるのは、γ\gamma が ∫∣x−y∣2 dγ\int|x-y|^2\,d\gamma を最小化するため、台の上の有限個の対の間でどのように組み替えても ∑i∣xi−yi∣2\sum_i|x_i-y_i|^2 を厳密に減らせないからである。ステップ3:μ≪\mu\ll ルベーグ測度であるため、アレクサンドロフの定理により φ\varphi は μ\mu-ほとんど至るところで微分可能であり、したがって γ\gamma はグラフ {(x,∇φ(x))}\{(x,\nabla\varphi(x))\} 上に集中する、すなわち γ=(id,∇φ)#μ\gamma=(\mathrm{id},\nabla\varphi)_{\#}\mu となる。これにより決定論的な最適写像 T=∇φT=\nabla\varphi が得られ、それはモンジュの問題も解くことになり、コストの狭義凸性により他のどんな最適写像も μ\mu-ほとんど至るところでこれと一致しなければならない。

この定理を使うトピック

ステップごとの証明

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

参考文献

  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