HDU3341AC自动机+dp

题意：N个子串，每个串价值为1，给你一个长串S，求重排后最大的价值（可重叠）；

明白题意后，很裸的自动机+DP，建完自动机状态表示成dp[2][500][20][20][20]，把其中3个必然小于20的元素来DP，然后TLE，40*500*20*20*20，有30组case，必然TLE；

/*
学习一：一种压缩方式，类似于进制表示，每个元素的个数确定后
num[],那么bas[0]=1,bas[1]=bas[i-1]*(num[i-1]+1);
例如num[]={2,1,3,2};bas[]={1,3,6,24};
编码 incode() ret+=num[i]*bas[i];
然后就是解码 decode();
这样把20^3*40的空间压缩成11^4的空间；
但是还是TLE；40*500*4*11^4*30;
我们发现状态的转移都是从长度为n-1转移到n;
而长度为n-1的编码一定是小于长度为n的编码的；
类似于1100->1101或1110数值都是变大的；
所以直接从编码小的转移到编码大的就可以了，
500*11^4*4*30  AC;
*/
//18:34
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include<map>
using namespace std;
const int N=500+10;
const int SZ=4;
int idx(char s){
    if (s=='A') return 0;
    if (s=='G') return 1;
    if (s=='T') return 2;
    return 3;
}
struct AC{
    int ch[N][SZ];
    int sz,f[N],match[N];
    queue<int> q;
    void init(){
        sz=1;memset(match,0,sizeof(match));
        memset(ch[0],0,sizeof(ch[0]));
    }
    void insert(char *s){
        int m=strlen(s),rt=0;
        for (int i=0;i<m;i++){
            int c=idx(s[i]);
            if (ch[rt][c]==0){
                ch[rt][c]=sz;
                memset(ch[sz],0,sizeof(ch[sz]));
                sz++;
            }
            rt=ch[rt][c];
        }
        match[rt]+=1;
    }
    void getfail(){
        while (!q.empty()) q.pop();
        for (int i=0;i<SZ;i++){
            int c=ch[0][i];
            if (c){ f[c]=0;q.push(c); }
        }
        while(!q.empty()){
            int u=q.front();q.pop();
            for (int i=0;i<SZ;i++){
                int r=ch[u][i];
                if(!r){
                    ch[u][i]=ch[f[u]][i];
                }else {
                    q.push(r);
                    int v=f[u];
                    f[r]=ch[v][i];
                    match[r]+=match[f[r]];
                }
            }
        }
    }
    void find(){
        cout<<"----------- "<<endl;
        for (int i=0;i<sz;i++){
            if (match[i]) cout<<i<<" "<<match[i]<<endl;
        }
        cout<<"+++++++++++ "<<endl;
    }
}H;
/////////////////////AC自动机部分///////////////////////
char T[50];
int up[10];
int bas[10];
void init(){
    scanf("%s",T);
    memset(up,0,sizeof(up));
    for (int i=0;T[i];i++){
        int c=idx(T[i]);
        up[c]++;
    }
    bas[0]=1;
    for (int i=1;i<=4;i++){
        bas[i]=bas[i-1]*(up[i-1]+1);
    }
  /*  cout<<"*** "<<endl;
    for (int i=0;i<=4;i++){
        cout<<bas[i]<<" ";
    }cout<<endl;
    cout<<"/// "<<endl;
*/
}
int dp[15000][N];
void decode(int *num,int k){
    for (int i=3;i>=0;i--){
        num[i]=k/bas[i];
        k=k%bas[i];
    }
}
void Max(int &x,int y){
    if (x==-1) {
        x=y; 
    }else {
        x=max(x,y);
    }
}
void DP(){
    memset(dp,-1,sizeof(dp));
    dp[0][0]=0;
    int num[10];
    for (int i=0;i<bas[4];i++){
       decode(num,i);
       for (int j=0;j<H.sz;j++){
            if (dp[i][j]!=-1){
                for (int k=0;k<SZ;k++){
                    int c=H.ch[j][k];
                    if (num[k]<up[k]){
                        int t=i+bas[k];
                        Max(dp[t][c],dp[i][j]+H.match[c]);
                    }
                }
            }
        }
    }

    int ret=-1,all=bas[4]-1;
    for (int i=0;i<H.sz;i++){
        if (dp[all][i]!=-1) Max(ret,dp[all][i]);
    }
    printf("%d\n",ret);
}
int n;
char p[15];
int main(){
    int cas=0;
    while (~scanf("%d",&n),n){
        H.init();
        for (int i=0;i<n;i++){
            scanf("%s",p);
            H.insert(p);
        }
        H.getfail();
        //
       // H.find();
        //
        init();
        printf("Case %d: ",++cas);
        DP();
    }
    return 0;
}