近似解法T-@niftyホームページサービス-ホームページ…可研报告.ppt

半正定値計画を用いた最大カット問題の .878 近似解法ver. C.22 東京大学大学院工学系研究科 計数工学専攻 松井知己 文献紹介 Improved Approximation Algorithms for Maximum Cut and Satisfiability Problems Using Semidefinite Programming, M. X. Goemans and D. P. Williamson, J. ACM, 42 (1995), pp.1115-1145. Approximation Algorithm ：近似解法 Maximum Cut Problem ：最大カット問題 Maximum Satisfiability Problem：最大充足可能性問題？ Semidefinite Programming ：半正定値計画 MAX CUT問題(問題例) カット重み=８　　　　　　　カット重み=１４ ???の定義 グラフG=（V,E) V：頂点集合 E：枝集合 MAX CUT問題(定義) MAX CUT問題 入力:無向グラフG=(V,E)，非負整数枝重みw:E→Z+． 問題：G中のカットで， 　　 重み最大となるものを求めよ． 頂点部分集合U?V に対するカット： δ(U)=｛e∈E |e はU中の頂点とU に無い頂点を結ぶ｝ カットδ(U)の重さ： 　w(U)= δ(U)中の枝重みの総和． MAX CUT問題：max { w(U ) | U?V } 　　　　　　　　　　　　　　　　　　NP‐困難[ Karp 72] MAX CUT 問題の重要性 0-1変数2次計画問題(制約無し） max{ xtQx + cｔx | x ∈{0,1}n } (Q:対称行列） 近似解法とは？ 近似解法：近似精度の存在する解法． 得られる解と最適解との距離の見積もりが存在． 発見的解法：近似精度は不明． 近似比率：　（最大化問題の場合） ＝（得られた解の目的関数値）／（最適値） α近似解法： 近似比率α以上の解を必ず求める解法． MAX CUTの1/2 近似解法 ●非常に単純なランダム化算法． １/２近似解法 MAX CUT問題：max { w(U) | U ?V } １/２近似解法： 各頂点を１/２の確率でUに入れる． 解析：各枝について，両端の内 　　ちょうど一方がUに入る確率は1/2． カットの重みの期待値ｚ は， ｚ　＝∑e ∈E w（e ） Pr[枝e 　がカットに入る] ＝∑e ∈E w（e ）(1/2) ＝(1/2) (枝重みの総和) ≧ (1/2)(最大カット重み) =(1/2)(最適値)． MAX CUTの近似解法の歴史 (1/2） [1976: Sahni and Gonzales] (1/2)+(1/2m) [1981: Vitany] (1/2)+((n -1)/4m)[1982: Poljak and Turzik] (1/2)+(1/2n) [1991: Haglin and Venkatesan] (1/2)+(1/2D) [1995: Hofmeister and Lefmann] n :頂点数 m :枝数 D :最大次数 （次数：頂点につながる枝数） 0.878 [1995: Goemans and Williamson] MAX CUTの0.878 近似解法 [Goemans and Williamson 95] ●MAX CUTを 　　非凸2次連続最適化問題に変形． ●半正定値計画緩和． ●ランダム化アルゴリズム． 球面上の配置問題（図） グラフGの頂点を球面S上に配置して， 　　枝重みと距離の積の総和を最大化する． max{∑w(i,j) ｄ(v(i )v( j)) | v(i )∈ S (? i ∈V )} 球面上の配置問題 入力:無向グラフG=(V,E)，非負枝重みw:E→ Z+， S：原点を中心とし，半径1/π の球の表面． 次元は任意． (大円の半円弧の長さ=１)． ?v,?v’∈ S , 　　ｄ(vv’)= v と v’を結ぶ大円弧の短い方の長さ． 球面上の配置問題 グラフGの頂点をS上に配置して， 　　枝重みと距離の積の総和を最大化する． max{∑w(i,j) ｄ(v(i )v( j)) | v(i )∈ S (?