SPOJ687 Repeats （后缀数组+lcp）

求循环节个数最大的子串的循环节个数。很绕吧。。基本和poj3693相同，此题要简单些。因为我先看了hzw的poj3693的题解。。所以我也写的是反过来做一遍求l。先入为主了呢。。网上还有很多人写的其实不需要反过来做一遍去求l，大家可以去借鉴下。。我觉得时间上差不多。。所以就懒得改了呢。。（思路有问题还是参考3693吧。。）

#include <cstdio>
#include <cstring>
#define N 50010
int tst,n,m,log[N],ans=1;
inline int min(int x,int y){return x<y?x:y;}
inline int swap(int &x,int &y){int t=x;x=y;y=t;}
inline int max(int x,int y){return x>y?x:y;}
struct SA{
	int s[N],n,m,st[N][18];
	int rank[N<<1],rank1[N],h[N],sa[N],tmp[N],count[N];
	void clear(){
		memset(s,0,sizeof(s));
	}
	void init(){
		m=2;
		memset(rank,0,sizeof(rank1));
		for(int i=1;i<=m;++i) count[i]=0;
		for(int i=1;i<=n;++i) count[s[i]]=1;
		for(int i=1;i<=m;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) rank[i]=count[s[i]];
		int k=0;
		for(int p=1;k!=n;p<<=1,m=k){
			for(int i=1;i<=m;++i) count[i]=0;
			for(int i=1;i<=n;++i) count[rank[i+p]]++;
			for(int i=1;i<=m;++i) count[i]+=count[i-1];
			for(int i=n;i>=1;--i) tmp[count[rank[i+p]]--]=i;
			for(int i=1;i<=m;++i) count[i]=0;
			for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;
			for(int i=1;i<=m;++i) count[i]+=count[i-1];
			for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];
			memcpy(rank1,rank,sizeof(rank1));
			rank[sa[1]]=k=1;
			for(int i=2;i<=n;++i){
				if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+p]!=rank1[sa[i-1]+p]) ++k;
				rank[sa[i]]=k;
			}
		}k=0;
		for(int i=1;i<=n;++i){
			if(rank[i]==1){h[1]=0;continue;}
			if(i==1||h[rank[i-1]]<=1) k=0;
			if(k) --k;
			while(s[i+k]==s[sa[rank[i]-1]+k]) ++k;
			h[rank[i]]=k;
		}
		for(int i=1;i<=n;++i) st[i][0]=h[i];
		for(int i=1;i<=log[n];++i)
			for(int j=1;j<=n;++j)
				if(j+(1<<i-1)<=n) st[j][i]=min(st[j][i-1],st[j+(1<<i-1)][i-1]);
				else break;
	}
	inline int lcp(int x,int y){
		int l=rank[x],r=rank[y];
		if(l>r) swap(l,r);++l;
		int t=log[r-l+1];
		return min(st[l][t],st[r-(1<<t)+1][t]);
	}
}c[2];
void calc(int x){
	for(int i=1;i+x<=n;i+=x){
		if(c[0].s[i]==c[0].s[i+x]){
			int r=c[0].lcp(i,i+x),l=c[1].lcp(n-i+2,n-(i+x)+2);
			ans=max(ans,(r+l)/x+1);
		}
	}
}
int main(){
//	freopen("a.in","r",stdin);
	scanf("%d",&tst);log[0]=-1;
	for(int i=1;i<=N-10;++i) log[i]=log[i>>1]+1;
	while(tst--){
		scanf("%d",&n);c[0].n=c[1].n=n;c[0].clear();c[1].clear();
		for(int i=1;i<=n;++i){
			char str[1];
			scanf("%s",str);
			c[0].s[i]=str[0]-'a'+1;
			c[1].s[n-i+1]=c[0].s[i];
		} 
		c[0].init();c[1].init();
		ans=1;
		for(int i=1;i<=n;++i) calc(i);
		printf("%d\n",ans);
	}
	return 0;
}