51nod 1113 矩阵连乘快速幂模板（对100000007取模）

Werky_blog

于 2017-07-26 23:38:30 发布

阅读量1.7k

点赞数

CC 4.0 BY-SA版权

分类专栏：快速幂

本文链接：https://blog.csdn.net/qq_33997572/article/details/76167666

快速幂专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种利用矩阵快速幂解决大规模矩阵乘法问题的方法。通过定义矩阵结构、初始化矩阵及实现矩阵乘法和快速幂运算，可以在较短时间内计算出大指数下的矩阵结果。适用于求解递归关系和矩阵乘法等问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目：点击打开链接

#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
#define M 1000000007
using namespace std;
const int maxn=105;
int N;
struct mat{
    ll m[maxn][maxn];
};
mat A,I;
void init(){
    for(int i=0;i<N;i++)
    for(int j=0;j<N;j++){
        scanf("%lld",&A.m[i][j]);
        A.m[i][j]%=M;
        I.m[i][j]=(i==j);
    }
}
mat multi(mat a,mat b){
    mat c;
    for(int i=0;i<N;i++)
    for(int j=0;j<N;j++){
            c.m[i][j]=0;
            for(int k=0;k<N;k++)
                c.m[i][j]+=a.m[i][k]*b.m[k][j]%M;
            c.m[i][j]%=M;
    }
    return c;
}
mat power(mat A,int k){
    mat ans=I,p=A;
    while(k){
        if(k&1){
            ans=multi(ans,p);
            k--;
        }
        k>>=1;
        p=multi(p,p);
    }
    return ans;
}
int main(){
    int k;
    scanf("%d%d",&N,&k);
    init();
    mat ans=power(A,k);
    for(int i=0;i<N;i++){
        printf("%lld",ans.m[i][0]);
        for(int j=1;j<N;j++)
            printf(" %lld",ans.m[i][j]);
        printf("\n");
    }
    return 0;
}