洛谷 P5357 【模板】AC自动机(二次加强版)

传送门


解题思路

来pyyz过掉的第一题祭
一开始在前来个模板上改了改,发现T飞了。

看一下题解发现,因为标记不会清空,所以会被卡到O(ST)。

于是考虑优化。

因为AC自动机本质上是一个有向无环图,所以考虑每次更新答案时只更新到当前节点,最后拓扑排序遍历一遍图统计答案即可。

复杂度降到了O(S+T)。

AC代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
using namespace std;
const int maxn=2e5+5;
string s[maxn];
int n,cnt,p[maxn],tot,ans[maxn],tr[maxn][30],fail[maxn],in[maxn],anss[maxn];
vector<int> num[maxn];
struct node{
	int v,next;
}e[maxn];
void insert(int u,int v){
	in[v]++;
	cnt++;
	e[cnt].v=v;
	e[cnt].next=p[u];
	p[u]=cnt;
}
void work(){
	queue<int> q; 
	for(int i=1;i<maxn;i++){
		if(in[i]==0) q.push(i);
	}
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int i=p[u];i!=-1;i=e[i].next){
			int v=e[i].v;
			ans[v]+=ans[u];
			in[v]--;
			if(in[v]==0) q.push(v);
		}
	}
}
void update(string s,int id){
	int now=0,len=s.length();
	for(int i=0;i<len;i++){
		int k=s[i]-'a';
		if(tr[now][k]) now=tr[now][k];
		else{
			tr[now][k]=++tot;
			now=tot;
		}
	}
	num[now].push_back(id);
}
void build(){
	queue<int> q;
	for(int i=0;i<26;i++){
		if(tr[0][i]) q.push(tr[0][i]);
	}
	while(!q.empty()){
		int now=q.front();q.pop();
		for(int i=0;i<26;i++){
			if(tr[now][i]){
				fail[tr[now][i]]=tr[fail[now]][i],q.push(tr[now][i]);
				insert(tr[now][i],fail[tr[now][i]]);
			}
			else tr[now][i]=tr[fail[now]][i];
		}
	}
}
void query(string s){
	int now=0,len=s.length();
	for(int i=0;i<len;i++){
		now=tr[now][s[i]-'a'];
		if(now){
			ans[now]++;
		}
	}
}
int main(){
	//freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	ios::sync_with_stdio(false);
	memset(p,-1,sizeof(p));
	cin>>n;
	for(int i=1;i<=n;i++) cin>>s[i],update(s[i],i);
	build();
	cin>>s[n+1];
	query(s[n+1]);
	work();
	for(int i=1;i<maxn;i++){
		int k=num[i].size();
		for(int j=0;j<k;j++) anss[num[i][j]]+=ans[i];
	}
	for(int i=1;i<=n;i++){
		cout<<anss[i]<<endl;
	}
	return 0;
}
posted @ 2021-10-11 15:59  尹昱钦  阅读(56)  评论(0编辑  收藏  举报