分治法的思想:将原问题分解为几个规模较小但类似于原问题的子问题,递归的求解这些子问题,然后再合并这些子问题的解来建立原问题的解。
分治法在每层递归是遵循的三个步骤:
(1)分解原问题为若干个子问题,这些子问题是原问题的规模较小的实例。
(2)解决这些子问题,队规的求解各个子问题,当子问题规模足够小的时候,直接求解。
(3)合并这些子问题的解构成原问题的解。
显然归并排序是一个非常经典规矩的分治法的例子,鉴于之前已经写过一篇关于归并排序的博文,这里不在使用归并排序作为例子。
注意分治法的每一层递归中的第一步分解,可能产生两个子问题(如归并排序、二分查找等),也可能产生多个子问题(如排列、组合等),产生两个子问题的时候当然比较容易理解,而产生多个子问题的时候需要使用环循罗列这些子问题。
下面就以排列和组合算法为例,介绍产生多个子问题的分治算法。
一、排列
问题:输入一个字符串,打印出该字符串中字符的所有排列。
分析:利用分治法的思想,
(1)先将原问题分解,假如输入的字符串长度是n,那么第一次选择可能是第一个字符、也可能是第二个、。。。也可能是第n个,但是不管是哪一个,只要选出第一个字符,就可以在剩下的n-1个字符里面继续选择一个了,所以需要将原问题分解为n个子问题(每个子问题为第一步选择的是i,然后再对除了i之外的字符进行全排列),到现在可以发现如果直接按照顺序分解之后,对除了i之外的字符进行全排列,不是那么容易实现递归,于是想到将每个元素(包括第一个元素)都与第一个元素交换,然后分解成的子问题就是先将每个元素与第一个元素交换并选出,然后对第二个到最后的所有元素全排列。注意每次个子问题考虑完之后需要将交换的元素换回。
(2)利用递归解决每个子问题
(3)当所有问题都解决的时候,子问题的解组合起来就是原问题的解了
如:输入字符串为abc ,排列函数为permutation()那么分解成的子问题为a+permutation(bc)、b+permutation(ac)、c+permutation(ab)
1 #include "stdafx.h" 2 #include<iostream> 3 using namespace std; 4 void print(char *str) 5 { 6 char *p=str; 7 while(*p) 8 { 9 cout<<*p<<' '; 10 p++; 11 } 12 } 13 void bianli(char *str,int begin,int length) 14 { 15 char temp; 16 int i; 17 if(begin==length-1) 18 { 19 print(str); 20 cout<<endl; 21 return ; 22 } 23 //可以选取某一个值(包括begin自己)与begin的位置交换,然后对剩下的字符全排列 24 //所以对于每一个位置要么选择先交换,然后递归,要么选择不交换(即交换两次) 25 for(i=begin;i<length;i++) 26 { 27 temp=str[begin]; 28 str[begin]=str[i]; 29 str[i]=temp; 30 31 bianli(str,begin+1,length); 32 33 temp=str[begin]; 34 str[begin]=str[i]; 35 str[i]=temp; 36 } 37 } 38 39 int _tmain(int argc, _TCHAR* argv[]) 40 { 41 char str[4]="123"; 42 bianli(str,0,3); 43 return 0; 44 }
二、组合
问题:找出从自然数1、2、3。。。n中任取r个元素的所有组合
分析:
1、分解:与排列不同,组合里每个元素在一种只出现一次,所以并不需要交换元素,而是每次从n个数中按照某种顺序取一个元素,然后考虑全面了即可,如每次取一个最大值,那么只要元素个数>k则是子问题的一种,剩下的思想和排列差不多。
1 #include<iostream> 2 using namespace std; 3 int a[100];//用于存放组合的结果 4 void zuhe(int n,int k) 5 { 6 for(int i=n;i>=k;i--)//顺序选取组合中最大的数 7 { 8 a[k]=i; 9 if(k>1) 10 { 11 zuhe(i-1,k-1); 12 } 13 else 14 { 15 for(int i=1;i<=a[0];i++) 16 { 17 cout<<a[i]<<" "; 18 } 19 cout<<endl; 20 } 21 } 22 } 23 int main() 24 { 25 int n,k; 26 cin>>n>>k; 27 a[0]=k; 28 zuhe(n,k); 29 return 0; 30 }