全排列函数、组合函数

1

1、求一个全排列函数：如p([1,2,3])输出：

         [123],[132],[213],[231],[321],[312].

2、求一个组合函数如p([1,2,3])输出：

         [1],[2],[3],[1,2],[2,3],[1,3],[1,2,3]

这两问可以用伪代码。

void swap(int *a, int *b) //交换函数
{int tmp;tmp =*a;*a=*b;*b=tmp;
}

以为很简单，手写全排序第一遍就错了。

func1( a[], n, m)if(m>n-1) return; endprint(a);for(i=m+1 to n-1)swap(a[m], a[i]);func1(a, n, m+1);swap(a[m], a[i]);end
end

调用：func1(a, n, 0)

发现以a[0]开头的只有一组，其他可以满足，因为i=m+1。func1的理解就是输出a，并将a[m]分别与a[m+1]..a[n]交换后的后n-m-1个元素的排列组合。当调用了func1(a, n, 0)时，以a[0]开头的只有一组。改为for(i=m to n)可以吗？

No！这时候，却多出很多。由于m>n-1才是终止条件。即要进入n-1层，每层有n-m+1个分支，于而每进入一层就会print(a)，这样就会输出的次数不等于n!了（一般情况下）

纠结了n久之后，再改，应该当m=n-1时才输出。Bingo！

func1( a[], n, m)
　　if(m==n-1) 
　　    print(a);
       return;
    end
 　 for(i=m+1 to n-1)
　　　　swap(a[m], a[i]);
　　　　func1(a, n, m+1);
　　　　swap(a[m], a[i]);
　　end
end

调用：func1(a, n, 0)

可是结果却和题目的不完全一致。因为是输出的位置的先后问题。换成尾递归即可。

void func1(int a[], int m, int n)
{        
    for (int j=m; j<n; j++)
    {
        swap(&a[m], &a[j]);
        func1(a, m+1, n);
        swap(&a[m], &a[j]);
    }
    if (m==n-1)
    {
        for (int i=0; i<n; i++)
        {
            cout<<a[i]<<" ";
        }
        cout<<endl;
        return;
    }
}

简单的题目也不见得简单啊~（其实个人水平问题。。。）

2、只考虑有最大有32个数的情况：

void func2(int a[], int n)
{
    int cnt=1;

    if (n>32 || n<=0)
    {
        return;
    }
    cnt = cnt<<n;

    unsigned int tmp;
    for (unsigned int i=0; i<=cnt; i++)
    {
        tmp = 1;
        for (int j=0; j<n; j++)
        {
            if (i & tmp)
            {
                cout<<a[j]<< " ";
            }
            tmp = tmp<<1;
        }
        cout<<endl;
    }
}

总结：

常用的排列组合算法，包括全排列算法，全组合算法，m个数选n个组合算法等。

2. 排列算法

常见的排列算法有：

(A)字典序法

(B)递增进位制数法

(C)递减进位制数法

(D)邻位对换法

(E)递归法

介绍常用的两种：

(1) 字典序法

对给定的字符集中的字符规定了一个先后关系，在此基础上按照顺序依次产生每个排列。

[例]字符集{1,2,3},较小的数字较先,这样按字典序生成的全排列是:123,132,213,231,312,321。

生成给定全排列的下一个排列 所谓一个的下一个就是这一个与下一个之间没有字典顺序中相邻的字符串。这就要求这一个与下一个有尽可能长的共同前缀，也即变化限制在尽可能短的后缀上。

算法思想：

设P是[1,n]的一个全排列。

P=P1P2…Pn=P1P2…Pj-1PjPj+1…Pk-1PkPk+1…Pn , j=max{i|Pi<Pi+1}, k=max{i|Pi>Pj} ,对换Pj，Pk，将Pj+1…Pk-1PjPk+1…Pn翻转， P’= P1P2…Pj-1PkPn…Pk+1PjPk-1…Pj+1即P的下一个

例子:839647521的下一个排列.

从最右开始,找到第一个比右边小的数字4(因为4<7，而7>5>2>1),再从最右开始,找到4右边比4大的数字5(因为4>2>1而4<5),交换4、5,此时5右边为7421,倒置为1247,即得下一个排列:839651247.用此方法写出全排列的非递归算法如下

该方法支持数据重复，且在C++ STL中被采用。

(2) 递归法

设一组数p = {r1, r2, r3, … ,rn}, 全排列为perm(p)，pn = p – {rn}。则perm(p) = r1perm(p1), r2perm(p2), r3perm(p3), … , rnperm(pn)。当n = 1时perm(p} = r1。

如：求{1, 2, 3, 4, 5}的全排列

1、首先看最后两个数4, 5。 它们的全排列为4 5和5 4, 即以4开头的5的全排列和以5开头的4的全排列。

由于一个数的全排列就是其本身，从而得到以上结果。

2、再看后三个数3, 4, 5。它们的全排列为3 4 5、3 5 4、 4 3 5、 4 5 3、 5 3 4、 5 4 3 六组数。

即以3开头的和4,5的全排列的组合、以4开头的和3,5的全排列的组合和以5开头的和3,4的全排列的组合.

以上的算法1的即是使用这个方法。

3组合算法

3.1 全组合

在此介绍二进制转化法，即，将每个组合与一个二进制数对应起来，枚举二进制的同时，枚举每个组合。如字符串：abcde

00000 <– –> null

00001<– –> e

00010 <– –> d

… …

11111 <– –> abcde

以上的算法2的即是使用这个方法。

3.2 从n中选m个数

(1) 递归

a. 首先从n个数中选取编号最大的数，然后在剩下的n-1个数里面选取m-1个数，直到从n-(m-1)个数中选取1个数为止。

b. 从n个数中选取编号次小的一个数，继续执行1步，直到当前可选编号最大的数为m。

下面是递归方法的实现：

/// 求从数组a[1..n]中任选m个元素的所有组合。

/// a[1..n]表示候选集，n为候选集大小，n>=m>0。

/// b[1..M]用来存储当前组合中的元素(这里存储的是元素下标)，

/// 常量M表示满足条件的一个组合中元素的个数，M=m，这两个参数仅用来输出结果。

void combine( int a[], int n, int m,  int b[], const int M )

{

  for(int i=n; i>=m; i--)   // 注意这里的循环范围

  {

    b[m-1] = i - 1;

    if (m > 1)

      combine(a,i-1,m-1,b,M);

    else                     // m == 1, 输出一个组合

    {

      for(int j=M-1; j>=0; j--)

      cout << a[b[j]] << " ";

      cout << endl;

    }

  }

}

void  func19(int m, int n) //非递归
{
    int *arr;
    int i, pos;

    if (m<n)
    {
        return ;
    }

    arr = new int[n];

    for (i=0; i<n; ++i)
    {
        arr[i] = i+1;
    }

    for (i=0; i<n; ++i)
    {
        cout<<arr[i]<<" ";
    }
    cout<<endl;

    pos = n-1;
    while(1)
    {
        if (arr[n-1] == m )
        {
            pos--;
        }
        else
        {
            pos = n-1;
        }

        arr[pos] ++;

        for (i=pos+1; i<n; i++)
        {
            arr[i] = arr[i-1]+1;
        }

        for (i=0; i<n; ++i)
        {
            cout<<arr[i]<<" ";
        }
        cout<<endl;

        if (arr[0] == m-n+1)
        {
            break;
        }
    }
}

 

(2) 01转换法

本程序的思路是开一个数组，其下标表示1到n个数，数组元素的值为1表示其代表的数被选中，为0则没选中。

首先初始化，将数组前n个元素置1，表示第一个组合为前n个数。

然后从左到右扫描数组元素值的“10”组合，找到第一个“10”组合后将其变为“01”组合，同时将其左边的所有“1”全部移动到数组的最左端。

当第一个“1”移动到数组的n-m的位置，即n个“1”全部移动到最右端时，就得到了最后一个组合。

例如求5中选3的组合：

1 1 1 0 0 //1,2,31 1 0 1 0 //1,2,41 0 1 1 0 //1,3,40 1 1 1 0 //2,3,41 1 0 0 1 //1,2,51 0 1 0 1 //1,3,50 1 1 0 1 //2,3,51 0 0 1 1 //1,4,50 1 0 1 1 //2,4,50 0 1 1 1 //3,4,5

4. 参考资料

(1) http://www.cnblogs.com/nokiaguy/archive/2008/05/11/1191914.html

(2) http://comic.sjtu.edu.cn/thucs/GD_jsj_025y/text/chapter01/sec05/default.htm

(3) http://blog.csdn.net/sharpdew/archive/2006/05/25/755074.aspx

(4) http://xiaomage.blog.51cto.com/293990/74094

以上总结源于：董的博客