hdu 1575

最新推荐文章于 2019-03-18 20:48:45 发布

korun

最新推荐文章于 2019-03-18 20:48:45 发布

阅读量644

点赞数

分类专栏：矩阵

本文链接：https://blog.csdn.net/u011560507/article/details/11619159

版权

矩阵专栏收录该内容

2 篇文章 0 订阅

订阅专栏

hdu 1575 Tr A;

简单的矩阵幂。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 15
int n , k;
struct Mat
{
    int val[maxn][maxn];
    void unit()  ///单位矩阵
    {
        for(int i=0; i<maxn; i++)
            val[i][i] = 1;
    }
    void zero()  ///零矩阵
    {
        memset(val,0,sizeof(val));
    }
} x;
int Ar(const Mat &a)
{
    int num = 0;
    for(int i=0; i<n; i++)
    {
        num += (a.val[i][i]);
        num%=9973;
    }
    return num;
}
Mat operator *(const Mat &a,const Mat &b)
{
    Mat tmp;
    tmp.zero();
     for(int k=0; k<n; k++)
     for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
            {
                tmp.val[i][j] += (a.val[i][k] * b.val[k][j]) % 9973;
                tmp.val[i][j]%=9973;  ///这个地方开始没写，竟然错了很多次。。
            }
    return tmp;
}
Mat operator ^(Mat x,int k)
{
        Mat tmp;
        tmp.zero();
        tmp.unit();
        while( k )
        {
            if(k%2==1) tmp = tmp*x;
            k/=2;
            x = x * x;
        }
        return tmp;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&k);
        x.zero();
        for(int i=0; i<n; i++)
            for(int j=0; j<n; j++)
            {
                scanf("%d",&x.val[i][j]);
            }
        Mat y = x^k;
        printf("%d\n",Ar(y));
    }
    return 0;
}