高速ソートと安定結婚問題

高速ソートと安定結婚問題 平成24年12月14日

高速ソート

単純なソートアルゴリズム 　　バブルソート　　Ｏ（ｎ２）　　　　選択ソート　　　Ｏ（ｎ２）　　　　挿入ソート　　　Ｏ（ｎ２）　　もっと高速なソートアルゴリズムは？計算量がＯ（ｎｌｏｇｎ）のソートアルゴリズムクイックソート　　ヒープソート　　マージソートこの中でクィックソートが最速（ランダムデータに対して．不得手な場合もある）

クイックソートの原理 考え方：分割統治法（ｄｉｖｉｄｅａｎｄｃｏｎｑｕｅｒ） • まず、配列中の適当な要素を選ぶ　　　→　これを枢軸（ｐｉｖｏｔ）と呼ぶ • 次のように要素を並べ替える　　　→　これを分割（ｐａｒｔｉｔｉｏｎ）というａ［０］　　・・・・・・・・・・　ａ［ｖ－１］ａ［ｖ］ａ［ｖ＋１］　・・・　ａ［ｎ－１］ａ［ｖ］より小（以下）ａ［ｖ］ａ［ｖ］より大（以上）

では、ａ［０］,・・・,ａ［ｖ－１］の部分と、　　　　ａ［ｖ＋１］,・・・,ａ［ｎ－１］の部分はどうするか？では、ａ［０］,・・・,ａ［ｖ－１］の部分と、　　　　ａ［ｖ＋１］,・・・,ａ［ｎ－１］の部分はどうするか？ ↓ ピボットの選択 →分割を再帰的に繰り返せばよい！

クィックソートのアルゴリズム概観 ①ソートするデータ区間のデータを得る ②区間内のデータ数が2以下の時2-1：データ数が1以下の時なにもしない2-2：データ数が２の時，必要なら交換 ③区間内のデータ数が3以上の時3-1:区間内でピボットを選ぶ（便宜的に右端）3-2:区間内で，ピボットより小さい区間（左半分）と大きい区間(右半分）に分割する3-3:3-2の両区間を再帰的処理する

/* 配列ａのうち a[l]～a[r] を整列する */ quick_sort ( int a[ ], int l, int r ) { if( 整列する要素が一つのみである ) return; 　　適当な要素 a[v] を枢軸にして、 a[v] より小さい要素を　a[l]～a[v-1] に集め、 a[v] より大きい要素を　a[v+1]～a[r] に集める quick_sort ( a, l, v-1 );/* 左の部分を整列する */ quick_sort ( a, v+1, r );/* 右の部分を整列する */ }

分割のアルゴリズム ０ａ．一番右の要素を枢軸に選ぶ０ｂ．ポインタｉを左端に、ポインタｊを枢軸の　　　すぐ左に設定する１．枢軸より大きな要素が見つかるまで　　ｉを右へスキャン２．枢軸より小さな要素が見つかるまで　　ｊを左へスキャン３．ｉが指す要素とｊが指す要素を交換する４．ｉとｊがぶつかるまで１～３を繰り返す５．ｊの指す要素と枢軸と交換する

ピボット i j ５５３７４２０１３８７４６３０前半に＜３０，後半に＞３０を置く．それに反するものを入れ替える j ５５３７４２０１３８７４６３０ i j １３３７４２０５５８７４６３０ j i １３３７４２０５５８７４６３０ j i １３３２０７４５５８７４６３０左スキャンと右スキャンが交錯したらピボットをそこに置く→＠ピボットの左にピボットより小さい値，右に大きい値があるはず j ｖ i １３３２０３０５５８７４６７４

再帰版クイックソート /* 配列a[l]~a[r]を分割する。枢軸の添え字を返す */ int partition( int a[], int l, int r ) { int i, j, pivot, t; i = l – 1; j = r; /* ポインタ i とｊを初期化する */ pivot = a[r]； /* いちばん右端の要素を枢軸にする */ for( ; ; ) { /* ポインタ i と j とがぶつかるまで繰り返す */ while( a[++i] < pivot ) ; /* ポインタ i を右へ進める */ while( i < j && pivot < a[j] ) ; /* ポインタ i を左へ進める */ if ( i >= j ) break; /* ポインタ i と j がぶつかったらループを抜ける */ t = a[i]: a［i］ = a[j]; a[j] = t; /* i の指す要素とｊの指す要素を交換する */ } t = a[i]: a［i］ = a[r]; a[r] = t; /* a[i] と枢軸を交換する */ return i; ｝

クイックソートの計算量 最適な分割：　ｐｉｖｏｔがほぼ中央に来る　→　ほぼ半分ずつに分割される　→　要素1までの分割数をｐとすると２ｐ =n ，よってp=log 2n回分割が必要　　　部分列の個数はｎ個．　　 →計算量は　Ｏ（ｎｌｏｇｎ）

最悪の分割：ピボットが端に来る　　枢軸枢軸より小さい要素ｋ－１個枢軸　　　↓ 分割のレベル数（深さ）はｎ枢軸より大きい要素０個　　　↓ 計算量は　Ｏ（ｎ２）まとめ：　クイックソートの計算量は　　最良　Ｏ（ｎｌｏｇｎ）　　平均も！　　最悪　Ｏ（ｎ２）

クイックソートの改善方法 • 部分配列の長さが１０程度以下になったら、　挿入ソートに切り替えるｋ１・ｎ２　＞？　ｋ２・ｎｌｏｇ２ｎ　　（ｋ１＜ｋ２）俗にquicker sort • 最悪の場合を回避する　部分配列の右端、中央、左端の３つのうちの　　中央値を枢軸に

再帰の深さを減らす →　左右の部分配列のうち、　　　　短いほうから先に処理をする　　　　　　→　再帰の深さは最大ｌｏｇ２ｎ程度 →　非再帰版　　　　　　→　再帰の処理をスタックを利用して　　　　　　　　　自分で行なう。

安定結婚問題

１２３４５ ① ② ③ ④ ⑤ 安定結婚問題フィーリングカップル5対5　（講義1回目）

安定結婚問題 ・最初の論文　→　[Gale & Shapley 1962] - アメリカの研修医配属問題がきっかけ。 - どんな例題にも、必ず安定マッチングが存在する。 - 安定マッチングを多項式時間で見つけることができる。 (Gale-Shapley アルゴリズム) ・様々な類似問題 - 安定ルームメイト問題 - Residents/Hospitals 問題・最近、様々な新種の問題・実世界での応用研修医配属 NRMP (National Resident Matching Program) CaRMS (Canadian Resident Matching Service) SPA (Scottish Pre-registration house officer Allocations) JRMP (Japanese Resident Matching Program)

不安定結婚とは • カップルでない男女双方ともに互いに対する好感度が自分のカップル相手よりも高い場合、そのカップルは持続せず解消してしまう。→不安定なカップル • 上記の状況が全く発生していない→安定なカップル

力任せで安定結婚問題を解くアルゴリズムコンピュータは速いから、すべての5組カップルを生成して（５！＝５＊４＊３＊２＊１＝１２０）、それが安定か不安定かを調べればいい？力任せで安定結婚問題を解くアルゴリズムコンピュータは速いから、すべての5組カップルを生成して（５！＝５＊４＊３＊２＊１＝１２０）、それが安定か不安定かを調べればいい？ N=111秒 N=１２　　　　　　４秒 N=１３　 2分52秒 N=１４　４2分３０秒 N=１５　１１時間・すべての組み合わせを作る計算量は　Ｎ！・安定性をチェックする計算量は　　　　　Ｎ＊Ｎ・よって総計算量は　　　　　　　　　　　　　Ｎ！＊Ｎ＊Ｎ・N=16のときの大よその実行時間は　　1週間

もっと速いアルゴリズムは？ • なぜ遅い？⇒すべての組合せを生成して，毎回浮気のチェックをしているので，　　入力データ数Nの壁に阻まれる生成検査法（Generate and Test)⇒一発で解が作れないか？⇒それは無理⇒でも無意味な組み合わせから始めるのはよくない⇒どのような組合せなら作成可能で無駄が少ないか？ ①理想的な状態（浮気が起こらない組合せ）から始めてみてはどうか？理想状態＝第一希望の異性とのカップル　　相手にパートナーがいるか否かを無視して、とりあえず強引に第一希望に申し込む。既にパートナーがいれば、その好感度を比較して、既存パートナーの好感度が高ければ、申込を破棄、低ければ、パートナーを解消して、新パートナーを作る。

もっと速いアルゴリズムは？ ②乗り換え発生時：カップル解消された人は困りますよね。その人は、一つ好感度が低い人に新たにカップルを申し入れて、先程と同様の処理をする。 ③カップル解消が起きて一人はじかれると、入れ替え（玉突き現象）が際限なく起こらない？ ⇒これは大丈夫。好感度が一つ低い相手を選ぶので、一度のカップル解消で最大入れ替え回数は？回

④⑤③①② ②③④①⑤ ④③②⑤① ③②⑤①④ ④②⑤③① ４１３２５１３５２４３５１２４４２１３５４３２１５下記の例で効率的アルゴリズムを実行してみると男子 × 女子

④⑤③①② ②③④①⑤ ④③②⑤① ③②⑤①④ ④②⑤③① ４１３２５１３５２４３５１２４４２１３５４３２１５フィーリングカップル5対5　の解答男子手順１男性１：④ 手順２男性１：④ 男性２：③ 手順３男性１：④ 男性２：③ 男性３：④→② 手順４男性１：④→③ 男性２：③→②→⑤ 男性３：④→② 男性４：④ 手順５男性１：④→③→② 男性２：③→②→⑤→① 男性３：④→②→⑤ 男性４：④ 男性５：②→③ × 女子

男性の希望リスト　　　　　　　　女性の希望リスト男性の希望リスト　　　　　　　　女性の希望リスト • 1: a c b d e a: 2 1 3 4 5 • 2: c a e b d b: 2 1 4 5 3 • 3: b a e d c c: 1 2 3 5 4 • 4: c b d e a d: 3 1 4 2 5 • 5: c d b e a e: 4 3 1 2 5 より形式的に安定結婚問題を考えると入力：男性N人　　　女性N人希望リスト N=5の例男性： 1,2,3,4,5 女性： a,b,c,d,e

1: a c b d e a: 2 1 3 4 5 2: c a e b d b: 2 1 4 5 3 3: b a e d c c: 1 2 3 5 4 4: c b d e a d: 3 1 4 2 5 5: c d b e a e: 4 3 1 2 5 出力：男女間の不安定マッチングの例 1 a 2b 3c 4d 5e 男性 1と女性 cはブロッキングペアブロッキングペアの存在しないマッチング: 安定マッチングこれは解の条件を満たしていない。安定結婚問題:　与えられた入力から、安定マッチングを見つける。

Theorem [Gale & Shapley 1962, Gusfield & Irving 1989] The Gale-Shapley algorithm finds a stable matching. The Gale-Shapley Algorithm [Gale & Shapley 1962] (Men-propose） × 1: a c b d e a: 2 1 3 4 5 2: c a e b d b: 2 1 4 5 3 3: b a e d c c: 1 2 3 5 4 4: c b d e a d: 3 1 4 2 5 5: c d b e a e: 4 3 1 2 5 × × × × × × ×

高速ソートと 安定結婚問題