POJ 1743 Musical Theme(二分+后缀数组)

题意：

求不重叠的最长相同变化的子串的最长长度，比如1 2 3 4 5 6 7 8 9 10,最长长度为5，因为子串1 2 3 4 5 和 6 7 8 9 10变化都一样的

思路：

男人八题之一，然而poj上a了都一万多了。。。。
思路在论文里面讲的很清楚了（其实应该是一个区间最大sa减去最小sa大于k，而不是大于等于，dicuss有人说了这个问题，我也是这么认为的，不过因为poj数据太水导致都能过）。这里讲一个坑点，就是表示前后的差值的时候会出现负数，要手动加上88，而因为都是差值且数值都是1-88，所以m最大也不会超过N

错误及反思：

数据很水，第一次交的有个地方其实不对但是也算我ac了

代码：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int N = 20100;int sa[N],rak[N],height[N];int arr[N];void construct(const int *s,int n,int m) {    static int t1[N],t2[N],c[N];    int *x = t1,*y = t2;    int i,j,k,p,l;    for (i = 0; i < m; ++ i) c[i] = 0;    for (i = 0; i < n; ++ i) c[x[i] = s[i]] ++;    for (i = 1; i < m; ++ i) c[i] += c[i - 1];    for (i = n - 1; i >= 0; -- i) sa[--c[x[i]]] = i;    for (k = 1; k <= n; k <<= 1) {        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 = 0; i < m; ++ i) c[i] = 0;        for (i = 0; i < n; ++ i) c[x[y[i]]] ++;        for (i = 1; i < m; ++ i) c[i] += c[i - 1];        for (i = n - 1; i >= 0; -- i) sa[--c[x[y[i]]]] = y[i];        std::swap(x,y);        p = 1; x[sa[0]] = 0;        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;    }    for (i = 0; i < n; ++ i) rak[sa[i]] = i;    for (i = 0,l = 0; i < n; ++ i) {        if (rak[i]) {            j = sa[rak[i] - 1];            while (s[i + l] == s[j + l]) l++;            height[rak[i]] = l;            if (l) l--;        }    }}bool judge(int x,int n){    for(int i=2;i<=n;i++)    {        int minn=n,maxn=-1;        while(height[i]>=x)        {            minn=min(minn,sa[i-1]);            minn=min(minn,sa[i]);            maxn=max(maxn,sa[i-1]);            maxn=max(maxn,sa[i]);            i++;        }        if(maxn-minn>x) return true;    }    return false;}int main(){    int n;    while(scanf("%d",&n)&&n)    {        for(int i=0;i<n;i++)            scanf("%d",&arr[i]);        for(int i=0;i<n-1;i++)            arr[i]=arr[i+1]-arr[i]+88;        construct(arr,n,20100);        int l=0,r=n-1;        while(l+1<r)        {            int m=(l+r)>>1;            if(judge(m,n-1)) l=m;            else r=m;        }        if(l+1<5) printf("0\n");        else printf("%d\n",l+1);    }}