注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

ydc的博客

 
 
 

日志

 
 

bzoj1031  

2013-02-05 13:17:15|  分类: bzoj |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

字符加密

如此裸的后缀数组!

字符串倍长之后构造后缀数组,然后输出

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 200010
using namespace std;
char a[MAXN];
int len,Sa[MAXN],rank[MAXN],height[MAXN];
int x[MAXN],y[MAXN],sum[MAXN];
void read()
{
    scanf("%s",a);
    len=strlen(a);
    for(int i=len;i<len+len;i++)
        a[i]=a[i-len];
    len<<=1;
}
void sort(int a[],int b[],int c[],int n,int m)
{
    for(int i=0;i<=m;i++)  sum[i]=0;
    for(int i=0;i<n;i++)  sum[c[a[i]]]++;
    for(int i=1;i<=m;i++)  sum[i]+=sum[i-1];
    for(int i=n-1;i>=0;i--)  b[--sum[c[a[i]]]]=a[i];
}
void make_Sa()
{
    for(int i=0;i<len;i++)
        x[i]=a[i],rank[i]=i;
    sort(rank,Sa,x,len,256);
    rank[Sa[0]]=1;
    for(int i=1;i<len;i++)
        rank[Sa[i]]=rank[Sa[i-1]]+(x[Sa[i]]!=x[Sa[i-1]]);
    for(int i=1;i<=len;i<<=1)
    {
        for(int j=0;j<len;j++)
            x[j]=rank[j],y[j]=j+i<len?rank[j+i]:0,Sa[j]=j;
        sort(Sa,rank,y,len,len),sort(rank,Sa,x,len,len);
        rank[Sa[0]]=1;
        for(int j=1;j<len;j++)
            rank[Sa[j]]=rank[Sa[j-1]]+(x[Sa[j]]!=x[Sa[j-1]]||y[Sa[j]]!=y[Sa[j-1]]);
        if(rank[Sa[len-1]]==len)  return ;
    }
}
void make_height()
{
    for(int i=0,j=0;i<len;i++)
    {
        if(j>0)  j--;
        if(rank[i]!=1)
            while(a[i+j]==a[Sa[rank[i]-2]+j])
                j++;
        height[rank[i]]=j;
    }
}
void print()
{
    for(int i=0;i<len;i++)
        if(Sa[i]<len/2)
            printf("%c",a[Sa[i]+len/2-1]);
    printf("\n");
}
int main()
{
    read();
    make_Sa();
    make_height();
    print();
    return 0;
}
  评论这张
 
阅读(544)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017