luoguP3390（矩阵快速幂模板题）

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

题意：矩阵快速幂模板题，思路和快速幂一致，只需提供矩阵的乘法即可。

AC代码：

#include<cstdio>
#include<cstring>
using namespace std;

typedef long long LL;

const int MOD=1e9+;
int n;
LL k;

struct Mat{
    LL m[][];
}a,e;

Mat mul(Mat& x,Mat& y){
    Mat res;
    memset(res.m,,sizeof(res.m));
    for(int i=;i<=n;++i)
        for(int j=;j<=n;++j)
            for(int l=;l<=n;++l){
                res.m[i][j]+=x.m[i][l]*y.m[l][j]%MOD;
                res.m[i][j]%=MOD;
            }
    return res;
}

Mat qpow(Mat& x,LL k){
    Mat ans=e;
    while(k){
        if(k&) ans=mul(ans,x);
        x=mul(x,x);
        k>>=;
    }
    return ans;
}

int main(){
    scanf("%d%lld",&n,&k);
    for(int i=;i<=n;++i)
        for(int j=;j<=n;++j)
            scanf("%lld",&a.m[i][j]);
    for(int i=;i<=n;++i)
        e.m[i][i]=;
    Mat ans=qpow(a,k);
    for(int i=;i<=n;++i){
        for(int j=;j<=n;++j)
            printf("%lld ",ans.m[i][j]);
        printf("\n");
    }
    return ;
}