[BJOI2017]魔法咒语(AC自动机+DP+矩阵快速幂)

文章目录

  • title
  • solution
  • code

title

solution

针对数据编程才是坠吊的!!!
观察数据,发现分隔数据的LLL跨度过大,没有衔接——推测很有可能是分数据做法


①:考虑L≤100L\le100L100的情况
可以暴力DPDPDP转移
dp[i][j]:dp[i][j]:dp[i][j]: 长度为iii指向自动机上jjj节点的方案数
枚举基本串kkk,写个find()find()find()函数暴力找jjj节点跳基本串kkk后指向的节点tototo
dp[i+len[k]][to]=(dp[i+len[k]][to]+dp[i][j])%moddp[i+len[k]][to]=(dp[i+len[k]][to]+dp[i][j])\%moddp[i+len[k]][to]=(dp[i+len[k]][to]+dp[i][j])%mod
O(L3)O(L^3)O(L3)完全可以,非常可以,很好很好


②:考虑基本串的长度只有1/21/21/2

思考一:
此时LLL急速飙升,一般上了1e7,1e8,1e91e7,1e8,1e91e7,1e8,1e9的数据点不是有规律周期可循,就是n,logn\sqrt n,lognn,logn
这里因为长度特别小,加上之前的做题经验,知道ACACAC自动机套矩阵快速幂的套路
于是自然而然就往矩阵加速方面靠了

思考二:
对于长度为lenlenlen的一个串sss,只有可能是由长度为s−1s-1s1的一个固定的串s′s's,或者长度为s−2s-2s2的一个固定的串s′′s''s转移而来
如果借用上面的dpdpdp思路,应该dp[len][to]=(dp[len−1][j]+dp[len−2][k])%moddp[len][to]=(dp[len-1][j]+dp[len-2][k])\%moddp[len][to]=(dp[len1][j]+dp[len2][k])%mod
这个狮子不仅代表的是转移方程,也在暗示每个状态只会与之前的最多两个状态挂钩
眼尖的犇犇更会发现有点像斐波拉契递推式——此时也会联想到矩阵加速求解


有长度为111的基本串,也有长度为222的基本串,矩阵要扩大两倍
也可以理解为把N×NN\times NN×N的一个矩阵看成一项fif_ifi,可以类比斐波拉契的转移推导过程

PS:下面的nnn仅代表一个字符,含义并不是题目里的nnn

初始矩阵OOO
(x0x1x2..xnx0′x1′..xn′)\begin{pmatrix} x_0\\ x_1\\ x_2\\ .\\ .\\ x_n\\ x_0'\\ x_1'\\ .\\ .\\ x_n' \end{pmatrix}x0x1x2..xnx0x1..xn
我们想让ta转移成
(x0′x1′x2′..xn′x0′′x1′′..xn′′)\begin{pmatrix} x_0'\\ x_1'\\ x_2'\\ .\\ .\\ x_n'\\ x_0''\\ x_1''\\ .\\ .\\ x_n'' \end{pmatrix}x0x1x2..xnx0x1..xn
在这里插入图片描述

code

#include <queue>
#include <cstdio>
#include <cstring>
using namespace std;
#define int long long
#define mod 1000000007
int n, m, L, tot;
char s[105];
int len[55];
int dp[105][105];
char S[55][105];queue < int > q;struct node {int trie[105][30], fail[105];bool End[105];void insert() {int now = 0, Len = strlen( s );for( int i = 0;i < Len;i ++ ) {int nxt = s[i] - 'a';if( ! trie[now][nxt] ) trie[now][nxt] = ++ tot;now = trie[now][nxt];}End[now] = 1;}void Fail() {memset( fail, 0, sizeof( fail ) );while( ! q.empty() ) q.pop();fail[0] = 0;for( int i = 0;i < 26;i ++ )if( trie[0][i] ) {fail[trie[0][i]] = 0;q.push( trie[0][i] );}while( ! q.empty() ) {int now = q.front(); q.pop();End[now] |= End[fail[now]];for( int i = 0;i < 26;i ++ ) {if( trie[now][i] ) {fail[trie[now][i]] = trie[fail[now]][i];q.push( trie[now][i] );}elsetrie[now][i] = trie[fail[now]][i];}}}int find( int now, int id ) {//求从now开始跳id串后指向的节点for( int i = 0;i < len[id];i ++ ) {now = trie[now][S[id][i] - 'a'];if( End[now] ) return -1;}return now;/*这种写法是错误的从End[now]处理的是从根节点到now这一条串的路上是否含有一个被禁止的子串从被禁止的子串开始往下跳不算for( int i = 0;i < len[id];i ++ ) {if( End[now] ) return -1;else now = trie[now][S[id][i] - 'a'];}if( End[now] ) return -1;else return now;*/}}AC;struct Matrix {int c[205][205];Matrix() {memset( c, 0, sizeof( c ) );}Matrix operator * ( const Matrix &p ) const {Matrix ans;for( int i = 0;i <= ( tot << 1 | 1 );i ++ )for( int k = 0;k <= ( tot << 1 | 1 );k ++ )if( c[i][k] )for( int j = 0;j <= ( tot << 1 | 1 );j ++ )ans.c[i][j] = ( ans.c[i][j] + c[i][k] * p.c[k][j] % mod ) % mod;return ans;}}O, V;Matrix qkpow( Matrix x, int y ) {Matrix ans;for( int i = 0;i <= ( tot << 1 | 1 );i ++ )ans.c[i][i] = 1;while( y ) {if( y & 1 ) ans = ans * x;x = x * x;y >>= 1;}return ans;
}
//dp[i][j]:长度为i指向自动机上j的方案数
void subtask1() {dp[0][0] = 1;for( int i = 0;i <= L;i ++ ) { for( int j = 0;j <= tot;j ++ ) {if( ! dp[i][j] ) continue;for( int k = 1;k <= n;k ++ ) {int to = AC.find( j, k ); if( ! ( ~ to ) || i + len[k] > L ) continue;else dp[i + len[k]][to] = ( dp[i + len[k]][to] + dp[i][j] ) % mod;}}}int ans = 0;for( int i = 0;i <= tot;i ++ ) ans = ( ans + dp[L][i] ) % mod;printf( "%lld\n", ans );
}void subtask2() {O.c[0][0] = 1;for( int i = 0;i <= tot;i ++ )V.c[i + tot + 1][i] = 1;for( int i = 0;i <= tot;i ++ ) {for( int j = 1;j <= n;j ++ ) {int to = AC.find( i, j );if( ! ( ~ to ) ) continue;if( len[j] & 1 ) {V.c[i + tot + 1][to + tot + 1] ++;if( ! i ) O.c[0][to + tot + 1] ++;}elseV.c[i][to + tot + 1] ++;}}O = O * qkpow( V, L );int ans = 0;for( int i = 0;i <= tot;i ++ )ans = ( ans + O.c[0][i] ) % mod;printf( "%lld\n", ans );
}signed main() {scanf( "%lld %lld %lld", &n, &m, &L );for( int i = 1;i <= n;i ++ )scanf( "%s", S[i] ), len[i] = strlen( S[i] );for( int i = 1;i <= m;i ++ ) {scanf( "%s", s );AC.insert();}AC.Fail();if( L <= 100 ) subtask1();else subtask2();return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317746.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

P6242-[模板]线段树3【吉司机线段树】

正题 题目链接:https://www.luogu.com.cn/problem/P6242 题目大意 给出一个长度为nnn的序列aaa&#xff0c;mmm次要求支持操作 区间加上一个值kkk区间所有aia_iai​变为min{ai,k}min\{a_i,k\}min{ai​,k}区间求和区间求最大值区间求历史最大值 1≤n,q≤51051\leq n,q\leq 5\…

洛谷P1912:诗人小G(二分栈、决策单调性)

二分栈&#xff0c;就是通过二分维护的栈 &#xff08;逃&#xff09; 解析 本题的决策单调性可以说是显然 但是本题是同维度&#xff08;其实只有一维&#xff09;自左向右转移&#xff0c;分治的写法是不能奏效的 所以我们使用决策点调性的另一种实现方法&#xff1a;二分栈…

[国家集训队]middle(二分+主席树[中位数思维题])

文章目录点击查看solutioncode点击查看 solution 简单口胡一下就跑 考虑二分答案ansansans 区间[x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d][x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d][x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d] 大于等于ansansans的设为111&#xff0c;小于ansans…

使用Roslyn脚本化C#代码,C#动态脚本实现方案

来源&#xff1a;https://www.cnblogs.com/7tiny/p/10279349.html【前言】Roslyn 是微软公司开源的 .NET 编译器。编译器支持 C# 和 Visual Basic 代码编译&#xff0c;并提供丰富的代码分析 API。Roslyn不仅仅可以直接编译输出&#xff0c;难能可贵的就是上述描述中的开放了编…

Moving On Gym - 102222F

Moving On Gym - 102222F 题意&#xff1a; 有 n 个城市&#xff0c;q 次询问. 给出每个城市的危险度 r 和 城市的邻接矩阵. 每次询问给出 u、v、w&#xff0c;求从 u 到 v 且不经过其他危险度超过 w 的城市的最短路. 题解&#xff1a; floyd 变形 我队友一开始想的是每次…

10.27模拟 总结

前言 220pts 100200100 按照gg的建议从《暴力模式》转化为《切题模式》 但是感觉本次有点难阿… 进入石头门困境… 再四道题中反复横跳 说实话心态是炸的 T1&#xff08;伪&#xff09;正解出来后才安了一些心 由于这次数据的强度确实不高 T1和T4的做法其实都很假 过掉也有运气…

[构造训练]CF1227G Not Same,CF1375H Set Merging,CF1364E X-OR

文章目录T1&#xff1a;CF1227G Not SamesolutioncodeT2&#xff1a;CF1364E X-ORsolutioncodeT3&#xff1a;CF1375H Set Mergingsolutioncode~~脑子是个好东西&#xff0c;希望人人都有构造真的不是个东西&#xff0c;看了一天视频&#xff0c;没有一道题会做~~ T1&#xff…

asp.net core 环境(Development、Staging 、Production)

一.在asp.net core中使用多个环境ASP.NET Core 配置是基于运行时环境, 使用环境变量。ASP.NET Core 在应用启动时读取环境变量ASPNETCORE_ENVIRONMENT&#xff0c;并将该值存储在 IHostingEnvironment.EnvironmentName 中。ASPNETCORE_ENVIRONMENT 可设置为任意值&#xff0c;但…

AT3913-XOR Tree【状压dp】

正题 题目链接:https://www.luogu.com.cn/problem/AT3913 题目大意 给出一棵有边权的树&#xff0c;你每次可以选择一条链让所有的边异或上同一个值&#xff0c;求最少的操作次数使得所有边的权值都为000。 2≤n≤105,0≤w<162\leq n\leq 10^5,0\leq w<162≤n≤105,0≤…

2018宁夏邀请赛

2018宁夏邀请赛 A Maximum Element In A Stack B Gym 102222B Rolling The Polygon C Gym 102222C Caesar Cipher D Gym 102222D Take Your Seat E Gym 102222E 2-3-4 Tree F Gym 102222F Moving On G Gym 102222G Factories H Gym 102222H Fight Against Monsters I Gym 10222…

洛谷P3349:小星星(容斥dp)

解析 先安利一波洛谷上我介绍如何用暴力日过去的博客 现在开始务正业 考虑把dp记录状态的一维s去掉 这样单次转移复杂度变成n3n^3n3 但是这样显然会算多啊&#xff01; 因为一个编号可能会用很多次 考虑容斥 设ansians_iansi​表示至少浪费了i个编号的答案 那么我们的答案显然…

ASP.NET Core 2.1 : 图解路由(2.1 or earler)

本文通过一张图来看一下路由的配置以及请求处理的机制。 一、概述路由主要有两个主要功能&#xff1a;将请求的URL与已定义的路由进行匹配&#xff0c;找到该URL对应的处理程序并传入该请求进行处理。根据已定义的路由生成URL这两个功能看起来这两个是相反的。A.路由的配置路由…

P2564 [SCOI2009]生日礼物

P2564 [SCOI2009]生日礼物 题意&#xff1a; n个彩珠&#xff0c;k个种类&#xff0c;分布在一个彩带上&#xff0c;现在要选取彩带的一部分&#xff0c;要求该部分包含所有种类的彩珠&#xff0c;且长度尽可能短&#xff0c;你能计算这个最短的长度吗&#xff1f; 1≤N≤100…

P6478-[NOI Online #2 提高组]游戏【dp,二项式反演】

正题 题目链接:https://www.luogu.com.cn/problem/P6478 题目大意 给出2m2m2m个点的一棵树&#xff0c;有mmm个白点mmm个黑点。 每个白点匹配一个黑点。 对于每个k∈[0,m]k\in[0,m]k∈[0,m]求恰好有kkk个匹配存在祖孙关系的方案。 答案对998244353998244353998244353取模。…

洛谷P6302:回家路线(斜率优化)

解析 《论什么是合理的实现》 本题dp的斜率式子还是不难 恶心在其他地方 由于不能时光倒流&#xff0c;新点必须在q时间后再插入 因此我们开一个堆来按找q升序排列&#xff0c;算完一个点就塞到堆里&#xff0c;每次把当前可以插入的点一起弹出来插入 一个重要的技巧是只存储火…

[数据结构专训][GXOI/GZOI2019]旧词,[hdu5118]GRE Words Once More!,[hdu6333]Problem B. Harvest of Apples

文章目录T1&#xff1a;[GXOI/GZOI2019]旧词solutioncodeT2&#xff1a;GRE Words Once More!solutioncodeT3&#xff1a;Problem B. Harvest of ApplessolutioncodeT1&#xff1a;[GXOI/GZOI2019]旧词 点击查看 solution 考虑k1k1k1的情况 由于dep[lca(x,y)]∣{z,zdep[lca(…

浅谈c#垃圾回收机制(GC)

写了一个window服务&#xff0c;循环更新sqlite记录&#xff0c;内存一点点稳步增长。三天后&#xff0c;内存溢出。于是&#xff0c;我从自己的代码入手&#xff0c;查找到底哪儿占用内存释放不掉&#xff0c;最终明确是调用servicestack.ormlite更新sqlite数据库造成的。至于…

动态规划专题复习

很久没做过dp了&#xff0c;都忘得差不多 普通dp 管道取珠 牛客网【每日一题】3月26日 合并回文子串 Music Problem 线性dp 树形dp 概率dp

Loj#6039-「雅礼集训 2017 Day5」珠宝【四边形不等式,dp】

正题 题目链接:https://loj.ac/p/6039 题目大意 有nnn个物品&#xff0c;第iii个费用为wiw_iwi​&#xff0c;价值为viv_ivi​&#xff0c;对于k∈[1,m]k\in[1,m]k∈[1,m]求费用为mmm时能获得的最大价值。 1≤n≤106,1≤m≤5104,1≤wi≤300,1≤vi≤1091\leq n\leq 10^6,1\leq…

洛谷P2497:基站建设(splay、斜率优化)

所谓splay斜率优化dp&#xff0c;就是利用splay和斜率对dp进行优化 &#xff08;逃&#xff09; 解析 在斜优的时候&#xff0c;有时我们会发现我们插入的点的横坐标并不单调 这个时候我们就无法利用单调队列维护凸包了 这时&#xff0c;我们就要请出今天的主角&#xff1a;s…