HDU 3065 病毒侵袭持续中（AC自动机）

题目链接：点击打开链接

题意：求n个串在一个很长的串中出现的次数。

思路： 典型的AC自动机模板题。

细节参见代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
typedef long long ll;
typedef long double ld;
const ld eps = 1e-9, PI = 3.1415926535897932384626433832795;
const int mod = 1000000000 + 7;
const int INF = 0x3f3f3f3f;
const int MAXS = 1000000*2 + 10;
// & 0x7FFFFFFF
const int seed = 131;
const int SIGMA_SIZE = 199;
const ll INF64 = ll(1e18);
const int maxn = 100;
const int MAXNODE = 1000*55 + 10;
int T,n,m;
int ch[MAXNODE][SIGMA_SIZE];
int f[MAXNODE]; // fail函数
int val[MAXNODE]; // 每个字符串的结尾结点都有一个非0的val
int last[MAXNODE]; // 输出链表的下一个结点
int cnt[MAXS];
int sz;
map<string, int> ms;
void init() {
sz = 1;
memset(ch[0], 0, sizeof(ch[0]));
memset(cnt, 0, sizeof(cnt[0])*(n+2));
ms.clear();
}
// 字符c的编号
int idx(char c) {
return c-'A';
}
// 插入字符串。v必须非0
void _insert(char *s, int v) {
int u = 0, n = strlen(s);
for(int i = 0; i < n; i++) {
int c = idx(s[i]);
if(!ch[u][c]) {
memset(ch[sz], 0, sizeof(ch[sz]));
val[sz] = 0;
ch[u][c] = sz++;
}
u = ch[u][c];
}
val[u] = v;
// ms[string(s)] = v;
}
// 递归打印以结点j结尾的所有字符串
void print(int j) {
if(j) {
cnt[val[j]]++;
print(last[j]);
}
}
// 在T中找模板
void _find(char* T) {
int n = strlen(T);
int j = 0; // 当前结点编号，初始为根结点
for(int i=0;i<n;i++) { // 文本串当前指针
int c = idx(T[i]);
while(j && !ch[j][c]) j = f[j]; // 顺着细边走，直到可以匹配
j = ch[j][c];
if(val[j]) print(j);
else if(last[j]) print(last[j]); // 找到了！
}
}
// 计算fail函数
int getFail() {
queue<int> q;
f[0] = 0;
// 初始化队列
for(int c = 0; c < SIGMA_SIZE; c++) {
int u = ch[0][c];
if(u) { f[u] = 0; q.push(u); last[u] = 0; }
}
// 按BFS顺序计算fail
while(!q.empty()) {
int r = q.front(); q.pop();
for(int c = 0; c < SIGMA_SIZE; c++) {
int u = ch[r][c];
if(!u) continue;
q.push(u);
int v = f[r];
while(v && !ch[v][c]) v = f[v];
f[u] = ch[v][c];
last[u] = val[f[u]] ? f[u] : last[f[u]];
}
}
}
char buf[1006][55];
char s[2000000 + 10];
int main() {
// ios::sync_with_stdio(false);
while(~scanf("%d",&n)) {
init();
for(int i=1;i<=n;i++) {
scanf("%s",buf[i]);
_insert(buf[i], i);
}
getFail();
scanf("%s",s);
_find(s);
for(int i=1;i<=n;i++) {
if(!cnt[i]) continue;
printf("%s: %d\n",buf[i],cnt[i]);
}
}
return 0;
}