L1-048 矩阵A乘以B (15 分)
输入格式:
输入先后给出两个矩阵A和B。对于每个矩阵,首先在一行中给出其行数R和列数C,随后R行,每行给出C个整数,以1个空格分隔,且行首尾没有多余的空格。输入保证两个矩阵的R和C都是正数,并且所有整数的绝对值不超过100。
输出格式:
若输入的两个矩阵的规模是匹配的,则按照输入的格式输出乘积矩阵AB,否则输出Error: Ca != Rb,其中Ca是A的列数,Rb是B的行数。
输入样例1:
2 3
1 2 3
4 5 6
3 4
7 8 9 0
-1 -2 -3 -4
5 6 7 8
输出样例1:
2 4
20 22 24 16
53 58 63 28
输入样例2:
3 2
38 26
43 -5
0 17
3 2
-11 57
99 68
81 72
输出样例2:
Error: 2 != 3
这一题对没有学过线性代数的朋友不是很友好(比如我)。我在理解题目上花了很多时间
所以咱们得先解释一下,什么是矩阵相乘!
1.矩阵相乘A*B简单来说,当前一个矩阵A的列数等于后面一个矩阵B的行数时,拿前面矩阵的行元素去乘后面一个矩阵的列元素,将所得的乘积求和。比如题目范例中A矩阵有两行,B矩阵有四列。
答案中的 20 22 24 16怎么来的呢?
即sum = 1 X 7+ 2 X(-1)+ 3 X 5 == 20
后面的元素以此类推,第二行的元素是A矩阵第二行与B矩阵相乘之后求和的结果,我们也可以得到:
新的矩阵行数与A相等,列数与B相等(ra,cb)
2.解决了什么是矩阵相乘,后面就到咱们的代码环节。既然是矩阵相乘,矩阵相乘的结果也是一个矩阵,所以我们需要三个二维数组,还有他们的行列元素
3.输入完后,我们第一步就要判断A的列与B的行是否相等(ca!=rb)因为矩阵相乘这个是前提,有这个前提才有相乘的资格。
4.到了我们的关键的计算环节了,从上面发现的规律我们可以知道,新的矩阵行数与A相等,列数与B相等,两层循环分别代表新矩阵行和列,这里很重要,先留意一下,之后我们会发现,只有这两层循环不够用,因为我们无法遍历整个A矩阵。
所以我们不得不使用第三层循环来代表第一个矩阵的列数,为什么不用第四层循环来代替B的行数呢?因为B的行数等于A的列数,当然我们还得有一个储存求和的值,这个值要放在第二层循环内,让每轮计算都可以从0开始计数。最后在第二层循环内,第三层循环外,将求到的和赋值给C的元素
5.最后的输出也有坑点!!!直接输出C矩阵会导致C矩阵每一行最后一个元素多一个空格符,无法AC!!所以要在输出每一行最后一个元素时,只输出它,不要多加空格符!!而且输出要放在else内!!
#include <iostream>
using namespace std;
int main() {
int a[100][100];
int b[100][100];
int c[100][100];
int i, j, k;
int ra, ca, rb, cb;
cin >> ra >> ca;
for (i = 0; i < ra; i++)
for (j = 0; j < ca; j++)
cin >> a[i][j];
cin >> rb >> cb;
for (i = 0; i < rb; i++)
for (j = 0; j < cb; j++)
cin >> b[i][j];
/*for (i = 0; i < ra; i++)
{
for (j = 0; j < ca; j++)
cout << a[i][j] << " ";
cout << endl;
}
cout << endl;
for (i = 0; i < rb; i++)
{
for (j = 0; j < cb; j++)
cout << b[i][j] << " ";
cout << endl;
}*/
if (ca != rb)//先判断ca和rb
printf("Error: %d != %d", ca, rb);
else
{
for (i = 0; i < ra; i++)
for (j = 0; j < cb; j++)
{
int sum = 0;//注意sum的位置
for (k = 0; k < ca; k++)
{
sum += a[i][k] * b[k][j];
}//题目要求的c[i][j]要放在二层内,三层外
c[i][j] = sum;
}
cout << ra << " " << cb << endl;
for (i = 0; i < ra; i++)//输出放在else内
{
for (j = 0; j < cb; j++)
if (j<cb-1)//到每一行最后一个时,只输出它
cout << c[i][j] << " ";
else
cout<<c[i][j];
cout << endl;
}
}
return 0;
}