【知识学习】快速幂+矩阵快速幂+斐波那契数列及其推广

已于 2022-04-10 11:09:17 修改
阅读量211 收藏
点赞数
分类专栏： 知识学习 文章标签： acm oi
于 2020-11-27 22:34:27 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/cls1277/article/details/110247996
版权

前

这一篇blog从有想法到写完所有代码总共花了好长好长时间

虽然最后可能在这篇博文里体现不出来，但是真的好难

快速幂和矩阵快速幂比较简单，稍微带过直接放代码

斐波那契数列才是这次的主角

快速幂

位运算比较快 用了>>和&

LL qpow(LL a , LL b) {
	if(b==0)
		return 1;
	LL c = qpow(a,b>>1)%k;
	if(b&1)
		return c*c%k*a%k;
	return c*c%k;
}

矩阵快速幂

重载一下乘法，再写几个结构体函数，其实其他和快速幂一样

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
typedef long long LL;
#define PI acos(-1)
#define fin freopen("data.txt","r",stdin)
#define INF 2147483647
#define eps 1e-7
#define L 100005
#define mod 76
#define Fo(i,a,b) for(LL i=(a),_=(b); i<=_; i++)
#define Ro(i,b,a) for(LL i=(b),_=(a); i>=_; i--)
inline LL read() {
	LL x=0,f=1;char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+(c^48ll),c=getchar();
    return x*f;
}

LL n , k;
struct Mat {
	LL num[105][105];
	Mat() {
		memset(num , 0 , sizeof(num));
	}
	void init() {
		Fo(i,1,n)
			num[i][i] = 1;
	}
	void print() {
		Fo(i,1,n) {
			Fo(j,1,n)
				printf("%d ",num[i][j]);
			printf("\n");
		}	
	}
};
Mat A , E;
Mat operator * (const Mat &x , const Mat &y) {
	Mat z;
	Fo(i,1,n)
		Fo(j,1,n)
			Fo(k,1,n)
				z.num[i][j]=(z.num[i][j]+x.num[i][k]*y.num[k][j]%mod)%mod;
	return z;
}

Mat qpow(Mat x , LL y) {
	if(y==0)
		return E;
	if(y==1)
		return x;
	Mat z = qpow(x , y>>1);
	if(y&1)
		return z*z*x;
	return z*z;
}

int main() {
	E.init();
	n=read();k=read();
	Fo(i,1,n)
		Fo(j,1,n)
			A.num[i][j]=read();
	qpow(A,k).print();
	return 0;
}

斐波那契数列及其推广

有几种形式：（目前遇到的，以下分别标号123）

f [ 1 ] = 1 , f [ 2 ] = 1 , f [ i ] = f [ i − 1 ] + f [ i − 2 ] f[1]=1,f[2]=1,f[i]=f[i-1]+f[i-2] f[1]=1,f[2]=1,f[i]=f[i−1]+f[i−2]
f [ 1 ] = 1 , f [ 2 ] = 1 , f [ i ] = A × f [ i − 1 ] + B × f [ i − 2 ] f[1]=1,f[2]=1,f[i]=A×f[i-1]+B×f[i-2] f[1]=1,f[2]=1,f[i]=A×f[i−1]+B×f[i−2]
f [ 1 ] = C , f [ 2 ] = D , f [ i ] = A × f [ i − 1 ] + B × f [ i − 2 ] f[1]=C,f[2]=D,f[i]=A×f[i-1]+B×f[i-2] f[1]=C,f[2]=D,f[i]=A×f[i−1]+B×f[i−2]

分别对应题目：
LG1720
HDU1005
WIT新生赛J题

第1种

就是常规的斐波那契数列，矩阵快速幂做法如下图

在这里插入图片描述
其实没必要非得用最后一行的公式，用倒数第二行的公式或许更能变通（比如第2种）

代码实现：

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;
#define PI acos(-1)
#define fin freopen("data.txt","r",stdin)
#define INF 2147483647
#define eps 1e-7
#define L 100005
#define mod 1000000007
#define Fo(i,a,b) for(LL i=(a),_=(b); i<=_; i++)
#define Ro(i,b,a) for(LL i=(b),_=(a); i>=_; i--)
#define Ms(a,b) memset((a),(b),sizeof(a))
inline LL read() {
	LL x=0,f=1;char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+(c^48ll),c=getchar();
    return x*f;
}

LL n;
struct Mat {
	LL num[6][6];
	Mat() {
		Ms(num,0);
	}
	void init() {
		num[1][1]=num[1][2]=num[2][1]=1;
		num[2][2]=0;
	}
	void E() {
		Fo(i,1,n)
			num[i][i]=1; 
	}
	LL fib() {
		return num[2][1];
	}
};
Mat E , A;

Mat operator * (const Mat &x , const Mat &y) {
	Mat z;
	Fo(i,1,2)
		Fo(j,1,2)
			Fo(k,1,2)
				z.num[i][j]=z.num[i][j]+x.num[i][k]*y.num[k][j];
	return z;
}

Mat qpow(Mat x , LL y) {
	if(!y) return E;
	if(y==1) return x;
	Mat z=qpow(x,y>>1);
	if(y&1)
		return z*z*x;
	return z*z;
}

int main() {
	n=read();
	E.E(); A.init();
	printf("%lld.00",qpow(A,n).fib());
	return 0;
}

第2种

就是第1种的变通（画图有点丑）在这里插入图片描述
代码实现

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;
#define PI acos(-1)
#define fin freopen("data.txt","r",stdin)
#define INF 2147483647
#define eps 1e-7
#define L 100005
#define mod 7
#define Fo(i,a,b) for(LL i=(a),_=(b); i<=_; i++)
#define Ro(i,b,a) for(LL i=(b),_=(a); i>=_; i--)
#define Ms(a,b) memset((a),(b),sizeof(a))
inline LL read() {
	LL x=0,f=1;char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+(c^48ll),c=getchar();
    return x*f;
}

LL A , B , n;
struct Mat {
	LL num[15][15];
	Mat() {
		Ms(num,0);
	}
	void init(LL a , LL b) {
		num[1][1]=a; num[1][2]=1;
		num[2][1]=b; num[2][2]=0;
	}
	void K() {
		num[1][1]=num[1][2]=1;
	}
	void E() {
		Fo(i,1,n)
			num[i][i]=1;
	}
	LL fib() {
		return num[1][1]%mod;
	}
};
Mat E , K , P;

Mat operator * (const Mat &x , const Mat &y) {
	Mat z;
	Fo(i,1,2)
		Fo(j,1,2)
			Fo(k,1,2)
				z.num[i][j]=(z.num[i][j]+x.num[i][k]*y.num[k][j]%mod)%mod;
	return z;
}

Mat qpow(Mat x , LL y) {
	if(!y)  return E;
	if(y==1) return x;
	Mat z=qpow(x,y>>1);
	if(y&1)
		return z*z*x;
	return z*z; 
}

int main() {
	E.E(); K.K();
	while(1) {
		A=read(); B=read(); n=read();
		if(A==0&&B==0&&n==0)
			break;
		if(n==1||n==2) {
			printf("1\n");	
			continue;		
		}
		P.init(A,B);
		printf("%lld\n",(K*qpow(P,n-2)).fib()%mod);
	}
	return 0;
}

第3种

这个坑点有点多，但是学到一种求大数的斐波那契数列的方法

就是先求底数矩阵的个位数次方，然后底数矩阵自己求10次方，再做十位依次向前

因为我的做法是求得n-1次方，所以有很大的坑点

所以 直接用n就行！最好别写n-1因为要判断0啥的有亿点点麻烦
在这里插入图片描述

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;
#define PI acos(-1)
#define fin freopen("data.txt","r",stdin)
#define INF 2147483647
#define eps 1e-7
#define L 100005
#define Fo(i,a,b) for(LL i=(a),_=(b); i<=_; i++)
#define Ro(i,b,a) for(LL i=(b),_=(a); i>=_; i--)
#define Ms(a,b) memset((a),(b),sizeof(a))
inline LL read() {
	LL x=0,f=1;char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+(c^48ll),c=getchar();
    return x*f;
}

string n;
LL A , B , C , D , mod , len , num[500005] , flag;

struct Mat {
	int num[10][10];
	Mat() {
		Ms(num,0);
	}
	void K(LL a , LL b) {
		num[1][1]=a; num[1][2]=b;
		num[1][3]=num[2][1]=num[3][3]=1;
	} 
	void P(LL f0 , LL f1 , LL f2) {
		num[1][1]=f2;
		num[1][2]=num[2][1]=f1;
		num[2][2]=f0;
		num[3][1]=num[3][2]=num[3][3]=7;
	}
	void E() {
		Fo(i,1,3)
			num[i][i]=1;
	}
	LL fib() {
		return num[1][1]%mod;
	}
};
Mat K , P , ans;

Mat operator * (const Mat &x , const Mat &y) {
	Mat z;
	Fo(i,1,3)
		Fo(j,1,3)
			Fo(k,1,3)
				z.num[i][j]=(z.num[i][j]+x.num[i][k]*y.num[k][j]%mod)%mod;
	return z;
}

Mat E;
Mat qpow(Mat x , LL y) {
	if(!y) return E;
	if(y==1) return x;
	Mat z=qpow(x,y>>1);
	if(y&1)
		return z*z*x;
	return z*z;
}

int main() {
	A=read(); B=read(); C=read(); D=read(); cin>>n; mod=read();
	E.E(); K.K(A,B); P.P(D-C,C,D); ans.E();
	len = n.length();
	Fo(i,0,len-1)
		num[i+1]=n[i]-'0'; 
	Ro(i,len,1) {
		if(i==len) {
			LL top=i;
			while(num[top]==0) {
				top--;
				flag = 1;
			}
			if(!flag) {
				ans=qpow(K,num[i]-1);
				K=qpow(K,10);			
			} else {
				num[top]--;
				Fo(j,top+1,len)
					num[j]=9;
				ans = ans * qpow(K,num[i]);
				K = qpow(K,10);					
			}
			continue;
		}
		ans = ans * qpow(K,num[i]);
		K = qpow(K,10);			
	}
	ans=ans*P;
	printf("%lld",ans.fib()%mod);
	return 0;
}

后

好了，这篇blog会一直更新，遇到新的类型就会更新上去，希望能对你对我有所帮助

cls1277
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
斐波那契数列最优解-矩阵快速幂
ZZY的博客
05-14 644
斐波那契数列（Fibonacci sequence），又称黄金分割数列、因数学家列昂纳多·斐波那契（Leonardoda Fibonacci）以兔子繁殖为例子而引入，故又称为“兔子数列”，指的是这样一个数列：1、1、2、3、5、8、13、21、34、……在数学上，斐波那契数列以如下被以递推的方法定义：F(1)=1，F(2)=1, F(n)=F(n - 1)+F(n - 2)（n ≥ 3，n ∈ N*）在现代物理、准晶体结构、化学等领域，斐波纳契数列都有直接的应用。 1、递归 ...
一文彻底搞懂快速幂(原理实现、矩阵快速幂)
03-19
在实际应用中，矩阵快速幂常用于解决线性递推问题，如斐波那契数列、最短路径问题等。它通过预处理递推矩阵，然后快速求解任意步数的递推值，大大提高了计算速度。 这篇"一文彻底搞懂快速幂(原理实现、矩阵快速幂)...
【剑指offer】面试题 10：斐波那契数列及其推广
Jinlong_Xu的博客
06-29 537
题目描述 大家都知道斐波那契数列，现在要求输入一个整数n，请你输出斐波那契数列的第n项。 n 时间限制：1秒 空间限制：32768K 热度指数：153571 思路 1、用递归来实现，但是递归存在调用栈溢出的风险。 调用栈溢出是因为递归的计算中，每次函数调用在内存栈中分配空间，而每个进程的栈的容量都是有限的。
斐波那契数列的几种编程实现及一般推广
weixin_30586085的博客
12-18 215
斐波那契数列是一列规律很简单、明显的数列，它的第0项是0，第1项是1，第2项是1，依此类推，之后每一项是之前两数的和。首几个数是： 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, 610, 987, 1597, 2584, 4181, 6765, 10946 ……（OEISA000045） 编程实现 实现它最容易想到的方法，可...
通过矩阵快速幂在O(logN)时间求斐波那契数列及其推广形式
andyL_05的博客
04-17 660
通过矩阵快速幂在O(logN)时间求斐波那契数列及其推广形式 1. 斐波那契数列 这个数列非常经典了，递推公式为： f(n)=f(n−1)+f(n−2)f(n)=f(n-1)+f(n-2)f(n)=f(n−1)+f(n−2) 前几项为： 1, 2, 3, 5, 8, 13, 21, 34, 55… 如果我们要求斐波那契数列的第n项，我们一般用到的是O(N)复杂度的解法，也即构造一个长度为n的数组，...
由斐波那契数列引述到矩阵快速幂技巧
Harrison的学习博客
03-28 472
求斐波那契数列矩阵乘法的方法 1）斐波那契数列的线性求解（O(N)）的方式非常好理解 2)同时利用线性代数，也可以改写出另一种表示： |F(N),F(N-1)| = |F(2),F(1)| * 某个二阶矩阵的N-2次方 3)求出这个二阶矩阵，进而最快求出这个二阶矩阵的N-2次方 斐波那契数列的线性求解 public static int f1(int n){ if(n<1){ return 0; } if(n==1 || n==
矩阵快速幂求解斐波那契
04-13
二. 矩阵相乘 若A为n×k矩阵，B为k×m矩阵，则它们的乘积AB(有时记做A·B)将是一个n×m矩阵。...快速幂的思路就是：设A为矩阵，求A的N次方，N很大。例如：A的9次方 A^9 = A*A*A*A*A*A*A*A*A 【一个一
快速幂算法详解：高效计算大数幂的利器.zip
03-24
例如，在解决斐波那契数列、线性递推关系、矩阵快速幂等问题时，都能看到它的身影。熟练掌握快速幂算法，不仅能够提升代码性能，还能在解决复杂问题时提供思路。 通过阅读“快速幂算法详解：高效计算大数幂的利器....
C&C++解决斐波那契数列（兔子数列）问题
10-18
使用C&C++解决斐波那契数列（兔子繁殖）问题，可输入月份，输出从第0个月到所输入月份的兔子数量。入门级，C与C++混用。
2017广西邀请赛 Covering（矩阵快速幂）
信仰.的博客
08-31 575
题意：给你4*n的矩形，用1*2的方格铺满。 题解：经典问题，我们可以先用公式找出前几项，然后求递推式即可，但是n很大，所以要快速幂搞搞。 我们可以由前几项推得递推公式为：f[i]=f[i-1]+5*f[i-2]+f[i-3]-f[i-4]; 因此我们可以构造矩阵为： 1 5 1 -1 1 0 0 0 0 1 0 0 0 0 1 0       然后矩阵快速幂求解就好啦。
矩阵快速幂 CodeVS 1250 Fibonacci数列
unknown
06-04 214
矩阵快速幂
ACM-数论-矩阵快速幂 HDU6030
CS_W98的博客
04-29 315
这里是题面 r:red b:blue 【题意】 有一串珍珠，长度为n（1e18） 每个珍珠要不染色成红色，要不染色成蓝色。 要求任何连续素数长度的珍珠，都必须是红色个数&amp;gt;=蓝色个数 让你求出有多少种对这串珍珠的染色方案。 一开始以为要打质数表，但是后来发现不需要，套路特别像斐波那契数列。 要求任何连续素数长度的珍珠，都必须是红色个数&amp;gt;=蓝色个数 1.用f[...
2018年湘潭大学校赛 G题 又见斐波拉契（矩阵快速幂）
HJ
05-08 435
https://www.nowcoder.com/acm/contest/105/G根据递推公式 f[i+1] = f[i] + f[i-1] + i^3 + i^2 + i + 1 f[i] = f[i] (i+1)^3 = i^3 + 3*i^2 + 3*i + 1 (i+1)^2 = i^2 + 2*i + 1 (i+1) = i + 1 1 = 1 因为构造矩阵就
快速幂+矩阵优化斐波那契数列(超详细教程)
weixin_46503238的博客
03-24 6785
文章目录前言一、快速幂二、矩阵优化斐波那契数列1.矩阵相关知识2.斐波那契数列用矩阵表示3.O(log2n)的斐波那契数列总结 前言 我们首先讲解快速幂,然后利用快速幂优化矩阵相乘,将O(n)算法变为O(log2n),紧接着我们用矩阵,优化斐波那契数列! 一、快速幂 通常我们算一个数(a)的n次幂,我们需要计算n次,也就是n个a相乘,这样难免太过缓慢,于是有了快速幂,即不需要n次操作就可以算出! 举例:计算A的9次幂 通过上述操作想必大家明白,快速幂的思想也就是二分的思想! 代码实现: #include
2024年公关服务行业分析报告.pptx
07-09
行业报告
惠州市建筑物矢量数据（Shp格式+带高度）.txt
最新发布
07-09
因文件较多，数据存放网盘，txt文件内包含下载链接及提取码，永久有效。失效会第一时间进行补充。
Matlab语言教程及案例
07-09
项目总结 本文详细介绍了Matlab的基础知识和一个完整的数据分析与可视化案例。内容包括： 1. **基本操作**：基本算术运算和矩阵操作。 2. **绘图**：二维和三维绘图。 3. **函数和脚本**：编写和运行Matlab函数和脚本。 4. **案例分析**：生成模拟数据，进行数据预处理、绘图和统计分析。 通过这个项目，读者可以掌握Matlab的基本使用方法，并了解如何在Matlab中进行数据分析和可视化。这对于初学者入门Matlab编程、提升数据处理和分析能力非常有帮助。希望本文能为读者提供有价值的参考，提升其Matlab编程和数据分析技能。
KB8825英文数据手册
07-09
三星KB8825英文数据手册

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

cls1277

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

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

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

打赏作者

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

抵扣说明：

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

余额充值