矩阵的乘法1055 Problem C: Matrix Problem (III) : Array Practice

问题 C: Matrix Problem (III) : Array Practice

时间限制: 1 Sec  内存限制: 4 MB
提交: 15450  解决: 5169
[提交][状态]

题目描述

求两个矩阵A、B的乘积C=AB。根据矩阵乘法的定义，只有A的列数和B的行数相同才能相乘。可以确保所有运算结果都在int类型的范围之内。

输入

输入数据为多个矩阵（最少2个），每个矩阵以两个正整数m和n开始，满足0<m,n<=100，接下来为一个m行n列的矩阵A。当输入的m和n均为0时，表示输入数据结束。

输出

对输入的矩阵两两相乘：第1个和第2个相乘、第1个和第2个相乘的结果和第3个相乘……按顺序输出矩阵相乘的结果：每行两个元素之间用一个空格分开，每行最后一个元素之后为一个换行，在下一行开始输出矩阵的下一行。

若前k个矩阵相乘的结果和第k+1个矩阵无法相乘（即不满足矩阵乘法定义），则输出“Not satisfied the definition of matrix multiplication!”。然后用第k+1个矩阵去和第k+2个矩阵相乘。最后一个矩阵只做乘数。

每两组输出之间用一个空行分隔开。

样例输入

2 3
1 1 1
1 1 1
3 3
1 2 3
4 5 6
7 8 9
3 1
0
0
0
0 0

样例输出

12 15 18
12 15 18

0
0

提示

矩阵的乘法就是一行乘以一列加起来做一个元素。

Append Code

解法如下：

#include <stdio.h>
#include <string.h>
int main() {
	int i, j, left1, left2, right1, right2;
	int a[101][101], b[101][101], des[101][101] = {0};
	scanf("%d %d", &left1, &right1);
	for (i = 0; i < left1; i++)
		for (j = 0; j < right1; j++)
			scanf("%d", &a[i][j]);
	scanf("%d %d", &left2, &right2);
	for (i = 0; i < left2; i++)
		for (j = 0; j < right2; j++)
			scanf("%d", &b[i][j]);
	if (right1 != left2) {
		printf("Not satisfied the definition of matrix multiplication!\n");
		left1 = left2;
		right1 = right2;
		for (i = 0; i < left1; i++)
			for (j = 0; j < right1; j++) {
				des[i][j] = b[i][j];
			}
	} else {
		for (i = 0; i < left1; i++)
			for (int j = 0; j < right2; j++)
				for (int k = 0; k < right1; k++) {
					des[i][j] += a[i][k] * b[k][j];
				}
		for (int i = 0; i < left1; i++) {
			for (int j = 0; j < right2; j++) {
				if (j != right2 - 1)
					printf("%d ", des[i][j]);
				else
					printf("%d", des[i][j]);
			}
			printf("\n");
		}
		right1 = right2;
	}

	while (1) {
		printf("\n");//两输出间的空行。
		int des2[101][101], c[101][101], k;
		memset(c, 0, sizeof(c));
		scanf("%d %d", &left2, &right2);
		if (left2 == 0 && right2 == 0)
			break;
		for (i = 0; i < left2; i++)
			for (j = 0; j < right2; j++)
				scanf("%d", &c[i][j]);
		if (right1 != left2) {
			printf("Not satisfied the definition of matrix multiplication!\n");
			left1 = left2;
			right1 = right2;
			memset(des, 0, sizeof(des));
			for (i = 0; i < left1; i++)
				for (j = 0; j < right1; j++) {
					des[i][j] = c[i][j];
				}
		} else {
			for (i = 0; i < left1; i++)
				for (int j = 0; j < right2; j++)
					for ( k = 0; k < right1; k++) {
						des2[i][j] += des[i][k] * c[k][j];
					}
			memset(des, 0, sizeof(des));
			for (int i = 0; i < left1; i++) {
				for (int j = 0; j < right2; j++) {
					if (j != right2 - 1)
						printf("%d ", des2[i][j]);
					else
						printf("%d", des2[i][j]);
					des[i][j] = des2[i][j];
				}
				printf("\n");
			}
			right1 = right2;
			memset(des2, 0, sizeof(des2));
		}
	}
	return 0;
}//所有的memset都是为了清零以防止数据干扰。

代码虽烂，但请勿喷。谁不是这样过来的呢。