2015 Multi-University Training Contest 八(hdu5384,AC自动机)
2015 Multi-University Training Contest 8(hdu5384,AC自动机)
题意:
给n个模式串,m个匹配串
求每个匹配串还有多少个模式串(可以重复,也就是说要将fail指针稍微改一下),AC自动机时间复杂度是O(n*m),刚开始一看时间10S,数据正好可以搞,1A。。。然后发现跑了0.3s就A了,就感觉奇怪了,一看时间限制是1s,我擦,那为什么还A了?。。。。。我擦,题目都说总共的模式串加起来不超过6*10^5…我擦,1s够了。。。。
请叫我四眼田鸡QAQ....
#include <bits/stdc++.h>
using namespace std;
const int kind = 26;
struct node
{
node *fail; //失败指针
node *next[kind]; //Tire每个节点的26个子节点(最多26个字母)
int count; //是否为该单词的最后一个节点
node() //构造函数初始化
{
fail=NULL;
count=0;
memset(next,NULL,sizeof(next));
}
}*q[500001]; //队列,方便用于bfs构造失败指针
char keyword[600000]; //输入的单词
string str[100005]; //模式串
int head,tail; //队列的头尾指针
void insert(char *temp,node *root)
{
node *p=root;
int i=0,index;
while(temp[i])
{
index=temp[i]-'a';
if(p->next[index]==NULL) p->next[index]=new node();
p=p->next[index];
i++;
}
p->count++;
}
void build_ac_automation(node *root)
{
int i;
root->fail=NULL;
q[head++]=root;
while(head!=tail)
{
node *temp=q[tail++];
node *p=NULL;
for(i=0; i<26; i++)
{
if(temp->next[i]!=NULL)
{
if(temp==root) temp->next[i]->fail=root;
else
{
p=temp->fail;
while(p!=NULL)
{
if(p->next[i]!=NULL)
{
temp->next[i]->fail=p->next[i];
break;
}
p=p->fail;
}
if(p==NULL) temp->next[i]->fail=root;
}
q[head++]=temp->next[i];
}
}
}
}
int query(int u,node *root)
{
int i=0,cnt=0,index,len=str[u].size();
node *p=root;
while(str[u][i])
{
index=str[u][i]-'a';
while(p->next[index]==NULL && p!=root) p=p->fail;
p=p->next[index];
p=(p==NULL)?root:p;
node *temp=p;
while(temp!=root && temp->count!=-1)
{
cnt+=temp->count;
//temp->count=-1;//如果把这句话删掉,a匹配aaa会得到3而不是1,根据加题意看加不加这句话
temp=temp->fail;
}
i++;
}
return cnt;
}
int main()
{
#ifdef xxz
freopen("in.txt","r",stdin);
#endif // xxz
int n,t,m;
cin>>t;
while(t--)
{
head=tail=0;
node *root=new node();
cin>>n>>m;
for(int i = 0; i < n; i++) cin>>str[i];
while(m--)
{
cin>>keyword;
insert(keyword,root);
}
build_ac_automation(root);//建立fail函数,相当于KMP里面的next函数
for(int i = 0; i < n; i++)
{
cout<<query(i,root)<<endl;
}
}
return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。