两个N*N矩阵的乘法

实现两个N*N矩阵的乘法，矩阵由一维数组表示。

[cpp] view plaincopy

1. void MatrixMultiplication(int * pArry1, int *pArry2, int * pDestArry, int len)  
2. {  
3.     int row ;  
4.     int col;  
5.     row= col = (int)sqrt((double)len);  
6.     if (row * col != len)  //矩阵不为方阵则提示矩阵无效  
7.     {  
8.         cout << "Matrix Invalidate!" << endl;  
9.     }  
10.     int i, j, k;  
11.     int sum;  
12.     for (i = 0; i < row; i++)  
13.     {  
14.         for (j = 0; j < col; j++)  
15.         {  
16.             sum = 0;  
17.             for (k = 0; k < row; k++)  
18.             {  
19.                 sum += pArry1[row * i + k] * pArry2[row * k + j];  
20.                //某行、某列所有对应元素相乘求和即为新矩阵某行某列的元素  
21.             }  
22.             pDestArry[row * i + j] = sum;  //矩阵相乘所得新矩阵        
23.         }  
24.     }  
25. }