1285. 单词 ac自动机 + fail树

传送门

文章目录

题意:

一篇论文由若干单词构成,且单词间是隔开的,给你 n n n个单词,要求你计算每个单词在论文中出现了多少次。

1 ≤ n ≤ 200 1\le n\le 200 1n200,所有单词总长不超过 1 e 6 1e6 1e6

思路:

考虑一种比较暴力的写法,我们将所有串仍 a c ac ac自动机里面,让后遍历每个单词,将其视为要匹配的串,每次都跑一遍,复杂度 O ( n 1 e 6 ) O(n1e6) O(n1e6),还带常数,给了 1 s 1s 1s显然过不去。

考虑每个单词在其他单词中出现的位置,都是一些前缀的后缀!所以我们可以跑每个前缀,让后暴跳 f a i l fail fail计算答案,这个过程显然可以用 f a i l fail fail树来优化,我们将所有前缀都打上标记,让后做树形 d p dp dp向上传递标记,这样到某个单词的时候,这个值就是他在所有单词中出现的次数。

或者可以不需要跑 f a i l fail fail树,直接利用队列中的拓扑结构转移也可。

// Problem: 单词
// Contest: AcWing
// URL: https://www.acwing.com/problem/content/description/1287/
// Memory Limit: 128 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int N=2000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;

int n;
int tr[N][26],idx,ed[N],fail[N];
int cnt[N],pos[N],ans[N],cnt1[N];
string s[N];
map<string,int>mp;
vector<int>v[N];

int newnode() {
	int u=++idx;
	for(int i=0;i<26;i++) tr[u][i]=0;
	return u; 
}

void insert(int id,string s) {
	int n=s.length(),p=0;
	for(int i=0;i<n;i++) {
		int u=s[i]-'a';
		if(!tr[p][u]) tr[p][u]=newnode();
		p=tr[p][u];
		cnt[p]++;
	}
	ed[p]=id;
}

void get_fail() {
	queue<int>q;
	for(int i=0;i<26;i++) if(tr[0][i]) q.push(tr[0][i]);
	while(q.size()) {
		int u=q.front(); q.pop(); 
		for(int i=0;i<26;i++) {
			int p=tr[u][i];
			if(p) fail[p]=tr[fail[u]][i],q.push(p);
			else tr[u][i]=tr[fail[u]][i];
		}
	}
}

void dfs(int u) {
	for(auto x:v[u]) dfs(x),cnt[u]+=cnt[x];
}

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);
	
	cin>>n;
	for(int i=1;i<=n;i++) {
		cin>>s[i];
		if(!mp.count(s[i])) mp[s[i]]=i;
		insert(mp[s[i]],s[i]);
	}
	get_fail();
	for(int i=1;i<=idx;i++) v[fail[i]].pb(i);
	dfs(0);
	// for(int i=idx;i>=1;i--) cnt[fail[v[i]]]+=cnt[v[i]];
	for(int i=1;i<=idx;i++) if(ed[i]) ans[ed[i]]+=cnt[i];
	for(int i=1;i<=n;i++) printf("%d\n",ans[mp[s[i]]]);
	


	return 0;
}
/*

*/









  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值