HDU 2457 DNA repair

题目链接：HDU-2457

题意：给出一些不合法的模式DNA串，给出一个原串，问最少需要修改多少个字符，使得原串中不包含非法串。

思路：沿着AC自动机dp，dp[i][j]表示前i个字符，当前为状态j的时候，需要修改的最少字符数。

代码：

#include<cstdio>
#include<cstdlib>
#include<vector>
#include<set>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<queue>
using namespace std;
typedef long long LL;
const LL MAXN=2e3+10;
const LL INF=0x3f3f3f3f;

struct AC
{
    int ch[MAXN][4];
    int val[MAXN];
    int chr[MAXN];
    int sz;
    AC() { init(); }
    void init() { sz=1; memset(ch[0],0,sizeof(ch[0])); memset(val,0,sizeof(val));}
    int idx(char c)
    {
        if(c=='A') return 0;
        if(c=='G') return 1;
        if(c=='C') return 2;
        if(c=='T') return 3;
        return 4;
    }
    void insert(char *s)
    {
        int u=0,n=strlen(s);
        for(int i=0;i<n;i++)
        {
            int c=idx(s[i]);
            if(!ch[u][c])
            {
                memset(ch[sz],0,sizeof(ch[sz]));
                val[sz]=0;
                ch[u][c]=sz++;
                chr[sz-1]=c;
            }
            u=ch[u][c];
        }
        val[u]++;
    }
    int f[MAXN];
    int last[MAXN];
    int find(char *T)
    {
        int ans=0;
        int n=strlen(T);
        int j=0;
        for(int i=0;i<n;i++)
        {
            int c=idx(T[i]);
            j=ch[j][c];
            if(val[j])
            {
                ans+=val[j];
                val[j]=0;
            }
            int k=j;
            while(last[k])
            {
                k=last[k];
                ans+=val[k];
                val[k]=0;
            }
        }
        return ans;
    }
    void getFail()
    {
        queue<int> q;
        f[0]=0;
        for(int c=0;c<4;c++)
        {
            int u=ch[0][c];
            if(u) {f[u]=0;q.push(u);last[u]=0;}
        }
        while(!q.empty())
        {
            int r=q.front();q.pop();
            for(int c=0;c<4;c++)
            {
                int u=ch[r][c];
                if(!u)
                {
                    ch[r][c]=ch[f[r]][c];
                    continue;
                }
                q.push(u);
                int v=f[r];
                while(v && !ch[v][c]) v=f[v];
                f[u]=ch[v][c];
                last[u]=val[f[u]]?f[u]:last[f[u]];
                if(!val[u]) val[u]=val[last[u]];
            }
        }
    }
};
AC T;
int f[MAXN][MAXN];
char s[MAXN];
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
    // freopen("out.txt","w",stdout);
#endif
    int n;
    int tt=0;
    while(scanf("%d",&n)!=EOF && n)
    {
        tt++;
        T.init();
        for(int i=1;i<=n;i++)
        {
            scanf("%s",s);
            T.insert(s);
        }
        T.getFail();
        scanf("%s",s);
        int l=strlen(s);
        memset(f,INF,sizeof(f));
        for(int k=0;k<4;k++) if(!T.val[T.ch[0][k]])
            f[0][T.ch[0][k]]=(k==T.idx(s[0]))?0:1;
        for(int i=0;i<l-1;i++)
            for(int j=0;j<T.sz;j++)
                if(!T.val[j] && f[i][j]!=INF)
                    for(int k=0;k<4;k++)
                        if(!T.val[T.ch[j][k]])
                            f[i+1][T.ch[j][k]]=min(f[i][j]+((k==T.idx(s[i+1]))?0:1),f[i+1][T.ch[j][k]]);
        int ans=INF;
        for(int j=0;j<T.sz;j++)
            if(!T.val[j])
                ans=min(ans,f[l-1][j]);
        if(ans>=INF) ans=-1;
        printf("Case %d: %d
",tt,ans);
    }
    return 0;
}