矩阵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];
这一句是关键，掌握了矩阵相乘的计算公式，则可得到