矩阵相乘的一维数组实现

欢迎访问我的新博客：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');
	}
}

这个程序的主要部分还是在矩阵相乘部分的处理吧。

（全文完）