poj 3691

AC自动机的失配函数的利用和DP，Trie图中的每个节点代表一种状态，即一个字符串的后缀为从根节点到本节点的字符串。就是说把字符串根据他的后缀分为不同的类别，AC自动机上的每一个节点代表一个类别。在一个字符串的后面再加一个字符，此字符串就会从Tire图中的一个节点转移到另一个节点，这个过程可通过失配函数的求解过程中求出。其中有的节点所代表的类是不合法的，要记录下来。然后就是dp,ans[i][j]表示母串的i长度的串要变成状态为j的串所需最小修改次数。初始化为无穷大，如果j不合法不求，并且j的下一个转移得到的节点也不从j转移得到，最后比较每个节点的ans[len][j],得到最小值。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const int maxn=50+10;
const int inf=0x3f3f3f3f;
int n;
int ch[maxn*20][4];
int val[maxn*20],tot,f[maxn*20];
int ans[1000+10][maxn*20];
int getV(char c)
{
    switch(c)
    {
    case 'A':return 0;
    case 'G':return 1;
    case 'C':return 2;
    case 'T':return 3;
    }
}
void Insert(char *s)
{
    int u=0,len=strlen(s),i,tv;
    for(i=0;i<len;i++)
    {
        tv=getV(s[i]);
        if(!ch[u][tv]) ch[u][tv]=tot++;
        u=ch[u][tv];
    }
    val[u]=1;
}
void getFail()
{
    queue<int> q;
    f[0]=0;
    int i,u;
    for(i=0;i<4;i++)
    {
        u=ch[0][i];
        if(u)
        {
            f[u]=0;
            q.push(u);
            //printf("%d %d\n",u,val[u]);
        }
    }
    int v,r;
    while(!q.empty())
    {
        r=q.front();q.pop();
        for(i=0;i<4;i++)
        {
            u=ch[r][i];
            if(!u) { ch[r][i]=ch[f[r]][i]; }
            else
            {
            if(val[r]) val[u]=1;//这部很关键
            q.push(u);
            v=f[r];
            while(v&&!ch[v][i]) v=f[v];
            f[u]=ch[v][i];
            if(val[f[u]]) val[u]=1;//这步很关键
            }
         //   printf("%d %d\n",u,val[u]);
        }
    }
}
int main()
{
    int cas=0;
    while(scanf("%d",&n)&&n)
    {
        int i,j;
        char s[1000+10];
        tot=1;
        memset(ch,0,sizeof(ch));
        memset(val,0,sizeof(val));
        for(i=0;i<n;i++)
        {
            scanf("%s",s);
            Insert(s);
        }
        scanf("%s",s);
        getFail();
        int len=strlen(s);
        for(i=0;i<=len;i++)
        {
            for(j=0;j<tot;j++)
                ans[i][j]=inf;
        }
        ans[0][0]=0;
        int k,v;
        for(i=0;i<len;i++)
        {
            for(j=0;j<tot;j++)
            {
                if(ans[i][j]<inf)
                {
                for(k=0;k<4;k++)
                {
                    if(!val[ch[j][k]])
                    {
                        v=ch[j][k];
                        ans[i+1][v]=min(ans[i+1][v],ans[i][j]+(getV(s[i])==k?0:1));
                    }
                }
            }
            }
        }
        int re=inf;
        for(i=0;i<tot;i++) re=min(re,ans[len][i]);
        printf("Case %d: %d\n",++cas,re<inf?re:-1);
    }
    return 0;
}