矩阵乘法
题目描述
分别给定
n
∗
p
n * p
n∗p 和
p
∗
m
p * m
p∗m 的两个矩阵
A
A
A 和
B
B
B,求
A
∗
B
A * B
A∗B。
1
≤
n
,
p
,
m
≤
500
1\le n,p,m\le 500
1≤n,p,m≤500,
−
1
0
9
≤
A
i
,
j
,
B
i
,
j
≤
1
0
9
-10^9\le A_{i,j}, B_{i,j}\le 10^9
−109≤Ai,j,Bi,j≤109
输入描述
第一行三个正整数
n
,
p
,
m
n,p,m
n,p,m,表示矩阵的长宽。
之后的
n
n
n 行,每行
p
p
p 个整数,表示矩阵
A
A
A。
之后的
p
p
p 行,每行
m
m
m 个整数,表示矩阵
B
B
B。
输出描述
输出 n n n 行,每行 m m m 个整数,表示矩阵 A ∗ B A * B A∗B,每个数模 1 0 9 + 7 10 ^ 9 + 7 109+7 输出。行末尾没有多余的空格或换行。
样例
输入
3 4 5
-2 -8 -9 8
-10 0 6 -8
-10 -6 6 9
4 -7 5 -5 9
10 -2 -10 5 5
-3 -7 -3 8 -2
-6 7 7 3 -2
输出
999999898 149 153 999999929 999999951
999999997 999999979 999999883 74 999999921
999999835 103 55 95 999999857
来源
MG
C语言实现
#include <stdio.h>
#include <string.h>
typedef long long int ll;
const unsigned int M = 1000000007;
ll a[1000][1000],b[1000][1000],c[1000][1000];
ll MM(ll n){
if(n<0)
return MM(n+M);
return n%M;
}
int main()
{
int n,p,m,i,j,k;
scanf("%d %d %d",&n,&p,&m);
memset(c,0,sizeof(c));
for(i=0;i<n;i++){
for(j=0;j<p;j++){
scanf("%lld",&a[i][j]);
}
}
for(i=0;i<p;i++){
for(j=0;j<m;j++){
scanf("%lld",&b[i][j]);
}
}
for(i=0;i<n;i++){
for(j=0;j<m;j++){
for(k=0;k<p;k++){
c[i][j]+=a[i][k]*b[k][j]%M;
c[i][j]=c[i][j]%M;
}
}
}
for(i=0;i<n;i++){
for(j=0;j<m-1;j++){
printf("%lld ",MM(c[i][j]));
}
printf("%lld\n",MM(c[i][m-1]));
}
return 0;
}