實現兩個矩陣相乘的C語言程式

程式功能：實現兩個矩陣相乘的C語言程式，並將其輸出

代碼如下：

 1 #include "stdafx.h"
 2 #include "windows.h"
 3 void Multi(int * left, int * right, int * result, int f1, int f2, int s1, int s2);
 4 
 5 int main()
 6 {
 7     int i, j;
 8     int a[4][3] = 
 9     {   {1,2,3},
10         {4,5,6},
11         {7,8,9},
12         {1,2,3}};
13     int b[3][3] =
14     { { 1,2,3 },
15     { 4,5,6 },
16     {7,8,9} };
17     //c矩陣為a矩陣與b矩陣的乘積
18     int c[4][3];
19     //計算兩個矩陣相乘
20     Multi(a[0], b[0],c[0], 4, 3, 3, 3);
21     //輸出c矩陣
22     for (i = 0; i < 4; i++)
23     {
24         for(j=0;j<3;j++)
25             printf("%4d  ", c[i][j]);
26         printf("\n");
27     }
28     system("pause");
29     return 0;
30 }
31 
32 //矩陣相乘演算法，最後四個參數是兩個相乘的矩陣的行數和列數
33 void Multi(int * left, int * right, int * result,int f1, int f2, int s1, int s2)
34 {
35     //int a[f1][f2], b[s1][s2], c[f1][s2];
36     //a[i][j]==a[i*f2+j]
37     //由矩陣相乘，要求f2=s1，以下用f2
38     for (int i = 0; i < f1; i++)    //i表示第i行
39     {
40         for (int j = 0; j < s2; j++) //j表示第j列
41         {
42             result[i*f2 + j] = 0;        //在這裡 result[i][j] = result[i*f2+j];
43             for(int p=0;p<f2;p++)
44             {
45                 result[i*f2 + j] +=left[i*f2+p]*right[p*f2+j];
46             }
47         }
48     }
49 }

 運行結果：