【AC自动机】hdu2222 hdu2896 hdu3065 zoj3430 poj2778 hdu2243

最新推荐文章于 2018-06-13 19:32:48 发布
最新推荐文章于 2018-06-13 19:32:48 发布
阅读量291 收藏
点赞数
分类专栏： 后缀数组、自动机 文章标签： ac自动机
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/acmer_cheer/article/details/50479899
版权

AC自动机用于多个模式串与多个母串的匹配。
第一步：根据模式串建立字典树

int len=strlen(w), r=root;
    for(int i=0;i<len;++i)
    {
    if(tree[r].ch[w[i]])r=tree[r].ch[w[i]];
        else r=tree[r].ch[w[i]]=++cnt;
    }
    ++tree[r].cnt;//cnt为在该节点结束的模式串的数量

第二步：计算每一个节点的fail指针（与kmp中next数组相似）。首先找到u的父亲的fail指针v。若v对应的儿子不为空，则u的fail指针指向v。否则访问v的fail指针。

void bfs()
{
    head=tail=0;
    q[++tail]=0;
    int p, tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<128;++i)//根据模式串的字符大小确定i的范围
        {
            if(tree[p].ch[i])
            {
                tmp=tree[p].ch[i];
                if(p)tree[tmp].fail=tree[tree[p].fail].ch[i];
                q[++tail]=tmp;
            }
            else tree[p].ch[i]=tree[tree[p].fail].ch[i];
        }
    }
}

第三步：在字典树中查找

void query()
{
    int len=strlen(w), ans=0, p=0, tmp;
    for(int i=0;i<len;++i)
    {
        p=tree[p].ch[w[i]-'a'];
        tmp=p;
        while(tree[tmp].pos)
        {
            ans+=tree[tmp].pos;
            tree[tmp].pos=0;
            tmp=tree[tmp].fail;
        }
    }
    printf("%d\n",ans);
}

模板题：
hdu2222
题目大意：问在一个母串中有多少个模式串出现
直接上模板

#include <iostream>
#include <cstdio>
#include <cstring>
#define MAXN 500005
#define MAXM 1000005
using namespace std;

int n;
char w[MAXM];

struct node
{
    int pos, ch[26], fail;
    inline void init()
    {
        fail=pos=0;
        memset(ch,0,sizeof ch);
    }
}tree[MAXN];
int cnt, root;

void add()
{
    int len=strlen(w), r=root;
    for(int i=0;i<len;++i)
    {
        w[i]-='a';
        if(tree[r].ch[w[i]])r=tree[r].ch[w[i]];
        else r=tree[r].ch[w[i]]=++cnt;
    }
    ++tree[r].pos;
}

int head, tail, q[MAXN];
void bfs()
{
    head=tail=0;
    q[++tail]=0;
    int p, tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<26;++i)
        {
            if(tree[p].ch[i])
            {
                tmp=tree[p].ch[i];
                if(p)tree[tmp].fail=tree[tree[p].fail].ch[i];
                q[++tail]=tmp;
            }
            else tree[p].ch[i]=tree[tree[p].fail].ch[i];
        }
    }
}

void query()
{
    int len=strlen(w), ans=0, p=0, tmp;
    for(int i=0;i<len;++i)
    {
        p=tree[p].ch[w[i]-'a'];
        tmp=p;
        while(tree[tmp].pos)
        {
            ans+=tree[tmp].pos;
            tree[tmp].pos=0;
            tmp=tree[tmp].fail;
        }
    }
    printf("%d\n",ans);
}

int main()
{
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
        cnt=root=0;
        scanf("%d",&n);
        for(int i=0;i<n;++i)
        {
            scanf("%s",w);
            add();
        }
        bfs();
        scanf("%s",w);
        query();
        for(int i=0;i<=cnt;++i)
            tree[i].init();
    }
    return 0;
}

hdu2896
一定要注意模式串字符的范围。蒟蒻在此处RE了几次才发现…
提供指针版
不过还是数组版的好调试一些 233

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define MAXN 100005
#define MAXM 10005
using namespace std;

int n, m;
bool vis[MAXN];
struct node
{
    int pos;
    node *ch[128], *fail;
    inline void init()
    {
        fail=0, pos=0;
        memset(ch,0,sizeof ch);
    }
}tree[MAXN], *cnt, *root;

char w[MAXM];
void add(int j)
{
    int len=strlen(w), v;
    node *r=root;
    for(int i=0;i<len;++i)
    {
        v=w[i];
        if(r->ch[v])r=r->ch[v];
        else r=r->ch[v]=++cnt;
    }
    r->pos=j;
    vis[j]=1;
}

node *q[MAXN];
int head, tail;
void bfs()
{
    head=tail=0;
    q[++tail]=root;
    node *p, *son, *tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<128;++i)
        {
            son=p->ch[i];
            if(son)
            {
                if(p==root)son->fail=p;
                else
                {
                    tmp=p->fail;
                    while(tmp)
                    {
                        if(tmp->ch[i])
                        {
                            son->fail=tmp->ch[i];
                            break;
                        }
                        tmp=tmp->fail;
                    }
                    if(!tmp)son->fail=root;
                }
                q[++tail]=son;
            }
        }
    }
}

int num[MAXM], tmp, tot;
void query(int j)
{
    int len=strlen(w), v;
    node *p=root, *temp;
    for(int i=0;i<len;++i)
    {
        v=w[i];
        while(!p->ch[v]&&p!=root)
            p=p->fail;
        p=p->ch[v];
        if(!p)p=root;
        temp=p;
        while(vis[temp->pos])
        {
            num[++tmp]=temp->pos;
            vis[temp->pos]=0;
            temp=temp->fail;
        }
    }
    if(tmp)
    {
        sort(num+1,num+tmp+1);
        printf("web %d:",j);
        for(int i=1;i<=tmp;++i)
        {
            printf(" %d",num[i]);
            vis[num[i]]=1;
        }
        ++tot, tmp=0;
        puts("");
    }
}

int main()
{
    while(~scanf("%d",&n))
    {
        root=cnt=tree, tot=0;
        for(int i=1;i<=n;++i)
        {
            scanf("%s",w);
            add(i);
        }
        bfs();
        scanf("%d",&m);
        for(int i=1;i<=m;++i)
        {
            scanf("%s",w);
            query(i);
        }
        printf("total: %d\n",tot);
        for(node *p=tree;p<=cnt;++p)
            p->init();
    }
    return 0;
}

hdu3065
这里就涉及到了计数的问题。
不进行标记。在之前的模板中都加了一个优化，即访问了一个节点就打上标记。

#include <iostream>
#include <cstdio>
#include <cstring>
#define MAXM 2000005
#define MAXN 50005
using namespace std;

int n, ans[1005];
char w[1005][55], s[MAXM];

struct node
{
    node *ch[26], *fail;
    int pos;
    void init()
    {
        fail=0, pos=0;
        memset(ch,0,sizeof ch);
    }
}tree[MAXN], *root, *cnt;

void add(node *r,int j)
{
    int len=strlen(w[j]), v;
    for(int i=0;i<len;++i)
    {
        v=w[j][i]-'A';
        if(r->ch[v])r=r->ch[v];
        else r=r->ch[v]=++cnt;
    }
    r->pos=j;
}

int head, tail;
node *q[MAXN];
void bfs()
{
    head=tail=0;
    q[++tail]=root;
    node *p, *son, *tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<26;++i)
        {
            son=p->ch[i];
            if(son)
            {
                if(p==root)son->fail=p;
                else
                {
                    tmp=p->fail;
                    while(tmp)
                    {
                        if(tmp->ch[i])
                        {
                            son->fail=tmp->ch[i];
                            break;
                        }
                        tmp=tmp->fail;
                    }
                    if(!tmp)son->fail=root;
                }
                q[++tail]=son;
            }
        }
    }
}

void query()
{
    int len=strlen(s), v;
    node *p=root, *temp;
    for(int i=0;i<len;++i)
    {
        if(s[i]<'A'||s[i]>'Z')
        {
            p=root;
            continue;
        }
        v=s[i]-'A';
        while(p)
        {
            if(p->ch[v])
            {
                p=p->ch[v];
                break;
            }
            p=p->fail;
        }
        if(!p)p=root;
        else
        {
            temp=p;
            while(temp)
            {
                if(temp->pos)++ans[temp->pos];
                temp=temp->fail;
            }
        }
    }
}

int main()
{
    while(~scanf("%d",&n))
    {
        root=cnt=tree;
        for(int i=1;i<=n;++i)
        {
            scanf("%s",w[i]);
            add(root,i);
        }
        bfs();
        scanf("%s",s);
        query();
        for(int i=1;i<=n;++i)
            if(ans[i])
            {
                printf("%s: %d\n",w[i],ans[i]);
                ans[i]=0;
            }
        for(node *p=tree;p<=cnt;++p)
            p->init();
    }
    return 0;

zoj3430
题目大意：给你一个加密规则：将所有字符写成二进制并串联起来。然后每6位数组成一个新的二进制数，再转化为十进制数，根据密码表翻译成字符。若len%3==1，就再加上=。若len%3==2，就加上==。
现在已知一些被加密后的模式串和一些加密后的母串。求每一个母串中出现了多少个模式串。

这道题巧妙地运用位运算可以很方便的还原字符串。要注意虽然密码串是从‘0’到‘z’，但原串是0到255。而且是多组数据，一定要清零。被坑惨了…

#include <iostream>
#include <cstdio>
#include <cstring>
#define MAXN 50050
using namespace std;

char w[MAXN];
bool vis[MAXN];
int n, temp[MAXN], cnt;

struct node
{
    int ch[256], fail, cnt;
    void init()
    {
        cnt=fail=0;
        memset(ch,0,sizeof ch);
    }
}tree[MAXN];

void add(int r)
{
    for(int i=1;i<=temp[0];++i)
    {
        if(tree[r].ch[temp[i]])r=tree[r].ch[temp[i]];
        else r=tree[r].ch[temp[i]]=++cnt;
    }
    tree[r].cnt=1;
}

int key[256];
void table()
{
    for(int i=0;i<26;++i)key[i+'A']=i;
    for(int i=0;i<26;++i)key[i+'a']=i+26;
    for(int i=0;i<10;++i)key[i+'0']=i+52;
    key['+']=62, key['/']=63;
}

void change(char s[])
{
    temp[0]=0;
    int len, x=0;
    for(len=strlen(s);s[len-1]=='=';--len);
    for(int i=0, tmp=0;i<len;++i)
    {
        x=(x<<6)|key[s[i]], tmp+=6;
        if(tmp>=8)
        {
            temp[++temp[0]]=(x>>(tmp-8))&255;
            tmp-=8;
        }
    }
}

int head, tail, q[MAXN];
void bfs()
{
    head=tail=0;
    q[++tail]=0;
    int p, tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<256;++i)
        {
            if(tree[p].ch[i])
            {
                tmp=tree[p].ch[i];
                if(p)tree[tmp].fail=tree[tree[p].fail].ch[i];
                q[++tail]=tmp;
            }
            else tree[p].ch[i]=tree[tree[p].fail].ch[i];
        }
    }
}

void query(int w[])
{
    int ans=0, p=0, tmp;
    for(int i=1;i<=w[0];++i)
    {
        p=tree[p].ch[w[i]];
        tmp=p;
        while(vis[tmp])
        {
            ans+=tree[tmp].cnt;
            vis[tmp]=0;
            tmp=tree[tmp].fail;
        }
    }
    printf("%d\n",ans);
}

int main()
{
    table();
    while(~scanf("%d",&n))
    {
        cnt=0;
        for(int i=1;i<=n;++i)
        {
            scanf("%s",w);
            change(w);
            add(0);
        }
        bfs();
        scanf("%d",&n);
        while(n--)
        {
            memset(vis,1,sizeof vis);
            scanf("%s",w);
            change(w);
            query(temp);
        }
        for(int i=0;i<=cnt;++i)
            tree[i].init();
        puts("");
    }
    return 0;
}

poj2778
题目大意：有 n 个病毒的DNA序，求长度为l的DNA序中不含病毒的个数。

不得不承认蒟蒻是看了题解的…
首先对这 n <script type="math/tex" id="MathJax-Element-5">n</script>个病毒DNA序构建AC自动机。自动集中的节点都是由有向边连接而成的，那么将那些打了标记的节点删除，剩下的就是一个图。问题就转化为从起点0开始，走N步有多少种方案。
这样就是一个矩阵加速的问题了。

#include <iostream>
#include <cstdio>
#include <cstring>
#define MAXN 110
#define LL long long int
#define mod 100000
using namespace std;

struct mat
{
    LL num[MAXN][MAXN], n;
    void init()
    {
        memset(num,0,sizeof num);
        n=0;
    }
    mat operator * (const mat &a)const
    {
        mat ans;
        ans.init();
        ans.n=n;
        for(int i=0;i<=n;++i)
            for(int j=0;j<=n;++j)
                for(int k=0;k<=n;++k)
                    ans.num[i][k]=(ans.num[i][k]+num[i][j]*a.num[j][k])%mod;
        return ans;
    }
}ans;

mat power(mat a,int pos)
{
    mat ans=a;
    while(pos)
    {
        if(pos&1)ans=ans*a;
        a=a*a;
        pos>>=1;
    }
    return ans;
}

inline int getid(char a)
{
    if(a=='A')return 0;
    if(a=='C')return 1;
    if(a=='G')return 2;
    return 3;
}

struct node
{
    int ch[5], fail, cnt;
    void init()
    {
        fail=cnt=0;
        memset(ch,0,sizeof ch);
    }
}tree[MAXN];

int cnt, root;
char w[MAXN];
void add()
{
    int len=strlen(w), r=root, v;
    for(int i=0;i<len;++i)
    {
        v=getid(w[i]);
        if(tree[r].ch[v])r=tree[r].ch[v];
        else r=tree[r].ch[v]=++cnt;
    }
    ++tree[r].cnt;
}

int q[MAXN], head, tail;
void bfs()
{
    head=tail=0;
    q[++tail]=0;
    int p, tmp;
    while(head<tail)
    {
        p=q[++head];
        for(int i=0;i<4;++i)
        {
            if(tree[p].ch[i])
            {
                tmp=tree[p].ch[i];
                if(p)
                {
                    tree[tmp].fail=tree[tree[p].fail].ch[i];
                    tree[tmp].cnt+=tree[tree[tree[p].fail].ch[i]].cnt;
                }
                q[++tail]=tmp;
            }
            else
                tree[p].ch[i]=tree[tree[p].fail].ch[i];
        }
    }
}

void build()
{
    ans.init();
    for(int i=0;i<=cnt;++i)
    {
        for(int j=0;j<4;++j)
        {
            if(tree[tree[i].ch[j]].cnt)continue;
            ++ans.num[i][tree[i].ch[j]];
        }
    }
    ans.n=cnt;
}

int n, m;
LL out;
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        root=cnt=0;
        for(int i=0;i<n;++i)
        {
            scanf("%s",w);
            add();
        }
        bfs();
        build();
        ans=power(ans,m-1);
        out=0;
        for(int i=0;i<=cnt;++i)out=(out+ans.num[0][i])%mod;
        printf("%d\n",out);
        for(int i=0;i<=cnt;++i)tree[i].init();
    }
    return 0;
}

hdu2243
这道题就是上面那道的加强版。
蒟蒻目前TLE中，期待持续更新…

RZXZ
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
hdu 3341(ac自动机+状态压缩)
10-29
hdu 3341(ac自动机+状态压缩) 题意：容易理解... 思路：首先一开始容易想到要用到dp，开设一个dp[41][41][41][41][501]的数组来解决，但是明显内存已经超出范围了，于是就想如何减少内存呢？只要知道A、T、C、G其中...
ac自动机.pptx
07-14
要学AC自动机需要自备两个前置技能:KMP和trie树(其实个人感觉不会kmp也行,失配指针的概念并不难) 其中,KMP是用于一对一的字符串匹配,而trie虽然能用于多模式匹配,但是每次匹配失败都需要进行回溯,如果模式串很长的话...
ZOJ 3430 Detect the Virus AC自动机
Azarath Metrion Zinthos
04-21 701
Detect the Virus Time Limit: 2 Seconds      Memory Limit: 65536 KB One day, Nobita found that his computer is extremely slow. After several hours' work, he finally found that it was a virus th
【ZOJ 3430 Detect the Virus】 AC自动机
lajiyuan的博客
06-13 241
ZOJ3430 题意就是给你一个加密的单词，加密方式就是先把字符按照ASCII表转为2进制，再截取6个为一段转为10进制，再去表中改为对应字符。 我们只需要把所有关键字和文本串都解码之后，就是一个最朴素的AC自动机了。 密文转换方法是向kuangbin大神学习的，我们通过思考可以发现，解码的过程就是把四个字符转换为3个字符的过程，所以我们只要对加密串进行一些二进制操作，就可以直接将24位二进...
zoj 3430 AC自动机
observadores的博客
04-24 226
题目链接：http://acm.hust.edu.cn/vjudge/problem/visitOriginUrl.action?id=16400 把输入的字符根据下表，转化成value数字，在转化成6位的二进制，载8个二进制一组划分，转化成0~255的数字（ASCII字符的十进制），然后就是AC自动机模版的事了。一个=表示转化得到的二进制末尾要去掉2个0 Value
ZOJ 3430 Detect the Virus ac自动机 超级精简的转化
MissZhou要努力
08-29 727
Description One day, Nobita found that his computer is extremely slow. After several hours' work, he finally found that it was a virus that made his poor computer slow and the virus was activated b
zoj3430 AC自动机模拟
JIBANCANYANG
05-06 912
/**********************jibancanyang************************** *Author* :jibancanyang *Created Time* : 五 5/ 6 16:14:59 2016 *File Name* : jy.cpp **Problem**: **Analyse**: 此题难就难在编码,还有题意问的
【AC自动机】hdu2222 hdu2896 hdu3065 zoj3430 poj2778
12-24 316
AC自动机用于多个模式串与多个母串的匹配。 第一步：根据模式串建立字典树int len=strlen(w), r=root; for(int i=0;i<len;++i) { if(tree[r].ch[w[i]])r=tree[r].ch[w[i]]; else r=tree[r].ch[w[i]]=++cnt; } ++tree[r].c
hdu与poj题目分类
qq_21577929的博客
08-07 6434
POJ初期: 一.基本算法: (1)枚举. (poj1753,poj2965) (2)贪心(poj1328,poj2109,poj2586) (3)递归和分治法. (4)递推. (5)构造法.(poj3295) (6)模拟法.(poj1068,poj2632,poj1573,poj2993,poj2996) 二.图算法:
AC自动机_AC自动机模板_
09-30
供信息学奥林匹克竞赛选手使用 AC自动机模板
高效Java敏感词过滤系统AC自动机算法源码，支持独立部署与集成注册中心
03-25
本项目是一款高效的Java敏感词过滤系统，基于AC自动机算法实现。系统支持独立部署，同时可便捷集成至注册中心，为各类项目提供敏感词过滤服务。包含文件共117个，其中主要构成如下： - Java源文件：49个 - Class...
中文高性能AC自动机代码
11-05
C语言实现，效率极高，实现了中文的关键字匹配，输出的格式为偏移量加上关键字（中文编码为GB2312）
keil5安装程序（包含keygen）
05-21
keil5安装程序（包含keygen）适用于单片机学习
grpcio-1.14.1-cp37-cp37m-manylinux1_i686.whl
最新发布
05-21
Python社区提供了大量的第三方库，如NumPy、Pandas和Requests，极大地丰富了Python的应用领域，从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径，而且为经验丰富的开发者提供了强大的工具，以高效率、高质量地完成复杂任务。例如，Matplotlib和Seaborn库在数据可视化领域内非常受欢迎，它们提供了广泛的工具和技术，可以创建高度定制化的图表和图形，帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
grpcio-1.14.0-cp36-cp36m-linux_armv7l.whl
05-21
Python社区提供了大量的第三方库，如NumPy、Pandas和Requests，极大地丰富了Python的应用领域，从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径，而且为经验丰富的开发者提供了强大的工具，以高效率、高质量地完成复杂任务。例如，Matplotlib和Seaborn库在数据可视化领域内非常受欢迎，它们提供了广泛的工具和技术，可以创建高度定制化的图表和图形，帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
CoreNLP一套Java核心自然语言处理工具，用于标记化、句子分词、NER分析、相互引用、情感分析等.zip
05-21
CoreNLP一套Java核心自然语言处理工具，用于标记化、句子分词、NER分析、相互引用、情感分析等
java+springboot图片处理Utils类文件
05-21
包含方法： 1、从图片的MultipartFile获取InputStream。 2、转换图片格式为jpg，保持原尺寸不变，并获取新的InputStream。 3、压缩图片文件大小,图片尺寸不变。 4、按照指定比例进行缩小和放大。 5、裁剪图片。 6、给某张图片绘制标记框，并输出至指定位置。 7、MultipartFile 转 File。
vscode-1.52.0.tar源码文件
05-21
vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件 vscode-1.52.0.tar源码文件vscode-1.52.0.tar源码文件
python ac自动机
09-02
Python AC自动机是一个用于字符串匹配的算法，它可以高效地在一段文本中查找多个预定义的模式。它的实现可以使用多种库，其中包括ac自动机python和ahocorasick-python。 ac自动机python是一个对标准的ac自动机算法...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值