【专题】矩阵快速幂

最新推荐文章于 2022-08-09 17:47:38 发布
最新推荐文章于 2022-08-09 17:47:38 发布
阅读量178 收藏
点赞数
分类专栏： ————数论相关———— # 矩阵快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Zy_Ming/article/details/89154502
版权

HDU-1757

思路：

这题是矩阵快速幂模板题。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int k, mod;
int num[10];
struct matrix{
	int m[10][10];
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	void init_1(){
		init_0();
		for(int i = 0; i < 10; ++i)
			m[i][i] = 1;
	} 
	void init(){
		init_0();
		m[0][0] = num[0], m[0][1] = num[1], m[0][2] = num[2], m[0][3] = num[3], m[0][4] = num[4], m[0][5] = num[5], m[0][6] = num[6], m[0][7] = num[7], m[0][8] = num[8], m[0][9] = num[9];
		m[1][0] = 1;
		m[2][1] = 1;
		m[3][2] = 1;
		m[4][3] = 1;
		m[5][4] = 1;
		m[6][5] = 1;
		m[7][6] = 1;
		m[8][7] = 1;
		m[9][8] = 1;
	}
	
	matrix operator*(const matrix &M)const{
		matrix temp;
		temp.init_0();
		for(int i = 0; i < 10; ++i){
			for(int j = 0; j < 10; ++j){
				for(int k = 0; k < 10; k++){
					temp.m[i][j] = (temp.m[i][j] + (m[i][k] * M.m[k][j])%mod) % mod;
				}
			}
		}
		return temp;
	}
};

matrix qpow(int t){
	matrix ans, p;
	ans.init_1();
//	cout << "ans init complete" << endl;
	p.init();
//	cout << "enter qpow" << endl;
	while(t){
		if(t&1)ans = ans * p;
		p = p * p;
		t >>= 1;
	}
	return ans;
}

int main(){
	int init[10] = {9,8,7,6,5,4,3,2,1,0};
	while(cin >> k >> mod){
		for(int i = 0; i < 10; ++i){
			cin >> num[i];
		} 
//		cout << "input complete" << endl;
		if(k>=10){
//			cout << "k>=10" << endl;
			matrix ans = qpow(k-9);
			ll a = 0;
			for(int i = 0; i < 10; ++i){
				a = a%mod + (ans.m[0][i]%mod) * (init[i]%mod);
			}
			cout << a%mod << endl;
		}
		else{
			cout << k%mod << endl;
		}
	}
	return 0;
}

 

HDU-1575

思路：

矩阵快速幂模板题。

求完快速幂后把主对角线的加起来就是答案。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 9973;
int t, n, k;
int a[11][11];
struct matrix{
	int m[11][11];
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	void init_1(){
		init_0();
		for(int i = 1; i <= n; ++i){
			m[i][i] = 1;
		}
	}
	void init(){
		for(int i = 1; i <= n; ++i){
			for(int j = 1; j <= n; ++j){
				m[i][j] = a[i][j];
			}
		}
	}
	matrix operator*(matrix &M)const{
		matrix temp;
		temp.init_0();
		for(int i = 1; i <= n ; ++i){
			for(int j = 1; j <= n; ++j){
				for(int k = 1; k <= n; ++k){
					temp.m[i][j] = (temp.m[i][j] + (m[i][k] * M.m[k][j])%mod)%mod;
				}
			}
		}
		return temp;
	} 
};

matrix qpow(ll k){
	matrix ans, p;
	ans.init_1();
	p.init();
	while(k){
		if(k&1) ans = ans * p;
		p = p * p;
		k >>= 1;
	}
	return ans;
}

int main(){
	cin >> t;
	while(t--){
		cin >> n >> k;
		for(int i = 1; i <= n; ++i){
			for(int j = 1; j <= n; ++j){
				cin >> a[i][j];
			}
		}
		matrix ans = qpow(k);
		ll sum = 0;
		for(int i = 1; i <= n; ++i){
			sum += ans.m[i][i];
		}
		cout << sum%mod << endl;
	}
	return 0;
}

 

HDU-2604

思路：

算出前面几个值后可以看出 n\geq 4 时，F(N)=F(n-1)+F(n-3)+F(N-4)。

根据公式就可以写出转移矩阵了。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int n, mod;
struct matrix{
	ll m[11][11];
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	void init_1(){
		init_0();
		for(int i = 1; i <= 4; ++i){
			m[i][i] = 1;
		}
	}
	void init(){
		init_0();
		m[1][1] = 1, m[1][3] = 1, m[1][4] = 1;
		m[2][1] = 1;
		m[3][2] = 1;
		m[4][3] = 1; 
	}
	matrix operator*(matrix &M)const{
		matrix temp;
		temp.init_0();
		for(int i = 1; i <= 4; ++i){
			for(int j = 1; j <= 4; ++j){
				for(int k = 1; k <= 4; ++k){
					temp.m[i][j] = (temp.m[i][j] + (m[i][k] * M.m[k][j])%mod)%mod;
				}
			}
		}
		return temp;
	} 
};

matrix qpow(ll k){
	matrix ans, p;
	ans.init_1();
	p.init();
	while(k){
		if(k&1) ans = ans * p;
		p = p * p;
		k >>= 1;
	}
	return ans;
}

int main(){
	int a[5] = {0, 1, 2, 4, 6};
	while(cin >> n >> mod){
		if(n<=3){
			cout << a[n+1]%mod << endl;;
		}
		else {
			ll sum = 0;
			matrix ans = qpow(n-3);
			for(int i = 1; i <= 4; ++i){
				sum += ans.m[1][i] * a[4-i+1] % mod;
			}
			cout << sum % mod << endl;
		}
	}
	return 0;
}

 

HDU-2256

思路：

代码：

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int mod = 1024;
int t;
ll n;

struct matrix{
	int m[3][3];
	
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	
	void init_1(){
		init_0();
		for(int i = 1; i <= 2; ++i){
			m[i][i] = 1;
		}
	}
	
	void init(){
		init_0();
		m[1][1] = 5, m[1][2] = 12;
		m[2][1] = 2, m[2][2] = 5;
	}
	
	matrix operator*(const matrix &M)const{
		matrix ans;
		ans.init_0();
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				for(int k = 1; k <= 2; ++k){
					ans.m[i][j] = (ans.m[i][j] + (m[i][k] * M.m[k][j]) % mod) % mod;
				}
			}
		}
		return ans;
	}
};
matrix a;
ll qpow(ll T){
	matrix ans, p;
	ans.init_1();
	p.init();
	while(T){
		if(T&1)ans = ans * p;
		p = p * p;
		T >>= 1;
	}
	ans = ans * a;
	return 2*ans.m[1][1] - 1;
}

int main(){
	scanf("%d", &t);
	a.m[1][1] = 5, a.m[2][1] = 2;
	while(t--){
		scanf("%lld", &n);
		if(n==1)printf("9\n");
		else printf("%lld\n", qpow(n-1)%1024);
	}
	return 0;
}

 

CodeForces - 185A

思路：

根据观察，一个上三角形生成 3 个小的上三角形；一个下三角形生成 1 个上三角形。

推出公式 F(N)=2*F(N-1)+4^{N-1}(N\geq 1)。

代码：

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int mod = 1000000007;
ll n;

struct matrix{
	ll m[3][3];
	
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	
	void init_1(){
		init_0();
		for(int i = 1; i <= 2; ++i){
			m[i][i] = 1;
		}
	}
	
	void init(){
		init_0();
		m[1][1] = 2, m[1][2] = 4;
		m[2][1] = 0, m[2][2] = 4;
	}
	
	matrix operator*(const matrix &M)const{
		matrix ans;
		ans.init_0();
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				for(int k = 1; k <= 2; ++k){
					ans.m[i][j] = (ans.m[i][j] + (m[i][k] * M.m[k][j]) % mod) % mod;
				}
			}
		}
		return ans;
	}
};
matrix a;
ll qpow(ll T){
	matrix ans, p;
	ans.init_1();
	p.init();
	while(T){
		if(T&1)ans = ans * p;
		p = p * p;
		T >>= 1;
	}
	ans = ans * a;
	return ans.m[1][1];
}

int main(){
	a.m[1][1] = 3, a.m[2][1] = 1;
	scanf("%I64d", &n);
	if(n==0)printf("1\n");
	else if(n==1)printf("3\n");
	else printf("%I64d\n", qpow(n-1)%mod);
	return 0;
}

 

HDU - 2842

思路：

首先，明确游戏规则：前 k 个环取走，第 k+1 还在时，可以取走第 k+2 个环。

例如，有 1、2、3、4、5、6 共 6 个环。

当前剩下 4、5、6，这时可以取走 5 号环。

即，要取走前 k 个环(F(K))，要取走前 K-2 个环(F(K-2))后，拿走第 K 个环，再把前 K-2 个环(F(K-2))放回去，再执行拿走前 K-1 个环(F(K-1))，最后再补上拿走第 1 个环的 1 次。

以上就可以得出递推式：F(K)=F(k-1)+2*F(k-2)+1 (K>2)。

之后就可以得出转移矩阵

\begin{bmatrix} 1& 2 &1 \\ 1& 0 &0 \\ 0& 0 &1 \end{bmatrix}

初始矩阵

\begin{bmatrix} 2& 0 &0 \\ 1& 0 &0 \\ 1& 0 &0 \end{bmatrix}

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 200907;
ll n;

struct matrix{
	ll m[5][5];
	
	void init_0(){
		memset(m, 0, sizeof(m));
	}
	void init_1(){
		init_0();
		for(int i = 1; i <= 3; ++i)m[i][i] = 1;
	}
	void init(){
		init_0();
		m[1][1] = 1, m[1][2] = 2, m[1][3] = 1;
		m[2][1] = 1;
		m[3][3] = 1;
	}
	matrix operator*(const matrix &M)const{
		matrix ans;
		ans.init_0();
		for(int i = 1; i <= 3; ++i)
			for(int j = 1; j <= 3; ++j)
				for(int k = 1; k <= 3; ++k)
					ans.m[i][j] = (ans.m[i][j] + (m[i][k] * M.m[k][j]) % mod) % mod;
		return ans;
	}
}; 
matrix a;
ll qpow(ll t){
	matrix ans, p;
	ans.init_1();
	p.init();
	while(t){
		if(t&1) ans = ans * p;
		p = p * p;
		t >>= 1;
	}
	ans = ans * a;
	return ans.m[1][1]%mod;
}

int main(){
	a.m[1][1] = 2, a.m[2][1] = 1, a.m[3][1] = 1;
	while(scanf("%lld", &n)!=EOF, n){
		if(n>=3){
			printf("%lld\n", qpow(n-2));
		}
		else{
			printf("%lld\n", n);
		}
	}
	return 0;
}

未完待续。。。

Zy_Ming
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
HDU 5015 233 Matrix （数论——矩阵快速幂）
SHUI
08-05 686
233 Matrix Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 0    Accepted Submission(s): 0 Problem Description In our daily life we
矩阵快速幂专题（矩阵快速幂入门、矩阵构造法、数论规律题）
林下的码路
08-04 2324
Link1：http://acm.hdu.edu.cn/showproblem.php?pid=1575 Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3628    Accepted Submissio
数论基础——循环节和矩阵快速幂的运用
多反思，多回顾，要坚持。
05-03 3970
首先我们来看一道基础题：     题目链接：HDU1005 Number Sequence     题目描述： Number Sequence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 147421    Acce
[算法]矩阵快速幂
wulalalawu的博客
01-03 147
注意：如果有n的k次方，那就把n的k次方写成k个n相乘，然后补k行k列就好了。 求(a_n = a_{n – 1} + 2 * a_{n – 2} + n^3) 首先我们需要知道立方和完全展开式：((a – b)^3 = a^3 – 3 a^2 b + 3 a b^2 – b^3) 然后配就行了((i + 1)^3 = i^3 + 3 i^2 + 3 i + 1)。 有一个要注意的地方是，最后要拿转移矩阵的n次方的第一行乘上初始列才是答案。 （推完之后发现这矩阵好像杨辉三角形，如果有更高次项应该能按杨辉三角
矩阵快速幂专题
weixin_43184669的博客
08-26 117
poj3233 经典题 因为k的值太大，如果直接算，肯定会超时。 因此采用二分的方法。 #include <cstdio> #include <cstring> using namespace std; int n,k,m; struct matrix { int mat[35][35]; }; matrix A; void init(matrix &a) ...
矩阵小专题（矩阵快速幂+矩阵加速）
bunny_1024的博客
08-09 340
矩阵小专题（矩阵快速幂+矩阵加速）
20191010 专题：矩阵快速幂
hongkongreporter
10-18 158
重载运算符 struct gg { int a[n][n]; friend operator * (const gg &a,const gg *b) { gg c; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { int num=0; for(int k=1;k<=n;k++) ...
矩阵快速幂专题【完结】
你那么爱笑，家里人知道吗？
09-03 375
第一题 hdu 1757 A Simple Math Problem 点击打开链接 思路：矩阵快速幂 分析: 1 最简单的矩阵快速幂的题目，直接利用矩阵求解即可 点击打开查看代码 第二题 hdu 1575 Tr A 点击打开hdu 1575 思路: 矩阵快速幂 分析: 1 题目给定一个n*n的矩阵要求矩阵的
MATLAB学习专题.
05-13
本学习专题旨在深入解析MATLAB的基础知识和常用功能，帮助初学者快速掌握这一工具。 1. MATLAB的特点： MATLAB 提供了交互式环境，支持矩阵运算，具有丰富的内置函数库，易于编程，并能与各种外部数据源进行接口。...
dp专题讲义1
08-03
3. 减少状态转移花费的时间：通过预处理或数据结构优化，加快每次状态转移的计算速度，例如使用矩阵快速幂、前缀和等技术。 动态规划的应用需要满足两个前提条件：最优化原理和无后效性。最优化原理指出，无论过去...
第一周数学专题题解1
08-03
8. **矩阵快速幂**：H题中，要求递推序列的第N项，使用矩阵快速幂可以高效求解，尤其在存在循环节的情况下。 9. **模运算与同余方程**：I题是基于模运算的同余方程求解，可以通过欧几里得算法和枚举法找出满足条件...
动态规划专题课件
12-06
例如，链式矩阵快速幂、记忆化搜索等高级技巧在解决复杂问题时能显著提高效率。 五、动态规划的优化技巧 1. 记忆化搜索：通过存储已经计算过的子问题结果，避免重复计算。 2. 自底向上：从简单的情况开始计算，逐渐...
2012版高考数学 3-2-1精品系列 专题16 矩阵与变换、行列式
08-19
4. **二阶矩阵与二元一次方程组**：矩阵可以用来表示线性方程组，系数矩阵的逆矩阵可以帮助我们快速求解方程组。线性方程组的解的存在性和唯一性是由矩阵的秩决定的。 5. **特征值与特征向量**：特征值和特征向量是...
【专题】莫比乌斯反演
Zy_Ming的博客
04-16 1483
问题引入： 给定整数 N 和 M。求满足且为质数的点对的个数。 数据范围： 接下来会见到以下内容： 莫比乌斯函数 莫比乌斯函数的线性筛 迪利克雷卷积介绍 莫比乌斯反演 整除分块 杜教筛介绍 莫比乌斯函数： 这里 else 是指：n有大于1的平方因子的情况，如，4、9、16等。 莫比乌斯函数的线性筛： 其实，莫比乌斯函数线性筛与普通的线性...
【LightOJ】1179 Josephus Problem - 约瑟夫环
Zy_Ming的博客
12-06 405
G - Josephus Problem  LightOJ - 1179  The historian Flavius Josephus relates how, in the Romano-Jewish conflict of 67 A.D., the Romans took the town of Jotapata which he was commanding. Escaping, Jo...
【51nod】1181 质数中的质数 - 质数筛法
Zy_Ming的博客
08-24 402
1181 质数中的质数（质数筛法）     题解： 易知，题目要我们求质数。 把求得的质数存放于数组 prime[MAXN] 中,同时用 tag[MAXN] 记录质数与非质数。 tag[i] == 0 为质数，tag[i] == 1 为非质数。 接着循环判断 prime[i] &gt;= n &amp;&amp; tag[i] == 0 即可。 为了提高效率，这里采用了线性筛法。 ...
【51nod】1008 N的阶乘 mod P - 同余定理
Zy_Ming的博客
08-24 334
1008 N的阶乘 mod P   题解： 这里用到了同余定理 (a+b)%m = (a%m + b%m) % m (a*b)%m = (a%m * b%m) % m  #include &lt;cstdio&gt; using namespace std; long long Solve(int n, int m){ if(n == 1) return 1...
【hdu】6470 Count - 矩阵快速幂
Zy_Ming的博客
03-29 268
Count Time Limit: 6000/3000 MS (Java/Others)Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 708Accepted Submission(s): 279 Problem Description Farmer John有n头奶牛. 某天奶牛想要数一数有多少头...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值