HDU4057 Rescue the Rabbit（AC自动机+状压DP）

题目大概是给几个DNA片段以及它们各自的权值，如果一个DNA包含某个片段那么它的价值就加上这个片段的权值，同时包含多个相同DNA片段也只加一次，问长度l的DNA可能的最大价值。

与HDU2825大同小异。

• dp[i][j][S]表示长度i（自动机转移i步）、后缀状态为自动机第j个结点、包含的DNA片段为集合S 的DNA最大价值
• dp[0][0][0]=0
• 我为人人转移，从dp[i][j][S]向ATCG四个方向更新dp[i+1][j'][S']

1. 注意的是要用滚动数组，不然要开上百兆数组。
2. 用滚动数组要注意初始化。
3. 转移过程有些计算可以预处理出来，存在数组里，虽然不预处理应该也不会TLE，题目时间给了10000ms。

#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
#define INF (1<<30)
int tn,ch[1111][4],fail[1111],flag[1111];
int idx[128];
void insert(char *s,int k){
    int x=0;
    for(int i=0; s[i]; ++i){
        int y=idx[s[i]];
        if(ch[x][y]==0) ch[x][y]=++tn;
        x=ch[x][y];
    }
    flag[x]|=1<<k;
}
void init(){
    memset(fail,0,sizeof(fail));
    queue<int> que;
    for(int i=0; i<4; ++i){
        if(ch[0][i]) que.push(ch[0][i]);
    }
    while(!que.empty()){
        int x=que.front(); que.pop();
        for(int i=0; i<4; ++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[2][1111][1<<10],VAL[1<<10];
int main(){
    idx['A']=0; idx['G']=1; idx['T']=2; idx['C']=3;
    int m,n,a;
    char str[111];
    int val[11];
    while(~scanf("%d%d",&m,&n)){
        tn=0;
        memset(ch,0,sizeof(ch));
        memset(flag,0,sizeof(flag));
        for(int i=0; i<m; ++i){
            scanf("%s%d",str,val+i);
            if(strlen(str)>100) continue;
            insert(str,i);
        }
        for(int i=1; i<(1<<m); ++i){
            for(int j=0; j<m; ++j){
                if((i>>j)&1) VAL[i]=VAL[i^(1<<j)]+val[j];
            }
        }
        init();
        for(int j=0; j<=tn; ++j){
            for(int k=0; k<(1<<m); ++k) d[0][j][k]=-INF;
        }
        d[0][0][0]=0;
        for(int i=0; i<n; ++i){
            int x=i&1;
            for(int j=0; j<=tn; ++j){
                for(int k=0; k<(1<<m); ++k) d[x^1][j][k]=-INF;
            }
            for(int j=0; j<=tn; ++j){
                for(int k=0; k<(1<<m); ++k){
                    if(d[x][j][k]==-INF) continue;
                    for(int y=0; y<4; ++y){
                        d[x^1][ch[j][y]][k|flag[ch[j][y]]]=max(d[x^1][ch[j][y]][k|flag[ch[j][y]]],d[x][j][k]+VAL[k^(k|flag[ch[j][y]])]);
                    }
                }
            }
        }
        int res=-INF;
        for(int i=0; i<=tn; ++i){
            for(int j=0; j<(1<<m); ++j) res=max(res,d[n&1][i][j]);
        }
        if(res<0) puts("No Rabbit after 2012!");
        else printf("%d
",res);
    }
    return 0;
}