CF666E Forensic Examination（后缀自动机+动态线段树）

题意

　　给你一个串 $S$ 以及一个字符串数组 $T[1..m]$ ， $q$ 次询问，每次问 $S$ 的子串 $S[p_l..p_r]$ 在 $T[l..r]$ 中的哪个串里的出现次数最多，并输出出现次数。如有多解输出最靠前的那一个。

题解

　　神仙题……虽然的确是好题……然而码量好大……好麻烦……因为一个sb错误调了一个下午

　　先膜一波zsy大佬……%%%

　　首先先把$T$给建出一个广义$SAM$。考虑每一个询问的$p_r$，如果在$SAM$上匹配到了一个节点，那么这个子串$S[p_l,p_r]$一定是这个节点的一个祖先（然后如果根本匹配不到的话直接跳过）

　　然后可以考虑倍增找到祖先，只要找到满足$len_u>=p_r-p_l+1$的最上面的节点就好了

　　然后考虑如何表示这个节点在$[l,r]$范围内在哪些串出现了多少次，以及出现次数的最大值。区间查询可以考虑用线段树。对于每一个$SAM$上的节点，我们对他建一棵线段树，表示有这个节点代表的字符串在$T$中的出现次数。然后通过倍增找到节点，在线段树上查询就是了

　　于是直接把串$S$放到$SAM$上跑，然后在每一个节点记录一下有哪些询问在这个点上，一波$dfs$的时候顺便合并父亲和儿子的出现次数的区间，并求出答案

　　说的好像很简单……然而真的码得有点想哭……

 //minamoto
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 using namespace std;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while((ch=getchar())<''||ch>'')
     (ch=='-')&&(flag=true);
     for(res=num;(ch=getchar())>=''&&ch<='';res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 char sr[<<],z[];int C=-,Z;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(int x,char ch){
     if(C><<)Ot();if(x<)sr[++C]=,x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]=ch;
 }
 const int N=5e5+;
 struct data{
     int x,y;
     inline bool operator <(const data &b)const
     {return x<b.x||x==b.x&&y>b.y;}
 }ans[N];
 int L[N*],R[N*];data V[N*];
 struct node{int l,r,pl,pr;}q[N];
 int n,m,Q,last=,cnt=,ch[N][],fa[][N],l[N],rt[N],tot=;
 int nxt[][N],head[][N];
 char s[N],t[N];
 void ins(int c){
     int p=last,np=++cnt;last=np,l[np]=l[p]+;
     for(;p&&!ch[p][c];p=fa[][p]) ch[p][c]=np;
     if(!p) fa[][np]=;
     else{
         int q=ch[p][c];
         if(l[q]==l[p]+) fa[][np]=q;
         else{
             int nq=++cnt;l[nq]=l[p]+;
             memcpy(ch[nq],ch[q],sizeof(ch[q]));
             fa[][nq]=fa[][q],fa[][q]=fa[][np]=nq;
             for(;ch[p][c]==q;p=fa[][p]) ch[p][c]=nq;
         }
     }
 }
 void modify(int &now,int l,int r,int x){
     now=++tot;
     if(l==r) return (void)(++V[now].x,V[now].y=l);
     int mid=l+r>>;
     if(x<=mid) modify(L[now],l,mid,x);
     else modify(R[now],mid+,r,x);
     V[now]=max(V[L[now]],V[R[now]]);
 }
 void merge(int &x,int y){
     if(!x||!y) return (void)(x|=y);
     if(!L[x]&&!R[x]) return (void)(V[x].x+=V[y].x);
     merge(L[x],L[y]),merge(R[x],R[y]);
     V[x]=max(V[L[x]],V[R[x]]);
 }
 data query(int p,int l,int r,int ql,int qr){
     if(ql<=l&&qr>=r) return V[p];
     int mid=l+r>>;
     if(qr<=mid) return query(L[p],l,mid,ql,qr);
     if(ql>mid) return query(R[p],mid+,r,ql,qr);
     return max(query(L[p],l,mid,ql,qr),query(R[p],mid+,r,ql,qr));
 }
 void dfs(int u){
     for(int i=head[][u];i;i=nxt[][i])
     dfs(i),merge(rt[u],rt[i]);
     for(int i=head[][u];i;i=nxt[][i])
     ans[i]=query(rt[u],,m,q[i].l,q[i].r);
 }
 int main(){
     scanf("%s",s+);n=strlen(s+);
     m=read();
     for(int i=;i<=m;++i){
         scanf("%s",t+);int len=strlen(t+);last=;
         for(int j=;j<=len;++j) ins(t[j]-'a'),modify(rt[last],,m,i);
     }
     Q=read();
     for(int i=;i<=Q;++i){
         int a=read(),b=read(),c=read(),d=read();
         q[i]=(node){a,b,c,d};
         nxt[][i]=head[][q[i].pr],head[][q[i].pr]=i;
     }
     for(int i=;i<=cnt;++i)
     nxt[][i]=head[][fa[][i]],head[][fa[][i]]=i;
     for(int j=;j<;++j)
     for(int i=;i<=cnt;++i)
     fa[j][i]=fa[j-][fa[j-][i]];
     for(int i=,now=,cnt=;i<=n;++i){
         while(now&&!ch[now][s[i]-'a']) now=fa[][now],cnt=l[now];
         if(!now){now=,cnt=;continue;}
         now=ch[now][s[i]-'a'],++cnt;
         for(int j=head[][i];j;j=nxt[][j]){
             int u=now;if(cnt<q[j].pr-q[j].pl+) continue;
             for(int k=;~k;--k)
             if(l[fa[k][u]]>=q[j].pr-q[j].pl+)
             u=fa[k][u];
             nxt[][j]=head[][u],head[][u]=j;
         }
     }
     dfs();
     for(int i=;i<=Q;++i){
         if(!ans[i].x) ans[i].y=q[i].l;
         print(ans[i].y,),print(ans[i].x,);
     }
     Ot();
     return ;
 }