高速ソート-yamaguti.comp.ae.keio.ac.jp.pptVIP

* * * * * * * * * * * * * * * * * * * * * 高速ソートと安定結婚問題 平成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 ５５ ３ ７４ ２０ １３ ８７ ４６ ３０ i j １３ ３ ７４ ２０ ５５ ８７ ４６ ３０ j １３ ３ ７４ ２０ ５５ ８７ ４６ ３０ i j １３ ３ ２０ ７４ ５５ ８７ ４６ ３０ １３ ３ ２０ ３０ ５５ ８７ ４６ ７４ i i j j ｖ 前半に＜３０，後半に＞３０を置く．それに反するものを入れ替える 左スキャンと右スキャンが交錯したらピボットをそこに置く→＠ピボットの左にピボットより小さい値，右に大きい値があるはず 再帰版クイックソート /* 配列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 の指す要素と ｊ の指す要素を交換する */