您的位置:首页 > 其它

hdu 3065 病毒侵袭持续中 ac自动机

2014-04-28 16:53 369 查看
#include <cstdio>
#include <cstring>
#include <queue>
#define MAXNODE 50000
#define SIGMA_SIZE 100
using namespace std;
int ch[MAXNODE][SIGMA_SIZE];
int f[MAXNODE];    // fail函数
int val[MAXNODE];  // 每个字符串的结尾结点都有一个非0的val
int last[MAXNODE]; // 输出链表的下一个结点
char web[2000100];
char vir[1010][60];
int cnt[MAXNODE];
int sz;

void init()
{
sz = 1;
memset(ch[0], 0, sizeof(ch[0]));
memset(cnt, 0, sizeof(cnt));
}

// 插入字符串。v必须非0
void insert(char *s, int v)
{
int u = 0, n = strlen(s);
for(int i = 0; i < n; i++)
{
int c = 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;
}

// 递归打印以结点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 = 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函数
void 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]];
}
}
}
int main()
{
int n;
int i,j;
while(scanf("%d",&n)!=EOF)
{
init();
for(i=1;i<=n;i++)
{
scanf("%s",vir[i]);
insert(vir[i],i);
}
getfail();
scanf("%s",web);
find(web);
for(i=1;i<=n;i++)
if(cnt[i]!=0)
printf("%s: %d\n",vir[i],cnt[i]);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: