POJ 1743 后缀数组（男人八题七）

题意

有一串音符，输出最长的两个不重叠的相似串的长度（如果小于5，则输出0）。相似的意思是两个串所有字符之差必须一样。

题解

很有意思的一道题。据说可以用HASH做，但是HASH看了半天都没看懂，只能用后缀数组了。
刚看到相似的时候我吓了一跳，不过如果转化成字符之差的话就舒服多了。我们可以用后一个字符减去前一个字符，这样可以得到一个差序列。如果两个差序列相等，则可以认为两个串相似。
对于这种找最长字符串的题，一般的套路就是二分。这里我们也二分一下长度就可以了。有一点需要注意的是，mx-mn>k，因为处理的是字符之差，如果等于的话，那么两个字符串的首尾便会重叠，所以不能取到=。

代码

#include <iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<string>
#include<set>
#include<map>
#include<bitset>
#include<stack>
#include<string>
#define UP(i,l,h) for(int i=l;i<h;i++)
#define DOWN(i,h,l) for(int i=h-1;i>=l;i--)
#define W(a) while(a)
#define MEM(a,b) memset(a,b,sizeof(a))
#define LL long long
#define INF 0x3f3f3f3f3f3f3f3f
#define MAXN 20050
#define MOD 1000000007
#define EPS 1e-3
using namespace std;

int a[MAXN];
int sa[MAXN],t[MAXN],t2[MAXN],c[MAXN],height[MAXN],xrank[MAXN],n;

void build_sa(int m) {
    int i,*x=t,*y=t2;
    UP(i,0,m) c[i]=0;
    UP(i,0,n) c[x[i]=a[i]]++;
    UP(i,1,m) c[i]+=c[i-1];
    DOWN(i,n,0) sa[--c[x[i]]]=i;
    for(int k=1; k<=n; k<<=1) {
        int p=0;
        UP(i,n-k,n) y[p++]=i;
        UP(i,0,n) if(sa[i]>=k) y[p++]=sa[i]-k;
        UP(i,0,m) c[i]=0;
        UP(i,0,n) c[x[y[i]]]++;
        UP(i,1,m) c[i]+=c[i-1];
        DOWN(i,n,0) sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=1;
        x[sa[0]]=0;
        UP(i,1,n) 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 i,j,k=0;
    UP(i,0,n) xrank[sa[i]]=i;
    UP(i,0,n-1) {
        if(k) k--;
        int j=sa[xrank[i]-1];
        W(a[i+k]==a[j+k]) k++;
        height[xrank[i]]=k;
//        cout<<i<<" "<<k<<endl;
    }
}

bool good(int k) {
    int mx=0,mi=INF;
//    cout<<k<<endl;
    UP(i,1,n) {
        if(height[i]>=k) {
            mi=min(mi,min(sa[i],sa[i-1]));
            mx=max(mx,max(sa[i],sa[i-1]));
            if(mx-mi>k) return true;
        }
    }
    return false;
}

int main() {
    W(~scanf("%d",&n)) {
        if(n==0) break;
        UP(i,0,n) scanf("%d",&a[i]);
        UP(i,0,n-1) {
            a[i]=a[i+1]-a[i]+88;
//            cout<<a[i]<<endl;
        }
        a[n]=0;
        build_sa(200);
        getHeight();
        int l=4,h=n-1,ans=-1;
        W(l<=h) {
            int m=(l+h)/2;
            if(good(m)) {
                l=m+1;
                ans=m;
            } else {
                h=m-1;
            }
        }
        printf("%d\n",ans+1);
    }
}