描述
斐波那契数列的定义如下:
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
题解
因为n范围巨大,所以只用递归求是不行的,可以使用矩阵原理。
代码
#include <iostream>
using namespace std;
#define mod(a, m) ((a) % (m) + (m)) % (m)
const int MOD = 1e9 + 9;
struct MATRIX
{
long long a[2][2];
};
MATRIX a;
long long f[2];
void ANS_Cf(MATRIX a)
{
f[0] = mod(a.a[0][0] + a.a[1][0], MOD);
f[1] = mod(a.a[0][1] + a.a[1][1], MOD);
return ;
}
MATRIX MATRIX_Cf(MATRIX a, MATRIX b)
{
MATRIX ans;
int k;
for (int i = 0; i < 2; i++)
{
for (int j = 0; j < 2; j++)
{
ans.a[i][j] = 0;
k = 0;
while (k < 2)
{
ans.a[i][j] += a.a[k][i] * b.a[j][k];
ans.a[i][j] = mod(ans.a[i][j], MOD);
++k;
}
}
}
return ans;
}
MATRIX MATRIX_Pow(MATRIX a, long long n)
{
MATRIX ans;
ans.a[0][0] = 1;
ans.a[1][1] = 1;
ans.a[0][1] = 0;
ans.a[1][0] = 0;
while (n)
{
if (n & 1)
{
ans = MATRIX_Cf(ans, a);
}
n = n >> 1;
a = MATRIX_Cf(a, a);
}
return ans;
}
int main()
{
long long n;
while (cin >> n)
{
if (n == 1)
{
cout << '1' << '\n';
continue;
}
a.a[0][0] = a.a[0][1] = a.a[1][0] = 1;
a.a[1][1] = 0;
a = MATRIX_Pow(a, n - 2);
ANS_Cf(a);
cout << f[0] << '\n';
}
return 0;
}