备战蓝桥杯---数学之矩阵快速幂基础

我们先不妨看一道题:

看见n的数据范围就知道直接按以前的递归写肯定狗带,那我们有什么其他的方法吗?

下面是分析:

我们就拿斐波那契数列试试手吧:

下面是AC代码,可以当作模板记:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int m,n,mod=1e9+7;
struct node{
	int m[100][100];
}ans,res;
node mul(node a,node b){
	node tmp;
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			tmp.m[i][j]=0;
		}
	}
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			for(int k=0;k<n;k++){
				tmp.m[i][j]=(tmp.m[i][j]+a.m[i][k]*b.m[k][j])%mod;
			}
		}
	}
	return tmp;
}
void quickpower(int m,int n){
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			if(i==j) ans.m[i][j]=1;
			else ans.m[i][j]=0;
		}
	}
	while(m){
		if(m&1) ans=mul(ans,res);
		res=mul(res,res);
		m=m>>1; 
	}
}
signed main(){
	cin>>m;
	m-=2;
	n=2;
	res.m[0][0]=0;
	res.m[1][0]=1;
	res.m[1][1]=1;
	res.m[0][1]=1;
	if(m<0) cout<<1;
	else {quickpower(m,n);
	cout<<(ans.m[1][0]+ans.m[1][1])%mod;}
}
相关推荐
怀澈1228 分钟前
高性能服务器模型之Reactor(单线程版本)
linux·服务器·网络·c++
chnming198731 分钟前
STL关联式容器之set
开发语言·c++
带多刺的玫瑰35 分钟前
Leecode刷题C语言之统计不是特殊数字的数字数量
java·c语言·算法
爱敲代码的憨仔36 分钟前
《线性代数的本质》
线性代数·算法·决策树
威桑42 分钟前
MinGW 与 MSVC 的区别与联系及相关特性分析
c++·mingw·msvc
熬夜学编程的小王1 小时前
【C++篇】深度解析 C++ List 容器:底层设计与实现揭秘
开发语言·数据结构·c++·stl·list
yigan_Eins1 小时前
【数论】莫比乌斯函数及其反演
c++·经验分享·算法
Mr.131 小时前
什么是 C++ 中的初始化列表?它的作用是什么?初始化列表和在构造函数体内赋值有什么区别?
开发语言·c++
阿史大杯茶1 小时前
AtCoder Beginner Contest 381(ABCDEF 题)视频讲解
数据结构·c++·算法
C++忠实粉丝1 小时前
计算机网络socket编程(3)_UDP网络编程实现简单聊天室
linux·网络·c++·网络协议·计算机网络·udp