分治法(归并排序为例)思想
什么是分治法:
在计算机科学中,分治法是建基于多项分支递归的一种很重要的算法范式。字面上的解释是“分而治之”,就是把一个复杂的问题分成两个或更多的相同或相似的子问题,直到最
后子问题可以简单的直接求解,原问题的解即子问题的解的合并。
这个技巧是很多高效算法的基础,如排序算法(快速排序、归并排序)、傅立叶变换(快速傅立叶变换)。
另一方面,理解及设计分治法算法的能力需要一定时间去掌握。正如以归纳法去证明一个理论,为了使递归能够推行,很多时候需要用一个较为概括或复杂的问题去取代原有
问题。而且并没有一个系统性的方法去适当地概括问题。
分治法这个名称有时亦会用于将问题简化为只有一个细问题的算法,例如用于在已排序的列中寻找其中一项的折半搜索算法(或是在数值分析中类似的勘根算法)。这些算法比一
般的分治算法更能有效地执行。其中,假如算法使用尾部递归的话,便能转换成简单的循环。但在这广义之下,所有使用递归或循环的算法均被视作“分治算法”。因此,有些作者考
虑“分治法”这个名称应只用于每个有最少两个子问题的算法。而只有一个子问题的曾被建议使用减治法这个名称。
分治算法通常以数学归纳法来验证。而它的计算成本则多数以解递回关系式来判定。
历史先例:
折半搜索算法——一个将原来问题连逐地拆细成大约一半大小的单一子问题的分治算法——拥有一段悠长历史。虽然算法在计算机上的清楚描述出现在1946年约翰莫齐利
(John Mauchly)的一篇文章里,然而利用已排序的物件序列去加快搜寻的构想早已在公元前200年的巴比伦尼亚出现。另一个单一子问题的分治算法是找出2个数的最大公因
子的辗转相除法(透过将数字化小至使子问题变得简单),于公元前数世纪已经出现。
一个早期有多个子问题的分治算法是高斯在1805年描述关于快速傅立叶变换的算法,尽管他没有量化地分析它的操作数目,而快速傅立叶变换直至在一世纪之后被重新发现之
前亦没有广泛流传。这个算法现在称为库利-图基快速傅里叶变换算法。
至于专门用于计算机之上而且正确地分析的分治算法早期例子,则可以数到约翰·冯·诺伊曼于1945年发明的归并排序。
另一个显著的例子是Anatolii Alexeevitch Karatsuba于1960年发明在步骤内将两个n位数相乘的算法。它反证了安德雷·柯尔莫哥洛夫于1956年认为这个乘法需
要步骤的猜想。
高德纳举了一个最初并没有涉及计算机的分治算法例子,就是一般邮局用于分发信件的方法:信件在主要邮局根据不同的地理范围而分到不同的袋里,每个袋亦在运送到地区邮
局时分到更小的袋里,如是者直至信件被派发为止。这个方法与早于1929年的打孔卡排序机所用的基数排序相类同。
分治法优点:
解决困难问题
算法效率
同步性
内存存取
修正控制
实现步骤:
循环递归
在每一层递归上都有三个步骤:
- 分解:将原问题分解为若干个规模较小,相对独立,与原问题形式相同的子问题。
- 解决:若子问题规模较小且易于解决时,则直接解。否则,递归地解决各子问题。
- 合并:将各子问题的解合并为原问题的解。
显堆栈
示例代码:
分治法在高级语言中主要的一个思想是递归,LISP语言中的体现出了极丰富的分治法。
以下是归并排序C语言的示例代码,输入参数中,需要排序的数组为array[],起始索引为first,终止索引为last。调用完成后,array[]中从first到last处于升序排列。
1 #include <stdlib.h> 2 #include <stdio.h> 3 4 void merge_sort(int* array , int first , int last); 5 void merge(int* array , int first , int mid , int last); 6 7 int main(void) 8 { 9 int i = 0; 10 int a[4] = {1,32,5,45}; 11 merge_sort(a,0,3); //函数调用 12 for (i = 0; i <= 3; i++) //显示结果 13 { 14 printf("%d\n",a[i]); 15 } 16 system("pause"); 17 return 0; 18 } 19 20 void merge_sort(int* array , int first , int last) 21 //将array【first-last】归并 22 { 23 int mid = 0; 24 if (first<last) 25 { 26 mid = (first+last)/2; //将array【first-last】平分为array【first-mid】和array【mid+1-last】 27 merge_sort(array,first,mid); //递归调用 28 merge_sort(array,mid+1,last); 29 merge(array,first,mid,last); //实现归并 30 } 31 } 32 33 void merge(int* array , int first , int mid , int last) 34 { 35 //将两个有序的子文件array[first..mid)和R[mid+1..last]归并成一个有序的 36 //子文件array[first..last] 37 int* arr = NULL; 38 int i = first , j =mid+1 ,lenth = last-first+1,p = 0 ; 39 arr = (int*)malloc(lenth*sizeof(int)); 40 if (arr == NULL) 41 { 42 printf("分配内存失败!"); 43 exit(0); //申请空间失败 44 } 45 46 while (i <= mid || j <= last) 47 { 48 if (i <= mid && j <= last) 49 { 50 arr[p++] = (array[i]<=array[j])?array[i++]:array[j++]; 51 //上一行代码可改写成如下: 52 /* if (array[i] >= array[j]) 53 { 54 arr[p] = array[j]; 55 p++; 56 j++; 57 } 58 else 59 { 60 arr[p] = array[i]; 61 p++; 62 i++; 63 }*/ 64 } 65 else if (i <= mid && j > last) //若第1个子文件非空,则复制剩余记录到arr中 66 { 67 arr[p++] = array[i++]; 68 } 69 else //若第2个子文件非空,则复制剩余记录到arr中 70 { 71 arr[p++] = array[j++]; 72 } 73 } 74 for (i = 0,j = first ; i < lenth ; i++,j++) //归并完成后将结果复制回array【first-last】 75 { 76 array[j] = arr[i]; 77 } 78 free(arr); 79 80 }
算法复杂度:
比较操作的次数介于和
。 赋值操作的次数是
。 归并算法的空间复杂度为:Θ (n)