POJ 1743 - Musical Theme 后缀数组模板+用后缀数组求不覆盖最长重复子序列

             题意:

                     给了N个数(1<=N<=20000)..现在问其中有最长的相同的子串是多长..首先要没有重合..再一个子串可以是一段全部加/减一个数和另一个子串比较..

             题解:

                     熟悉模板...

Program:

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<set>
#include <stack>
#include<queue>
#include<algorithm>
#include<cmath>
#define eps 1e-4
#define oo 1000000007
#define MAXN 20105
#define MAXM 100005<<1
#define ll long long
#define pi acos(-1.0) 
using namespace std;    
//----------------后缀数组模板-------------------- 
//sa[]是存好的后缀排序起点 
//rank[]是存着每个后缀的排序后位置(后面重新统计得到)
//hight[]存的是排好序后前后两个后缀最长公共前缀 
//后缀数组中的所有元素是非负数 
const int maxn = 20105;  
int rank[maxn],wb[maxn],wv[maxn],wss[maxn];
int n;    
bool 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)  // n串的总长+1,m字符集的个数 
{   
      int i,j,p,*x=rank,*y=wb,*t;  
      for(i=0;i<m;i++) wss[i]=0;  
      for(i=0;i<n;i++) wss[x[i]=r[i]]++;  
      for(i=1;i<m;i++) wss[i]+=wss[i-1];  
      for(i=n-1;i>=0;i--) sa[--wss[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++) wss[i]=0;  
            for(i=0;i<n;i++) wss[wv[i]]++;  
            for(i=1;i<m;i++) wss[i]+=wss[i-1];  
            for(i=n-1;i>=0;i--) sa[--wss[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++;  
      }   
      return;  
}  
int height[maxn];  
void calheight(int *r,int *sa,int n)  // 求height,n为正整数 
{  
      int i,j,k=0;  
      for(i=1;i<=n;i++) rank[sa[i]]=i; //得到的sa[]下标是从1开始的
      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++);
      return;  
}  
//----------------后缀数组模板-------------------- 
int s[maxn],sa[maxn];
bool ok(int len,int n)
{
      int i,left,right;
      left=right=0;
      for (i=0;i<n;i++)
         if (height[i]<len) left=right=sa[i];
           else
           {
                  left=min(left,sa[i]);
                  right=max(right,sa[i]);
                  if (right-left>=len) return true;  
           }
      return false; 
}
int main()
{         
      int n,i,l,r,mid;   
      while (~scanf("%d",&n) && n)
      {
               n--;
               for(i=0;i<=n;i++) scanf("%d",&s[i]); 
               for(i=0;i<n;i++)  s[i]=s[i+1]-s[i]+100;
               s[n]=0;
               da(s,sa,n+1,200);
               calheight(s,sa,n);
               l=0,r=n;
               while (r-l>1)
               {
                       mid=r+l>>1;
                       if (ok(mid,n)) l=mid;
                                 else r=mid;
               }
               l++;
               if (l<5) l=0;
               printf("%d\n",l);
      }
      return 0;
}