poj2774 sa模版

学习地址：http://blog.csdn.net/yxuanwkeith/article/details/50636898

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=;
char s[maxn],s1[maxn],s2[maxn];
int ans,tax[maxn],tp[maxn],sa[maxn],ra[maxn],H[maxn],n,m;
void rsort(){
    for(int i=;i<=m;++i)tax[i]=;
    for(int i=;i<=n;++i)tax[ra[tp[i]]]++;
    for(int i=;i<=m;++i)tax[i]+=tax[i-];
    for(int i=n;i>=;--i)sa[tax[ra[tp[i]]]--]=tp[i];
}
int cmp(int x,int y,int w){return tp[x]==tp[y]&&tp[x+w]==tp[y+w];}
void pre(){
    for(int i=;i<=n;++i)ra[i]=s[i]-'a'+,tp[i]=i;
    m=;rsort();
    for(int w=,p=,i;p<n;w+=w,m=p){
    //p==n是说明已经得到各个后缀之间的大小关系，故可跳出循环
    //并且最后p一定能==n因为各个后缀长度不一样，一定能排出一个序；
        for(p=,i=n-w+;i<=n;++i)tp[++p]=i;
        for(i=;i<=n;++i)if(sa[i]>w)tp[++p]=sa[i]-w;
        rsort();swap(ra,tp);ra[sa[]]=p=;
    //这里换一下才方便由这一次的ra得到下一次的ra;
        for(i=;i<=n;++i)ra[sa[i]]=cmp(sa[i],sa[i-],w)?p:++p;
    }
    int j,k=;
    for(int i=;i<=n;H[ra[i++]]=k)
        for(k=k?k-:k,j=sa[ra[i]-];s[i+k]==s[j+k];++k);
}
int main(){
    scanf("%s%s",s1+,s2+);
    int len1=strlen(s1+),len2=strlen(s2+);
    for(int i=;i<=len1;++i)s[i]=s1[i];
    s[len1+]='a'+;
    for(int i=;i<=len2;++i)s[i+len1+]=s2[i];
    n=len1+len2+;
    pre();ans=;
    //for(int i=1;i<=n;++i)cout<<sa[i]<<' '<<H[i]<<endl;
    for(int i=;i<=n;++i){
        if((sa[i-]<=len1)!=(sa[i]<=len1)&&ans<H[i])ans=H[i];
    }
    cout<<ans;
    system("pause");
    return ;
}
/*
yeshowmuchiloveyoumydearmotherreallyicannotbelieveit
yeaphowmuchiloveyoumydearmother
*/