题意:
给出n个数字的串,将这个串后一个减去前一个,所得的串尽心操作。问所得的串最长不重复子串的长度。
题解:
这种题目很明显没办法直接求得长度,那么可以考虑二分答案判断是否满足条件。我们可以二分得到k,判断是否存在这样的k。对于如何判断,我们可以将后缀根据height的值进行分组,每组的后缀之间的最长公共前缀必须至少为k,这样之后只要找出每组中最大SA值和最小SA的值的差即可。好好理解下,这种做法,值得细细考虑。
#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=200005;
int rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],cnt[maxn];
int a[maxn];
void Swap(int*& x,int*& y)
{
int *temp=x;
x=y;
y=temp;
}
int lcp(int t[],int a,int b,int l)
{
return t[a]==t[b]&&t[a+l]==t[b+l];
}
void build_SA(int s[],int len,int up){
int *k1=t1,*k2=t2,*r=t3;
for(int i=0;i<up;i++)cnt[i]=0;
for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;
for(int d=1,p=0;p<len;d<<=1){
for(int i=len-d;i<len;i++)k2[p++]=i;
for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
for(int i=0;i<len;i++)r[i]=k1[k2[i]];
for(int i=0;i<up;i++)cnt[i]=0;
for(int i=0;i<len;i++)cnt[r[i]]++;
for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];
Swap(k1,k2);
k1[SA[0]]=0;
p=1;
for(int i=1;i<len;i++){
k1[SA[i]]= lcp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
}
if(p>=len)return ;///这句话切记要添加,忘了,检查了好几个小时。ORZ!
up=p;
p=0;
}
}
void get_height(int s[],int len)
{
for(int i=1;i<=len;i++)rank[SA[i]]=i;
for(int i=0,p=0;i<len;i++){
int j=SA[rank[i]-1];
while(s[i+p]==s[j+p])p++;
height[rank[i]]=p;
if(p)p--;
}
}
bool ok(int tag,int len)
{
int i=1;
int mx,mi;
while(1)
{
while(i<len&&height[i]<tag)i++;
if(i>=len)break;
mx=SA[i-1];
mi=SA[i-1];
while(i<len&&height[i]>=tag){
mx=max(mx,SA[i]);
mi=min(mi,SA[i]);
i++;
}
if(mx-mi>=tag)return true;
}
return false;
}
int solve(int len)
{
int l=0,r=len,mid;
int ans=0;
while(l<=r){
mid=(l+r)>>1;
if(ok(mid,len)){
l=mid+1;
ans=mid;
}else{
r=mid-1;
}
}
return ans;
}
int main(){
int n;
while(scanf("%d",&n)!=EOF){
if(n==0)break;
for(int i=0;i<n;i++)
scanf("%d",&a[i]);
n--;
for(int i=0;i<n;i++)
a[i]=a[i+1]-a[i]+90;
a[n]=0;
build_SA(a,n+1,200);
get_height(a,n);
int ans=solve(n);
if(ans>=4)ans++;
else ans=0;
printf("%d\n",ans);
}
return 0;
}//main
/**
*/