poj1743--Musical Theme（后缀数组）

题意：求一列数字中走向相同的两个字序列，长度要求大于5

题解：相邻数字求差，原题就变成求相同的长度大于4的子串。

［存疑：在保证两个子串不相交时觉得限定条件应该是大于x，但是wa了＝ ＝ 不是很理解］

/**************************************

Problem: 1743        User: G_lory

Memory: 1392K        Time: 204MS

Language: G++        Result: Accepted

**************************************/

//后缀数组

#include <stdio.h>

#include <cstring>

#include <iostream>

#include <algorithm>

using namespace std;

const int N = int(2e5)+;

int cmp(int *r,int a,int b,int l){

    return (r[a]==r[b]) && (r[a+l]==r[b+l]);

}

// 用于比较第一关键字与第二关键字,

// 比较特殊的地方是,预处理的时候,r[n]=0(小于前面出现过的字符)

int wa[N],wb[N],wss[N],wv[N];

int sa[N];        // 排第几的是谁 0~n-1

int rk[N],     // 谁排第几

    height[N];    // 排名相邻的两个后缀的最长公共前缀长度:suffix(sa[i-1])和(sa[i]) 的最长公共前缀，

int a[N];

void DA(int *r,int *sa,int n,int m){                    // 此处N比输入的N要多1，为人工添加的一个字符，用于避免CMP时越界

    int i,j,p,*x=wa,*y=wb,*t;

    for(i=;i<m;i++) wss[i]=;

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

    for(i=;i<m;i++) wss[i]+=wss[i-];

    for(i=n-;i>=;i--) sa[--wss[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++) wss[i]=;

        for(i=;i<n;i++) wss[wv[i]]++;

        for(i=;i<m;i++) wss[i]+=wss[i-];

        for(i=n-;i>=;i--) sa[--wss[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++;

    }

}

void calheight(int *r,int *sa,int n){                    // 此处N为实际长度

    int i,j,k=;

    for(i=;i<=n;i++) rk[sa[i]]=i;

    for(i=;i<n; height[rk[i++]] = k )

    for(k?k--:,j=sa[rk[i]-]; r[i+k]==r[j+k]; k++);

}

bool ok(int x, int n)

{

    int minn, maxn;

    minn = maxn = sa[];

    for (int i = ; i <= n; ++i)

    {

        if (height[i] >= x) {

            minn = min(minn, sa[i]);

            maxn = max(maxn, sa[i]);

        } else {

            minn = maxn = sa[i];

        }

        if (maxn - minn >= x) return true;

    }

    return false;

}

int main(int argc, char const *argv[])

{

    //freopen("in", "r", stdin);

    int n;

    while (cin >> n && n) {

        for (int i = ; i < n; ++i) scanf("%d", a+i);

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

        a[--n] = ;

        DA(a, sa, n+, );

        calheight(a, sa, n);

        int l = , r = n;

        while (l +  < r) {

            int mid = (l+r) >> ;

            if (ok(mid, n)) l = mid;

            else r = mid;

        }

        ++l;

        printf("%d\n", l <  ?  : l);

    }

    return ;

}