poj 1743

若序列ai,ai+1,ai+2,ai+k-1,和aj,aj+1,aj+2,..aj+k-1,满足要求，那ai+1-ai==aj+1-aj,以此类推。所以对原序列进行处理，得到有n-1个元素的序列b，bi代表ai+1-ai的值，那么原问题就转化为了不可重叠最长重复子串。二分结果值k，看是否存在两个长度为k的子串相同，并且不重叠。也就是看是否存在一对后缀，这两个后缀的LCP大于等于k，并且编号差值>=k。把根据height数组分组，得到的每组任意两个后缀的LCP>=k，并且是封闭的,任意一对后缀满足LCP大于等于k的后缀都分别属于这几组，即分了类，然后看每一类中是否至少存在一对。即求出每一组中，sa[i]最大值和最小值的差值，看是否>=k，若存在这样的一个分组，则当前k成立，二分下一个k。

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    using namespace std;
    const int maxn=20000+10;
    int num[maxn];
    int sa[maxn],t[maxn],t2[maxn],c[maxn],height[maxn],rank[maxn];
    int n;
    void build_sa(int m,int n,int *s)
    {
        int i,*x=t,*y=t2;
        for(i=1;i<m;i++) c[i]=0;
        for(i=0;i<n;i++) c[x[i]=s[i]]++;
        for(i=2;i<m;i++) c[i]+=c[i-1];
        for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
        for(int k=1;k<=n;k<<=1)
        {
            int 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;
            for(i=1;i<m;i++) c[i]=0;
            for(i=0;i<n;i++) c[x[y[i]]]++;
            for(i=2;i<m;i++) c[i]+=c[i-1];
            for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
            swap(x,y);
            p=2;x[sa[0]]=1;
            for(i=1;i<n;i++) x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
            if(p>n) break;
            m=p;
        }
    }
    void getHeight(int n)
    {
        int i,j,k=0;
        for(i=0;i<n;i++) rank[sa[i]]=i;
        for(i=0;i<n;i++)
        {
            if(k) k--;
            if(rank[i])
            {
                int j=sa[rank[i]-1];
                while(num[i+k]==num[j+k]) k++;
                height[rank[i]]=k;
            }
        }
    }
    int bsearch(int n)
    {
        int low=3,high=n+1,mid,i,maxv,minv;
        int flag;
        while(high>low)
        {
            flag=0;
            mid=low+(high-low)/2;
            maxv=minv=sa[0];
            for(i=1;i<n;i++)
            {
                if(height[i]>=mid)
                {
                    if(sa[i]>maxv) maxv=sa[i];
                    if(sa[i]<minv) minv=sa[i];
                }
                else
                {
                    if((maxv-minv)>=mid)
                    {
                        flag=1;
                        break;
                    }
                    maxv=minv=sa[i];
                }
            }
            if((maxv-minv)>=mid) flag=1;
            if(flag) low=mid+1;
            else high=mid;
        }
        return (low-1)>3?low:0;
    }
    int main()
    {
        //freopen("1.txt","r",stdin);
        while(scanf("%d",&n)&&n)
        {
            int i,tem,v;
            scanf("%d",&tem);
            for(i=0;i<n-1;i++)
            {
                scanf("%d",&v);
                num[i]=v-tem+88;
                tem=v;
            }
            build_sa(88+88,n-1,num);
            getHeight(n-1);
            int ans=bsearch(n-1);
            printf("%d\n",ans);
        }
        return 0;
    }