归并排序

转自： http://blog.csdn.net/morewindows/article/details/6678165/

http://blog.csdn.net/u014082714/article/details/44132111

归并排序是建立在归并操作上的一种有效的排序算法。该算法是采用分治法（Divide and Conquer）的一个非常典型的应用。

首先考虑下如何将将二个有序数列合并。这个非常简单，只要从比较二个数列的第一个数，谁小就先取谁，取了后就在对应数列中删除这个数。然后再进行比较，如果有数列为空，那直接将另一个数列的数据依次取出即可。

[cpp]  view plain  copy

1. //将有序数组a[]和b[]合并到c[]中  
2. void MemeryArray(int a[], int n, int b[], int m, int c[])  
3. {  
4.     int i, j, k;  
5.   
6.     i = j = k = 0;  
7.     while (i < n && j < m)  
8.     {  
9.         if (a[i] < b[j])  
10.             c[k++] = a[i++];  
11.         else  
12.             c[k++] = b[j++];   
13.     }  
14.   
15.     while (i < n)  
16.         c[k++] = a[i++];  
17.   
18.     while (j < m)  
19.         c[k++] = b[j++];  
20. }  

可以看出合并有序数列的效率是比较高的，可以达到O(n)。

解决了上面的合并有序数列问题，再来看归并排序，其的基本思路就是将数组分成二组A，B，如果这二组组内的数据都是有序的，那么就可以很方便的将这二组数据进行排序。如何让这二组组内数据有序了？

可以将A，B组各自再分成二组。依次类推，当分出来的小组只有一个数据时，可以认为这个小组组内已经达到了有序，然后再合并相邻的二个小组就可以了。这样通过先递归的分解数列，再合并数列就完成了归并排序。

[cpp]  view plain  copy

1. //将有二个有序数列a[first...mid]和a[mid...last]合并。  
2. void mergearray(int a[], int first, int mid, int last, int temp[])  
3. {  
4.     int i = first, j = mid + 1;  
5.     int m = mid,   n = last;  
6.     int k = 0;  
7.       
8.     while (i <= m && j <= n)  
9.     {  
10.         if (a[i] <= a[j])  
11.             temp[k++] = a[i++];  
12.         else  
13.             temp[k++] = a[j++];  
14.     }  
15.       
16.     while (i <= m)  
17.         temp[k++] = a[i++];  
18.       
19.     while (j <= n)  
20.         temp[k++] = a[j++];  
21.       
22.     for (i = 0; i < k; i++)  
23.         a[first + i] = temp[i];  
24. }  
25. void mergesort(int a[], int first, int last, int temp[])  
26. {  
27.     if (first < last)  
28.     {  
29.         int mid = (first + last) / 2;  
30.         mergesort(a, first, mid, temp);    //左边有序  
31.         mergesort(a, mid + 1, last, temp); //右边有序  
32.         mergearray(a, first, mid, last, temp); //再将二个有序数列合并  
33.     }  
34. }  
35.   
36. bool MergeSort(int a[], int n)  
37. {  
38.     int *p = new int[n];  
39.     if (p == NULL)  
40.         return false;  
41.     mergesort(a, 0, n - 1, p);  
42.     delete[] p;  
43.     return true;  
44. }  

 

归并排序的效率是比较高的，设数列长为N，将数列分开成小数列一共要logN步，每步都是一个合并有序数列的过程，时间复杂度可以记为O(N)，故一共为O(N*logN)。因为归并排序每次都是在相邻的数据中进行操作，所以归并排序在O(N*logN)的几种排序方法（快速排序，归并排序，希尔排序，堆排序）也是效率比较高的。

 

在本人电脑上对冒泡排序，直接插入排序，归并排序及直接使用系统的qsort()进行比较（均在Release版本下）

对20000个随机数据进行测试：

对50000个随机数据进行测试：

再对200000个随机数据进行测试：

 

注：有的书上是在mergearray()合并有序数列时分配临时数组，但是过多的new操作会非常费时。因此作了下小小的变化。只在MergeSort()中new一个临时数组。后面的操作都共用这一个临时数组。

-----------------------------

归并排序的核心思想是将两个已经排序的序列合并成一个序列，那如何得到两个已经排序的序列呢？我们知道， 如果一个序列只有一个元素，那该序列是已经排序的，这样我们就可以利用分治的思想，将未排序的序列划分成更小的序列，只到我们可以很方便的对小序列进行排序（比如划分到序列只有一个元素， 或者序列很小可以方便的使用其它排序算法进行排序），然后再将小序列逐次合并，得到最后的排序结果。

[cpp]  view plain  copy

 print ?

1. <span style="font-size:18px;">#include<stdio.h>  
2. #include<stdlib.h>  
3. #include<string.h>  
4. int num =0;  
5. void merge(int in[], int out[], int l, int m, int r)  
6. {  
7.     int i = l;  
8.     int k  = m + 1;  
9.   
10.     while (l <= m && k <= r)   
11.     {  
12.         if (in[l] <= in[k]) {  
13.         out[i++] = in[l++];  
14.     }   
15.     else   
16.     {  
17.         out[i++] = in[k++];  
18.         /* 因为此时in[l...m]已经排序，如果in[l] 与 in[k] 逆序, 
19.        则in[l+1], in[l+2]...in[m] 都与in[k]逆序, 共 m-l+1 对*/  
20.         num += m - l + 1; /* calculate the inversion number */  
21.     }  
22. }  
23.     while (l <= m)   
24.     {  
25.         out[i++] = in[l++];  
26.     }  
27.     while (k <= r)  
28.     {  
29.         out[i++] = in[k++];  
30.     }  
31. }  
32. /* 
33. * @brief 递归将序列划分为只有一个元素的子序列, 然后逐次对子序列进行合并 
34. */  
35. void m_sort(int in[], int out[], int l, int r)  
36. {  
37. /* 仅有一个元素, 已排序, 递归结束 */  
38.     if (l >= r)   
39.     {  
40.         return;  
41.     }     
42. /* 计算 l 和 r 的中间值, 防止溢出 */  
43.     int m = (l & r) + ((l ^ r) >> 1);  
44. /* note that in and out are swapped */  
45.     m_sort(out, in, l, m);  
46.     m_sort(out, in, m + 1, r);  
47.   
48.     merge(in, out, l, m, r);  
49. }  
50. /* 
51. * @brief merge sort 
52. * 统一申请空间, 避免反复申请释放 
53. */  
54. int merge_sort(int a[], int n)  
55. {  
56.     int *b = (int *)malloc(n * sizeof(int));  
57.     if (b)   
58.     {  
59.         memcpy(b, a, n * sizeof(int));  
60.         m_sort(b, a, 0, n - 1);  
61.         free(b);  
62.         return 0;  
63.     }  
64.   
65.     return -1;  
66. }  
67. int main()  
68. {  
69.     int a[9]={2,1,3,8,5,7,4,6,10};  
70.     merge_sort(a,9);  
71.     for(int i=0;i<9;i++)  
72.     printf("%d   ",a[i]);  
73.     printf("\n");  
74.     return 0;  
75. }</span>  

这是一个递归算法，这个算法的理解其实可以借助下面这个图：

对于原始的数组2,1,3,8,5,7,6,4,10，在整个过程执行的是顺序是途中红色编号1-20。虽然我们描述中说的是程序先分解，再归并，但实际过程是一边分解一边归并，前半部分分先排好序，后半部分再拍好，最后整个归并为一个完整的序列，途中的merge过程它所在层的两个序列的merge过程：下图展示了每个merge过程对作用于数组的哪部分（红色）。

整个过程就像一个动态的树，执行顺序就是对树的先序遍历顺序。

 

最后实验验证，对10000个倒序的数组进行排序，直接插入排序需1024ms，二归并排序只需20ms。