UVA 11149.Power of Matrix-矩阵快速幂倍增
Power of Matrix
UVA - 11149
代码:
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 #define maxn 44 7 #define mod 10 8 int n; 9 struct matrix{ 10 int f[maxn][maxn]; 11 }; 12 matrix sum(matrix a,matrix b){ 13 int i,j; 14 matrix s; 15 for(i=0;i<n;i++) 16 for(j=0;j<n;j++) 17 s.f[i][j]=(a.f[i][j]+b.f[i][j])%mod; 18 return s; 19 } 20 matrix mul(matrix a,matrix b){ 21 int i,j,k; 22 matrix s; 23 memset(s.f,0,sizeof(s.f)); 24 for(i=0;i<n;i++) 25 for(j=0;j<n;j++) 26 for(k=0;k<n;k++) 27 s.f[i][j]=(s.f[i][j]+a.f[i][k]*b.f[k][j])%mod; 28 return s; 29 } 30 matrix pows(matrix e,int b){ 31 matrix s,a; 32 int i,j,k; 33 a=e; 34 for(i=0;i<n;i++) 35 for(j=0;j<n;j++) 36 if(i==j)s.f[i][j]=1; 37 else s.f[i][j]=0; 38 while(b){ 39 if(b&1) 40 s=mul(s,a); 41 a=mul(a,a); 42 b=b>>1; 43 } 44 return s; 45 } 46 matrix work(matrix e,int k){ 47 matrix s,a,b; 48 if(k==1) 49 return e; 50 a=work(e,k/2); 51 s=sum(a,mul(a,pows(e,k/2))); 52 if(k&1) 53 s=sum(s,pows(e,k)); 54 return s; 55 } 56 int main(){ 57 int k; 58 while( cin>>n>>k){ 59 if(n==0) 60 break; 61 int i,j,a; 62 matrix e; 63 for(i=0;i<n;i++) 64 { 65 for(j=0;j<n;j++) 66 { 67 cin>>a; 68 e.f[i][j]=a%10; 69 } 70 } 71 e=work(e,k); 72 for(i=0;i<n;i++) 73 { 74 cout<<e.f[i][0]; 75 for(j=1;j<n;j++) 76 cout<<" "<<e.f[i][j]; 77 cout<<endl; 78 } 79 cout<<endl; 80 } 81 return 0; 82 }
转载于:https://www.cnblogs.com/ZERO-/p/9692425.html
UVA 11149.Power of Matrix-矩阵快速幂倍增相关推荐
- BZOJ.4180.字符串计数(后缀自动机 二分 矩阵快速幂/倍增Floyd)
题目链接 先考虑 假设S确定,使构造S操作次数最小的方案应是:对T建SAM,S在SAM上匹配,如果有S的转移就转移,否则操作数++,回到根节点继续匹配S.即每次操作一定是一次极大匹配. 简单证明:假设 ...
- POJ 3233 Matrix Power Serie (矩阵快速幂)
Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + - + Ak. ...
- HDU - 5015 233 Matrix(矩阵快速幂)
题目链接:点击查看 题目大意:初始化:第一行依次为233,2333,23333....第一列依次为a0,a1,a2....(题目中会给出),再给出递推公式:,求矩阵中第n行m列的数字是多少 题目分析: ...
- Matrix Power Series(矩阵快速幂,二分法/构造矩阵套矩阵)
传送门 题意: 思路: 方法一:(二分) 分析: 奇偶分类讨论: 1.代码: #include <iostream> #include <cstdio> #include &l ...
- I-Matrix Power Series POJ - 3233 矩阵快速幂+分治
I-Matrix Power Series POJ - 3233 矩阵快速幂+分治 Problem Description Given a n × n matrix A and a positive ...
- UVA 10870 Recurrences 矩阵快速幂
题意:给你递推公式要求矩阵快速幂的第n项 解题思路: 经典的矩阵快速幂递推第n项目做法. 解题代码: 1 // File Name: temp.cpp 2 // Author: darkdream 3 ...
- 【矩阵快速幂】233 Matrix HDU - 5015
Think: 1知识点:矩阵快速幂 2题意:定义一个233矩阵,第一行形如233, 23333, 233333, 233-(即a[0][1] = 233, a[0][2] = 2333,a[0][3] ...
- 【做题】SRM701 Div1 Hard - FibonacciStringSum——数学和式&矩阵快速幂
原文链接 https://www.cnblogs.com/cly-none/p/SRM701Div1C.html 题意:定义"Fibonacci string"为没有连续1的01串 ...
- H - Fibonacci POJ - 3070 (矩阵快速幂)
H - Fibonacci POJ - 3070 (矩阵快速幂) Description In the Fibonacci integer sequence, F0 = 0, F1 = 1, and ...
最新文章
- DesignPattern(四)结构型模式(下)
- SortedDictionaryTKey, TValue 类 表示根据键进行排序的键/值对的集合。
- 如何解决win8系统下卸载软件出现错误代码为2502和2503的问题
- 探索 YOLO v3 实现细节 - 第6篇 预测 (完结)
- LUA表与函数的深入理解
- JAVA面向对象程序设计(第二版) 袁绍欣 第六章答案
- 226. Invert Binary Tree
- MySQL索引的一些问题
- 【机器学习】李宏毅 - 02 回归
- AB 罗克韦尔PLC(1769-ERC)相关与 CK-UR12-E01 关于 EtherNetIP 通信 PLC 配置指南
- python怎样遍历列表中数字_Python基础-列表操作(2):列表的遍历和数字列表
- 计算机机房线槽,镀锌线槽机房布线的使用方式
- 研发平台演进:从工具链到开放平台
- java图片识别(文字、字母、数字)
- 为什么宇宙会将最大速度限制在光速
- 5G建设即将收官,杀手级应用还在路上
- android安卓最新版本,安卓X86最新版本发布,目前已支持大部分安卓程序
- html5录音+百度语音实现语音识别
- android版本 6手机版下载,granboard安卓版下载-GranBoard 安卓版v6.1.2-PC6安卓网
- js获取控制浏览器,屏幕,对象宽度和高度,文件大小
热门文章
- 【转】DCT变换的透彻解析
- 业务层勿用继承,不要为了方便舍弃了性能。TʌT不好意思我错了
- 1 用存储过程实现分页,除了上一页,下一页,第一页,和末页外还要有go按钮,以及go到那里的文本框。另外还要在Lable显示“当前x页,一共y页”。注意验证控件的使用和 链接存储过程的内容。...
- BREW做的第一个程序--Hello world!
- 26 JSX深度剖析与使用技巧
- 网络七层协议之物理层
- C++编程经验总结1
- KVM虚拟机共享存储动态迁移与冷迁移
- BZOJ 4066 简单题 ——KD-Tree套替罪羊树
- java文件处理之压缩,分割