矩阵3-n阶矩阵相乘

步骤：

1、两个矩阵的乘法

设有矩阵A33 * B33 = C33，根据矩阵乘法定律，有：

• c00 = a00*b00 + a01b10 + a02*b20;
• c01 = a00*b01 + a01b11 + a02*b21;
• ...

 从而可以得到两个矩阵相乘的算法。其中m1、m2表示乘数，m是积。

void MatrixMutiply(int n, int **m1, int **m2, int **m) {
    //n阶矩阵相乘
    int tmp;

    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            tmp = 0;
            for (int k = 0; k < n; ++k) {
                tmp += m1[i][k] * m2[k][j];
            }
            m[i][j] = tmp;
        }
    }
}

2、求n次幂

有了两个矩阵的运算方法，就可以求n次幂了。两个矩阵呢相乘的结果保存在res，待一次运算结束后，再把res赋给下一次要相乘的矩阵之一，直到得出x次。

void PowMatrix(int x, int n, int **m, int **res) {
    //求n阶矩阵的x次幂，结果用res接收
    int **p1 = m, **p2 = m;
    while (x > 0) {
        MatrixMutiply(n, p1, p2, res);
        **p2 = **res;
        x--;
    }

}

3、调用

至此算法的主要逻辑已经完成，调用时主要需要注意二阶数组的定义和传参问题。 

int **CreateMatrix(int n) {
    // 创建n阶矩阵
    int **matrix = (int **) malloc(sizeof(int) * n * n);
    for (int i = 0; i < n; ++i) {
        matrix[i] = (int *) malloc(sizeof(int) * n);
    }

    return matrix;
}

void PrintMatrix(int **m, int n) {
    // 打印n阶矩阵
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            printf("%d ", m[i][j]);
        }
        printf("\n");
    }
    printf("\n");
}

int main() {
    int n = 3;

    //创建初始矩阵
    int **p = CreateMatrix(n);
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            p[i][j] = i + j + 1;
        }
    }

    //打印初始矩阵
    PrintMatrix(p, n);

    //定义接收结果的矩阵
    int **res = CreateMatrix(n);

    //求x次幂
    int x = 1;
    PowMatrix(x, n, p, res);

    //打印结果
    PrintMatrix(res, n);

    free(p);
    free(res);
    return 0;
}