hdu_2457_DNA repair(AC自动机+DP)

题目连接：hdu_2457_DNA repair

题意：

给你N个字符串，最后再给你一个要匹配的串，问你最少修改多少次，使得这个串不出现之前给的N的字符串

题解：

刚学AC自动机，切这题还真不知道怎么来DP，然后看了一下题解，需要在失败指针那里做文章，这里我们要将trie的每一个节点当作一个状态，然后设dp[i][j]表示考虑到第i个字符，j这个trie节点时的最小修改次数，为什么要这样考虑，因为AC自动机在匹配失败的时候会转向其他的节点，所以这里我们要考虑每一个节点的状态，然后当前节点的子节点如果不存在也要处理一下，就指向这个节点的fail指针，这样我们在后面dp的时候才能保证匹配失败的时候回到fail节点

 #include<cstdio>
 #include<cstring>
 #define F(i,a,b) for(int i=a;i<=b;i++)

 inline void up(int &x,int y){if(x>y)x=y;}
 const int AC_N=,inf=1e8;
 struct AC_automation{
     int tr[AC_N][],cnt[AC_N],Q[AC_N],fail[AC_N],tot;
     int gt(char x){
         if(x=='A')return ;
         if(x=='G')return ;
         if(x=='C')return ;
         if(x=='T')return ;
     }
     void nw(){cnt[++tot]=;memset(tr[tot],-,sizeof(tr[tot]));}
     void init(){tot=-,fail[]=-,nw();}
     void insert(char *s,int x=){
         for(int len=strlen(s),i=,w;i<len;x=tr[x][w],i++)
             if(tr[x][w=gt(s[i])]==-)nw(),tr[x][w]=tot;
         cnt[x]=;//串尾标记
     }
     void build(int head=,int tail=){
         for(Q[++tail]=;head<=tail;){
             for(int i=,x=Q[head++],p=-;i<=;i++)if(~tr[x][i]){
                 if(x==)fail[tr[][i]]=;
                 else{
                     for(p=fail[x],fail[tr[x][i]]=;~p;p=fail[p])
                     if(~tr[p][i]){fail[tr[x][i]]=tr[p][i];break;}
                 }//如果他失败指针指向的节点的子节点为危险DNA那么这点的子节点也不能取
                 if(cnt[fail[tr[x][i]]])cnt[tr[x][i]]=;
                 Q[++tail]=tr[x][i];
             }else if(x==)tr[][i]=;//不存在的节点指向失败指针的位置
             else tr[x][i]=tr[fail[x]][i];
         }
     }
     int dp[][];
     int ask(char *s){
         int len=strlen(s),ans=inf;
         F(i,,len)F(j,,tot)dp[i][j]=inf;
         dp[][]=;
         F(i,,len)F(j,,tot){
             if(cnt[j]||dp[i-][j]==inf)continue;
             F(k,,){
                 int nxt=tr[j][k];
                 if(cnt[nxt])continue;
                 up(dp[i][nxt],dp[i-][j]+(gt(s[i-])!=k));
             }
         }
         F(i,,tot)up(ans,dp[len][i]);
         return ans==inf?-:ans;
     }
 }AC;

 char buf[];
 int main(){
     int n,ic=;
     while(~scanf("%d",&n),n){
         AC.init();
         F(i,,n)scanf("%s",buf),AC.insert(buf);
         AC.build();
         scanf("%s",buf);
         printf("Case %d: %d\n",ic++,AC.ask(buf));
     }
     return ;
 }