矩阵乘法+\(AC\)自动机

是道很不错的题了

首先是前六十分,就是一个\(AC\)自动机上的套路\(dp\),设\(dp[i][j]\)表示匹配出的长度为\(i\)在自动机上位置为\(j\)的方案数,转移的话就枚举下一个单词选择哪个放到自动机上一波匹配就好了

后面\(40\)分强行变成了另外一道题,\(L\)变成了\(1e8\),一看就是矩乘的复杂度了

但是单词的长度都非常小,于是转移\(dp[i][j]\)的时候只需要从\(dp[i-1][]\)和\(dp[i-2][]\)里转移,发现这非常像斐波那契的转移,于是提前在\(ac\)机上的每个位置都处理一下对应的转移之后矩乘就好了

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#define re register
#define LL long long
#define maxn 205
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
const LL mod=1e9+7;
char S[maxn];
int fail[maxn],flag[maxn],son[maxn][26];
char T[55][maxn],len[maxn];
int n,m,L,cnt;
inline void ins()
{
scanf("%s",S+1);
int len=strlen(S+1),now=0;
for(re int i=1;i<=len;i++)
{if(!son[now][S[i]-'a']) son[now][S[i]-'a']=++cnt;now=son[now][S[i]-'a'];}
flag[now]=1;
}
inline void Build()
{
std::queue<int> q;
for(re int i=0;i<26;i++) if(son[0][i]) q.push(son[0][i]);
while(!q.empty())
{
int k=q.front();q.pop();
flag[k]|=flag[fail[k]];
for(re int i=0;i<26;i++)
if(son[k][i]) fail[son[k][i]]=son[fail[k]][i],q.push(son[k][i]);
else son[k][i]=son[fail[k]][i];
}
}
namespace solve1
{
int dp[maxn][maxn];
inline int query(int x,int y)
{
int now=x;
for(re int i=1;i<=len[y];i++)
{
if(flag[now]) return -1;
now=son[now][T[y][i]-'a'];
}
if(flag[now]) return -1;
return now;
}
inline void work()
{
dp[0][0]=1;
for(re int i=0;i<L;i++)
for(re int j=0;j<=cnt;j++)
for(re int k=1;k<=n;k++)
{
if(i+len[k]>L) continue;
if(!dp[i][j]) continue;
int v=query(j,k);
if(v==-1) continue;
dp[i+len[k]][v]=(dp[i+len[k]][v]+dp[i][j])%mod;
}
int ans=0;
for(re int i=0;i<=cnt;i++) ans=(ans+dp[L][i])%mod;
printf("%d\n",ans);
}
}
namespace solve2
{
LL ans[maxn][maxn],a[maxn][maxn];
int M;
inline void did_a()
{
LL mid[maxn][maxn];
for(re int i=0;i<=M;i++)
for(re int j=0;j<=M;j++) mid[i][j]=a[i][j],a[i][j]=0;
for(re int k=0;k<=M;k++)
for(re int i=0;i<=M;i++)
for(re int j=0;j<=M;j++)
{a[i][j]+=((mid[i][k]*mid[k][j])%mod);if(a[i][j]>mod) a[i][j]%=mod;}
}
inline void did_ans()
{
LL mid[maxn][maxn];
for(re int i=0;i<=M;i++)
for(re int j=0;j<=M;j++) mid[i][j]=ans[i][j],ans[i][j]=0;
for(re int k=0;k<=M;k++)
for(re int i=0;i<=M;i++)
for(re int j=0;j<=M;j++)
{ans[i][j]+=((a[i][k]*mid[k][j])%mod);if(ans[i][j]>mod) ans[i][j]%=mod;}
}
inline void quick(int b){while(b) {if(b&1) did_ans();b>>=1;did_a();}}
inline void work()
{
M=cnt+cnt+1;
for(re int i=0;i<=cnt;i++)
{
if(flag[i]) continue;
for(re int j=1;j<=n;j++)
if(len[j]==1)
{
int v=son[i][T[j][1]-'a'];
if(!flag[v]) a[v+cnt+1][i+cnt+1]++;
}
else if(len[j]==2)
{
int v=son[i][T[j][1]-'a'];
int vv=son[v][T[j][2]-'a'];
if(flag[v]||flag[vv]) continue;
a[vv+cnt+1][i]++;
}
}
for(re int j=cnt+1;j<=M;j++) a[j-cnt-1][j]++;
for(re int i=0;i<=M;i++) ans[i][i]=1;
quick(L);
LL Ans=0;
for(re int i=cnt+1;i<=M;i++) Ans=(ans[i][cnt+1]+Ans)%mod;
printf("%lld\n",Ans);
}
}
int main()
{
scanf("%d%d%d",&n,&m,&L);
for(re int i=1;i<=n;i++) scanf("%s",T[i]+1),len[i]=strlen(T[i]+1);
for(re int i=1;i<=m;i++) ins();
Build();
if(L<=100) solve1::work();
else solve2::work();
return 0;
}

【[BJOI2017]魔法咒语】的更多相关文章

  1. bzoj4861 / P3715 [BJOI2017]魔法咒语

    P3715 [BJOI2017]魔法咒语 AC自动机+dp+矩阵乘法 常规思路是按基本串建立AC自动机 然鹅这题是按禁忌串建立AC自动机 对后缀是禁忌的点以及它的失配点做上标记$(a[i].ed)$, ...

  2. [BJOI2017]魔法咒语 --- AC自动机 + 矩阵优化

    bzoj 4860   LOJ2180   洛谷P3175 [BJOI2017]魔法咒语 题目描述: Chandra 是一个魔法天才. 从一岁时接受火之教会洗礼之后,Chandra 就显示出对火元素无 ...

  3. Luogu-3250 [BJOI2017]魔法咒语(AC自动机,矩阵快速幂)

    Luogu-3250 [BJOI2017]魔法咒语(AC自动机,矩阵快速幂) 题目链接 题解: 多串匹配问题,很容易想到是AC自动机 先构建忌讳词语的AC自动机,构建时顺便记录一下这个点以及它的所有后 ...

  4. P3715 [BJOI2017]魔法咒语

    P3715 [BJOI2017]魔法咒语 用基本词汇组成\(L\)长度的单词,其中不能包含禁忌词汇 用禁忌词汇建强大的\(tire\)图 解决: 分类讨论,\(L<=100\)用普通dp暴力在\ ...

  5. [BZOJ4861][BJOI2017]魔法咒语(AC自动机+矩阵优化DP)

    4861: [Beijing2017]魔法咒语 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 217  Solved: 105[Submit][Sta ...

  6. [BJOI2017]魔法咒语

    Description Chandra 是一个魔法天才. 从一岁时接受火之教会洗礼之后, Chandra 就显示出对火元素无与伦比的亲和力,轻而易举地学会种种晦涩难解的法术.这也多亏 Chandra ...

  7. 【题解】AC自动机题解合集

    最近貌似大家都在搞字符串?很长一段时间都没有写博客了……还是补一补坑吧. 感觉AC自动机真的非常优美了,通过在trie树上建立fail指针可以轻松解决多模匹配的问题.实际上在AC自动机上的匹配可以看做 ...

  8. AHOI2018训练日程(3.10~4.12)

    (总计:共90题) 3.10~3.16:17题 3.17~3.23:6题 3.24~3.30:17题 3.31~4.6:21题 4.7~4.12:29题 ZJOI&&FJOI(6题) ...

  9. Trie图(AC自动机)总结

    AC自动机构建完成后,某个节点沿着Fail链向上能从长到短走到自己的所有后缀.一般的,遍历主串进行匹配,就是在Trie图上定向移动的过程. 构造(一遍 BFS) void build_AC() { ; ...

随机推荐

  1. [javaSE] 基本类型(String相关)

    字符串是一个特殊的对象 字符串一旦初始化就不可以被改变 获取字符串的长度 调用String对象的length()方法,返回int长度 获取某个索引位置的字符 调用String对象的charAt()方法 ...

  2. kubernetes学习资源

    参考文章: 1.kubernetes学习资源 1. <Kubernetes与云原生应用>系列之Kubernetes的系统架构与设计理念 2.[docker专业介绍的网站dockerinfo ...

  3. confirm(message)

    如果用户点击确定按钮,则 confirm() 返回 true.如果点击取消按钮,则 confirm() 返回 false. 在用户点击确定按钮或取消按钮把对话框关闭之前,它将阻止用户对浏览器的所有输入 ...

  4. c语言中函数的形参test(int *&a)?

    今天在看一段c代码的时候看到一个函数的形参是(int *&a)居然是这个东西,这让我好生疑惑啊,不知道用这么多的地址符号用意何在呢?传址么? 那也不必这样,只需要用指针完全能够达到这样的效果啊 ...

  5. jQuery源码学习笔记一

    学习jQuery源码,我主要是通过妙味视频上学习的.这里将所有的源码分析,还有一些自己弄懂过程中的方法及示例整理出来,供大家参考. 我用的jquery v2.0.3版本. var rootjQuery ...

  6. 浅析javascript高阶函数

    什么是高阶函数:在数学和计算机科学中,高阶函数是至少满足下列一个条件的函数: 1. 接受一个或多个函数作为输入: 2. 输出一个函数.在数学中它们也叫做算子(运算符)或泛函.微积分中的导数就是常见的例 ...

  7. react解决roadhog buildDll 【转】

    本地删了 node module 目录,重新安装的时候,提示 找了找,可如下解决 ------- 转自: https://www.cnblogs.com/huhanhaha/p/7605722.htm ...

  8. ASP.Net与JSP如何共享Session值

     思路: ASP.NET中序列化Session以二进制数据保存到数据库,然后由JSP读取数据库中的二进制数据反序列化成Session对象,再强制转化成JAVA的Session对象登录的ASPX文件 ...

  9. redis list 查询、下标查询、删除、裁剪、压入弹出、队列实现

    查询  lrange list 0 1 // 注意0和1之间是空格:这个命令和pop命令不一样,不会删除里面的值lrange list 0 -1 // 所有的 下标查询 lpush person zs ...

  10. Windows 安装mkvirtualenv虚拟python环境

    pip install virtualenvwrapper-win mkvirtualenv --python=python.exe 新python虚拟环境名称 使用方法 所有的命令可使用:virtu ...