矩阵运算
Neo__Z
这个作者很懒,什么都没留下…
展开
-
NOIP模拟 Matrix 矩阵运算
题目描述给出两个 N×N 的矩阵 A、B,矩阵每行每列标号 0~N-1 。 定义这两个矩阵的乘积 AB 为(AB)i,j=∑k=0n−1Ai,jBi,j(AB)_{i,j}=\sum\limits_{k=0}^{n-1}A_{i,j}B_{i,j}现在要在这两个矩阵上依次进行 Q 次修改操作,两种操作描述如下: A i j K ,将 Ai,j 的值修改为 K 。 B i j K ,将 Bi,j原创 2017-07-22 21:54:13 · 500 阅读 · 0 评论 -
bzoj5298: [Cqoi2018]交错序列【二项式定理+动态规划+矩阵快速幂】
传送门解题思路:xayb=(n−y)ayb=∑i=0a(ai)ni(−1)a−iya+b−ixayb=(n−y)ayb=∑i=0a(ai)ni(−1)a−iya+b−ix^ay^b=(n-y)^ay^b=\sum\limits_{i=0}^a\binom{a}{i}n^i(-1)^{a-i}y^{a+b-i}设f[k][i][0/1]f[k][i][0/1]f[k][i][0/1]表示...原创 2018-05-03 14:12:40 · 305 阅读 · 0 评论 -
BJ模拟 排列【状压dp+矩阵快速幂】
题目大意:求长度为nnn,且满足对于任意i,|pi−i|≤ki,|pi−i|≤ki,|p_i-i|\le k的排列数。n≤1e9,k≤5n≤1e9,k≤5n\le 1e9,k\le 5解题思路:由于kkk很小,很容易想到记录一个2k+12k+12k+1位的状态表示当前位还有那些位置可以用,再用矩阵快速幂优化,但还是会TLE。注意到我们每一步转移恰好增加一个占用,至多减少一个占用,所...原创 2018-05-02 09:03:12 · 209 阅读 · 0 评论 -
bzoj4944: [Noi2017]泳池【概率dp+特征多项式】
题目大意:有一个1001×n1001×n1001×n的的网格,每个格子有qqq的概率是安全的,1−q1−q1−q的概率是危险的。 定义一个矩形是合法的当且仅当: 1.这个矩形中每个格子都是安全的 2.必须紧贴网格的下边界 问你最大的合法子矩形大小恰好为kkk的概率是多少。解题思路:首先求恰好为kkk的概率一般转化为求≤k≤k\le k的概率减去≤k−1≤k−1\le k-1...原创 2018-04-09 20:05:48 · 659 阅读 · 0 评论 -
BJ模拟 数列【特征多项式and生成函数】
题目大意:有数列:fm,n=⎧⎩⎨⎪⎪an,n=1...m∑k=1m(a−1)fm,n−k−1(1)(1)fm,n={an,n=1...m∑k=1m(a−1)fm,n−k−1\begin{align}f_{m,n}=\begin{cases}a^n,n=1...m\\\sum\limits_{k=1}^m(a-1)f_{m,n-k-1}\end{cases}\end{align}...原创 2018-03-15 10:53:03 · 449 阅读 · 1 评论 -
bzoj4161: Shlw loves matrixI【特征多项式优化常系数齐次线性递推模板】
Description给定数列 {hn}前k项,其后每一项满足 hn = a1*h(n-1) + a2*h(n-2) + … + ak*h(n-k) 其中 a1,a2…ak 为给定数列。请计算 h(n),并将结果对 1000000007 取模输出。Input第 1 行包含两个整数 n,k 第 2 行包含 k 个整数 a1,a2…ak 第 3 行包含 k 个整数h[0],h[...原创 2018-03-14 21:00:35 · 990 阅读 · 0 评论 -
hdu2243考研路茫茫——单词情结 【AC自动机+动态规划+矩阵快速幂】
这道题就是bzoj1030的翻版,那道题题解见这里。还是用总方案数减去一个单词都不包含的方案数,只不过节点很少,不到30个,就可以用矩阵快速幂优化。不过注意要求的是长度小于L的,所以矩阵多开一维记录前缀和,总方案数也变成了260+261+262+……+26L" role="presentation" style="position: relative;">260+2原创 2018-01-25 22:51:57 · 284 阅读 · 0 评论 -
bzoj2510 弱题 【期望dp+循环矩阵快速幂】
解题思路:很容易想到每一轮的dp方程: a[i]=(1−1m)a[i]+1ma[i−1=0?n:i−1]a[i]=(1-\frac{1}{m})a[i]+\frac{1}{m}a[i-1=0?n:i-1]考虑矩阵快速幂,可以O(n3log2k)O(n^3log_2k)求解,但仍会TLE。注意到转移矩阵是一个循环矩阵,即每一行都是上一行平移一格得到(列也是如此),而循环矩阵的性质有:原创 2018-01-05 18:00:48 · 328 阅读 · 0 评论 -
hdu5564 Clarke and digits【数位dp+矩阵快速幂】
题目大意:求长度在[l,r][l,r]之间的能被7整除且相邻数位之和不为kk的正整数有多少个。 1≤l,r≤1091\le l,r \le10^9解题思路:设f[x][i][j]f[x][i][j]表示第x位为j,模7余数为i的方案数,那么易得转移方程:f[x+1][(i∗10+k)f[x+1][(i*10+k)%7][k]+=f[x][i][j]7][k]+=f[x][i][j]。原创 2018-01-04 21:28:48 · 274 阅读 · 0 评论 -
51nod 1195 斐波那契数列的循环节【斐波那契数列&&二次剩余&&欧拉判定准则】
解题思路:先说明一下结论在下都不会证明,囧……。对于一个正整数n,我们求Fib数模n的循环节的长度的方法如下: (1)将n分解质因数,即n=pk11pk22……pkmm(1)将n分解质因数,即n=p_1^{k_1}p_2^{k_2}……p_m^{k_m} (2)分别计算Fib数模pkii的循环节的长度,假设是x1,x2,……xm(2)分别计算Fib数模p_i^{k_i}的循环节的长度,假设是x_原创 2017-12-08 19:30:21 · 566 阅读 · 0 评论 -
loj#2325. 「清华集训 2017」小 Y 和恐怖的奴隶主
传送门解题思路:设 p[i][j][k][l]p[i][j][k][l]p[i][j][k][l] 表示 iii 回合奴隶主 1、2、31、2、31、2、3 血剩余情况为 j、k、lj、k、lj 、k 、l 的概率,那么对答案的贡献就是 p[i][j][k][l]j+k+l+1p[i][j][k][l]j+k+l+1\frac{p[i][j][k][l]}{j+k+l+1} 。先预...原创 2018-05-24 20:54:22 · 588 阅读 · 0 评论