暑假集训 || 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 ...

随机推荐

Java的四大基础特性
Java的四大基础特性一.抽象父类为子类提供一些属性和行为,子类根据业务需求实现具体的行为. 抽象类使用abstract进行修饰,子类要实现所有的父类抽象方法否则子类也是抽象类. 二.封装把对象 ...
3、HTML的body内标签1
一.特殊符号的表示 #代指空格 < #代指,< > #代指,> ...... #这玩意有很多,记也记不完,用的时候查一下即可: 二.p和br标签 <p>< ...
设计模式——模板模式（Template Pattern）
在读Spring源码的时候,发现Spring代码中运用了大量的模板模式,比如根据文件系统目录加载配置文件(FileSystemXmlApplicationContext),类路径加载配置文件(Clas ...
（水题）洛谷 - P1149 - 火柴棒等式
https://www.luogu.org/problemnew/show/P1149 一开始还分类重复了.在非0的dfs中居然赋值了0,脑残得一笔. 其实就按 $lead0$ 分类就好了, $lea ...
c# KeyDown KeyPress 函数中event 的 Handled属性
很奇怪的 KeyDown中的 Handled.true 只能使 Keys.Back 这类失效, 如果要使比如数字失效,必须设置一个变量 _bHandled = true 然后在紧接着会触发的 KeyP ...
poj 1733 Parity game【hash+带权并查集】
hash一下然后用带权并查集做模2下的前缀和 #include<iostream> #include<cstdio> #include<map> #include& ...
NOIP2017 赛后总结
NOIP2017 确实,一场很深刻的考试结束了. 现在也已经搞了两周的学科了,在补之前两个月的学科的内容. 距离11.12已经过去12天了. 姓名准考证号 math complexity park ...
influxdb数据库的安装部署（windows）二
1.把下载好的文件解压到目录中 2.以管理员身份运行cmd进入influxdb-1.7.6-1文件夹 3.生成influxdb数据的默认配置Config文件 4.修改Config配置信息 5.用配置好 ...
hdu 3461 Code Lock 并查集(有点难想到)★★
#include<stdio.h> #include<math.h> ]; int count; #define mod 1000000007 int find(int x) ...
UOJ #35 后缀排序哈希做法
题面 http://uoj.ac/problem/35 题解后缀数组当然可以这里用哈希做首先排序的问题在哪里在于比较两个后缀的复杂度是O(length)的但是我们可以通过找LCP来优化比较 ...

暑假集训 || AC自动机

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

随机推荐

热门专题