博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
全排列(1) ----字典序
阅读量:4141 次
发布时间:2019-05-25

本文共 2227 字,大约阅读时间需要 7 分钟。

引言

对一个给定数据进行全排列,在各种场合经常会用到。组合数学中,生成全排列的方法有很多,卢开澄老师的《组合数学》中就介绍了三种:序数法,字典序法,临位互换法等。其中以字典序法由于算法简单,并且使用的时候可以依照当前状态获取下一个状态,直到所有排列全部完成,方便在程序中随要随用,应用比较广泛,STL中的Next_permutation也是使用此法。

 

算法定义

首先看什么叫字典序,顾名思义就是按照字典的顺序(a-z, 1-9)。以字典序为基础,我们可以得出任意两个数字串的大小。比如 "1" < "12"<"13"。 就是按每个数字位逐个比较的结果。对于一个数字串,“123456789”, 可以知道最小的串是 从小到大的有序串“123456789”,而最大的串是从大到小的有序串“*987654321”。这样对于“123456789”的所有排列,将他们排序,即可以得到按照字典序排序的所有排列的有序集合。

如此,当我们知道当前的排列时,要获取下一个排列时,就可以范围有序集合中的下一个数(恰好比他大的)。比如,当前的排列时“123456879”, 那么恰好比他大的下一个排列就是“123456897”。 当当前的排列时最大的时候,说明所有的排列都找完了。所有排列是一个递增的过程。

于是可以有下面计算下一个排列的算法:

设P是1~n的一个全排列:p=p1p2......pn=p1p2......pj-1pjpj+1......pk-1pkpk+1......pn

  1)从排列的右端开始,找出第一个比右边数字小的数字的序号j(j从左端开始计算),即 j=max{i|pi<pi+1}
  2)在pj的右边的数字中,找出所有比pj大的数中最小的数字pk,即 k=max{i|pi>pj}(右边的数从右至左是递增的,因此k是所有大于pj的数字中序号最大者)
  3)对换pi,pk

  4)再将pj+1......pk-1pkpk+1......pn倒转得到排列p'=p1p2.....pj-1pjpn.....pk+1pkpk-1.....pj+1,这就是排列p的下一个排列。

       例如839647521是数字1~9的一个排列。从它生成下一个排列的步骤如下:
  自右至左找出排列中第一个比右边数字小的数字4  839647521
  在该数字后的数字中找出比4大的数中最小的一个5  839647521
  将5与4交换 839657421
  将7421倒转 839651247
  所以839647521的下一个排列是839651247。

字典序算法实现

#include <iostream>

using namespace std;

void showPermutation( int *p, int n );

bool nextPermutation( int *p, int n );
void swap( int &a, int &b );

int main( int argc, char* argv[] )

{
    int iCount;
    int *pNumber;

    if( argc < 2 ){

        cout << "Usage: " << argv[0] << " number" <<endl;
        cout << "\t" << "number is a plus integer less than 100!" <<endl;
        return 0;
    }

    sscanf( argv[1], "%d", &iCount );

    cout << "Permutation of " << iCount << " numbers :" <<endl;

    pNumber = new int[iCount];

    for( int i=0; i<iCount; i++ ){
        pNumber[i] = i+1;
    }

    showPermutation( pNumber, iCount );

    while( nextPermutation( pNumber, iCount ) ) {

        showPermutation( pNumber, iCount );
        //sleep(1);
    }

    if( pNumber ){

        delete[] pNumber;
    }
    return 0;
}

void showPermutation( int *p, int n )

//打印输出当前序列
{
    for( int i=0; i<n; i++ ){
        cout << p[i] << " ";
    }
    cout << endl;
}

bool nextPermutation( int *p, int n )

//按字典序计算下一个序列
{
    
    int i=0, j=0, k=0;

    i = n-2;

    while( i>=0 && p[i]>p[i+1] ){
        i--;
    }

    if( i<0 ) return false;

    j=n-1;

    while(p[j]<p[i])

        j--;

    swap( p[i], p[j] );

    for( j=n-1, k=i+1; k<j; k++, j-- ){

        swap( p[j], p[k] );
    }

    return true;

}

void swap( int &a, int &b )

{
    int tmp = a;
    a = b;
    b = tmp;
}

转载地址:http://ezrvi.baihongyu.com/

你可能感兴趣的文章
RQP-DEF-0177
查看>>
MySQL字段类型的选择与MySQL的查询效率
查看>>
Java的Properties配置文件用法【续】
查看>>
JAVA操作properties文件的代码实例
查看>>
IPS开发手记【一】
查看>>
Java通用字符处理类
查看>>
文件上传时生成“日期+随机数”式文件名前缀的Java代码
查看>>
Java代码检查工具Checkstyle常见输出结果
查看>>
北京十大情人分手圣地
查看>>
Android自动关机代码
查看>>
Android中启动其他Activity并返回结果
查看>>
2009年33所高校被暂停或被限制招生
查看>>
GlassFish 部署及应用入门
查看>>
iWatch报错: Authorization request cancled
查看>>
iWatch报错: Authorizationsession time out
查看>>
X-code7 beta error: warning: Is a directory
查看>>
Error: An App ID with identifier "*****" is not avaliable. Please enter a different string.
查看>>
Xcode 报错: Extra argument in call
查看>>
iTunes Connect 上传APP报错: Communication error. please use diagnostic mode to check connectivity.
查看>>
#import <Cocoa/Cocoa.h> 报错 Lexical or Preprocessor Issue 'Cocoa/Cocoa.h' file not found
查看>>