[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,一经查实,立即删除!

相关文章

[国家集训队]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;难能可贵的就是上述描述中的开放了编…

[构造训练]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;但…

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

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

[数据结构专训][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数据库造成的。至于…

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

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

MediatR 知多少

引言首先不用查字典了&#xff0c;词典查无此词。猜测是作者笔误将Mediator写成MediatR了。废话少说&#xff0c;转入正题。先来简单了解下这个开源项目MediatR&#xff08;作者Jimmy Bogard&#xff0c;也是开源项目AutoMapper的创建者&#xff0c;在此表示膜拜&#xff09;&a…

网络分析(带权并查集)

网络分析 题意&#xff1a; 有n个节点&#xff0c;一开始彼此独立&#xff0c;有两个操作&#xff0c;第一个操作时是连接两个节点&#xff0c;第二个操作是对一个节点x&#xff0c;&#xff08;在进行第二个操作时&#xff0c;与该点相连的点也会x&#xff09; 问每个节点的…

良心发现,时隔一年再回首莫比乌斯反演(公式性质证明+题目练习)

文章目录莫比乌斯反演引入公式性质模板公式证明莫比乌斯函数前缀和题目练习完全平方数[HAOI2011]ProblembYY的GCD[SDOI2014]数表[国家集训队]Crash的数字表格/JZPTAB[SDOI2015]约数个数和寒假疫情期间跟着lmm学了一遍&#xff0c;完全是懵逼到底状态&#xff0c;以至于后面考到…

平面切分

平面切分 问题描述 题解&#xff1a; 我对这种题极其非常不擅长。。。 另外吐槽为什么acwing的数据卡的这么死&#xff0c;蓝桥杯官网数据那么水 其实题目很简单&#xff0c;如果只有一个直线&#xff0c;那么就是两部分&#xff0c;如果是两个直线&#xff0c;这两个直线不相…

新起点!新征程!微软技术俱乐部(苏州)成立大会暨微软技术交流会

2019年1月19日&#xff0c;苏州微软将举办苏州史上最盛大的开发者聚会&#xff0c;微软技术俱乐部成立大会暨微软技术交流会。超越苹果&#xff0c;登顶世界市值第一&#xff01;云与AI两大技术支柱支撑起的微软帝国&#xff0c;正向万亿美元俱乐部挺进&#xff01;微软的改变我…

跨平台、跨语言应用开发工具,Elements 介绍

目录1&#xff0c;Elements 介绍2&#xff0c;Elements 版本3&#xff0c;Elements 能干嘛4&#xff0c;Elements IDES5&#xff0c;Elements 工具1&#xff0c;Elements 介绍RemObjects Elements&#xff0c;是多平台移动项目开发工具&#xff0c;是一款可以帮助开发人员在 不…

后缀自动机(SAM)构造实现过程演示+习题集锦

文章目录后缀自动机算法实现过程模板习题洛谷后缀自动机模板题品酒大会[HEOI2015]最短不公共子串字符串蒟蒻写这篇blogblogblog主要是存一下&#xff0c;后缀自动机的详细搭建过程&#xff0c;方便以后复习 具体的某些证明&#xff0c;为什么这么做&#xff0c;正确性劈里啪啦一…

使用logdashboard进行可视化的日志追踪

本文源码在Github可以找到下载LogDashboard如果你还不了解LogDashboard请看这里。 LogDashboard 1.1版本支持请求追踪,虽然目前版本还没有发布。不过这个功能可以先睹为快效果图下载项目首先我们可以在 https://github.com/liangshiw/LogDashboard/tree/master/samples/Request…

ML.NET 0.9特性简介

ML.NET 0.9已于上周发布&#xff0c;距离上次0.8版本的发布只有一个多月&#xff0c;此次增加的新特性主要包括特征贡献计算&#xff0c;模型可解释性增强&#xff0c;ONNX转换对GPU的支持&#xff0c;Visual Studio ML.NET项目模板预览&#xff0c;以及API改进。特征贡献计算特…

[学习笔记] 初次见面,请多关照 (公式推导+题集)——杜教筛

筛积性函数的前缀和常见积性函数公式推导狄利克雷卷积杜教筛实现常见卷积习题集Sum神犇和蒟蒻简单的数学题常见积性函数 μ\muμφφφd(n)d(n)d(n)&#xff1a;nnn的约数个数σ(n)σ(n)σ(n)&#xff1a;nnn的约数和ϵ(n)ϵ(n)ϵ(n)&#xff1a;单位元函数&#xff0c;e(n)[n1…

AOP 还在配置吗改用打标签模式吧!

为什么我喜欢打标签来配置AOP1. 配置多很混乱&#xff0c;代码里面很难分辨出来哪些是AOP容器(比如属性注入)2. 对于代码生成器生成的代码里面还需要手动加到配置里面3. 连java spring现在都是清一色的注解来代替xml&#xff0c;这个就是趋势所在我基于Autofac开发了一个基于标…

DZY Loves Math IV(杜教筛)

文章目录titlesolutioncodetitle solution 这道题是多么的妙啊&#xff0c;完全不是我能推出来的式子呢&#xff01; 观察数据范围&#xff0c;有点奇怪欸&#xff0c;在暗示我&#xff1f;&#xff1f; 考虑暴力枚举nnn S(n,m)∑i1mφ(ni)S(n,m)\sum_{i1}^mφ(n\times i)S…