HDU3341 Lost's revenge（AC自动机+DP）

题目是给一个DNA重新排列使其包含最多的数论基因。

考虑到内存大概就只能这么表示状态：

dp[i][A][C][G][T]，表示包含各碱基个数为ACGT且当前后缀状态为自动机第i的结点的字符串最多的数论基因数

其中ACGT可以hash成一个整数（a*C*G*T+c*G*T+g*T+T），这样用二维数组就行了，而第二维最多也就11*11*11*11个。

接下来转移依然是我为人人型，我是丢进一个队列，用队列来更新状态的值。

这题果然挺卡常数的，只好手写队列，最后4500msAC，还是差点超时，代码也搞得挺乱的。

 #include<cstdio>
 #include<cstring>
 using namespace std;
 int tn,ch[][],fail[],flag[];
 int idx[];
 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];
 }
 int que[];
 void init(){
     memset(fail,,sizeof(fail));
     int front=,rear=;
     for(int i=; i<; ++i){
         if(ch[][i]) que[rear++]=ch[][i];
     }
     while(front!=rear){
         int x=que[front++];
         for(int i=; i<; ++i){
             if(ch[x][i]) que[rear++]=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 quex[*],quey[*];
 int main(){
     idx['A']=; idx['C']=; idx['G']=; idx['T']=;
     char str[];
     int n,cse=;
     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);
         int times[]={};
         for(int i=; str[i]; ++i){
             ++times[idx[str[i]]];
         }
         int tcal0=(times[]+)*(times[]+)*(times[]+);
         int tcal1=(times[]+)*(times[]+);
         int tcal2=(times[]+);
         memset(d,-,sizeof(d));
         d[][]=;
         int front=,rear=,x,y,ny,cnt[];
         quex[rear]=; quey[rear]=; ++rear;
         while(front!=rear){
             x=quex[front]; y=quey[front]; ++front;
             ny=y;
             cnt[]=ny/tcal0;
             ny-=cnt[]*tcal0;
             cnt[]=ny/tcal1;
             ny-=cnt[]*tcal1;
             cnt[]=ny/tcal2;
             cnt[]=ny-cnt[]*tcal2;
             for(int i=; i<; ++i){
                 if(cnt[i]>=times[i]) continue;
                 ++cnt[i];
                 ny=cnt[]*tcal0+cnt[]*tcal1+cnt[]*tcal2+cnt[];
                 if(d[ch[x][i]][ny]==-){
                     d[ch[x][i]][ny]=d[x][y]+flag[ch[x][i]];
                     quex[rear]=ch[x][i]; quey[rear]=ny; ++rear;
                 }else if(d[ch[x][i]][ny]<d[x][y]+flag[ch[x][i]]){
                     d[ch[x][i]][ny]=d[x][y]+flag[ch[x][i]];
                 }
                 --cnt[i];
             }
         }
         y=times[]*tcal0+times[]*tcal1+times[]*tcal2+times[];
         int res=;
         for(int x=; x<=tn; ++x){
             if(res<d[x][y]) res=d[x][y];
         }
         printf("Case %d: %d\n",++cse,res);
     }
     return ;
 }