BZOJ4556 HEOI2016字符串

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

By：大奕哥

 #include<bits/stdc++.h>
 using namespace std;
 const int N=2e5+;
 int rt[N],num,n,q,aa,bb,cc,dd,cnt,head[N],fa[N][];
 char s[N];
 struct tree{
     int l,r;
 }t[N*];
 void change(int &x,int l,int r,int p)
 {
     if(!x)x=++num;
     if(l==r)return;
     int mid=l+r>>;
     if(p>mid)change(t[x].r,mid+,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 ;
     if(l==L&&r==R)return ;
     int mid=l+r>>;
     if(mid<L)return query(t[x].r,mid+,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+,r,mid+,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=;i<=;++i)fa[x][i]=fa[fa[x][i-]][i-];
     for(int i=head[x];i;i=e[i].nex)
     {
         int y=e[i].to;
         fa[y][]=x;dfs(y);
     }
 }
 struct SAM
 {
     int cnt,root,last,c[N][],f[N],dp[N],tt[N],pos[N],p[N],r[N],id[N],l[N];
     SAM(){cnt=;last=root=++cnt;}
     void add(int x,int dd)
     {
         int now=last,a=++cnt;last=a;id[a]=dd;
         l[a]=l[now]+;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]+)f[a]=q;
             else{
                 int b=++cnt;id[b]=dd;
                 l[b]=l[now]+;
                 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=;i<=cnt;++i)tt[l[i]]++;
         for(int i=;i<=n;++i)tt[i]+=tt[i-];
         for(int i=;i<=cnt;++i)r[tt[l[i]]--]=i;
         return;
     }
     bool judge(int x,int ll)
     {
         for(int i=;i>=;--i)
         if(l[fa[x][i]]>=ll)x=fa[x][i];
         return query(rt[x],,n,aa+ll-,bb);
     }
     void work()
     {
         for(int i=cnt;i>=;--i)
         {
             int x=r[i];
             change(rt[x],,n,id[x]);
             merge(rt[f[x]],rt[x]);
             addedge(f[x],x);
         }dfs();int ans=;
         for(int i=;i<=q;++i)
         {

             scanf("%d%d%d%d",&aa,&bb,&cc,&dd);
             int ll=,rr=min(bb-aa+,dd-cc+);ans=;
             aa=n-aa+;bb=n-bb+;cc=n-cc+;dd=n-dd+;
             swap(aa,bb);swap(cc,dd);int x=p[dd];
             while(ll<=rr)
             {
                 int mid=ll+rr>>;
                 if(judge(x,mid))ll=mid+,ans=mid;
                 else rr=mid-;
             }
             printf("%d\n",ans);
         }
         return;
     }
 }A;
 int main()
 {
     scanf("%d%d",&n,&q);
     scanf("%s",s+);
     reverse(s+,s++n);
     for(int i=;i<=n;++i)A.add(s[i]-'a',i);
     A.sort();A.work();
     return ;
 }