poj 3691

ac自动机+dp 自动机上的节点来作为状态

dp[i][j]表示长度为i状态为j至少需要转换多少个字符

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const int root=0,N=4;
char a[1111];
int dp[1111][1111],c[1111][10];
struct
{
    int next[5],tmp,fail;
}trie[1111];
int lon;
void trieini()
{
    memset(trie,0,sizeof(trie));
    lon=0;
}

void insert(char s[])
{
    int t=root;
    int n=strlen(s+1);
    for(int i=1;i<=n;i++)
    {
        if(trie[t].next[s[i]-'A']==0)
        trie[t].next[s[i]-'A']=++lon;
        t=trie[t].next[s[i]-'A'];
        if(i==n)
        trie[t].tmp++;
    }
}

void getfail()
{
    queue <int> q;
    q.push(root);
    while(!q.empty())
    {
        int t=q.front();
        q.pop();
        for(int i=0;i<N;i++)
        if(trie[t].next[i])
        {
            int u=trie[t].next[i];
            int tmp=trie[t].fail;
            while(tmp!=root&&trie[tmp].next[i]==0)
            tmp=trie[tmp].fail;
            trie[u].fail=trie[tmp].next[i];
            if(t==root) trie[u].fail=root;
            trie[u].tmp+=trie[trie[u].fail].tmp;
            q.push(u);
        }
    }
}

void find()
{
    for(int i=0;i<=lon;i++)
    {
        for(int k=0;k<N;k++)
        {
            int t=i;
            while(t&&trie[t].next[k]==0) t=trie[t].fail;
            if(trie[t].next[k])
            c[i][k]=trie[t].next[k];
            else
            c[i][k]=root;
        }
    }
}

void work(char a[])
{
    int n=strlen(a+1);
    for(int i=1;i<=n;i++)
    {
        if(a[i]=='T')
        a[i]='B';
        else if(a[i]=='G')
        a[i]='D';
    }
}

int main()
{
    int n;
    int tcase=0;
    while(scanf("%d",&n),n)
    {
        trieini();
        for(int i=1;i<=n;i++)
        {
            char s[22];
            scanf("%s",s+1);
            work(s);
            insert(s);
        }
        getfail();
        find();
        scanf("%s",a+1);
        work(a);
        memset(dp,50,sizeof(dp));
        dp[0][0]=0;
        int m=strlen(a+1);
        for(int k=0;k<m;k++)
        for(int i=0;i<=lon;i++)
        if(trie[i].tmp==0)
        {
            for(int p=0;p<N;p++)
            {
                int tmp=p+'A'!=a[k+1];
                dp[k+1][c[i][p]]=min(dp[k+1][c[i][p]],dp[k][i]+tmp);
            }
        }
        int ans=1111;

        for(int i=0;i<=lon;i++)
        if(trie[i].tmp==0)
        {
            ans=min(ans,dp[m][i]);
//            printf("%d %d %d\n",m,i,dp[m][i]);
        }

        printf("Case %d: ",++tcase);
        if(ans==1111)
        printf("-1\n");
        else
        printf("%d\n",ans);
    }
    return 0;
}