题目
开始后缀数组啦,开心,什么都看不懂,但只要理解好几个数组的意思,抄起模板来随便过,呜呜呜,真的理解不了国家集训队的论文了,神仙高中生啊,照着看了几遍了还是只会抄代码
楼教主男人八题之一(我不配做男人啦)
找序列的最长不重叠子串
题解
太经典了,论文中写的比我的思考清晰多了(什么你不知道是那篇论文?你还学什么后缀数组)
代码
主函数的格式有点小问题,Devc太垃圾了,没有一键缩进,一键校正代码的功能(但是用起来是真)
#include<iostream>
#include<algorithm>
#include<cstdio>
#define N 40005
using namespace std;
const int maxn=40005;
#define INF (1<<30)
int t1[N],t2[N],sa[N],h[N],rk[N],c[10*N],a[N];
char s[N];
int m,n;
void calcsa(int n,int m){
int *x=t1,*y=t2,p=0,f=0;
for(int i=1;i<=m;i++)c[i]=0;
for(int i=1;i<=n;i++)c[x[i]=a[i]]++;
for(int i=1;i<=m;i++)c[i]+=c[i-1];
for(int i=n;i;i--)sa[c[x[i]]--]=i;
for(int i=1;i<=n&&p<=n;i<<=1){
p=0;
for(int j=n-i+1;j<=n;j++)y[++p]=j;
for(int j=1;j<=n;j++)if(sa[j]>i)y[++p]=sa[j]-i;
for(int j=1;j<=m;j++)c[j]=0;
for(int j=1;j<=n;j++)c[x[y[j]]]++;
for(int i=1;i<=m;i++)c[i]+=c[i-1];
for(int j=n;j;j--)sa[c[x[y[j]]]--]=y[j];
swap(x,y);x[sa[1]]=1;p=2;
for(int j=2;j<=n;j++)
x[sa[j]]=y[sa[j]]==y[sa[j-1]]&&y[sa[j]+i]==y[sa[j-1]+i]?p-1:p++;
m=p;
}
for(int i=1;i<=n;i++)rk[sa[i]]=i;
for(int i=1;i<=n;i++){
int j=sa[rk[i]-1];
if(f)f--;while(a[i+f]==a[j+f])f++;
h[rk[i]]=f;
}
}
int jiancha(int k)
{
int minm=INF,maxm=0;
for(int i=2;i<=n;i++)
{
//cout<<h[i]<<endl;
if(h[i]<k)
{
minm=INF;
maxm=0;
}
else
{
minm=min(minm,min(sa[i],sa[i-1]));
maxm=max(maxm,max(sa[i],sa[i-1]));
if(maxm-minm>=k)
{
return 1;
}
}
}
return 0;
}
int b[maxn];
int main()
{
int x1;
while(~scanf("%d",&n)&&n)
{
for(int i=1;i<=n;i++)
{
scanf("%d",&b[i]);
}
for(int i=2;i<=n;i++)
{
a[i-1]=b[i]-b[i-1]+88;
}
a[n]=0;
calcsa(n,88*2);
/*for(int i=1;i<=n-1;i++)
{
cout<<a[rk[i]]<<' ';
}
cout<<endl;
for(int i=1;i<=n-1;i++)
{
cout<<' '<<h[i]<<' ';
}
cout<<endl;
for(int i=1;i<=n-1;i++)
{
cout<<' '<<sa[i]<<' ';
}
cout<<endl;*/
int lf=0,rg=n;
while(lf<rg)
{
//cout<<lf<<endl;
int mid=(lf+rg+1)/2;
if(jiancha(mid))
{
//cout<<lf<<endl;
lf=mid;
}
else
{
rg=mid-1;
}
}
if(lf+1<5)
{
cout<<'0'<<endl;
}
else
{
cout<<lf+1<<endl;
}
}
}