矩阵相乘的一维数组实现
欢迎访问我的新博客: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'); } }
这个程序的主要部分还是在矩阵相乘部分的处理吧。
(全文完)