BZOJ4556 HEOI2016字符串

没错，又是这题，使用后缀自动机，反向建树，主席树维护right集合。

By：大奕哥

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
int rt[N],num,n,q,aa,bb,cc,dd,cnt,head[N],fa[N][19];
char s[N];
struct tree{
    int l,r;
}t[N*50];
void change(int &x,int l,int r,int p)
{
    if(!x)x=++num;
    if(l==r)return;
    int mid=l+r>>1;
    if(p>mid)change(t[x].r,mid+1,r,p);
    else change(t[x].l,l,mid,p);
}
void merge(int &x,int y)
{
    if(!x||!y){x=x+y;return;}
    ++num;t[num]=t[x];x=num;
    merge(t[x].l,t[y].l);
    merge(t[x].r,t[y].r);
}
bool query(int x,int l,int r,int L,int R)
{
    if(!x)return 0;
    if(l==L&&r==R)return 1;
    int mid=l+r>>1;
    if(mid<L)return query(t[x].r,mid+1,r,L,R);
    else if(mid>R)return query(t[x].l,l,mid,L,R);
    else return query(t[x].l,l,mid,L,mid)|query(t[x].r,mid+1,r,mid+1,R); 
}
struct edge{
    int to,nex;
}e[N];
void addedge(int x,int y)
{
    e[++cnt].to=y;e[cnt].nex=head[x];head[x]=cnt;
}
void dfs(int x)
{
    for(int i=1;i<=18;++i)fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int i=head[x];i;i=e[i].nex)
    {
        int y=e[i].to;
        fa[y][0]=x;dfs(y);
    }
}
struct SAM
{
    int cnt,root,last,c[N][26],f[N],dp[N],tt[N],pos[N],p[N],r[N],id[N],l[N];
    SAM(){cnt=0;last=root=++cnt;}
    void add(int x,int dd)
    {
        int now=last,a=++cnt;last=a;id[a]=dd;
        l[a]=l[now]+1;p[dd]=a;
        for(;now&&!c[now][x];now=f[now])c[now][x]=a;
        if(!now)f[a]=root;
        else{
            int q=c[now][x];
            if(l[q]==l[now]+1)f[a]=q;
            else{
                int b=++cnt;id[b]=dd;
                l[b]=l[now]+1;
                f[b]=f[q];
                f[a]=f[q]=b;
                memcpy(c[b],c[q],sizeof(c[q]));
                for(;now&&c[now][x]==q;now=f[now])c[now][x]=b;
            }
        }
        return;
    }
    void sort(){
        for(int i=1;i<=cnt;++i)tt[l[i]]++;
        for(int i=1;i<=n;++i)tt[i]+=tt[i-1];
        for(int i=1;i<=cnt;++i)r[tt[l[i]]--]=i;
        return;
    }
    bool judge(int x,int ll)
    {
        for(int i=18;i>=0;--i)
        if(l[fa[x][i]]>=ll)x=fa[x][i];
        return query(rt[x],1,n,aa+ll-1,bb);
    }
    void work()
    {
        for(int i=cnt;i>=2;--i)
        {
            int x=r[i];
            change(rt[x],1,n,id[x]);
            merge(rt[f[x]],rt[x]);
            addedge(f[x],x);
        }dfs(1);int ans=0;
        for(int i=1;i<=q;++i)
        {
            
            scanf("%d%d%d%d",&aa,&bb,&cc,&dd);
            int ll=1,rr=min(bb-aa+1,dd-cc+1);ans=0;
            aa=n-aa+1;bb=n-bb+1;cc=n-cc+1;dd=n-dd+1;
            swap(aa,bb);swap(cc,dd);int x=p[dd];
            while(ll<=rr)
            {
                int mid=ll+rr>>1;
                if(judge(x,mid))ll=mid+1,ans=mid;
                else rr=mid-1;
            }
            printf("%d
",ans);
        }
        return;
    }
}A;
int main()
{
    scanf("%d%d",&n,&q);
    scanf("%s",s+1);
    reverse(s+1,s+1+n);
    for(int i=1;i<=n;++i)A.add(s[i]-'a',i);
    A.sort();A.work();
    return 0;
}