您的位置:首页 > 其它

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;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签:  dp