注意事项: 恶补一波线性代数,矩阵的零次幂竟然是单位矩阵,太伤心了┭┮﹏┭┮
#include<iostream> #include<algorithm> #include<cmath> #include<string> #include<cstring> using namespace std; int a[1001][1001]; int b[1001][1001]; int c[1001][1001]; int main() { int n, m; cin >> n >> m; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) cin >> a[i][j]; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) b[i][j] = a[i][j]; if (m == 0) { for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) { if (i == j) a[i][j] = 1; else a[i][j] = 0; } } else { for (int i = 1; i < m; i++) { memset(c, 0, sizeof(c)); for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) { for (int l = 0; l < n; l++) { c[j][k] += a[j][l] * b[l][k]; } } } for (int j = 0; j < n; j++) for (int k = 0; k < n; k++) a[j][k] = c[j][k]; } } for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { cout << a[i][j]; if (j != n) cout << " "; } cout << endl; } return 0; }
0.0分
0 人评分
C语言程序设计教程(第三版)课后习题1.6 (C语言代码)浏览:589 |
三角形 (C++代码)递推浏览:772 |
简单的a+b (C语言代码)浏览:550 |
printf基础练习2 (C语言代码)浏览:515 |
C语言程序设计教程(第三版)课后习题10.1 (C语言代码)浏览:774 |
C语言程序设计教程(第三版)课后习题1.5 (C语言代码)浏览:555 |
C二级辅导-公约公倍 (C语言代码)浏览:1313 |
1202题解浏览:639 |
简单的a+b (C语言代码)浏览:576 |
C二级辅导-阶乘数列 (C语言代码)浏览:663 |