POJ3070：Fibonacci——题解

http://poj.org/problem?id=3070

题目大意：求Fibonacci数列第n项，对10000取模。

矩阵乘法板子题……实在不知道写什么了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
ll n,m=;
struct node{
    ll g[][];
}f,res;
void buildI(node &x){//构造单位矩阵
    for(int i=;i<=;i++){
        for(int j=;j<=;j++){
            if(i==j)x.g[i][j]=1LL;
            else x.g[i][j]=0LL;
        }
    }
    return;
}
void multi(node &x,node &y,node &z){//z=x*y
    memset(z.g,,sizeof(z.g));
    for(int i=;i<=;i++){
        for(int j=;j<=;j++){
            if(x.g[i][j]){
                for(int k=;k<=;k++){
                    z.g[i][k]+=x.g[i][j]%m*y.g[j][k]%m;
                    z.g[i][k]%=m;
                }
            }
        }
    }
    return;
}
void qpow(ll k){
    buildI(res);
    node tmp=f,t;
    while(k!=){
        if(k&){
            multi(res,tmp,t);
            res=t;
        }
        multi(tmp,tmp,t);
        tmp=t;
        k>>=;
    }
    return;
}
ll solve(){
    if(n==)return 0LL;
    if(n<=)return 1LL;
    qpow(n-);
    ll ret=res.g[][]%m+res.g[][]%m;
    return ret%m;
}
int main(){
    while(scanf("%lld",&n)!=EOF){
        if(n==-)break;
        f.g[][]=;f.g[][]=;
        f.g[][]=;f.g[][]=;
        ll res=solve();
        printf("%lld\n",res);
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++