矩阵k次幂 采用三重循环

#include<iostream>
using namespace std;

int main()
{
    int n,k;
    long long a[][],b[][],c[][];
    while(cin>>n>>k)
    {
        int x=k;

        for(int i=;i<n;i++)
            for(int j=;j<n;j++)
                {
                    cin>>a[i][j];
                    b[i][j]=a[i][j];
                    c[i][j]=;
                }
    k=k-;
      while(k--)
      {

         for(int i=;i<n;i++)
            for(int j=;j<n;j++)
                for(int m=;m<n;m++)

                 {
             c[i][j]=a[i][m]*b[m][j]+c[i][j];
             }
               if(k>=)
       {   

            for(int i=;i<n;i++)
            for(int j=;j<n;j++)
            {
                a[i][j]=c[i][j];
                c[i][j]=;
            }

        }

       }
       if(x==)
       {
         for(int i=;i<n;i++)
                 {
                    for(int j=;j<n-;j++)
                cout<<b[i][j]<<' ';
             cout<<b[i][n-]<<endl;
         }
       }
       else
       {
                  for(int i=;i<n;i++)
                 {
                    for(int j=;j<n-;j++)

                cout<<c[i][j]<<' ';
             cout<<c[i][n-]<<endl;
         }

     }
     cout<<endl;
    }
    return ;
}

 c[i][j]=a[i][m]*b[m][j]+c[i][j];
这一句是关键，掌握了矩阵相乘的计算公式，则可得到