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;
}