[JSOI2007]文本生成器

本文介绍使用AC自动机解决字符串匹配问题的方法,并通过一道具体题目详细阐述如何计算不包含特定模式的字符串数量。文章提供了完整的C++代码实现。

这道题是很裸的trie图(AC自动机)

王赟的论文里面讲过一道题,求长度为n的字符串中不出现给定字符串的方案数,跟这道题差不多

这道题dp方程多加一维,f[i][j][k]表示长度为i,trie图位置在j,是否经过了至少一个匹配串的方案数

trie图很久没写了……今天复习了一下……


用的是未改进的trie图,因为字符集比较小……


PS:吐槽这题的数据很水……我在trie图上加边的过程一开始写错了也能AC……


//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i)
#define drep(i,a,b) for(int i=a,tt=b;i>=tt;--i)
#define erep(i,e,x) for(int i=x;i;i=e[i].next)
#define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define read() (strtol(ipos,&ipos,10))
#define sqr(x) ((x)*(x))
#define pb push_back
#define PS system("pause");
typedef long long ll;
typedef pair<int,int> pii;
const int oo=~0U>>1;
const double inf=1e20;
const double eps=1e-6;
string name="",in=".in",out=".out";
//Var
int n,m;
int fa[6008],q[6008],dan[6008],e[6008][27],ans,mod=10007,tot=1;
int f[108][6008][2];
void Init()
{
	scanf("%d%d\n",&n,&m);char s[108];
	rep(i,1,n)
	{
		gets(s+1);int k=1;
		rep(j,1,strlen(s+1))
		{
			s[j]-='A'-1;
			if(!e[k][s[j]])e[k][s[j]]=++tot;
			k=e[k][s[j]];
		}
		dan[k]=true;
	}
}
void Work()
{
	int x;q[1]=1;
	for(int l=1,r=1;l<=r;l++)
	{
		x=fa[q[l]];
		dan[q[l]]=dan[q[l]]||(x&&dan[x]);
		rep(i,1,26)
			if(e[q[l]][i])
			{
				x=fa[q[l]];while(x&&!e[x][i])x=fa[x];
				if(x&&e[x][i])fa[e[q[l]][i]]=e[x][i];else fa[e[q[l]][i]]=1;
				q[++r]=e[q[l]][i];
			}
			else
			{
				x=fa[q[l]];while(x&&!e[x][i])x=fa[x];
				if(x&&e[x][i])e[q[l]][i]=e[x][i];else e[q[l]][i]=1;
			}
	}
	f[0][1][0]=1;
	rep(i,1,m)rep(j,1,tot)rep(k,1,26)rep(l,0,1)
	{
		if(dan[e[j][k]])x=1;else x=l;
		f[i][e[j][k]][x]=(f[i][e[j][k]][x]+f[i-1][j][l])%mod;
	}
	rep(i,1,tot)ans=(ans+f[m][i][1])%mod;
	printf("%d\n",ans);
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
	return 0;
}


# P4052 [JSOI2007] 文本生成器 ## 题目描述 JSOI 交给队员 ZYX 一个任务,编制一个称之为“文本生成器”的电脑软件:该软件的使用者是一些低幼人群,他们现在使用的是 GW 文本生成器 v6 版。 该软件可以随机生成一些文章——总是生成一篇长度固定且完全随机的文章。 也就是说,生成的文章中每个字符都是完全随机的。如果一篇文章中至少包含使用者们了解的一个单词,那么我们说这篇文章是可读的(我们称文章 $s$ 包含单词 $t$,当且仅当单词 $t$ 是文章 $s$ 的子串)。但是,即使按照这样的标准,使用者现在使用的 GW 文本生成器 v6 版所生成的文章也是几乎完全可读的。ZYX 需要指出 GW 文本生成器 v6 生成的所有文本中,可读文本的数量,以便能够成功获得 v7 更新版。你能帮助他吗? 答案对 $10^4 + 7$ 取模。 ## 输入格式 第一行有两个整数,分别表示使用者了解的单词总数 $n$ 和生成的文章长度 $m$。 接下来 $n$ 行,每行一个字符串 $s_i$,表示一个使用者了解的单词。 ## 输出格式 输出一行一个整数表示答案对 $10^4 + 7$ 取模的结果。 ## 输入输出样例 #1 ### 输入 #1 ``` 2 2 A B ``` ### 输出 #1 ``` 100 ``` ## 说明/提示 #### 数据规模与约定 对于全部的测试点,保证: - $1 \leq n \leq 60$,$1 \leq m \leq 100$。 - $1 \leq |s_i| \leq 100$,其中 $|s_i|$ 表示字符串 $s_i$ 的长度。 - $s_i$ 中只含大写英文字母。 为什么WA了(过了样例): ```cpp #include <bits/stdc++.h> using namespace std; const int MOD = 1e4 + 7; const int MAX_N = 100; const int LENGTH = 10050; int n, m, trie[LENGTH][26], tot; int fail[LENGTH], dp[MAX_N][LENGTH]; bool vis[LENGTH]; void insert(string s) { int cur = 0; for (char ch : s) { int i = ch - 'A'; if (!trie[cur][i]) trie[cur][i] = ++tot; cur = trie[cur][i]; } vis[cur] = true; } void build_fail() { queue<int> q; for (int i = 0; i < 26; i++) if (trie[0][i]) q.push(trie[0][i]); while (!q.empty()) { int cur = q.front(); q.pop(); vis[cur] |= vis[fail[cur]]; for (int i = 0; i < 26; i++) if (trie[cur][i]) { fail[trie[cur][i]] = trie[fail[cur]][i]; q.push(trie[cur][i]); } else trie[cur][i] = trie[fail[cur]][i]; } } int main() { cin >> n >> m; for (int i = 1; i <= m; i++) { string s; cin >> s; insert(s); } build_fail(); dp[0][0] = 1; for (int i = 0; i < n; i++) for (int j = 0; j <= tot; j++) for (int k = 0; k < 26; k++) if (!vis[trie[j][k]]) dp[i + 1][trie[j][k]] = (dp[i + 1][trie[j][k]] + dp[i][j]) % MOD; int ans = 1; for (int i = 1; i <= n; i++) ans = ans * 26 % MOD; for (int i = 0; i <= tot; i++) ans = (ans - dp[n][i]) % MOD; cout << (ans + MOD) % MOD << '\n'; return 0; } ```
最新发布
08-02
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值