【学习笔记】CF1037H Security

仰望星空的蚂蚁

于 2023-03-18 10:56:12 发布

阅读量97

点赞数

分类专栏： codeforces 文章标签：深度优先图论算法

本文链接：https://blog.csdn.net/cqbzlydd/article/details/129508748

版权

codeforces 专栏收录该内容

6 篇文章 0 订阅

订阅专栏

本文介绍了一种基于SAM（SuffixAutomaton）的数据结构解决字符串相关问题的方法，特别是关于字典序和最长公共前后缀（LCS）的查询。文章提到，在固定右边界的情况下，问题可以通过离线处理简化，避免了使用线段树合并的复杂性。代码示例展示了如何构建和操作SAM来高效地处理这些查询。

摘要由CSDN通过智能技术生成

非常基础的 $S A M$ 题目。

假设给你了 $S_1$ ，那么显然非常好求，随便怎么搞都行。

字典序这玩意也非常好翻译啊，直接枚举前缀 $L CS$ 就完了。

然后就是实现细节。乍一看这玩意要线段树合并，还要可持久化啊 ~~直接把我劝退了~~ 。不过如果固定了 $r$ 就非常好做，因此直接一波离线即可。

~~很久没打这种数据结构了~~

一遍过了，就很离谱。要是考场上有这手感还愁打不出来数据结构？？？

#include<bits/stdc++.h>
#define pb push_back
#define fi first
#define se second
#define ll long long
using namespace std;
const int N=4e5+5;
int Q,n,tot,last,L[N],R[N],ed[N];
int num,dfn[N],sz[N],id[N];
string s;
string q[200005];
string ans[200005];
vector<int>G[N];
vector<int>g[N];
struct node{
	int link,len,to[26];
}t[N];
void extend(char c){
	int cur=++tot;
	t[cur].len=t[last].len+1;
	int p=last;
	while(p!=-1&&!t[p].to[c-'a']){
		t[p].to[c-'a']=cur;
		p=t[p].link;
	}
	if(p!=-1) {
		int q=t[p].to[c-'a'];
		if(t[p].len+1==t[q].len){
			t[cur].link=q;
		}
		else{
			int clone=++tot;
			t[clone].len=t[p].len+1;
			t[clone].link=t[q].link;
			for(int i=0;i<26;i++)t[clone].to[i]=t[q].to[i];
			while(p!=-1&&t[p].to[c-'a']==q){
				t[p].to[c-'a']=clone;
				p=t[p].link;
			}
			t[q].link=t[cur].link=clone;
		}
	}
	last=cur;
} 
void dfs(int u){
	sz[u]=1,dfn[u]=++num;
	for(auto v:g[u]){
		dfs(v),sz[u]+=sz[v];
	}
}
int T[N<<2];
void ins(int p,int l,int r,int x,int y){
	if(l==r){
		T[p]=y;
		return;
	}
	int mid=l+r>>1;
	x<=mid?ins(p<<1,l,mid,x,y):ins(p<<1|1,mid+1,r,x,y);
	T[p]=max(T[p<<1],T[p<<1|1]);
}
int qry(int p,int l,int r,int ql,int qr){
	if(ql<=l&&r<=qr){
		return T[p];
	}
	int mid=l+r>>1;
	if(qr<=mid)return qry(p<<1,l,mid,ql,qr);
	if(mid<ql)return qry(p<<1|1,mid+1,r,ql,qr);
	return max(qry(p<<1,l,mid,ql,qr),qry(p<<1|1,mid+1,r,ql,qr));
}
int ask(int x,int y){
	if(x&&qry(1,1,num,dfn[x],dfn[x]+sz[x]-1)>=y){
		return 1;
	}
	return 0;
}
void work(int j){
	int len=q[j].size();
	int it=0,k=0;
	while(k<len&&t[it].to[q[j][k]-'a']){
		id[k]=it,it=t[it].to[q[j][k]-'a'],k++;
	}
	id[k]=it;
	for(;~k;k--){
		int down=(k==len)?0:q[j][k]-'a'+1;
		for(int l=down;l<26;l++){
			if(ask(t[id[k]].to[l],L[j]+k)){
				for(int i=0;i<k;i++){
					ans[j]+=q[j][i];
				}
				ans[j]+=('a'+l);
				return;
			}
		}
	}
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0); 
	cin>>s>>Q,n=s.size();t[0].link=-1;memset(T,-1,sizeof T);
	for(int i=0;i<n;i++)extend(s[i]),ed[i]=last;
	for(int i=1;i<=Q;i++){
		cin>>L[i]>>R[i]>>q[i],L[i]--,R[i]--;
		G[R[i]].pb(i);
	}
	for(int i=1;i<=tot;i++){
		g[t[i].link].pb(i);
	}
	dfs(0);
	for(int i=0;i<n;i++){
		ins(1,1,num,dfn[ed[i]],i);
		for(auto j:G[i]){
			work(j);
		}
	}
	for(int i=1;i<=Q;i++){
		if(!ans[i].size()){
			cout<<"-1"<<"\n";
		}
		else {
			cout<<ans[i]<<"\n";
		}
	}
}