矩阵的乘法C++

无聊中自己写了个矩阵乘法。
已知A，B两个矩阵相乘为C矩阵。



代码如下

#include<iostream>
using namespace std;
#define ROW 100
#define COL 100
int main()
{
	int a[ROW][COL],b[ROW][COL],c[ROW][COL];
	int i,j,k,m,n,q,p;
	cin>>m>>n;//输入第一个矩阵行数列数
	cin>>p>>q;//输入第二个矩阵行数列数
	for(i = 0; i < m; ++i)
		for(j = 0; j < n; ++j)
		cin>>a[i][j];//输入第一个矩阵
	for(i = 0; i < p; ++i)
		for(j = 0; j < q; ++j)
		cin>>b[i][j];//输入第二个矩阵
	for(i = 0; i < m; ++i)//矩阵相乘
	{
		for(j = 0; j < q; ++j)
		{
			c[i][j]=0;
			for(k = 0; k < n; ++k)
			c[i][j] += a[i][k] * b[k][j]; 
		}	
	}
	for(i = 0; i < m; ++i)//输出相乘后的矩阵
	{
		for(j = 0; j < q; ++j)
		cout<<c[i][j]<<" ";
		cout<<endl;
	}
	return 0;
 } 

运行结果如下