C语言-矩阵乘法

代码

#include<stdio.h>
#include<string.h>

int a[100][100];//初始矩阵
int c[100][100]= {0}; //上一次乘积
int d[100][100]= {0}; //每次结果
int n;
int b;
int m;

int main()
{
    scanf ("%d %d %d",&n,&b,&m);
    void nyuu(int n);
    nyuu(n);//录入矩阵
    int i,j,p;
    for (i=0; i<n; i++)
    {
        for (j=0; j<n; j++)
        {
            c[i][j]=a[i][j];
        }
    }
    if (b==0)
    {
        for (i=0; i<n; i++)
            for (j=0; j<n; j++)
            {
                if (i==j)
                    d[i][j]=1;
            }
    }//0次幂
    if (b==1)
    {
        for (i=0; i<n; i++)
        {
            int j;
            for (j=0; j<n; j++)
            {
                printf ("%d ",a[i][j]%m);
            }
            printf("\n");
        }
        system("pause");
        return 0;
    }//1次幂
    if (b>1)
    {
        for(p=0; p<b-1; p++) //进行b-1次乘法计算
        {
            for (i=0; i<n; i++) //把c清零
                for (j=0; j<n; j++)
                    d[i][j]=0;
            for (i=0; i<n; i++) //一次乘法
            {
                for (j=0; j<n; j++) //将矩阵乘法的结果的同一行的每个元素分别录入到c矩阵的对应位置
                {
                    int k;
                    for (k=0; k<n; k++)
                        d[i][j]+=c[i][k]*a[k][j];
                }
            }
            for (i=0; i<n; i++)
            {
                for (j=0; j<n; j++)
                {
                    c[i][j]=d[i][j];
                }
            }//到此完成一次乘法
        }
    }
    for (i=0; i<n; i++)
    {
        int j;
        for (j=0; j<n; j++)
        {
            printf ("%d ",d[i][j]%m);
        }
        printf("\n");
    }
    system("pause");
    return 0;
}

void nyuu(int n)//录入矩阵的函数
{
    int i;
    for (i=0; i<n; i++)
    {
        int j;
        for (j=0; j<n; j++)
        {
            scanf ("%d",&a[i][j]);
        }
    }
}