POJ1699 Best Sequence（AC自动机+状压DP）

题目，求包含所有的给定的n个DNA片段的序列的最短长度。

AC自动机上的DP题。

• dp[S][u]表示已经包含的DNA片段集合为S，且当前后缀状态是自动机第u个结点的最短长度
• dp[0][0]=0
• 我为人人+队列轻松转移。。

 #include<cstdio>
 #include<cstring>
 #include<queue>
 #include<algorithm>
 using namespace std;
 #define INF (1<<30)
 int tn,ch[][],fail[],flag[];
 int idx[];
 void insert(char *s,int k){
     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]|=(<<k);
 }
 void getfail(){
     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];
         }
     }
 }
 struct Node{
     int s,x;
     Node(int _s,int _x):s(_s),x(_x){}
 };
 int d[<<][],vis[<<][];
 int main(){
     idx['A']=; idx['C']=; idx['G']=; idx['T']=;
     int t,n;
     char str[];
     scanf("%d",&t);
     while(t--){
         tn=;
         memset(ch,,sizeof(ch));
         memset(flag,,sizeof(flag));
         scanf("%d",&n);
         for(int i=; i<n; ++i){
             scanf("%s",str);
             insert(str,i);
         }
         getfail();
         for(int i=; i<(<<n); ++i){
             for(int j=; j<=tn; ++j) d[i][j]=INF;
         }
         d[][]=;
         memset(vis,,sizeof(vis));
         vis[][]=;
         queue<Node> que;
         que.push(Node(,));
         while(!que.empty()){
             int s=que.front().s,x=que.front().x; que.pop();
             for(int i=; i<; ++i){
                 int ns=s|flag[ch[x][i]],nx=ch[x][i];
                 if(d[ns][nx]>d[s][x]+){
                     d[ns][nx]=d[s][x]+;
                     if(!vis[ns][nx]){
                         vis[ns][nx]=;
                         que.push(Node(ns,nx));
                     }
                 }
             }
             vis[s][x]=;
         }
         int res=INF;
         for(int i=; i<=tn; ++i){
             res=min(res,d[(<<n)-][i]);
         }
         printf("%d\n",res);
     }
     return ;
 }