矩阵乘法（简单基础）

 矩阵乘法 
问题描述
给定一个N阶矩阵A，输出A的M次幂（M是非负整数）
例如：
　　A =
　　1 2
　　3 4
　　A的2次幂
　　7 10
　　15 22

输入格式
第一行是一个正整数N、M（1<=N<=30, 0<=M<=5），表示矩阵A的阶数和要求的幂数
接下来N行，每行N个绝对值不超过10的非负整数，描述矩阵A的值

输出格式
输出共N行，每行N个整数，表示A的M次幂所对应的矩阵。相邻的数之间用一个空格隔开

样例输入
2 2
1 2
3 4

样例输出
7 10

15 22

详见代码。

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

int n;

void fun(int b[35][35],int a[35][35])
{
    int c[35][35];
    memset(c,0,sizeof(c));
    for (int i=1; i<=n; i++) //行
    {
        for (int j=1; j<=n; j++) //列
        {
            for (int k=1; k<=n; k++) //矩阵中的每一个数
            {
                c[i][j]+=a[i][k]*b[k][j];
            }
        }
    }
    memcpy(b,c,sizeof(c));
}

int main()
{
    int a[35][35],b[35][35];
    int m;
    while (~scanf("%d%d",&n,&m))
    {
        memset(b,0,sizeof(b));
        for (int i=1; i<=n; i++)
        {
            for (int j=1; j<=n; j++)
                scanf("%d",&a[i][j]);
        }
        memcpy(b,a,sizeof(a));
        for (int i=1; i<m; i++)
        {
            fun(b,a);//b乘以a，把结果存在b内
        }
        for (int i=1; i<=n; i++)
        {
            for (int j=1; j<n; j++)
                printf ("%d ",b[i][j]);
            printf ("%d",b[i][n]);
            printf ("\n");
        }
    }
    return 0;
}