矩阵相乘的一维数组实现

106 阅读 0 评论 70 点赞

我是靠谱客的博主朴素苗条，这篇文章主要介绍矩阵相乘的一维数组实现，现在分享给大家，希望可以做个参考。

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
# 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');
	}
}