4556: [Tjoi2016&Heoi2016]字符串 后缀数组+主席树+二分

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

题解：

先搞出后缀数组的一堆东西，二分答案，转为判定性问题。然后二分出答案可以取的rank区间，这个可以用RMQ$O(1)$询问区间最小值辅助完成，然后在$[a,b-mid+1]$区间内找是否有符合条件的rank，这个可以用主席树完成，要注意二分时的各种上下界、细节。
什么？你说这不是后缀数组？反正我过了

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<LL,int>
const int Maxn=100010;
const int inf=2147483647;
const unsigned int base=233;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*f;
}
int sa[Maxn],len,q,height[Maxn],rank[Maxn],mn[Maxn][17],Log[Maxn];
unsigned int Hash[Maxn],Pow[Maxn];
char s[Maxn];
bool same(int x,int y,int l)
{
    if(!l)return true;
    if(x<y)swap(x,y);
    return((Hash[y]-Hash[y+l])==(Hash[x]-Hash[x+l])*Pow[x-y]);
}
bool cmp(int x,int y)
{
    int mx=len-max(x,y)+1;
    if(same(x,y,mx))return x>y;
    int l=1,r=mx-1;
    while(l<=r)
    {
        int mid=l+r>>1;
        if(same(x,y,mid))l=mid+1;
        else r=mid-1;
    }
    return s[x+l-1]<s[y+l-1];
}
void get_height()
{
    for(int i=1;i<=len;i++)rank[sa[i]]=i;
    for(int i=1;i<=len;i++)
    {
        int j=sa[rank[i]-1],k=height[rank[i-1]];
        if(k)k--;
        while(s[i+k]==s[j+k])k++;
        height[rank[i]]=k;
    }
}
int query_min(int l,int r)
{
    if(l>r)return inf;
    int t=Log[r-l+1];
    return min(mn[l][t],mn[r-(1<<t)+1][t]);
}
int root[Maxn],lc[Maxn*20],rc[Maxn*20],c[Maxn*20],tot=0;
void insert(int &u,int l,int r,int p)
{
    if(!u)u=++tot;
    c[u]++;
    if(l==r)return;
    int mid=l+r>>1;
    if(p<=mid)insert(lc[u],l,mid,p);
    else insert(rc[u],mid+1,r,p);
}
void merge(int &u1,int u2)
{
    if(!u1){u1=u2;return;}
    if(!u2)return;
    c[u1]+=c[u2];
    merge(lc[u1],lc[u2]);
    merge(rc[u1],rc[u2]);
}
int query(int r1,int r2,int l,int r,int ql,int qr)
{
    if(c[r1]==c[r2])return 0;
    if(ql==l&&qr==r)return c[r1]-c[r2];
    int mid=l+r>>1;
    if(qr<=mid)return query(lc[r1],lc[r2],l,mid,ql,qr);
    else if(ql>mid)return query(rc[r1],rc[r2],mid+1,r,ql,qr);
    else
    {
        if(query(lc[r1],lc[r2],l,mid,ql,mid))return 1;
        return query(rc[r1],rc[r2],mid+1,r,mid+1,qr);
    }
}
int main()
{
    len=read(),q=read();
    scanf("%s",s+1);
    Log[1]=0;
    for(int i=2;i<=len;i++)Log[i]=Log[i>>1]+1;
    for(int i=1;i<=len;i++)sa[i]=i;
    Pow[0]=1;
    for(int i=1;i<=len;i++)Pow[i]=Pow[i-1]*base;
    Hash[len+1]=0;
    for(int i=len;i;i--)Hash[i]=Hash[i+1]+Pow[len-i+1]*s[i];
    sort(sa+1,sa+1+len,cmp);
    get_height();
    for(int i=1;i<=len;i++)mn[i][0]=height[i];
    for(int i=1;(1<<i)<=len;i++)
    for(int j=1;j+(1<<i)-1<=len;j++)
    mn[j][i]=min(mn[j][i-1],mn[j+(1<<(i-1))][i-1]);
    root[0]=0;
    for(int i=1;i<=len;i++)insert(root[i],1,len,rank[i]),merge(root[i],root[i-1]);
    while(q--)
    {
        int a=read(),b=read(),c=read(),d=read();
        int l=0,r=min(b-a+1,d-c+1);
        while(l<=r)
        {
            int mid=l+r>>1;
            if(!mid)
            {
                l=mid+1;
                continue;
            }
            int ll,rr,L,R;
            /*二分左端点*/
            ll=0,rr=rank[c];
            while(ll<=rr)
            {
                int mmid=ll+rr>>1;
                if(query_min(mmid+1,rank[c])>=mid)rr=mmid-1;
                else ll=mmid+1;
            }
            L=rr+1;
            /*二分右端点*/
            ll=rank[c]+1,rr=len;
            while(ll<=rr)
            {
                int mmid=ll+rr>>1;
                if(query_min(L+1,mmid)>=mid)ll=mmid+1;
                else rr=mmid-1;
            }
            R=ll-1;
            if(L<=R)
            {
                if(query(root[b-mid+1],root[a-1],1,len,L,R))l=mid+1;
                else r=mid-1;
            }
            else r=mid-1;
        }
        printf("%d\n",l-1);
    }
}