算法导论-矩阵相乘(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