矩阵快速幂刷题小记

最新推荐文章于 2024-06-16 11:27:34 发布
最新推荐文章于 2024-06-16 11:27:34 发布
阅读量113 收藏
点赞数 3
分类专栏： 数论 文章标签： 矩阵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Sun_QH_/article/details/113868717
版权

矩阵快速幂刷题小记

矩阵快速幂模板

传送门： P3390 【模板】矩阵快速幂
Code:

#include<bits/stdc++.h>
using namespace std;
const long long p=1000000007;
inline long long Read(){
	long long dx=0,fh=1;
	char c=getchar();
	while(c>'9'||c<'0'){
		if(c=='-')	fh=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		dx=dx*10+c-'0';
		c=getchar();
	}
	return dx*fh;
}
long long k,n;
struct Matricx{
	long long a[101][101];
}A;
Matricx mul(Matricx x,Matricx y){
	Matricx z;
	for(long long i=1;i<=n;++i)
		for(long long j=1;j<=n;++j){
			z.a[i][j]=0;
			for(long long u=1;u<=n;++u)
				z.a[i][j]=(z.a[i][j]+(x.a[i][u]*y.a[u][j])%p)%p;
		}
	return z;
}
Matricx quick_pow(Matricx x,long long b){
	Matricx r,base=x;
	for(int i=1;i<=n;++i)	for(int j=1;j<=n;++j)	r.a[i][j]=0;
	for(int i=1;i<=n;++i)	r.a[i][i]=1;
	while(b){
		if(b&1)	r=mul(r,base);
		base=mul(base,base);
		b>>=1;
	}
	return r;
}
int main(){
	n=Read();k=Read();
	for(long long i=1;i<=n;++i)
		for(long long j=1;j<=n;++j)
			A.a[i][j]=Read();
	Matricx B=quick_pow(A,k);
	for(long long i=1;i<=n;++i){
		for(long long j=1;j<=n;++j)
			printf("%lld ",B.a[i][j]);
		putchar('\n');
	}	
	return 0; 
}

反思与总结：
刚定义的矩阵一定要先全部赋值为0

斐波那契数列

传送门:P1962 斐波那契数列
Code:

#include<bits/stdc++.h>
using namespace std;
const long long p=1000000007;
inline long long Read(){
	long long dx=0,fh=1;
	char c=getchar();
	while(c>'9'||c<'0'){
		if(c=='-')	fh=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		dx=dx*10+c-'0';
		c=getchar();
	}
	return dx*fh;
}
long long k;
struct Matricx{
	long long a[3][3];
}A;
Matricx mul(Matricx x,Matricx y){
	Matricx z;
	for(long long i=1;i<=2;++i)
		for(long long j=1;j<=2;++j){
			z.a[i][j]=0;
			for(long long u=1;u<=2;++u)
				z.a[i][j]=(z.a[i][j]+(x.a[i][u]*y.a[u][j])%p)%p;
		}
	return z;
}
Matricx quick_pow(Matricx x,long long b){
	Matricx r,base=x;
	for(int i=1;i<=2;++i)	for(int j=1;j<=2;++j)	r.a[i][j]=0;
	for(int i=1;i<=2;++i)	r.a[i][i]=1;
	while(b){
		if(b&1)	r=mul(r,base);
		base=mul(base,base);
		b>>=1;
	}
	return r;
}
int main(){
	k=Read();
	if(k==1||k==2)	printf("1");
	else{
		A.a[1][1]=1,A.a[1][2]=1,A.a[2][1]=1,A.a[2][2]=0;
		Matricx B=quick_pow(A,k-2);
		printf("%lld",(B.a[1][1]+B.a[1][2])%p);
	}
	return 0; 
}

反思与总结：

  • 要特判n=1或n=2的情况
  • 记得mod p

广义斐波那契数列

传送门: P1349 广义斐波那契数列
题目分析: 简单推一下:
{ f [ n ] = p × f [ n − 1 ] + q × f [ n − 2 ] f [ n − 1 ] = 1 × f [ n − 1 ] + 0 × f [ n − 2 ]    ⟹    [ f [ n ] f [ n − 1 ] ] = [ p q 1 0 ] ∗ [ f [ n − 1 ] f [ n − 2 ] ] \begin{cases} f[n]=p×f[n-1]+q×f[n-2] \\f[n-1]=1×f[n-1]+0×f[n-2] \end{cases}\implies\begin{bmatrix} f[n] \\ f[n-1] \end{bmatrix} =\begin{bmatrix} p & q \\ 1 & 0 \end{bmatrix}*\begin{bmatrix} f[n-1] \\ f[n-2] \end{bmatrix} {f[n]=p×f[n−1]+q×f[n−2]f[n−1]=1×f[n−1]+0×f[n−2]​⟹[f[n]f[n−1]​]=[p1​q0​]∗[f[n−1]f[n−2]​]
最后别忘了 f [ 1 ] = a 1 f[1]=a1 f[1]=a1, f [ 2 ] = a 2 f[2]=a2 f[2]=a2
Code:

#include<bits/stdc++.h>
using namespace std;
inline long long Read(){
	long long dx=0,fh=1;
	char c=getchar();
	while(c>'9'||c<'0'){
		if(c=='-')	fh=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		dx=dx*10+c-'0';
		c=getchar();
	}
	return dx*fh;
}
long long n,p,a2,a1,q,m;
struct Matricx{
	long long a[3][3];
}A;
Matricx mul(Matricx x,Matricx y){
	Matricx z;
	for(long long i=1;i<=2;++i)
		for(long long j=1;j<=2;++j){
			z.a[i][j]=0;
			for(long long u=1;u<=2;++u)
				z.a[i][j]=(z.a[i][j]+(x.a[i][u]*y.a[u][j])%m)%m;
		}
	return z;
}
Matricx quick_pow(Matricx x,long long b){
	Matricx r,base=x;
	for(int i=1;i<=2;++i)	for(int j=1;j<=2;++j)	r.a[i][j]=0;
	for(int i=1;i<=2;++i)	r.a[i][i]=1;
	while(b){
		if(b&1)	r=mul(r,base);
		base=mul(base,base);
		b>>=1;
	}
	return r;
}
int main(){
	p=Read();q=Read();a1=Read();a2=Read();n=Read();m=Read();
	if(n==1)	printf("%lld",a1);
	else if(n==2)	printf("%lld",a2);
	else{
		A.a[1][1]=p,A.a[1][2]=q,A.a[2][1]=1,A.a[2][2]=0;
		Matricx B=quick_pow(A,n-2);
		printf("%lld",((B.a[1][1]*a2)%m+(B.a[1][2]*a1)%m)%m);
	}
	return 0; 
}

【模板】矩阵加速(数列)

传送门:P1939 【模板】矩阵加速(数列)
算法分析:
简单分析一下:
{ f [ n ] = 1 × f [ n − 1 ] + 0 × f [ n − 2 ] + 1 × f [ n − 3 ] f [ n − 1 ] = 1 × f [ n − 1 ] + 0 × f [ n − 2 ] + 0 × f [ n − 3 ] f [ n − 2 ] = 0 × f [ n − 1 ] + 1 × f [ n − 2 ] + 0 × f [ n − 3 ]    ⟹    [ f [ n ] f [ n − 1 ] f [ n − 2 ] ] = [ 1 0 1 1 0 0 0 1 0 ] ∗ [ f [ n − 1 ] f [ n − 2 ] f [ n − 3 ] ]    ⟹    [ f [ n ] f [ n − 1 ] f [ n − 2 ] ] = [ 1 0 1 1 0 0 0 1 0 ] n − 3 × [ f [ 3 ] f [ 2 ] f [ 1 ] ]    ⟹    [ f [ n ] f [ n − 1 ] f [ n − 2 ] ] = [ 1 0 1 1 0 0 0 1 0 ] n − 3 × [ 1 1 1 ] \begin{cases} f[n]=1×f[n-1]+0×f[n-2]+1×f[n-3] \\f[n-1]=1×f[n-1]+0×f[n-2]+0×f[n-3] \\f[n-2]=0×f[n-1]+1×f[n-2]+0×f[n-3] \end{cases} \implies \begin{bmatrix} f[n]\\f[n-1]\\f[n-2] \end{bmatrix} =\begin{bmatrix} 1&0&1 \\ 1&0&0 \\0&1&0 \end{bmatrix} *\begin{bmatrix} f[n-1] \\ f[n-2]\\f[n-3] \end{bmatrix} \implies \begin{bmatrix} f[n]\\f[n-1]\\f[n-2] \end{bmatrix} =\begin{bmatrix} 1&0&1 \\ 1&0&0 \\0&1&0 \end{bmatrix}^{n-3}×\begin{bmatrix}f[3]\\f[2]\\f[1] \end{bmatrix}\implies \begin{bmatrix} f[n]\\f[n-1]\\f[n-2] \end{bmatrix} =\begin{bmatrix} 1&0&1 \\ 1&0&0 \\0&1&0 \end{bmatrix}^{n-3}×\begin{bmatrix}1\\1\\1 \end{bmatrix} ⎩⎪⎨⎪⎧​f[n]=1×f[n−1]+0×f[n−2]+1×f[n−3]f[n−1]=1×f[n−1]+0×f[n−2]+0×f[n−3]f[n−2]=0×f[n−1]+1×f[n−2]+0×f[n−3]​⟹⎣⎡​f[n]f[n−1]f[n−2]​⎦⎤​=⎣⎡​110​001​100​⎦⎤​∗⎣⎡​f[n−1]f[n−2]f[n−3]​⎦⎤​⟹⎣⎡​f[n]f[n−1]f[n−2]​⎦⎤​=⎣⎡​110​001​100​⎦⎤​n−3×⎣⎡​f[3]f[2]f[1]​⎦⎤​⟹⎣⎡​f[n]f[n−1]f[n−2]​⎦⎤​=⎣⎡​110​001​100​⎦⎤​n−3×⎣⎡​111​⎦⎤​
设 [ 1 0 1 1 0 0 0 1 0 ] \begin{bmatrix} 1&0&1 \\ 1&0&0 \\0&1&0 \end{bmatrix} ⎣⎡​110​001​100​⎦⎤​为 A A A，则 f [ n ] = A [ 1 ] [ 1 ] + A [ 1 ] [ 2 ] + A [ 1 ] [ 3 ] f[n]=A[1][1]+A[1][2]+A[1][3] f[n]=A[1][1]+A[1][2]+A[1][3]
Code:

#include<bits/stdc++.h>
using namespace std;
const long long p=1000000007;
inline long long Read(){
	long long dx=0,fh=1;
	char c=getchar();
	while(c>'9'||c<'0'){
		if(c=='-')	fh=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		dx=dx*10+c-'0';
		c=getchar();
	}
	return dx*fh;
}
long long k,T;
struct Matricx{
	long long a[4][4];
}A;
Matricx mul(Matricx x,Matricx y){
	Matricx z;
	for(long long i=1;i<=3;++i)
		for(long long j=1;j<=3;++j){
			z.a[i][j]=0;
			for(long long u=1;u<=3;++u)
				z.a[i][j]=(z.a[i][j]+(x.a[i][u]*y.a[u][j])%p)%p;
		}
	return z;
}
Matricx quick_pow(Matricx x,long long b){
	Matricx r,base=x;
	for(int i=1;i<=3;++i)	for(int j=1;j<=3;++j)	r.a[i][j]=0;
	for(int i=1;i<=3;++i)	r.a[i][i]=1;
	while(b){
		if(b&1)	r=mul(r,base);
		base=mul(base,base);
		b>>=1;
	}
	return r;
}
int main(){
	A.a[1][1]=1,A.a[1][2]=0,A.a[1][3]=1;
	A.a[2][1]=1,A.a[2][2]=0,A.a[2][3]=0;
	A.a[3][1]=0,A.a[3][2]=1,A.a[3][3]=0;
	T=Read();
	for(;T;--T){
		k=Read();
		if(k==1||k==2||k==3)	printf("1\n");
		else{
			Matricx B=quick_pow(A,k-3);
			printf("%lld\n",(B.a[1][1]+B.a[1][2]+B.a[1][3])%p);
		}
	}	
	return 0; 
}
Sun_QH_
关注
  • 3
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
牛客网 矩阵幂（矩阵快速幂、北邮机试）
sunlanchang的博客
02-20 310
题目描述 给定一个n*n的矩阵，求该矩阵的k次幂，即P^k。 输入描述 第一行：两个整数n（2&amp;amp;lt;=n&amp;amp;lt;=10）n（2&amp;amp;lt;=n&amp;amp;lt;=10）n（2&amp;lt;=n&amp;lt;=10）、k（1&amp;amp;lt;=k&amp;amp;lt;=5）k（1&amp;amp;lt;=k&amp;amp;lt;=5）k（1&amp;lt;=
矩阵快速幂刷题系列
weixin_30703911的博客
08-28 163
来源自http://blog.csdn.net/chenguolinblog/article/details/10309423 hdu 1575 Tr A Time Limit: 1000/1000 MS (Java/Others)Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5587Accepte...
矩阵快速幂：数列考验
qq_44654498的博客
05-24 104
Vjudge:矩阵快速幂板子题
矩阵幂
Lassi的博客
03-19 154
题目描述： 给定一个n*n的矩阵，求该矩阵的k次幂，即P^k。 输入描述： 第一行：两个整数n（2<=n<=10）、k（1<=k<=5），两个数字之间用一个空格隔开，含义如上所示。 接下来有n行，每行n个正整数，其中，第i行第j个整数表示矩阵中第i行第j列的矩阵元素Pij且（0<=Pij<=10）。另外，数据保证最后结果不会超过10^8。 输出描述： 对于每组测试数据，输出其结果。格式为： n行n列个整数，每行数之间用空格隔开，注意，每行最后一个数后面不应该有多余的空格。
矩阵快速幂-Fast Matrix Calculation HDU - 4965
njuptACMcxk的博客
02-03 163
矩阵快速幂-Fast Matrix Calculation HDU - 4965 题目： One day, Alice and Bob felt bored again, Bob knows Alice is a girl who loves math and is just learning something about matrix, so he decided to make a craz...
python进行爬虫小记
01-15
Python爬虫技术是一种用于自动化网页数据抓取的编程方法，尤其适合初学者快速入门。Python在爬虫领域具有显著优势，因为其拥有丰富的第三方库，如requests、lxml和parsel等，使得编写爬虫代码变得简洁高效。此外，...
Beatles小记
03-02
这篇小记主要处于两方面考虑：首先，希望打破一提到海量数据分析，就只有hadoop基础上的一系列工具，更多的时候很多企业需要的是更轻量的设计（办喜酒杀猪杀鸡未必都要用一把刀），因此将开放平台基础分析组件重构...
Scrum过程实践小记
02-26
严格来说，不能算是真正的scrum实践，但实践敏捷的过程本身也是一种“敏捷方法”，所以就算是“敏捷实践之敏捷开发方法－scrum过程”吧。1.Scrum团队（5-7个人的小项目组）。　2.Backlog:急待完成的一系列任务，包括...
快速幂+矩阵快速幂（总结+例题）
zhoupingqi2017的博客
04-02 153
转载链接：https://blog.csdn.net/red_red_red/article/details/90208713
算法之矩阵快速幂
zkc_home
03-15 365
矩阵快速幂 矩阵的快速幂是用来高效地计算矩阵的高次方的。将朴素的o（n）的时间复杂度，降到log（n）。 这里先对原理（主要运用了矩阵乘法的结合律）做下简单形象的介绍： 一般一个矩阵的n次方，我们会通过连乘n-1次来得到它的n次幂。 但做下简单的改进就能减少连乘的次数，方法如下： 把n个矩阵进行两两分组，比如：A*A*A*A*A*A =&gt; (A*A)(A*A)(A*A) 这样...
矩阵快速幂总结
weixin_30484739的博客
07-23 88
依然主要还是自用 首先矩阵一条性质的概述和证明 概述:对于一个临接矩阵$G$来说，它自乘$G^k$次方中$G[i][j]$含义为从i走到j走k步方案数. 证明:比较麻烦,我们设f[i]表示从1走到其他点方案数,那么根据矩阵递推优化,k次就是转移了k次,那么每次走一步,k次就是走了k步 最好手模一下,然后自己随便写一个数据验证一下. 然后矩阵可以优化递推,比如优化菲波纳契数列的递...
矩阵快速幂（附模板）
Running Snail
08-30 379
求解矩阵 A 的 N 次方，我们可以类比整数快速幂，写一个矩阵的结构体，用一个matmul函数来定义矩阵的乘法，具体实现过程与整数快速幂类似（整数快速幂） 模板 struct mat { ll m[maxn][maxn]; }unit; void init() { for(int i=1;i<maxn;i++) unit.m[i][i]=1; } mat matmul(mat a,mat b)//ans=矩阵a*矩阵b { mat ans; ll tmp
快速幂、矩阵快速幂及其模板
whx1003__的博客
05-18 163
题目描述 给定333个整数a,b,pa,b,pa,b,p，求ab mod pa^b\ mod\ pab mod p的值。 输入输出格式 输入格式： 333个整数a,b,pa,b,pa,b,p。 输出格式： ab mod pa^b\ mod\ pab mod p的值。 输入输出样例 输入样例#1： 2 3 7 输出样...
简单计数(递推式 矩阵快速幂加速)
ccsu_deer
05-30 642
链接：https://ac.nowcoder.com/acm/contest/879/A 来源：牛客网 简单计数 时间限制：C/C++ 1秒，其他语言2秒 空间限制：C/C++ 524288K，其他语言1048576K 64bit IO Format: %lld 题目描述 thewindandrainisinmywayandnevergoingawaythewindandrainisi...
快速幂算法及矩阵快速幂
winycg的博客
03-06 886
Description 给定三个数A, B, K, 求 A的B次方除以K的余数 。 Input 输入只有一行，为三个正整数A(1 Output 一个整数，(A ^ B) % K 的值。 Sample Input 11 100 7 Sample Output 4 #include using namespace std; int main() { long long
pinocchio计算MPC控制状态方程AB矩阵
xiaoc100200的博客
06-16 63
Pinocchio是一个用于机器人动力学和运动学计算的库，广泛应用于机器人控制和仿真中。为了在MPC控制中使用Pinocchio计算状态方程的线性化矩阵 A 和 B
cryptography-0.9-cp34-none-win32.whl
06-18
Python库是一组预先编写的代码模块，旨在帮助开发者实现特定的编程任务，无需从零开始编写代码。这些库可以包括各种功能，如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库，如NumPy、Pandas和Requests，极大地丰富了Python的应用领域，从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径，而且为经验丰富的开发者提供了强大的工具，以高效率、高质量地完成复杂任务。例如，Matplotlib和Seaborn库在数据可视化领域内非常受欢迎，它们提供了广泛的工具和技术，可以创建高度定制化的图表和图形，帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
antdpro Demo
最新发布
06-18
antdpro Demo
protobuff使用小记
05-31
下面是使用Protobuf的一些小记： 1. 定义消息格式 首先，需要定义消息格式，以便Protobuf可以将数据序列化和反序列化。消息格式定义在.proto文件中，使用protobuf语言编写。例如，下面是一个简单的消息格式定义： ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值