给m个单词,由这m个单词组成的一个新单词(两个单词可以重叠包含)长度为n,且新单词中包含的基本单词数目不少于k个。问这样的新单词共有多少个?
m很小,用二进制表示新单词中包含基本单词的情况。
用m个单词建立AC自动机,可以求出所有单词之间相互包含的情况,AC自动机的后缀特性(每个结点的失配边指向新结点,新节点到trie树根的字符串是当前节点字符串的后缀)。
动态规划:
f(i, j, k):长度为i的单词,在trie树中第j个结点处,包含基本单词的情况为k(二进制),的总方案数。
转移方程:
在当前字符串后边再加上一个字母
f(i+1, u, k|val[u]) += f(i, j, k)
u是j结点的子结点(或者是失配边指向的结点),k|val[u]表示加上一个新字母后包含基本单词的情况。
#include <cstdio> #include <cstring> #include <algorithm> #include <queue> using namespace std; struct AC_Automata { #define N 102 #define M 26 int ch[N][M], val[N], last[N], f[N], sz; void clear() { sz = 1; memset(ch[0], 0, sizeof(ch[0])); } int idx(char c) { return c-'a'; } void insert(char s[], int v) { int u = 0; for (int i=0; s[i]; i++) { int c = idx(s[i]); if (!ch[u][c]) { memset(ch[sz], 0, sizeof(ch[sz])); val[sz] = 0; ch[u][c] = sz++; } u = ch[u][c]; } val[u] = 1<<v; } void build() { queue<int> q; f[0] = 0; for (int c=0; c<M; c++) { int u = ch[0][c]; if (u) { f[u] = last[u] = 0; q.push(u); } } while (!q.empty()) { int r = q.front(); q.pop(); for (int c=0; c<M; c++) { int u = ch[r][c]; if (!u) { ch[r][c] = ch[f[r]][c]; val[r] = val[r] | val[f[r]]; //从根到当前结点组成的字符串后缀包含基本单词的情况(传递) continue; } q.push(u); f[u] = ch[f[r]][c]; last[u] = val[f[u]] ? f[u] : last[f[u]]; } } } } ac; int c[1030]; ///记录每个状态含有几个串 int f[27][102][1025]; void init() { memset(c, 0, sizeof(c)); for (int i=0; i<1024; i++) { for (int j=0; j<10; j++) if ((1<<j) & i) c[i]++; } } void solve(int n, int m, int p) { #define Mod 20090717 int u; for (int i=0; i<=n; i++) for (int j=0; j<ac.sz; j++) for (int k=0; k<(1<<m); k++) f[i][j][k] = 0; f[0][0][0] = 1; for (int i=0; i<n; i++) for (int j=0; j<ac.sz; j++) { for (int k=0; k<(1<<m); k++) { if (f[i][j][k] == 0) continue; for (int jj=0; jj<M; jj++) { u = ac.ch[j][jj]; f[i+1][u][ac.val[u]|k] = (f[i+1][u][ac.val[u]|k] + f[i][j][k]) % Mod; } } } int ans = 0; for (int i=0; i<ac.sz; i++) for (int j=0; j<(1<<m); j++) if (c[j] >= p) ans = (ans + f[n][i][j]) % Mod; printf("%d\n", ans); } int main() { int n, m, k; char s[12]; init(); while (scanf("%d%d%d", &n, &m, &k) == 3) { if (n == 0 && m == 0 && k == 0) break; ac.clear(); for (int i=0; i<m; i++) { scanf(" %s", s); ac.insert(s, i); } ac.build(); solve(n, m, k); } return 0; }