Luogu3390 【模板】矩阵快速幂

原题链接：https://www.luogu.org/problemnew/show/P3390

【模板】矩阵快速幂

题目背景

矩阵快速幂

题目描述

给定n*n的矩阵A，求A^k

输入输出格式

输入格式：

第一行，n,k

第2至n+1行，每行n个数，第i+1行第j个数表示矩阵第i行第j列的元素

输出格式：

输出A^k

共n行，每行n个数，第i行第j个数表示矩阵第i行第j列的元素，每个元素模10^9+7

输入输出样例

输入样例#1：

2 1
1 1
1 1

输出样例#1：

1 1
1 1

说明

n<=100, k<=10^12, |矩阵元素|<=1000 算法：矩阵快速幂

题解

我就学习学习非递归的矩阵快速幂。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=105,mod=1e9+7;
ll n,k;
struct matrix{ll sq[M][M];matrix(){memset(sq,0,sizeof(sq));}}mat,one;
matrix operator *(matrix a,matrix b)
{
    matrix c;
    for(int i=1;i<=n;++i)for(int j=1;j<=n;++j)for(int k=1;k<=n;++k)
    (c.sq[i][j]+=a.sq[i][k]*b.sq[k][j]%mod)%=mod;
    return c;
}
matrix power(matrix x,ll p)
{
    matrix ans=one;
    for(;p;p>>=1,x=x*x)if(p&1)ans=ans*x;
    return ans;
}
void in()
{
    scanf("%lld%lld",&n,&k);
    for(int i=1;i<=n;++i)for(int j=1;j<=n;++j)scanf("%lld",&mat.sq[i][j]);
}
void ac()
{
    for(int i=1;i<=n;++i)one.sq[i][i]=1;
    mat=power(mat,k);
    for(int i=1;i<=n;++i,putchar(10))for(int j=1;j<=n;++j)printf("%lld ",mat.sq[i][j]);
}
int main(){in();ac();}