YBT 2.4 AC自动机

其实这个专题NOIP几乎不考

AC自动机，就是能让题自动AC的东西，是不是十分神奇

对的，就是这么神奇

AC自动机是解决多模式串与文本串匹配的问题

是KMP+Trie树的结合，也是一个毒瘤算法

Keywords Search

link

此题是AC自动机的板子，刚才说过AC自动机是解决匹配问题的，这道题便如此

板子题，记录结尾即可

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

inline int read()

{

    int f=,ans=;char c;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}

    return f*ans;

}

int t,n;

int bo[],tot,ch[][];

char str[];

void insert(char *str)

{

    int u=;

    int len=strlen(str);

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

    {

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

        if(ch[u][c]==) ch[u][c]=++tot;

        u=ch[u][c];

    }

    bo[u]++;

    return;

}

int nex[],que[];

void bfs()

{

    que[]=,nex[]=;

    int head=,tail=;

    while(head<=tail)

    {

        int t=que[head]; head++;

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

        {

            if(ch[t][i]==) ch[t][i]=ch[nex[t]][i];

            else{

                que[++tail]=ch[t][i];

                int v=nex[t];

                nex[ch[t][i]]=ch[v][i];

            }

        }

    }

}

int ans;

void find(char *str)

{

    int u=,k;

    int len=strlen(str);

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

    {

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

        k=ch[u][c],u=ch[u][c];

        while(k>)

        {

            ans+=bo[k];

            bo[k]=;

            k=nex[k];

        }

    }

    return;

}

int main()

{

    t=read();

    while(t--)

    {

        ans=;

        tot=;

        memset(ch,,sizeof(ch));

        memset(bo,,sizeof(bo));

        for(int i=;i<=;i++) ch[][i]=;

        n=read();

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

        {

            scanf("%s",str);

            insert(str);

        }

        bfs();

        scanf("%s",str);

        find(str);

        cout<<ans<<endl;

    }

}

玄武密码

link

此题求模式串在文本串前缀最大匹配长度

所以让文本串在AC自动机上走一遍，用flag记录经过的点

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

inline int read()

{

    int f=,ans=;char c;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}

    return f*ans;

}

int n,m,le[];

char cx(char u)

{

    if(u=='E') return '';

    if(u=='S') return '';

    if(u=='W') return '';

    if(u=='N') return '';

}

int fa[],ch[][],tot;

bool flag[];

char str[];

int nex[];

int que[];

char s[];

int son[];

void insert(char *str,int an)

{

    int u=;

    int len=strlen(str);

    le[an]=len;

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

    {

        int c=str[i]-'';

        if(ch[u][c]==) ch[u][c]=++tot,fa[tot]=u;

        u=ch[u][c];

    }

    son[an]=u;

    return;

}

void bfs()

{

    nex[]=,que[]=;

    int head=,tail=;

    while(head<=tail)

    {

        int x=que[head];head++;

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

        {

            if(ch[x][i]==) ch[x][i]=ch[nex[x]][i];

            else{

                que[++tail]=ch[x][i];

                nex[ch[x][i]]=ch[nex[x]][i];

            }

        }

    }

    return;

}

int ls;

void find(char *s)

{

    int u=,k;

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

    {

        int c=s[i]-'';

        k=ch[u][c];

        while(k>)

        {

            if(flag[k]) break;

            flag[k]=;

            k=nex[k];

        }

        u=ch[u][c];

    }

    return;

}

int work(int an)

{

    int u=son[an];

    for(int i=le[an];i>=;i--)

    {

        if(flag[u]) return i;

        u=fa[u];

    }

    return ;

}

int main()

{

    memset(ch,,sizeof(ch));

    tot=;

    n=read(),m=read();

    scanf("%s",s);

    for(int i=;i<;i++) ch[][i]=;

    ls=strlen(s);

    for(int i=;i<ls;i++) s[i]=cx(s[i]);

    for(int i=;i<;i++) ch[][i]=;

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

    {

        scanf("%s",str);

        int l=strlen(str);

        for(int j=;j<l;j++) str[j]=cx(str[j]);

        insert(str,i);

    }

    bfs();

    find(s);

    for(int i=;i<=m;i++) printf("%d\n",work(i));//cout<<work(i)<<endl;

}

Censoring

link

也是让文本串在AC自动机上跑，然后记录一下位置

#include<iostream>

#include<cstdio>

#include<cstring>

#include<stack>

#include<algorithm>

#include<cmath>

using namespace std;

inline int read()

{

    int f=,ans=;char c;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}

    return f*ans;

}

char str[];

int n;

int ch[][],tot;

char s[];

int que[],nex[];

int bo[];

void insert(char *s)

{

    int len=strlen(s);

    int u=;

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

    {

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

        if(ch[u][c]==) ch[u][c]=++tot;

        u=ch[u][c];

    }

    bo[u]=len;

    return;

}

void bfs()

{

    int head=,tail=;

    que[]=,nex[]=;

    while(head<=tail)

    {

        int x=que[head];head++;

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

        {

            if(ch[x][i]==) ch[x][i]=ch[nex[x]][i];

            else{

                que[++tail]=ch[x][i];

                nex[ch[x][i]]=ch[nex[x]][i];

            }

        }

    }

    return;

}

int stk1[],stk2[];

void find(char *str)

{

    int u=;

    int tot=;

    int len=strlen(str);

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

    {

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

        u=ch[u][c];

        stk1[++tot]=u;

        stk2[tot]=i;

        while(bo[u])

        {

            tot-=bo[u];

            u= stk1[tot];

        }

    }

    for(int i=;i<=tot;i++) cout<<str[stk2[i]];

    cout<<endl;

}

int main()

{

    tot=;

    for(int i=;i<;i++) ch[][i]=;

    scanf("%s",str);

    n=;

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

    {

        scanf("%s",s);

        insert(s);

    }

    bfs();

    find(str);

}

单词(待补此坑)

最短母串(待补此坑)

病毒

link

此题是不让文本串在AC自动机上跑到每个模式串的末尾

所以想这个问题，如果说你有无穷长度的文本串

在AC自动机上跑着，但是不让他经过末尾

所以说这个图是有环的，AC自动机不是因为节省时间直接相连nex吧，所以好像这个叫做fail图

所以看一下在AC自动机中不经过每个字符串尾端是否有环即可

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<queue>

using namespace std;

inline int read()

{

    int f=,ans=;char c;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}

    return f*ans;

}

char s[];

int tot,st[];

int ch[][];

int n;

int quee[],nex[];

void insert(char *s)

{

    int len=strlen(s),u=;

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

    {

        int c=s[i]-'';

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

        u=ch[u][c];

    }

    st[u]=;

    return;

}

void bfs()

{

    int head=,tail=;

    quee[]=,nex[]=;

    while(head<=tail)

    {

        int x=quee[head];head++;

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

        {

            if(ch[x][i]==) ch[x][i]=ch[nex[x]][i];

            else{

                quee[++tail]=ch[x][i];

                nex[ch[x][i]]=ch[nex[x]][i];

                if(st[nex[ch[x][i]]]) st[ch[x][i]]++;

            }

        }

    }

    return;

}

int vis1[],vis2[];

bool flag;

void dfs(int f)

{

    if(flag) return;

    vis2[f]=;

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

    {

        int s=ch[f][i];

        if(flag) return;

        if(st[s]) continue;

        if(vis2[s])

        {

            flag=;

            cout<<"TAK";

            return;

        }

        if(flag) return;

        if(vis1[s]) continue;

        vis1[s]=;

        dfs(s);

    }

    vis2[f]=;

}

int main()

{

    tot=;

    n=read();

    for(int i=;i<=;i++) ch[][i]=;

    for(int i=;i<=n;i++) scanf("%s",s),insert(s);

    bfs();

    flag=false;

    dfs();

    if(flag==) cout<<"NIE";

    return ;

}

/*

2

1

0

*/

文本生成器

link

第一次在AC自动机上跑dp

想一个问题，如果这个有m长度的文本串符合条件的话，他至少是一个模式的母串吧

你怎么去统计呢，有可能是一个，两个，三个，四个~~~

在数学上有一种想法是退而求进，所以可以把这道题转换成有多少个无法模式串匹配的情况

这时候求的就是0了，然后去拿总可能-这种可能便是答案

为什么说这一道题是dp呢，因为满足dp的思想了

dp[i]为长度为i有多少个

但是无法往前推

所以再加一维，加什么呢

很容易想到是访问在哪个节点吧

所以dp模型就已经定下来了

dp[i][j]表示现在长度为i，节点是j共有多少个串

AC自动机有什么用呢，其实就是建图有用

将字符串转换成一张图

不让他走末尾节点

所以现在就有一个要注意的坑点了

***当你找到一点x，若此点的nex是一个字符串的结尾，所以这个x其实就是一个病毒了

其实就是每人对nex理解了

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#define mod 10007

using namespace std;

inline long long read()

{

    long long f=,ans=;char c;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}

    return f*ans;

}

long long n,m,tot,u,ch[][],bo[];

char s[];

void insert(char *s)

{

    long long len=strlen(s),u=;

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

    {

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

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

        u=ch[u][c];

    }bo[u]=;

    return;

}

long long que[],nex[];

void bfs()

{

    long long head=,tail=;

    que[]=,nex[]=;

    while(head<=tail)

    {

        long long xx=que[head];head++;

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

        {

            if(ch[xx][i]==) ch[xx][i]=ch[nex[xx]][i];

            else{

                que[++tail]=ch[xx][i];

                nex[ch[xx][i]]=ch[nex[xx]][i];

                bo[ch[xx][i]]|=bo[ch[nex[xx]][i]];

            }

        }

    }

    return;

}

long long dp[][];

int main()

{

    for(long long i=;i<;i++) ch[][i]=;

    tot=;

    n=read(),m=read();

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

    {

        scanf("%s",s);

        insert(s);

    }

    bfs();

    dp[][]=;

    for(long long i=;i<=m;i++)

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

            for(long long z=;z<;z++)

                if(!bo[ch[j][z]]) dp[i][ch[j][z]]+=dp[i-][j],dp[i][ch[j][z]]%=mod;

    long long ans=;

    for(long long i=;i<=tot;i++) ans+=dp[m][i],ans%=mod;

    long long sum=;

    for(long long i=;i<=m;i++) sum*=,sum%=mod;

    cout<<(sum-ans+mod)%mod;

}

YBT 2.4 AC自动机的更多相关文章

021（Keywords Search）（AC自动机）
题目:http://ybt.ssoier.cn:8088/problem_show.php?pid=1479 题目思路:一道AC自动机的模板题备注:还不会字典树和KMP的尽早回去重修如果让你在一篇 ...
基于trie树做一个ac自动机
基于trie树做一个ac自动机 #!/usr/bin/python # -*- coding: utf-8 -*- class Node: def __init__(self): self.value ...
AC自动机-算法详解
What's Aho-Corasick automaton? 一种多模式串匹配算法,该算法在1975年产生于贝尔实验室,是著名的多模式匹配算法之一. 简单的说,KMP用来在一篇文章中匹配一个模式串:但 ...
python爬虫学习(11) —— 也写个AC自动机
0. 写在前面本文记录了一个AC自动机的诞生! 之前看过有人用C++写过AC自动机,也有用C#写的,还有一个用nodejs写的.. C# 逆袭--自制日刷千题的AC自动机攻克HDU OJ HDU 自 ...
BZOJ 2434: [Noi2011]阿狸的打字机 [AC自动机 Fail树树状数组 DFS序]
2434: [Noi2011]阿狸的打字机 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 2545 Solved: 1419[Submit][Sta ...
BZOJ 3172: [Tjoi2013]单词 [AC自动机 Fail树]
3172: [Tjoi2013]单词 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 3198 Solved: 1532[Submit][Status ...
BZOJ 1212: [HNOI2004]L语言 [AC自动机 DP]
1212: [HNOI2004]L语言 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1367 Solved: 598[Submit][Status ...
[AC自动机]【学习笔记】
Keywords Search Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)To ...
AC自动机 HDU 3065
大概就是裸的AC自动机了 #include<stdio.h> #include<algorithm> #include<string.h> #include< ...

随机推荐

VIN码识别/车架号OCR识别：快速占领汽车后市场数据入口
大数据时代,企业在数据入口方面的竞争越来越激烈,这种对于入口级的大数据“争夺战”,让很多企业在数据挖掘和收集的技术方面开始加快更新速度. 在当前IT行业激烈竞争环境之下,对于入口产品的控制成为了大数据 ...
Qt-第一个QML程序-4-自定义按钮类，动画，状态
上篇中,我们写到了自己定义了一个按钮,但是呢,按照这样的写法,要写一个程序出来,那要累死了,所以,qml给我的感觉就是各种随便调用,所以了,可以自己写一个自己Button的qml,这样在以后用到了,就 ...
对网页进行截图（selenium）
import os def insert_img(driver,file_name): #获取当前路径,并转换为字符串 base_dir=str(os.path.dirname(__file__)) ...
第一篇 HTML基础
浏览网页,就是上网,上网的本质就是下载内容. 浏览器是个解释器,用来执行HTML.css.JS代码的. HTML,CSS, JavaScript 号称网络三剑客. 1. 浏览器发送一个域名给服务端 2 ...
python基础之变量和简单数据类型
1.1 变量的命名和使用规范变量名可以包含数字.字母.下划线,但是不能以数字开头. 变量名不能包含空格,可使用下划线来分割其中的单词. 不要将Python关键字和函数名用作变量名. 变量名应既简短又 ...
Python元组与列表的区别和联系？
1. 元组和列表比较相似,不过它们之间也有着不同: (1)列表:一个大仓库,你可以随时往里边添加和删除任何东西. (2)元组:封闭的列表,一旦定义,就不可改变(不能添加.删除或修改). 2. 什么情 ...
进度条加载与案例优化对比——python使用perf_count方法实现
本章我们将讨论python3 perf_counter()的用法及它的实际应用我从中选取两个python基于rquests库的爬虫实例代码源文件进行举例 Python3 perf_counter() ...
Java学习 · 初识面向对象基础一
面向对象基础 1.1面向过程与面向对象的区别面向过程和面向对象二者都是思考问题的方式,再简单的事物时,可以线性思考时使用面向过程,但当事物较为复杂时,只能使用面向对象设计.但二者并不是对立的,在解决 ...
[C++] Variables and Basic Types
Getting Started compile C++ program source $ g++ -o prog grog1.cc run C++ program $ ./prog The libra ...
[LeetCode] 53. Maximum Subarray 解题思路
Find the contiguous subarray within an array (containing at least one number) which has the largest ...

YBT 2.4 AC自动机

YBT 2.4 AC自动机的更多相关文章

随机推荐

热门专题