码迷,mamicode.com
首页 > 其他好文 > 详细

(矩阵快速幂)51NOD 1242斐波那契数列的第N项

时间:2018-11-26 00:09:57      阅读:170      评论:0      收藏:0      [点我收藏+]

标签:等于   mod   class   快速   inpu   icon   lock   text   tle   

斐波那契数列的定义如下:
 
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的结果即可。
 

输入

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

输出

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

输入样例

11

输出样例

89
解:由于斐波那契数列的第N(N>2)项等于N-1个{{1,1},{1,1}}矩阵相乘后的第一项。
  由于这种矩阵形式上的特殊性(对称,乘法可交换),我们可以借助快速幂的思想可以快速求解这个答案。
 1 #include <stdio.h>
 2 
 3 #define MOD 1000000009
 4 
 5 int main()
 6 {
 7     long long n;
 8     while (scanf_s("%lld", &n) != EOF)
 9     {
10         long long a[2][2] = { 1,0,0,1 }, tmp[2][2] = { 1,1,1,0 };
11         if (n < 2)printf("%d\n", n);
12         else
13         {
14             --n;
15             while (n)
16             {
17                 if (n % 2)
18                 {
19                     int q, w, e;
20                     q = (tmp[0][0] * a[0][0] + tmp[0][1] * a[1][0]) % MOD;
21                     w = (tmp[0][0] * a[0][1] + tmp[0][1] * a[1][1]) % MOD;
22                     e = (tmp[1][0] * a[0][1] + tmp[1][1] * a[1][1]) % MOD;
23                     a[0][0] = q;
24                     a[0][1] = a[1][0] = w;
25                     a[1][1] = e;
26                 }
27                 int q, w, e;
28                 q = (tmp[0][0] * tmp[0][0] + tmp[0][1] * tmp[1][0]) % MOD;
29                 w = (tmp[0][0] * tmp[0][1] + tmp[0][1] * tmp[1][1]) % MOD;
30                 e = (tmp[1][0] * tmp[0][1] + tmp[1][1] * tmp[1][1]) % MOD;
31                 tmp[0][0] = q;
32                 tmp[0][1] = tmp[1][0] = w;
33                 tmp[1][1] = e;
34                 n >>= 1;
35             }
36             printf("%d\n", a[0][0]);
37         }
38     }
39 }

 

(矩阵快速幂)51NOD 1242斐波那契数列的第N项

标签:等于   mod   class   快速   inpu   icon   lock   text   tle   

原文地址:https://www.cnblogs.com/Ekalos-blog/p/10017725.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!