矩阵乘法

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<ctype.h>
#define N 100
int a[N][N],b[N][N];//a,b矩阵
int c[N][N];//a,b的乘积
int main()
{
 int p,q,r,s;//a,b的行列
 int i,j,k;
 
 scanf("%d%d",&p,&q);
 for(i=1;i<=p;i++)
 for(j=1;j<=q;j++)
 scanf("%d",&a[i][j]);
 
 scanf("%d%d",&r,&s);
 for(i=1;i<=r;i++)
 for(j=1;j<=s;j++)
 scanf("%d",&b[i][j]);
 
 for(i=1;i<=p;i++)
 {
  for(j=1;j<=s;j++)
  {
   for(k=1;k<=q;k++)
   c[i][j]+=a[i][k]*b[k][j];//a矩阵的i行与b矩阵在j列对应相乘再相加
   if(j!=1) printf(" ");
   printf("%d",c[i][j]);
  }
  puts("");
 }
 return 0;
}