传说中的楼教主男人必做八题。。。
这题啊。。。凭借我当时对后缀数组的理解,那是肯定做不出来的。。。
要用到二分判定答案以及height数组很重要的一个性质,分块法。
题目抽象一下后就变成了:求一个序列中最长的不重叠子串。
因为是不重叠的,所以不能直接遍历一遍height数组来求最大值。
height数组有着这样的一个性质,我总是把他理解成传递性。如果在区间[l,r)中,height值均大于等于k,那么,对sa数组中[l,r]对应的后缀中任意两个,他们的lcp均大于等于k。这个还是蛮容易理解的。
所以这题我们可以二分答案,每一次二分答案的时候,需要遍历一遍height数组,看是否存在这样的一个区间,rank值满足要求。
PS:我的程序中height数组遍历是要用两个for循环的,其实是可以变成1个的。
#include "cstdio"
#include "cstring"
#include "cmath"
#define max 200005
int a[max];
int sa[max], rank[max], height[max];
int wa[max], wb[max], wv[max], wd[max];
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 n, int m){// 倍增算法 r为待匹配数组 n为总长度 m为字符范围
int i, j, p, *x = wa, *y = wb, *t;
for(i = 0; i < m; i ++) wd[i] = 0;
for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++;
for(i = 1; i < m; i ++) wd[i] += wd[i-1];
for(i = n-1; i >= 0; i --) sa[-- wd[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 ++) wd[i] = 0;
for(i = 0; i < n; i ++) wd[wv[i]] ++;
for(i = 1; i < m; i ++) wd[i] += wd[i-1];
for(i = n-1; i >= 0; i --) sa[-- wd[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 ++;
}
}
}
void calHeight(int *r, int n){ // 求height数组。
int i, j, k = 0;
for(i = 1; i <= n; i ++) rank[sa[i]] = i;
for(i = 0; i < n; height[rank[i ++]] = k){
for(k ? k -- : 0, j = sa[rank[i]-1]; r[i+k] == r[j+k]; k ++);
}
}
int cal(int x){
return x>0 ? x:-x;
}
int judge(int mid,int n){
for(int i=2;i<=n;i++){
if(height[i]<mid) continue;
for(int j=i-1;j>=2;j--){
if(cal(sa[i]-sa[j])>=mid){
return 1;
}
if(height[j]<mid) break;
}
}
return 0;
}
int main(){
int n;
while(scanf("%d",&n)!=EOF&&n!=0){
for(int i=0;i<n;i++) scanf("%d",&a[i]);
if(n<10){
printf("0\n");
continue;
}
for(int i=0;i<n-1;i++) a[i]=a[i+1]-a[i]+90;
a[n-1]=0,n--;
//printf("%d\n",n);
//for(int i=0;i<n;i++) printf("%d ",a[i]);
//printf("\n");
da(a,n+1,200);
calHeight(a,n);
int ans=0;
int left=3,right=n,mid;
while(right>=left){
mid=(left+right)/2;
if(judge(mid,n)){
left=mid+1;
ans=mid;
}else{
right=mid-1;
}
}
if(ans<4){
printf("0\n");
}else{
printf("%d\n",ans+1);
}
}
return 0;
}