矩阵快速幂

最新推荐文章于 2023-06-22 16:37:05 发布
最新推荐文章于 2023-06-22 16:37:05 发布
阅读量261 收藏 1
点赞数 3
分类专栏： algorithm and data structure
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45530271/article/details/104183623
版权

目录

Learn

· 矩阵乘法

矩阵a[ i行 * k列 ] x 矩阵b[ k行 * j列 ] = 矩阵c[ i行 * j列 ]
（①第一个矩阵的列数必须等于第二个矩阵的行数，②结果矩阵的行数等于第一个矩阵的行数，列数等于第二个矩阵的列数）
c[ i ][ j ] = a[ i ][ 1 ] * b[ 1 ][ j ] + a[ i ][ 2 ] * b[ 2 ][ j ] +…+ a[ i ][ k ] * b[ k ][ j ]
在这里插入图片描述

· 矩阵快速幂

类比整数的快速幂，
[N * N的矩阵]1 × [N * N的矩阵]1 = [N * N的矩阵]2
[N * N的矩阵]2 × [N * N的矩阵]1 = [N * N的矩阵]3
所以矩阵也能快速幂，用struct封装矩阵，重载 * ^ 运算符即可实现。

struct Matrix{
	ll m[N][N];
	int size=3; //矩阵大小
	Matrix(){//根据题目进行初始化  
		memset(m,0,sizeof(m));
		m[1][1]=1;m[1][2]=2;m[1][3]=1;
		m[2][1]=1;m[2][2]=0;m[2][3]=0;
		m[3][1]=0;m[3][2]=0;m[3][3]=1;
	}
	
	void clear(){
		memset(m,0,sizeof(m));
		for(int i=1;i<=size;i++)
			m[i][i]=1;
	}
	
	void display(){
		cout<<"Matrix's begin:"<<endl;
		for(int i=1;i<=size;i++)
			for(int j=1;j<=size;j++)
				if(j<size) cout<<m[i][j]<<" ";
				else cout<<m[i][j]<<endl;
		cout<<"Matrix's end:"<<endl;
	}
	
	friend Matrix operator*(Matrix a,Matrix b){
		Matrix ans;
		for(int i=1;i<=size;i++) //size*size为矩阵大小 
			for(int j=1;j<=size;j++){
				ans.m[i][j]=0;
				for(int k=1;k<=size;k++)
					ans.m[i][j]=(ans.m[i][j]+a.m[i][k]*b.m[k][j])%mod;
			}
		return ans;
	}
	
	friend Matrix operator^(Matrix base,ll k){
		Matrix ans;
		ans.clear();
		while(k){
			if(k&1) ans=ans*base;
			base=base*base;
			k>>=1;
		}
		return ans;
	}
};

//使用方法：
//Matrix base;
//base=base^n;

Practice

· [HDU4990] Reading comprehension

题目链接

题意：给出1<=n,m<=1e9，序列a定义为
a[0]=0，a[i为奇数]=a[i-1]*2+1，a[i为偶数]=a[i-1]*2，计算第n项模m的结果

思路：a[0]=0，a[1]=1，a[2]=2，a[3]=5，a[4]=10，a[5]=21，a[6]=42，a[7]=85，a[8]=170
看出a(i)=a(i-1)+2*a(i-2)+1，(i>2)
所以可以如下构造矩阵，用矩阵快速幂求解

在这里插入图片描述

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<map>
#include<algorithm>
#include<queue>
#include<math.h>
#include<set>
#include<vector>
using namespace std;
#define Inf 0x7fffffff
typedef long long ll;
const int N=100007;
ll n,mod;
int size=3;
struct Matrix{
	ll m[15][15];
	
	Matrix(){
		memset(m,0,sizeof(m));
		m[1][1]=1;m[1][2]=2;m[1][3]=1;
		m[2][1]=1;m[2][2]=0;m[2][3]=0;
		m[3][1]=0;m[3][2]=0;m[3][3]=1;
	}
	
	void clear(){
		memset(m,0,sizeof(m));
		for(int i=1;i<=size;i++)
			m[i][i]=1;
	}
	
	void display(){
		cout<<"Matrix's begin:"<<endl;
		for(int i=1;i<=size;i++)
			for(int j=1;j<=size;j++)
				if(j<size) cout<<m[i][j]<<" ";
				else cout<<m[i][j]<<endl;
		cout<<"Matrix's end:"<<endl;
	}
	
	friend Matrix operator*(Matrix a,Matrix b){
		Matrix ans;
		for(int i=1;i<=size;i++) //size*size为矩阵大小 
			for(int j=1;j<=size;j++){
				ans.m[i][j]=0;
				for(int k=1;k<=size;k++)
					ans.m[i][j]+=a.m[i][k]*b.m[k][j]%mod;
			}
		return ans;
	}
	
	friend Matrix operator^(Matrix base,ll k){
		Matrix ans;
		ans.clear();
		while(k){
			if(k&1) ans=ans*base;
			base=base*base;
			k>>=1;
		}
		return ans;
	}
};

int main(){
	while(~scanf("%lld%lld",&n,&mod)){
		Matrix base;
		base=base^(n-1);
		ll ans=0;
		ll a[4];a[1]=1;a[2]=0;a[3]=1;
		for(int i=1;i<=size;i++)		
			ans=(ans+a[i]*base.m[1][i])%mod;
		cout<<ans<<"\n";
	}
}

· [牛客寒假集训营1J题] u’s的影响力

题目链接

题意：给出1<=n,x,y,a,b<=1e12，f(1)=x，f(2)=y，f(i) = f(i-1) * f(i-2) * ab
计算 f(n)%1000000007的值

思路：
在这里插入图片描述
需要注意的一点是，矩阵里的数最终都是要作为指数，因为模数mod是质数，所以对作为指数的数要%(mod-1)，最后算快速幂时再quick_pow(x,e+mod-1) [欧拉降幂]

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<map>
#include<algorithm>
#include<queue>
#include<math.h>
#include<set>
#include<vector>
using namespace std;
#define Inf 0x7fffffff
typedef long long ll;
const ll mod=1000000007;
int size=3;
ll n,x,y,a,b,arr[4];
struct Matrix{
    ll m[15][15];
     
    Matrix(){//根据题目进行初始化 
        memset(m,0,sizeof(m));
        m[1][1]=1;m[1][2]=1;m[1][3]=1;
        m[2][1]=1;m[2][2]=0;m[2][3]=0;
        m[3][1]=0;m[3][2]=0;m[3][3]=1;
    }
     
    void clear(){
        memset(m,0,sizeof(m));
        for(int i=1;i<=size;i++)
            m[i][i]=1;
    }
     
    void display(){
        cout<<"Matrix's begin:"<<endl;
        for(int i=1;i<=size;i++)
            for(int j=1;j<=size;j++)
                if(j<size) cout<<m[i][j]<<" ";
                else cout<<m[i][j]<<endl;
        cout<<"Matrix's end:"<<endl;
    }
     
    friend Matrix operator*(Matrix a,Matrix b){
        Matrix ans;
        for(int i=1;i<=size;i++) //size*size为矩阵大小
            for(int j=1;j<=size;j++){
                ans.m[i][j]=0;
                for(int k=1;k<=size;k++)
                    ans.m[i][j]=(ans.m[i][j]+a.m[i][k]*b.m[k][j])%(mod-1);
            }
        return ans;
    }
     
    friend Matrix operator^(Matrix base,ll k){
        Matrix ans;
        ans.clear();
        while(k){
            if(k&1) ans=ans*base;
            base=base*base;
            k>>=1;
        }
        return ans;
    }
};
ll qpow(ll x,ll e){
    ll re=1;
    while(e){
        if(e&1){
            re=re*x%mod;
        }
        x=x*x%mod;
        e>>=1;
    }
    return re%mod;
}
 
int main(){
    cin>>n>>x>>y>>a>>b;
    x%=mod;y%=mod;a%=mod;b=b%(mod-1);
    if(n==1){
        cout<<x;
        return 0;
    }
    if(n==2){
        cout<<y;
        return 0;
    }
    Matrix base;
    base=base^(n-2);
    ll ans=1;
    ans=(ans*qpow(y,base.m[1][1]+mod-1))%mod;
    ans=(ans*qpow(x,base.m[1][2]+mod-1))%mod;
    ll t=qpow(a,b+mod-1);
    ans=(ans*qpow(t,base.m[1][3]+mod-1))%mod;
    cout<<ans%mod;
}

· [Poj3613] Cow Relays

题目链接

题意：给出一个k (<=2*t)个点，t (<=100)条边的无向连通图，求s点到e点经过n (<=1e6)条边的最短距离

思路：
思考一下，如果一个矩阵，表示走k条边后，一张图的点与点的最短路径，(a,b)表示从a到b的最短路径，然后我们把它与自己，按照矩阵乘法的格式“相乘”，把其中的乘改为取min，ans.m[i][j] = min(ans.m[i][j],a.m[i][k]+b.m[k][j])【类似于Floyd】。这样得到的是走k+k条边的矩阵，同理，一个走了a条边的矩阵×一个走了b条边的矩阵，得到的是走了a+b条边的矩阵。
所以对点离散化，用矩阵快速幂的方法，修改一下opereator重载的 * 式子求basen即可。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<map>
#include<algorithm>
#include<queue>
#include<math.h>
#include<set>
#include<vector>
using namespace std;
typedef long long ll;
const int N=207;
int m,s,t,l[N],r[N];
ll w[N],n,Inf;
int U[N],cntu,size;
struct Matrix{
	ll m[N][N];
	
	Matrix(){//根据题目进行初始化  
		memset(m,63,sizeof(m));
	}
	
	void clear(){
		memset(m,0,sizeof(m));
		for(int i=1;i<=size;i++)
			m[i][i]=1;
	}
	
	void display(){
		cout<<"Matrix's begin:"<<endl;
		for(int i=1;i<=size;i++)
			for(int j=1;j<=size;j++)
				if(j<size) cout<<m[i][j]<<" ";
				else cout<<m[i][j]<<endl;
		cout<<"Matrix's end:"<<endl;
	}
	
	friend Matrix operator*(Matrix a,Matrix b){
		Matrix ans;
		for(int i=1;i<=size;i++) //size*size为矩阵大小 
			for(int j=1;j<=size;j++)
				for(int k=1;k<=size;k++){
					if(ans.m[i][j]==Inf)ans.m[i][j]=a.m[i][k]+b.m[k][j];
					else ans.m[i][j]=min(ans.m[i][j],a.m[i][k]+b.m[k][j]);
				}
		return ans;
	}
	
	friend Matrix operator^(Matrix base,ll k){
		k--;
		Matrix ans=base;
		while(k){
			if(k&1) ans=ans*base;
			base=base*base;
			k>>=1;
		}
		return ans;
	}
};


int main(){
    cin>>n>>m>>s>>t;
    for(int i=1;i<=m;i++){
		scanf("%lld%d%d",&w[i],&l[i],&r[i]);
		U[++cntu]=l[i];
		U[++cntu]=r[i];
	}
    sort(U+1,U+1+cntu);
    size=cntu=unique(U+1,U+1+cntu)-U-1;
    Matrix base;
	Inf=base.m[0][0];
	for(int i=1;i<=m;i++){
    	l[i]=lower_bound(U+1,U+1+cntu,l[i])-U;
    	r[i]=lower_bound(U+1,U+1+cntu,r[i])-U;
   		base.m[l[i]][r[i]]=base.m[r[i]][l[i]]=w[i];
    }
    s=lower_bound(U+1,U+1+cntu,s)-U;
    t=lower_bound(U+1,U+1+cntu,t)-U;
    base=base^n;
    cout<<base.m[s][t];
}
linkscx
关注
  • 3
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
矩阵快速幂详解
Gene_I must wait for the sunrise
06-01 1149
前言 如果你对这篇文章可感兴趣，可以点击「【访客必读 - 指引页】一文囊括主页内所有高质量博客」，查看完整博客分类与对应链接。 设想这样一个场景，面试官给了你一道算法题，你很快确定这是一道递推问题，并给出了 O(n) 的解题方法，然而面试官却继续问：“还能继续优化吗？” 这样类似的场景并不少见，因为算法不仅追求「正确」，还追求「效率」，而这也正是优化方法的意义。 本文即将介绍的「矩阵快速幂」便是一种常见的优化递推的方法，能将 O(n) 的递推过程加速至 O(log(n))，使效率显著提升。 矩阵运算 首先我
【蓝桥杯】简单数论2——快速幂&矩阵快速幂
荷叶田田的博客
01-22 3160
蓝桥杯之简单数论：快速幂和矩阵快速幂
矩阵快速幂算法
果光的博客
03-10 3805
矩阵快速幂算法、利用矩阵快速幂计算斐波那契数列 矩阵快速幂的本质就是 矩阵+快速幂，思路没什么变化。 快速幂的思路见之前的 快速幂介绍，这里就不多说了。 至于矩阵的话，知道矩阵是啥，怎么算乘法就可以了。 相关矩阵基础 注：已经了解矩阵的可以跳过本部分内容。 简单来说 n∗mn*mn∗m 的矩阵就是个 nnn 行 mmm 列 的大方块。 矩阵要进行乘法运算必须保证第一个矩阵的列数=第二个矩阵的行数。
快速幂 矩阵快速幂
issey的博客
10-22 5871
前言：好像没啥好写的，链接可能还没有更新完 快速幂 快速幂，用于解决当n很大时的情况。通常与取模同时应用。 用最笨的方法求,即。时间复杂度为,而快速幂（附带取模），可以将时间复杂度降低为。 利用倍增的思想，例如，等于，又等于,即,如果n为奇数，,那么换成代码就是： ll Powermod(ll a,ll b) { ll ans=1; a=a%mod; while(b>0) { if(b%2==1) ans=(ans*a)%...
矩阵快速幂及其应用
weixin_43667715的博客
06-22 400
在学习矩阵快速幂之前，我们必须先了解学习一下矩阵乘法的事情。
矩阵快速幂总结
热门推荐
努力
11-25 4万+
基础知识：(会基础的直接看应用部分)(1)矩阵乘法简单的说矩阵就是二维数组，数存在里面，矩阵乘法的规则:A*B=C其中c[i][j]为A的第i行与B的第j列对应乘积的和，即:代码:const int N=100; int c[N][N]; void multi(int a[][N],int b[][N],int n)//n是矩阵大小，n&lt;N { memset(c,0,sizeof c...
矩阵快速幂加速递推
qq_66500237的博客
05-16 229
可以看到，我们用matrixpower函数来计算矩阵的n次方，这里用了一种二分法的思想，将n拆成二进制后，每一位对应一个矩阵的乘积。矩阵快速幂是一种优化递推的算法，可以用来求解递推式的第n项。将递推式转化为矩阵形式：[F(n), F(n-1)] = [F(n-1), F(n-2)] * [1, 1;首先，斐波那契数列的递推式为：F(n) = F(n-1) + F(n-2)，其中F(0)=0，F(1)=1。1, 0]]的n次方，来得到第n项斐波那契数。以下是一个斐波那契数列的例子，用矩阵快速幂来求解第n项。
矩阵快速幂模板以及其应用矩阵加速递推
M3ng@L的博客
08-27 847
计算这样的递推次数很多的式子，我们需要一个更简化更快捷的计算方式，也就是矩阵快速幂。我们只需要构造最小限度能解决问题的矩阵就好，无论是矩阵初始值还是转移矩阵，所以还是需要具体问题具体分析地来构造矩阵。的矩阵，其实列向量还是行向量都可以，只要转移矩阵也是对应的即可）那么根据转移矩阵表现出来的与原数列的关系，我们可以这样来看。的矩阵初始值浪费了计算空间，所以构造其他形式的初始矩阵。的初始矩阵的话，对应数列的递推关系是不完整的；得到的新矩阵的第一项数据即是斐波那契数列的第。把矩阵乘法展开，得到。......
C++题解：矩阵快速幂 求 斐波那契数列
Keven_11的博客
08-18 1253
一般先构造最后表示答案的矩阵，比如这里我们构造一个 1×2 的矩阵里边填上最开始的两项[a1​,a2​] ，接下来我们考虑把它往后推到 [a2​,a3​] ，那么我们需要构造一个矩阵，让 [a1​,a2​] 乘这个矩阵得到 [a2​,a3​]。依次看后边矩阵的每一项，看它可以由前边矩阵里边的每一项怎么凑出来，会发现 a2​=a2​,a3​=a1​+a2​，所以构造一个 2×2 的矩阵 A=[0,1，1,1​] ，这样就可以按照矩阵乘法递推，然后用矩阵快速幂解决问题了。这里n的范围极大，只能用矩阵快速幂。..
yxy版c++教程 浅谈矩阵快速幂
01-09
矩阵快速幂与C++实现 矩阵快速幂是一种高效的算法，用于计算矩阵的幂运算。该算法基于二进制表示法，通过将幂运算分解为小幂运算，从而大大减少了计算次数。在C++中，矩阵快速幂可以通过模板实现，提高了计算效率。...
矩阵快速幂_矩阵快速幂_
09-30
矩阵快速幂是一种高效的算法，主要用于解决大整数的幂运算问题。在计算机科学，特别是算法竞赛和数值计算中，它有着广泛的应用。本算法的核心思想是利用矩阵乘法的结合律和分配律来减少计算次数，从而提高计算效率。...
一文彻底搞懂快速幂(原理实现、矩阵快速幂)
03-19
在本文中，我们将深入探讨快速幂的原理，以及其在矩阵运算中的应用——矩阵快速幂。 快速幂算法的核心在于二进制分解。对于一个大整数n，我们可以通过将其转换为二进制表示，然后逐位处理，来计算a^n。例如，如果n=...
矩阵快速幂求解斐波那契
04-13
矩阵快速幂 幂又称乘方。表示一个数字乘若干次的形式,如n个a相乘的幂为a^n ，或称a^n为a的n次幂。a称为幂的底数，n称为幂的指数。 快速幂的思路就是：设A为矩阵，求A的N次方，N很大。例如：A的9次方 A^9 = A*A*A...
【202309】毕马威-中国第六届领先汽车科技50_109页.pdf
07-28
【202309】毕马威-中国第六届领先汽车科技50_109页.pdf
毕业设计—流媒体视频直播服务器-可执行内含文档代码-可执行内含文档代码.zip
07-28
毕业设计—流媒体视频直播服务器-可执行内含文档代码-可执行内含文档代码.zip
【202206】中国汽车潮流观察报告-汽车之家研究院_54页.pdf
07-28
【202206】中国汽车潮流观察报告-汽车之家研究院_54页.pdf
tour-project【程序员VIP专用】.zip
最新发布
07-28
【项目简介】 HTML5+CSS实现雪花旅行社
矩阵快速幂详解：从概念到应用
"yxy版c++教程 浅谈矩阵快速幂" 本文主要介绍了矩阵快速幂这一高效的算法，它是快速幂思想在矩阵运算上的应用。快速幂算法的基本原理是利用二进制分解，大大减少计算一个数的幂次所需的乘法次数。在数论中，快速幂...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

linkscx

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值