HDU2457 DNA repair（AC自动机+DP）

题目一串DNA最少需要修改几个基因使其不包含一些致病DNA片段。

这道题应该是AC自动机+DP的入门题了，有POJ2778基础不难写出来。

dp[i][j]表示原DNA前i位（在AC自动机上转移i步）且后缀状态为AC自动机结点j的最少需要修改的基因数

转移我为人人型，从dp[i][j]向ATCG四个方向转移到dp[i+1][j']，如果结点被标记包含致病基因就不能转移。

 #include<cstdio>
 #include<cstring>
 #include<queue>
 using namespace std;
 #define INF (1<<30)
 int tn,ch[][],fail[],idx[];
 bool flag[];
 void insert(char *s){
     int x=;
     for(int i=; s[i]; ++i){
         int y=idx[s[i]];
         if(ch[x][y]==) ch[x][y]=++tn;
         x=ch[x][y];
     }
     flag[x]=;
 }
 void init(){
     memset(fail,,sizeof(fail));
     queue<int> que;
     for(int i=; i<; ++i){
         if(ch[][i]) que.push(ch[][i]);
     }
     while(!que.empty()){
         int x=que.front(); que.pop();
         for(int i=; i<; ++i){
             if(ch[x][i]) que.push(ch[x][i]),fail[ch[x][i]]=ch[fail[x]][i],flag[ch[x][i]]|=flag[ch[fail[x]][i]];
             else ch[x][i]=ch[fail[x]][i];
         }
     }
 }
 int d[][];
 int main(){
     idx['A']=; idx['G']=; idx['C']=; idx['T']=;
     char str[];
     int n,t=;
     while(~scanf("%d",&n) && n){
         tn=;
         memset(ch,,sizeof(ch));
         memset(flag,,sizeof(flag));
         while(n--){
             scanf("%s",str);
             insert(str);
         }
         init();
         scanf("%s",str+);
         n=strlen(str+);
         for(int i=; i<=n; ++i){
             for(int j=; j<=tn; ++j) d[i][j]=INF;
         }
         d[][]=;
         for(int i=; i<n; ++i){
             for(int j=; j<=tn; ++j){
                 if(d[i][j]==INF || flag[j]) continue;
                 for(int k=; k<; ++k){
                     if(flag[ch[j][k]]) continue;
                     if(idx[str[i+]]==k) d[i+][ch[j][k]]=min(d[i+][ch[j][k]],d[i][j]);
                     else d[i+][ch[j][k]]=min(d[i+][ch[j][k]],d[i][j]+);
                 }
             }
         }
         int res=INF;
         for(int i=; i<=tn; ++i) res=min(res,d[n][i]);
         if(res==INF) printf("Case %d: %d\n",++t,-);
         else printf("Case %d: %d\n",++t,res);
     }
     return ;
 }