排列组合算法小结(未完)

做GCJ 2011 Qualification Round中的C题遇到的，以前居然没有接触过，囧。。。

组合算法有三个，这里给出了递归和回溯的方法，我照搬了一下，不过改成了C实现。。。

回溯法貌似挺慢的，做题时有个大数据的test，结果跑了6，7分钟都没跑完，这里又给出了一个高效的排列组合算法，我把组合实现了，排列我觉得不是很高效，暂时没写。。。

这里有全排列的递归算法，我给搬过来了，做swap时本想用异或实现的，可是注意到i，j可能指向同一个值，结果异或之后变成了0，无奈改成普通的交换方式

 

---

#include <stdio.h> //组合递归算法 void combine1(int a[], int i, int b[], int j, int n, int m, int M) { // b size == m if(m <= 0) { for(j = 0; j < M; j++) printf("%d ", a[b[j]]); printf("/n"); return; } for(; i <= n - m; i++) { b[j] = i; combine1(a, i+1, b, j+1, n, m-1, M); } } //组合回溯算法 void combine2(int a[], int b[], int n, int m) { //b size == m + 1 int i, k, flag; for(i = 0; i <= m; i++) b[i] = i - 1; k = m; flag = 1; while(b[0] == -1) { if(flag) { for(i = 1; i <= m; i++) { printf("%d ", a[b[i]]); } printf("/n"); flag = 0; } b[k]++; if(b[k] == n) { b[k--] = 0; continue; } if(k < m) { b[++k] = b[k-1]; continue; } if(k == m) flag = 1; } } //高效组合算法 void combine3(int a[], int b[], int n, int m) { int i, j, k; int q = 0; for(i = 0; i < n; i++) { //b[] = 1 1 1 0 0 if(i < m) b[i] = 1; else b[i] = 0; } while(!q) { for(i = 0; i < n; i++) { if(b[i]) printf("%d ", a[i]); } printf("/n"); q = 1; for(i = 0; i < n-1; i++) { if(b[i] == 1 && b[i+1] == 0) { // <1, 0> -> <0, 1> q = 0; b[i] = 0; b[i+1] = 1; for(j = 0, k = i-1; j < k; j++, k--) { if(b[k]) { b[k] = b[j]; b[j] = 1; } } break; } } } } //全排列递归算法 void swap(int a[], int i, int j) { //i == j so a[i] = a[j] = 0, that's really bad... //a[i] ^= a[j]; //a[j] ^= a[i]; //a[i] ^= a[j]; //a[i] = a[i] ^ a[j]; //a[j] = a[i] ^ a[j]; //a[i] = a[i] ^ a[j]; int tmp; tmp = a[i]; a[i] = a[j]; a[j] = tmp; } void permutation1(int a[], int s, int n) { int i; if(s >= n) { for(i = 0; i < n; i++) printf("%d ", a[i]); printf("/n"); return; } for(i = s; i < n; i++) { swap(a, i, s); permutation1(a, s+1, n); swap(a, i, s); } } //高效全排列算法 void permutation2(int a[], int n) { printf("permutation for efficiency, to be continued.../n"); } int main() { int a1[5] = {1, 2, 3, 4, 5}; int a2[7] = {1, 2, 3, 4, 5, 6, 7}; int a3[3] = {1, 2, 3}; int b1[6], b2[8], b3[5]; combine1(a1, 0, b1, 0, 5, 3, 3); //combine1(a2, 0, b2, 0, 7, 4, 4); printf("/n/n"); combine2(a1, b1, 5, 3); printf("/n/n"); combine3(a1, b3, 5, 3); printf("/n/n"); permutation1(a3, 0, 3); printf("/n/n"); permutation2(a3, 3); printf("/n/n"); }  

---

还不太完整，有机会补全。。。