51nod - 1242 斐波那契数列的第N项

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

 

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

 思路：

构造矩阵来根据F(n-1), F(n-2)得到F(n), F(n-1) 。然后利用矩阵快速幂来计算。

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
#include <stdio.h>
using namespace std;
 
typedef long long int ll;
const int MOD = 1e9 + 9;
struct Node
{
	ll v[2][2];
};
 
ll n;
 
Node mul(Node a, Node b)
{
	Node result;
	result.v[0][0] = ((a.v[0][0] * b.v[0][0]) % MOD + (a.v[0][1] * b.v[1][0]) % MOD) % MOD;
	result.v[0][1] = ((a.v[0][0] * b.v[0][1]) % MOD + (a.v[0][1] * b.v[1][1]) % MOD) % MOD;
	result.v[1][0] = ((a.v[1][0] * b.v[0][0]) % MOD + (a.v[1][1] * b.v[1][0]) % MOD) % MOD;
	result.v[1][1] = ((a.v[1][0] * b.v[0][1]) % MOD + (a.v[1][1] * b.v[1][1]) % MOD) % MOD;
 
	return result;
}
 
Node pow(ll n)
{
	Node base = {1, 1, 1, 0};
	Node result = {1, 0, 0, 0};
	while (n)
	{
		if (n & 1)
		{
			result = mul(result, base);
		}
 
		n >>= 1;
		base = mul(base, base);
	}
 
	return result;
}
 
int main()
{
	cin >> n;
	if (n == 1)
	{
		cout << 1 << endl;
		return 0;
	}

	Node a = pow(n - 1);
	ll result = a.v[0][0];
	cout << result << endl;
 
    return 0;
}