strassen矩阵相乘算法C++代码.docxVIP

Strassen 矩阵相乘算法代码 如果喜欢，请下载使用，谢谢。 #include tchar.h #include iostream #include ctime #include Windows.h using namespace std; templatetypename T class Strassen_class { public: void ADD(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize); void SUB(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize); void MUL(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize);//朴素算法实现 void FillMatrix(T** MatrixA, T** MatrixB, int length);//A,B矩阵赋值 void PrintMatrix(T **MatrixA, int MatrixSize);//打印矩阵 void Strassen(int N, T **MatrixA, T **MatrixB, T **MatrixC);//Strassen算法实现 }; templatetypename T void Strassen_classT::ADD(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) { for (int i = 0; i MatrixSize; i++) { for (int j = 0; j MatrixSize; j++) { MatrixResult[i][j] = MatrixA[i][j] + MatrixB[i][j]; } } } templatetypename T void Strassen_classT::SUB(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) { for (int i = 0; i MatrixSize; i++) { for (int j = 0; j MatrixSize; j++) { MatrixResult[i][j] = MatrixA[i][j] - MatrixB[i][j]; } } } templatetypename T void Strassen_classT::MUL(T** MatrixA, T** MatrixB, T** MatrixResult, int MatrixSize) { for (int i = 0; iMatrixSize; i++) { for (int j = 0; jMatrixSize; j++) { MatrixResult[i][j] = 0; for (int k = 0; kMatrixSize; k++) { MatrixResult[i][j] = MatrixResult[i][j] + MatrixA[i][k] * MatrixB[k][j]; } } } } /* c++使用二维数组，申请动态内存方法 申请 int **A; A = new int *[desired_array_row]; for ( int i = 0; i desired_array_row; i++) A[i] = new int [desired_column_size]; 释放 for ( int i = 0; i your_array_row; i++) delete [] A[i]; delete[] A; */ templatetypename T void Strassen_classT::Strassen(int N, T **MatrixA, T **MatrixB, T **MatrixC) { int HalfSize = N / 2; int newSize = N / 2; if (N = 64) //分治门槛，小于这个值时不再进行递归计算，而是采用常规矩阵计算方法 { MUL(MatrixA, MatrixB, MatrixC, N); } else { T** A11; T** A12; T** A21; T** A22; T** B11; T** B12; T** B21; T** B2