hdu1575之矩阵快速幂入门
来源:互联网 发布:java linux 时间戳 编辑:程序博客网 时间:2024/06/02 21:07
Tr A
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1887 Accepted Submission(s): 1395
Problem Description
A为一个方阵,则Tr A表示A的迹(就是主对角线上各项的和),现要求Tr(A^k)%9973。
Input
数据的第一行是一个T,表示有T组数据。
每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行,每行有n个数据,每个数据的范围是[0,9],表示方阵A的内容。
每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行,每行有n个数据,每个数据的范围是[0,9],表示方阵A的内容。
Output
对应每组数据,输出Tr(A^k)%9973。
Sample Input
22 21 00 13 999999991 2 34 5 67 8 9
Sample Output
22686
#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<string>#include<queue>#include<algorithm>#include<map>#include<iomanip>#define INF 99999999using namespace std;const int MAX=10;const int mod=9973;int s[MAX][MAX],sum[MAX][MAX];int n;void MatrixMult(int a[MAX][MAX],int b[MAX][MAX]){int c[MAX][MAX]={0};for(int i=0;i<n;++i){for(int j=0;j<n;++j){for(int k=0;k<n;++k){c[i][j]+=a[i][k]*b[k][j];}}}for(int i=0;i<n;++i){for(int j=0;j<n;++j)a[i][j]=c[i][j]%mod;}}int Matrix(int k){memset(sum,0,sizeof sum);for(int i=0;i<n;++i)sum[i][i]=1;//初始化矩阵使sum*a=a while(k){if(k&1)MatrixMult(sum,s);MatrixMult(s,s);k>>=1;}int ans=0;for(int i=0;i<n;++i)ans=(ans+sum[i][i])%mod;return ans;}int main(){int t,k;cin>>t;while(t--){cin>>n>>k;for(int i=0;i<n;++i){for(int j=0;j<n;++j)cin>>s[i][j];}cout<<Matrix(k)<<endl;}return 0;}
- hdu1575之矩阵快速幂入门
- HDU1575 矩阵快速幂入门
- hdu1575(矩阵快速幂入门题)
- 快速矩阵幂HDU1575
- 矩阵快速幂 HDU1575
- HDU1575[矩阵快速幂]模板
- hdu1575 TrA 矩阵快速幂
- hdu1575 Tr A (矩阵快速幂)
- hdu1575(矩阵乘法快速幂)
- Hdu1575 - Tr A - 矩阵快速幂
- HDU1575-Tr A(矩阵快速幂)
- hdu1575 Tr A(矩阵快速幂)
- HDU1575 Tr A 【矩阵快速幂】
- HDU1575 Tr A(矩阵快速幂)
- hdu1575 Tr A 矩阵快速幂
- hdu1575 Tr A(矩阵快速幂)
- hdu1575 Tr A (矩阵快速幂)
- [矩阵快速幂]hdu1575 Tr A
- LINUX内核经典面试题
- OpenVPN莫名其妙断线的问题及其解决-confirm
- HDOJ 1824 - Let's go home 简单构图2-set
- 吃牛肉有哪些好处?
- 敌兵布阵
- hdu1575之矩阵快速幂入门
- 软件工程之开发模型及其选择
- mmc线性0-1规划问题
- 畅通工程-并查集-hdu 1232
- hdu4099(斐波纳契数列+字典树)
- Linux系统如何平滑生效NAT-BUGFIX
- 在世界坐标系中放置物体SetTransform函数
- quagga 中 命令框架
- 如何学习IOS