矩阵加速(数列)

最新推荐文章于 2024-05-26 17:31:07 发布
最新推荐文章于 2024-05-26 17:31:07 发布
阅读量401 收藏 2
点赞数
分类专栏： Daily algorithm
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/CC_1012/article/details/91443638
版权
本文介绍了如何使用矩阵加速技术来高效地计算斐波那契数列及其广义形式的问题。具体包括P1962斐波那契数列模运算，P1349广义斐波那契数列的余数计算，以及P1939矩阵加速模板的应用。通过矩阵乘法，可以解决大范围内的数列计算问题。
摘要由CSDN通过智能技术生成

P1962 斐波那契数列

题目背景

大家都知道，斐波那契数列是满足如下性质的一个数列：

• f(1) = 1

• f(2) = 1

• f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数)

题目描述

请你求出 f(n) mod 1000000007 的值。

输入输出格式

输入格式：

 

·第 1 行：一个整数 n

 

输出格式：

 

第 1 行： f(n) mod 1000000007 的值

 

题解

/*
落谷:p1962 
题意：求斐波那契数列f(n)mod1e9+7 n在long long范围内 
*/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL mod=1e9+7;
LL n,p;
struct node
{
    LL m[3][3];
    node(){
        memset(m,0,sizeof(m));
    }
    friend node operator *(node x,node y){
        node t;
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                t.m[i][j]=0;
            }
        }
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                for(int k=1;k<=2;k++){
                	//尽量这样mod规范些 
                    t.m[i][j]+=((x.m[i][k]%mod)*(y.m[k][j]%mod))%mod;
                }
            }
        }
        return t;
    }
}a,ans;
node binary(node x,LL y)
{
    while(y){
        if(y&1) ans=ans*x;
        x=x*x;
        y>>=1;
    }
    return ans;
    
}
void Init()
{
    a.m[1][1]=1;
    a.m[1][2]=1;
    a.m[2][1]=1;
    a.m[2][2]=0;
    
    ans.m[1][1]=1;
    ans.m[1][2]=1;
}
int main()
{
    Init();
    LL n;scanf("%lld",&n);
    if(n<=2) printf("1\n");
    else{
        node t=binary(a,n-2);
        printf("%lld\n",t.m[1][1]%mod);
    } 
    return 0;
}

 

 

P1349 广义斐波那契数列

题目描述

广义的斐波那契数列是指形如an=p\times a_{n-1}+q\times a_{n-2}an=p×an−1​+q×an−2​的数列。今给定数列的两系数pp和qq，以及数列的最前两项a_1a1​和a_2a2​，另给出两个整数nn和mm，试求数列的第nn项a_nan​除以mm的余数。

输入输出格式

输入格式：

 

输入包含一行6个整数。依次是pp,qq,a_1a1​,a_2a2​,nn,mm，其中在pp,qq,a_1a1​,a_2a2​整数范围内，nn和mm在长整数范围内。

 

输出格式：

 

输出包含一行一个整数，即a_nan​除以mm的余数。

 

输入输出样例

输入样例#1： 复制

1 1 1 1 10 7

输出样例#1： 复制

6

说明

数列第10项是55，除以7的余数为6。

 

和第一个的区别是要改下基本矩阵  想办法相乘为第n项就可以

p 1

1 0    这样的一个矩阵


#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
LL mod;
LL p,q,a1,a2,n,m;
struct node
{
    LL m[3][3];
    node(){
        memset(m,0,sizeof(m));
    }
    friend node operator *(node x,node y){
        node t;
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                t.m[i][j]=0;
            }
        }
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                for(int k=1;k<=2;k++){
                    t.m[i][j]+=(x.m[i][k]%mod*y.m[k][j]%mod)%mod;
                }
            }
        }
        return t;
    }
}a,e,b,ans;
node binary(node x,LL y)
{
    ;
    while(y){
        if(y&1) ans=ans*x;
        x=x*x;
        y>>=1;
    }
    return ans;
    
}
void Init()
{
    a.m[1][1]=p;
    a.m[1][2]=1;
    a.m[2][1]=q;
    a.m[2][2]=0;
    
    ans.m[1][1]=a2;
    ans.m[1][2]=a1;
    mod=m;
}
int main()
{
    
    scanf("%lld %lld %lld %lld %lld %lld",&p,&q,&a1,&a2,&n,&m);
    Init();
    if(n==1) printf("%lld\n",a1);
    else if(n==2) printf("%lld\n",a2);
    else{
        node t=binary(a,n-2);
        printf("%lld\n",t.m[1][1]%mod);
    } 
    return 0;
}

 

 

P1939 【模板】矩阵加速（数列）

题目描述

a[1]=a[2]=a[3]=1

a[x]=a[x-3]+a[x-1] (x>3)

求a数列的第n项对1000000007（10^9+7）取余的值。

输入输出格式

输入格式：

 

第一行一个整数T，表示询问个数。

以下T行，每行一个正整数n。

 

输出格式：

 

每行输出一个非负整数表示答案。

 

输入输出样例

输入样例#1： 复制

3
6
8
10

输出样例#1： 复制

4
9
19

说明

对于30%的数据 n<=100；

对于60%的数据 n<=2*10^7；

对于100%的数据 T<=100，n<=2*10^9；

 

 

 

 

 

推导过程

 

 

/*
落谷：p3390
给定n*n的矩阵A，求A^k %p 
*/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL mod=1e9+7;
LL p,q,a1,a2,n,m;
struct node
{
	LL m[4][4];
	node(){
		memset(m,0,sizeof(m));
	}
	friend node operator *(node x,node y){
		node t;
		for(int i=1;i<=3;i++){
			for(int j=1;j<=3;j++){
				t.m[i][j]=0;
			}
		}
		for(int i=1;i<=3;i++){
			for(int j=1;j<=3;j++){
				for(int k=1;k<=3;k++){
					t.m[i][j]+=(x.m[i][k]%mod*y.m[k][j]%mod)%mod;
				}
			}
		}
		return t;
	}
}a,ans;
node binary(node x,LL y)
{
	while(y){
		if(y&1) ans=ans*x;
		x=x*x;
		y>>=1;
	}
	return ans;
	
}
void Init()
{
	memset(a.m,0,sizeof(a.m));
	memset(ans.m,0,sizeof(ans.m));
	a.m[1][1]=1;
	a.m[1][2]=1;
	a.m[2][3]=1;
	a.m[3][1]=1;
	
	ans.m[1][1]=1;
	ans.m[1][2]=1;
	ans.m[1][3]=1;
	
}
int main()
{
	int T;scanf("%d",&T);
	while(T--){
		Init();
		LL n;
		scanf("%lld",&n);
		if(n<=3) printf("1\n");
		else{
			node t=binary(a,n-3);
			printf("%lld\n",t.m[1][1]%mod);
		}
	} 
	return 0;
}

例题4：和P1349 广义斐波那契数列一样但是n的范围为：

 

十进制拆分不过就是b是10的倍数的话，我们让a每次变为原来的10次幂，b变为原来的1/10

他要是给n是一个16进制串，那咱们就搞个16进制拆分快速幂

 

 
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e6+100;
LL mod;
LL p,q,a1,a2,len;
char s[N];
struct node
{
    LL m[3][3];
    node(){
        memset(m,0,sizeof(m));
    }
    friend node operator *(node x,node y){
        node t;
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                t.m[i][j]=0;
            }
        }
        for(int i=1;i<=2;i++){
            for(int j=1;j<=2;j++){
                for(int k=1;k<=2;k++){
                    t.m[i][j]+=(x.m[i][k]%mod*y.m[k][j]%mod)%mod;
                }
            }
        }
        return t;
    }
}a,e,b,ans;
node binary(node x,LL y)
{
    node ret;
    for(int i=1;i<=2;i++) ret.m[i][i]=1;
    while(y){
        if(y&1) ret=ret*x;
        x=x*x;
        y>>=1;
    }
    return ret;
}
node solve(node A)
{
    node ret;
    for(int i=1;i<=2;i++) ret.m[i][i]=1;
    for(int i=len;i>=1;i--){
        int x=s[i]-'0';
        ret=ret*binary(A,x);
        A=binary(A,10);
    }
    return ret;
}
void Init()
{
    a.m[1][1]=p;
    a.m[1][2]=1;
    a.m[2][1]=q;
    a.m[2][2]=0;
    
    ans.m[1][1]=a2;
    ans.m[1][2]=a1;
}
int main()
{
    scanf("%lld %lld %lld %lld",&a1,&a2,&p,&q);
    scanf("%s",s+1);
    scanf("%lld",&mod);
    len=strlen(s+1);
    Init();
    node t=ans*solve(a);
    printf("%lld\n",t.m[1][2]%mod);
    return 0;
}

 

「已注销」
关注
专栏目录
[矩阵加速] 数列
烟尘的博客
04-02 339
文章目录题目描述解题思路 题目描述 数列 解题思路 首先说一下矩阵乘法： 假设矩阵AAA和矩阵BBB做乘法得到矩阵CCC，则必须满足矩阵AAA的列数等于矩阵BBB的行数 即如果矩阵AAA有nnn行mmm列，则矩阵BBB有mmm行ppp列，矩阵CCC就有nnn行ppp列 那么矩阵乘法怎么做呢？ 公式：C[i][j]=∑i=1n∑j=1p∑k=1mA[i][k]∗B[k][j]C[i][j] = \...
浅谈矩阵加速——以时间复杂度为O(log n）的算法实现裴波那契数列第n项及前n之和使用矩阵加速法的优化求法
C20200521的博客
04-01 743
首先请连矩阵乘法乘法都还没有了解的同学简单看一下这篇博客： https://blog.csdn.net/weixin_44049566/article/details/88945949 首先直接暴力求使用O(n)的时间复杂度肯定是不行的，所以我们应该使用更优的时间复杂度。 设f(n)为裴波那契数列第n项。让我们来构造两个矩阵： 和. 现在我们不妨将两个矩阵相乘，化简过后可以得到：,...
P1939 【模板】矩阵加速（数列）
hjx
08-19 380
luogu 传送门 矩阵乘法快速幂，注意记录过程量； 还有就是如果前面已经算过了，直接赋值就好了。#include<iostream> #include<cstdio> #include<algorithm> #include<string> #include<cstring> #define LL long long #define MOD 1000000007 using namespac
【模板】矩阵加速（数列）
weixin_30608131的博客
04-18 97
洛谷 1939 a[1]=a[2]=a[3]=1 a[x]=a[x-3]+a[x-1] (x>3) 求a数列的第n项对1000000007（10^9+7）取余的值。 1 #include<cstdio> 2 #include<algorithm> 3 #include<cstring> 4 #define LL long lo...
矩阵加速(递推式)(数列)
姬小野的博客
07-26 1334
P1939 【模板】矩阵加速（数列） 题目描述 a[1]=a[2]=a[3]=1 a[x]=a[x-3]+a[x-1] (x&gt;3) 求a数列的第n项对1000000007（10^9+7）取余的值。 输入输出格式 输入格式： 第一行一个整数T，表示询问个数。 以下T行，每行一个正整数n。 输出格式： 每行输出一个非负整数表示答案。 输入输出样例 输入样例#1： 3 ...
洛谷 1939 【模板】矩阵加速（数列）
苦鬼·。
08-08 1026
题目描述 a[1]=a[2]=a[3]=1 a[x]=a[x-3]+a[x-1] (x>3) 求a数列的第n项对1000000007（10^9+7）取余的值。 输入输出格式 输入格式： 第一行一个整数T，表示询问个数。 以下T行，每行一个正整数n。 输出格式： 每行输出一个非负整数表示答案。 输入输出样例 输入样例#1： 3 6 8 10 输出样例#1： 4 9
python 使用矩阵加速计算斐波那契数列第N项
05-21
使用python实现的的计算线性齐次递推数列第N项的算法
矩阵加速技术在数列与算法中的应用探析
"浅谈矩阵加速.pdf" 矩阵加速是一种在计算领域中用于优化线性问题的高效算法，尤其在处理线性递推关系时表现出色。矩阵加速的核心是矩阵快速幂技术，它允许我们以较低的时间复杂度求解矩阵的幂。这种技术广泛应用于...
【洛谷P1939】[模板]矩阵加速(数列)【矩阵乘法】
绫
12-15 255
题目描述 linklinklink 已知一个数列 a，它满足： ax=1(x<=1,2,3)a_x=1(x<=1,2,3)ax​=1(x<=1,2,3) ax=ax−1+ax−3(x>=4)a_x=a_{x-1}+a_{x-3}(x>=4)ax​=ax−1​+ax−3​(x>=4) 求 aaa 数列的第 nnn 项对109+710^9+7109+7取余的值。 输入格式 第一行一个整数 TTT，表示询问个数。 以下 TTT 行，每行一个正整数 nnn。 输出格式 每行输出一
矩阵加速（数列）
LELElele01的博客
11-12 203
矩阵加速（数列） 矩阵优化是个很不错的工具，但struct的坑太多了！好难改. 首先我们知道初始状态 [f[3]f[2]f[1]]\begin{bmatrix}f[3]&\\f[2]&\\f[1]\end{bmatrix}⎣⎡​f[3]f[2]f[1]​​⎦⎤​ [f[i]f[i−1]f[i−2]]\begin{bmatrix}f[i]&\\f[i-1]&\\f[i...
洛谷 P1939 【模板】矩阵加速（数列）：优化递推式的方法——矩阵快速幂
Karshilov的博客
12-15 1748
在大多数情况下，O(n)的效率都是值得骄傲的，然而，有时候并不是，比如如何在一秒钟内算出一个递推式的第1e9项，很明显O(n)不行了。 然而常数级又不太现实，除非你的数学非常好，这题又比较简单，你推了一个特征方程的通项公式…… 所以考虑log的做法：矩阵快速幂 如果你还不知道矩阵快速幂是什么，请走这边：传送门 对于这道题，嗯，模板嘛，已经告诉你了式子，就只需要考虑矩阵了，对于整个过程，我们
矩阵加速 （线性递推）
最新发布
qijingbo_的博客
05-26 1389
矩阵加速算法
星星之火OIer：矩阵加速例题（一）——TR的数列
qq_43890190的博客
04-09 285
在之前的讲解中 我们已经初步了解了矩阵加速 在这一讲中，我们就来更深入地了解一下矩阵加速 斐波那契数列——番外番外篇——TR的数列 题目描述： 这道题 首先确定原数组：： 然后我们定义 那么我们看怎么用和推出：： 又因为：： 所以：： 自此，我们就可以定出原数组：： 即 然后就可以推出这道题的加速数组：： 到需要再加即原数组的第三项 所...
[矩阵乘法]数列
正汰的学习笔记
10-12 1458
嘿嘿嘿…….矩阵乘法！！ https://yunpan.cn/cvxtdIUiCITNk 访问密码 38b3 细细~具体就是看PPT反正不是我写的~ OJ 【题目描述】 a[1]=a[2]=a[3]=1 a[x]=a[x-3]+a[x-1] (x>3) 求a数列的第n项对1000000007（10^9+7）取余的值。 【输入格式】 第一行一个整数T，表示询问个数。
C++矩阵加速例题斐波拉契数列变式—————TR的数列
C2020lax的博客
04-04 308
题目描述： TR非常喜欢数学，经常一个人拿出草稿纸研究奇奇怪怪的数学问题，最近，他突然对数列产生了兴趣，他找到一个数列，类似于斐波拉契，即：Tn=1*f1+2*f2+3*f3+……+n*fn (fn为斐波拉契的第n项值) 现在TR想请你帮忙求Tn%m的值 输入： 两个用空格隔开的整数n和m 1≤n,m≤-1 输出： Tn mod m的值 输入样例： 5 5 输出样例： ...
斐波那契数列（数论+矩阵加速）
Stockholm_Sun的博客
10-11 1138
斐波那契数列输入一个long long范围的一个数x，要求输出斐波那契数列第x项。思路显然这个题目，用O（N）的朴素做法是肯定不行的。所以我们要寻求一种新的方法。 接下来介绍一种计算：矩阵加法。这样的一种运算具有满足结合律的性质。对这个有兴趣的同学可以参阅人教版数学选修4-2的第一章和第二章（第一章带你入门，第二章教你矩阵乘法）。 大概矩阵加速就是和快速幂类似的，下面代码里面的操作可以教会我们其
HDU 3306
weixin_30437847的博客
09-01 198
先转一些 http://www.cnblogs.com/frog112111/archive/2013/05/19/3087648.html Fibonacci数列：F(0)=1 , F(1)=1 , F(n)=F(n-1)+F(n-2) 我们以前快速求Fibonacci数列第n项的方法是 构造常系数矩阵 （一） Fibonacci数列f[n]=f[n-1]+f[n-2],f[1]=...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值