poj 1743 Musical Theme

主要是使用后缀数组，参考《后缀数组----处理字符串的有力工具》

sa[i]表示排名第i的后缀是哪个

可以得到height[i]数组：表示sa[i]和sa[i-1]的子串的最大长度。

然后怎么找到最大的解呢，首先就是二分答案k，然后按照顺序按照hight[i]>=k分组，每一组里面j个height里面总共j+1个sa，如果里面最小和最大的差>=k，那么就找到了长度为k并且不重叠的公共字符串。

不过题目要求的是差值一定的公共子串，如何做呢？

先预处理输入的数据，i=0..n-2,r[i] = data[i+1]-data[i]+88，r[n-1] = 0;

为什么这样处理呢

假设data[j]-data[i]=data[j+1]-data[i+1]=……data[j+k]=data[i+k]

移一下项就知道了data[i+1]-data[i] = data[j+1]-data[j]

#include<iostream>

#include<cstdio>
#include<cstring>
#include<cstring>
using namespace std;
const int MAXN = 21000;
const int MAXM = 90;
int r[MAXN];
int sa[MAXN];
int ran[MAXN];
int height[MAXN];
int wa[MAXN],wb[MAXN];
int wv[MAXN],wws[MAXN];
int low,high;
int cmp(int *r,int a,int b,int l){
return r[a]==r[b] && r[a+l]==r[b+l];
}
void da(int *r,int *sa,int n,int m){
int i,j,p,*x = wa,*y = wb,*t;
for(i=0;i<m;i++)wws[i] = 0;
for(i = 0;i<n;i++)wws[x[i]=r[i]]++;
for(i = 1;i<m;i++)wws[i]+=wws[i-1];
for(i=n-1;i>=0;i--)sa[--wws[x[i]]] = i;
for(j=1,p=1;p<n;j*=2,m=p){
for(p=0,i=n-j;i<n;i++)y[p++]=i;
for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
for(i=0;i<n;i++)wv[i]=x[y[i]];
for(i=0;i<m;i++)wws[i]=0;
for(i=0;i<n;i++)wws[wv[i]]++;
for(i=1;i<m;i++)wws[i]+=wws[i-1];
for(i=n-1;i>=0;i--)sa[--wws[wv[i]]]=y[i];
for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
}
return ;
}
void calheight(int *r,int *sa,int n){
int i,j,k = 0;
r[n] = 0;
for(i = 0;i<n;i++)ran[sa[i]] = i;
for(i = 0;i<n;height[ran[i++]] = k,k>high?high = k:0){
if(ran[i]==0){
k = 0;
continue;
}
for(k?k--:0,j = sa[ran[i]-1];r[i+k]==r[j+k];){
k++;
}
}
return ;
}
int n;
bool binary(int mid){
int min_,max_;
for(int i = 1;i<n;i++){
if(height[i]>=mid){
min_ = max_ = sa[i-1];
while(i<n && height[i]>=mid){
if(sa[i]<min_)min_ = sa[i];
if(sa[i]>max_)max_ = sa[i];
if(max_-min_>=mid)return true;
i++;
}
if(i<n)i--;

}
}
return false;
}
int main(){

while(scanf("%d",&n),n){
int data[MAXN];
for(int i = 0;i<n;i++)scanf("%d",&data[i]);
for(int i = 0;i<n-1;i++)r[i] = data[i+1]-data[i]+88;
r[n-1] = 0;
da(r,sa,n,200);
///for(int i = 0;i<n;i++)cout <<sa[i]<<endl;
low = high = 0;
calheight(r,sa,n);
//for(int i = 0;i<n;i++)printf("height[%d]=%d\n",i,height[i]);
///low = 0;
int mid;
int ans = 0;
while(low<=high){
mid = (low+high)/2;
if(binary(mid)){
ans = mid;
low = mid+1;
}else {
high = mid-1;
}
}
if(ans+1>=5)printf("%d\n",ans+1);
else printf("0\n");
}


return 0;
}