标签:cdn 串匹配 code 过程 指针 font clu ++ add
首先,看清楚了,这是AC自动机不是自动AC机
引用AC自动机模板题上的一句话:
ovo
在学习AC自动机之前,应该先掌握一下两个前置技能:
AC自动机,通俗地讲,就是在Trie上跑KMP。AC自动机利用Trie的性质和KMP的思想,可以实现字符串的多模匹配。KMP是单模匹配,而它与AC自动机最大的区别就在fail指针的求法,其余思想基本相同。
所谓多模匹配,即给出若干个模式串和一个文本串,需要查找这些模式串在文本串中出现的情况。
AC自动机的操作分为三步:
既然是要利用Trie,自然要先建立一棵Trie了。本文以she,he,say,her,shr五个字符串为例建立一棵Trie:
其中,root为根节点,绿色节点表示该节点为某个单词的结尾,也就是结尾标记。
AC自动机的建树方法与Trie完全相同,在这里就不再赘述。
建树代码:
void add(string s)
{
int p=0;
for(int i=0;i<s.size();i++)
{
if(!ac[p].son[s[i]-‘a‘])
ac[p].son[s[i]-‘a‘]=++tot;
p=ac[p].son[s[i]-‘a‘];
}
ac[p].end++;
}
求fail指针是AC自动机最精髓的地方,也是最大的难点。不过,在掌握了KMP算法之后,理解起来也不难。
AC自动机中fail指针的作用与KMP中next数组的作用相同,都是在当前字符串失配时跳转到它指向的位置继续匹配。而AC自动机之所以能够进行多模匹配,就是因为fail指针。
在AC自动机中,fail指针用BFS来求。
步骤:
仍然以she,he,say,her,shr五个字符串为例,如图:
队列的状态是这样的:
h s
s e
e a h
a h r
h r y
r y e r
这样讲可能有点乱,请结合图和队列状态理解,不会难。
在实际实现过程中,若一直重复以上的第4、5步,时间复杂度难免会高。这里有一个巧妙的方法:当发现一个节点x没有某一个儿子s时,直接将s作为指针指向x的fail指针与s相同的这个儿子。这样实际上就是在模拟第4、5步反复查找的过程,这个指针会从上到下传递下来。因为当我们将根节点的编号设为0时,若一个节点没有儿子,就相当于这个儿子作为指针指向了根节点。这样可以更加方便地实现第4、5步。
求fail指针代码:
void build()
{
for(int i=0;i<26;i++)
if(ac[0].son[i])
{
ac[ac[0].son[i]].fail=0;
q.push(ac[0].son[i]);
}//将与根相连的节点的fail指针指向根节点并将它们入队
while(q.size())
{
int now=q.front();
q.pop();//取出队头
for(int i=0;i<26;i++)
if(ac[now].son[i])
{
ac[ac[now].son[i]].fail=ac[ac[now].fail].son[i];
q.push(ac[now].son[i]);
}
else
ac[now].son[i]=ac[ac[now].fail].son[i];
}//重复第4、5步
}
字符串匹配的思想与KMP基本相同,实现方式与Trie上查找字符串类似。将文本串从头到尾一位一位在Trie上查找,对于每一个节点,若没有被遍历过,沿着它的fail指针走,直到根节点或一个已遍历过的点。对于路径上的所有点,将答案加上它的结尾标记,即当前节点为几个字符串的结尾,然后将其结尾标记改为-1,以显示其已遍历过。
还是以这个图为例:
若文本串为yasherhs,则:
为什么走到一个已遍历过的点也要停止呢?因为若一个节点已被遍历,则沿着它的fail指针走直到根节点的这条路径上的所有节点一定都被遍历过了,若在走一遍则属于浪费时间,因此直接停止即可。
字符串匹配代码:
int get()
{
int p=0,ans=0;
for(int i=0;i<f.size();i++)
{
p=ac[p].son[f[i]-‘a‘];
for(int j=p;j && ac[j].end!=-1;j=ac[j].fail)
{
ans+=ac[j].end;
ac[j].end=-1;
}
}
return ans;
}
最后奉上完整代码:
#include<iostream>
#include<string>
#include<queue>
using namespace std;
const int N=1e6;
int tot=0,n;
string f;
queue<int> q;
struct T
{
int end=0,fail=0,son[26];分别表示结尾标记,fail指针和儿子节点
}ac[N];
void add(string s)
{
int p=0;
for(int i=0;i<s.size();i++)
{
if(!ac[p].son[s[i]-‘a‘])
ac[p].son[s[i]-‘a‘]=++tot;
p=ac[p].son[s[i]-‘a‘];
}
ac[p].end++;
}//建树
void build()
{
for(int i=0;i<26;i++)
if(ac[0].son[i])
{
ac[ac[0].son[i]].fail=0;
q.push(ac[0].son[i]);
}
while(q.size())
{
int now=q.front();
q.pop();
for(int i=0;i<26;i++)
if(ac[now].son[i])
{
ac[ac[now].son[i]].fail=ac[ac[now].fail].son[i];
q.push(ac[now].son[i]);
}
else
ac[now].son[i]=ac[ac[now].fail].son[i];
}
}//求fail指针
int get()
{
int p=0,ans=0;
for(int i=0;i<f.size();i++)
{
p=ac[p].son[f[i]-‘a‘];
for(int j=p;j && ac[j].end!=-1;j=ac[j].fail)
{
ans+=ac[j].end;
ac[j].end=-1;
}
}
return ans;
}//字符串匹配
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
string s;
cin>>s;
add(s);
}//输入模式串并插入Trie
ac[0].fail=0;//将根节点的fail指针指向自己,其实这步可以不要,因为默认就是0
build();//求fail指针
cin>>f;
cout<<get()<<endl;//输入文本串并匹配,直接输出答案
return 0;
}
习题:
声明:本文部分内容参考了一些大佬的博客
参考资料:
2019.5.2 于厦门外国语学校石狮分校
标签:cdn 串匹配 code 过程 指针 font clu ++ add
原文地址:https://www.cnblogs.com/TEoS/p/11384548.html