P2536 [AHOI2005]病毒检测（记忆化搜索）

题目描述

科学家们在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片段不是病毒。

输入输出格式

输入格式：

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

输出格式：

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

输入输出样例

输入样例#1： 复制

A*G?C
    3
AGTC
AGTGTC
AGTGC

输出样例#1： 复制

1

说明

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

---

---

---

在做这个提之前，你要知道一个事情，就是rna必须从头开始匹配 。。。。。。

然而开始不知道这个，在想题解为什么不建fail  qwq

模板到了*,模板走串不走 或者 模板不走串走

模板到了agut ，能匹配就走

模板到了？ ，都走一步

vis[x][y] 表示模板到了第x位，到了树上的y节点，记忆一波

唯有吸氧才能救这个码 嘿嘿

---

---

---

#include<iostream>
#include<cstring>
#include<cstdio>
#include<bitset>

#include<queue>
using namespace std;
const int N=500+5;
const int L=1000+5;
const int M=N*N;
template<class T>inline void read(T &num){
    char ch;
    while(!isdigit(ch=getchar()));
    num=ch-'0';
    while(isdigit(ch=getchar()))num=num*10+ch-'0';
}

int mmp[300];//愤怒之情溢于数组
int tri[M][4],is_end[M],fla[M],totn=1,n,ans,res[M],valu[L];
char str[L],s[N];
inline void Insert(char s[]){
    int len=strlen(s),x=0;
    for(int p=0;p<len;++p){
        int ch=mmp[s[p]];
        if(tri[x][ch]==0){
            tri[x][ch]=++totn;
        }
        x=tri[x][ch];
    }
    is_end[x]++;
}

struct Data{
    int first,second;
    Data(const int a=0,const int b=0){
        first=a,second=b;
    }
};

queue<Data>que;
bitset<L> vis[M];
inline void Push(const int x,const int y){
    if(vis[x][y])return;
    que.push(Data(x,y));
    vis[x][y]=true;
}

inline void bfs(){
    int len=strlen(str),x,p;
    Push(0,0);
    while(que.size()){
        x=que.front().first,p=que.front().second;que.pop();
        if(p==len){
            if(fla[x]==0){
                ans-=is_end[x];
                fla[x]=1;
            }
            continue;
        }
        int ch=mmp[str[p]];
        if(ch<4&&tri[x][ch])Push(tri[x][ch],p+1);
        else if(ch==4){
            Push(x,p+1);
            for(int i=0;i<4;++i)if(tri[x][i])Push(tri[x][i],p);
        }
        else if(ch==5)for(int i=0;i<4;++i)if(tri[x][i])Push(tri[x][i],p+1);
    }
}

int main(){
    mmp['G']=1,mmp['C']=2,mmp['T']=3,mmp['*']=4,mmp['?']=5;
    scanf("%s",str);
    read(n);
    ans=n;
    for(int i=1;i<=n;++i){
        scanf("%s",s);
        Insert(s);
    }
    bfs();
    printf("%d
",ans);
    return 0;
}