poj 1743 二分答案+后缀数组 求不重叠的最长重复子串

题意：给出一串序列，求最长的theme长度

（theme：完全重叠的子序列，如1 2 3和1 2 3  or  子序列中每个元素对应的差相等，如1 2 3和7 8 9）

要是没有差相等这个条件那就好办多了，直接裸题。

一开始想了个2B方法，后来发现真心2B啊蛤蛤蛤

 for i= to  do
 {
     for j= to length
     {
         r2[j]=r[j]+i;
         if (r2[j]>）    r2[i]-=;
     }
     把新序列r2连接到原序列r的后面
     process[r+r2]
     .....
     求for i=->88中得到的最大答案
 }

正解：http://bbezxcy.iteye.com/blog/1407395

对原序列相邻两元素作差就行了= =

比如：序列8 8 8 15 24 3 3 3

作差得到0 0 7 9 -21 0 0

对新序列再求最长不重叠重复子串，得到的结果+1，就是原序列的最长不重叠重复子串

求最长不重叠重复子串：参考IOI2009论文

Code：

 #include "stdio.h"
 #include "string.h"
 #define maxn 20010

 int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
 int rank[maxn],height[maxn];
 int r[maxn],sa[maxn],ans[maxn];
 int n;

 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=; i<m; i++) ws[i]=;
     for(i=; i<n; i++) ws[x[i]=r[i]]++;
     for(i=; i<m; i++) ws[i]+=ws[i-];
     for(i=n-; i>=; i--) sa[--ws[x[i]]]=i;
     for(j=,p=; p<n; j*=,m=p)
     {

         for(p=,i=n-j; i<n; i++) y[p++]=i;
         for(i=; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;
         for(i=; i<n; i++) wv[i]=x[y[i]];
         for(i=; i<m; i++) ws[i]=;
         for(i=; i<n; i++) ws[wv[i]]++;
         for(i=; i<m; i++) ws[i]+=ws[i-];
         for(i=n-; i>=; i--) sa[--ws[wv[i]]]=y[i];
         for(t=x,x=y,y=t,p=,x[sa[]]=,i=; i<n; i++)
             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
     }
     return;
 }

 void calheight(int *r,int *sa,int n)
 {
     int i,j,k=;
     for(i=; i<=n; i++) rank[sa[i]]=i;
     for(i=; i<n; height[rank[i++]]=k)
         for(k?k--:,j=sa[rank[i]-]; r[i+k]==r[j+k]; k++);
     return;
 }

 bool calc(int x)
 {
     int nm=;
     for (int i=;i<=n;i++)
         if (height[i]<x)
         {
             nm++;
             ans[nm]=i;
         }
     nm++;
     ans[nm]=n+;
     for (int i=;i<nm;i++)
     {
         int tl=ans[i],tr=ans[i+]-;
         int tn=maxn,tx=;
         for (int j=tl;j<=tr;j++)
         {
             if (tn>sa[j])   tn=sa[j];
             if (tx<sa[j])   tx=sa[j];
         }
         if (tx-tn>=x)   return true;
     }
     return false;
 }

 //da(r,sa,n+1,128);
 //calheight(r,sa,n);
 int main()
 {
 //    freopen("in.txt","r",stdin);
     while (~scanf("%d",&n))
     {
         if (n!=)
         {
             memset(sa,,sizeof(sa));
             memset(height,,sizeof(height));
             memset(ans,,sizeof(ans));

             for (int i=; i<n; i++)
                 scanf("%d",&r[i]);
             r[n]=;

 //          for (int i=0;i<=n;i++)
 //              printf("%d ",r[i]);
 //          printf("\n %d\n",n);

             for (int i=;i<n;i++)
                 r[i-]=r[i]-r[i-];
             n--;
             r[n]=;

             for (int i=;i<n;i++)
                 r[i]+=;

 //          for (int i=0;i<=n;i++)
 //              printf("%d ",r[i]);
 //          printf("\n %d\n",n);

             da(r,sa,n+,);
             calheight(r,sa,n);

 //          for (int i=0; i<=n; i++)
 //              printf("%d  %d\n",sa[i],height[i]);
 //          printf("\n");

             int l=maxn,r=,res=;
             for (int i=;i<=n;i++)
             {
                 if (height[i]<l)   l=height[i];
                 if (height[i]>r)   r=height[i];
             }
             while (r>=l)
             {
                 int mid=(l+r)/;
                 if (calc(mid))
                 {
                     l=mid+;
                     res=mid;
                 }
                 else
                     r=mid-;
             }
             res++;
             if (res<)  res=;
             printf("%d\n",res);
         }
     }
     return ;
 }

本题出自USACO5.1 theme，据说是男人八题之一.....orz