1、直接嵌套循环法

最常用的计算矩阵乘积的方法，使用三个嵌套的for循环。

#include <stdio.h>

void multiplyMatrices(int rows1, int cols1, int cols2, int mat1[rows1][cols1], int mat2[cols1][cols2], int result[rows1][cols2]) {
    for (int i = 0; i < rows1; i++) {
        for (int j = 0; j < cols2; j++) {
            result[i][j] = 0;
            for (int k = 0; k < cols1; k++) {
                result[i][j] += mat1[i][k] * mat2[k][j];
            }
        }
    }
}

int main() {
    int mat1[2][2] = {{1, 2}, {3, 4}};
    int mat2[2][2] = {{5, 6}, {7, 8}};
    int result[2][2];

    multiplyMatrices(2, 2, 2, mat1, mat2, result);

    for (int i = 0; i < 2; i++) {
        for (int j = 0; j < 2; j++) {
            printf("%d ", result[i][j]);
        }
        printf("\n");
    }

    return 0;
}

2、通过指针进行矩阵乘法

使用指针来遍历矩阵元素进行乘法操作。

#include <stdio.h>

void multiplyMatrices(int rows1, int cols1, int cols2, int *mat1, int *mat2, int *result) {
    for (int i = 0; i < rows1; i++) {
        for (int j = 0; j < cols2; j++) {
            *(result + i * cols2 + j) = 0;
            for (int k = 0; k < cols1; k++) {
                *(result + i * cols2 + j) += *(mat1 + i * cols1 + k) * *(mat2 + k * cols2 + j);
            }
        }
    }
}

int main() {
    int mat1[2][2] = {{1, 2}, {3, 4}};
    int mat2[2][2] = {{5, 6}, {7, 8}};
    int result[2][2];

    multiplyMatrices(2, 2, 2, (int *)mat1, (int *)mat2, (int *)result);

    for (int i = 0; i < 2; i++) {
        for (int j = 0; j < 2; j++) {
            printf("%d ", result[i][j]);
        }
        printf("\n");
    }

    return 0;
}

3、递归法

虽然递归法不是矩阵乘法的常用方法，但可以使用递归将矩阵乘法分解成较小的子问题来解决。

#include <stdio.h>

void multiplyMatricesRecursive(int rows1, int cols1, int cols2, int mat1[rows1][cols1], int mat2[cols1][cols2], int result[rows1][cols2], int i, int j, int k) {
    if (i >= rows1) return;

    if (j >= cols2) {
        multiplyMatricesRecursive(rows1, cols1, cols2, mat1, mat2, result, i + 1, 0, 0);
        return;
    }

    if (k < cols1) {
        result[i][j] += mat1[i][k] * mat2[k][j];
        multiplyMatricesRecursive(rows1, cols1, cols2, mat1, mat2, result, i, j, k + 1);
    } else {
        multiplyMatricesRecursive(rows1, cols1, cols2, mat1, mat2, result, i, j + 1, 0);
    }
}

int main() {
    int rows1 = 2, cols1 = 2, cols2 = 2;
    int mat1[2][2] = {{1, 2}, {3, 4}};
    int mat2[2][2] = {{5, 6}, {7, 8}};
    int result[2][2] = {0};

    multiplyMatricesRecursive(rows1, cols1, cols2, mat1, mat2, result, 0, 0, 0);

    for (int i = 0; i < rows1; i++) {
        for (int j = 0; j < cols2; j++) {
            printf("%d ", result[i][j]);
        }
        printf("\n");
    }

    return 0;
}