51nod 1242

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

上面的公式就是n项的值  

首先用到了快速幂 ，不会快速幂的去百度百科看  

然后就是矩阵的乘法

#include <iostream>
#include <cstring>
#include <stdio.h>
#include <cmath>
using namespace std;
const long long mod=1000000009;
int main()
{
    long long n;
    while(cin>>n)
    {

        long long d[4][4];
        long long p[4][4];
        long long s[4][4];
        int a[100];
        memset(d,0,sizeof(d));
        memset(p,0,sizeof(p));
        memset(s,0,sizeof(s));
        memset(a,0,sizeof(a));
        long long w=n;
        int j=0,i=0;
        while(w>0)
        {
            if(w%2==1)
            {
                a[j]=i;
                //cout<<a[j]<<' '<<j<<endl;
                j++;
            }
            i++;
            w=w/2;
        }
       // cout<<i<<endl;

        int l=j;


        for(i=0;i<l;i++)
        {
                d[1][1]=1;
                d[1][2]=1;
                d[2][1]=1;
                d[2][2]=0;
                for(j=0;j<a[i];j++)
                {
                    p[1][1]=d[1][1]*d[1][1] %mod+d[1][2]*d[2][1] %mod; p[1][2]=d[1][1]*d[1][2] %mod+d[1][2]*d[2][2] %mod;
                    p[2][1]=d[2][1]*d[1][1] %mod+d[2][2]*d[2][1] %mod; p[2][2]=d[2][1]*d[1][2] %mod+d[2][2]*d[2][2] %mod;
                    d[1][1]=p[1][1]%mod;
                    d[1][2]=p[1][2]%mod;
                    d[2][1]=p[2][1]%mod;
                    d[2][2]=p[2][2]%mod;
                }
                /*
            cout<<"sssssss"<<endl;
            cout<<s[1][1]<<' '<<s[1][2]<<endl;
            cout<<s[2][1]<<' '<<s[2][2]<<endl;
            cout<<"ddddd"<<endl;
            cout<<d[1][1]<<' '<<d[1][2]<<endl;
            cout<<d[2][1]<<' '<<d[2][2]<<endl;
            */
            if(i==0)
            {
                s[1][1]=d[1][1];
                s[1][2]=d[1][2];
                s[2][1]=d[2][1];
                s[2][2]=d[2][2];
            }
            else
            {
                p[1][1]=s[1][1]*d[1][1] %mod+s[1][2]*d[2][1] %mod;p[1][2]=s[1][1]*d[1][2] %mod+s[1][2]*d[2][2] %mod;
                p[2][1]=s[2][1]*d[1][1] %mod+s[2][2]*d[2][1] %mod;p[2][2]=s[2][1]*d[1][2] %mod+s[2][2]*d[2][2] %mod;
                    s[1][1]=p[1][1]%mod;
                    s[1][2]=p[1][2]%mod;
                    s[2][1]=p[2][1]%mod;
                    s[2][2]=p[2][2]%mod;
            }
            /*

            cout<<"sssssss"<<endl;
            cout<<s[1][1]<<' '<<s[1][2]<<endl;
            cout<<s[2][1]<<' '<<s[2][2]<<endl;
            cout<<"ddddd"<<endl;
            cout<<d[1][1]<<' '<<d[1][2]<<endl;
            cout<<d[2][1]<<' '<<d[2][2]<<endl;
            */


        }
        cout<<s[2][1]%mod<<endl;
    }
}

会快速幂和矩阵相乘就是挺水的一个题