bzoj1009矩阵快速面+kmp

【bzoj1009矩阵快速面+kmp】的更多相关文章

bzoj1009矩阵快速面+kmp

其实kmp真的很次要,求长度为20的kmp感觉真的有点杀鸡用牛刀这题思路相当明确:一看题就是数位dp,一看n的大小就是矩阵矩阵的构造用m*m比较方便,本来想写1*m的矩阵乘m*m的,但是感觉想起来太麻烦就偷懒,没想到1A了 log的速度的确可以,87ms贼快,好久没见这么短的运行时间了 #include <cstdio> int n,m,mod,k; char ch; ],ne[]; ][],t[][],zy[][]; void mul(bool b) { ;i<m;i++) ;j&…

BZOJ1009 矩阵快速幂+DP+KMP

Problem 1009. -- [HNOI2008]GT考试 1009: [HNOI2008]GT考试 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 3773 Solved: 2314[Submit][Status][Discuss] Description 阿申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字. 他的不吉利数学A1A2...Am(0<=Ai<=9)…

BZOJ 1009: [HNOI2008]GT考试( dp + 矩阵快速幂 + kmp )

写了一个早上...就因为把长度为m的也算进去了... dp(i, j)表示准考证号前i个字符匹配了不吉利数字前j个的方案数. kmp预处理, 然后对于j进行枚举, 对数字0~9也枚举算出f(i, j)表示dp(x-1, j)对dp(x, i)的贡献.然后用矩阵快速幂就可以了. 时间复杂度O(M3logN + M) ------------------------------------------------------------------- #include<bits/stdc++.h>…

[HNOI2008] GT考试(DP+矩阵快速幂+KMP)

题目链接:https://www.luogu.org/problemnew/show/P3193#sub 题目描述阿申准备报名参加 GT 考试,准考证号为 N 位数 X1,X2…Xn(0 <= Xi <= 9) ,他不希望准考证号上出现不吉利的数字. 他的不吉利数学 A1,A2…Am(0≤Ai≤9) 有 M 位,不出现是指 X1,X2…Xn 中没有恰好一段等于 A1,A2…Am ,A1 输入输出格式输入格式: 第一行输入N,M,K.接下来一行输入M位的数. 输出格式:…

[BZOJ1009][HNOI2008]GT考试 DP+矩阵快速幂+KMP

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1009 我们令$dp(i,j)$表示已经填了$i$位,而且后缀与不幸运数字匹配了$j$位,那么转移方程就是$dp(i,j)=dp(i,k)*a(j,k)$,其中$a(j,k)$表示从$j$位可以转移到$k$位的方案数,这个可以利用kmp的ne数组来计数求得.因为DP方程是线性的,所以可以用矩阵优化. #include<cstdio> #include<cstring> #in…

BZOJ1009: [HNOI2008]GT考试矩阵快速幂+kmp+dp

这个题你发现打暴力的话可以记忆化搜素加剪枝,那么意味着可以递推,我们搜的话就是1010^9我们就往下匹配遇到匹配成功就return,那么我们可以想一下什么决定了状态,我们考虑kmp的过程,对于我们目前匹配到的距离,下一次在匹配时不会用他之后的字符,那么只要我们知道匹配到的距离和已匹配长度就行了,那么我们考虑状态的转移,我们由于要像kmp那样匹配于是我们只要知道在匹配到k位时往下走一个数时匹配到哪,算出a[k][j](在k时到j的方案数),那么新的f[i][j]=∑f[i-1][k]*a[k][j…