矩阵快速幂(运算符重载)
2022/7/31 6:22:57
本文主要是介绍矩阵快速幂(运算符重载),对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
https://www.luogu.com.cn/problem/P3390
- 把*重载成矩阵的乘法
- 再用普通的快速幂就行
- (AC代码是copy的,实在debug不出了)
#include <algorithm> #include <iostream> #include <cstring> #include <cstdio> #include <cctype> #define ll long long #define gc() getchar() #define maxn 105 #define mo 1000000007 using namespace std; inline ll read() { ll a = 0; int f = 0; char p = gc(); while (!isdigit(p)) { f |= p == '-'; p = gc(); } while (isdigit(p)) { a = (a << 3) + (a << 1) + (p ^ 48); p = gc(); } return f ? -a : a; } int n; struct ahaha { ll a[maxn][maxn]; ahaha() { memset(a, 0, sizeof a); } inline void build() { //建造单位矩阵 for (int i = 1; i <= n; ++i) a[i][i] = 1; } } a; ahaha operator*(const ahaha &x, const ahaha &y) { ahaha z; for (int k = 1; k <= n; ++k) for (int i = 1; i <= n; ++i) for (int j = 1; j <= n; ++j) z.a[i][j] = (z.a[i][j] + x.a[i][k] * y.a[k][j] % mo) % mo; return z; } ll k; inline void init() { n = read(); k = read(); for (int i = 1; i <= n; ++i) for (int j = 1; j <= n; ++j) a.a[i][j] = read(); } int main() { init(); ahaha ans; ans.build(); do { if (k & 1) ans = ans * a; a = a * a; k >>= 1; } while (k); for (int i = 1; i <= n; putchar('\n'), ++i) for (int j = 1; j <= n; ++j) printf("%d ", ans.a[i][j]); return 0; }
这篇关于矩阵快速幂(运算符重载)的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-13TiDB + ES:转转业财系统亿级数据存储优化实践
- 2024-05-09“2024鸿蒙零基础快速实战-仿抖音App开发(ArkTS版)”实战课程已上线
- 2024-05-09聊聊如何通过arthas-tunnel-server来远程管理所有需要arthas监控的应用
- 2024-05-09log4j2这么配就对了
- 2024-05-09nginx修改Content-Type
- 2024-05-09Redis多数据源,看这篇就够了
- 2024-05-09Google Chrome驱动程序 124.0.6367.62(正式版本)去哪下载?
- 2024-05-09有没有大佬知道这种数据应该怎么抓取呀?
- 2024-05-09这种运行结果里的10.100000001,怎么能最快改成10.1?
- 2024-05-09企业src漏洞挖掘-有意思的命令执行