矩阵快速幂hdu 1575

本文主要是介绍矩阵快速幂hdu 1575，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

矩阵快速幂就是求方阵A的n次方，一定得是方阵.复杂度为logn，对角线相加就是矩阵的值
主要是套模板,下面给出我收藏的模板

struct mat{int m[maxn][maxn];
}unit;
//矩阵乘法
mat operator * (mat a,mat b)
{mat ret;ll x;for(int i=0;i<n;i++)for(int j=0;j<n;j++){x=0;for(int k=0;k<n;k++)x+=mod((ll)a.m[i][k]*b.m[k][j]);ret.m[i][j]=mod(x);}return ret;
}
void init_unit()
{for(int i=0;i<maxn;i++)unit.m[i][i]=1;return ;
}
mat pow_mat(mat a,ll n)
{mat ret=unit;while(n){if(n&1){ret= ret*a;}n>>=1;a=a*a;}return ret;
}

#include<iostream>
#include<algorithm>
#include<iomanip>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<string>
typedef long long ll;
using namespace std;
#define mod(x) ((x)%9973)
int n;
int t;
int x;
int k;
const int maxn = 11;
struct mat{int m[maxn][maxn];
}unit;
//矩阵乘法
mat operator * (mat a,mat b)
{mat ret;ll x;for(int i=0;i<n;i++)for(int j=0;j<n;j++){x=0;for(int k=0;k<n;k++)x+=mod((ll)a.m[i][k]*b.m[k][j]);ret.m[i][j]=mod(x);}return ret;
}
void init_unit()
{for(int i=0;i<maxn;i++)unit.m[i][i]=1;return ;
}
mat pow_mat(mat a,ll n)
{mat ret=unit;while(n){if(n&1){ret= ret*a;}n>>=1;a=a*a;}return ret;
}
int main()
{cin>>t;init_unit();while(t--){mat a;cin>>n>>k;for(int i=0;i<n;i++)for(int j=0;j<n;j++){cin>>a.m[i][j];}a = pow_mat(a, k);ll ans=0;for(int i=0;i<n;i++)ans+=a.m[i][i];ans%=9973;cout<<ans<<endl;}return 0;
}

这篇关于矩阵快速幂hdu 1575的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---