poj Musical Theme 后缀数组

二分长度 然后判断它是否存在且设两个指针存储它可能出现最大位置和最小位置
第一发后缀数组

#include<cstdio>
#include<iostream>
#include<algorithm>
#define maxn 20005
using namespace std;
int wa[maxn];//第一关键字
int wb[maxn];//第二关键字
int wv[maxn];//储存按第二关键字排好后的第一关键字
int wscnt[maxn];//计数排序用的辅助数组
int sa[maxn],ran[maxn],height[maxn];//后缀数组，排名数组，最长公共前缀数组
int s[maxn],n;//本题数据
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void mda(int *r,int *sa,int n,int m)//n数组长度，m基数排序记录数的最大范围
{
     int i,j,p,*x = wa,*y = wb,*t; // x 合并后的价值，y此时各个区间的价值
     for(i=0;i<m;i++)wscnt[i] = 0;
     for(i=0;i<n;i++)wscnt[x[i] = r[i]]++;
     for(i=1;i<m;i++)wscnt[i]+=wscnt[i-1];
     for(i=n-1;i>=0;i--)sa[--wscnt[x[i]]] = i;

     for(int 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++)wscnt[i] = 0;
         for(i=0;i<n;i++)wscnt[wv[i]]++;
         for(i=1;i<m;i++)wscnt[i]+=wscnt[i-1];
         for(i=n-1;i>=0;i--)sa[--wscnt[wv[i]]] = y[i];

         for(t=x,x = y,y = t,p=1,x[sa[0]]=0,i=1;i<n;i++)
             x[sa[i]] = cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
void getHeight(int *r,int n)//n不保存最后的0
{
     int i,j,k = 0;
     for(int i=0;i<=n;i++) ran[sa[i]]  = i;
     for(int i=0;i<n;i++)
     {
         if(k)k--;
         else k = 0;
         j = sa[ran[i]-1];
         while(r[i+k]==r[j+k])k++;
         height[ran[i]] = k;
     }
}
bool ok(int hi)
{
     int i,maxi,mini;
     maxi = mini = sa[0];
     for(int i=0;i<=n;i++)
     {
        if(height[i]>=hi)
        {
           maxi = max(maxi,sa[i]);
           mini = min(mini,sa[i]);
           if(maxi-mini>hi)return 1;
           continue ;
        }
        if(maxi-mini>hi)return 1;
        maxi = mini = sa[i];
     }
     return 0;
}
int main()
{
    while(scanf("%d",&n)&&n)
    {
         int ans = 0;
         for(int i=0;i<n;i++)scanf("%d",&s[i]);
         for(int i=0;i<n-1;i++)s[i] = s[i+1]-s[i]+100;
         s[--n] = 0;
         mda(s,sa,n+1,200);
         getHeight(s,n);
         int l = 0,r = n;
         while(l<=r)
         {
             int mid = (l+r)/2;
             if(ok(mid)){
                ans = mid;
                l = mid+1;
             }
             else r = mid-1;
         }
         ans++;
         if(ans<=4)printf("0\n");
         else printf("%d\n",ans);
    }
    return 0;
}