zoj 3228 Searching the String (ac自动机)

题意:

给出一个目标串,已经一些单词,每个单词有一个限制,0表示找到目标串中包含重叠部分出现的次数,1相反。

可以有多个相同的单词。

题解:

如果直接kmp 复杂度O(n^2)超时。那么可以考虑ac自动,将单词存入字典树整体匹配。最蛋疼的是多个单词重复是允许的。那么我们要存对应编号的单词在字典树中的位置这样输出处理都方便只要处理字典树上的对应节点即可。对于不可以重叠的单词的计数,添加一个数组标记上次出现的位置,然后通过判断是否重叠来累加。


#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const ll MOD=20090717;
const int maxn=100005;
const int SIZE=maxn*6;
const int alph=26;
char str[maxn],ss[10];
int ans[SIZE][2];
int last[SIZE];
int pos[maxn],type[maxn];

struct AC
{
    int next[SIZE][alph],fail[SIZE],length[SIZE];
    int root,cnt;

    void Init()
    {
        cnt=0;
        root=newNode();
    }

    int newNode()
    {
        for(int i=0;i<alph;i++)
            next[cnt][i]=-1;
        length[cnt++]=0;
        return cnt-1;
    }

    void Insert(char buff[],int id)
    {
        int now=root;
        int len=strlen(buff);
        for(int i=0,k;i<len;i++)
        {
            k=buff[i]-'a';
            if(next[now][k]==-1)
            {
                next[now][k]=newNode();
                length[next[now][k]]=i+1;
            }
            now=next[now][k];
        }
        pos[id]=now;
    }

    void build()
    {
        queue<int>Q;
        fail[root]=root;
        int now=root;
        for(int i=0;i<alph;i++)
        {
            if(next[now][i]==-1)
                next[now][i]=root;
            else
            {
                fail[next[now][i]]=root;
                Q.push(next[now][i]);
            }
        }
        while(!Q.empty())
        {
            now=Q.front();
            Q.pop();
            for(int i=0;i<alph;i++)
            {
                if(next[now][i]==-1)
                    next[now][i]=next[fail[now]][i];
                else
                {
                    fail[next[now][i]]=next[fail[now]][i];
                    Q.push(next[now][i]);
                }
            }
        }
    }

    void Search(char buff[])
    {
        int len=strlen(buff);
        int now=root;
        memset(ans,0,sizeof ans);
        memset(last,-1,sizeof last);
        for(int i=0;i<len;i++)
        {
            int k=buff[i]-'a';
            now=next[now][k];
            int temp=now;
            while(temp!=root)
            {
                ans[temp][0]++;
                if(i-last[temp]>=length[temp])
                {
                    ans[temp][1]++;
                    last[temp]=i;
                }
                temp=fail[temp];
            }
        }
    }

};
AC ac;

int main()
{
    int n,cas=1,f;
    while(scanf("%s",&str)!=EOF)
    {
        scanf("%d",&n);
        ac.Init();
        for(int i=1;i<=n;i++)
        {
            scanf("%d %s",&type[i],ss);
            ac.Insert(ss,i);
        }
        ac.build();
        ac.Search(str);
        printf("Case %d\n",cas++);
        for(int i=1;i<=n;i++)
            printf("%d\n",ans[pos[i]][type[i]]);
        puts("");
    }
    return 0;
}
/**

*/




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值