51nodoj 1113 矩阵快速幂
题目连接
基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题
收藏
关注
给出一个N * N的矩阵,其中的元素均为正整数。求这个矩阵的M次方。由于M次方的计算结果太大,只需要输出每个元素Mod (10^9 + 7)的结果。
Input
第1行:2个数N和M,中间用空格分隔。N为矩阵的大小,M为M次方。(2 <= N <= 100, 1 <= M <= 10^9)
第2 - N + 1行:每行N个数,对应N * N矩阵中的1行。(0 <= N[i] <= 10^9)
Output
共N行,每行N个数,对应M次方Mod (10^9 + 7)的结果。
Input示例
2 3
1 1
1 1
Output示例
4 4
4 4
这个题是用operator重载一下矩阵的乘法运算,然后按照快速幂的方法来求解
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
#define M 110
#define LL long long
#define MOD 1000000007struct node
{LL m[M][M];
} unit;int n;
void init()//单位矩阵
{for(int i=0; i<M; i++){unit.m[i][i] = 1;}
}
node operator * (node t1, node t2)//重载一下乘法
{node res;LL 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 += (t1.m[i][k] * t2.m[k][j]) % MOD;x %= MOD;}res.m[i][j] = x;}}return res;
}
node quick_pow(node x, int m)
{node res = unit;//单位矩阵while(m){if(m&1){res = res * x;}x = x * x;m >>= 1;}return res;
}
int main()
{int m;init();node a = unit;scanf("%d%d", &n, &m);for(int i=0; i<n; i++){for(int j=0; j<n; j++){scanf("%d", &a.m[i][j]);}}a = quick_pow(a, m);for(int i=0; i<n; i++){for(int j=0; j<n; j++){printf("%d%c", a.m[i][j], j==n-1 ? '\n' : ' ');}}return 0;
}
51nodoj 1113 矩阵快速幂相关推荐
- 51nod 1113 矩阵快速幂
1113 矩阵快速幂 基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大 ...
- 51nod 1113 矩阵快速幂 模板题
1113 矩阵快速幂 基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 收藏 关注 给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计 ...
- 51nod 1113 矩阵快速幂【裸题】【内含黑科技】
1113 矩阵快速幂 基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大 ...
- 51nod 1113 矩阵快速幂 (矩阵的n次方)
给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大,只需要输出每个元素Mod (10^9 + 7)的结果. Input 第1行:2个数N和M,中间用空格分隔. ...
- 51nod 1113 矩阵快速幂 模板
给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大,只需要输出每个元素Mod (10^9 + 7)的结果. Input 第1行:2个数N和M,中间用空格分隔. ...
- 51nod-1113 矩阵快速幂
原题链接 1113 矩阵快速幂 基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 收藏 关注 给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方 ...
- 矩阵快速幂(51nod)
1113 矩阵快速幂 基准时间限制:3 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 给出一个N * N的矩阵,其中的元素均为正整数.求这个矩阵的M次方.由于M次方的计算结果太大, ...
- 矩阵快速幂+构造方法
与快速幂一样,可以将递推式通过二进制的方式来进行优化,这个学了快速幂就是十分容易理解 大概的板子如下: struct mat///自己定义大小的矩阵 {ll m[11][11]; }; mat mul ...
- 【做题】SRM701 Div1 Hard - FibonacciStringSum——数学和式&矩阵快速幂
原文链接 https://www.cnblogs.com/cly-none/p/SRM701Div1C.html 题意:定义"Fibonacci string"为没有连续1的01串 ...
最新文章
- [原创]java WEB学习笔记94:Hibernate学习之路---session 的管理,Session 对象的生命周期与本地线程绑定...
- Apache JMeter 字体、字体大小修改
- 绑定事件的几个方法总结
- Java 使用SAX解析XML文档
- HDU 1027 G - Can you answer these queries?
- @程序员,如何在五分钟内构建个人网盘服务?| 技术头条
- java构造方法可以是私有的吗_Java中的构造函数可以私有吗?
- 1.1.15 word调整文字与下划线之间的间距
- 祥福别园GhostXP_SP3纯净自选版_2013.04[NTFS格式]
- 力软敏捷开发框架源码7.0.6旗舰版
- 手机app性能测试简介了解
- 详解推荐系统的算法与应用
- 一篇想要获取积分的博
- 迅雷和BT有什么区别?迅雷是不是不毁硬盘?速度快吗?
- [Java安全]—fastjson漏洞利用
- html在浏览器显示图片,html - 在所有Web浏览器中显示TIFF图像
- Win10自动更新后,输入法故障问题
- 拓臻生物任命Senthil Sundaram为首席执行官、Mark Vignola, Ph.D.为首席财务官、Erin Quirk, M.D.为总裁
- Matlab:实现光波偏振态仿真
- 免费获取论文的网站,获取免费国外书籍的网站,SCI-Hub科研论文全文下载可用网址