pku1743

这个地址上的题解已经好详细了

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;  

}