實現兩個N*N矩陣的乘法，矩陣由一維數組表示

實現兩個N*N矩陣的乘法，矩陣由一維數組表示。

先介紹一下矩陣的加法：

      Add( rows,          ( i= ;i<rows;i++         ( j=;j<cols;j++           result[i][j]=mat1[i][j]+      } 

若兩個矩陣要做乘法運：只有在一個矩陣的行數與另一個矩陣的列數相同時，才能做兩個矩陣的乘法。

如何得到矩陣的轉置：

矩陣的轉置也是一個矩陣，原始矩陣中的行轉變為轉置矩陣的列。例如，有下述一個3×3矩陣：

1 2 3
6 7 8
4 5 9

那麼它的轉置矩陣為：

1 6 4
2 7 5
3 8 9

讓我們從程序員的角度仔細地考察一下這一現象。假設原始數組為M，轉置矩陣為MT。那麼M[1][0]＝6，在轉置矩陣中我們發現MT [0][1]＝6。因此，我們能夠得到程序化的結論：轉置一個矩陣實際上就是對換下標變量。用技術術語講：

1. MT[Row][Column] = M[Column][Row]; 

下面是得到轉置矩陣的C代碼：

1. void show_transpose(float mat[][10],int row,int col)    
2. {    
3.    int i,j;    
4.    for(i=0;i<row;i++)    
5.    {    
6.       for(j=0;j<col;j++)    
7.          printf("%f\t",mat[j][i]);    
8.       printf("\n");    
9.    }    
10. }   

這個方法顯示了矩陣的轉置。

1. #include<iostream>  
2. using namespace std;  
3. #define size 2  
4.   
5. int multi(int *a , int *b , int N)  
6. {  
7.     int i , j , k , temp;  
8.     int *c = (int*)malloc(N * N * sizeof(int));  
9.   
10.     for(i = 0 ; i < N ; i++)  
11.     {  
12.         for(j = 0 ; j < N ; j++)  
13.         {  
14.             temp = i * N + j;  
15.             *(c + temp) = 0;  
16.             for(k = 0 ; k < N ; k++)  
17.             {  
18.                 *(c + temp) += a[i * N + k] * b[k * N + j];  
19.             }  
20.             cout<<*(c + temp)<<" ";  
21.         }  
22.     }  
23.     return *c;  
24. }  
25.   
26. int main()  
27. {  
28.     int a[size * size] = {2 , 1 , 4 , 3};  
29.     int b[size * size] = {1 , -1 , 3 , 2};  
30.     multi(a , b , size);  
31.   
32.     return 0;  
33. }