NYOJ 623 A*B Problem II

A*B Problem II

时间限制：1000 ms  |  内存限制：65535 KB

描述
ACM的C++同学有好多作业要做，最头痛莫过于线性代数了，因为每次做到矩阵相乘的时候，大量的乘法都会把他搞乱，所以他想请你写个程序帮他检验一下计算结果是否正确。

输入有多组测试数据，每行给出一组m,n,k（0<m,n,k<=50）。m,n,k表示两个矩阵的大小,其中：

矩阵A：m行n列。

矩阵B：n行k列。

接下来给出m*n个数表示矩阵A和n*k个数表示矩阵B，对于每个数s，0<=s<=1000。

当m，n，k同时为0时结束。
输出
计算两个矩阵的乘积并输出。

样例输入

2 1 3
1
2
1 2 3
2 2 3
1 2
3 4
1 0 1
0 1 0
0 0 0

样例输出

1 2 3
2 4 6
1 2 1
3 4 3

 思路：当矩阵A的行与矩阵的B的列相等时，两个矩阵才能相乘，其结果为矩阵A的行矩阵B的列

       一个m×n的矩阵a(m,n)左乘一个n×p的矩阵b(n,p)，会得到一个m×p的矩阵c(m,p)。

       矩阵在计算机中也就相当于二维数组，计算算法很简单，就是3个for循环。

#include<stdio.h>
#include<string.h>
#include<algorithm>

using namespace std;

int main()
{
int n,m,k;
while(~scanf("%d%d%d",&m,&n,&k))
{
if(m==0&&n==0&&k==0)
break;
int i,j,t;
int a[51][51] = {0},b[51][51] = {0},c[51][51] = {0};
for(i = 1; i <= m; i ++)
for(j = 1; j <= n; j ++)
scanf("%d",&a[i][j]);
for(i = 1; i <= n; i++)
for(j = 1; j <= k; j++)
scanf("%d",&b[i][j]);
for(i = 1; i <= m; i++)
for(j = 1; j <= k ; j++)
for(t = 1; t <= n; t++)
c[i][j] += a[i][t]*b[t][j];  //c数组用来储存新的来的矩阵

for(i = 1; i <= m ; i ++)
{
for(j = 1; j <= k ; j ++)
printf("%d ",c[i][j]);
printf("\n");
}

}  return 0;
}