[BZOJ3940]:[Usaco2015 Feb]Censoring（AC自动机）

题目传送门

---

题目描述

FJ把杂志上所有的文章摘抄了下来并把它变成了一个长度不超过10⁵的字符串S。他有一个包含n个单词的列表，列表里的n个单词
记为t₁…t_N。他希望从S中删除这些单词。
FJ每次在S中找到最早出现的列表中的单词(最早出现指该单词的开始位置最小)，然后从S中删除这个单词。他重复这个操作直到S中
没有列表里的单词为止。（注意删除一个单词后可能会导致S中出现另一个列表中的单词）
FJ注意到列表中的单词不会出现一个单词是另一个单词子串的情况，这意味着每个列表中的单词在S中出现的开始位置是互不相同的
请帮助FJ完成这些操作并输出最后的S。

---

输入格式

第一行包含一个字符串S。
第二行包含一个整数N接下来的N行，每行包含一个字符串，第i行的字符串是t_i​​。

---

输出格式

一行，输出操作后的S。

---

样例

样例输入：

begintheescapexecutionatthebreakofdawn
2
escape
execution

样例输出：

beginthatthebreakofdawn

---

数据范围与提示

n≤2000

---

题解

不难想到这样一道题：

给出两个字符串S和T，每次从前往后找到S的一个子串A=T并将其删除，空缺位依次向前补齐，重复上述操作多次，直到S串中不含T串。输出最终的S串。（原题见[BZOJ3942]:[Usaco2015 Feb]Censoring）

很显然是一道KMP水题，大体思路就是求T的nxt数组，将S依次入栈，跑KMP，在此就不做过多的赘述。

那么这时候出现了多个串，怎么办呢？

暴力hash时间复杂度O(n²)，显然不能接受。（虽说我有个队友直接卡常A掉了叭~）

显然这道题与BZOJ3942的区别在于这道题是多模式串匹配，而那道题是单模式串匹配，不难想到AC自动机。

在此不会AC自动机的大佬可以看我这篇博客：AC自动机讲解+[HDU2222]Keywords Search（AC自动机）。

那么我们再结合BZOJ3942的思想，首先Trie树中end数组的意义定义为在当前点结束的串的长度是多少，然后匹配的时候开一个栈，边匹配边压栈，发现匹配到的当前点end[p]不为0（即为有串在当前点结束），则需要暴力将栈中后end[p]位弹出（即为将栈中的这个单词弹出）。需要注意的一点是，弹栈之后要将当前匹配的节点恢复到栈顶元素的位置，然后从下一个元素开始继续匹配。

---

代码时刻

#include<bits/stdc++.h>
using namespace std;
int cnt,sum,nxt[100001],trie[100001][30],end[100001],que[100001];
char s[100001],str[100001];
int ans[100001],ani[100001];
int n;
void insert(char *str)
{
	int len=strlen(str);
	int p=1;
	for(int k=0;k<len;k++)
	{
		int ch=str[k]-'a';
		if(!trie[p][ch])
			trie[p][ch]=++cnt;
		p=trie[p][ch];
	}
	end[p]=len;//end数组记录在当前点结束的串的长度
}
void _doudou()
{
	for(int i=0;i<26;i++)
		trie[0][i]=1;
	que[1]=1;
	for(int head=1,tail=1;head<=tail;head++)
		for(int i=0;i<26;i++)
			if(!trie[que[head]][i])trie[que[head]][i]=trie[nxt[que[head]]][i];
			else
			{
				que[++tail]=trie[que[head]][i];
				nxt[trie[que[head]][i]]=trie[nxt[que[head]]][i];
			}
}
void find(char *str)
{
	int len=strlen(str),p=1;
	for(int i=0;i<len;i++)
	{
		int ch=str[i]-'a';
		p=trie[p][ch];
		ans[++sum]=p;
		ani[sum]=i;//记录栈中每一个点在原串中的位置
		if(end[p])//发现有单词在这个点结束
		{
			sum-=end[p];//暴力弹栈
			p=ans[sum];
		}
	}
}
int main()
{
	scanf("%s",str);
	cnt=1;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%s",s);
		insert(s);
	}
	_doudou();
	find(str);
	for(int i=1;i<=sum;i++)//输出
		cout<<str[ani[i]];
	return 0;
}

cpp

---

下面是hash暴力O(n²)过掉这道题的队友，疯狂register卡常，想打正解的大佬就不要往下看啦哈哈

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<map>
#define ull unsigned long long
using namespace std;
const int p=131;
char s[100050];
char s1[100050];
ull h[2005];
ull hstack[100050];
ull base[100050];
int istack[100050];
int l[2005];
int top,n,len;
int main(){
	scanf("%s",s+1);
	len=strlen(s+1);
	scanf("%d",&n);
	base[0]=1;
	for(register int i=1;i<=n;i++){
		scanf("%s",s1+1);
		l[i]=strlen(s1+1);
		for(int j=1;j<=l[i];j++)
			h[i]=h[i]*p+s1[j]-'a'+1;
	}
	for(register int i=1;i<=100050;i++)
		base[i]=base[i-1]*p;
	for(register int i=1;i<=len;i++){
		istack[++top]=i;
		hstack[top]=hstack[top-1]*p+s[i]-'a'+1;
		for(register int j=1;j<=n;j++)
			if(hstack[top]-hstack[top-l[j]]*base[l[j]]==h[j]){
				top-=l[j];break;
			}
	}
	for(register int i=1;i<=top;i++)
		putchar(s[istack[i]]);
	return 0;
}

---

rp++