Codevs1732 Fibonacci数列2

最新推荐文章于 2018-02-13 21:53:03 发布

Leo_CT

最新推荐文章于 2018-02-13 21:53:03 发布

阅读量286

点赞数

分类专栏： OI-数论文章标签：矩阵乘法快速幂数论

本文链接：https://blog.csdn.net/Chorolop/article/details/72834177

版权

OI-数论专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题目

题目描述 Description
在“1250 Fibonacci数列”中，我们求出了第n个Fibonacci数列的值。但是1250中，n<=109。现在，你的任务仍然是求出第n个Fibonacci数列的值，但是注意：n为整数，且1 <= n <= 100000000000000

输入描述 Input Description
输入有多组数据，每组数据占一行，为一个整数n(1 <= n <= 100000000000000)

输出描述 Output Description
输出若干行。每行输出第（对应的输入的）n个Fibonacci数（考虑到数会很大，mod 1000000007）

样例输入 Sample Input
3
4
5

样例输出 Sample Output
2
3
5

数据范围及提示 Data Size & Hint
1 <= n <= 100000000000000

代码

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

struct M{
    long long mat[5][5];
    void INIT(){
        memset(mat,0,sizeof(mat));
    }

    M operator * (M B) const{
        M C;
        for(int i = 1;i <= 2;i++){
            for(int j = 1;j <= 2;j++){
                C.mat[i][j] = 0;
                for(int k = 1;k <= 2;k++)
                    C.mat[i][j] = (C.mat[i][j]+(mat[i][k]*B.mat[k][j]))%1000000007;
            }
        }
        return C;
    }
};

M ksm(M A,long long n){
    M B = A;
    while(n){
        if(n&1)
            B = B * A;
        A = A*A;
        n >>= 1;
    }
    return B;
}

int main(){
    long long q;
    while(cin >> q)
    {
        if(q == 1 || q == 2){
            cout << 1 << endl;
            continue;
        }
//      scanf("%d",&q);
        M A;
        A.mat[1][1] = A.mat[1][2] = A.mat[2][1] = 1;
        A.mat[2][2] = 0;
        A = ksm(A,q-2);
        printf("%lld\n",A.mat[1][1]%1000000007);
    }
    return 0;
}