hdu 3117 斐波那契数列

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3117

只会求后四位的，前四位也是看别人的博客学的。。。。

http://blog.csdn.net/xieqinghuang/article/details/7789908

#include<iostream>
#include<cmath>
const int m=10000;
const double t=(1+sqrt(5.0))/2;
using namespace std;
int n;
int f[40]={0,1,1,2,3,5,8,13,21,34,55};

struct Matrix{
    int map[2][2];
};

Matrix mat;

Matrix Mul(Matrix &a,Matrix &b){
    Matrix c;
    for(int i=0;i<2;i++){
        for(int j=0;j<2;j++){
            c.map[i][j]=0;
            for(int k=0;k<2;k++){
                c.map[i][j]+=a.map[i][k]*b.map[k][j];
                c.map[i][j]%=m;
            }
        }
    }
    return c;
}

Matrix Pow(int n){
    if(n==1)return mat;
    else if(n&1){
        return Mul(mat,Pow(n-1));
    }else {
        Matrix temp=Pow(n>>1);
        return Mul(temp,temp);
    }
}

int main(){
    while(scanf("%d",&n)!=EOF){
        if(n<40){
            for(int i=11;i<40;i++){
                f[i]=f[i-1]+f[i-2];
            }
            printf("%d\n",f[n]);
            continue;
        }
        mat.map[0][0]=1;
        mat.map[0][1]=1;
        mat.map[1][0]=1;
        mat.map[1][1]=0;
        mat=Pow(n);
        double ans=log10(1.0/sqrt(5.0))+n*log10(t);
        ans=ans-int(ans);//减去整数部分，得到小数
        ans=pow(10.0,ans);
        ans*=1000;
        printf("%d...",(int)ans);
        int temp=mat.map[0][1];
        if(temp<10){
            printf("000%d\n",temp);
        }else if(temp<100){
            printf("00%d\n",temp);
        }else if(temp<1000){
            printf("0%d\n",temp);
        }else {
            printf("%d\n",temp);
        }
    }
    return 0;
}