欢迎访问我的新博客:http://www.milkcu.com/blog/
原文地址:http://www.milkcu.com/blog/archives/1369108860.html
矩阵相乘,听起来很简单,首先想到的是二维数组实现。那用一维数组如何实现呢?其实多维数组在本质上也是一维数组。
源代码如下,一个小小的程序,却总是因为某个字母而出错。
# include <stdio.h>
# define MAX 1000
int main(void)
{
int a[MAX], b[MAX];
int ma, na;
int mb, nb; //na == mb
//matrix a
scanf("%d%d", &ma, &na);
for(int i = 0; i < ma; i++) {
for(int j = 0; j < na; j++) {
scanf("%d", &a[i * na + j]);
}
}
printf("----- ----- ----- -----\n");
for(int i = 0; i < ma; i++) {
for(int j = 0; j < na; j++) {
printf("%5d", a[i * na + j]);
}
putchar('\n');
}
printf("----- ----- ----- -----\n");
//matrix b
scanf("%d%d", &mb, &nb);
for(int i = 0; i < mb; i++) {
for(int j = 0; j < nb; j++) {
scanf("%d", &b[i * nb + j]);
}
}
printf("----- ----- ----- -----\n");
for(int i = 0; i < mb; i++) {
for(int j = 0; j < nb; j++) {
printf("%5d", b[i * nb + j]);
}
putchar('\n');
}
printf("----- ----- ----- -----\n");
//matrix a*b
for(int i = 0; i < ma; i++) {
for(int j = 0; j < nb; j++) {
int s = 0;
for(int x = 0; x < ma; x++) {
for(int y = 0; y < na; y++) {
//s += a[i][x] * b[y][j];
s += a[i * na + x] * b[y * nb + j];
//printf("+%d+", s);
}
}
printf("%5d", s);
}
putchar('\n');
}
}
这个程序的主要部分还是在矩阵相乘部分的处理吧。
(全文完)