【BZOJ1030】[JSOI2007]文本生成器【AC自动机】【计数DP】

http://www.lydsy.com/JudgeOnline/problem.php?id=1030

AC自动机经典题型,数据不大,不用矩阵快速幂优化。

/* Footprints In The Blood Soaked Snow */
#include <cstdio>
#include <cstring>
#define cls(a, x) memset(a, x, sizeof(a))

using namespace std;

const int maxn = 6005, maxm = 105, p = 10007, maxq = 10000;

int dp[maxm][maxn], q[maxq];

struct _acm {
	int son[maxn][26], fail[maxn], acmcnt;
	bool flag[maxn];

	void init() {
		cls(son, 0); for(int i = 0; i < maxn; i++) fail[i] = flag[i] = 0;
		acmcnt = 0;
	}

	void insert(char *s) {
		int now = 0, len = strlen(s);
		for(int i = 0; i < len; i++) {
			int &pos = son[now][s[i] - 'A'];
			if(!pos) pos = ++acmcnt;
			now = pos;
		}
		flag[now] = 1;
	}

	void getfail() {
		int h = 0, t = 0;
		for(int i = 0; i < 26; i++) if(son[0][i]) q[t++] = son[0][i];
		while(h != t) {
			int u = q[h++];
			for(int i = 0; i < 26; i++)
				if(!son[u][i]) son[u][i] = son[fail[u]][i];
				else {
					fail[q[t++] = son[u][i]] = son[fail[u]][i];
					flag[son[u][i]] |= flag[fail[son[u][i]]];
				}
		}
	}
} acm;

inline int qpow(int a, int n) {
	int ans = 1;
	for(int t = a; n; n >>= 1, t = t * t % p) if(n & 1) ans = ans * t % p;
	return ans;
}

char str[maxm];

int main() {
	int n, m; scanf("%d%d", &n, &m);
	for(int i = 1; i <= n; i++) {
		scanf("%s", str);
		acm.insert(str);
	}
	acm.getfail();

	dp[0][0] = 1;
	for(int i = 1; i <= m; i++) for(int j = 0; j <= acm.acmcnt; j++) for(int k = 0; k < 26; k++)
		if(!acm.flag[acm.son[j][k]]) dp[i][acm.son[j][k]] = (dp[i][acm.son[j][k]] + dp[i - 1][j]) % p;

	int ans = qpow(26, m);
	for(int i = 0; i <= acm.acmcnt; i++) ans = (ans - dp[m][i] + p) % p;
	printf("%d\n", ans);
	return 0;
}
 


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值