題目背景

矩陣快速冪

題目描述

給定n*n的矩陣A，求A^k

輸入輸出格式

輸入格式：

?

第一行，n,k

第2至n+1行，每行n個(gè)數(shù)，第i+1行第j個(gè)數(shù)表示矩陣第i行第j列的元素

?

輸出格式：

?

輸出A^k

共n行，每行n個(gè)數(shù)，第i行第j個(gè)數(shù)表示矩陣第i行第j列的元素，每個(gè)元素模10^9+7

?

輸入輸出樣例

輸入樣例#1：
2 1 1 1 1 1 輸出樣例#1：
1 1 1 1

說(shuō)明

n<=100, k<=10^12, |矩陣元素|<=1000

?

//上板子！ #include<iostream> #include<cstdio> #define ll long long #define mod 1000000007using namespace std; ll n,m; struct node {ll a[101][101]; }ans,base;ll init() {ll x=0,f=1;char c=getchar();while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; }node mul(node a,node b) {node res;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){res.a[i][j]=0;for(int k=1;k<=n;k++)res.a[i][j]=(res.a[i][j]+a.a[i][k]*b.a[k][j])%mod;}return res; }node qw(node a,ll k) {node res=a;while(k){if(k&1) a=mul(a,res);res=mul(res,res);k>>=1;}return a; }int main() {n=init();m=init();for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){ans.a[i][j]=init();}m--;ans=qw(ans,m);for (int i=1;i<=n;i++){for (int j=1;j<n;j++) printf("%d ",ans.a[i][j]);printf("%d\n",ans.a[i][n]);} }

?

算法：矩陣快速冪

轉(zhuǎn)載于:https://www.cnblogs.com/L-Memory/p/6366930.html

總結(jié)

以上是生活随笔為你收集整理的P3390矩阵快速幂的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。