矩阵(快速幂/高斯消元/子矩阵)
文章平均质量分 61
矩阵相关
Code92007
No Saturday , no Sunday , no holiday .
展开
-
Atcoder TUPC 2023(東北大学プログラミングコンテスト 2023)E. And DNA(矩阵快速幂+拆位讨论)
2. 如果m=1,由于1+(1&1)=2,0+(0&x)=0,所以不能有三个1相邻,不能有两个0相邻。对于i∈[1,n],均满足a[i]+(a[i-1]&a[i+1])=m,求这样可能的数组的方案数。①如果是1,只能是1+(1&1)=2,此时还给倒数第二位贡献了1,需要再递归求m/2-1的方案。长为n(3<=n<=1e9)的数组,第i个数ai在[0,m](m<=1e9)之间。特别地,认为a[0]=a[n],a[n+1]=a[1],即这个数组是个环形的数组。有f(m)=f(1)+f((m-1)/2)原创 2024-03-15 20:49:43 · 419 阅读 · 0 评论 -
AtCoder Beginner Contest 258 Ex. Odd Steps(带限制的矩阵快速幂)
给定前缀和值s和长为n(n<=1e5)的序列A,第i个数ai(1<=a1<=a2<=...<=an<=s<=1e18)2. 给出矩阵快速幂转移式,并且有初始向量f[0]=1,s[-1]=0,s[-2]=0(边界条件)3. 有前缀和不能为ai的限制,所以在n次矩阵快速幂转移时,手动给当前f[i]的值置为0。求满足以下条件的序列X的数量,答案对998244353取模。3. 记X序列的前缀和序列为Y,Y中不包含给定序列a中的值。2. X序列中的数的和为s(s<=1e18)1. X序列中的每个元素都是正奇数。原创 2024-02-04 02:44:59 · 489 阅读 · 0 评论 -
Educational Codeforces Round 157 (Rated for Div. 2) F. Fancy Arrays(容斥+组合数学)
第一个数字选0,后面每个数都有2k+1种选择方式,最后把最小值往上平移到[0,x+k-1]之间。t(t<=50)组样例,每次给定n(1<=n<=1e9),x(1<=x<=40),由于长度为1时对应的[0,x-1]的向量均为1,所以将每一行的和从答案中减掉即可。注意到x<=40,所以可以dp[i][j]表示长为i的数组最后一个是j的方案数。转移时,只要abs(j1-j2)<=k,就可以从j1转移到j2,即长为n的数列,使用的值均在[0,x-1]的方案数,的字眼,首先想到容斥,用总的减不满足的,原创 2023-11-06 01:14:33 · 444 阅读 · 0 评论 -
Harbour.Space Scholarship Contest 2023-2024 (Div. 1 + Div. 2) H. Asterism Stream(期望dp+矩阵快速幂/生成函数)
你可能会说,f(2)的时候,0.25f(0)+0.25f(0)和0.5f(0)不是一样的么...题解1. 矩阵快速幂的比较难想,没遇到过矩阵快速幂维护x/2的,感觉也是一个典。但是,n稍微大一点的时候,n-1的第二行和第三行就不一样了,这种做法具有普适性。为了说明与ctz(二进制尾0的个数/最低位1的位置)的关系,举俩小点的例子。实际大概2的60次方到1e18,维护一个60*60的矩阵,需要用几行,即n和n-1相比,除以2的k次方什么时候不变,矩阵快速幂,以倍增的思想,求出处理出前n-1个转移矩阵积,原创 2023-09-04 06:12:02 · 76 阅读 · 0 评论 -
2019 ICPC Asia Nanchang Regional J. Summon(polya定理+矩阵快速幂优化dp)
题目n(4<=n<=1e5)个珠子的项链,项链由四种颜色0123组成,旋转相同时视为同种方案,m(0<=m<=256)个限制,第i次给出四个整数a b c d(0<=a,b,c,d<=3)表示顺时针看这串项链时,abcd段不能出现,求合法的方案数%998244353思路来源https://www.cnblogs.com/zxytxdy/p/12582065.html题解poj2888原题既视感,高度相似的一个题,只是变二维为四维了.原创 2020-06-11 16:29:50 · 497 阅读 · 0 评论 -
CCPC-Wannafly & Comet OJ 夏季欢乐赛(2019)E.飞行棋(期望dp+矩阵快速幂)
题目飞行棋的规则如下:1、每名玩家有一个棋子,每个回合可以掷一次骰子。2、如果使用的骰子为k面,则这k面上的点数分别为 1,2,3,…,k,且掷得每种点数的概率均为。3、如果当前回合掷得的点数为Q,则玩家控制的棋子前进Q步。4、若当前棋子的位置到终点的距离d < Q,则棋子先行动d步到终点,再倒退Q-d步。(即到终点的距离变为Q−d)5、某一回...原创 2019-08-06 12:34:04 · 358 阅读 · 0 评论 -
hdu6395 Sequence(矩阵快速幂+数论分块)
题目思路来源https://blog.csdn.net/Originum/article/details/82111066题解对于每个P/n向下取整相同值的地方用矩阵快速幂,[l,r]的值都为v==P/i,那么求r=P/(P/i),r+1就得到下一个l,是为数论分块g1,g2为last Fn-1,Fn-2的值,f1,f2为now Fn-1,Fn-2的值,类似滚动dp的存...原创 2019-07-02 17:51:33 · 147 阅读 · 0 评论 -
Codeforces Round #566 (Div. 2) E.Product Oriented Recurrence(矩阵快速幂)
题目思路来源涛神+贤神题解先构造出g(x)的系数矩阵,计算出cnum,num[1],num[2]和num[3]在要算g(x)的时候,偷梁换柱成f(x)计算,加上费马小定理降幂即可题解有左右同乘,令,则代码#include<iostream>#include<cstdio>#include<cstring>#in...原创 2019-07-01 19:51:20 · 200 阅读 · 0 评论 -
hdu6470 Count(矩阵快速幂)
题目第一头奶牛为1号,第二头奶牛为2号,第三头奶牛之后,假如当前奶牛是第n头,那么他的编号就是2倍的第n-2头奶牛的编号加上第n-1头奶牛的编号再加上自己当前的n的三次方为自己的编号.现在Farmer John想知道,第n头奶牛的编号是多少,估计答案会很大,你只要输出答案对于123456789取模.第一行输入一个T,表示有T组样例接下来T行,每行有一个正整数n,表示有n头奶牛 (n&g...原创 2019-03-18 15:56:25 · 330 阅读 · 0 评论 -
hdu5015 233 Matrix(矩阵快速幂)
题目给你一个矩阵,第一行是233,233,2333,然后会给你a[1][0],a[2][0],…,a[n][0]且a[i][j]=a[i-1][j]+a[i][j-1]%mod求a[n][m]思路来源https://www.cnblogs.com/whatbeg/p/3971994.html题解把第一行都左移一位就可以对系数矩阵快速幂了然后最后a[n]这一行...原创 2019-03-09 22:27:29 · 151 阅读 · 0 评论 -
poj3233 Matrix Power Series(矩阵快速幂 求前1-k次幂矩阵和)
题目输入n*n的矩阵A,求A的前1-k次幂矩阵和的矩阵S输出S各元素%m思路来源挑战程序竞赛 P205题解构造一个左上角是A(输入矩阵),右上角是0(零矩阵)左下角和右下角都是I(单位矩阵)的矩阵然后根据,就可以求了把构造的矩阵快速幂一下,则=构造矩阵的k次幂的左下角矩阵*I其实之前见过这个题……这不是要做hdu2243么整理一下……代码...原创 2019-02-15 11:51:42 · 283 阅读 · 1 评论 -
hdu6030 Happy Necklace(矩阵快速幂+递推式)
思路来源https://www.cnblogs.com/bolderic/p/7216739.html题意给你一个长为n的珠子串,要求你给这条串涂色,要求任意素数长度的串中,红色数不少于蓝色数,只有红蓝两色。问长为n的串的涂色方案数。题解这个递推式还是挺好想的,就是好久不用矩阵快速幂一下愣住了。dp写久了小范围的用不了矩阵快速幂的还是把自己写傻了QAQ注意到...原创 2018-11-05 22:02:28 · 179 阅读 · 0 评论