矩阵乘法 （C++代码）

参考代码:

#include<iostream>
using namespace std;

#define MX 50
int ia[MX][MX] = {};
int ib[MX][MX] = {};
int ab[MX][MX] = {};
size_t m, n;

void matrix_mul(int arr[][MX],int arr1[][MX],int arr2[][MX])
{
	for (size_t i = 0; i != m; ++i) {
		for (size_t j = 0; j != m; ++j) {
			for (size_t k = 0; k != m; ++k)
				arr[i][j] += arr1[i][k] * arr2[k][j];
		}
	}
}
void matrix_mi(int n)
{
	if (n == 0) {
		for (size_t i = 0; i != m; ++i) {
			for (size_t j = 0; j != m; ++j) {
				ab[i][j] = ib[i][j];
			}
		}
	}
	else {
		while (n != 0) {
			for (size_t i = 0; i != m; ++i) {
				for (size_t j = 0; j != m; ++j) {
					ab[i][j] = 0;
				}
			}
			matrix_mul(ab, ia, ib);
			for (size_t i = 0; i != m; ++i) {
				for (size_t j = 0; j != m; ++j) {
					ib[i][j] = ab[i][j];
				}
			}
			n--;
		}
	}
}
int main()
{
	cin >> m >> n;
	for (size_t i = 0; i != m; ++i) {
		for (size_t j = 0; j != m; ++j) {
			cin >> ia[i][j];
			if (i == j) {
				ib[i][j] = 1;
			}
		}
	}
	matrix_mi(n);
	for (size_t i = 0; i != m; ++i) {
		for (size_t j = 0; j != m; ++j)
			cout << ab[i][j] << " ";
		cout << endl;
	}
	return 0;
}