矩阵快速幂刷题小记

最新推荐文章于 2024-06-18 16:45:27 发布
最新推荐文章于 2024-06-18 16:45:27 发布
阅读量115 收藏
点赞数 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;=
python进行爬虫小记
01-15
Python爬虫技术是一种用于自动化网页数据抓取的编程方法，尤其适合初学者快速入门。Python在爬虫领域具有显著优势，因为其拥有丰富的第三方库，如requests、lxml和parsel等，使得编写爬虫代码变得简洁高效。此外，...
矩阵快速幂刷题系列
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 155
题目描述： 给定一个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...
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
矩阵补全学习笔记
jacke121的专栏
06-17 250
矩阵补全学习笔记
短视频矩阵系统：高效运营，解决多账号管理难题
最新发布
qq_66975719的博客
06-18 592
这款短视频矩阵系统以其全面的功能和高效的运营模式，真正实现了多账号管理的便捷与高效，无疑是当下短视频运营者的得力助手。它不仅提升了视频创作的效率，更在保障内容质量的同时，助力运营者轻松应对多平台挑战，实现短视频内容的最大化传播。这样的系统，确实值得每一位追求高效的短视频运营者推荐与尝试。
重构大学数学基础_week05_雅各比矩阵与雅各比行列式
https://github.com/foxpup11?tab=repositories
06-17 461
https://www.bilibili.com/video/BV11a4y1L7BV?p=69&vd_source=7dace3632125a1ef7fd32c285eb2fbac
protobuff使用小记
05-31
下面是使用Protobuf的一些小记： 1. 定义消息格式 首先，需要定义消息格式，以便Protobuf可以将数据序列化和反序列化。消息格式定义在.proto文件中，使用protobuf语言编写。例如，下面是一个简单的消息格式定义： ...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值