POJ 3691 TRIE图+DP

题意：

已知一个DNA串和一些病毒DNA序列，求出最少改变DNA串中多少个字符，能使得串中不包含任意一个病毒序列。

题解：

嗯，和上一个trie图一样，把病毒建成trie，只要母串不能再trie上走到危险节点即可（危险节点就是病毒dna序列的终止的那个被标有fg的节点）。

然后trie图上的dp，要走到危险节点的时候，枚举转移即可。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>

#define M 1111
#define N 111
#define INF 0x3f3f3f3f

using namespace std;

struct TR
{
    int son[4]; int f; bool fg;
}tr[N*N];

int map[N],cnt,q[N*N],n,m,cas;
char str[N],a[M];
int dp[M][N*N];

inline void insert(char *s)
{
    int len=strlen(s+1);
    int now=1;
    for(int i=1;i<=len;i++)
    {
        if(!tr[now].son[map[s[i]]]) tr[now].son[map[s[i]]]=++cnt;
        now=tr[now].son[map[s[i]]];
    }
    tr[now].fg=true;
}

inline void build()
{
    int h=1,t=2,sta,now,tmp;
    q[1]=1;
    while(h<t)
    {
        sta=q[h++];
        for(int i=0;i<4;i++)
        {
            now=tr[sta].son[i];
            if(sta==1) tmp=1;
            else tmp=tr[tr[sta].f].son[i];
            if(now==0) tr[sta].son[i]=tmp;
            else
            {
                tr[now].f=tmp; tr[now].fg|=tr[tmp].fg;
                q[t++]=now;
            }
        }
    }
}

inline void read()
{
    memset(tr,0,sizeof tr);
    tr[cnt=1].f=1;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",str+1);
        insert(str);
    }
    scanf("%s",a+1);
    m=strlen(a+1);
    build();
}

inline void go()
{
    memset(dp,0x3f,sizeof dp);
    dp[0][1]=0;
    for(int i=0;i<m;i++)
        for(int j=1;j<=cnt;j++)
        {
            for(int k=0;k<4;k++)
                if(!tr[tr[j].son[k]].fg) dp[i+1][tr[j].son[k]]=min(dp[i+1][tr[j].son[k]],dp[i][j]+1);
            if(!tr[tr[j].son[map[a[i+1]]]].fg) dp[i+1][tr[j].son[map[a[i+1]]]]=min(dp[i+1][tr[j].son[map[a[i+1]]]],dp[i][j]);
        }
    int ans=INF;
    for(int i=1;i<=cnt;i++) ans=min(ans,dp[m][i]);
    if(ans==INF) ans=-1;
    printf("Case %d: %d\n",++cas,ans);
}

int main()
{
    map['A']=0; map['G']=1; map['C']=2; map['T']=3;
    while(scanf("%d",&n),n) read(),go();
    return 0;
}