【CF666E】Forensic Examination（后缀自动机，线段树合并）
【CF666E】Forensic Examination（后缀自动机，线段树合并）
题面
翻译：
给定一个串\(S\)和若干个串\(T_i\)
每次询问\(S[pl..pr]\)在\(T_l..T_r\)中出现的最多次数，以及出现次数最多的那个串的编号。
题解
好题啊。
我们对于所有的\(T\)串构建出广义后缀自动机之后
把\(S\)丢到\(SAM\)上匹配，对于每组询问，
相当于回答\(S[pl..pr]\)所代表的节点的\(right\)集合所代表的串的众数是哪个串，以及这个众数出现的次数。
考虑如何维护\(right\)集合关于每个\(T\)的出现次数
我们对于每个节点开一棵线段树，线段树的值域是\(1..m\)，用来维护每个\(T\)串的出现次数。
每个节点维护的值就是\(T\)串的最多出现次数。
这样就挺好做了，\(right\)集合只需要线段树合并就可以快速求得。
对于每个询问挂链，然后\(dfs\)一遍\(parent\)树的时候，顺便在线段树上查询一下区间最大就行了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 1111111
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
struct SAM
{
struct Node
{
int son[26];
int ff,len;
}t[MAX];
int last,tot;
void init(){last=tot=1;}
void extend(int c)
{
int p=last,np=++tot;last=np;
t[np].len=t[p].len+1;
while(p&&!t[p].son[c])t[p].son[c]=np,p=t[p].ff;
if(!p)t[np].ff=1;
else
{
int q=t[p].son[c];
if(t[q].len==t[p].len+1)t[np].ff=q;
else
{
int nq=++tot;
t[nq]=t[q];
t[nq].len=t[p].len+1;
t[q].ff=t[np].ff=nq;
while(p&&t[p].son[c]==q)t[p].son[c]=nq,p=t[p].ff;
}
}
}
}SAM;
int n,m,f[22][MAX];
struct Data{int v,p;}ans[MAX];
bool operator<(Data a,Data b){return (a.v<b.v)||(a.v==b.v&&a.p>b.p);}
struct SegmentTree{int ls,rs;Data v;}t[MAX<<4];
int Tot,rt[MAX];
struct query{int l,r,pl,pr;}q[MAX];
char S[MAX],T[MAX];
void Modify(int &x,int l,int r,int p)
{
if(!x)x=++Tot;
if(l==r){t[x].v.v++;t[x].v.p=p;return;}
int mid=(l+r)>>1;
if(p<=mid)Modify(t[x].ls,l,mid,p);
else Modify(t[x].rs,mid+1,r,p);
t[x].v=max(t[t[x].ls].v,t[t[x].rs].v);
}
void Merge(int &x,int y)
{
if(!x||!y){x|=y;return;}
if(!t[x].ls&&!t[x].rs){t[x].v.v+=t[y].v.v;return;}
Merge(t[x].ls,t[y].ls);Merge(t[x].rs,t[y].rs);
t[x].v=max(t[t[x].ls].v,t[t[x].rs].v);
}
Data Query(int x,int l,int r,int L,int R)
{
if(L==l&&r==R)return t[x].v;
int mid=(l+r)>>1;
if(R<=mid)return Query(t[x].ls,l,mid,L,R);
if(L>mid)return Query(t[x].rs,mid+1,r,L,R);
return max(Query(t[x].ls,l,mid,L,mid),Query(t[x].rs,mid+1,r,mid+1,R));
}
struct Link
{
struct Line{int v,next;}e[MAX];
int h[MAX],cnt;
void Add(int u,int v){e[++cnt]=(Line){v,h[u]};h[u]=cnt;}
}Par,Qy,Aw;
void dfs(int u)
{
for(int i=Par.h[u];i;i=Par.e[i].next)
dfs(Par.e[i].v),Merge(rt[u],rt[Par.e[i].v]);
for(int i=Aw.h[u];i;i=Aw.e[i].next)
ans[Aw.e[i].v]=Query(rt[u],1,m,q[Aw.e[i].v].l,q[Aw.e[i].v].r);
}
int main()
{
scanf("%s",S+1);n=strlen(S+1);
m=read();SAM.init();
for(int i=1;i<=m;++i)
{
SAM.last=1;scanf("%s",T+1);
for(int j=1,l=strlen(T+1);j<=l;++j)
SAM.extend(T[j]-97),Modify(rt[SAM.last],1,m,i);
}
int Q=read();
for(int i=1;i<=Q;++i)
{
q[i]=(query){read(),read(),read(),read()};
Qy.Add(q[i].pr,i);
}
for(int i=2;i<=SAM.tot;++i)Par.Add(f[0][i]=SAM.t[i].ff,i);
for(int i=1;i<22;++i)
for(int j=1;j<=SAM.tot;++j)f[i][j]=f[i-1][f[i-1][j]];
for(int i=1,nw=1,len=0;i<=n;++i)
{
int c=S[i]-97;
while(nw&&!SAM.t[nw].son[c])nw=SAM.t[nw].ff,len=SAM.t[nw].len;
if(!nw){nw=1,len=0;continue;}
nw=SAM.t[nw].son[c];len+=1;
for(int j=Qy.h[i];j;j=Qy.e[j].next)
{
int v=Qy.e[j].v,u=nw;
if(len<q[v].pr-q[v].pl+1)continue;
for(int k=21;~k;--k)
if(SAM.t[f[k][u]].len>=q[v].pr-q[v].pl+1)u=f[k][u];
Aw.Add(u,v);
}
}
dfs(1);
for(int i=1;i<=Q;++i)
{
if(!ans[i].v)ans[i].p=q[i].l;
printf("%d %d\n",ans[i].p,ans[i].v);
}
return 0;
}
