poj2774 Long Long Message（后缀数组）

来源：互联网发布：淘宝ulzzang短发女模特编辑：程序博客网时间：2024/05/29 23:23

很简单的后缀数组，求两个串的最长公共子串。

把这两个串，用特殊字符隔开连在一起，更新合法的最大公共前缀即可。

#include <cstdio>#include <cstring>#define N 200005char a[N];int rank[N<<1],sa[N],st[256],rank1[N],count[N],tmp[N],ans=0,h[N];inline int max(int x,int y){return x>y?x:y;}int main(){//freopen("a.in","r",stdin);scanf("%s",a+1);int len=strlen(a+1);a[++len]='z'+1;scanf("%s",a+len+1);int n=strlen(a+1);memset(st,0,sizeof(st));memset(rank,0,sizeof(rank));memset(rank1,0,sizeof(rank1));for(int i=1;i<=n;++i) st[a[i]]=1;for(int i=1;i<=255;++i) st[i]+=st[i-1];for(int i=n;i>=1;--i) rank[i]=st[a[i]];int k=0;for(int l=1;k!=n;l<<=1){memset(count,0,sizeof(count));for(int i=1;i<=n;++i) count[rank[i+l]]++;for(int i=1;i<=n;++i) count[i]+=count[i-1];for(int i=n;i>=1;--i) tmp[count[rank[i+l]]--]=i;memset(count,0,sizeof(count));for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;for(int i=1;i<=n;++i) count[i]+=count[i-1];for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];memcpy(rank1,rank,sizeof(rank)>>1);rank[sa[1]]=k=1;for(int i=2;i<=n;++i){if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+l]!=rank1[sa[i-1]+l]) ++k;rank[sa[i]]=k;}}k=0;for(int i=1;i<=n;++i){//i--第i个后缀 if(rank[i]==1){h[1]=0;continue;}if(i==1||h[rank[i-1]]<=1) k=0;if(k) --k;while(a[i+k]==a[sa[rank[i]-1]+k]) ++k;h[rank[i]]=k;}for(int i=2;i<=n;++i){//i--排名第i if(!(sa[i]>len^sa[i-1]<=len)) ans=max(ans,h[i]);}printf("%d\n",ans);return 0;}

阅读全文

0 0