矩阵的快速幂
2016-05-20 14:28
369 查看
原博客地址:http://blog.csdn.net/mark_meego/article/details/8654619讲矩阵快速幂的http://blog.csdn.net/y990041769/article/details/39694583 http://www.cnblogs.com/vongang/archive/2012/04/01/2429015.html题目描述: 给定一个n*n的矩阵,求该矩阵的k次幂,即P^k。
输入:
输入包含多组测试数据。
数据的第一行为一个整数T(0<T<=10),表示要求矩阵的个数。
接下来有T组测试数据,每组数据格式如下:
第一行:两个整数n(2<=n<=10)、k(1<=k<=5),两个数字之间用一个空格隔开,含义如上所示。
接下来有n行,每行n个正整数,其中,第i行第j个整数表示矩阵中第i行第j列的矩阵元素Pij且(0<=Pij<=10)。另外,数据保证最后结果不会超过10^8。
输出:
对于每组测试数据,输出其结果。格式为:
n行n列个整数,每行数之间用空格隔开,注意,每行最后一个数后面不应该有多余的空格。
样例输入:
样例输出:
[/code]
输入:
输入包含多组测试数据。
数据的第一行为一个整数T(0<T<=10),表示要求矩阵的个数。
接下来有T组测试数据,每组数据格式如下:
第一行:两个整数n(2<=n<=10)、k(1<=k<=5),两个数字之间用一个空格隔开,含义如上所示。
接下来有n行,每行n个正整数,其中,第i行第j个整数表示矩阵中第i行第j列的矩阵元素Pij且(0<=Pij<=10)。另外,数据保证最后结果不会超过10^8。
输出:
对于每组测试数据,输出其结果。格式为:
n行n列个整数,每行数之间用空格隔开,注意,每行最后一个数后面不应该有多余的空格。
样例输入:
3 2 2 9 8 9 3 3 3 4 8 4 9 3 0 3 5 7 5 2 4 0 3 0 1 0 0 5 8 5 8 9 8 5 3 9 6 1 7 8 7 2 5 7 3
样例输出:
153 96 108 81 1216 1248 708 1089 927 504 1161 1151 739 47 29 41 22 16 147 103 73 116 94 162 108 153 168 126 163 67 112 158 122 152 93 93 111 97
由于时间比较仓促,忽略题目的一些细节要求。
[code]<pre class="cpp" name="code">#include <stdio.h> #include <cstring> #define T 10 struct matrix { int m[T][T]; }res,origin[T]; void calmatrix(matrix ori[],int q[][2],int matrix_count); matrix multiply(matrix x,matrix y,int n); int main(void) { int matrix_count; int matrix_n_k[T][2]; do { scanf("%d",&matrix_count); if(matrix_count>0) { for(int i=0;i<matrix_count;i++) { scanf("%d %d",&matrix_n_k[i][0],&matrix_n_k[i][1]); for(int j=0;j<matrix_n_k[i][0];j++) for(int k=0;k<matrix_n_k[i][0];k++) scanf("%d",&origin[i].m[j][k]); } calmatrix(origin,matrix_n_k,matrix_count); } }while(matrix_count!=0); printf("finish\n"); return 0; } void calmatrix(matrix ori[],int q[][2],int matrix_count) { int a,b; for(int count=0;count<matrix_count;count++) { int n=q[count][1]; memset(res.m,0,sizeof(res.m)); for( a=0;a<q[count][0];a++) res.m[a][a]=1; //将res.m初始化为单位矩阵 //这里用到二进制思想,如A^156,而156(10)=10011100(2) ,也就有A^156=>(A^4)*(A^8)*(A^16)*(A^128) while(n) { if(n&1) { res=multiply(res,ori[count],q[count][0]); } ori[count]=multiply(ori[count],ori[count],q[count][0]); n>>=1; } for( a=0;a<q[count][0];a++) { for( b=0;b<q[count][0];b++) printf("%d ",res.m[a][b]); printf("\n"); } } } matrix multiply(matrix x,matrix y,int n) { int a,b,c; matrix temp; memset(temp.m,0,sizeof(temp.m)); for(a=0;a<n;a++) for( b=0;b<n;b++) for( c=0;c<n;c++) temp.m[a][b]+=x.m[a][c]*y.m[c][b]; return temp; }
[/code]
相关文章推荐
- 怎么写inf cab makecab
- 一步一步 搞定RSA(公钥、私钥)
- iOS学习之—— UIVisualEffectView
- gnocchi-采样数据存储流程分析(002)--数据的异步统计
- JavaWeb项目使用BlazeDS整合Flex进行RemoteObject数据交互
- Linux下使用yum安装软件
- 实验3:白盒测试
- linux基础命令学习(三)文件搜索 find
- jquery.cookie.js插件使用方法
- 关于接入新浪微博sdk项目崩溃的解决办法
- Android View的加载过程
- 在线网站技术分析工具
- android studio 2.2下载地址
- jsp放在web-inf下的注意事项
- HTML+CSS基础(三):认识html标签
- [JavaScript] +function(){} 作用
- [转]PHP--优化PHP代码建议
- Error running app: Instant Run requires 'Tools | Android | Enable ADB integration' to be enabled.
- Linux下时间戳格式和标准时间格式的转换
- Android替代HashMap的SparseArray