组合枚举方法与常见优化

ACM中常见的有两种问题:

 

1.按字典顺序输出全排列

 

如: abc 可产生

 

abc

acb

bac

bca

cab

cba

 

2.m字母选n个字母，按字典顺序输出

 

如: abc 中3选2

ab

ac

bc

ba

ca

cb

 

测试1:

 

首先说全排列的产生方式,中学的思想就是n!(n的阶乘),那么我们就可以用递归的方式来实现.

 

1.排序

2.枚举

   模拟dance_link的思想，使其达到N!

   枚举过程中，要解决字串中有重复字符问题

 

#include<stdio.h> #include<string.h> #include<stdlib.h> #define MAX_LEN 16 int link[MAX_LEN]; char path[MAX_LEN]; int cmp(void const *a,void const *b) { return (*(char *)a)-*((char *)b); } void dfs(int inc,int len,char *str) { int p,back; if(inc == len) { printf("%s/n",path); return; } //有点像仿dance_link的做法,使程序的运行次数,刚好是N! for(p = 0;link[p] > 0;) { back = link[p]; //保存节点 path[inc] = str[back-1]; //保存路径 link[p] = link[back]; //删除节点 dfs(inc+1,len,str); link[p] = back; //恢复节点 //解决字串中有多个相同字符的全排列问题,如字串aabb的全排列 //也算是剪枝的作用 p = link[p]; while(link[p]>0 && str[back-1] == str[link[p]-1]) p = link[p]; } } int main() { int len,i; char temp[MAX_LEN]; while(scanf("%s",temp)!=EOF) { len = strlen(temp); //模仿dance_link,初始化链表 for(i = 0;i<len; i++) link[i] = i+1; link[len] = 0; //先排序 qsort(temp,len,sizeof(char),cmp); memset(path,0,sizeof(char)*(MAX_LEN)); dfs(0,len,temp); } return 0; }

 

 

 

测试2:

 

m字母选n个字母，按字典顺序输出,问题就容易解决了，程序跟上边一样，

 

只是 dfs(inc,len,str) 中的len = n就可以了;

 

 

与网上一段流行程序作个比较:

 

#include <stdio.h> #include<string.h> #include<ctime> int n = 0; void swap(char *a, char *b) { int m; m = *a; *a = *b; *b = m; } void perm(char list[], int k, int m) { int i; if(k > m) { n++; } else { for(i = k; i <= m; i++) { swap(&list[k], &list[i]); perm(list, k + 1, m); swap(&list[k], &list[i]); } } } int main() { clock_t t1,t2; char list[18]; while(scanf("%s",list)!=EOF) { n = 0; t1 = clock(); perm(list, 0, strlen(list)-1); t2 = clock(); printf("%lf total:%d/n", double(t2-t1)/CLOCKS_PER_SEC,n); } return 0; }

 

 

测试3:

 

STL的写法:

 

#include<ctime> #include <iostream> #include <algorithm> #include <iterator> int main() { using namespace std; int n =0; char a[] = "abcdefghijk"; clock_t t1,t2; t1 = clock(); do { n++; //copy( a,a+11,ostream_iterator<char>(cout," ") ); //cout << endl; }while( next_permutation(a,a+11) ); t2 = clock(); printf("%lf total:%d/n", double(t2-t1)/CLOCKS_PER_SEC,n); return 0; }

 

 

测试环境: CPU AMD X2 4200+ 2.2G

 

测试字串:abcdefghijk

 

测试1 (我优化的程序)测试时间: 4.531s

 

测试2 流行程序测试时间: 7.859s

 

测试3 STL 程序测试时间: 9.375s

 

4.531/7.859 = 0.57....

4.531/9.375 = 0.48....

 

提高将近40%~50多的速度