算法提高 矩阵相乘
时间限制:1.0s 内存限制:256.0MB
问题描述
小明最近在为线性代数而头疼,线性代数确实很抽象(也很无聊),可惜他的老师正在讲这矩阵乘法这一段内容。
当然,小明上课打瞌睡也没问题,但线性代数的习题可是很可怕的。
小明希望你来帮他完成这个任务。
现在给你一个ai行aj列的矩阵和一个bi行bj列的矩阵,
要你求出他们相乘的积(当然也是矩阵)。
(输入数据保证aj=bi,不需要判断)
当然,小明上课打瞌睡也没问题,但线性代数的习题可是很可怕的。
小明希望你来帮他完成这个任务。
现在给你一个ai行aj列的矩阵和一个bi行bj列的矩阵,
要你求出他们相乘的积(当然也是矩阵)。
(输入数据保证aj=bi,不需要判断)
输入格式
输入文件共有ai+bi+2行,并且输入的所有数为整数(long long范围内)。
第1行:ai 和 aj
第2~ai+2行:矩阵a的所有元素
第ai+3行:bi 和 bj
第ai+3~ai+bi+3行:矩阵b的所有元素
第1行:ai 和 aj
第2~ai+2行:矩阵a的所有元素
第ai+3行:bi 和 bj
第ai+3~ai+bi+3行:矩阵b的所有元素
输出格式
输出矩阵a和矩阵b的积(矩阵c)
(ai行bj列)
(ai行bj列)
样例输入
2 2
12 23
45 56
2 2
78 89
45 56
12 23
45 56
2 2
78 89
45 56
样例输出
1971 2356
6030 7141
6030 7141
分析:很简单的矩阵相乘,线代的基础公式套用一下,注意不要搞错变量位置即可。
代码:
#include<stdio.h>
#include<limits.h>
#include<string.h>
#define MAXLEN 1000
long long a[MAXLEN][MAXLEN];
long long b[MAXLEN][MAXLEN];
long long c[MAXLEN][MAXLEN];
int main()
{
long long m1, n1;
long long m2, n2;
memset(a, 0, sizeof(a));
memset(b, 0, sizeof(b));
memset(c, 0, sizeof(c));
scanf("%lld %lld", &m1, &n1);
for(int i = 0; i < m1; i++)
for(int j = 0; j < n1; j++)
scanf("%lld", &a[i][j]);
scanf("%lld %lld", &m2, &n2);
for(int i = 0; i < m2; i++)
for(int j = 0; j < n2; j++)
scanf("%lld", &b[i][j]);
for(int i = 0; i < m1; i++)
for(int j = 0; j < n2; j++)
for(int k = 0; k < n1; k++)
c[i][j] += a[i][k] * b[k][j];
for(int i = 0; i < m1; i++)
{
for(int j = 0; j < n2; j++)
printf("%d ", c[i][j]);
printf("\n");
}
return 0;
}