bzoj1030(ac自动机+dp)
2016-08-06 17:21
447 查看
设f[i][j]表示走i步,现在在j号节点的路径条数。
最后把所有f[m][i]累和就是不可读的串,一般问含有任意一个的方案数,就可以用补集的思想,总方案数减去不含有任何一个的方案数,剩下的就一定含有任意一个了!
注意一个节点的fail节点,及表示它后缀的节点,如果是单词节点,那么这个点可以在getfail时直接标记为单词节点,这样就处理出树上所有不可到达的节点,实际求不包含任意一个,就是在自动机上走m次,不经过任意一个单词节点
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
using namespace std;
const int mod=10007;
int ch[10005][30],tot,val[10005],fail[10005],last[10005],f[105][10005];
int n,m;
char s[105];
inline int idx(char ch) {return ch-'A';}
void insert(char *s)
{
int k=0,len=strlen(s);
for (int i=0;i<len;i++)
{
int id=idx(s[i]);
if (!ch[k][id]) ch[k][id]=++tot;
k=ch[k][id];
}
val[k]=1;
}
void getfail()
{
queue<int> q;
for (int i=0;i<26;i++)
if (ch[0][i]) q.push(ch[0][i]);
while (!q.empty())
{
int v,u=q.front();q.pop();
for (int i=0;i<26;i++)
if (v=ch[u][i])
{
int j=fail[u];
while (ch[j][i]==0&&j) j=fail[j];
j=ch[j][i];
fail[v]=j;
if (val[fail[v]]) val[v]=1;//如果他的后缀是不可行的,他也是不可行的,由此标记节点是否可行,感觉ac自动机这个位置很灵活
q.push(v);
}
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
{
scanf("%s",s);
insert(s);
}
getfail();
int sum=1,ans=0;
for (int i=1;i<=m;i++) sum=sum*26%mod;
f[0][0]=1;//最开始,走了0步,走到了0位置的方案数
for (int i=0;i<m;i++)//
{
for (int j=0;j<=tot;j++)
if (f[i][j]&&!val[j])//这个位置不能是单词末尾,且必须到达过
for (int k=0;k<26;k++)//枚举下一个字母是
{
int l=j;
while (!ch[l][k]&&l) l=fail[l];//沿失配边走,找到第一个可以放k这个字母的位置
l=ch[l][k];//如果到了0这个位置,没有k儿子,那么还会是0位置,更新到它的方案数,表示加上k这个字母,匹配数为0
if (!val[l]) f[i+1][l]=(f[i+1][l]+f[i][j])%mod;//如果该点不是单词节点,即加上新的字母,原文本不会出现单词,就进行转移,否则该情况不合法
}
}
for (int i=0;i<=tot;i++) if (!val[i]) ans=(ans+f[m][i])%mod;//刚开始忘了还有0,就是到最后匹配了0的位置
printf("%d",((sum-ans)%mod+mod)%mod);
return 0;
}
最后把所有f[m][i]累和就是不可读的串,一般问含有任意一个的方案数,就可以用补集的思想,总方案数减去不含有任何一个的方案数,剩下的就一定含有任意一个了!
注意一个节点的fail节点,及表示它后缀的节点,如果是单词节点,那么这个点可以在getfail时直接标记为单词节点,这样就处理出树上所有不可到达的节点,实际求不包含任意一个,就是在自动机上走m次,不经过任意一个单词节点
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
using namespace std;
const int mod=10007;
int ch[10005][30],tot,val[10005],fail[10005],last[10005],f[105][10005];
int n,m;
char s[105];
inline int idx(char ch) {return ch-'A';}
void insert(char *s)
{
int k=0,len=strlen(s);
for (int i=0;i<len;i++)
{
int id=idx(s[i]);
if (!ch[k][id]) ch[k][id]=++tot;
k=ch[k][id];
}
val[k]=1;
}
void getfail()
{
queue<int> q;
for (int i=0;i<26;i++)
if (ch[0][i]) q.push(ch[0][i]);
while (!q.empty())
{
int v,u=q.front();q.pop();
for (int i=0;i<26;i++)
if (v=ch[u][i])
{
int j=fail[u];
while (ch[j][i]==0&&j) j=fail[j];
j=ch[j][i];
fail[v]=j;
if (val[fail[v]]) val[v]=1;//如果他的后缀是不可行的,他也是不可行的,由此标记节点是否可行,感觉ac自动机这个位置很灵活
q.push(v);
}
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
{
scanf("%s",s);
insert(s);
}
getfail();
int sum=1,ans=0;
for (int i=1;i<=m;i++) sum=sum*26%mod;
f[0][0]=1;//最开始,走了0步,走到了0位置的方案数
for (int i=0;i<m;i++)//
{
for (int j=0;j<=tot;j++)
if (f[i][j]&&!val[j])//这个位置不能是单词末尾,且必须到达过
for (int k=0;k<26;k++)//枚举下一个字母是
{
int l=j;
while (!ch[l][k]&&l) l=fail[l];//沿失配边走,找到第一个可以放k这个字母的位置
l=ch[l][k];//如果到了0这个位置,没有k儿子,那么还会是0位置,更新到它的方案数,表示加上k这个字母,匹配数为0
if (!val[l]) f[i+1][l]=(f[i+1][l]+f[i][j])%mod;//如果该点不是单词节点,即加上新的字母,原文本不会出现单词,就进行转移,否则该情况不合法
}
}
for (int i=0;i<=tot;i++) if (!val[i]) ans=(ans+f[m][i])%mod;//刚开始忘了还有0,就是到最后匹配了0的位置
printf("%d",((sum-ans)%mod+mod)%mod);
return 0;
}
相关文章推荐
- bzoj1030: [JSOI2007]文本生成器(AC自动机+Dp)
- [BZOJ1030][JSOI2007]文本生成器(AC自动机+dp)
- bzoj 1030 [JSOI2007]文本生成器(AC自动机+DP)
- BZOJ1030 [JSOI2007]文本生成器(AC自动机+dp)
- bzoj 1030: [JSOI2007]文本生成器 (ac自动机上的dp)
- [BZOJ1030][JSOI2007]文本生成器(AC自动机+dp)
- [BZOJ 1030][JSOI2007]文本生成器(AC自动机+DP)
- BZOJ1030 [JSOI2007]文本生成器 补全AC自动机+简单DP
- bzoj1030 文本生成器(AC自动机+dp)
- 【bzoj 1030】文本生成器(AC自动机+DP)
- bzoj 1030: [JSOI2007]文本生成器(AC自动机+DP)
- 【bzoj 1030】 [JSOI2007]文本生成器(AC自动机+dp)
- BZOJ_1030_[JSOI2007]_文本生成器_(AC自动机+DP)
- bzoj 3796: Mushroom追妹纸 AC自动机+后缀自动机+dp
- bzoj 1212: [HNOI2004]L语言(AC自动机+DP)
- [BZOJ3530][SDOI2014]数数(AC自动机+数位dp)
- bzoj 1030: [JSOI2007]文本生成器 (ac自己主动机上的dp)
- bzoj1559 [JSOI2009]密码(AC自动机+状压DP)
- [bzoj1030][JSOI2007]文本生成器——AC自动机
- [BZOJ3530][Sdoi2014]数数(AC自动机+数位DP)