计算矩阵的乘积

问题：

　　编写程序，计算两矩阵a与b的乘积，并输出结果

分析：

　　根据两矩阵相乘的规则，矩阵a的列数应与矩阵b的行数相同。若a是4行3列矩阵，b是3行2列矩阵，设乘积为c，则c是4行两列矩阵，由于数学上的矩阵恰好和C语言上二维数组相对应，因此本例说明三个二维数组变量，其所分配的存储空间恰好可以存放a、b、c三个矩阵所有的元素值，为简便计，不妨设矩阵元素为整型数值，矩阵a、b的元素值均由键盘输入，为增强直观性，希望所输入的数据以及输出的结果数据形式均与矩阵的书写顺序相对应，，即各矩阵按行列左右分开、上下对齐。

//实现任意两个矩阵的乘积，并输出结果 
#include<stdio.h>
int main(){
    int i,j,k,s;
    int M,P,X,Y;
    printf("Enter your first of matrix a rows and cols："); 
    scanf("%d %d",&M,&P); 
    printf("Enter your second of matrix a rows and cols：");
    scanf("%d %d",&X,&Y); 
    int a[M+1][P+1],b[X+1][Y+1],c[M+1][Y+1];
    printf("Enter the elements of matrix a(%d*%d)：
",M,P);
    for(i=1;i<=M;i++)
        for(j=1;j<=P;j++)
            scanf("%d",&a[i][j]);
    printf("
");
    printf("Enter the elements of matrix b(%d*%d)：
",X,Y);
    for(j=1;j<=X;j++)
        for(k=1;k<=Y;k++){
            scanf("%d",&b[j][k]);
        }
    printf("

");
    printf("The result of a*b is
");
    for(i=1;i<=M;i++){
        for(k=1;k<=Y;k++){
            s=0;
            for(j=1;j<=X;j++)
                s+=a[i][j]*b[j][k];
            c[i][k]=s;
            printf("%6d",s); 
        }
        printf("
");
    }        
    return 0;
}