1966: [Ahoi2005]VIRUS 病毒检测

1966: [Ahoi2005]VIRUS 病毒检测

Time Limit: 3 Sec   Memory Limit: 64 MB
Submit: 360   Solved: 143
[ Submit][ Status][ Discuss]

Description

科学家们在Samuel星球上的探险仍在继续。非常幸运的,在Samuel星球的南极附近,探险机器人发现了一个巨大的冰湖!机器人在这个冰湖中搜集到了许多RNA片段运回了实验基地。科学家们经过几个昼夜的研究,发现这些RNA片段中有许多是未知的病毒!每个RNA片段都是由A、C、T、G组成的序列。科学家们也总结出了Samuel星球上的“病毒模版片段”。一个模版片段是由A、C、T、G的序列加上通配符 * 和 ? 来表示。其中 * 的意思是可以匹配上0个或任意多个字符,而 ? 的意思是匹配上任意一个字母。如果一个RNA片段能够和“病毒模版片段”相匹配,那么这个RNA片段就是未知的病毒。例如,假设“病毒模版片段”为A*G?C。RNA片段:AGTC,AGTGTC都是未知的病毒,而RNA片段AGTGC则不是病毒。由于,机器人搜集的这些RNA片段中除去病毒的其他部分都具有非常高的研究价值。所以科学家们希望能够分辨出其中哪些RNA片段不是病毒,并将不是病毒的RNA片段运回宇宙空间站继续进行研究。科学家将这项任务交给了小联。现在请你为小联编写程序统计哪些RNA片段不是病毒。

Input

第一行有一个字符串,由A、C、T、G、*、? 组成。表示“病毒模版片段”。“病毒模版片段”的长度不超过1000。第二行有一个整数N(0<N<500),表示机器人搜集到的RNA片段的数目。随后的N行,每一行有一个字符串,由A、C、T、G组成,表示一个RNA片段。每个RNA片段的长度不超过500。注意:“病毒模版片段”和RNA片段的长度都至少为1。

Output

只有一行输出,为整数M,即不是病毒的RNA片段的数目。

Sample Input

A*G?C
3
AGTC
AGTGTC
AGTGC

Sample Output

1

HINT

输入中的RNA片段AGTGC不是病毒。

Source

[ Submit][ Status][ Discuss]




瞎逼dp下

f[i][j]:主串匹配到第i位,当前串匹配到第j位,是否可行

转移随意脑补,,,,(WA了好多次是什么鬼?!)


听说AC自动机也可以呀

话说,,,我这种傻逼写法要特判最后一位是否为'*'

#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<bitset>
#include<algorithm>
#include<cstring>
#include<map>
#include<stack>
#include<set>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 1010;
const int maxm = 505;

int n,m,cnt,len,ans,f[maxn][maxm];
char a[maxn],b[maxm];
bool not_cho[maxn];

int main()
{
	#ifdef DMC
		freopen("virus8.in","r",stdin);
	#endif
	
	scanf("%s",a + 1); int tot = strlen(a + 1);
	for (int i = 2; i <= tot; i++)
		if (a[i] == '*' && a[i-1] == a[i])
			not_cho[i] = 1;
	for (int i = 1; i <= tot; i++)
		if (!not_cho[i])
			a[++n] = a[i];
	cin >> m;	
	for (int l = 1; l <= m; l++) {
		scanf("%s",b + 1); len = strlen(b + 1);
		f[0][0] = ++cnt;
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= len; j++) {
				if (a[i-1] == '*') {
					if (a[i] == '?') {
						if (f[i-2][j-1] == cnt) f[i][j] = cnt;
					}
					else if (a[i] == b[j] && f[i-2][j-1] == cnt) f[i][j] = cnt;
				}
				if (a[i] == '*') {
					if (f[i][j-1] == cnt) f[i][j] = cnt;
					if (f[i-1][j-1] == cnt) f[i][j] = cnt;
				}
				else if (a[i] == '?') {
					if (f[i-1][j-1] == cnt) f[i][j] = cnt;
				}
				else if (a[i] == b[j] && f[i-1][j-1] == cnt) f[i][j] = cnt;
			} 
		if (f[n][len] == cnt) ++ans;
		else if (a[n] == '*' && f[n-1][len] == cnt) ++ans;
	}
	cout << m - ans;
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值