poj 1743 Musical Theme (后缀数组+二分答案)

题意：

给出n个数字的串，将这个串后一个减去前一个，所得的串尽心操作。问所得的串最长不重复子串的长度。

题解：

这种题目很明显没办法直接求得长度，那么可以考虑二分答案判断是否满足条件。我们可以二分得到k，判断是否存在这样的k。对于如何判断，我们可以将后缀根据height的值进行分组，每组的后缀之间的最长公共前缀必须至少为k，这样之后只要找出每组中最大SA值和最小SA的值的差即可。好好理解下，这种做法，值得细细考虑。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=200005;
int rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],cnt[maxn];
int a[maxn];

void Swap(int*& x,int*& y)
{
    int *temp=x;
    x=y;
    y=temp;
}

int lcp(int t[],int a,int b,int l)
{
    return t[a]==t[b]&&t[a+l]==t[b+l];
}

void build_SA(int s[],int len,int up){

    int *k1=t1,*k2=t2,*r=t3;
    for(int i=0;i<up;i++)cnt[i]=0;
    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;

    for(int d=1,p=0;p<len;d<<=1){

        for(int i=len-d;i<len;i++)k2[p++]=i;
        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
        for(int i=0;i<len;i++)r[i]=k1[k2[i]];

        for(int i=0;i<up;i++)cnt[i]=0;
        for(int i=0;i<len;i++)cnt[r[i]]++;
        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];

        Swap(k1,k2);
        k1[SA[0]]=0;
        p=1;
        for(int i=1;i<len;i++){
                k1[SA[i]]= lcp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
        }
        if(p>=len)return ;///这句话切记要添加，忘了，检查了好几个小时。ORZ!
        up=p;
        p=0;
    }
}

void get_height(int s[],int len)
{
    for(int i=1;i<=len;i++)rank[SA[i]]=i;
    for(int i=0,p=0;i<len;i++){

        int j=SA[rank[i]-1];
        while(s[i+p]==s[j+p])p++;
        height[rank[i]]=p;
        if(p)p--;
    }
}

bool ok(int tag,int len)
{
    int i=1;
    int mx,mi;
    while(1)
    {
        while(i<len&&height[i]<tag)i++;
        if(i>=len)break;
        mx=SA[i-1];
        mi=SA[i-1];
        while(i<len&&height[i]>=tag){

            mx=max(mx,SA[i]);
            mi=min(mi,SA[i]);
            i++;
        }
        if(mx-mi>=tag)return true;
    }
    return false;
}

int solve(int len)
{
    int l=0,r=len,mid;
    int ans=0;
    while(l<=r){

        mid=(l+r)>>1;
        if(ok(mid,len)){
            l=mid+1;
            ans=mid;
        }else{
            r=mid-1;
        }
    }
    return ans;
}

int main(){

    int n;
    while(scanf("%d",&n)!=EOF){

        if(n==0)break;
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]);
        n--;
        for(int i=0;i<n;i++)
            a[i]=a[i+1]-a[i]+90;
        a[n]=0;
        build_SA(a,n+1,200);
        get_height(a,n);
        int ans=solve(n);
        if(ans>=4)ans++;
        else ans=0;
        printf("%d\n",ans);
    }
    return 0;
}//main
/**

*/