矩阵快速幂
2021/12/4 6:20:30
本文主要是介绍矩阵快速幂,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
- 初始化单元矩阵
- 输入矩阵数值
- 类似于快速幂求解矩阵乘法
- 重载运算法*,用来表示矩阵乘法
#include<bits/stdc++.h> using namespace std; const int maxn=110; const int mod=1e9+7; struct mat{ int m[maxn][maxn]; }unit; int n; void init_unit() { for(int i=0;i<maxn;i++) unit.m[i][i]=1; return ; } mat operator *(mat a,mat b) { mat ret; long long x; for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ x=0; for(int k=0;k<n;k++) x=(x%mod+((long long)a.m[i][k]*b.m[k][j])%mod)%mod; ret.m[i][j]=x%mod; } } return ret; } mat pow_mat(mat a,long long y) { mat res=unit; while(y) { if(y&1) res=res*a; y>>=1; a=a*a; } return res; } int main() { long long x; init_unit(); while(cin>>n>>x) { mat a; for(int i=0;i<n;i++) for(int j=0;j<n;j++) cin>>a.m[i][j]; a=pow_mat(a,x); for(int i=0;i<n;i++) for(int j=0;j<n;j++) if(j+1==n) cout<<a.m[i][j]<<endl; else cout<<a.m[i][j]<<" "; } return 0; }
这篇关于矩阵快速幂的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-09-28AI给的和自己写的Python代码,都无法改变输入框的内容,替换也不行
- 2024-09-27Sentinel配置限流资料:新手入门教程
- 2024-09-27Sentinel配置限流资料详解
- 2024-09-27Sentinel限流资料:新手入门教程
- 2024-09-26Sentinel限流资料入门详解
- 2024-09-26Springboot框架资料:初学者入门教程
- 2024-09-26Springboot框架资料详解:新手入门教程
- 2024-09-26Springboot企业级开发资料:新手入门指南
- 2024-09-26SpringBoot企业级开发资料新手指南
- 2024-09-26Springboot微服务资料入门教程