后缀数组 poj. 1743. Musial Theme

题意感觉有点难懂。但就是问你，给你一个数字串，然后问有没有相同子串（不能重合且长度大于五），或者不用相同的，全部加上一个或减去一个数之后相等的就好。
先是二分长度，然后通过height数组来判断是否可行

#include <iostream>#include <cstring>#include <cstdio>#include <cmath>#include <algorithm>#include<string>using namespace std;const int  maxn =4e4+10;int n;struct Suffix{    int s[maxn];    int SA[maxn],ranked[maxn],height[maxn];    int x[maxn],y[maxn],buc[maxn],len,m;    void init(int *str)    {        //len=strlen(str);        len=n;        for(int i=0;i<len;i++) s[i]=(int)str[i];        m=200;//    }    void build()    {        for (int i = 0; i < m; i++) buc[i] = 0; // buc 是一个桶        for (int i = 0; i < len; i++) buc[x[i] = s[i]]++;        for (int i = 1; i < m; i++)             buc[i] += buc[i - 1];        for (int i = len - 1; i >= 0; i--)            SA[--buc[x[i]]] = i;        for (int k = 1; k <= len; k <<= 1) { // k 倍增          int p = 0;          //对第二关键字排序，y[i]:第i大的第二关键字是谁          // 后缀 len - k 及之后的所有后缀第二关键字最小。为0          for (int i = len - 1; i >= len - k; i--)            y[p++] = i;          for (int i = 0; i < len; i++)                if (SA[i] >= k)                    y[p++] = SA[i] - k;         //总体来排个序，求出SA          for (int i = 0; i < m; i++) buc[i] = 0;          for (int i = 0; i < len; i++)                buc[x[y[i]]]++;          for (int i = 1; i < m; i++)                buc[i] += buc[i - 1];          for (int i = len - 1; i >= 0; i--)                SA[--buc[x[y[i]]]] = y[i];          swap(x, y);          p = 1; x[SA[0]] = 0;          // 重新计算每个一元的名次。则x数组里存的是总体的顺序          for (int i = 1; i < len; i++) {            if (y[SA[i - 1]] == y[SA[i]] && y[SA[i - 1] + k] == y[SA[i] + k])              x[SA[i]] = p - 1;            else x[SA[i]] = p++;          }          if (p >= len) break; // 每个后缀的名次已经完全不同，不需要继续倍增              m = p; // 更新名次的最大值。        }    }    void LCP()    {        for(int i=0;i<len;i++)            ranked[SA[i]]=i;        int k=0;        for(int i=0;i<len;i++)        {            if(ranked[i]==0) {height[0]=0;continue;}            if(k) k--;            int j=SA[ranked[i]-1];            while(s[i+k]==s[j+k]&&i+k<len&&j+k<len) k++;            height[ranked[i]]=k;        }    } /*  int LCS(char *s1)    {        int l=strlen(s1);        s1[l]=' ';        scanf("%s",s1+l+1);        init(s1);        build();        LCP();        int ans=0;        for(int i=1;i<len;i++)        if((SA[i-1]<l&&SA[i]>l)||(SA[i-1]>l&&SA[i]<l))                ans=max(ans,height[i]);        return ans;    }*/};Suffix a;int ok(int mid){    int i=1;    while(1){    for(;i<n;i++)        if(a.height[i]>=mid)            break;    if(i>=n) break;    int minn=a.SA[i-1];    int maxx=a.SA[i-1];    while(i<n&&a.height[i]>=mid)    {        minn=min(a.SA[i],minn);        maxx=max(a.SA[i],maxx);        i++;    }    if(maxx-minn>=mid) return 1;    }    return 0;}int num[maxn];int str[maxn];int main(){    while(scanf("%d",&n)!=EOF&&n)    {       memset(str,0,sizeof(str));      memset(num,0,sizeof(num));        for(int i=0;i<n;i++)            scanf("%d",&num[i]);        if(n<10) printf("0\n");        else        {            for(int i=1;i<n;i++)                str[i-1]=num[i]-num[i-1]+100;             n--;         //   str[n-1]=0;            a.init(str);            a.build();            a.LCP();            int l=4,r=n/2,res=0;            while(l<=r)            {                int mid=l+(r-l)/2;                if(ok(mid)){                    res=mid;                    l=mid+1;                }                else r=mid-1;            }            if(res>=4)            printf("%d\n",res+1);            else printf("0\n");        }    }    return 0;}