HDU 2825 Wireless Password （AC自动机 + 状态压缩DP）

题目链接：Wireless Password

解析：给 m 个单词构成的集合，统计所有长度为 n 的串中，包含至少 k 个单词的方案数。

AC自动机 + 状态压缩DP。

DP[i][j][k]：长度为i的字符串匹配到状态j且包含k个magic word的可能字符串个数。

AC代码:

#include <algorithm>
#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>
using namespace std;

const int MOD = 20090717;
int n, m, k;
int dp[30][110][1<<10];
int num[5000];

struct Trie{
    int next[110][26], fail[110], end[110];
    int root, L;
    int newnode(){
        for(int i=0; i<26; i++) next[L][i] = -1;
        end[L++] = 0;
        return L-1;
    }
    void init(){
        L = 0;
        root = newnode();
    }
    void insert(char buf[], int id){
        int len = strlen(buf);
        int now = root;
        for(int i=0; i<len; i++){
            if(next[now][buf[i] - 'a'] == -1)
                next[now][buf[i] - 'a'] = newnode();
            now = next[now][buf[i] - 'a'];
        }
        end[now] |= (1<<id);             //注意
    }
    void build(){
        queue<int> Q;
        fail[root] = root;
        for(int i=0; i<26; i++){
            if(next[root][i] == -1)
                next[root][i] = root;
            else{
                fail[next[root][i]] = root;
                Q.push(next[root][i]);
            }
        }
        while(!Q.empty()){
            int now = Q.front();
            Q.pop();
            end[now] |= end[fail[now]];
            for(int i=0; i<26; i++){
                if(next[now][i] == -1)
                    next[now][i] = next[fail[now]][i];
                else{
                    fail[next[now][i]] = next[fail[now]][i];
                    Q.push(next[now][i]);
                }
            }
        }
    }
    int solve(){
//        memset(dp, 0, sizeof(dp));          //memset初始化DP数组，比下面的慢
        for(int i=0; i<=n; i++)
            for(int j=0; j<L; j++)
                for(int p=0; p<(1<<m); p++)
                    dp[i][j][p] = 0;
        dp[0][0][0] = 1;
        for(int i=0; i<n; i++)
            for(int j=0; j<L; j++)
                for(int p=0; p<(1<<m); p++)
                    if(dp[i][j][p] > 0){
                        for(int x=0; x<26; x++){
                            int newi = i+1;
                            int newj = next[j][x];
                            int newp = (p | end[newj]);
                            dp[newi][newj][newp] += dp[i][j][p];
                            dp[newi][newj][newp] %= MOD;
                        }
                    }
        int ans = 0;
        for(int p=0; p<(1<<m); p++){
            if(num[p] < k) continue;
            for(int i=0; i<L; i++)
                ans = (ans + dp[n][i][p]) % MOD;
        }
        return ans;
    }
};

char buf[20];
Trie ac;
int main(){
//    freopen("in.txt", "r", stdin);

    for(int i=0; i<(1<<10); i++){
        num[i] = 0;
        for(int j=0; j<10; j++)
            if(i & (1<<j)) num[i] ++;
    }
    while(scanf("%d%d%d", &n, &m, &k) == 3){
        if(n == 0 && m == 0 && k == 0) break;
        ac.init();
        for(int i=0; i<m; i++){
            scanf("%s", buf);
            ac.insert(buf, i);
        }
        ac.build();
        printf("%d\n", ac.solve());
    }
    return 0;
}