HDU 4057 Rescue the Rabbit ( AC自动机 + 状态压缩DP )

模板来自notonlysuccess.

模式串只有10个，并且重复出现的分值不累加，因此很容易想到状态压缩。

将模式串加入AC自动机，最多有10*100个状态。

dp[i][j][k]：串长为i，在Trie图上的状态为j，已经包含的模式串为k（二进制表示，第x位为1代表已经包含第x个串）。

dp[i][j][k]为true或false代表该状态是否可达。

沿着Trie图中的边走进行DP，时间复杂度O( 100*1000*1024 );

最后枚举一下串长为L的所有可达状态，最大值即为结果。

PS.内存有限，需要用滚动数组

PS2.若使用静态队列的话，队列大小最好开状态数的二倍或者更多

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>

using namespace std;

const int MAX_NODE = ;
const int CHILD_NUM = ;
const int MAXN = ;
const int INF =  << ;

struct ACAutomaton
{
    int chd[MAX_NODE][CHILD_NUM]; //每个节点的儿子,即当前节点的状态转移
    int val[MAX_NODE];            //记录题目给的关键数据
    int fail[MAX_NODE];           //传说中的fail指针
    int Q[MAX_NODE<<];           //队列,用于广度优先计算fail指针
    int ID[];                  //字母对应的ID
    int sz;                       //已使用节点个数

    //初始化,计算字母对应的儿子ID,如:'a'->0 ... 'z'->25
    void Initialize()
    {
        fail[] = ;
        ID['A'] = ;
        ID['G'] = ;
        ID['C'] = ;
        ID['T'] = ;
        return;
    }
    //重新建树需先Reset
    void Reset()
    {
        memset(chd[] ,  , sizeof(chd[]));
        val[] = ;
        sz = ;
    }
    //将权值为key的字符串a插入到trie中
    void Insert(char *a,int key)
    {
        int p = ;
        for ( ; *a ; a ++)
        {
            int c = ID[*a];
            if (!chd[p][c])
            {
                memset(chd[sz] ,  , sizeof(chd[sz]));
                val[sz] = ;
                chd[p][c] = sz ++;
            }
            p = chd[p][c];
        }
        val[p] = key;
    }
    //建立AC自动机,确定每个节点的权值以及状态转移
    void Construct()
    {
        int *s = Q , *e = Q;
        for (int i =  ; i < CHILD_NUM ; i ++)
        {
            if (chd[][i])
            {
                fail[ chd[][i] ] = ;
                *e ++ = chd[][i];
            }
        }
        while (s != e)
        {
            int u = *s++;
            for (int i =  ; i < CHILD_NUM ; i ++)
            {
                int &v = chd[u][i];
                if (v)
                {
                    *e ++ = v;
                    fail[v] = chd[ fail[u] ][i];
                    //以下一行代码要根据题目所给val的含义来写
                    val[v] |= val[ fail[v] ];
                }
                else
                {
                    v = chd[ fail[u] ][i];
                }
            }
        }
    }
} AC;

int N, L;
int cost[MAXN];
bool dp[][MAX_NODE][ ( << ) +  ];

int getVal( int S )
{
    int res = ;
    for ( int j = ; j < N; ++j )
        if ( S & (  << j ) )
            res += cost[j];
    return res;
}

int main()
{
    char gene[];
    AC.Initialize();
    while ( scanf( "%d%d", &N, &L ) ==  )
    {
        AC.Reset();
        for ( int i = ; i < N; ++i )
        {
            scanf( "%s%d", gene, &cost[i] );
            AC.Insert( gene,  << i );
        }
        AC.Construct();

        int all =  << N;

        memset( dp[], false, sizeof(dp[]) );
        int pre = , cur = ;
        dp[][][] = true;
        for ( int i = ; i <= L; ++i )
        {
            memset( dp[cur], , sizeof(dp[cur]) );
            for ( int j = ; j < AC.sz; ++j )
            {
                for ( int k = ; k < ; ++k )
                {
                    int next = AC.chd[j][k];
                    for ( int S = ; S < all; ++S )
                    {
                        if ( dp[pre][j][S] )
                            dp[cur][next][ S|AC.val[next] ] = true;
                    }
                }
            }
            pre ^= ;
            cur ^= ;
        }

        int ans = -;
        for ( int i = ; i < all; ++i )
        {
            for ( int j = ; j < AC.sz; ++j )
            {
                if ( dp[pre][j][i] )
                {
                    ans = max( ans, getVal(i) );
                    break;
                }
            }
        }

        if ( ans <  ) puts("No Rabbit after 2012!");
        else printf( "%d\n", ans );
    }
    return ;
}