矩阵快速幂--理解

原理：

斐波那契数列定义

用矩阵形式来描述递推关系：

以上递推式可以写成矩阵形式：

观察结果

例如F(2) = 1, F(3) = 2...

 

而实际上f(n)就是矩阵的左上角元素就是斐波那契数列的第n项的值

例题：

1242 斐波那契数列的第N项 

基准时间限制：1 秒 空间限制：131072 KB 分值: 0 难度：基础题

 

斐波那契数列的定义如下：

F(0) = 0

F(1) = 1

F(n) = F(n - 1) + F(n - 2) (n >= 2)

(1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...)

给出n，求F(n)，由于结果很大，输出F(n) % 1000000009的结果即可。

Input

输入1个数n(1 <= n <= 10^18)。
Output

输出F(n) % 1000000009的结果。
Input示例

11
Output示例

89

 代码实现：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9+9;
const int N=2;
ll i,j,k;
ll tmp[N][N],res[N][N];
void multi(ll a[][N],ll b[][N],int n) {
	
	memset(tmp,0,sizeof(tmp));
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			for(k=0;k<n;k++){
				tmp[i][j]+=(a[i][k]*b[k][j])%mod;
			}
		}
	}
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			a[i][j]=tmp[i][j];
		}
	}
}

void pow(ll a[][N],ll m,int n){
	memset(res,0,sizeof(res));
	for(i=0;i<n;i++) res[i][i]=1;
	while(m){
	if(m&1)       //m为奇数 
	multi(res,a,n);
	multi(a,a,n);
	m/=2;
  }
}

int main()
 {
 	ll m;
 	ll a[N][N];
 	int n;
 	while(scanf("%lld",&m)!=EOF){
 		 n=2;
		 a[0][0]=1,a[0][1]=1,a[1][0]=1,a[1][1]=0;
		 pow(a,m-1,n);
		 printf("%lld\n",res[0][0]);
	 }
	return 0;
}