poj1743 后缀数组入门题

版权声明：欢迎转载，转载请注明出处，如有错误，还望指出，谢谢 博客地址：https://blog.csdn.net/lanyanzhiji123asd https://blog.csdn.net/lanyanzhiji123asd/article/details/89462677

题目链接：http://poj.org/problem?id=1743

//poj1743

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<set>
#include<queue>
using namespace std;
#define ll long long
const int ding=2e4+5;
const int INF=0x3f3f3f3f;
int sa[ding],x[ding],y[ding],c[ding],r[ding],rank[ding],height[ding];
int n,m;
void get_SA()
{
	int i,k;
	memset(c,0,sizeof(c));//清空桶；
	for(i=0;i<n;i++)	c[x[i]=r[i]]++; 
	for(i=1;i<m;i++)	c[i]+=c[i-1];
	for(i=n-1;i>=0;i--)	
	{
		sa[--c[x[i]]]=i;
	}
	int p=1;
	for(k=1;k<n;k<<=1,m=p)
	{
		p=0;
		for(i=n-k;i<n;i++)	y[p++]=i;//????????????????????????????
		for(i=0;i<n;i++)	if(sa[i]>=k) y[p++]=sa[i]-k;//大于等于k 
		//sa是第一关键字的排名， 然后再根据第二关键字的排名来确定排名， 
		memset(c,0,sizeof(c));
		for(i=0;i<n;i++)	c[x[y[i]]]++;//统计当前字符串 
		for(i=1;i<m;i++)	c[i]+=c[i-1];//统计个数 
		for(i=n-1;i>=0;i--)		sa[--c[x[y[i]]]]=y[i];//重置sa数组
		memcpy(y,x,sizeof(x));
		memset(x,0,sizeof(x));
		p=1;//从1开始编号 
		x[sa[0]]=0;
		//因为上面已经重置了x的值  所以下面是从1开始循环，然后是和前一个比较sa[i]和sa[i-1] 比较 
		for(i=1;i<n;i++) 
			x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k])?p-1:p++;
		
		//x数组里面记住的是x[i]=j  下标为i的后缀，排名为j
		//y数组的内容是y[i]=j 下标为j 的后缀排名为i 
		if(p>=n)break;
	}
}

void get_height()
{
	int i,k=0,j;
	for(i=0;i<n;i++) rank[sa[i]]=i;
	for(i=0;i<n-1;i++)
	{//为什么循环到n-1，因为我们在原始字符串后面加了a，所以这个后缀一定是最小的。 
		if(k)k--;
		for(j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
		height[rank[i]]=k; 
	}
		
} 

int check(int x)
{
	int minn=INF,maxn=-INF,i;
	for(i=2;i<n;i++)
	{
		if(height[i]>=x)
		{
			minn=min(minn,min(sa[i],sa[i-1]));
			maxn=max(maxn,max(sa[i],sa[i-1]));
			if(maxn-minn>x)return 1;
		}
		else
		{
			minn=INF;
			maxn=-INF;
		}
	}
	return 0;
}

int main()
{
 	int i,L,R,ans;
	
	while(~scanf("%d",&n)&&n){
		m=180;
		for(i=0;i<n;i++)
		{
			scanf("%d",&r[i]);
		}
		for(i=0;i<n-1;i++)
			r[i]=r[i+1]-r[i]+89;
		r[n-1]=0;
		get_SA(); 
		get_height();
		L=0,R=(n-1)>>1;
		while(L<R)
		{
			int mid=(L+R+1)>>1;
			if(check(mid))
			{
				L=mid;
				
			}
			else R=mid-1;
		}
		if(L>=4)printf("%d\n",L+1);
		else printf("0\n");
	}
	

 return 0;
}