斐波那契数列-矩阵乘法
来源:互联网 发布:阿里云计算机有限公司 编辑:程序博客网 时间:2024/05/24 06:31
先简单介绍一下矩阵乘法求斐波那契数列的原理
问题的求解就变成的解决,而幂的求可用二分法来求。
int a[2][2];
mul2(result, base);
f(n) 是第n项的值。
f(1)= 1; f(2) =1;
f(n)= f(n-1) + (n-2)
下面的介绍是我从网上查到了,收益匪浅。
分两步推导:
二分法可用递归和非递归来求:
下面是代码:
定义矩阵
struct matrix //定义2*2的矩阵
{
};
void mul(matrix& x,matrix& y) // 矩阵乘法 x = x*y
{
matrix ans;
ans.a[0][0]=x.a[0][0]*y.a[0][0]+x.a[0][1]*y.a[1][0];
ans.a[0][1]=x.a[0][0]*y.a[0][1]+x.a[0][1]*y.a[1][1];
ans.a[1][0]=x.a[1][0]*y.a[0][0]+x.a[1][1]*y.a[1][0];
ans.a[1][1]=x.a[1][0]*y.a[0][1]+x.a[1][1]*y.a[1][1];
x = ans;
}
//下面这个函数是我从网上查到的
matrix power(matrix x,long e)
{
matrix ans,tmp;
if(e==0) //指数为0的时候,返回单位矩阵
{
ans.a[0][0]=1;
ans.a[0][1]=0;
ans.a[1][0]=0;
ans.a[1][1]=1;
return ans;
}
if( e==1 )
return x;
tmp=power(x,e>>1); //X的e/2次方
ans=mul(tmp,tmp);
if( e&1 ) //如果e为奇数,还要再乘以X
ans=mul(ans,x);
return ans;
}
非递归代码:
matrix result = {{1,0,0,1}};
matrix base ={{1,1,1,0}};
if(n&1) result = base;
n>>=1;
while(n)
{
mul2(base, base);
if(n&1)
{
}
n>>=1;
}
<<result.a[0][1]就是第 n项 的值。
0 0
- 斐波那契数列-矩阵乘法
- SSL 1529 斐波那契数列Ⅱ 矩阵乘法
- SSL 1530 斐波那契数列Ⅲ 矩阵乘法
- SSL 1531 斐波那契数列Ⅳ 矩阵乘法
- 洛谷1349 广义斐波那契数列 【矩阵乘法】
- Codevs 1574 广义斐波那契数列(矩阵乘法)
- 【矩阵乘法】广义斐波那契数列
- 矩阵乘法解fibonacci斐波那契数列
- 斐波那契数列三种实现+矩阵乘法+矩阵cimi
- 矩阵乘法 求斐波那契数列
- 矩阵乘法 求斐波那契数列
- 矩阵乘法求斐波那契数列
- 矩阵乘法求斐波那契数列
- HDOJ 4549 - M斐波那契数列 费马小定理,矩阵乘法
- [NOI2011]兔农(斐波那契数列+乘法逆+矩阵加速)
- 51nod--1242 斐波那契数列第N项 (矩阵乘法优化)
- 斐波那契数列的第 n 项 mod 1000000007(矩阵乘法)
- hdu 4549 M斐波那契数列(矩阵乘法+降幂公式)
- SDUTOJ 2712 5-2 派生类的构造函数
- Update视图数据(oracle)
- 排序算法
- VS2012 快捷键
- Run Loop详解
- 斐波那契数列-矩阵乘法
- UIImageView userInteractionEnabled
- Android的px、dp和sp
- 【诚邀各位开发者参加】2014-10-25 深圳OGEngine手游开发者交流会议
- 命令行执行java程序
- find命令【linux】
- [Leetcode] Remove Duplicates from Sorted Array II
- 程序员需要人脉
- 为何需要静音检测