bzoj4556 [Tjoi2016&Heoi2016]字符串

http://www.elijahqi.win/archives/2947
Description

佳媛姐姐过生日的时候，她的小伙伴从某东上买了一个生日礼物。生日礼物放在一个神奇的箱子中。箱子外边写了
一个长为n的字符串s，和m个问题。佳媛姐姐必须正确回答这m个问题，才能打开箱子拿到礼物，升职加薪，出任CE
O，嫁给高富帅，走上人生巅峰。每个问题均有a,b,c,d四个参数，问你子串s[a..b]的所有子串和s[c..d]的最长公
共前缀的长度的最大值是多少？佳媛姐姐并不擅长做这样的问题，所以她向你求助，你该如何帮助她呢？
Input

输入的第一行有两个正整数n,m，分别表示字符串的长度和询问的个数。接下来一行是一个长为n的字符串。接下来
m行，每行有4个数a,b,c,d，表示询问s[a..b]的所有子串和s[c..d]的最长公共前缀的最大值。1<=n,m<=100,000,
字符串中仅有小写英文字母，a<=b,c<=d,1<=a,b,c,d<=n

Output

对于每一次询问，输出答案。

Sample Input

5 5
aaaaa
1 1 1 5
1 5 1 1
2 3 2 3
2 4 2 3
2 3 2 4
Sample Output

1
1
2
2
2
看了一会题目 发现不会做 看了眼标签 好大概yy了一个sa+主席树的做法 主要是二分答案没想到
首先做SA 然后针对每个rk数组建立主席树可以想象为横坐标是下标纵坐标是rk里面的数
那么每次做的时候相当于先二分一个答案 然后求一下c串在sa数组中往前往后能延伸多远求出这个L,R之后在主席树rt[a-1],rt[b-len+1]之间判断下L,R区间是否存在后缀即可

#include<cstdio>
#include<cctype>
#include<cstring> 
#include<algorithm>
#define N 100010
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int out_len=1<<16;
char obuf[out_len],*oh=obuf;
inline void write_char(char c){
    if (oh==obuf+out_len) fwrite(obuf,1,out_len,stdout),oh=obuf;
    *oh++=c;
}
template<class T>
inline void W(T x){
    static int buf[30],cnt;
    if (!x) write_char('0');else{
        if (x<0) write_char('-'),x=-x;
        for (cnt=0;x;x/=10) buf[++cnt]=x%10+48;
        while(cnt) write_char(buf[cnt--]);
    }
}
inline void flush(){fwrite(obuf,1,oh-obuf,stdout);}
char s[N];int mn[N][20];
int n,q,rk[N<<1],rk1[N<<1],m,height[N],cnt[N],tmp[N],sa[N],k,rt[N],num,Log[N];
struct node{
    int left,right,s;
}tree[N*20];
inline int lcp(int l,int r){
    ++l;int t=Log[r-l+1];
    return min(mn[l][t],mn[r-(1<<t)+1][t]);
}
inline void insert1(int &x,int l,int r,int p){
    tree[++num]=tree[x];x=num;++tree[x].s;
    if(l==r) return;int mid=l+r>>1;
    if (p<=mid) insert1(tree[x].left,l,mid,p);else insert1(tree[x].right,mid+1,r,p);
}int L,R;
inline int query(int rt1,int rt2,int l,int r,int l1,int r1){
    if (l1<=l&&r1>=r) return tree[rt2].s-tree[rt1].s;int mid=l+r>>1,tmp=0;
    if (l1<=mid) tmp=query(tree[rt1].left,tree[rt2].left,l,mid,l1,r1);if(tmp) return tmp;
    if (r1>mid) tmp=query(tree[rt1].right,tree[rt2].right,mid+1,r,l1,r1);return tmp;
}
inline void find(int x,int len){
    L=x;R=x;static int l,r,mid;
    if (x>1){
        l=1,r=x-1;
        while(l<=r){
            mid=l+r>>1;
            if (lcp(mid,x)>=len) r=mid-1;
            else l=mid+1;
        }L=l;
    }
    if (x<n){
        l=x+1,r=n;
        while(l<=r){
            mid=l+r>>1;
            if (lcp(x,mid)>=len) l=mid+1;else r=mid-1;
        }R=r;
    }
}
inline bool check(int a,int b,int c,int md){
    find(rk[c],md);return query(rt[a-1],rt[b-md+1],1,n,L,R);
}
int main(){
    freopen("bzoj4556.in","r",stdin);
    n=read();q=read();char ch=gc();while(ch<'a'||ch>'z') ch=gc();Log[0]=-1;
    for (int i=1;i<=n;++i) s[i]=ch,ch=gc(),Log[i]=Log[i>>1]+1;
    for (int i=1;i<=n;++i) cnt[s[i]]=1;
    for (int i=1;i<=255;++i) cnt[i]+=cnt[i-1];
    for (int i=1;i<=n;++i) rk[i]=cnt[s[i]];m=30;
    for (int p=1;k!=n;p<<=1,m=k){static int i;
        for (int i=1;i<=m;++i) cnt[i]=0;
        for (i=1;i<=n;++i) ++cnt[rk[p+i]];
        for (i=1;i<=m;++i) cnt[i]+=cnt[i-1];
        for (i=n;i;--i) tmp[cnt[rk[i+p]]--]=i;
        for (int i=1;i<=m;++i) cnt[i]=0;
        for (i=1;i<=n;++i) ++cnt[rk[i]];
        for (i=1;i<=m;++i) cnt[i]+=cnt[i-1];
        for (i=n;i;--i) sa[cnt[rk[tmp[i]]]--]=tmp[i]; 
        memcpy(rk1,rk,sizeof(rk)>>1);k=rk[sa[1]]=1;
        for (i=2;i<=n;++i){
            if (rk1[sa[i]]!=rk1[sa[i-1]]||rk1[sa[i]+p]!=rk1[sa[i-1]+p]) ++k;
            rk[sa[i]]=k;
        }
    }k=0;
    for (int i=1;i<=n;++i){
        if(rk[i]==1) continue;k=k==0?0:k-1;
        while(s[i+k]==s[sa[rk[i]-1]+k]) ++k;
        height[rk[i]]=k;
    }
    for (int i=1;i<=n;++i) mn[i][0]=height[i];
    for (int j=1;j<=Log[n];++j)
        for (int i=1;i+(1<<j)<=n+1;++i)
            mn[i][j]=min(mn[i][j-1],mn[i+(1<<j-1)][j-1]);
    for (int i=1;i<=n;++i) rt[i]=rt[i-1],insert1(rt[i],1,n,rk[i]);
    for (int i=1;i<=q;++i){static int a,b,c,d,l,r;
        a=read();b=read();c=read();d=read();l=1;r=min(b-a+1,d-c+1);
        while(l<=r){
            static int mid;mid=l+r>>1;
            if (check(a,b,c,mid)) l=mid+1;else r=mid-1;
        }
        W(r);write_char('\n');
    }flush();
    return 0;
}