暑假集训 || AC自动机

HDU 2222

题意：给n个模式串和一个字符串，求有多少个模式串在这个字符串中出现

思路：裸题，注意数组开的大小

#include <iostream>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <algorithm>

#include <queue>

using namespace std;

typedef long long LL;

const int SZ = ;

char keystr[];

int tot;

int fail[SZ], ch[SZ][], sum[SZ];

char str[];

void insert(char s[])

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'a';

        if(!ch[p][c]) ch[p][c] = ++tot;

        p = ch[p][c];

    }

    sum[p]++;

}

queue<int> q;

void getfail()

{

    fail[] = ;

    q.push();

    while(q.size())

    {

        int u = q.front();

        q.pop();

        for(int i = ; i < ; i++)

        {

            if(ch[u][i])

            {

                if(u == ) fail[ch[u][i]] = ;

                else

                {

                    int v = fail[u];

                    while(v && ch[v][i] == ) v = fail[v];

                    fail[ch[u][i]] = ch[v][i];

                }

                q.push(ch[u][i]);

            }

            else ch[u][i] = ch[fail[u]][i];

        }

    }

}

int match(char s[])

{

    int ans = , p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'a';

        while(p && ch[p][c] == ) p = fail[p];

        p = ch[p][c];

        int v = p;

        while(sum[v] > )

        {

            ans += sum[v];

            sum[v] = ;

            v = fail[v];

        }

    }

    return ans;

}

void init()

{

    tot = ;

    memset(sum, , sizeof(sum));

    memset(ch, , sizeof(ch));

    memset(fail, , sizeof(fail));

}

int main()

{

    int T;

    scanf("%d", &T);

    while(T--)

    {

        int n;

        scanf("%d", &n);

        init();

        for(int i = ; i < n; i++)

        {

            scanf("%s", keystr);

            insert(keystr);

        }

        getfail();

        scanf("%s", str);

        printf("%d\n", match(str));

    }

    return ;

}

HDU 3065

题意：给n个模式串，一个字符串，输出在这个字符串中出现的模式串的出现次数

思路：裸题，记录idx[p] = id;//以p结尾的是第id个模式串

#include <iostream>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <algorithm>

#include <queue>

using namespace std;

typedef long long LL;

typedef long double LD;

const int SZ = ;

char keystr[][];

int tot;

int fail[SZ];

int ch[SZ][], idx[SZ], cnt[SZ];

char str[];

void insert(char s[], int id)

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'A';

        if(!ch[p][c]) ch[p][c] = ++tot;

        p = ch[p][c];

    }

    idx[p] = id;//以p结尾的是第id个模式串

}

queue<int> q;

void getfail()

{

    fail[] = ;

    q.push();

    while(q.size())

    {

        int u = q.front();

        q.pop();

        for(int i = ; i < ; i++)

        {

            if(ch[u][i])

            {

                if(u == ) fail[ch[u][i]] = ;

                else

                {

                    int v = fail[u];

                    while(v && ch[v][i] == ) v = fail[v];

                    fail[ch[u][i]] = ch[v][i];

                }

                q.push(ch[u][i]);

            }

            else ch[u][i] = ch[fail[u]][i];

        }

    }

}

void match(char s[])

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'A';

        if(s[i] < 'A' || s[i] > 'Z') c = ;

        while(p && ch[p][c] == ) p = fail[p];

        p = ch[p][c];

        int v = p;

        while(v)

        {

            if(idx[v] > ) cnt[idx[v]]++;

            v = fail[v];

        }

    }

}

void init()

{

    tot = ;

    memset(cnt, , sizeof(cnt));

    memset(ch, , sizeof(ch));

    memset(fail, , sizeof(fail));

    memset(idx, , sizeof(idx));

}

int main()

{

    int n;

    while(~scanf("%d", &n))

    {

        init();

        for(int i = ; i <= n; i++)

        {

            scanf("%s", keystr[i]);

            insert(keystr[i], i);

        }

        getfail();

        scanf("%s", str);

        match(str);

        for(int i = ; i <= n; i++)

            if(cnt[i]) printf("%s: %d\n", keystr[i], cnt[i]);

    }

    return ;

}

ZOJ 3228

题意：给n个模式串，一个字符串，分别求每个模式串在字符串中出现的次数，其中输入中0表示可以覆盖着出现，1表示不能

思路：如果没有1那种情况就是裸题，1的情况下

用last[i]记录Trie节点i在上一次匹配时所对应的字符在文本串中的位置。

用pos[i]记录Trie节点i所对应的字符在模式串中的位置。

没有重叠的判断 —— 当前字符位置 - last[当前节点] >= pos[当前节点]。

abababac - aba

pos[1] = 1 pos[2] = 2 pos[3] = 3

i = 4时，last[3] = 2 pos[3] = 3 而i - last[3] < pos[3] 所以不行

3再往前到1，last[1] = 2 pos[1] = 1 可以了，然后last[1] = 4，节点1匹配到了文本串中的第4位

太高端了。。。。

卡数组大小，记得开n*len的

#include <iostream>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <algorithm>

#include <queue>

using namespace std;

typedef long long LL;

typedef long double LD;

const int SZ = ;

char keystr[];

int tot;

int fail[SZ];

int ch[SZ][], idx[SZ], cnt[][SZ];

int last[SZ], pos[SZ], typ[SZ];

char str[];

void insert(char s[], int id)

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'a';

        if(!ch[p][c]) ch[p][c] = ++tot;

        p = ch[p][c];

        pos[p] = i+;

    }

    idx[id] = p;

}

queue<int> q;

void getfail()

{

    fail[] = ;

    q.push();

    while(q.size())

    {

        int u = q.front();

        q.pop();

        for(int i = ; i < ; i++)

        {

            if(ch[u][i])

            {

                if(u == ) fail[ch[u][i]] = ;

                else

                {

                    int v = fail[u];

                    while(v && ch[v][i] == ) v = fail[v];

                    fail[ch[u][i]] = ch[v][i];

                }

                q.push(ch[u][i]);

            }

            else ch[u][i] = ch[fail[u]][i];

        }

    }

}

void match(char s[])

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = s[i] - 'a';

        p = ch[p][c];

        int v = p;

        while(v)

        {

            cnt[][v]++;

            if(i - last[v] >= pos[v])

            {

                cnt[][v]++;

                last[v] = i;

            }

            v = fail[v];

        }

    }

}

void init()

{

    tot = ;

    memset(cnt, , sizeof(cnt));

    memset(ch, , sizeof(ch));

    memset(last, -, sizeof(last));

    memset(pos, , sizeof(pos));

}

int main()

{

    int n, tt = ;

    while(~scanf("%s", str))

    {

        init();

        scanf("%d", &n);

        for(int i = ; i < n; i++)

        {

            scanf("%d %s", &typ[i], keystr);

            insert(keystr, i);

        }

        getfail();

        match(str);

        printf("Case %d\n", ++tt);

        for(int i = ; i < n; i++)

            printf("%d\n", cnt[typ[i]][idx[i]]);

        printf("\n");

    }

    return ;

}

HDU 2457 AC自动机+DP

题意：给n种病毒序列（只含ACGT），给一串基因，为最少修改多少个碱基（？）能使得其中不含病毒，如果不能做到则输出-1

思路：考虑神仙DP

f[i][j] 表示文本串的前 i 个字符，处于 Trie 图的节点 j 且不经过终止节点（危险节点）的最少修改个数

每个点的是否是终止节点：如果它或它顺着 fail 函数能达到的点中有一个是终止节点，那么它也是终止节点

这个讲的蛮明白的：https://blog.csdn.net/human_ck/article/details/6577142

转移方程：dp[i+1][ch[p][j]] = min(dp[i+1][ch[p][j]], dp[i][p] + (j == c ? 0 : 1)); //j表示枚举把这个非危险节点的点修改成什么，如果和原来相同则操作数不变，如果不同则+1

#include <iostream>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cstring>

#include <algorithm>

#include <queue>

using namespace std;

typedef long long LL;

const int SZ = ;

const int INF = ;

char keystr[];

int tot;

int fail[SZ];

int ch[SZ][], dp[][SZ];

char str[];

bool flag[SZ];

int tran(char c)

{

    if(c == 'A') return ;

    if(c == 'G') return ;

    if(c == 'C') return ;

    if(c == 'T') return ;

}

void insert(char s[])

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i < len; i++)

    {

        int c = tran(s[i]);

        if(!ch[p][c]) ch[p][c] = ++tot;

        p = ch[p][c];

    }

    flag[p] = true;

}

queue<int> q;

void getfail()

{

    fail[] = ;

    q.push();

    while(q.size())

    {

        int u = q.front();

        q.pop();

        for(int i = ; i <= ; i++)

        {

            if(ch[u][i])

            {

                if(u == ) fail[ch[u][i]] = ;

                else

                {

                    int v = fail[u];

                    while(v && ch[v][i] == ) v = fail[v];

                    fail[ch[u][i]] = ch[v][i];

                }

                if(flag[fail[ch[u][i]]]) flag[ch[u][i]] = true;

                q.push(ch[u][i]);

            }

            else ch[u][i] = ch[fail[u]][i];

        }

    }

}

void match(char s[])

{

    int p = ;

    int len = strlen(s);

    for(int i = ; i <= len; i++)

        for(int j = ; j <= tot; j++)

            dp[i][j] = INF;

    dp[][] = ;

    for(int i = ; i < len; i++)

    {

        int c = tran(s[i]);

        for(p = ; p <= tot; p++)

        {

            if(dp[i][p] == INF) continue;

            for(int j = ; j <= ; j++)

            {

                if(flag[ch[p][j]]) continue;

                dp[i+][ch[p][j]] =  min(dp[i+][ch[p][j]], dp[i][p] + (j == c ?  : ));

            }

        }

    }

}

void init()

{

    tot = ;

    memset(ch, , sizeof(ch));

    memset(fail, , sizeof(fail));

    memset(flag, false, sizeof(flag));

}

int main()

{

    int n, tt = ;

    while(scanf("%d", &n) && n)

    {

        init();

        for(int i = ; i <= n; i++)

        {

            scanf("%s", keystr);

            insert(keystr);

        }

        getfail();

        scanf("%s", str);

        int len = strlen(str);

        match(str);

        int minn = INF;

        for(int i = ; i <= tot; i++)

            if(!flag[i]) minn = min(minn, dp[len][i]);

        if(minn == INF) minn = -;

        printf("Case %d: %d\n", ++tt, minn);

    }

    return ;

}

OTZ

暑假集训 || AC自动机的更多相关文章

【集训第二天·翻水的老师】--ac自动机+splay树
今天是第二天集训.(其实已经是第三天了,只是昨天并没有机会来写总结,现在补上) 上午大家心情都很愉快,因为老师讲了splay树和ac自动机. 但到了下午,我们的教练竟然跑出去耍了(excuse me? ...
【暑假】[实用数据结构] AC自动机
Aho-Corasick自动机算法: <功能> AC自动机用于解决文本一个而模板有多个的问题. AC自动机可以成功将多模板匹配,匹配意味着算法可以找到每一个模板在文本中出现的位置. & ...
暑假集训(2)第七弹 -----今年暑假不AC（hdu2037)
J - 今年暑假不AC Crawling in process... Crawling failed Time Limit:1000MS Memory Limit:32768KB 64 ...
今年暑假不AC - HZNU寒假集训
今年暑假不AC "今年暑假不AC?" "是的." "那你干什么呢?" "看世界杯呀,笨蛋!" "@#$%^&a ...
2015UESTC 暑假集训总结
day1: 考微观经济学去了…… day2: 一开始就看了看一道题目最短的B题,拍了半小时交了上去wa了感觉自己一定是自己想错了,于是去拍大家都过的A题,十分钟拍完交上去就A了然后B题写了一发暴力 ...
字符串处理-AC自动机
估计在OJ上刷过题的都会对AC自动机这个名词很感兴趣,同样,记得去年ACM暑期集训的时候,在最后讲到字符串部分,听说了这个算法的名字之后就对于它心向往之,AC正好是Accept的简称,字面意义上的理解 ...
「LuoguP3808」【模板】AC自动机（简单版）
题目背景通过套取数据而直接“打表”过题者,是作弊行为,发现即棕名. 这是一道简单的AC自动机模板题. 用于检测正确性以及算法常数. 为了防止卡OJ,在保证正确的基础上只有两组数据,请不要恶意提交. ...
【模版 Luogu P3808/P3796/P5357】AC自动机（简论）
浙江集训Day9,没有出任何实质性成果,只好把昨天打完的板子记一下. 该博客基于luogu的三道模版题.只有一个大致的讲解,主要提供代码给自己参考. ------------------------- ...
基于trie树做一个ac自动机
基于trie树做一个ac自动机 #!/usr/bin/python # -*- coding: utf-8 -*- class Node: def __init__(self): self.value ...

随机推荐

CodeForces 719B Anatoly and Cockroaches (水题贪心)
题意:给定一个序列,让你用最少的操作把它变成交替的,操作有两种,任意交换两种,再就是把一种变成另一种. 析:贪心,策略是分别从br开始和rb开始然后取最优,先交换,交换是最优的,不行再变色. 代码如下 ...
android 带checkbox的List
可实现点击内容即可选中 http://blog.csdn.net/harvic880925/article/details/40475367
转-sql之left join、right join、inner join的区别
left join(左联接) 返回包括左表中的所有记录和右表中联结字段相等的记录 right join(右联接) 返回包括右表中的所有记录和左表中联结字段相等的记录inner join(等值连接) 只 ...
165 Compare Version Numbers 比较版本号
比较两个版本号 version1 和 version2.如果 version1 大于 version2 返回 1,如果 version1 小于 version2 返回 -1, 除此以外返回 0.您可 ...
python staticmethod&classmethod
python中的这两种方法都通过修饰器来完成静态方法: 不需要传递类对象或者类的实例可以通过类的实例.方法名a().foo()或者类名.方法a.foo()名来访问当子类继承父类时,且实例化子类时 ...
Git团队协作 - 新feature的开发过程
新feature的开发过程建议使用SmartGit,以下是命令行操作 git checkout -b dev (对于没有分支的人)新建dev分支 git pull origin dev拉取最新数据 ...
Outlook读取奇妙清单Wunderlist日历失败的解决办法
错误: Outlook.com日历订阅奇妙清单的日历链接时报错 This calendar wasn't updated because of a problem with the publisher ...
PL/SQL笔记(1)-流程控制,循环,异常,块
流程控制 1.If,then,else,elsif(不是elseif) ' then null; endif; 2.Case 简单case表达式: 搜索型Case表达式: 3.goto语句 begin ...
Redis学习笔记1-安装配置
一.Redis安装 Redis官网:http://www.redis.io/download 注意:版本号2.4,2.6,2.8等偶数结尾为稳定版,2.5等为非稳定版本,生成环境应该使用稳定版下载解 ...
SpringBoot 2.x （6）：使用Filter、Servlet、Listener
Filter的使用: 1.常见的默认Filter有: characterEncodingFilter hiddenHttpMethodFilter httpPutFormContentFilter r ...

暑假集训 || AC自动机

暑假集训 || AC自动机的更多相关文章

随机推荐

热门专题