ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

矩阵快速幂(大佬们看看我的qwq)

2022-01-22 14:02:12  阅读:169  来源: 互联网

标签:node int ll 矩阵 ans 大佬 qwq 105


在学习矩阵快速幂之前,首先要了解什么是矩阵

上百度百科:

矩阵是一个按照长方阵列排列的复数或实数集合

形如 n * m m * p的两个矩阵可以进行矩阵乘法,也就是前一个矩阵的行数等于后一个的列数

在相乘时要将对应的行列的对应元素相乘后 相加,得到一个n * p 的矩阵

用公式就是这样

 

 因为太懒所以盗一张大佬的图

那么矩阵快速幂就显而易见,不过是把数换成了矩阵

附上代码和传送门

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int mod = 1e9 +7;
ll n,k;
struct node{
    ll a[105][105];
}s,ans;
node mul(node x,node y){
    node z;
    memset(z.a,0,sizeof(z.a));
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j){
            for(int k=1;k<=n;++k){
                z.a[i][j] += x.a[i][k] * y.a[k][j];
                z.a[i][j] %= mod;
            }
        }
    }
    return z;
} 
void q(ll k){
    while(k){
        if(k&1)ans = mul(ans,s);
        s = mul(s,s);
        k >>= 1;
    }
}
int main(){
    memset(ans.a,0,sizeof(ans.a));
    
    scanf("%lld%lld",&n,&k);
    
    for(int i=1;i<=n;++i){
        ans.a[i][i] = 1;
        for(int j=1;j<=n;++j){
            scanf("%lld",&s.a[i][j]);
        }
    }
    q(k);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j)printf("%lld ",ans.a[i][j]);
        putchar('\n');
    }
}

谢谢观看~

 

标签:node,int,ll,矩阵,ans,大佬,qwq,105
来源: https://www.cnblogs.com/wsxxs/p/15833544.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有