POJ 1743 Musical Theme

　　感觉最近好混乱......各种OJ都刷一点，感觉不太好......尤其是这种英文题

　　这道题一开始还没有看懂。听了ljh大犇的解释后终于明白了。下面我为英语和我一样的人翻译一下题面：

　　输入n个数。求最长的两端不交叉的序列，使它们的变化量相同，即相邻两位的差相同。当这个长度小于5时输出0。

　　好了，这道题是不是和最长不重叠子串有点像？把原数组差分后，简直一模一样。我们先构出原数组的后缀数组，然后二分一个答案（答案显然具有单调性是吧？），扫一遍height数组看一看是否存在后缀i，j，使得lcp(i,j)>=x && abs(len(j)-len(i))>x 。找到了当前答案（x）就是可行的。一路判断下去就可以辣！不过这道题要注意：原数组差分之后只有n-1个！

　　下面贴代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout);
#define maxn 200010

using namespace std;
typedef long long llg;

int n,ht[maxn],rk[maxn],sa[maxn],c[maxn],a[maxn];

int getint(){
    int w=0;bool q=0;
    char c=getchar();
    while((c>'9'||c<'0')&&c!='-') c=getchar();
    if(c=='-') q=1,c=getchar();
    while(c>='0'&&c<='9') w=w*10+c-'0',c=getchar();
    return q?-w:w;
}

void build_sa(int m){
    int i,*x=rk,*y=ht;
    for(i=1;i<=m;i++) c[i]=0;
    for(i=1;i<=n;i++) c[x[i]=a[i]]++;
    for(i=1;i<=m;i++) c[i]+=c[i-1];
    for(i=n;i;i--) sa[c[x[i]]--]=i;
    for(int k=1,p;k<=n;k<<=1){
        p=0;
        for(i=n-k+1;i<=n;i++) y[++p]=i;
        for(i=1;i<=n;i++) if(sa[i]>k) y[++p]=sa[i]-k;
        for(i=1;i<=m;i++) c[i]=0;
        for(i=1;i<=n;i++) c[x[y[i]]]++;
        for(i=1;i<=m;i++) c[i]+=c[i-1];
        for(i=n;i;i--) sa[c[x[y[i]]]--]=y[i];
        swap(x,y); x[sa[1]]=1; p=1;
        for(i=2;i<=n;i++)
            x[sa[i]]=(y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k])?p:++p;
        if(p==n) break; m=p;
    }
    for(i=1;i<=n;i++) rk[sa[i]]=i;
    for(int i=1,j,k=0;i<=n;i++){
        if(k) k--;
        j=sa[rk[i]-1];
        while(a[i+k]==a[j+k]) k++;
        ht[rk[i]]=k;
    }
}

bool pd(int x){
    for(int i=2,ma,mi;i<=n;){
        while(ht[i]<x && i<=n) i++;
        ma=mi=sa[i-1];
        while(ht[i]>=x && i<=n){
            ma=max(ma,sa[i]);
            mi=min(mi,sa[i]);
            i++;
        }
        if(ma-mi>x) return 1;
    }
    return 0;
}

int main(){
    File("a");
    n=getint();
    while(n){
        for(int i=0;i<n;i++) a[i]=getint();
        for(int i=n-1;i;i--) a[i]=a[i]-a[i-1]+100;
        build_sa(200);
        int l=0,r=n+1,mid;
        while(l!=r){
            mid=l+r>>1;
            if(pd(mid)) l=mid+1;
            else r=mid;
        }
        if(l<=4) printf("0
");
        else printf("%d
",l);
        n=getint();
    }
    return 0;
}