《離散凸解析》のソースを表示

'''【りさんとつかいせき (discrete convex analysis) 】'''

　離散的な集合(例えば整数格子点の集合)の上で定義された関数の構造を，凸解析 ([6]) とマトロイド理論 ([1, 7, 8]) の両方の視点から考察する方法論を，[[離散凸解析]] (discrete convex analysis) ([1, 3, 4, 5])と呼ぶ．より一般的には，解析的な視点と組合せ論的な視点の両方から「組合せ論的な凸性」という構造を考察する方法論を指す．離散最適化 ([2])，システム解析 ([5])，数理経済学 ([5, 7]) などへの応用がある．

　整数格子点上で定義され整数値をとる関数<math>f: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ \pm\infty \}\, </math>を考える(<math>V\, </math>は有限集合である)．<math>{{\rm dom\,}} f = \{ x \in {\mathbf Z}^{V} \mid -\infty < f(x) < +\infty \}\, </math>を<math>f\, </math>の実効定義域と呼び，以下では，<math>{{\rm dom\,}} f \not= \emptyset\, </math>であるような関数だけを考える．<math>i \in V\, </math>に対してその特性ベクトルを<math>\chi_{i} \ (\in \{ 0,1 \}^{V})\, </math> と表わす．ベクトル<math>x \in {\mathbf Z}^{V}\, </math>に対して，<math>{{\rm supp}^{+}}(x) = \{ i \in V \mid x_{i} > 0 \}\, </math>, <math>{{\rm supp}^{-}}(x) = \{ i \in V \mid x_{i} < 0 \}\, </math>とおく．関数<math>f: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ +\infty \}\, </math>が交換公理:

<table align="center">
<tr>
<td>任意の <math>x, y \in {{\rm dom\,}} f\, </math> と任意の <math>i \in {{\rm supp}^{+}}(x-y)\, </math> に対して, ある<br> <math>j \in {{\rm supp}^{-}}(x-y)\, </math> が存在して
</td>
</tr>
</table>


<center>
<math>f(x)+f(y) \geq f(x-\chi_{i}+\chi_{j}) + f(y+\chi_{i}-\chi_{j})\, </math> 
</center>


を満たすとき，[[M凸関数]] (M-convex function) という．M凸関数<math>f\, </math>の実効定義域<math>{{\rm dom\,}} f\, </math>は整基多面体(に含まれる整数格子点)である．

　関数g:<math> {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ +\infty \}\, </math>が2条件:
　

<table align="center">
<tr>
<td><math>g(p) + g(q) \geq g(p \vee q) + g(p \wedge q)
\qquad ( p, q \in {\mathbf Z}^{V})  ,\, </math>
</td>
</tr>
<tr>
<td><math>\exists r \in {\mathbf Z}, \forall p \in {\mathbf Z}^{V}: \ 
 g(p+{\mathbf 1}) = g(p) + r  ,\, </math>
</td>
</tr>
</table>


を満たすとき，[[L凸関数]] (L-convex function) という．ここで，<math>p \vee q, p \wedge q\, </math>は，それぞれ，成分毎に最大値, 最小値をとって得られるベクトル(すなわち，<math>(p \vee q)_{i} = \max(p_{i}, q_{i})\, </math>, <math>(p \wedge q)_{i} = \min(p_{i}, q_{i}))\, </math>を表し，<math>{\mathbf 1}=(1,1,\ldots,1) \in {\mathbf Z}^{V}\, </math>である．L凸関数<math>g\, </math>の実効定義域<math>{{\rm dom\,}} g\, </math>は<math>\vee\, </math>, <math>\wedge\, </math>に関して<math>{\mathbf Z}^{V}\, </math>の部分束を成す．また，正斉次L凸関数は，劣モジュラ集合関数と同一視することができる．

　一般に関数<math>h: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ \pm\infty \}\, </math> の凸共役<math>h^{\bullet}\, </math>，凹共役<math>h^{\circ}\, </math>を


<center>
<math>\begin{array}{lll}
 h^{\bullet}(p) 
 &=& \sup\{  \langle p, x \rangle - h(x) \mid x \in {\mathbf Z}^{V} \}
\qquad ( p \in {\mathbf Z}^{V}) ,
\\
 h^{\circ}(p) 
 &=& \inf\{  \langle p, x \rangle - h(x) \mid x \in {\mathbf Z}^{V} \}
\qquad ( p \in {\mathbf Z}^{V})
\end{array}
\, </math>
</center>


と定義する．ここで，<math>\textstyle \langle p, x \rangle = \sum_{i \in V} p_{i}x_{i}\, </math>である．この対応<math>h \mapsto h^{\bullet}\, </math>, <math>h \mapsto h^{\circ}\, </math> を(凸，凹)離散フェンシェル・ルジャンドル(Fenchel-Legendre)変換と呼ぶ．M凸関数とL凸関数は離散フェンシェル・ルジャンドル変換に関して共役関係にあり，対応<math>f \mapsto f^{\bullet}\, </math>, <math>g \mapsto g^{\bullet}\, </math>はM凸関数<math>f\, </math>とL凸関数<math>g\, </math>の間の１対１対応を与える．すなわち，M凸関数<math>f\, </math>とL凸関数<math>g\, </math>に対して，<math>f^{\bullet}\, </math>はL凸関数, <math>g^{\bullet}\, </math>はM凸関数で，<math>(f^{\bullet})^{\bullet}=f\, </math>, <math>(g^{\bullet})^{\bullet}=g\, </math>が成り立つ(共役性定理)．

　M凸関数やL凸関数に対して，[[離散分離定理]] (discrete separation theorem) や[[フェンシェル型双対定理]] (Fenchel-type duality theorem) に象徴されるような離散双対性が成り立つ．

　M凸関数に関する離散分離定理(M分離定理)を述べる:


:[M分離定理] <math>f: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ +\infty \}\, </math> をM凸関数，<math>g: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ -\infty \}\, </math>をM凹関数とし, <math>{{\rm dom\,}} f \cap {{\rm dom\,}} g \not= \emptyset\, </math>または<math>{{\rm dom\,}} f^{\bullet} \cap {{\rm dom\,}} g^{\circ} \not= \emptyset\, </math>であると仮定する．このとき, <math>f(x) \geq g(x)\, </math> <math>(\forall \ x \in {\mathbf Z}^{V})\, </math>ならば, ある<math>\alpha \in {\mathbf Z}\, </math>, <math>p \in {\mathbf Z}^{V}\, </math>が存在して 


<center>
<math>f(x) \geq \alpha + \langle p, x \rangle  \geq g(x)   \qquad  (\forall \ x \in {\mathbf Z}^{V})\, </math>
</center>


　が成り立つ(<math>g\, </math>がM凹関数とは<math>-g\, </math>がM凸関数のことである)．

ここで，<math>p\, </math>が整数ベクトルに選べることが離散性の反映である．上の主張で，<math>f\, </math>をL凸関数，<math>g\, </math>をL凹関数に置き換えたものも成立する(L分離定理)．L分離定理は，その特殊ケースとして，劣モジュラ集合関数に関するA. Frankの離散分離定理 ([1] 参照) を含んでいる．

　M分離定理とL分離定理は互いに共役の関係にあるが，次に述べるフェンシェル型双対定理は自己共役の形になっている．


:[フェンシェル型双対定理]   <math>f: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ +\infty \}\, </math> をM凸関数，<math>g: {\mathbf Z}^{V} \to {\mathbf Z} \cup \{ -\infty \}\, </math>をM凹関数とし, <math>{{\rm dom\,}} f \cap {{\rm dom\,}} g \not= \emptyset\, </math>または<math>{{\rm dom\,}} f^{\bullet} \cap {{\rm dom\,}} g^{\circ} \not= \emptyset\, </math>であると仮定する．このとき, 


<center>
<math>\inf\{ f(x) - g(x) \mid x \in {\mathbf Z}^{V}  \}
 = \sup\{ g^{\circ}(p) - f^{\bullet}(p) \mid p \in {\mathbf Z}^{V} \}\, </math>
</center>


　が成り立つ．さらに，この両辺が有限値なら，<math>\inf\, </math>を達成する<math>x \in {{\rm dom\,}} f \cap {{\rm dom\,}} g</math>と<math>\sup\, </math>を達成する<math>p \in {{\rm dom\,}} f^{\bullet} \cap {{\rm dom\,}} g^{\circ}\, </math>が存在する．

上の定理は，非線形整数計画に関する強双対性を主張しており，その本質的な部分は，<math>\inf\, </math>と<math>\sup\, </math>をとる範囲をそれぞれ整数ベクトルに限ってよいという主張にある．

　M凸関数，L凸関数に関する種々の問題に対して効率的なアルゴリズムが開発されている．これに関しては [4, 5] を参照されたい．



----
'''参考文献'''

[1] S. Fujishige, ''Submodular Functions and Optimization'',2nd ed.,  Elsevier, 2005.

[2] N. Katoh and T. Ibaraki, "Resource allocation problems," in ''Handbook of Combinatorial Optimization, Vol.2'', D. -Z. Du and P. M. Pardalos, eds., Kluwer, 159-260, 1998.

[3] K. Murota, "Discrete convex analysis", ''Mathematical Programming'', '''83''' (1998), 313-371.

[4] 室田一雄, 「離散凸解析」, 『共立出版』，2001.

[5] K.Murota, ''Discrete Convex Analysis'', SIAM, 2003.

[6] R. T. Rockafellar, ''Convex Analysis'', Princeton University Press, 1970. 

[7] A.Tamura, "Applications of discrete convex analysis to mathmatical economics", ''Publ. RIMS'', '''40''' (2004) , 1015-1037.

[8] D. J. A. Welsh, ''Matroid Theory''，Academic Press, 1976.

[9] N. White, ed., ''Theory of Matroids'', Cambridge University Press, 1986.