矩阵（加法&&乘法&&快速幂）

矩阵相乘：

#include<cstdio>int main(){    int m1,m2,n1,n2;    int t;    int a[105][105];    int b[105][105];    int c[105][105];    scanf("%d",&t);    while(t--)    {        scanf("%d%d",&m1,&n1);        for(int i=0;i<m1;i++)            for(int j=0;j<n1;j++)                scanf("%d",&a[i][j]);         scanf("%d%d",&m2,&n2);        for(int i=0;i<m2;i++)            for(int j=0;j<n2;j++)                scanf("%d",&b[i][j]);         if(n1!=m2)//矩阵相乘的条件：第一个矩阵的列等于第二个矩阵的行         printf("NO\n");        else        {            int sum=0;            printf("YES\n");            for(int i=0;i<m1;i++)            {                for(int j=0;j<n2;j++)                {                    for(int k=0;k<n1;k++)                    {                         sum=sum+a[i][k]*b[k][j];                    }                    c[i][j]=sum;                    sum=0;                }            }               for(int i=0;i<m1;i++)                for(int j=0;j<n2;j++)                    printf("%d%c",c[i][j],j==n2-1?'\n':' ');            }       }    return 0; }

矩阵相加：

#include<cstdio>int main(){    int a[55][55];    int b[55][55];    int c[55][55];    int m,n;    scanf("%d%d",&m,&n);    for(int i=0;i<m;i++)        for(int j=0;j<n;j++)            scanf("%d",&a[i][j]);    for(int i=0;i<m;i++)        for(int j=0;j<n;j++)            scanf("%d",&b[i][j]);    for(int i=0;i<m;i++)        for(int j=0;j<n;j++)            c[i][j]=a[i][j]+b[i][j];    for(int i=0;i<m;i++)        for(int j=0;j<n;j++)            printf("%d%c",c[i][j],j==n-1?'\n':' ');     return 0;}

矩阵快速幂：

#include<cstdio>struct Node{    int a[100][100];}m,n;int M,N;Node mtmul(Node p,Node q)//矩阵乘法 {    int sum=0;    for(int i=0;i<M;i++)    {        for(int j=0;j<N;j++)        {            for(int k=0;k<M;k++)            {                sum+=p.a[i][k]*q.a[k][j];            }            n.a[i][j]=sum;            sum=0;        }    }    return n;}Node mtPow(Node A,int n)//参考快速幂算法 {    if(n==1)    return A;    A=mtPow(A,n/2);    if(n%2==0)        return mtmul(A,A);    else         return mtmul(mtmul(A,A),m); }int main(){    scanf("%d%d",&M,&N);        for(int i=0;i<M;i++)            for(int j=0;j<N;j++)            scanf("%d",&m.a[i][j]);        //mtmul(m,m);        //mtPow(m,3);           for(int i=0;i<M;i++)            for(int j=0;j<N;j++)                printf("%d%c",n.a[i][j],j==N-1?'\n':' ');       return 0;}