-
Notifications
You must be signed in to change notification settings - Fork 2
/
ACA0.cpp
105 lines (94 loc) · 1.86 KB
/
ACA0.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
/*
ID: mfs6174
PROG: AC auto
LANG: C++
*/
#include<iostream>
#include<fstream>
#include<string>
#include<sstream>
#include<cstring>
#include<algorithm>
#include<map>
#include<vector>
#define NODE 10000 //最大可能出现的节点数
#define CH 26 //分支数
using namespace std;
ifstream inf("ti.in");
//ofstream ouf("ti.out");
//freopen("ti.i","r",stdin);
const int maxlongint=2147483647;
int i,j,k,n,t;
int cc,c;//cc是统计使用了多少个节点
int zh[NODE][CH];//自动机机体
int shu[NODE];//相应节点的数据域
int fail[NODE];//失败指针,貌似是指向[彻底]失败的位置
int sn[300];//每个字符的代号,无效字符是0
int q[NODE];//队列
void init() //每次都要先执行
{
fail[0]=0;
memset(zh[0],0,sizeof(zh[0]));
cc=0;
}
void ins(char *s, int d) //建立trie
{
int p=0;
for(;*s;s++)
{
int t=sn[*a];
if(!zh[p][t])
{
cc++;
memset(zh[cc],0,sizeof(zh[cc]));
shu[cc] = 0;
zh[p][t] =cc;
}
p =zh[p][t];
}
shu[p]=d;
}
void acinit()//自动机初始化,执行完以后zh里就是goto或fail的位置
{
int *s=q,*e=q,i;
for (i=0;i<=CH;i++) //先把紧邻root的fail设成0
if (zh[0][i])
{
fail[zh[0][i]]=0;
*e++=zh[0][i];
}
while (s!=e)
{
int p=*s++;
for (i=0;i<=CH;i++)
{
if (zh[p][i])//如果goto存在
{
int v=zh[p][i];
*e++=v;
fail[v]=zh[fail[p]][i];//fail直接设成彻底fail后应该转移的位置(最长后缀位置)
}
else
zh[p][i]=zh[fail[p]][i];//goto不存在,转移到fail位置
}
}
}
void com(char *s) //查找每个模式串出现的次数
{
int p=0;
for (;*s;s++)
{
p=zh[p][sn[*s]];//转移
int t=p;
while (t)//跟随fail找符合的模式
{
if (shu[t])
ci[shu[t]]++;
t=fail[t];
}
}
}
int main()
{
return 0;
}