全排列算法整理 本文整理了全排列算法的:
- 递归实现
- 含重复元素的递归实现
- 深搜实现(字典序)
- 含重复元素的深搜实现(字典序)
- 非递归字典序实现——下一个排列(可含重复元素)
设R={r1,r2,...,rn}R = { r 1 , r 2 , . . . , r n }是要进行排列的nn个元素, Ri=R?{ri}R i = R ? { r i } 。集合XX中元素的全排列记为 Perm(X)P e r m ( X ) 。 (ri)Perm(X)( r i ) P e r m ( X )表示在全排列Perm(X)P e r m ( X )的每一个排列钱加上前缀rir i得到的排列。 RR的全排列可归纳定义如下:依次递归定义,可设计产生Perm(R)P e r m ( R )的递归算法代码如下:
当n=1n = 1时, Perm(R)=(r)P e r m ( R ) = ( r ),其中rr是集合RR中唯一的元素;
当n>1n > 1时, Perm(R)P e r m ( R )由(r1)Perm(R1),(r2)Perm(R2),...,(rn)Perm(Rn)( r 1 ) P e r m ( R 1 ) , ( r 2 ) P e r m ( R 2 ) , . . . , ( r n ) P e r m ( R n )构成。
#include
#include
using namespace std;
template inline void Swap(Type &a, Type &b) {
Type temp = a;
a = b;
b = temp;
}
template
void Perm(Type list[], int k, int n) {
if(k == n) {
for(int i = 0;
i <= n;
i++) cout << list[i];
cout << endl;
}
else {
for(int i = k;
i <= n;
i++) {
Swap(list[k], list[i]);
Perm(list, k+1, n);
Swap(list[k], list[i]);
}
}
}int main() {
char s[] = "abcd";
int len = strlen(s);
Perm(s, 0, len-1);
return 0;
}
算法 Perm(list,k,m)P e r m ( l i s t , k , m )递归地产生所有前缀是 list[0:k?1]l i s t [ 0 : k ? 1 ],且后缀是 list[k:m]l i s t [ k : m ]的全排列的所有排列。函数调用 Perm(list,0,n?1)P e r m ( l i s t , 0 , n ? 1 )则产生 list[0:n?1]l i s t [ 0 : n ? 1 ]的全排列。
在一般情况下, k
数字122,第一个数字1和第二个数字2交换,可以得到2 12,这将生成以2为前缀,以12的全排列为后缀的所有排列。
第一个数字1与第三个数字2交换,可以得到2 21, 这将生成以2为前缀,以21的全排列为后缀的所有排列。
而12和21的全排列是相同的,所以上述两种交换得到的排列一定是重复的。
去重的全排列就是从第一个数字起每个数分别与它后面非重复出现的数字交换。
即第i个数与第j个数交换时,要求[i,j)中没有与第j个数相等的数。
#include
#include
using namespace std;
template inline void Swap(Type &a, Type &b)
{
Type temp = a;
a = b;
b = temp;
}
/*
* 去重的全排列就是从第一个数字起每个数分别与它后面非重复出现的数字交换。
* 第i个数与第j个数交换时,要求[i,j)中没有与第j个数相等的数。
*/
bool isSwap(char *list, int begin, int end)
{
for(int i = begin;
i < end;
i++) {
if(list[i] == list[end])
return false;
}
return true;
}
void Perm(char *list, int k, int n)
{
if(k == n) {
for(int i = 0;
i <= n;
i++) cout << list[i];
cout << endl;
}
else {
for(int i = k;
i <= n;
i++) {
if(isSwap(list, k, i)) {
Swap(list[k], list[i]);
Perm(list, k+1, n);
Swap(list[k], list[i]);
}
}
}
}int main()
{
char s[] = "1232";
int len = strlen(s);
Perm(s, 0, len-1);
return 0;
}
深搜实现 深搜也可以实现全排列,用一个visit[]数组记录当前那些元素已经被选,继续搜索没有被选的元素。
1. 数的全排列 给定数字n,输出数字1~n的全排列。比如给定数字3,输出:1 2 3 ,1 3 2 ,2 1 3 ,2 3 1 ,3 1 2 ,3 2 1
代码:
#include
#include
#define MAXN 256
using namespace std;
int n;
//for dfs_num()
int perm[MAXN];
//全排列
int visit[MAXN];
void dfs_num(int step);
//数的全排列,step为当前选择元素数量int main()
{
n = 4;
memset(visit, 0, sizeof(visit));
memset(perm, 0, sizeof(perm));
dfs_num(1);
return 0;
}
/*
* 数的全排列
*/
void dfs_num(int step)
{
if(step == n+1) {//step到n+1,即已经选满n个数,则深搜结束,得到全排列
for(int i = 1;
i <= n;
i++) {
cout << s[perm[i]-1] << " ";
}
cout << endl;
}
else {
//从1~n中选一个未被选中的元素,然后继续深搜。深搜结束后回溯。
for(int i = 1;
i <= n;
i++) {
if(!visit[i]) {//如果i未被标记
perm[step] = i;
//记录
visit[i] = 1;
//标记
dfs_num(step+1);
//深搜下一步
visit[i] = 0;
//回溯
}
}
}
}
2.字符的全排列 上面的数的全排列稍微改动一下就是字符的全排列。我们可以生成数字0~n-1的全排列,这个数的全排列即字符数组的下标的全排列。输出的时候按下标输出字符即可。visit[]记录的还是下标。
代码:
#include
#include
#define MAXN 256
using namespace std;
string s;
// for dfs_ch()
int len;
//s.length
char perm_ch[MAXN];
int visit[MAXN];
void dfs_ch(int step);
//通过记录下标实现的不含重复元素的全排列
void dfs_ch1(int step);
//通过记录元素实现的不含重复元素的字符的全排列void dfs_ch2(int step);
//通过记录下标实现的含有重复元素的全排列
void dfs_ch3(int step);
//通过记录元素出现个数实现的含重复元素的全排列int main()
{
s = "abcd";
len = s.length();
memset(visit, 0, sizeof(visit));
memset(perm, 0, sizeof(perm));
dfs_ch(0);
return 0;
}
/*
* 通过记录下标的不含重复元素的字符的全排列
*/
void dfs_ch(int step)
{
if(step == len) {
for(int i = 0;
i < len;
i++) {
cout << s[perm[i]] << " ";
}
cout << endl;
}
else {
for(int i = 0;
i < len;
i++) {
if(!visit[i]) {
perm[step] = i;
visit[i] = 1;
dfs_ch(step+1);
visit[i] = 0;
}
}
}
}
也可以用一个256个元素的visit[]数组记录当前元素是否被选中过,而不是记录当前的下标。
/*
* 通过记录元素的不含重复元素的字符的全排列
*/
void dfs_ch1(int step) {
if(step == len) {
for(int i = 0;
i < len;
i++) {
cout << perm_ch[i] << " ";
}
cout << endl;
}
else {
for(int i = 0;
i < len;
i++) {
if(!visit[(int)s[i]]) {
perm_ch[step] = s[i];
//perm_ch[]是一个char数组,用于记录排列
visit[(int)s[i]] = 1;
dfs_ch1(step+1);
visit[(int)s[i]] = 0;
}
}
}
}
含重复元素的深搜实现 上面第二种字符的全排列方式要求字符串s[]中不允许出现重复元素,否则该算法将无法满足终止条件。
如果用上面第一种算法,按照下标的全排列来生成字符串122的全排列,全排列的解答树如下:
文章图片
我们注意到橘黄色的树显示的都是已经出现过的排列,即我们需要剪掉的树枝。
我们希望在搜索子树的时候,不要出现重复的子树,否则一定会出现重复的排列。
如上图的第一层的第三棵子树,其根节点为2,与第二棵子树的根节点相同,又由于第二棵子树和第三棵拥有同样的祖先,所以在这两个子树下搜索时visit[]数组是完全相同的,所以这两个子树也必定完全相同。
所以,在第一层的搜索中,当发现第三个节点与前面已经出现过的一个节点相同时(即第二个节点),就不该继续搜索此树,应将其剪掉。
如果字符串s[]已经有序,则每次搜索新的子树的时候,只需要判断当前搜索到的子树的根节点是否与上一个搜索到的子树的根节点相同。用一个变量pre记录上一个搜索过的节点的下标,如果当前搜索的节点s[i]==s[pre],则跳过该次搜索,即剪掉这棵子树。
代码:
/*
* 通过记录下标实现的含重复元素的全排列
*/
void dfs_ch2(int step)
{
if(step == len) {
for(int i = 0;
i < len;
i++) {
cout << s[perm[i]] << " ";
}
cout << endl;
}
else {
int pre = 0;
//记录前一个选中的根节点
for(int i = 0;
i < len;
i++) {
if(!i || s[pre] != s[i]) { //与前一个选中的子树的根节点不重复
if(!visit[i]) {
perm[step] = i;
visit[i] = 1;
dfs_ch2(step+1);
visit[i] = 0;
pre = i;
//记录pre
}
}
}
}
}
如果用visit[]记录元素的方式, 我们可以用一个count[]数组记录每个元素出现的次数,表明该元素可以在全排列中选中的次数,将原本的 if(!visit[(int)s[i]])可以改为if(visit[(int)s[i]] < count[(int)s[i]])。
如果字符串s[]已经有序,则每次只需要搜索s[]的第一个元素和“与前一个元素不相同”的元素。即每次搜索子树的时候,在 for(int i = 0; i < len; i++) 和起后的大括号之前加上if(!i || s[i] != s[i-1])。这样就保证了重复元素绝对不会出现在同一层(s[]已经有序的情况下)。在 for(int i = 0; i < len; i++) 和起后的大括号之前加上if(!i || s[i] != s[i-1]),问题就解决了。
代码:
/*
* 通过记录元素出现个数实现的含重复元素的全排列
*/
void dfs_ch3(int step) {
if(step == len) {
for(int i = 0;
i < len;
i++) {
cout << perm_ch[i] << " ";
}
cout << endl;
}
else {
for(int i = 0;
i < len;
i++) {
if(!i || s[i] != s[i-1]) {
if(visit[(int)s[i]] < count[(int)s[i]]) {
perm_ch[step] = s[i];
visit[(int)s[i]]++;
dfs_ch3(step+1);
visit[(int)s[i]]--;
}
}
}
}
}
非递归字典序实现 非递归的全排列即不断地求当前字符串在字典序下的“下一个排列”。如”1234”的下一个排列就是”1243”,“1243”的下一个排列是“1324” …… 直到求得“4321”,找不到比当前排列字典序更大的排列,算法终止。
如何计算字符串的下一个排列了?来考虑”926520”这个字符串,我们从后向前找第一双相邻的递增数字,”20”、”52”都是非递增的,”26 “即满足要求,称前一个数字2为替换数,替换数的下标称为替换点,再从后面找一个比替换数大的最小数(这个数必然存在),0、2都不行,5可以,将5和2交换得到”956220”,然后再将替换点后的字符串”6220”颠倒即得到”950226”。
对于像”4321”这种已经是最“大”的排列,采用STL中的处理方法,将字符串整个颠倒得到最“小”的排列”1234”并返回false。
用更严谨一点的描述:
设P是1~n的一个全排列: p=p1p2......pn=p1p2......pj?1pjpj+1......pk?1pkpk+1......pnp = p 1 p 2 . . . . . . p n = p 1 p 2 . . . . . . p j ? 1 p j p j + 1 . . . . . . p k ? 1 p k p k + 1 . . . . . . p n关于这个next_permutation算法的证明,参考:http://blog.csdn.net/lohiaufung/article/details/52901564
1. 从排列的右端开始,找出第一个比右边数字小的数字的序号j(j从左端开始计算), j=max{i|pi2. 在 pjp j 的右边的数字中,找出所有比pj大的数中最小的数字pk,即k=maxi|pi>pjk = m a x i | p i > p j (右边的数从右至左是递增的,因此k是所有大于 pjp j 的数字中序号最大者)
3. 对换pi,pk
4. 再将 pj+1......pk?1pkpk+1......pnp j + 1 . . . . . . p k ? 1 p k p k + 1 . . . . . . p n 倒转得到排列 p′=p1p2.....pj?1pjpn.....pk+1pkpk?1.....pj+1p ′ = p 1 p 2 . . . . . p j ? 1 p j p n . . . . . p k + 1 p k p k ? 1 . . . . . p j + 1 ,这就是排列p的下一个排列。
还有:http://blog.fly2best.com/algo/2013/06/21/next_permutation/
代码:
#include
#include
using namespace std;
template inline void Swap(Type &a, Type &b){
Type temp = a;
a = b;
b = temp;
}void Reverse(char *i, char *j) {
while(i < j) {
Swap(*i, *j);
i++;
j--;
}
}bool Next_permutation(char *s) {
int len = strlen(s);
if(!len) return false;
char *pEnd, *p, *q, *pFind;
pEnd = s+len-1;
p = pEnd;
while(p != s) {
q = p;
p--;
if(*p < *q) {
pFind = pEnd;
while(*pFind <= *p) pFind--;
Swap(*p, *pFind);
Reverse(q, pEnd);
return true;
}
}
Reverse(s, pEnd);
return false;
}
int main()
{
char s[] = "1234";
do{
cout << s << endl;
}while(Next_permutation(s));
return 0;
}
【分治|全排列算法整理】本文参考了:
1. http://blog.csdn.net/morewindows/article/details/7370155/
2. http://blog.csdn.net/cpfeed/article/details/7376132
3. http://blog.csdn.net/lohiaufung/article/details/52901564
4. http://blog.fly2best.com/algo/2013/06/21/next_permutation/
推荐阅读
- 数论|hdu 5322 Hope(分治+NTT)
- codeforces|Codeforces Global Round 10 D. Omkar and Bed Wars(思维,分块)
- hdu|【hdu 5354】Bipartite Graph【分治 并查集】
- codeforce514 D. Nature Reserve 凸函数三分
- Pow(x, n) (计算x的n次方)
- 模板|LOJ #2340. 「WC2018」州区划分(FMT子集卷积)
- 点分治|点分治时间复杂度
- LeetCode 973 TopK问题 分治算法