有关矩阵的运算请参考:http://www2.edu-edu.com.cn/lesson_crs78/self/j_0022/soft/ch0605.html
//矩阵乘法A*B
/*
矩阵相乘的条件:前面矩阵的列数等于后面矩阵的行数(对应这里的na=mb)
A:ma*na B:mb*nb
C=A*B: ma*nb
*/
#include<stdio.h>
#include<stdlib.h>
int main()
{
int ma,na,mb,nb,i,j,k;
int t=0;
int **A,**B,**C;
int m,n,p;
printf("(Matrix A)Input ma and na\n");
scanf("%d%d",&ma,&na);
A=(int**)malloc(ma*sizeof(int*));
for(i=0;i<ma;i++)
A[i]=(int*)malloc(na*sizeof(int));//Initial A(为矩阵A分配空间)
/*----------------------------------------------*/
printf("Input A:\n");
for(i=0;i<ma;i++)
for(j=0;j<na;j++)
scanf("%d",&A[i][j]);//Input A
/*----------------------------------------------*/
printf("(Matrix B)Input mb and nb\n");
scanf("%d%d",&mb,&nb);
B=(int**)malloc(mb*sizeof(int*));
for(i=0;i<mb;i++)
B[i]=(int*)malloc(nb*sizeof(int));//Initial B(为矩阵B分配空间)
/*----------------------------------------------*/
printf("Input B:\n");
for(i=0;i<mb;i++)
for(j=0;j<nb;j++)
scanf("%d",&B[i][j]);//Input B
/*----------------------------------------------*/
printf("matrix A:\n");
for(i=0;i<ma;i++)
{
for(j=0;j<na;j++)
printf("%d ",A[i][j]);
printf("\n");
}//Output A
/*----------------------------------------------*/
printf("matrix B:\n");
for(i=0;i<mb;i++)
{
for(j=0;j<nb;j++)
printf("%d ",B[i][j]);
printf("\n");
}//Output B
/*----------------------------------------------*/
/*----------------------------------------------*/
m=ma;
n=nb;
p=mb;//p保存A的列数和B的行数
printf("m:%d n:%d\n",m,n);
C=(int**)malloc(m*sizeof(int*));
for(i=0;i<m;i++)
C[i]=(int*)malloc(n*sizeof(int));//Initial C=A*B(为矩阵C分配空间)
/*----------------------------------------------*/
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
C[i][j]=0;
for(k=0;k<p;k++)
C[i][j]+=A[i][k]*B[k][j];//计算矩阵相乘的要点(循环内累加)
}
}//Calculate C=A*B
/*----------------------------------------------*/
printf("matrix C=A*B:\n");
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
printf("%d ",C[i][j]);
printf("\n");
}//Output C=A*B
return 0;
}