【后缀数组】不可重叠最长重复子串

hz2016评测《《点击访问caioj《《点击访问这题阴森可怕，两个后缀可以不相同？哦天哪，怎么判断呢。于是，某个聪明的人发现，后缀中两两字符的ANSI值的差是相等的，于是我们可以维护一个由差值组成的后缀数组。所以呢，我们就可以用后缀数组解这道题了，但是这一题，我们还要引入一个新的概念height数组和h数组。height数组的定义就是排名为 i 与 i-1 的最长公共前缀。h就是代表排名为 i  的后缀与其前一个后缀的最长前缀
#include<map>#include<queue>#include<cmath>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#define Maxchar 28192 //20000 +2^13#define Maxs 30#define mes(x,y) memset(x,y,sizeof(x));#define mpy(x,y) memcpy(x,y,sizeof(x))#define INF 2147483647using namespace std; int n,a[Maxchar+1],tt[Maxchar+1],Rank[Maxchar+1],sa1[Maxchar+1],Rank2[Maxchar+1],sa2[Maxchar+1],Rsort[Maxchar+1];void get_sa(int n,int m){    memcpy(Rank,a,sizeof(Rank));    memset(Rsort,0,sizeof(Rsort));    for(int i=1;i<=n;i++)Rsort[Rank[i]]++;    for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];    for(int i=n;i>=1;i--)sa1[Rsort[Rank[i]]--]=i;    int ln=1,p=0;    while(p<n){        int k=0;        for(int i=n-ln+1;i<=n;i++)sa2[++k]=i;        for(int i=1;i<=n;i++)if(sa1[i]-ln>0)sa2[++k]=sa1[i]-ln;        memset(Rsort,0,sizeof(Rsort));        for(int i=1;i<=n;i++)Rsort[Rank[i]]++;        for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];         for(int i=n;i>=1;i--)sa1[Rsort[Rank[sa2[i]]]--]=sa2[i];        for(int i=1;i<=n;i++)tt[i]=Rank[i];        p=1;Rank[sa1[1]]=1;        for(int i=2;i<=n;i++){            if(tt[sa1[i]]!=tt[sa1[i-1]]||tt[sa1[i]+ln]!=tt[sa1[i-1]+ln])p++;            Rank[sa1[i]]=p;        }         m=p;ln*=2;    }} //上面是基础的后缀数组模版。 int height[Maxchar+1];void get_he(int n){    int j,k=0;    for(int i=2;i<=n;i++){        j=sa1[Rank[i]-1];        if(k!=0)k--;         while(a[j+k]==a[i+k])k++;//匹配最初前缀         //h数组的定义就是后缀 i 和后缀 i-1的最长前缀         height[Rank[i]]=k;//记录两个后缀的最长前缀         //最后记录下排名为 i 的后缀与其前一个后缀的最长前缀     }}//就是新的h数组，这个才是后缀数组的真谛。 bool check(int k,int n){    for(int i=2;i<=n;i++){        if(height[i]>=k){            for(int j=i-1;j>=2;j--){                if(abs(sa1[i]-sa1[j])>=k)return true;//满足直接true                 if(height[j]<k)break;//二分假如现在条件不满足就跳出             }        }    }    return false;}int main(){while(scanf("%d",&n)!=EOF){if(n==0)break;for(int i=1;i<=n;i++)scanf("%d",&a[i]);int mmax=-9999999;        for(int i=1;i<n;i++){             a[i]=a[i+1]-a[i]+88;            if(mmax<a[i])mmax=a[i];            //取maxx是为了减少后缀数组的基数排序时间复杂度         }        a[n]=0;n--;        get_sa(n,mmax);get_he(n);        //然后通过二分枚举答案         int l=1,r=n,ans=1;        while(l<=r){            int mid=(l+r)/2;            if(check(mid,n)==true){//二分枚举长度                 ans=mid;l=mid+1;            }            else r=mid-1;        }        if(ans<4)printf("0\n");         else printf("%d\n",ans+1);}    return 0;}
 

查看原文：http://hz2016.tk/blog/?p=29