hdu_3341_Lost's revenge(AC自动机+状态hashDP)

题目链接：hdu_3341_Lost's revenge

题意：

有n个模式串，一个标准串，现在让标准串重组，使得包含最多的模式串，可重叠，问重组后最多包含多少模式串

题解：

显然是AC自动机上的状态压缩DP，不过如果直接开40⁴*500的数组显示开不下，所以这样要将状态hash一下，然后再DP，因为这个字母加起来为40，所以状态数不超过15000，所以可以接受。

 #include<bits/stdc++.h>
 #define F(i,a,b) for(int i=a;i<=b;i++)
 using namespace std;

 const int AC_N=*,tyn=;//数量乘串长,类型数
 struct AC_automation{
     int tr[AC_N][tyn],cnt[AC_N],Q[AC_N],fail[AC_N],tot;
     inline int getid(char x){
         if(x=='A')return ;
         if(x=='C')return ;
         if(x=='G')return ;
         return ;
     }
     void nw(){cnt[++tot]=,fail[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=getid(s[i])])nw(),tr[x][w]=tot;
         cnt[x]++;//串尾标记
     }
     void build(int head=,int tail=){
         for(int i=;i<tyn;i++)if(tr[][i])Q[++tail]=tr[][i];
         while(head<=tail)for(int x=Q[head++],i=;i<tyn;i++)
             if(tr[x][i])fail[tr[x][i]]=tr[fail[x]][i],Q[++tail]=tr[x][i],cnt[tr[x][i]]+=cnt[tr[fail[x]][i]];
             else tr[x][i]=tr[fail[x]][i];
     }

 }AC;

 inline void up(int &a,int b){if(a<b)a=b;}

 char s[];
 int n,len,cnt,ic=,na,nc,ng,nt;
 int hsh[][][][],dp[][];

 void fuck()
 {
     int ed=,ans=;
     na=nc=ng=nt=;
     for(int i=;s[i]!='\0';i++)if(s[i]=='A')na++;
     else if(s[i]=='C')nc++;
     else if(s[i]=='G')ng++;
     else nt++;
     F(i,,na)F(j,,nc)F(k,,ng)F(l,,nt)hsh[i][j][k][l]=ed++;
     memset(dp,-,sizeof(dp)),dp[][]=;
     F(i,,na)F(j,,nc)F(k,,ng)F(l,,nt)if(i+j+k+l!=)
     {
         int x1=hsh[i][j][k][l],x2;
         F(ii,,AC.tot)F(jj,,)
         {
             if(jj==&&i->=)x2=hsh[i-][j][k][l];
             else if(jj==&&j->=)x2=hsh[i][j-][k][l];
             else if(jj==&&k->=)x2=hsh[i][j][k-][l];
             else if(jj==&&l->=)x2=hsh[i][j][k][l-];
             else continue;
             if(dp[x2][ii]==-)continue;
             int now=AC.tr[ii][jj];
             up(dp[x1][now],dp[x2][ii]+AC.cnt[now]);
             up(ans,dp[x1][now]);
         }
     }
     printf("Case %d: %d\n",++ic,ans);
 }

 int main()
 {
     while(scanf("%d",&n),n)
     {
         AC.init();
         F(i,,n)scanf("%s",s),AC.insert(s);
         AC.build(),scanf("%s",s),fuck();
     }
     return ;
 }