博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
数据结构与算法系列----AC自己主动机
阅读量:6696 次
发布时间:2019-06-25

本文共 4097 字,大约阅读时间需要 13 分钟。

一:概念

首先简要介绍一下AC自己主动机:Aho-Corasick automation,该算法在1975年产生于贝尔实验室,是著名的多模匹配算法之中的一个。一个常见的样例就是给出n个单词,再给出一段文章(长度是m),让你找出有多少个单词在文章里出现过。

要搞懂AC自己主动机。先得有字典树Trie的基础知识(也有人说需要KMP的知识,我认为暂且不要理会这个。

可是在看这篇文章之前,Trie字典树,你是必需要先搞懂,假设你还不理解Trie,请參考)。

与其它字符匹配不同,KMP算法是单模式串的字符匹配算法,AC自己主动机是多模式串的字符匹配算法。匹配时间复杂度是O(N)。线性复杂度!

二:算法过程(三步走)

举个样例,假如如今给出5个模式串:say she shr he her 

主串是:yasherhs

如今问你,这5个模式串有几个出如今主串里的?

OK,如今就拿这个样例来完毕这个算法的过程。

第一步:构建Trie树,这非常easy的了。

构建好后。出现下图:

第二步:构建失败指针

构建失败指针是AC自己主动机的核心所在,玩转了它也就玩转了AC自己主动机,失败指针就是。当我的主串在trie树中进行匹配的时候,假设当前节点不能再继续进行匹配。那么我们就会走到当前节点的fail节点继续进行匹配。

构造失败指针的过程概括起来就一句话:对于root的儿子节点。fail指针直接指向root,其它的全部节点(用到了BFS和队列),设这个节点上的字母为C。沿着它父亲的失败指针走。直到走到一个节点,它的儿子中也有字母为C的节点。然后把当前节点的失败指针指向那个字母为C的节点。

假设一直走到了root都没找到,那就把失败指针指向root。

构建好后,例如以下图:

针对图中红线的”h。e“这两个节点。我们想起了什么呢?对”her“中的”e“来说,e到root距离的n个字符恰好与”she“中的e向上的n个字符相等

第三步:模式匹配

匹配过程分两种情况:

(1)  当前字符匹配成功,表示从当前节点沿着树边有一条路径能够到达目标字符,此时仅仅需沿该路径走向下一个节点继续匹配就可以,目标字符串指针移向下个字符继续匹配;

(2)  当前字符不匹配,则去当前节点失败指针所指向的字符继续匹配,匹配过程随着指针指向root结束。反复这2个过程中的随意一个。直到模式串走到结尾为止。

注意:主串全部字符在匹配完后都必需要走fail节点来结束自己的旅途,相当于一个回旋,这样做的目的防止包括节点被忽略掉。

见下图,比方:我匹配到了"she",必定会匹配到该字符串的后缀”he",要想在程序中匹配到,则必须节点要走失败指针来结束自己的旅途。

三:完整代码

#include
#include
#define MAX 26//如果仅仅出现26个小写英文字母#define ROW 4#define COLUMN 10using namespace std;char pattern[ROW][COLUMN] = { "nihao","hao","hs","hsr" };char *s = "sdmfhsgnshejfgnihaofhsrnihao";struct Node{ int index;//存储模式串的下标 char x; Node *parent; Node *next[MAX]; Node *fail; Node() { index = -1;//pattern数组下标从0開始,-1代表该节点不是单词结尾 fail = nullptr; parent = nullptr; for (int i = 0; i < MAX; i++) next[i] = nullptr; }};class ACTree{public: Node *root; ACTree() { root = new Node; root->fail = root; } void Add(const char *ch, int index); //第一步 void NodeToQueue(Node *node, queue
&q); // void BuildFailPointer(); //第二步 void ACSearch(const char *s); //第三步};int main() { ACTree tree; for (int i = 0; i < ROW; i++) tree.Add(pattern[i], i); tree.BuildFailPointer(); cout << "待匹配字符串为(依次5个一组的输出):\n"; for (int i = 1; i <= strlen(s); i++) { cout << s[i]; if (i % 5 == 0) cout << " "; } cout << endl << endl; cout << "匹配结果例如以下:\n"; cout << "位置\t" << "编号\t" << "模式\n"; tree.ACSearch(s); return 0;}void ACTree::Add(const char *ch,int index){ int len = strlen(ch); if (len == 0) return; Node *p = root; for (int i = 0; i < len; i++) { int k = ch[i] - 'a'; if (p->next[k] == nullptr) { p->next[k] = new Node; p->next[k]->parent = p; p->next[k]->x = ch[i]; } p = p->next[k]; } p->index = index;//注意,在此保证输入的模式串不反复,否则index会被覆盖}void ACTree::NodeToQueue(Node *node, queue
&q){ if (node != nullptr) { for (int i = 0; i < MAX; i++) { if (node->next[i]) q.push(node->next[i]);//不知道这是干嘛的??想想BFS层次遍历的那些事 } }}void ACTree::BuildFailPointer(){ queue
q; for (int i = 0; i < MAX; i++) { if (root->next[i]) { NodeToQueue(root->next[i], q);//注意。切不可写成q.push(root->next[i]); root->next[i]->fail = root; } } Node *parent, *p; char ch; while (!q.empty()) { p = q.front(); ch = p->x; parent = p->parent; q.pop(); NodeToQueue(p, q); while (1) { if (parent->fail->next[ch - 'a'] != nullptr) { p->fail = parent->fail->next[ch - 'a']; break; } else { if (parent->fail == root) { p->fail = root; break; } else parent = parent->fail->parent; } } }}void ACTree::ACSearch(const char *s){ int len = strlen(s); if (len == 0) return; Node *p = root; int i = 0; while (i < len) { int k = s[i] - 'a'; if (p->next[k] != nullptr) { p = p->next[k]; if (p->index != -1) cout << i - strlen(pattern[p->index]) + 1 << "\t" << p->index << "\t" << pattern[p->index] << endl; i++; } else { if (p == root) i++; else { p = p->fail; if (p->index != -1) cout << i - strlen(pattern[p->index]) + 1 << "\t" << p->index << "\t" << pattern[p->index] << endl; } } } while (p != root) { p = p->fail; if(p->index!=-1) cout << i - strlen(pattern[p->index]) + 1 << "\t" << p->index << "\t" << pattern[p->index] << endl; }}
四:数据測试

返回文件夹---->

你可能感兴趣的文章
scala(7)-----IF...ELSE 语句
查看>>
dede后台反应特别慢-转
查看>>
2015年1月25日
查看>>
ubuntu查看系统版本和内核版本
查看>>
GZFramework错误(升级修改)日志
查看>>
临时忽略某些文件
查看>>
Python:while语句,for语句
查看>>
Spark:读取hdfs的文本文件
查看>>
CF961G Partitions(第二类斯特林数)
查看>>
hibernate延迟加载(get和load的区别)
查看>>
TYVJ P1081 最近距离 Label:这不是分治!!!
查看>>
为什么要使用索引?
查看>>
pycharm 基础教程
查看>>
appium界面运行过程(结合日志截图分析)
查看>>
WCF REST 工作总结
查看>>
Java - 数组排序 -- 浅析稳定性与复杂度
查看>>
bzoj3689
查看>>
Dreamweaver 制作图片热点之后,点击热点部分会有个提示框,怎么去掉
查看>>
Codeforces Round #545 (Div. 1) 简要题解
查看>>
购物商城---页面缓存oscached
查看>>