这个地址上的题解已经好详细了
http://wenku.baidu.com/view/87f5b7f80242a8956bece4df.html
不过本人觉得,加上这部分题解的话,会更好理解
若在假设重复子串的长度最多为L的限制下有解, 则对于任意一个比L小的限制L'<L, 也一定有解. 这就说明存在解的连续性, 这样就可以用二分查找答案长度L.给出一个关于LCP的定理LCP(SA[i], SA[j]) = RMQ(Height[i+1..j]). 由此, 若存在k, 满足Height[k] < L, 则对于所有i, j 满足i < k < j, 有LCP(SA[i], SA[j]) < L. 即公共长度至少为L的两个后缀, 不会跨过一个小于L的Height低谷k, 所以我们可以得到一些由这些低谷划分开的连续的段.
在某段内, 若存在i, j 满足SA[i]+L<SA[j], 则存在一个长度至少为L的2个相同不交迭子串. 实现时只要记录在每段内, 最大和最小的SA值即可.
--Amber大牛,《男人不容易系列Solution》
#include <iostream>
using namespace std;
#define MAXN 20010
int a[MAXN],b[MAXN],array[4][MAXN],*sa,*nsa,*rank1,*nrank,height[MAXN],n;
void make_sa(){
int i,k;
sa=array[0];
nsa=array[1];
rank1=array[2];
nrank=array[3];
memset(b,0,sizeof(b));
for(i=0;i<n;i++)
b[a[i]]++;
for(i=1;i<=256;i++)
b[i]+=b[i-1];
for(i=n-1;i>=0;i--)
sa[--b[a[i]]]=i;
for(rank1[sa[0]]=0,i=1;i<n;i++){
rank1[sa[i]]=rank1[sa[i-1]];
if(a[sa[i]]!=a[sa[i-1]])
rank1[sa[i]]++;
}
for(k=1;k<n && rank1[sa[n-1]]<n-1;k*=2){
for(i=0;i<n;i++)
b[rank1[sa[i]]]=i;
for(i=n-1;i>=0;i--)
if(sa[i]-k>=0)
nsa[b[rank1[sa[i]-k]]--]=sa[i]-k;
for(i=n-k;i<n;i++)
nsa[b[rank1[i]]--]=i;
for(nrank[nsa[0]]=0,i=1;i<n;i++){
nrank[nsa[i]]=nrank[nsa[i-1]];
if(rank1[nsa[i]]!=rank1[nsa[i-1]] || rank1[nsa[i]+k]!=rank1[nsa[i-1]+k])
nrank[nsa[i]]++;
}
int *t=sa;sa=nsa;nsa=t;
t=rank1;rank1=nrank;nrank=t;
}
}
void cal_height(){
int i,j,k;
for(k=0,i=0;i<n;i++){
if(rank1[i]==0)
height[rank1[i]]=0;
else{
for(j=sa[rank1[i]-1];a[i+k]==a[j+k];k++);
height[rank1[i]]=k;
if(k>0)
k--;
}
}
}
bool OK(int len){
int i,mn,mx;
mn=n;
mx=0;
for(i=1;i<n;i++){
if(height[i]<len){
mn=n;
mx=0;
}
else{
if(sa[i]>mx)
mx=sa[i];
if(sa[i]<mn)
mn=sa[i];
if(sa[i-1]>mx)
mx=sa[i-1];
if(sa[i-1]<mn)
mn=sa[i-1];
if(mx-mn>=len)
return true;
}
}
return false;
}
int main(){
int i,ans,l,r,len;
while(scanf("%d",&n) && n){
for(i=0;i<n;i++)
scanf("%d",&a[i]);
for(i=1;i<n;i++)
a[i-1]=a[i]-a[i-1]+88;
n--;
a[n++]=0;
make_sa();
cal_height();
ans=0;
l=0;r=n-1;
while(l<=r){
len=(l+r)/2;
if(OK(len)){
ans=len;
l=len+1;
}
else
r=len-1;
}
if(ans>=4)
printf("%d\n",ans+1);
else printf("0\n");
}
return 0;
}
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
int const N= 20100;
int wa[N], wb[N], ws[N], wv[N];
int rank[N], sa[N], height[N], r[N], n, k;
int cmp( int* r, int a, int b, int L ){
return r[a]== r[b] && r[a+ L]== r[b+ L];
}
void da( int* r, int* sa, int n, int m ){
int i, j, p, *x= wa, *y= wb, *t;
for( i= 0; i< m; ++i ) ws[i]= 0;
for( i= 0; i< n; ++i ) ws[ x[i]= r[i] ]++;
for( i= 1; i< m; ++i ) ws[i]+= ws[i-1];
for( i= n- 1; i>= 0; i-- ) sa[ --ws[ x[i] ] ]= i;
for( j= 1, p= 1; p< n; j*= 2, m= p ){
for( p= 0, i= n- j; i< n; ++i ) y[p++]= i;
for( i= 0; i< n; ++i )
if( sa[i]>= j ) y[p++]= sa[i]- j;
for( i= 0; i< n; ++i ) wv[i]= x[y[i]];
for( i= 0; i< m; ++i ) ws[i]= 0;
for( i= 0; i< n; ++i ) ws[ wv[i] ]++;
for( i= 1; i< m; ++i ) ws[i]+= ws[i-1];
for( i= n- 1; i>= 0; i-- ) sa[ --ws[ wv[i] ] ]= y[i];
t= x, x= y, y= t, p= 1; x[ sa[0] ]= 0;
for( i= 1; i< n; ++i )
x[ sa[i] ]= cmp( y, sa[i-1], sa[i], j )? p- 1: p++;
}
}
void callheight( int* r, int*sa, int n ){
int i, j, k= 0;
for( i= 1; i<= n; ++i ) rank[ sa[i] ]= i;
for( i= 0; i< n; height[ rank[i++] ]= k )
for( k?k--:0, j= sa[ rank[i]- 1]; r[i+k]== r[j+k]; k++ );
}
bool OK(int len)
{
int mx=0,mi=n;
for(int i=1;i<n;i++)
{
if(height[i]<len)
{
mx=0;mi=n;
}
else
{
if(sa[i]>mx)
mx=sa[i];
if(sa[i-1]>mx)
mx=sa[i-1];
if(sa[i]<mi)
mi=sa[i];
if(sa[i-1]<mi)
mi=sa[i-1];
if(mx-mi>=len)
return true;
}
}
return false;
}
int main(){
int i,ans,l,len,right;
while(scanf("%d",&n) && n){
for(i=0;i<n;i++)
scanf("%d",&r[i]);
for(i=1;i<n;i++)
r[i-1]=r[i]-r[i-1]+88;
n--; r[n++]=0;
da(r,sa,n+1,180);//最后一个参数表示r[]数组的最大值,数组ws[]的大小随最后一个参数而变
callheight(r,sa,n);
ans=0;
l=0;right=n-1;
while(l<=right){
len=(l+right)/2;
if(OK(len)){
ans=len;
l=len+1;
}
else
right=len-1;
}
if(ans>=4)
printf("%d\n",ans+1);
else printf("0\n");
}
return 0;
}