矩阵的算法

(1)矩阵的加、减运算
C(i,j)=a(i,j)+b(i,j) 加法
C(i,j)=a(i,j)-b(i,j) 减法
(2)矩阵相乘
(矩阵A有M*L个元素，矩阵B有L*N个元素，则矩阵C=A*B有M*N个元素)。矩阵C中任一元素 (i=1,2,…,m; j=1,2,…,n)
#define M 2
#define L 4
#define N 3
void mv(int a[M][L], int b[L][N], int c[M][N])
{
  int i, j, k;
  for(i=0; i<M; i++)
    for(j=0; j<N; j++)
    {
      c[i][j]=0;
      for(k=0; k<L; k++)
      c[i][j]+=a[i][k]*b[k][j];
    }
}
void main()
{
   int a[M][L]={{1,2,3,4},{1,1,1,1}};
   int b[L][N]={{1,1,1},{1,2,1},{2,2,1},{2,3,1}}, c[M][N];
   int i, j;
   mv(a,b,c);
   for(i=0; i<M; i++)
   {
     for(j=0; j<N; j++)
        printf("%4d", c[i][j]);
     printf("/n");
   }
}
（3）矩阵传置
例:有二维数组a(5,5)，要对它实现转置，可用下面两种方式：
#define N 3
void ch1(int a[N][N])
{
   int i, j, t;
   for(i=0; i<N; i++)
     for(j=i+1; j<N; j++)
     {
       t=a[i][j];
       a[i][j]=a[j][i];
       a[j][i]=t;
      }
}
void ch2(int a[N][N])
{
    int i, j, t;
    for(i=1; i<N; i++)
       for(j= 0; j<i; j++)
       {
          t=a[i][j];
          a[i][j]=a[j][i];
          a[j][i]=t;
        }
}
void main()
{
   int a[N][N]={{1,2,3},{4,5,6},{7,8,9}}, i, j;
   ch1(a); /*或ch2(a);*/
   for(i=0; i<N; i++)
   {
     for(j=0; j<N; j++)
         printf("%4d", a[i][j]);
     printf("/n");
   }
}
（4）求二维数组中最小元素及其所在的行和列
基本思路同一维数组，可用下面程序段实现（以二维数组a[3][4]为例）：
‘变量max中存放最大值，row,column存放最大值所在行列号
#define N 4
#define M 3
void min(int a[M][N])
{
    int min, row, column, i, j;
     min=a[0][0];
      row=0;
      column=0;
      for(i=0; i<M; i++)
     for(j=0; j<N; j++)
if(a[i][j]<min)
{ min=a[i][j];
row=i;
column=j;
}
printf("Min=%d/nAt Row%d,Column%d/n", min, row, column);
}
main()
{ int a[M][N]={{1,23,45,-5},{5,6,-7,6},{0,33,8,15}};
min(a);
}