bzoj1030(ac自动机+dp)

来源:互联网 发布:大学网络宣传部 编辑:程序博客网 时间:2024/04/27 20:10

设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这个字母,匹配数为0if (!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;}


0 1
原创粉丝点击