Poj 1743Musical Theme|后缀数组|二分答案

首先差分一下

话说为什么我花了好久才弄懂题意，然后才知道为什么差分……

然后就是二分答案，给height分组judge

1A赞（调了一早晨的sb错误）

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define T 511111
using namespace std;
int sc()
{
	int i=0;char c=getchar();
	while(c>'9'||c<'0')c=getchar();
	while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
	return i;
}
int s[T];
int t1[T],t2[T],cc[888888],sa[T],rank[T],height[T];
int len;
bool cmp(int *y,int a,int b,int k)
{
	int a1=y[a],b1=y[b];
	int a2=a+k>=len?-1:y[a+k];
	int b2=b+k>=len?-1:y[b+k];
	return a1==b1&&a2==b2;
}
void make_sa()
{
	int *x=t1,*y=t2,m=333;
	for(int i=0;i<m;i++)cc[i]=0; 
	for(int i=0;i<len;i++) ++cc[x[i]=s[i]+100];
	for(int i=1;i<m;i++)cc[i]+=cc[i-1];
	for(int i=len-1;~i;i--)sa[--cc[x[i]]]=i;
	for(int k=1;k<=len;k<<=1)
	{
		int p=0;
		for(int i=len-k;i<len;i++)y[p++]=i;
		for(int i=0;i<len;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
		for(int i=0;i<m;i++)cc[i]=0;
		for(int i=0;i<len;i++)++cc[x[y[i]]];
		for(int i=1;i<m;i++)cc[i]+=cc[i-1];
		for(int i=len-1;~i;i--)sa[--cc[x[y[i]]]]=y[i];
		swap(x,y),m=1,x[sa[0]]=0;
		for(int i=1;i<len;i++)x[sa[i]]=cmp(y,sa[i],sa[i-1],k)?m-1:m++;
		
		if(m>=len)return;
	}
}		
void make_height()
{
	int k=0;
	for(int i=0;i<len;i++)rank[sa[i]]=i;
	for(int i=0;i<len;i++)
	{
		if(!rank[i])continue;
		int j=sa[rank[i]-1];
		if(k)k--;
		while(s[i+k]==s[j+k])k++;
		height[rank[i]]=k;
	}
}
bool jud(int x)
{
	int now=0,l,r,lt;
	x--;
	while(now<len)
	{
		l=sa[now],r=sa[now];lt=now; now++;
		while(height[now]>=x&&now<len)
		{
			l=min(l,sa[now]);
			r=max(r,sa[now]);
			now++;
		}
		if(now!=lt+1)now--;
		if(r-l>=x)return 1;
	}
	return 0;
}	
int main()
{
	len=sc();
	while(len)
	{
		for(int i=0;i<len;i++)s[i]=sc();
		for(int i=len-1;i;i--)s[i]-=s[i-1];
		if(len<10)
		{
			puts("0");
			len=sc();
			continue;
		}
		make_sa(),make_height();
		int l=5,r=len/2,ans=0;
		while(l<=r)
		{
			int mid=l+r>>1;
			if(jud(mid))ans=mid,l=mid+1;
			else r=mid-1;
		}
		printf("%d\n",ans);
		len=sc();
	}
	return 0;
}