kuangbin ac自动机 常用打印

最新推荐文章于 2021-09-21 21:23:12 发布
最新推荐文章于 2021-09-21 21:23:12 发布
阅读量85 收藏
点赞数 1
分类专栏： 题解打印
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45891413/article/details/111713098
版权

第一行给一个整数，表示有几个案例。
每个案例先给出1个整数，表示有几个短字符串。
在每个案例的最后再给出长字符。
每个案例输出一行整数，表示短字符串里有几个能在长字符里找到。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<map>
#include<stack>
#include<set>
#include<queue>
#include<vector>
#include<stdlib.h>
 
#include<iostream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cstdio>
#include<cstdlib>
 
#define inf 0x3f3f3f3f
#define cha 1e-6
#define ll long long    
using namespace std;
const int maxn=1e6+6;
int trie[maxn][26];
int cntword[maxn];
int fail[maxn];//失配指针
int cnt=0;
void insert(string s){
    int root=0;
    for(int i=0;i<s.size();i++){
        int x=s[i]-'a';

        if(trie[root][x]==0){
            trie[root][x]=++cnt;
        }
        root=trie[root][x];
    }
    cntword[root]++;
}

void getfail(){
    queue<int>qu;
    for(int i=0;i<26;i++){
        if(trie[0][i]){
            fail[trie[0][i]]=0;
            qu.push(trie[0][i]);
        }
    }
    while(!qu.empty()){
        int now=qu.front();
        qu.pop();
        for(int i=0;i<26;i++){
            if(trie[now][i]){
                fail[trie[now][i]]=trie[fail[now]][i];
                qu.push(trie[now][i]);
            }
            else{
                trie[now][i]=trie[fail[now]][i];
            }
        }
    }
}

int query(string s){
    int now=0,ans=0;
    for(int i=0;i<s.size();i++){
        now=trie[now][s[i]-'a'];
        for(int j=now;j&&cntword[j]!=-1;j=fail[j]){//已经访问过
            ans+=cntword[j];
            cntword[j]=-1;
        }
    }
    return ans;
}

int main(){
    ios::sync_with_stdio(false);
    int n;
    int t;
    cin>>t;
    while(t--){
        memset(trie,0,sizeof(trie));
        memset(cntword,0,sizeof(cntword));
        cin>>n;
        string s;
        for(int i=0;i<n;i++){
            cin>>s;
            insert(s);
        }
        fail[0]=0;
        getfail();
        cin>>s;
        cout<<query(s)<<endl;
    }
    return 0;

}

第一行，一个整数N（1<=N<=500），表示病毒特征码的个数。
接下来N行，每行表示一个病毒特征码，特征码字符串长度在20—200之间。
每个病毒都有一个编号，依此为1—N。
不同编号的病毒特征码不会相同。
在这之后一行，有一个整数M（1<=M<=1000），表示网站数。
接下来M行，每行表示一个网站源码，源码字符串长度在7000—10000之间。
每个网站都有一个编号，依此为1—M。
以上字符串中字符都是ASCII码可见字符（不包括回车）。
Output
依次按如下格式输出按网站编号从小到大输出，带病毒的网站编号和包含病毒编号，每行一个含毒网站信息。
web 网站编号: 病毒编号 病毒编号 …
冒号后有一个空格，病毒编号按从小到大排列，两个病毒编号之间用一个空格隔开，如果一个网站包含病毒，病毒数不会超过3个。
最后一行输出统计信息，如下格式
total: 带病毒网站数
冒号后有一个空格。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<map>
#include<stack>
#include<set>
#include<queue>
#include<vector>
#include<stdlib.h>
 
#include<iostream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cstdio>
#include<cstdlib>
 
#define inf 0x3f3f3f3f
#define cha 1e-6
#define ll long long    
using namespace std;
const int maxn=1e5+5;
int trie[maxn][130];
int cntword[maxn];
int fail[maxn];//失配指针
int tag[maxn],s1[maxn];
int vis[maxn];
int cnt,ans,number,d,tol;

void insert(string s){
    int root=0;
    for(int i=0;i<s.size();i++){
        int x=s[i];
        if(trie[root][x]==0){
            trie[root][x]=++cnt;
        }
        root=trie[root][x];
    }
    cntword[root]++;
    tag[root]=++number;
}

void getfail(){
    queue<int>qu;
    for(int i=0;i<130;i++){
        if(trie[0][i]){
            fail[trie[0][i]]=0;
            qu.push(trie[0][i]);
        }
    }
    while(!qu.empty()){
        int now=qu.front();
        qu.pop();
        for(int i=0;i<130;i++){
            if(trie[now][i]){
                fail[trie[now][i]]=trie[fail[now]][i];
                qu.push(trie[now][i]);
            }
            else{
                trie[now][i]=trie[fail[now]][i];
            }
        }
    }
}

void query(string s){
    int now=0;
    for(int i=0;i<s.size();i++){
        now=trie[now][s[i]];
        for(int j=now;j&&vis[j]!=1;j=fail[j]){
             ans+=cntword[j];
             if(cntword[j]>0) s1[d++]=tag[j];
             vis[j]=1;
        }
    }
}

int main(){
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    string s;
    for(int i=1;i<=n;i++){
        cin>>s;
        insert(s);
    }
    getfail();
    int t;
    cin>>t;
    for(int i=1;i<=t;i++){
        memset(vis,0,sizeof(vis));
        memset(s1,0,sizeof(s1));
        ans=0,d=0;
        cin>>s;
        query(s);
        if(ans){
            tol++;
            printf("web %d:",i);
            sort(s1,s1+d);
            for(int j=0;j<d;j++){
                printf(" %d",s1[j]);
            }
            cout<<endl;
        }
    }    
    printf("total: %d\n",tol);
    return 0;

}

第一行，一个整数N（1<=N<=1000），表示病毒特征码的个数。
接下来N行，每行表示一个病毒特征码，特征码字符串长度在1—50之间，并且只包含“英文大写字符”。任意两个病毒特征码，不会完全相同。
在这之后一行，表示“万恶之源”网站源码，源码字符串长度在2000000之内。字符串中字符都是ASCII码可见字符（不包括回车）。
Output
按以下格式每行一个，输出每个病毒出现次数。未出现的病毒不需要输出。
病毒特征码: 出现次数
冒号后有一个空格，按病毒特征码的输入顺序进行输出。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<map>
#include<stack>
#include<set>
#include<queue>
#include<vector>
#include<stdlib.h>
 
#include<iostream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cstdio>
#include<cstdlib>
 
#define inf 0x3f3f3f3f
#define cha 1e-6
#define ll long long    
using namespace std;
const int maxn=50004;
int trie[maxn][130];
int cnt[maxn];
int fail[maxn];//失配指针
int tag[maxn],s1[maxn];
int vis[maxn],num[maxn];
char p[1003][55],s[2000006];
int ans,number,d,tol,pos;

void insert(int idx){
    int lp=strlen(p[idx]),root=0;
    for(int i=0;i<lp;i++){
        int x=p[idx][i];
        if(trie[root][x]==0){
            trie[root][x]=++pos;
        }
        root=trie[root][x];
    }
    num[root]=idx;
}

void getfail(){
    queue<int>qu;
    for(int i=0;i<130;i++){
        if(trie[0][i]){
            fail[trie[0][i]]=0;
            qu.push(trie[0][i]);
        }
    }
    while(!qu.empty()){
        int now=qu.front();
        qu.pop();
        for(int i=0;i<130;i++){
            if(trie[now][i]){
                fail[trie[now][i]]=trie[fail[now]][i];
                qu.push(trie[now][i]);
            }
            else{
                trie[now][i]=trie[fail[now]][i];
            }
        }
    }
}

void query(){
    int now=0,ls=strlen(s);
    for(int i=0;i<ls;i++){
        now=trie[now][s[i]];
        for(int j=now;j;j=fail[j]){
            cnt[num[j]]++;
        }
    }
}

int main(){
    ios::sync_with_stdio(false);
    int n;
    while(cin>>n){
        memset(trie,0,sizeof(trie));
        memset(cnt,0,sizeof(cnt));
        memset(num,0,sizeof(num));
        pos=0;
        for(int i=1;i<=n;i++){
            cin>>p[i];
            insert(i);
        }
        getfail();
        cin>>s;
        query();
        for(int i=1;i<=n;i++){
            if(cnt[i])
                printf("%s: %d\n",p[i],cnt[i]);
        }
    }
    return 0;

}

First line contains two integer m (0 <= m <= 10), n (1 <= n <=2000000000). Here, m is the number of genetic disease segment, and n is the length of sequences.

Next m lines each line contain a DNA genetic disease segment, and length of these segments is not larger than 10.
Output
An integer, the number of DNA sequences, mod 100000.
Sample Input
4 3
AT
AC
AG
AA
Sample Output
36

#include<cstdio>
#include<cstring>
#include<queue>
#include<map>
using namespace std;
const int maxn = 102;
const int mod = 100000;
int trie[maxn][4], fail[maxn], tail[maxn];
int  n, m, pos;
char s[15];
map<char, int>idx;
void insert() {
	int ls = strlen(s), u = 0;
	for (int i = 0; i < ls; i++) {
		int v = idx[s[i]];
		if (trie[u][v] == 0)
			trie[u][v] = ++pos;
		u = trie[u][v];
	}
	tail[u] = 1;
}
void getfail() {
	queue<int>q;
	for (int i = 0; i < 4; i++) {
		if (trie[0][i]) {
			fail[trie[0][i]] = 0;
			q.push(trie[0][i]);
		}
	}
	while (!q.empty()) {
		int cur = q.front();
		q.pop();
		for (int i = 0; i < 4; i++) {
			if (trie[cur][i]) {
				fail[trie[cur][i]] = trie[fail[cur]][i];
				q.push(trie[cur][i]);
			}
			else trie[cur][i] = trie[fail[cur]][i];
			tail[trie[cur][i]] |= tail[trie[fail[cur]][i]]; //注意是或，只要包含病毒就不行
		}
	}
}
struct matrix {
	long long a[maxn][maxn];
	matrix() {
		memset(a, 0, sizeof(a));
	}
};
matrix  operator*(const matrix& x, const matrix& y) {
	matrix  m;
	for (int i = 0; i <= pos; ++i) 
		for (int j = 0; j <= pos; ++j) 
			for (int k = 0; k <= pos; ++k) 
				m.a[i][j] = (m.a[i][j] + x.a[i][k] * y.a[k][j]) % mod;
	return m;
}
matrix fastm(matrix a, int n) {
	matrix res;
	for (int i = 0; i <= pos; ++i) res.a[i][i] = 1;
	while (n) {
		if (n & 1) res = res * a;
		a = a * a;
		n >>= 1;
	}
	return res;
}
int main() {
	idx['A'] = 0, idx['C'] = 1;
	idx['T'] = 2, idx['G'] = 3;
	while (~scanf("%d%d", &m, &n)) {
		pos = 0;
		memset(trie, 0, sizeof(trie));
		memset(tail, 0, sizeof(tail));
		while (m--) {
			scanf("%s", s);
			insert();
		}
		getfail();
		matrix x;
		for (int i = 0; i <= pos; ++i)  //构建初始矩阵
			if (!tail[i]) //如果本身不含病毒
				for (int j = 0; j < 4; ++j)
					if (!tail[trie[i][j]]) //其子节点也不含病毒
						x.a[i][trie[i][j]]++; //那么节点i到该子节点是可行方案+1
		x = fastm(x, n);
		int ans = 0;
		for (int i = 0; i <= pos; ++i)	//统计
			ans = (ans + x.a[0][i]) % mod;
		printf("%d\n", ans);
	}
	return 0;
}
我不是手机
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
【kuangbin带你飞】专题十七AC自动机
cjlu zx的博客
11-25 183
文章目录hdu3247 - Resource Archiver思路代码 hdu3247 - Resource Archiver 思路 n的数据范围只有10，一开始想的是设置状态dp[i][j]dp[i][j]dp[i][j]表示第一维当前选择好字符串的状态为i，那么第一维最对为(1<<10)，此时AC自动机上的状态为j。转移方程即为dp[i∣vis[p]][p]=min(dp[i][j]+1)dp[i|vis[p]][p]=min(dp[i][j]+1)dp[i∣vis[p]][p]=min(d
kuangbin 模板
09-08
kuangbin 模板
kuangbin_ac自动机
weixin_43442908的博客
09-21 201
前言 kuangbinkuangbinkuangbin的acacac自动机还是比较基础的，大多数题多是状态机模型的dpdpdp，其中有两道需要用到矩阵优化来记录路径数，然后是关于在tire图上找最短路问题，还有些要状压等等，用于字符串匹配的题很少，废话不多说，先大致分个类吧。 状态机dpdpdp 前言： 状态机dpdpdp​主要为建立AC自动机，开一维AC自动机的状态，在AC自动机上计数或求最大值等等。通常用记忆化能写，但记忆化不能压空间，也不能优化（比如：矩阵优化），这是只能递推写，递推时要注意设初始状
kuangbin专题十七AC自动机总结
qq_34921856的博客
04-17 284
这个专题写的我头皮发麻，出现了好多小bug耗费了我好多时间，但总体看不算太难，只要把思路缕清就行了。 AC自动机的题目有两类，一类是字符串找子串个数的，另一类则是建立状态，然后进行dp或者矩阵快速幂。 B - 病毒侵袭 这题不算太难，但有一个坑点，就是字符都是ASCII码可见字符。 算上空格, 从32到126共95个可见字符. 不算上空格则为94个. 所以我们读入需要用gets，而字...
kuangbin专题十七 HDU2222 AC自动机模板题
Start_to_crazy的博客
01-25 254
题解： AC自动机模板题，下面是我学习的两个博客： http://blog.csdn.net/creatorx/article/details/71100840 http://blog.csdn.net/niushuai666/article/details/7002823#reply 在学习AC自动机之前最好先学习一下字典树和KMP,前者用来搭建AC自动机的
AC自动机模板
qq_34921856的博客
07-29 227
AC自动机判母串中出现几次子串的个数(hdu2222) #include&lt;iostream&gt; #include&lt;queue&gt; using namespace std; const int maxn = 500005; const int Max = 26; struct Trie { int next[maxn][Max]; int fail[m...
AC自动机模板->kuangbin
我的博客
07-29 602
http://www.cnblogs.com/kuangbin/p/3157240.html
ac自动机模板－kuangbin
lanxuan365的专栏
09-04 678
hdu2222 #include #include #include #include #include using namespace std; const int MAX=500005; struct Trie { int next[MAX][26],fail[MAX],end[MAX]; int root,L; int newnode()
AC自动机+状压DP HDU 4057
A_Bo的博客
05-01 182
题目意思有n个字符串，每个字符串都有对应的权值，那么我们要组建一个长度为m的字符串，并且要是这个字符串的权值最大，如果最大权值小于零，输出"No Rabbit after 2012!"，否则输出最大的权值解决思路这里我们选择相应的字符串，要在Trie树上选择，并且我们还要进行模式串的匹配（匹配最大的数值），那么我们就一定会用到AC自动机，并且我们这里要确定一定长度的字符串使得权值最大，那么我们就会...
AC自动机
SIOFive的专栏
08-19 1055
AC自动机 AC自动机是KMP和Trie的结合，主要处理多模板串匹配问题。下面推荐一个博客，有助于学习AC自动机。 NOTONLYSUCCESS  这里还有一个Kuangbin开的比赛，大家也可以做一下，加深对算法的理解。 下面是比赛中的题目，采用了notonlysuccess的模板。 HDU 2222 Keywords Search 题意：最裸的模板题，给定一些模
ac自动机模板
ccsu_deer
02-27 275
#include <bits/stdc++.h> #ifdef LOCAL #define debug(x) cout<<#x<<" = "<<(x)<<endl; #else #define debug(x) 1; #endif #define lson id<<1,l,mid #define rson id<&...
kuangbin的ACM模板（2018.7更新）
07-25
kuangbin说他之前的老模板有一些错误，然后这个月分享给了我们新的模板。本来这个东西是公开的免费的，但是上传资源积分至少1个，不支持免费分享。大家如果没有积分的话可以给我私信，我不忙的时候一定会发的~~~欢迎...
kuangbin的ACM模板（新）.pdf
11-15
8 1.4 AC 自动机 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.5 后缀数组 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 1.5.1 DA . ...
kuangbin的ACM模板（新）_模板_kuangbin_kuangbinacm_ACM_
09-29
《kuangbin的ACM模板（新）》是面向ACM（国际大学生程序设计竞赛）参赛者的一份宝贵的参考资料，由kuangbin分享并命名为"kuangbinacm"。这份模板聚焦于算法和数据结构的实践应用，尤其是动态规划和计算几何两大核心...
kuangbin的ACM模板（新）.zip_kuangbin ACM算法模板
09-19
《kuangbin的ACM模板（新）》是一个专门针对ACM（国际大学生程序设计竞赛，International Collegiate Programming Contest）比赛设计的算法模板，由知名程序员kuangbin提供。这个压缩包包含了一份kuangbin的ACM算法...
pillow_heif-0.17.0-pp39-pypy39_pp73-macosx_14_0_arm64.whl
07-27
基本介绍 名称与起源：Pillow，原名为PIL（Python Imaging Library），但PIL只支持Python 2版本。随着Python 3的普及，Pillow作为PIL的一个分支出现，兼容Python 3并提供更多的功能和改进。 主要功能：Pillow支持多种图像格式的打开、保存、显示以及基本的图像操作和处理，如裁剪、缩放、旋转、翻转、滤镜应用等。 跨平台性：Pillow库可以在不同的操作系统上运行，包括Windows、Linux和MacOS等。 主要功能模块 Pillow库包含多个功能模块，每个模块都提供了特定的图像处理功能。以下是一些常用的模块： Image：用于处理图像文件，提供打开、保存、调整大小、旋转、裁剪、滤镜等功能。 ImageDraw：提供在图像上绘制各种形状（如线条、矩形、圆形）和文本的功能。 ImageFont：用于加载和使用TrueType字体文件，以便在图像上绘制文本时设置字体样式、大小和颜色。 ImageFilter：提供各种滤镜效果，如模糊、锐化、边缘增强等，用于图像增强、特效处理和图像识别等应用。 ImageEnhance：用于调整图像的亮度、对比度、颜色饱和度等参数，使图像更加清晰、明亮或具有特定的调色效果。 高级功能 除了基本的图像处理功能外，Pillow还支持一些高级功能，如色彩空间转换、直方图均衡化等。这些功能可以帮助用户进行更复杂的图像处理和分析。
cykooz.resizer-3.0.0-cp310-cp310-macosx_11_0_arm64.whl
07-27
python库，可以直接安装使用，再也不怕自己pip下载，出现等待超时的现象。欢迎下载使用，使用过程中如有使用问题，请及时与我沟通交流，帮你解决！ python库，可以直接安装使用，再也不怕自己pip下载，出现等待超时的现象。欢迎下载使用，使用过程中如有使用问题，请及时与我沟通交流，帮你解决！ python库，可以直接安装使用，再也不怕自己pip下载，出现等待超时的现象。欢迎下载使用，使用过程中如有使用问题，请及时与我沟通交流，帮你解决！ python库，可以直接安装使用，再也不怕自己pip下载，出现等待超时的现象。欢迎下载使用，使用过程中如有使用问题，请及时与我沟通交流，帮你解决！ python库，可以直接安装使用，再也不怕自己pip下载，出现等待超时的现象。欢迎下载使用，使用过程中如有使用问题，请及时与我沟通交流，帮你解决！
[毕设]Delphi题库管理与试卷自动生成系统.zip
最新发布
07-27
[毕设]Delphi题库管理与试卷自动生成系统
131 种水果、蔬菜和坚果的 90380张的高质量图像数据集+图像识别与定位+该数据集相关的研究论文
07-27
<项目介绍> 包含131 种水果、蔬菜和坚果的 90380张的高质量图像数据集+该数据集相关的研究论文。可供深度学习图像识别与定位。 - 训练集大小：67700 个图像 测试集大小：22888 个图像 文件名格式：image_index_100.jpg（例如 32_100.jpg）或 r_image_index_100.jpg（例如 r_32_100.jpg）或 r2_image_index_100.jpg 或 r3_image_index_100.jpg。 “r”代表旋转水果。 “r2”表示水果绕第三轴旋转。 “100”来自图像尺寸（100x100 像素）。 水果是如何拍摄的: 1:使用 C920 摄像机拍摄水果。最好的网络摄像头之一。 2:在水果后面，我们放了一张白纸作为背景。 3:将水果和蔬菜种植在低速电机（3转/分钟）的轴中，并录制一段20秒的短片。 下载后请首先打开README.md文件（如有），仅供学习参考, 切勿用于商业用途。 --------
ACM模板：字符串处理与数学算法详解
"ACM常用算法模板文档由kuangbin提供，收录了丰富的字符串处理、数学及相关数论算法。该模板主要分为两个部分：字符串处理和数学。 在字符串处理章节，核心算法包括： 1. KMP (Knuth-Morris-Pratt)算法，用于在文本...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值