话不多说,上题:
如果A是一个x行y列的矩阵,B是一个y行z列的矩阵,把A和B相乘,其结果就是另一个x行z列的矩阵C。C矩阵中的每个元素就是由每个元素对应的行数所对应的A矩阵的行的每个元素乘以B矩阵的对于的列数的每个元素之积 , 然后求和(A与B的行 相同);
例:
(1)、第一种对于一个int * 类型的指针对一维数组操作(对于cols 和 rank这是define的一个常量):
void matrix_multiply(int* m1, int* m2, int* r, int x, int y, int z)
{
int k, mid = 0;
for (i = 0; i < x; i++)
{
for (j = 0; j < z; j++)
{
for (k = 0; k < y; k++)
mid += *(m1 + k) * *(m2 + j + rank * k);
*(r + j) = mid;
mid = 0;
}
m1 += cols;
r += rank;
}
}
解释:
对于前两个循环就是对于C数组的每个元素进行赋值 这里就是 * (r + j)操作;
这里的第二个循环就是求值 mid 就集合 每个元素乘积之和;
m1 指针的改变就是,指向下一行,r也是;
2.第二种的解法就是使用int * []
void matrix_multiply(int(*m1)[cols], int(*m2)[rank], int(*r)[rank])
{
for (i = 0; i < rows; i++)
{
for (j = 0; j < rank; j++)
r[i][j] = m1[i][0] * m2[0][j] + m1[i][1] * m2[1][j];
}
}
3.然后就是对于主函数的:
#include <stdio.h>
#define rows 3
#define cols 2
#define rank 4
int i, j;
int main()
{
int a[rows][cols] ;
int b[cols][rank] ;
int c[rows][rank] ;
int x, y, z;
void matrix_multiply(int* m1, int* m2, int* r, int x, int y, int z);
scanf_s("%d%d%d", &x, &y, &z);
for (i = 0; i < x; i++)
{
for (j = 0; j < y; j++)
scanf_s("%d", &a[i][j]);
}
for (i = 0; i < y; i++)
{
for (j = 0; j < z; j++)
scanf_s("%d", &b[i][j]);
}
matrix_multiply(a, b, c, x, y, z);
for (i = 0; i < x; i++)
{
for (j = 0; j < z; j++)
printf("%d ", c[i][j]);
printf("\n");
}
return 0;
}
这个是对应第一种解法的,对于第二种要进行调节,该题来自于《C与指针》;