有点难理解。。。
参考代码:
#include <iostream> #include <cstdio> using namespace std; const int N=2; int m; struct Matrix{ int m[N][N]; }matrix; //定义表示矩阵的结构体 Matrix I={1,0, 0,1}; //单位矩阵 Matrix p; inline Matrix multiply(Matrix a,Matrix b){ //inline表示内联函数 实现矩阵乘法 int i,j,k; Matrix re; for(i=0;i<N;i++) for(j=0;j<N;j++){ re.m[i][j]=0; for(k=0;k<N;k++) re.m[i][j]+=a.m[i][k]*b.m[k][j]; re.m[i][j]%=m; } return re; } inline Matrix quick_pow(int n){ Matrix re=p,b=I; while(n>0){ if(n&1){ //当是奇数的时候 b=multiply(b,re); } n=n>>1; //向右移一位 re=multiply(re,re); } return b; } int main() { int b; while(cin>>b>>m){ for(int i=0;i<2;i++) for(int j=0;j<2;j++) cin>>p.m[i][j]; matrix=quick_pow(b); for(int i=0;i<2;i++){ for(int j=0;j<2;j++) cout<<matrix.m[i][j]%m<<" "; //坑 如果m=1且输出单位矩阵就要mod成0 cout<<endl; } } return 0; }
0.0分
1 人评分
C语言程序设计教程(第三版)课后习题5.4 (C语言代码)浏览:1282 |
求组合数 (C语言代码)浏览:1152 |
【排队买票】 (C语言代码)浏览:898 |
WU-蓝桥杯算法提高VIP-Quadratic Equation (C++代码)浏览:1741 |
wu-淘淘的名单 (C++代码)浏览:1329 |
哥德巴赫曾猜测 (C语言代码)浏览:2313 |
C语言程序设计教程(第三版)课后习题6.3 (C语言代码)from DQM浏览:652 |
1014题解浏览:484 |
有关字符,字符串的输入输出函数说明浏览:477 |
C二级辅导-等差数列 (C语言代码)浏览:808 |