矩阵快速幂及斐波那契数列模板

本篇博客先给出矩阵快速幂以及利用矩阵快速幂求斐波那契数列的模板，讲解待更新……

const int N=10;
int tmp[N][N];
void multi(int a[][N],int b[][N],int n)
{memset(tmp,0,sizeof tmp);for(int i=0;i<n;i++)for(int j=0;j<n;j++)for(int k=0;k<n;k++)tmp[i][j]+=a[i][k]*b[k][j];for(int i=0;i<n;i++)for(int j=0;j<n;j++)a[i][j]=tmp[i][j];
}
int res[N][N];
void Pow(int a[][N],int n)
{memset(res,0,sizeof res);//n是幂，N是矩阵大小for(int i=0;i<N;i++) res[i][i]=1;while(n){if(n&1)multi(res,a,N);//res=res*a;复制直接在multi里面实现了；multi(a,a,N);//a=a*an>>=1;}
}

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2;
const int MOD=1000000009;
struct mat
{ll a[N][N];
};
mat mat_mul(mat x,mat y)
{mat res;memset(res.a,0,sizeof(res.a));for(int i=0;i<2;i++)for(int j=0;j<2;j++)for(int k=0;k<2;k++)res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;return res;
}
void mat_pow(ll n)
{mat c,res;c.a[0][0]=c.a[0][1]=c.a[1][0]=1;c.a[1][1]=0;memset(res.a,0,sizeof(res.a));for(int i=0;i<2;i++) res.a[i][i]=1;while(n){if(n&1) res=mat_mul(res,c);c=mat_mul(c,c);n=n>>1;}printf("%I64d\n",res.a[0][1]);
}
int main()
{ll n;scanf("%lld",&n);mat_pow(n);return 0;
}

转载于:https://www.cnblogs.com/aerer/p/9930934.html

矩阵快速幂及斐波那契数列模板相关推荐

矩阵快速幂求斐波那契数列（初学整理）
参考文章: http://blog.csdn.net/u013795055/article/details/38599321 http://blog.csdn.net/g_congratulation ...
A - 【深蓝】教科书般的亵渎 (Gym - 102072A )（矩阵快速幂+计算斐波那契数列）
"现在给大家展示一波教科书般的亵渎" "扭了扭了" 炉石传说是一款比较火的卡牌游戏,不同的卡牌间可以打出许多不可思议的配合. 炉石传说的战斗大部分由随从来完成. ...
牛客18987 粉嘤花之恋（矩阵快速幂、斐波那契数列）
链接:https://ac.nowcoder.com/acm/problem/18987 来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 32768K,其他语言65536 ...
矩阵快速幂及斐波那契数列
1.通用模板:矩阵快速幂前几天写过一篇矩阵快速幂的一种解法,但是当时只是初学,所以用了当时我更好理解的方法去做.那个解法需要调用两个很相像的函数,而且容易造成紊乱,也不利于矩阵快速幂的推广,所以在观 ...
2017.3.25 矩阵快速幂求斐波那契数列第n项
对于矩阵快速幂只要知道矩阵取模.乘法原理就完全可以手推口诀:行列被计算的行列的交点是结果对应的位置: 剩下的就是推矩阵: 其实根据矩阵的方程意义就很好推了: 码: #includ ...
利用矩阵快速幂求斐波那契数列
我们知道如果用记忆化搜索逐项递推可以将复杂度降低到O(n),但是对于更大规模的输入,这个算法效率还是不够高,那么我们考虑更高效的算法: 二阶递推:f(n+2)=(1 1) f(n+1) f(n+1) ...
快速幂求解斐波那契数列
斐波那契数列斐波那契数列是很经典也很简单的一条题目.其满足: F n = { 1 ( n ≤ 2 ) F n − 1 + F n − 2 ( n ≥ 3 ) F_{n}= \begin{cases} ...
ZYH的斐波那契数列【线段树动态开点+矩阵快速幂求斐波那契】
描述 ZYH最近研究数列研究得入迷啦! 现在有一个斐波拉契数列(f[1]=f[2]=1,对于n>2有f[n]=f[n-1]+f[n-2]), 但是斐波拉契数列太简单啦,于是ZYH把它改成了斐波拉 ...
加强版斐波那契数列（矩阵快速幂）
关于快速幂的讲解可以参见我的上一篇博客<快速幂> 题目链接:又见斐波那契题目描述这是一个加强版的斐波那契数列. 给定递推式求F(n)的值,由于这个值可能太大,请对10 9+7取模. ...

矩阵快速幂及斐波那契数列模板

矩阵快速幂及斐波那契数列模板相关推荐

最新文章

热门文章