电子科大POJ "任意阶矩阵相乘"

    　　　任意阶矩阵的乘法

  Time Limit: 3000/1000MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others)

C-source:

<span style="color:#333333;">#include<stdio.h>
#include<malloc.h>
int main(void)
{
 int m,g,n;
 int i,j,k;

 int number;
 double **a=NULL,**b=NULL,**result=NULL;
 printf("Please input number:
");
 scanf("%d",&number);

while(number!=0)
{
    number--;

    printf("Please input m,g,n:
");
    scanf("%d%d%d",&m,&g,&n);


    a=(double**)malloc(m*sizeof(double*));
    for(i=0;i<m;i++)
        a[i]=(double*)malloc(g*sizeof(double));

    printf("Please input the first jz:
");
    for(i=0;i<m;i++)
        for(j=0;j<g;j++)
        scanf("%lf",&a[i][j]);

    b=(double**)malloc(g*sizeof(double *));
    for(i=0;i<g;i++)
        b[i]=(double*)malloc(n*sizeof(double));

    printf("Please input the second jz
");
    for(i=0;i<g;i++)
        for(j=0;j<n;j++)
            scanf("%lf",&b[i][j]);

    result=(double**)malloc(m*sizeof(double*));
        for(i=0;i<m;i++)
            result[i]=(double*)malloc(n*sizeof(double));

    for(i=0;i<m;i++)
        for(j=0;j<n;j++)
            result[i][j]=0;

    for(i=0;i<m;i++)
        for(j=0;j<n;j++)
            for(k=0;k<g;k++)
                result[i][j]+=a[i][k]*b[k][j];

    printf("answer:
");
    for(i=0;i<m;i++)
    {
        for(j=0;j<n;j++)
            printf("%5g",result[i][j]);
            printf("
");
    }
}
    for(i=0;i<m;i++)
    {
        free(a[i]);
        a[i]=NULL;
        free(result[i]);
        result[i]=NULL;
    }
    free(a);
    a=NULL;
    free(result);
    result=NULL;


    for(i=0;i<g;i++)
    {   free(b[i]);
        b[i]=NULL;
    }
    free(b);
    b=NULL;

    return 0;
}