练习赛一 C 斐波那契进阶

C 斐波那契进阶

时间限制：1000ms   内存限制：65536kb

通过率：56/60 (93.33%)    正确率：56/104 (53.85%)

题目描述

AlvinZH听说大家都掌握了斐波那契数列，先简单了解一下，斐波那契整数序列中，有F(0)=0，F(1)=1，F(n)=F(n-1)+F(n-2) （n≥2）。这样我们可以快速地求出斐波那契数列的前几项：0,1,1,2,3,5,8,13,21,34,...

你知道吗？斐波那契序列还有另外一种定义：

AlvinZH不满足于此，他想知道当n很大很大时，F(n)有多大，你来告诉他吧！（结果对10007取模）

输入

输入将由多组测试数据组成，以EOF结尾。

每组数据只含一行，为题中描述的整数n（0≤n≤10^9）。

输出

对于每组数据，输出一行，为F(n)对10007取模的值。

输入样例

0
1
1000000000

输出样例

0
1
7300

解析：

快速幂求解斐波那契数列。附上两个链接：

1. 求解斐波那契数列复杂度分析

2. 矩阵快速幂求斐波那契数列（初学整理）

代码：

#include<cstdio>
#define M 10007
struct Matrix{
    long long v[3][3];
}mat;

Matrix mtMul(Matrix A, Matrix B)
{
    Matrix C;
    C.v[0][0]=(A.v[0][0]*B.v[0][0]+A.v[0][1]*B.v[1][0])%M;
    C.v[0][1]=(A.v[0][0]*B.v[0][1]+A.v[0][1]*B.v[1][1])%M;
    C.v[1][0]=(A.v[1][0]*B.v[0][0]+A.v[1][1]*B.v[1][0])%M;
    C.v[1][1]=(A.v[1][0]*B.v[0][1]+A.v[1][1]*B.v[1][1])%M;
    return C;
}

Matrix mtPow(Matrix A, int k)
{
    if(k ==1)
        return A;
    A=mtPow(A,k/2);
    if(k%2== 0)
    {
        return mtMul(A,A);
    }
    else{
        return mtMul(mtMul(A, A), mat);
    }
}

int main(){
    int n;
    mat.v[0][0]=1,mat.v[0][1]=1;
    mat.v[1][0]=1,mat.v[1][1]=0;
    while(scanf("%d",&n)!=EOF)
    {
        Matrix A = mtPow(mat,n);
        printf("%d\n",A.v[0][1]);
    }
    return 0;
}