poj 1147 Binary codes BWT压缩算法

来源:互联网 发布:js for in循环 编辑:程序博客网 时间:2024/04/20 19:19

题意:一个长度为N的01序列,会有N个不同的轮换(当然,字符相同,其中也可能会有相同的),将这N个不同轮换按字典序排

序,取排序后的每个轮换的最后一排,组成一个序列。题中给出压缩后的序列,求原始序列,输出的是字典序最小的那个序列。

思路:这题基于一个性质:在已经排序好的矩阵中,对于首位相同的两行,经过左移一位的操作后,形成的新的两行的先后次序不发

生改变。即:设i行在j行前面,i行左移一位变成p行,j左移一位后变成q行,p还是在q的前面。已知最后一列,那么我们可以知道一行

有几个零(cnt0个)几个一(cnt1个),那么我们自然能得到第一列,前cnt0个为零,之后cnt1个为一。由上述性质,第一列第i个零所在的

p行一定是转移到最后一列第i个零所在的q行,且q行的首位是是p行首位的后一位。那么我们利用性质即可得到行转移的数组next。

最后沿着next进行行转移,依次输出行的首位即为第一行,详见代码:

/*********************************************************  file name: poj1147.cpp  author : kereo  create time:  2015年02月21日 星期六 18时08分02秒*********************************************************/#include<iostream>#include<cstdio>#include<cstring>#include<queue>#include<set>#include<map>#include<vector>#include<stack>#include<cmath>#include<string>#include<algorithm>using namespace std;typedef long long ll;const int sigma_size=26;const int N=100+50;const int MAXN=100000+50;const int inf=0x3fffffff;const double eps=1e-8;const int mod=1000000000+7;#define L(x) (x<<1)#define R(x) (x<<1|1)#define PII pair<int, int>#define mk(x,y) make_pair((x),(y))int n,cnt0,cnt1;int a[MAXN],num0[MAXN],num1[MAXN],next[MAXN];int main(){    while(~scanf("%d",&n)){        cnt0=cnt1=0;        for(int i=1;i<=n;i++){            scanf("%d",&a[i]);            if(a[i])                cnt1++,num1[cnt1]=i;            else                 cnt0++,num0[cnt0]=i;        }        for(int i=1;i<=n;i++){            if(i<=cnt0)                next[i]=num0[i];            else                 next[i]=num1[i-cnt0];        }        int k=1,ans;        for(int i=1;i<=n;i++){            if(i == 1)                printf("%d",ans=k<=cnt0 ? 0 : 1);            else                printf(" %d",ans=k<=cnt0 ? 0 : 1);            k=next[k];        }        printf("\n");    }return 0;}


0 0
原创粉丝点击