欢迎使用CSDN-markdown编辑器

算法导论-矩阵相乘（SQUARE-MATRIX-MULTIPLY） c++

算法小白慢慢入门中~

这段代码主要问题是对二级指针和二维数组的理解上的问题，具体参考二维数组和二级指针
二级指针是指向指针的指针，二维数组其实依然是一维数组，只是每个数组成员是一个数组，所以最终组合成二维数组，故用指针指向二维数组时，不能使用二级指针，正确做法如下：

int a[2][2];
int (*p)[2] = a;

这个算法的核心思想是将n*n的矩阵转化为4个小矩阵，然后递归成1*1求乘积，但这仅适用于n是2的m次幂.
这其实类似于快排策略，不过将一维变为二维，将排序变为乘积。做法都是要对半划分，然后子集递归操作。
具体代码如下，可能写法不太好，希望各位大神指出，蟹蟹

#include<iostream>  
using namespace std;


int n;


void caculate(int **aa, int **bb, int **cc, int x)
{
    int(*a)[2] = (int(*)[2]) aa;
    int(*b)[2] = (int(*)[2])bb;
    int(*c)[2] = (int(*)[2])cc;
    if (x == 1)
    {
        c[0][0] += a[0][0] * b[0][0];
        return;
    }
    int _x = x / 2;
    caculate((int **)a, (int **)b, (int **)c, _x);
    caculate((int **)a, (int **)b + _x, (int **)c + _x, _x);


    caculate((int **)a + _x, (int **)b + n * _x, (int **)c, _x);
    caculate((int **)a + _x, (int **)b + n * _x + _x, (int **)c + _x, _x);


    caculate((int **)a + n * _x, (int **)b, (int **)c + n * _x, _x);
    caculate((int **)a + n * _x, (int **)b + _x, (int **)c + n * _x + _x, _x);


    caculate((int **)a + n * _x + _x, (int **)b + n * _x, (int **)c + n * _x, _x);
    caculate((int **)a + n * _x + _x, (int **)b + _x + n * _x, (int **)c + n * _x + _x, _x);
}


int main()
{
    n = 2;
    int a[2][2] = { { 1, 2 }, { 3, 4 } }, b[2][2] = { { 1, 2 }, { 3, 4 } }, c[2][2] = { {0,0},{0,0} };
    caculate((int **)a, (int **)b, (int **)c, 2);
    system("pause");
    return 0;
}eet/stackedit