HDU 2457 DNA repair (AC自动机+DP)

题意:给N个串，一个大串，要求在最小的改变代价下，得到一个不含上述n个串的大串。

思路:dp，f[i][j]代表大串中第i位，AC自动机上第j位的最小代价。

 #include<algorithm>
 #include<cstdio>
 #include<cmath>
 #include<cstring>
 #include<iostream>
 #include<queue>
 std::queue<int>Q;
 int ch[][],end[],f[][],sz,root,fail[],n,Tcase;
 char s[];
 int newnode(){
     for (int i=;i<=;i++)
      ch[sz][i]=-;
     end[sz]=;
     sz++;
     return sz-;
 }
 void clear(){
     sz=;root=newnode();
 }
 int num(char c){
     if (c=='A') return ;
     else
     if (c=='T') return ;
     else
     if (c=='G') return ;
     else
     if (c=='C') return ;
     return ;
 }
 void insert(){
     scanf("%s",s);
     int len=strlen(s);
     int now=root;
     for (int i=;i<len;i++){
         if (ch[now][num(s[i])]==-) ch[now][num(s[i])]=newnode();
         now=ch[now][num(s[i])];
     }
     end[now]=;
 }
 void build(){
     int now=root;
     std::queue<int>Q;
     fail[root]=root;
     for (int i=;i<=;i++)
      if (ch[root][i]==-) ch[root][i]=root;
      else fail[ch[root][i]]=root,Q.push(ch[root][i]);
     while (!Q.empty()){
         int now=Q.front();Q.pop();
         if (end[fail[now]]==) end[now]|=;
         for (int i=;i<=;i++){
             if (ch[now][i]==-) ch[now][i]=ch[fail[now]][i];
             else fail[ch[now][i]]=ch[fail[now]][i],Q.push(ch[now][i]);
         }
     }
 }
 void dp(){
     scanf("%s",s);
     int len=strlen(s);
     for (int i=;i<=len;i++)
      for (int j=;j<=sz;j++)
       f[i][j]=;
     f[][root]=;
     for (int i=;i<len;i++)
      for (int j=;j<sz;j++)
       if (f[i][j]<){
             for (int k=;k<=;k++){
                 int now=ch[j][k];
                 if (end[now]) continue;
                 int tmp;
                 if (k==num(s[i])) tmp=f[i][j];
                 else tmp=f[i][j]+;
                 f[i+][now]=std::min(f[i+][now],tmp);
             }
       }
     int ans=;
     for (int i=;i<sz;i++)
      ans=std::min(ans,f[len][i]);
     if (ans==) ans=-;
     printf("%d\n",ans);
 }
 int main(){
     while (scanf("%d",&n)!=EOF){
         if (n==) return ;
         Tcase++;
         printf("Case %d: ",Tcase);
         clear();
         for (int i=;i<=n;i++)
          insert();
         build();
         dp();
     }
 }