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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值