入门训练 Fibonacci数列

问题描述
Fibonacci数列的递推公式为：Fn=Fn-1+Fn-2，其中F1=F2=1。当n比较大时，Fn也非常大，现在我们想知道，Fn除以10007的余数是多少。
输入格式
输入包含一个整数n。
输出格式
输出一行，包含一个整数，表示Fn除以10007的余数。
说明：在本题中，答案是要求Fn除以10007的余数，因此我们只要能算出这个余数即可，而不需要先计算出Fn的准确值，再将计算的结果除以10007取余数，直接计算余数往往比先算出原数再取余简单。
样例输入10
样例输出55
样例输入22
样例输出7704
数据规模与约定1 <= n <= 1,000,000。

对于这道题，由于数据较大，所以用递归或者for循环，肯定会超时的.网上有很多优化的递归办法这里就不写了。这里提供一种快速的矩阵快速幂解决办法。在存在递推和形如 Fibonacci数列的式子的时候可以直接拿来用。只要会够造矩阵，使用矩阵快速幂就不是问题。
f[n-2],f[n-1]]*A = [f[n-1],f[n]]怎样构造A呢？
首先构造f[n-1].通过等式左边的式子可以观察到f[n-2]*0+f[n-1]*1=f[n-1],所以二维矩阵A的第0列元素为（0 1）。由于f[n]=f[n-2]+f[n-1]得f[n]=f[n-2]*1+f[n-1]*1这样就得到就矩阵的第二列元素（1，1）。这样就构造好了矩阵A既代码main中的temp。

#include<stdio.h>
#include<string.h>
unsigned long long int n;
struct fib{
    long long a[2][2];
};
fib multiply(fib x, fib y){
    fib temp;
    int i,j,k;
    memset(temp.a,0,sizeof(temp.a));
    for(i=0;i<2;i++){
        for(j=0;j<2;j++){
            for(k=0;k<2;k++){
                temp.a[i][j]=(temp.a[i][j]+x.a[i][k]*y.a[k][j]%10007)%10007;
            }
        }
    }
    return temp;
}
int main(void){
    int i,j;
    fib temp={0,1,1,1};//[f[n-2],f[n-1]]*A = [f[n-1],f[n]]从这个式子我们可以构造出A既temp 
    fib ans;            
    long long f[2]={1,1},sum=0;
    scanf("%llu",&n);
    memset(ans.a,0,sizeof(ans.a));
    for(i=0;i<2;i++){
        ans.a[i][i]=1;//单位矩阵，根据矩阵性质单位矩阵乘以任何一个矩阵，矩阵不变 
    }
    if(n==1 || n==2){
        printf("%d\n",1);
        return 0;
    }
    n=n-2;
    while(n){//while循环是矩阵快速幂的模板 
        if(n & 1){
            ans=multiply(ans,temp);
        }
        n= n>> 1;
        temp=multiply(temp,temp);
    }
    for(i=0;i<2;i++){
        sum=(sum+f[i]*ans.a[i][1]%10007)%10007;//这个式子的构造可以在草稿纸上写一个二维矩阵推导 
    }
    printf("%d",sum);
    return 0;

}