bzoj 4556: [Tjoi2016&Heoi2016]字符串

题记

讲真…这题看了好久啊好久啊好久..久到我都..不想说什么了
嗯这题..正解..还..有待..学习【捂脸】
（挖坑）

---

题解

bzoj4556
后缀数组..+暴力【捂脸】
这个暴力啊..就..从c开始的后缀的排名..往前、后查
因为..i开始的后缀与c开始的后缀的lcp=i与c之间height的最小值
因此..就..依次往前、后推好啦，如果..此时的height值小于答案，就退出循环

---

代码

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 100010
char str[N];
int n,m,cnt[N],rank[2*N],rank1[2*N],tmp[N],sa[N],height[N];
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while('0'<=ch && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline void work(int a,int b,int c,int d){
    int ans=0,len=n;
    if(a<=c && c<=b) ans=min(b-c+1,d-c+1);
    for(int i=rank[c];i>=2;i--){
        len=min(len,height[i]); 
        if(len<=ans) break;
        if(a<=sa[i-1] && sa[i-1]<=b) ans=max(ans,min(len,b-sa[i-1]+1));
    }len=n;
    for(int i=rank[c]+1;i<=n;i++){
        len=min(len,height[i]); 
        if(len<=ans) break;
        if(a<=sa[i] && sa[i]<=b) ans=max(ans,min(len,b-sa[i]+1));
    }print(min(d-c+1,ans));
}
int main(){
    n=read();m=read();
    scanf("%s",str+1);
    memset(cnt,0,sizeof(cnt));
    for(int i=1;i<=n;i++) cnt[str[i]]=1;
    for(int i=1;i<=260;i++) cnt[i]+=cnt[i-1];
    for(int i=1;i<=n;i++) rank[i]=cnt[str[i]];
    int k=1;
    for(int len=1;k!=n;len*=2){
        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=n;i++) cnt[rank[i+len]]++;
        for(int i=1;i<=n;i++) cnt[i]+=cnt[i-1];
        for(int i=n;i>=1;i--) tmp[cnt[rank[i+len]]--]=i;
        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=n;i++) cnt[rank[tmp[i]]]++;
        for(int i=1;i<=n;i++) cnt[i]+=cnt[i-1];
        for(int i=n;i>=1;i--) sa[cnt[rank[tmp[i]]]--]=tmp[i];
        memcpy(rank1,rank,sizeof(rank));
        k=1;rank[sa[1]]=1;
        for(int i=2;i<=n;i++){
            if(rank1[sa[i]]!=rank1[sa[i-1]] || rank1[sa[i]+len]!=rank1[sa[i-1]+len]) k++;
            rank[sa[i]]=k;
        }
    }k=0;
    for(int i=1;i<=n;i++){
        if(rank[i]==1) k=0;
        else{
            if(k>0) k--;
            while(str[i+k]==str[sa[rank[i]-1]+k]) k++;
        }height[rank[i]]=k;
    }
    for(int i=1;i<=m;i++){
        int a=read(),b=read(),c=read(),d=read();
        work(a,b,c,d); 
    }
    return 0;
}