Codeforces 666E Forensic Examination SAM+权值线段树
Posted tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了Codeforces 666E Forensic Examination SAM+权值线段树相关的知识，希望对你有一定的参考价值。
第一次做这种$SAM$带权值线段树合并的题 然而$zjq$神犇看完题一顿狂码就做出来了 $Orz$
首先把所有串当成一个串建$SAM$ 我们对$SAM$上每个点 建一棵权值线段树 每个叶子节点表示一个匹配串能到达这个点的子串个数 这样我们对最后的$SAM$的权值线段树按$parent$树合并 询问的时候找到对应的$SAM$上的权值线段树直接查询就好了
具体的操作看代码吧~
#include<bits/stdc++.h>
using namespace std;
#define FO(x) {freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);}
#define pa pair<int,int>
#define mod 1000000007
#define ll long long
#define mk make_pair
#define pb push_back
#define fi first
#define se second
#define cl(x) memset(x,0,sizeof x)
#ifdef Devil_Gary
#define bug(x) cout<<(#x)<<" "<<(x)<<endl
#define debug(...) fprintf(stderr, __VA_ARGS__)
#else
#define bug(x)
#define debug(...)
#endif
const int INF = 0x7fffffff;
const int N=1230010;
/*
char *TT,*mo,but[(1<<15)+2];
#define getchar() ((TT==mo&&(mo=(TT=but)+fread(but,1,1<<15,stdin),TT==mo))?-1:*TT++)//*/
inline int read(){
int x=0,rev=0,ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')rev=1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return rev?-x:x;
}
pa operator * (const pa&A,const pa&B){return A.fi==B.fi?(A.se<B.se?A:B):(A.fi>B.fi?A:B);}
pa operator + (const pa&A,const pa&B){return mk(A.fi+B.fi,A.se);}
pa t[N*23];
char s[N];
int n,m,Q,c[N][27],par[N],len[N],fa[N][21],ls[N*23],rs[N*23],h[N],rt=1,sz=1,last=1,root[N],w[N],b[N],id;
void extend(int x){
int np=++sz,p=last;
len[np]=len[p]+1,last=np;
for(;p&&!c[p][x];p=par[p]) c[p][x]=np;
if(!p) par[np]=rt;
else{
int q=c[p][x];
if(len[p]+1==len[q]) par[np]=q;
else{
int nq=++sz;
len[nq]=len[p]+1,par[nq]=par[q],par[q]=par[np]=nq;
memcpy(c[nq],c[q],sizeof c[nq]);
for(;p&&c[p][x]==q;p=par[p]) c[p][x]=nq;
}
}
}
void modify(int&pos,int l,int r,int x){
if(!pos) pos=++id;
if(l==r){
t[pos].fi++,t[pos].se=l;
// cout<<t[pos].fi<<" "<<t[pos].se<<endl;
return;
}
int mid=l+r>>1;
if(x<=mid) modify(ls[pos],l,mid,x);
else modify(rs[pos],mid+1,r,x);
t[pos]=t[ls[pos]]*t[rs[pos]];
}
inline int merge(int x,int y,int l,int r)
{
if (!x||!y) return x|y;
int z=++id;
if (l==r) {
t[z]=t[x]+t[y];
return z;
}
int mid=(l+r)>>1;
ls[z]=merge(ls[x],ls[y],l,mid),rs[z]=merge(rs[x],rs[y],mid+1,r);
t[z]=t[ls[z]]*t[rs[z]];
return z;
}
pa Query(int pos,int l,int r,int nl,int nr){
if(!pos) return mk(0,0);
if(nl<=l&&r<=nr) return t[pos];
int mid=l+r>>1;
if(nr<=mid) return Query(ls[pos],l,mid,nl,nr);
else if(nl>mid) return Query(rs[pos],mid+1,r,nl,nr);
else return Query(ls[pos],l,mid,nl,nr)*Query(rs[pos],mid+1,r,nl,nr);
}
void topsort(){
for(int i=1;i<=sz;i++) w[len[i]]++;
for(int i=1;i<=sz;i++) w[i]+=w[i-1];
for(int i=1;i<=sz;i++) b[w[len[i]]--]=i;
for(int i=sz;i;i--){
int j=b[i];
if(par[j]) root[par[j]]=merge(root[par[j]],root[j],1,n);
}
}
void get_fa(){
for(int i=1;i<=sz;i++) fa[i][0]=par[i];
for(int i=1;i<=20;i++) for(int j=1;j<=sz;j++) fa[j][i]=fa[fa[j][i-1]][i-1];
}
int get(int l,int r){
int L=r-l+1,x=h[r];
for(int i=20;~i;i--) if(len[fa[x][i]]>=L) x=fa[x][i];
// -len[par[x]]-len[par[fa[x][i]]]
return x;
}
int main(){
#ifdef Devil_Gary
freopen("in.txt","r",stdin);
#endif
scanf("%s",s+1),m=strlen(s+1);
for(int i=1;i<=m;i++) extend(s[i]-'a'),h[i]=last;
n=read();
for(int i=1;i<=n;i++){
scanf("%s",s+1),m=strlen(s+1),extend(26);
for(int j=1;j<=m;j++) extend(s[j]-'a'),modify(root[last],1,n,i);
}
// bug(sz);
topsort(),get_fa();
for(Q=read();Q;Q--){
int l=read(),r=read(),pl=read(),pr=read(),x=get(pl,pr);
pa ans=Query(root[x],1,n,l,r);
if(!ans.fi) printf("%d 0\n",l);
else printf("%d %d\n",ans.se,ans.fi);
}
}
以上是关于Codeforces 666E Forensic Examination SAM+权值线段树的主要内容，如果未能解决你的问题，请参考以下文章
[CF 666E] Forensic Examination